Probing neutron-hidden neutron transitions with the MURMUR experiment

MURMUR is a new passing-through-walls neutron experiment designed to constrain neutron/hidden neutron transitions allowed in the context of braneworld scenarios or mirror matter models. A nuclear reactor can act as a hidden neutron source, such that neutrons travel through a hidden world or sector. Hidden neutrons can propagate out of the nuclear core and far beyond the biological shielding. However, hidden neutrons can weakly interact with usual matter, making possible for their detection in the context of low-noise measurements. In the present work, the novelty rests on a better background discrimination and the use of a mass of a material - here lead - able to enhance regeneration of hidden neutrons into visible ones to improve detection. The input of this new setup is studied using both modelizations and experiments, thanks to tests currently performed with the experiment at the BR2 research nuclear reactor (SCK$\cdot$CEN, Mol, Belgium). A new limit on the neutron swapping probability p has been derived thanks to the measurements taken during the BR2 Cycle 02/2019A: $p<4.0 \ \times 10^{-10}$ at 95% CL. This constraint is better than the bound from the previous passing-through-wall neutron experiment made at ILL in 2015, despite BR2 is less efficient to generate hidden neutrons by a factor 7.4, thus raising the interest of such experiment using regenerating materials.

MURMUR is a new passing-through-walls neutron experiment designed to constrain neutronhidden neutron transitions allowed in the context of braneworld scenarios or mirror matter models. A nuclear reactor can act as a hidden neutron source, such that neutrons travel through a hidden world or sector. Hidden neutrons can propagate out of the nuclear core and far beyond the biological shielding. However, hidden neutrons can weakly interact with usual matter, making possible for their detection in the context of low-noise measurements. In the present work, the novelty rests on a better background discrimination and the use of a mass of a material -here lead -able to enhance regeneration of hidden neutrons into visible ones to improve detection. The input of this new setup is studied using both modelizations and experiments, thanks to tests currently performed with the experiment at the BR2 research nuclear reactor (SCK·CEN, Mol, Belgium). A new limit on the neutron swapping probability p has been derived thanks to the measurements taken during the BR2 Cycle 02/2019A: p < 4.0 × 10 −10 at 95% CL. This constraint is better than the bound from the previous passing-through-wall neutron experiment made at ILL in 2015, despite BR2 is less efficient to generate hidden neutrons by a factor 7.4, thus raising the interest of such experiment using regenerating materials.

I. INTRODUCTION
In many extensions of both the Standard Model (SM) of particle physics and the ΛCDM cosmological model, the visible sector of particles containing the usual baryonic matter is supplemented by a "hidden" sector. These models aim at addressing the shortcomings of the standard models, in particular the lack of a microscopic description of Dark Matter, or they aim at exploring the possible low-energy phenomenology associated with quantum gravity. A hidden sector could take the form of a duplication of the SM particle content, in which case each particle would have a hidden mirror partner [1][2][3]. Alternatively, a hidden sector could take a geometrical meaning in the braneworld hypothesis: if our visible world is confined on a brane inside a bulk of more than 4 spacetime dimensions, the bulk could also contain a parallel brane where hidden particles could travel [4][5][6][7][8][9]. Here, we consider the generic subclass of hidden sector models for which the fermions could exist in a hidden state as well as in the visible state. In particular the neutron n would have a hidden state n ′ . The hidden state is sterile as it does not interact via the SM interactions with the visible matter. However, the hidden state could mix with the visible state through mass mixing (in the theoretical context of mirror matter [1]), or through kinetic or geometric mixing (in the context of the parallel brane model [10]). The mixing induces n → n ′ (or n ′ → n) transitions. The low energy phenomenology (see [10] for a derivation in the context of the parallel brane model) is obtained by considering the evolution of the four-state quantum problem with the following Hamiltonian: Here, E v and E h describe the energy of the visible and hidden states, respectively, and the mixing is described by ǫ. Those are in principle 2 × 2 Hermitian matrices acting on the spin states, which we assume diagonal to simplify the following discussion without loss of generality. In the regime of quasi-degenerate states, i.e. |E v − E h | < ǫ, neutrons would undergo slow oscillations to the hidden state with large amplitude. Several measurements searching for resonant disappearance of ultracold neutrons have been reported [11][12][13][14][15][16][17], setting already stringent limits in the quasi-degenerate regime. On the contrary, in the regime of nondegenerate states, i.e. |E v − E h | ≫ ǫ, the probability for a visible neutron to be found in the hidden state oscillates rapidly, with a mean value of [18]: The non-degenerate regime can be probed through regeneration effects, i.e. n → n ′ → n transitions. First, a visible neutron would convert into a hidden neutron by scattering on nuclei in a converter (in practice the moderator material of a large neutron source). Indeed, a collision with a visible nucleus acts as a quantum measurement process and the neutron wave-function collapses either in the visible state with a probability (1 − p) or in the hidden state with a probability p. The cross section for the n → n ′ transition at each collision is given by [18]: where σ s (n → n) is the usual elastic scattering cross section. After the transition, the hidden neutron would propagate freely and escape the source. Outside the source, hidden neutrons could convert back to visible neutrons by scattering on nuclei in a regenerator, by using the reverse process n ′ → n, which has the same weak cross section as the direct process n → n ′ given by Eq. (3). Experiments based on n → n ′ → n transitions [18][19][20] are called neutron passing-through-wall experiments [18,19]. Some of us (M.S.,G.P.,J.L.,O.M.,G.T.) have conducted the first measurement of this kind at the Institut Laue Langevin (ILL, Grenoble, France) [18,19]. No significant regeneration of neutrons was observed, which set the limit p < 4.6 × 10 −10 (95% C.L.) [19]. The present paper describes a new passing-through-wall experiment performed at the BR2 reactor (SCK·CEN, Mol, Belgium). The novel detector, which we call MURMUR, separates the material to regenerate neutrons through n ′ → n transition from the material to detect the visible neutrons and is also designed to improve the background rejection. The paper is organized as follows. In section 2 we elaborate on the concept of passing-though-walls neutron experiments and recall how to evaluate the expected flux of regenerated neutrons at the detector location. Section 3 describes the BR2 reactor and the evaluation of the hidden neutron flux. Section 4 describes the MURMUR detector. In section 5 we discuss the neutron background and mitigation, and in section 6 we present the results of the first run of measurement.

II. GENERAL FRAMEWORK OF A PASSING-THROUGH-WALLS NEUTRON EXPERIMENT
A neutron has a nonzero probability p to convert from a visible state to a hidden one and reciprocally [18]. It is possible to induce disappearance (reappearance) toward (from) a hidden state thanks to nuclei which possess a high scattering cross section. The general setup of a passing-through-wall neutron experiment is constituted by a high neutron flux source, a converter, a regenerator and a neutron shielding wall, as shown in Fig. 1. The converter, made of a material with a high scattering cross section, is assumed to generate a flux of hidden neutrons in an hidden state from the incident high neutron flux. Then, a shielding wall makes possible to stop visible neutrons from the source, while hidden neutrons would interact very weakly with the visible matter of the wall and are thus free to pass through it. Finally, a certain proportion of the hidden neutron flux could be regenerated into a visible one thanks to a regenerator, also constituted by a material with a high scattering cross section and with a low absorption cross section. After having swapped back in a visible state, these regenerated neutrons could be detected thanks to a neutron detector, with a neutron detection efficiency ξ.
A nuclear reactor with a high neutron flux, which undergoes many elastic collisions with the moderator of the core, would be then a very intense source of hidden neutrons. In the first experiment of this kind in 2015 [19], there was no regenerator distinct from detector. A 33.5 cm 3 3 He proportional counter at 4 atm simultaneously played the role of regenerator and of neutron detector. In the present new experiment -called MURMUR -the regeneration of hidden neutrons and the detection of regenerated visible neutrons are separated. This improved concept of measurement results into an increase of the efficiency to detect hidden neutrons. Indeed, while the neutron detector remains the . Improved concept of a passing-through-wall neutron experiment. Such an experiment is constituted by a converter and a regenerator, both with high scattering cross sections, separated by a shielding wall necessary to shield the rest of the experiment from neutron background. The converter makes possible to transform one part of a visible neutron flux Φv into a hidden one Φ h , with an efficiency proportional to p, while the regenerator enhances their reconversion into visible neutrons, with an efficiency also proportional to p, to be detected by a neutron detector separated from the regenerator. ΓR is the regenerated neutron rate in the regenerator. A certain proportion ξ of the regenerated neutrons are detected by the neutron detector with a neutron counting rate ΓDet.
same than in the previous experiment, a low-cost lead matrix now surrounds the detector thus enhancing the neutron regeneration as detailed hereafter. Let us assume that our Universe is endowed with two mutually invisible sectors. One considers a neutron flux Φ v inside a nuclear core. Our aim is to determine the intensity of the hidden neutron flux Φ h , in the vacuum sector of the core neighbourhood. For a given volume element of the reactor, since each neutron collision allows for a hidden neutron conversion with probability p, the hidden neutron source is proportional to the macroscopic scattering cross section Σ s of the reactor element (first converter in Fig. 1) and to the local neutron flux Φ v . More specifically, we get the source term corresponding to the number of hidden neutrons per unit volume and unit of time (see reference [18]): Equation (4) is derived by using a matrix density approach and a two-state Universe description, as shown in ref. [18]. From this source term, we deduce the hidden neutron flux Φ h at the position r by integrating over the reactor volume V : The relation (5) shows that hidden neutron signal should decay as 1/D 2 when the distance D between the reactor and the detector increases and is valid as neutrons propagate in the vacuum of the hidden sector. Φ v is derived from ab initio computations using MCNP as explained later. We now want to build a detector which could detect the hidden neutron flux. The present concept is based on the use of a "regenerator" (see Fig. 1) able to convert the hidden neutron flux Φ h into a visible neutron flux Φ ′ v in order to detect it. The calculations of the equations introduced in this section are detailed in A using an approach similar to this of Ref. [18]. It can be shown that the visible neutron flux Φ ′ v regenerated by the regenerator follows: where D is now the diffusion coefficient and Σ a the absorption macroscopic cross section of the regenerator medium. The expression of the regenerated visible source term is given by (see A): where Φ h is the hidden neutron flux given by equation (5) and Σ ′ s the macroscopic scattering cross section of the regenerator. Lead has a density much higher than 3 He, while having a macroscopic scattering cross section (≈ 0.4 cm −1 ) close to the absorption cross section of 3 He (≈ 0.53 cm −1 ) and a low absorption cross section (≈ 5.6 × 10 −3 cm −1 ). Thanks to its high atomic mass, lead has also the advantage to thermalize very poorly epithermal and fast neutrons. As a result, lead is an excellent candidate as regenerating medium. However, epithermal and fast neutrons coming from the astrophysical and reactor backgrounds -although weakly thermalized in lead -still lead to a residual source of noise for the 3 He proportional counter. This is the reason why one deals with the background mitigation in section V. Moreover, regenerated neutrons can escape out of the lead matrix by scattering. Thus, only part of the regenerated neutrons can be detected by the 3 He detector located in the lead matrix (see Fig. 1), resulting then into a regenerated neutron detection efficiency ξ < 1. Equation (6) can be propagated thanks to MCNP6 numerical computations, as detailed in section VI, to obtain ξ in order to calculate the regenerated flux inside lead from the measurements of the 3 He detector. It is noteworthy that hidden neutrons can also been regenerated directly in the 3 He proportional counter [18], where a regeneration is here also a detection (ξ = 1).

III. THE BR2 NUCLEAR CORE AS VISIBLE NEUTRON/HIDDEN NEUTRON CONVERTER
The Belgian Reactor (BR2) is a highly heterogeneous high flux engineering research thermal reactor operated by the Belgian Nuclear Research Centre (SCK·CEN) in Mol, Belgium. This tank-in-pool reactor is cooled by light water in a compact high-enriched uranium (HEU) core (93% 235 92 U ). The moderator as well as the reflector of the reactor is a hyperboloid matrix of beryllium surrounding the HEU core (see Fig. 2 a). The beryllium matrix is constituted of a large number of skew irregular hexagonal prisms, which form together a twisted hyperbolic bundle around the central 200 mm channel H1 containing beryllium plug. Due to the particular shape of the beryllium matrix, there is no fixed core configuration, which allows an important flexibility for many different core arrangements. The beryllium matrix holds 79 reactor cylindrical channels: 64 standard channels (⊘84 mm), 10 reflector channels (⊘50 mm) and 5 large channels (⊘200 mm). All channels can contain fuel components, control rods, reflector plugs or experiments. The reactor can be operated at a power of 40 − 100 MW, usually about 150 to 180 full power days per year. The thermal neutron flux is about 1 × 10 15 cm −2 s −1 and the fast neutron (defined with an energy higher than 0.1 MeV) flux 6.0 × 10 14 cm −2 s −1 at the power of 60 MW. At the moment, the configuration being used for runs of the experiment is an extended core loaded with about 33 fuel elements and operated at power levels 55 MW to 100 MW with 6 or 8 shim-safety control rods, as seen in Fig. 2 b which gives the geometry of the mid-plan of the BR2 for the Cycle 02/2019A. The moderator, i.e the beryllium matrix, appears in yellow and the control rods in red. There are 9 control rods, 8 of which are hafnium and one is cadmium. A 3D model of the BR2 core was developed by the SCK·CEN team using the latest version of the Monte Carlo transport code MCNP6 [21,22]. The model is a complete 3D description of the BR2 hyperboloid nuclear core formed by twisted and inclined reactor channels (see Fig. 3). Each channel, fuel component, beryllium plug and experimental device is represented separately, with its own position and inclination. The MCNP6 simulation is coupled with CINDER90 for the isotope evolution part [21,22]. The credibility of the MCNP model was demonstrated by multiple comparisons of code predictions with available experimental data, such as control rod worth's, neutron fluxes, gamma heating, linear power and fission rates, reactivity effects, etc [22]. The thermal neutron flux (integrated neutron flux for neutron energies below 0.5 eV) of the Cycle 02/2019A (see the reactor configuration in Fig. 2 b) has been numerically calculated at the center of each hexagon of the mid-plan using MCNP model described just above. The results are summarized in Fig. 4. Thermal neutron flux in a hexagon vary with the Z coordinate. The thermal neutron flux in function of Z is considered constant between −20 cm and +4 cm and then assumed to drops linearly with Z. The hidden neutron can be calculated from equation (5). As seen in Fig. 2 b, the core has a complex geometry. One considers here regular hexagonal prisms as elementary core volumes dV , regularly spaced in height and where the thermal neutron flux is constant. One also considers that hexagonal channels are straight and not twisted and inclined (we consider this as negligible since the inclination angle is less than 10°). The numerical integration is made over the moderator regions of the core, i.e. over the beryllium matrix. Control rod regions, where the flux is the most important, are filled with beryllium in order to be conservative (the insertion of hafnium or cadmium which have a high scattering cross section in the core varies with the elapsed time of the cycle). The rest of the core's volume which is not control rods or the moderator are replaced by aluminium (which has a proper weak scattering cross section). The elementary volumes are calculated with the proper surface in the reference frame in Fig. 4. The light water pool surrounding the core behaves also as a visible neutrons/hidden neutrons converter. However, the neutron flux at this location being much weaker than in beryllium, its contribution can be neglected. . Beryllium -the moderator of the reactor -appears in yellow, control rods made of hafnium or cadmium are in red and fuel components in pink and cyan. There are 9 control rods and only the 9th control rod is made of cadmium. It is noteworthy that the core configuration can vary from one run to another. Indeed, the different channels can be filled with fuel, beryllium or samples as required.

IV. MURMUR EXPERIMENT: LEAD AS HIDDEN NEUTRON RENEGERATOR
The MURMUR experiment [20] is installed at 6.18 meters from the centre of the BR2 nuclear reactor of the SCK·CEN at Mol (Belgium) as shown in Fig. 5 where the global layout of the experiment near the nuclear core is described. The nuclear core is mainly constituted of beryllium and is surrounded by a light water tank for reactor cooling. A light water pool and an heavy concrete wall constitute the biological shield, which guarantees the security of the BR2 site against radiations. The low-noise MURMUR detector itself is detailed in Fig. 6. The beryllium matrix which plays the role of moderator for the BR2 nuclear core acts as a visible neutrons/hidden neutrons converter. The regenerator which regenerates hidden neutrons into visible ones is a 23.6 × 17.7 × 11.5 cm 3 matrix of lead, constituted of lead blocks with a full volume of 4389.4 cm 3 and a total mass of 49.82 kg. While having a weak absorption cross section, lead has a high macroscopic scattering cross section close to the macroscopic absorption cross section of 3 He (0.53 cm −1 ) used for the first experiment at the ILL in 2015 [18]. Thus, the hidden neutrons can be regenerated into visible neutrons in a large volume of lead, while being weakly absorbed. Regenerated visible neutrons are free to propagate in lead up to an cylindrical 3 He proportional counter of 33.5 cm 3 with a gas pressure of 4 atm located inside the matrix of lead, making possible the detection of hidden neutrons. Hidden neutrons can also  The BR2 nuclear core (1) provides the visible neutron flux Φv. A small part of this visible neutron flux is expected to be converted into hidden neutron flux Φ h , which can freely escape the light water pool (2) and the concrete wall (3) of the reactor. In a second configuration of the experiment, a wall made of paraffin (4) has been added in order to reduce the neutron background coming from the nuclear core. In this configuration, the experiment (5) is located at 6.18 meters from the centre of the BR2 nuclear core.
be regenerated directly in the 3 He counter while being detected at the same time. In addition, due to heavy nuclei, lead thermalizes very poorly neutrons compared with other materials as graphite or heavy water. This property is important, as epithermal or fast neutrons could succeed to pass through the wall shielding and would thermalize in the regenerator, representing then a possible source of noise in the 3 He proportional counter. All the regenerated neutrons in lead are not detected by the 3 He proportional counter, only a certain proportion ξ of them are captured by the 3 He gas. This proportion is numerically calculated in section VI thanks to MCNP6. A 4π boron carbide (B 4 C) box of 3.6 cm of thickness is used as shielding wall in order to isolate the lead and the 3 He proportional counter from the external environment. A paraffin shielding of 40 cm of thickness has also been added between the core concrete wall and the experiment in order to slow down epithermal and fast neutrons coming from the nuclear core. A plane plastic organic scintillator EJ200 of 50 × 50 cm 2 is finally placed at the top of the setup in order to veto the data acquisition. Indeed, as discussed in section V, cosmic muons can induce neutrons in lead, becoming a source of noise. The presence of the scintillator on the top of the matrix of lead allows both to monitor the cosmic particle flux and to prevent the data acquisition from the muon-induced neutrons. Signals from both detectors of the experiment ( 3 He counter and scintillator) are connected to a DT5790 digitizer of the CAEN firmware [23] including a flash ADC (analog-to-digital converter). The acquisition is triggered with a threshold on the signal amplitude and for reach triggered event, the charge is integrated. To allow particle discrimination, two different time windows are used for the integration. Indeed, for a given charge, neutron-captures which induce nuclear recoils (proton and tritium), make shorter tracks than gamma rays which produce electrons, by photoelectric and Compton effects. These differences makes possible to discriminate them according to the shape of their voltage pulses. The Pulse Shape Discrimination (P SD) is a number calculated online according to the expression: where Q L and Q S are respectively the long and the short integrated charges. Time windows have been chosen empirically in order to optimize the discrimination between γ rays and neutrons (see section V B). The long gate is 20 µs and the short gate is 1.5 µs.
The software creates two output files every hour, one for each detector ( 3 He and scintillator), which contain information about each measured pulse accessible for readout at any time by the user: the long charge, the short charge, the P SD and the treatment time of the pulse by the software in ms. An anti-coincidence system is also implemented in the software in order to veto the 3 He proportional counter. This veto is used to overcome the muons-induced neutrons in lead, possible source of noise of the experiment. Each time an event is measured in the 3 He proportional counter in a time window of 8 µs after an event detection in the scintillator, it is automatically suppressed by the software and does not appear in the list event output file of the 3 He detector.

V. BACKGROUND MITIGATION
The MURMUR detector must reach a very low rate of background: the lower the background, the better the constrain on the swapping probability p. This section describes the different possible sources of background and characterizes the efficiency of the experimental setup to deal with them, in the light of preliminary measurements.

A. Background description
The sources of background can be classify in two categories. The first one is related to the 3 He detector in the form of γ rays and intern α particles. The second one gathers the visible neutrons coming from the reactor, or be created in the atmosphere, in the reactor building and surrounding materials by spallation of cosmic muons or (α, n) reactions. As mentioned above, γ rays can interact in the detector by photoelectric effect or Compton scattering. Given the path length of electrons, γ rays are at low energy. However, a certain proportion of γ ray events, for example, events which have a longitudinal track, can overlap in the region of the neutron event energy spectrum, leading to an overestimate of the neutron rate in the 3 He detector. The α-particle emission from the actinides present as contaminants in counter materials is a well known issue of low-background experiments [24,25]. Their energy can vary from 0 to several MeV and overlaps in the region of the neutron capture products in the energy spectrum of 3 He counters. The α signal can thus not be distinguish from a neutron signal. The mean α emission rate depends on the nature of the material and is more important in aluminium (∼ 6 × 10 −4 − 3 × 10 −2 s −1 ) than in steel (∼ 10 −4 s −1 ) [25]. The 3 He proportional counter used for the experiment is made in steel but in reason of the very weak counting rate specific of passing-through-wall neutron experiment (∼ 10 −4 s −1 ), α emissions can be a source for an undesirable background noise. For this reason, measurements of the natural radioactivity of the 3 He proportional counter have been carried out with a HPGe detector in a low-noise system made of one ton of lead surrounded by plastic scintillators for the mitigation of the muon background in the HPGe detector [26,27]. We did not measure any difference of signal for the main γ rays coming from actinides with and without the 3 He counter inside the system, which means that the counter contains a very weak quantity of actinides. The intern α background of the counter can thus be neglected. Cosmogenic neutrons, which are mainly induced by muon spallation, constitute an important source of background for experiment looking for rare events, as hidden neutrons. In the case of MURMUR, cosmogenic neutrons can thus be induced in the surrounding of the experiment [28,29] or directly in the matrix of lead [30]. From the muon-induced neutron yield [28,29] and from the muon spectrum [31], one can roughly estimate that 0.5 neutron per second is produced in 50 kg of lead. These neutrons having a mean energy of 8 M eV , the probability they would be captured by the 3 He is weak. However, some of the lowest-energy neutrons produced by spallation could be thermalized enough in the lead to become thermal neutrons. Finally, all materials surrounding the experiment, as concrete and shielding, have the ability to produce neutrons under the effect of the cosmic muon flux. These fast neutrons could be energetic enough to pass through the shielding, loss enough energy inside it to reach the thermal energy and have good chance to be detected by the 3 He proportional counter. Finally, (α, n) reactions can also be a source of background due to the presence of the large quantity of concrete as biological shielding [32][33][34]. Figure 7 shows the P SD as a function of the measured charge from the 3 He counter taken for 15 days together with the reconstructed energy spectrum of the 3 He counter. It is clearly observed that the neutron spectrum is well reconstructed as we can observe the full energy peak at 764 keV as well as the p and 3 H escapes at 191 and 573 keV respectively. We can conclude that the events due to neutrons are clearly discriminated using a P SD selection in addition to the energy cut (see red box in Fig. 7). Table I (1-a), is in agreement with the measurement which was done to prepare the ILL experiment [19]. The addition of lead (1-b) increases the neutron signal in the 3 He detector of one order of magnitude due to muon spallation since the same measurement carried out using the active veto (1-c) comes back to the initial value. Muons coming mostly from the zenith [31], one has chosen to put a scintillator on the top only of the boron carbide box. The sides of the experiment are thus not protected by the veto, even though they profit from the umbrella effect of the edges of the scintillator. which the muons can interact by spallation. Therefore, we expect less muons, as seen with configurations (2-b) and (4-b) where the rates are due to muon-induced neutrons in the lead, but more neutrons as seen in configuration (4-a). The distribution of arrival direction is also expected to be less vertical. This explains the lower impact of the active veto. While the rates is reduced by a factor of 8 in Namur, it remains almost unchanged at BR2 (configurations (2-c) and (4-c)). During the ON period, we can notice a higher rate of neutrons in the 3 He counter (3-c) and also in the veto (3-d), coming from the reactor. This reactor-induced background justifies the addition of the paraffin block which allows, by thermalizing high energy neutrons, to reduce significantly the counting rates in both detectors (5-c and 5-d), With this configuration, there is almost no difference between OFF and ON periods. One can conclude that the low-noise efficiency of the MURMUR experiment is limited by the cosmic rays background which induces neutrons both in the environment and in the matrix of lead by the sides.

D. Background subtraction
Some sources of background still remain, as muon-induced neutrons which are not filtered by the veto and neutrons coming from (α, n) reactions. In the aim of removing their contributions, reactor ON and reactor OFF measurements will be realized thanks to BR2's shut-down periods. The neutron counting rate obtained during the OFF period of the reactor could be thus subtracted from the neutron counting rate of the reactor ON measurements. From the signal difference, we thus obtain an upper estimation of the neutron flux induced by the reactor, whatever its origin, possibly due to hidden neutron generation and reappearance. This flux makes possible to fix a new constraint on the swapping probability p thanks to numerical computations. Long acquisition times of about 30 days have been done in order to have enough statistic and thus to reduce the uncertainty on the measurements.

VI. ANALYSIS IN TERM OF NEUTRON-HIDDEN NEUTRON TRANSITIONS
This section describes the method used to derive the constraint on the neutron swapping probability p thanks to the MURMUR experiment measurements. From Eqs. (5) and (7), the regenerated neutron rate in lead is a function of p: where is numerically calculated from the visible neutron flux Φ v induced in the BR2 nuclear core and calculated in section III.
Only a certain proportion of hidden neutrons regenerated in lead will be detected by the 3 He proportional counter. Monte Carlo simulations are required to study the propagation of regenerated neutrons in lead and evaluate the proportion of them which are captured by the 3 He gas. In this purpose, we used MCNP6. The detection efficiency ξ depends on the neutron energy. Indeed, neutron capture cross section decreases with increasing energies. In order to ensure the relevance of the simulations, MCNP6 results have been confronted to experimental data. For this purpose, a calibration experiment has been carried out with an 3 He detector and a neutron source of 1.79(±0.04) × 10 3 Bq supplemented by a 50 × 40 × 10 cm 3 high density polyethylene block (HDPE) used to thermalize fast neutrons coming from the source. MCNP6 simulations have been performed with the whole setup and the addition of a thermal treatment card for the HDPE. The number of neutronic captures counted with MCNP6 has been compared to the number of detected neutrons after noise subtraction. Table II summarizes the results. The detection efficiency of regenerated neutrons has been calculated with MCNP6 against the neutron energy according to the experimental setup geometry. Figure 8 shows MCNP6 results for the regenerated neutron detection efficiency ξ against thermal neutron energy. In the thermal energy region (between 0.01 and 0.5 eV), the detection efficiency varies from 2.6 to 1 %. These values are then convoluted with the thermal energy spectrum of the BR2 [35] in order to obtain a mean detection efficiency of 1.64 ± 0.17%, where the uncertainty of 0.17% is considered in order to be conservative with table II. MCNP6 results have been confirmed thanks to Geant4 numerical computations with a discrepancy between the two software lower than 10 %. Table III gives the measurements used to derive a new constraint on the swapping probability p. 760 hours of measurement were considered during the Cycle 02/2019A of the BR2 nuclear reactor in April 2019. Two shut-down periods were considered to do the noise subtraction, respectively before and after the Cycle 02/2019A, i.e. 831 hours  of shut-down in December 2018/January 2019 and 998 hours in May/June 2019. The mean atmospheric pressure during these measurements are closed to each other and the counting rate in the 3 He counter during the two shut-down periods are statistically compatible, which makes possible to use the mean value of these two measurements to do the shut-down subtraction. The resulting neutron counting rate Γ Det in the 3 He proportional counter is then given by: which is better than reference [19] by a factor 4 and where confidence intervals have been calculated according to the method described in reference [36]. The numerical calculation of the regenerated neutron rate in lead from Eqs. (9) and (10) leads to: The BR2 visible neutron flux Φ v has been calculated as mentioned in section III and can be seen in Fig. 4. The numerical calculation method of the hidden neutron flux Φ h is described in section III. The value of the beryllium neutron scattering cross section is 6.2 b, which is a convoluted value according to the thermal neutron energy spectrum of the BR2 [35] (cross sections for different neutron energies come from END/F VII library). The length of the lead matrix (24 cm) is not enough to induce a significant variation of the hidden neutron flux along its propagation pathway.
The flux is then calculated at the center of the lead block and considered constant everywhere in lead. The neutron scattering cross section of the lead being constant in the the thermal energy range, it does not require to be averaged over the thermal neutron energy spectrum. The regeneration of hidden neutrons directly inside the 3 He counter is also considered. In this case, a regeneration is a detection, so ξ = 1. The numerical calculation of the regenerated hidden neutron in the 3 He counter by using a formula similar to Eq. (10), but where the macroscopic cross section of lead is replaced by the neutron capture cross section of the 3 He (0.53 cm −1 ) and the volume by the 33.5 cm 3 of gas, leads to: d3 He = 8.6 × 10 13 s −1 .
The constraint on the swapping probability p is then given by: Equations (11), (12), (13), and (14) finally lead to the following constraint on the neutron swapping probability: This constraint is slightly better than the one obtained in 2016 (p < 4.6 × 10 −10 ) [19] with the first experiment at the ILL (Grenoble , France) [18]. It must be underlined that the expected efficiency for the BR2 reactor to produce hidden neutron flux is weaker than the ILL research reactor. Indeed, the d3 He parameter of the ILL experiment was 6.4 × 10 14 s −1 , i.e. 7.4 times the value at BR2, which is mainly due to the larger effective volume of the reactor moderator at the ILL. The addition of the lead matrix as regenerator makes possible to increase the detected regenerated neutron flux by a factor equal to (V P b Σ P b ξ + V3 He Σ3 He )/(V3 He Σ3 He ) = 2.5 with respect to the first experiment at the ILL in 2016 [19]. While the BR2 nuclear reactor presents a lower ability to produce hidden neutrons, MURMUR makes possible to give a constraint better than the previous one found at the ILL [19] thanks to a lower counting rate and a regenerator made of lead. This result is thus promising for future low-noise passing-through-walls neutron experiments using regenerating materials as lead. Improvements of the MURMUR experimental setup are under consideration in order to reduce the background and make possible to get a better constraint on p (the constraint on p is as small as the neutron counting rate is, as expressed by Eq. (14)), while continuing to benefit from the long-term location provided by the BR2. MURMUR is assumed to stay in function at BR2 until 2021, which makes possible to test many other experimental configurations.

VII. CONCLUSIONS
In the last five years, some concerns have focused on the realization of passing-through-wall neutron experiments able to constrain the existence of hidden sectors thanks to disappearance/reappearance of neutrons in the surroundings of a nuclear core. MURMUR is a new and improved passing-through-wall neutron experiment, installed for now near the BR2 reactor of SCK·CEN at Mol (Belgium). Many improvements with respect to the first passing-through-wall neutron experiment carried out at the ILL (Grenoble, France) in 2015 have been achieved. The major ones are a lower neutron counting rate thanks to a background substraction, a pulse shape discrimination and the use of 50 kg of lead as a hidden neutron regenerator separated from the neutron detector. This makes possible to obtain a new constraint on the neutron swapping probability p better than the last one found at the ILL [19].
If we insert Eq. (A21) into Eqs.(A19) and (A20), we obtain: where we have set p = 2Ω 2 /η 2 [10]. Let us now consider an initial local neutron flux Φ 0 = u 0 v where u 0 is the initial neutron density. The neutron fluxes in each sector are given by Φ v,h = Φ 0 ρ v,h . Then, from Eq. (A24) we deduce: Since the neutron flux Φ h = u h v, where u h is the local neutron density in the hidden sector, then: Since u h is now local, Eq. (A26) must be supplemented by a divergence term ∇ · j v to account for the local behaviour of the neutron current j v = u v v. We then deduce the continuity equation for neutrons in the second sector: which is the continuity equation endowed with a source term S ′ v : Considering j v = −D∇Φ v in order to take into account the diffusion of the regenerated visible neutrons, we deduce the expected diffusion equation in the steady state regime (i.e. ∂ t u v = 0): Thanks to the quickly oscillating terms related to cos(ηt) in Eq. (A28), and since the expression must be averaged against time assuming a statistic distribution of the time origin ϕ/η such that cos(ηt + ϕ) (ϕ is different for each neutron), we get: where σ tot ≈ σ s for a regenerator with a low absorption cross section.