# Status of the GRANIT Facility

Authors: Damien Roulier; Francis Vezzu; Stefan Baeßler; Benoît Clément; Daniel Morton; Valery V. Nesvizhevsky; Guillaume Pignol; Dominique Rebreyend
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

The GRANIT facility is a follow-up project, which is motivated by the recent discovery of gravitational quantum states of ultracold neutrons. The goal of the project is to approach the ultimate accuracy in measuring parameters of such quantum states and also to apply this phenomenon and related experimental techniques to a broad range of applications in particle physics as well as in surface and nanoscience studies. We overview the current status of this facility, the recent test measurements, and the nearest prospects.


---

## Body

## 1. Introduction

The GRANIT facility [1, 2] is a follow-up project, which is motivated by the recent discovery of gravitational quantum states of ultracold neutrons (UCNs) [3–5]. The main goal of the project is to realize the resonance spectroscopy of those quantum states, with the prospect of achieving an unprecedented sensitivity. Such precision measurements would address in particular searches for extra short-range fundamental forces [6–9]. Also the phenomenon of gravitational quantum states and related experimental techniques could be applied to a broad range of other applications in particle physics as well as in surface and nanoscience studies [10].GRANIT is located at the level C of the Institut Laue-Langevin (ILL) in Grenoble, at the H172A beamline as shown in Figure1. The instrument comprises an ultracold neutron source based on the production of UCNs in superfluid helium and a spectrometer installed in an ISO 5 class clean room, pictured in Figure 2. A monochromatic neutron beam (0.89 nm wavelength) is extracted from a white cold neutron beam with a monochromator [11] and guided towards the superfluid helium bath where UCNs are produced [2, 12, 13]. An extraction guide allows us then to transfer those UCNs to the spectrometer.Figure 1
The GRANIT instrument at Level C of ILL, Grenoble.Figure 2
The spectrometer in the cleanroom. Measurement and control instruments are installed on a flat massive granite table with the mass of 550 kg, all inside a welded aluminum vacuum chamber with passivated wall surfaces. All this is placed in a clean controlled area (ISO 5) in order to protect the delicate optical elements.In the spectrometer, UCNs are first stored in an intermediate storage volume. To exit this volume, neutrons must go through an extraction slit of height ≈100μm, a compromise between the total UCN flux and the UCN phase-space density. Then neutrons will bounce over high quality mirrors very close to the surface and the method of resonance spectroscopy [14] will be applied. The resonance can be induced by a vibration of the bottom mirror, as used by the QBounce collaboration [15] or an oscillating magnetic field gradient as in the GRANIT spectrometer [16].In 2013 we have performed extensive tests of the various components of the facility and connected for the first time the source to the spectrometer. In this paper we will present the characterization of the whole UCN chain: the 0.89 nm neutron beam, the cryogenic production volume, the extraction guides, and the mirror assembly. We also present a confrontation of the measurements to Monte-Carlo simulations. Finally we present the first UCN flux measurement in the GRANIT spectrometer.


## 2. 0.89 nm Neutrons Beam

In superfluid helium, cold neutrons with the wavelength ofλ=0.89 nm can be converted into UCNs through resonant phonon excitation [17]. Therefore, the UCN production rate in the source will depend directly on the neutron flux at this precise wavelength (dΦ/dλ)|0.89nm.


### 2.1. Monochromator Adjustment

The monochromator [11] is composed of 18 intercalated stage-2 KC24 crystals, with a lattice constant of d=8.74 Å. According to the Bragg formula for the first order reflection, 2dsinθ=λ, the outgoing beamline angle corresponding to λ=0.89 nm is found to be 2θ=61.2 degrees, defining the geometry of the installation downstream.The position and orientation of the monochromator can be adjusted remotely with five parameters: rotation, two tilt angles, and two translation axes. These parameters are optimized by maximizing the UCN flux out of the source. The most critical parameter is the rotation of the monochromator, for which the neutron count rate varies by 80% of maximum a few degrees away from the optimal position. For the other parameters, within their whole range, the count rate varies at most by 40%. The result of the scan is shown in Figure3. We checked that maximum UCN flux coincides with the maximum cold neutron flux, indicating that the setup is aligned correctly.Figure 3
Neutron rate as a function of the monochromator rotation angle, for UCNs (dots) and cold neutrons (squares).


### 2.2. Differential Neutron Flux

We have characterized the wavelength distribution in the beam using the time of flight (TOF) technique. Two measurements were performed: the first over a flight length of78±1 cm, the second over 28±1 cm. The spectrum of the first measurement is presented in Figure 4, together with a fit of the peaks. We obtained for the first order peak the central wavelength λ = 0.879(11) nm. The uncertainty is dominated by the error on the flight length, which is itself defined by the uncertainty of knowledge of the position in the gaseous detector where the reaction occured. The width of the peak, σ=0.022 nm, is compatible with the expected TOF resolution.Figure 4
Long-distance (78 cm) TOF spectrum. It is assumed that the three main peaks arise from the three first orders of Bragg diffraction which fixes the wavelengthsλ/2 and λ/3 for the second and third order peaks. There are 21 other free parameters to account for the offset of the X-axis, the amplitudes and width of the Gaussian peaks, as well as ad hoc background description.The relative intensities of the peak shown in Figure4 do not reflect the relative intensities of the corresponding neutron fluxes due to at least three important effects. First, the detector efficiency is proportional to λ according to the 1/v law. Second, losses in air, estimated to be 13% for 0.89 nm neutrons, are also proportional to λ. Finally, since different wavelength components of the beam have different angular divergencies, the geometrical acceptance of the setup is also wavelength dependent. To suppress the latter, a second TOF measurement with a shorter flight length of 28±1 cm was performed.We show in Figure5 the TOF spectrum recorded at the short distance. After correcting for the detector efficiency, we determined the relative contribution of the first and second order diffraction peaks to the total flux as(1)p0.89=0.74,p0.445=0.26.Figure 5
Short-distance (28 cm) TOF spectrum. The analysis is similar to that for Figure4. The dashed line is here to guide the eye, and the corresponding data points were not included in the fit.


### 2.3. Integral Neutron Flux

To determine the integral neutron flux of the H172A beamline, we have used the standard gold foil activation technique that allows us to measure the capture flux density defined by(2)Φgold=∫‍dϕdλλλthdλ,where λth=0.18 nm is the wavelength of thermal neutrons. In July 2012 with a reactor power of 48 MW, we have obtained the following capture flux density:(3)Φgold=6.5×108cm-2s-1.Knowing the proportion of each diffracted peakpλ, neutron flux in the 0.89 peak can be calculated:(4)Φ[0.89nm]=Φgold×0.18nm0.89nm×∑λ∈peaks‍pλp0.89×λ0.89nm-1.From the differential neutron flux results (equation (1)), we estimate Φ[0.89nm]=(1.1±0.1)×108cm-2s-1. Although the intrinsic width of the peak is unknown, it cannot be larger than the measured TOF width σ=0.02 nm. Assuming a Gaussian shape for dΦ/dλ, we deduce a worst case scenario estimate (with our limit on σ) of the differential neutron flux at 0.89 nm:(5)dΦdλ0.89nm=Φ[0.89nm]×1σ2π.The actual flux inside the UCN conversion volume is further reduced by two effects. First, the walls of the conversion volume (0.85 mm thick aluminum and 1 mm thick beryllium) attenuates the beam by a factor of 0.895. Then, because of the angular divergence of the beam of 30 mrad, the fraction of the cold beam interacting with the source is 0.72 (estimated with a Monte-Carlo simulation).Thus, the effective 0.89 nm differential flux inside the source is(6)dΦdλ0.89nmeff=(1.3±0.2)×109cm-2s-1nm-1with the reactor power of 48 MW. As the conversion rate in BeO vessel can be estimated [2], the volumic production rate is(7)P=(4.97±0.38)×10-9nm·cm-1×dΦdλ0.89nmeff.We deduce that we produce P×V=32000 UCN/s in our 5-liter vessel with the reactor power of 48 MW.


## 2.1. Monochromator Adjustment

The monochromator [11] is composed of 18 intercalated stage-2 KC24 crystals, with a lattice constant of d=8.74 Å. According to the Bragg formula for the first order reflection, 2dsinθ=λ, the outgoing beamline angle corresponding to λ=0.89 nm is found to be 2θ=61.2 degrees, defining the geometry of the installation downstream.The position and orientation of the monochromator can be adjusted remotely with five parameters: rotation, two tilt angles, and two translation axes. These parameters are optimized by maximizing the UCN flux out of the source. The most critical parameter is the rotation of the monochromator, for which the neutron count rate varies by 80% of maximum a few degrees away from the optimal position. For the other parameters, within their whole range, the count rate varies at most by 40%. The result of the scan is shown in Figure3. We checked that maximum UCN flux coincides with the maximum cold neutron flux, indicating that the setup is aligned correctly.Figure 3
Neutron rate as a function of the monochromator rotation angle, for UCNs (dots) and cold neutrons (squares).


## 2.2. Differential Neutron Flux

We have characterized the wavelength distribution in the beam using the time of flight (TOF) technique. Two measurements were performed: the first over a flight length of78±1 cm, the second over 28±1 cm. The spectrum of the first measurement is presented in Figure 4, together with a fit of the peaks. We obtained for the first order peak the central wavelength λ = 0.879(11) nm. The uncertainty is dominated by the error on the flight length, which is itself defined by the uncertainty of knowledge of the position in the gaseous detector where the reaction occured. The width of the peak, σ=0.022 nm, is compatible with the expected TOF resolution.Figure 4
Long-distance (78 cm) TOF spectrum. It is assumed that the three main peaks arise from the three first orders of Bragg diffraction which fixes the wavelengthsλ/2 and λ/3 for the second and third order peaks. There are 21 other free parameters to account for the offset of the X-axis, the amplitudes and width of the Gaussian peaks, as well as ad hoc background description.The relative intensities of the peak shown in Figure4 do not reflect the relative intensities of the corresponding neutron fluxes due to at least three important effects. First, the detector efficiency is proportional to λ according to the 1/v law. Second, losses in air, estimated to be 13% for 0.89 nm neutrons, are also proportional to λ. Finally, since different wavelength components of the beam have different angular divergencies, the geometrical acceptance of the setup is also wavelength dependent. To suppress the latter, a second TOF measurement with a shorter flight length of 28±1 cm was performed.We show in Figure5 the TOF spectrum recorded at the short distance. After correcting for the detector efficiency, we determined the relative contribution of the first and second order diffraction peaks to the total flux as(1)p0.89=0.74,p0.445=0.26.Figure 5
Short-distance (28 cm) TOF spectrum. The analysis is similar to that for Figure4. The dashed line is here to guide the eye, and the corresponding data points were not included in the fit.


## 2.3. Integral Neutron Flux

To determine the integral neutron flux of the H172A beamline, we have used the standard gold foil activation technique that allows us to measure the capture flux density defined by(2)Φgold=∫‍dϕdλλλthdλ,where λth=0.18 nm is the wavelength of thermal neutrons. In July 2012 with a reactor power of 48 MW, we have obtained the following capture flux density:(3)Φgold=6.5×108cm-2s-1.Knowing the proportion of each diffracted peakpλ, neutron flux in the 0.89 peak can be calculated:(4)Φ[0.89nm]=Φgold×0.18nm0.89nm×∑λ∈peaks‍pλp0.89×λ0.89nm-1.From the differential neutron flux results (equation (1)), we estimate Φ[0.89nm]=(1.1±0.1)×108cm-2s-1. Although the intrinsic width of the peak is unknown, it cannot be larger than the measured TOF width σ=0.02 nm. Assuming a Gaussian shape for dΦ/dλ, we deduce a worst case scenario estimate (with our limit on σ) of the differential neutron flux at 0.89 nm:(5)dΦdλ0.89nm=Φ[0.89nm]×1σ2π.The actual flux inside the UCN conversion volume is further reduced by two effects. First, the walls of the conversion volume (0.85 mm thick aluminum and 1 mm thick beryllium) attenuates the beam by a factor of 0.895. Then, because of the angular divergence of the beam of 30 mrad, the fraction of the cold beam interacting with the source is 0.72 (estimated with a Monte-Carlo simulation).Thus, the effective 0.89 nm differential flux inside the source is(6)dΦdλ0.89nmeff=(1.3±0.2)×109cm-2s-1nm-1with the reactor power of 48 MW. As the conversion rate in BeO vessel can be estimated [2], the volumic production rate is(7)P=(4.97±0.38)×10-9nm·cm-1×dΦdλ0.89nmeff.We deduce that we produce P×V=32000 UCN/s in our 5-liter vessel with the reactor power of 48 MW.


## 3. UCN Source

The GRANIT superthermal UCN source is an evolution of the SUN1 apparatus [12, 13]. The UCN conversion volume consists of a vessel made out of BeO/Be filled in with superfluid 4He, where 0.89 nm neutrons are downscattered to the UCN energy range by resonant phonon excitation. This rectangular volume of the size 7×7×100 cm3 is placed in continuity with the neutron guide, which is also of squared section (7×7 cm2). The conversion volume is encased in a cryostat that allows cooling the volume down to below 0.8 K. A first commissioning of SUN1 at its definitive position was reported in [12, 13], where a large diameter extraction guide from the source to a UCN detector was installed. In order to avoid diluting of UCNs in the phase-space density, a more elaborate extraction guide assembly, with smaller diameter in particular, was designed and built. In the following we report the necessary modifications to the source as well as the current performances.To determine these performances, a temporary set-up was installed. A box (filled in with Argon to reduce the UCN losses) was connected to the UCN extraction guide outside the clean room. That set-up allowed us to perform our first measurements with UCNs.


### 3.1. Cryostat

With the assembly of the new extraction presented in Figure6, which connects the coldest part to the ambient temperature, first tests showed a loss of cryogenic power. In the first configuration it was not possible to cool down the UCN volume below the temperature of 1 K; thus we had to increase the cryogenic power of the refrigerator. A first temporary set-up using a LN2 cooling system for the thermal screens (around the inner parts of the extraction) was installed in order to validate that increasing power could solve the problem. This set-up permitted us to reach the temperature of 0.74 K. In the second configuration this LN2 cooling system was replaced with a Sumitomo cold head 150 W @ 77 K. The outlet box of the cryostat and some thermal screens were also replaced. The current cryogenic system is now adapted for our configuration. Figure 7 shows a cool down and filling of the conversion volume.Figure 6
Extraction guides from the source to the spectrometer. The extraction guides are composed of several tubular elements, which are thin foils of stainless steel inserted inside tubes. The design of the guides allows compensating for the misalignment between the source and the spectrometer.Figure 7
A successful cooldown in 100 hours. The temperature of several areas of the cryostat is constantly monitored. It takes 48 hours from the beginning to start4He liquefaction, which helps cooling the conversion volume. 48 hours later, the conversion volume is cold enough for keeping liquid 4He. 52 additional hours are needed to fill in this volume. Then, He-II can be cooled down to the temperature of 0.74 K.


### 3.2. Separation Windows

In several occurrences, we need to physically separate volumes at different pressures keeping UCN-transparency. This is the case for any gaseous UCN detector. Also, the spectrometer vacuum has to be separated from the extraction vacuum to avoid reheating in case of spectrometer openings. Any material separation is a cause of extra UCN losses due to quantum reflection and absorption.The first choice for a material was aluminum that has both low optical potential (54 neV, corresponding to a critical velocity of 3.2 m/s) and small neutronic absorption cross-section (~102 barn for 5 m/s neutrons). However, soft UCNs (with the velocity lower than 3.2 m/s) are supposed to represent a nonnegligible fraction of UCNs extractable out of the source. A simple calculation of the transmission through two windows (one in the extraction and one for the gaseous detector) considering an isotropic angular distribution of UCNs, illustrated in Figure8, shows that better materials exist for our application.Figure 8
An analytical calculation of the probability of transmission of UCNs through two foils (15μm each for Ti, V, and TiAl alloy, 30 μm each for Al) for an isotropic UCN gas as a function of UCN energy.We conclude that titanium would be better-suited as long as soft UCNs are available, and pure titanium windows are quite easy to set up. However, the thickness of the windows is far more critical for titanium than for aluminum (because of the absorption cross-section for neutrons). Thus, depending on mechanical constraints, the windows should be as thin as possible. Figure9 shows the extraction window between the extraction and the spectrometer that was designed for GRANIT.Figure 9
Vacuum-separating window in the extraction guides. The foil must guarantee the vacuum tightness but also it should resist to an overpressure of 1.5 bar max in both directions in case of eventual fast heating of the source. Such overpressure is limited using two safety valves. The window is made with 2 aluminum half pieces milled by electrocutting reinforced by letting ribs 0.25 mm thick. Between these two parts a titanium foil 0.015 mm thick and 2 o’rings are inserted. The assembly is bolted. The window assembly had been tested with success under vacuum and under 5 bar pressure.


### 3.3. UCN Count Rate versus Temperature

A measurement of the UCN count rate as a function of He-II temperature was realized with a3He counter. The extraction window was made of a 30 μm aluminum foil, as well as the detector entrance window. The result is presented in Figure 10.Figure 10
UCN count rate versus the temperature of He-II. The cold neutron beam constantly passes through the source and the UCN valve is opened periodically.Two regimes appear. WhenT≥1.2 K, the UCN flux is highly dependent on temperature because of the upscattering of UCNs in He-II. At lower temperature (T≤0.9 K), the UCN flux saturates, because the main losses are caused by absorption of UCNs in the beryllium walls of the conversion volume. If this effect was an order of magnitude smaller, the transition would occur at a lower temperature and the saturation UCN density would be nearly an order of magnitude higher.When the UCN valve of the source is open, radiative heat on He-II causes an important increase in temperature (0.045 K/min). For this reason, the valve should not be opened for longer than 10–15 seconds in order to have a reliable and stable measurement, as well as to be able to cool back down He-II (~10 min). Thus, the source can operate in a pulsed regime at a temperature below 0.9 K (stable as long as the opening time of the valve is short), or in a continuous mode at a higher but always stable temperature (~1.3 K), but with a smaller UCN flux.


### 3.4. Source and Extraction Characteristic Times

We define two characteristic times for this system. The emptying timeτemptying is the average time necessary to extract UCNs from the apparatus. The storage time τstorage is the average lifetime of UCNs in the isolated source.The emptying time is measured by accumulating UCNs in the source during 2 min, then closing the cold beam shutter and waiting for a few seconds. The UCNs are then released through the extraction to a UCN detector, with a differential UCN count measurement. The data are fitted with a single exponential as illustrated in Figure11. The results for different wainting times are summarized in Table 1. One can notice that the longer the waiting time, the longer the emptying time. We conclude that the softer UCNs are stored longer in the source, and that a sufficiently long time must be chosen to integrate the number of neutrons extracted without introducing a bias between measurements with different waiting times.Table 1
Emptying time of UCNs in the source and extraction with two titanium windows for several waiting times.
Waiting time (s)
τ
emptying (s)
0
12.9
±
0.4
50
18.6
±
0.6
100
23.8
±
1.4Figure 11
Emptying time of UCNs in the source and extraction with Ti windows.The storage time is measured in a similar way, but the total number of extracted neutrons is counted, and the procedure is repeated for different waiting times. This time is obtained by fitting the data with a double exponential, as shown in Figure12, and calculating the weighted geometric mean of the two decay constants. The two exponentials allow accounting for several UCN populations with different velocities and storage times where a single exponential is not sufficient.Figure 12
Storage time of UCNs in the source with Ti windows.The storage time measurement was performed for several combinations of extraction and detector windows materials. The results are summarized in Table2. Having at least one aluminum window cuts out the soft UCNs, resulting in similar storage times for Al-Al and Ti-Al configurations. Using only titanium allows recovering the soft UCNs. The noticable increase of the storage time is expected as slower UCNs have less collisions on the source walls. The choice of titanium windows is therefore justified and gives us access to a nonnegligible fraction of the UCN velocity spectrum.Table 2
Storage time in different configurations.
Extraction
Detector
τ
storage (s) @ 0.81 K
Al
Al
21.0
±
0.4
Ti
Al
21.3
±
0.4
Ti
Ti
30.4
±
0.7The results show thatτemptying<τstorage: the source can work in an accumulation mode where UCNs are accumulated in the source and then released all at the same time towards the spectrometer.


### 3.5. UCN Velocity Spectrum

The velocity distribution of UCNs can be determined with a free fall experiment, as described in Figure13. The height of free fall for a UCN with a defined horizontal velocity is(8)h=g2dfreefallvUCN2.Figure 13
Set-up for free fall measurement: in an argon-filled box, the UCN beam is collimated with a 12 cm long absorber-mirror system forming a slit of height 500μm. A detector with a 15 μm thick titanium window is placed at the distance of 26 cm from the slit, and record the neutron count as a function of the fixed free fall height.Assuming we have spatially isotropic distribution of UCNs within our collimation system (this assumption is valid with high accuracy for any broad angular distribution of UCNs in the extraction system), the measurement of horizontal velocity of the collimated UCNs is representative of the velocity distribution of the UCNs in the extraction.The recorded UCNs counts as a function of the free fall height are shown in Figure14. In order to obtain a velocity spectrum for this measurement, we make a Monte-Carlo simulation of the experiment to fit the initial velocity spectra to the data. Two shapes were assumed for the initial spectrum: a Gaussian distribution and an asymetric triangular distribution. Both fitted correctly the data, and the obtained results for the mean and RMS of the distribution are(9)μGauss=5.1m·s-1,σGauss=1.6m·s-1,μTriangle=5.3m·s-1,σTriangle=1.4m·s-1.The asymetric triangular distribution is zero below 2.2 m/s, maximum at 4.5 m/s, and zero above 9.0 m/s. In both cases, we obtain a coarse UCN velocity spectrum and notice a quite wide distribution around the mean value.Figure 14
Number of UCNs detected for 30 min versus the free fall height. During the measurement, the temperature of He-II was 1.35 K. Since the UCN spectrum is defined by the storage and transmission properties of the whole system, and the temperature of He-II affects only the total count rates, the exact value of the temperature is of no importance. This measurement was fitted with a Monte-Carlo simulation, assuming an initial triangular or Gaussian spectrum.Because of the low statistics (some hundreds of UCNs counted for 30 minutes), the valve was always open (thus the temperature was 1.35 K). In this configuration, there is no accumulation of soft UCNs in the source, thus a higher mean velocity than what we could have. This method of measurement for the velocity spectrum of UCNs, though lacking in precision, was successful. The precision could be improved by designing a UCN detector with a dedicated geometry or a large position-sensitive detector and a vacuum-tight environment. The result itself is in agreement with what we expected and confirms that the source/extraction system we use is well-suited for the GRANIT spectrometer.


## 3.1. Cryostat

With the assembly of the new extraction presented in Figure6, which connects the coldest part to the ambient temperature, first tests showed a loss of cryogenic power. In the first configuration it was not possible to cool down the UCN volume below the temperature of 1 K; thus we had to increase the cryogenic power of the refrigerator. A first temporary set-up using a LN2 cooling system for the thermal screens (around the inner parts of the extraction) was installed in order to validate that increasing power could solve the problem. This set-up permitted us to reach the temperature of 0.74 K. In the second configuration this LN2 cooling system was replaced with a Sumitomo cold head 150 W @ 77 K. The outlet box of the cryostat and some thermal screens were also replaced. The current cryogenic system is now adapted for our configuration. Figure 7 shows a cool down and filling of the conversion volume.Figure 6
Extraction guides from the source to the spectrometer. The extraction guides are composed of several tubular elements, which are thin foils of stainless steel inserted inside tubes. The design of the guides allows compensating for the misalignment between the source and the spectrometer.Figure 7
A successful cooldown in 100 hours. The temperature of several areas of the cryostat is constantly monitored. It takes 48 hours from the beginning to start4He liquefaction, which helps cooling the conversion volume. 48 hours later, the conversion volume is cold enough for keeping liquid 4He. 52 additional hours are needed to fill in this volume. Then, He-II can be cooled down to the temperature of 0.74 K.


## 3.2. Separation Windows

In several occurrences, we need to physically separate volumes at different pressures keeping UCN-transparency. This is the case for any gaseous UCN detector. Also, the spectrometer vacuum has to be separated from the extraction vacuum to avoid reheating in case of spectrometer openings. Any material separation is a cause of extra UCN losses due to quantum reflection and absorption.The first choice for a material was aluminum that has both low optical potential (54 neV, corresponding to a critical velocity of 3.2 m/s) and small neutronic absorption cross-section (~102 barn for 5 m/s neutrons). However, soft UCNs (with the velocity lower than 3.2 m/s) are supposed to represent a nonnegligible fraction of UCNs extractable out of the source. A simple calculation of the transmission through two windows (one in the extraction and one for the gaseous detector) considering an isotropic angular distribution of UCNs, illustrated in Figure8, shows that better materials exist for our application.Figure 8
An analytical calculation of the probability of transmission of UCNs through two foils (15μm each for Ti, V, and TiAl alloy, 30 μm each for Al) for an isotropic UCN gas as a function of UCN energy.We conclude that titanium would be better-suited as long as soft UCNs are available, and pure titanium windows are quite easy to set up. However, the thickness of the windows is far more critical for titanium than for aluminum (because of the absorption cross-section for neutrons). Thus, depending on mechanical constraints, the windows should be as thin as possible. Figure9 shows the extraction window between the extraction and the spectrometer that was designed for GRANIT.Figure 9
Vacuum-separating window in the extraction guides. The foil must guarantee the vacuum tightness but also it should resist to an overpressure of 1.5 bar max in both directions in case of eventual fast heating of the source. Such overpressure is limited using two safety valves. The window is made with 2 aluminum half pieces milled by electrocutting reinforced by letting ribs 0.25 mm thick. Between these two parts a titanium foil 0.015 mm thick and 2 o’rings are inserted. The assembly is bolted. The window assembly had been tested with success under vacuum and under 5 bar pressure.


## 3.3. UCN Count Rate versus Temperature

A measurement of the UCN count rate as a function of He-II temperature was realized with a3He counter. The extraction window was made of a 30 μm aluminum foil, as well as the detector entrance window. The result is presented in Figure 10.Figure 10
UCN count rate versus the temperature of He-II. The cold neutron beam constantly passes through the source and the UCN valve is opened periodically.Two regimes appear. WhenT≥1.2 K, the UCN flux is highly dependent on temperature because of the upscattering of UCNs in He-II. At lower temperature (T≤0.9 K), the UCN flux saturates, because the main losses are caused by absorption of UCNs in the beryllium walls of the conversion volume. If this effect was an order of magnitude smaller, the transition would occur at a lower temperature and the saturation UCN density would be nearly an order of magnitude higher.When the UCN valve of the source is open, radiative heat on He-II causes an important increase in temperature (0.045 K/min). For this reason, the valve should not be opened for longer than 10–15 seconds in order to have a reliable and stable measurement, as well as to be able to cool back down He-II (~10 min). Thus, the source can operate in a pulsed regime at a temperature below 0.9 K (stable as long as the opening time of the valve is short), or in a continuous mode at a higher but always stable temperature (~1.3 K), but with a smaller UCN flux.


## 3.4. Source and Extraction Characteristic Times

We define two characteristic times for this system. The emptying timeτemptying is the average time necessary to extract UCNs from the apparatus. The storage time τstorage is the average lifetime of UCNs in the isolated source.The emptying time is measured by accumulating UCNs in the source during 2 min, then closing the cold beam shutter and waiting for a few seconds. The UCNs are then released through the extraction to a UCN detector, with a differential UCN count measurement. The data are fitted with a single exponential as illustrated in Figure11. The results for different wainting times are summarized in Table 1. One can notice that the longer the waiting time, the longer the emptying time. We conclude that the softer UCNs are stored longer in the source, and that a sufficiently long time must be chosen to integrate the number of neutrons extracted without introducing a bias between measurements with different waiting times.Table 1
Emptying time of UCNs in the source and extraction with two titanium windows for several waiting times.
Waiting time (s)
τ
emptying (s)
0
12.9
±
0.4
50
18.6
±
0.6
100
23.8
±
1.4Figure 11
Emptying time of UCNs in the source and extraction with Ti windows.The storage time is measured in a similar way, but the total number of extracted neutrons is counted, and the procedure is repeated for different waiting times. This time is obtained by fitting the data with a double exponential, as shown in Figure12, and calculating the weighted geometric mean of the two decay constants. The two exponentials allow accounting for several UCN populations with different velocities and storage times where a single exponential is not sufficient.Figure 12
Storage time of UCNs in the source with Ti windows.The storage time measurement was performed for several combinations of extraction and detector windows materials. The results are summarized in Table2. Having at least one aluminum window cuts out the soft UCNs, resulting in similar storage times for Al-Al and Ti-Al configurations. Using only titanium allows recovering the soft UCNs. The noticable increase of the storage time is expected as slower UCNs have less collisions on the source walls. The choice of titanium windows is therefore justified and gives us access to a nonnegligible fraction of the UCN velocity spectrum.Table 2
Storage time in different configurations.
Extraction
Detector
τ
storage (s) @ 0.81 K
Al
Al
21.0
±
0.4
Ti
Al
21.3
±
0.4
Ti
Ti
30.4
±
0.7The results show thatτemptying<τstorage: the source can work in an accumulation mode where UCNs are accumulated in the source and then released all at the same time towards the spectrometer.


## 3.5. UCN Velocity Spectrum

The velocity distribution of UCNs can be determined with a free fall experiment, as described in Figure13. The height of free fall for a UCN with a defined horizontal velocity is(8)h=g2dfreefallvUCN2.Figure 13
Set-up for free fall measurement: in an argon-filled box, the UCN beam is collimated with a 12 cm long absorber-mirror system forming a slit of height 500μm. A detector with a 15 μm thick titanium window is placed at the distance of 26 cm from the slit, and record the neutron count as a function of the fixed free fall height.Assuming we have spatially isotropic distribution of UCNs within our collimation system (this assumption is valid with high accuracy for any broad angular distribution of UCNs in the extraction system), the measurement of horizontal velocity of the collimated UCNs is representative of the velocity distribution of the UCNs in the extraction.The recorded UCNs counts as a function of the free fall height are shown in Figure14. In order to obtain a velocity spectrum for this measurement, we make a Monte-Carlo simulation of the experiment to fit the initial velocity spectra to the data. Two shapes were assumed for the initial spectrum: a Gaussian distribution and an asymetric triangular distribution. Both fitted correctly the data, and the obtained results for the mean and RMS of the distribution are(9)μGauss=5.1m·s-1,σGauss=1.6m·s-1,μTriangle=5.3m·s-1,σTriangle=1.4m·s-1.The asymetric triangular distribution is zero below 2.2 m/s, maximum at 4.5 m/s, and zero above 9.0 m/s. In both cases, we obtain a coarse UCN velocity spectrum and notice a quite wide distribution around the mean value.Figure 14
Number of UCNs detected for 30 min versus the free fall height. During the measurement, the temperature of He-II was 1.35 K. Since the UCN spectrum is defined by the storage and transmission properties of the whole system, and the temperature of He-II affects only the total count rates, the exact value of the temperature is of no importance. This measurement was fitted with a Monte-Carlo simulation, assuming an initial triangular or Gaussian spectrum.Because of the low statistics (some hundreds of UCNs counted for 30 minutes), the valve was always open (thus the temperature was 1.35 K). In this configuration, there is no accumulation of soft UCNs in the source, thus a higher mean velocity than what we could have. This method of measurement for the velocity spectrum of UCNs, though lacking in precision, was successful. The precision could be improved by designing a UCN detector with a dedicated geometry or a large position-sensitive detector and a vacuum-tight environment. The result itself is in agreement with what we expected and confirms that the source/extraction system we use is well-suited for the GRANIT spectrometer.


## 4. Spectrometer

### 4.1. Description

UCNs are transported through the extraction guide from the He-II source to the spectrometer as shown in Figure6. The extraction guides and the extraction window are the same as the previous configuration. The guides connect the source to an intermediate storage volume depicted in Figure 15. Optical elements, at the heart of the spectrometer, are connected at the exit of the storage volume. They are sets of silica pieces with different coatings and roughness states, providing different conditions to the transport of UCNs. They are described in Figures 16 and 17. All these parts have been provided by the SESO Company in Aix en Provence (France), and Diamond-Like Carbon (DLC) coatings were applied by Advanced Material Laboratory (IN2P3, France).Figure 15
The DLC-coated storage volume. The inner sizes are40×40×340 mm3.Figure 16
Schematics and pictures of the optical elements. The extraction mirror and the scatterer have the same sizes (300×100×20 mm3) and each consists of 3 pieces of 100×100×20 mm3 glued together. The mean roughness of the extraction mirror (floor) is very low (0.5 nm) to allow for specular bouncing, whereas the mean roughness of the scatterer (ceiling) is intentionally high (5.6 μm) to allow UCNs with too high vertical velocity to be diffused back. The transport mirror (300×250×70 mm3) has a surface mean roughness of 0.5 nm and a planarity of 80 nm. The absorber mirror is 300×90×30 mm3.Figure 17
The optical elements on the granit table. The extraction mirrors assembly and the transport mirror are placed on two separate adjustable supports. Their adjustment could be done with3+3 micrometric screws. To adjust the height and the orientation of the surface of the transport mirror with a great accuracy, we use 3 piezo-electric elements. The distance between the absorber and the transport mirror is adjustable as well using 3 piezo-electric elements. The piezos are driven from the control computer with a Labview application.The intermediate storage volume, made of several aluminum parts coated with DLC, is designed to randomize the UCNs trajectories. In order to close the storage volume, a nickel-coated butterfly valve is used, preventing the UCNs from returning into the source. The exit of the intermediate storage volume is closed by the first optical elements of the spectrometer, forming a slit and allowing only UCNs with negligible vertical velocity to be transmitted. The slit is composed of an extraction mirror and a scatterer. The scatterer is placed above the extraction mirror at a height between 50 and 200μm. Both are coated with DLC in order to provide reflection of UCNs with a broad velocity range from the surfaces and thus to assure proper operation of the so-called semidiffusive slit [18, 19].The second part of the spectrometer is a transport mirror, aligned with the slit (or positionned a few micrometers above/below the slit, allowing the selection of quantum states). The magnetic excitation array [16], described in Figure 18, can be placed above it. To induce resonant transitions between quantum states, a periodic magnetic field gradient will be produced with an array of wires located above the transport mirror. The third part of the spectrometer is an absorber, placed above and at the end of the transport mirror at an adjustable height to filter the quantum states.Figure 18
The magnetic excitation array is made of 4 modules of 32 wires each (1 mm2) constituting a 128-wire array. The wires are spaced by 0.25 mm.


### 4.2. Cooling of the Magnetic Excitation Array

The magnetic excitation array was built and tested before set-up. The current needed (5 A in average) to generate a sufficient magnetic field gradient produces 50 W of power that has to be dissipated without affecting too much the transport mirror. The wire system must be cooled down. Circulation of cooled gaseous N2 will be used for that purpose. To test the cooling in conditions as close as in the experimental set-up, the magnetic excitation array was placed in a vacuum chamber and connected to a power supply (5 A, 8 V). The heat exchanger was placed in a LN2 tank at a distance of 9 m from the magnetic excitation array, and a circulating N2 circuit linked the heat exchanger to the magnetic excitation array. Results of that test were very good, and we could easily stabilize the temperature of the wires at 15°C during 12 h without human intervention. The magnetic excitation array equipped with the cooling system will be installed during autumn 2014.


### 4.3. Preliminary Results

In summer 2013, we conducted the first tests of the full extraction chain. During this cycle, the source temperature could not be cooled lower than 1.35 K. Removing the transport mirror and the absorber, and placing the detector at the exit of the extraction slit of height 127μm, we measured the first UCNs in the GRANIT spectrometer. The 3He detector was equipped with a titanium window. In this configuration, we measured (10.9±1.5)×10-3 UCN/s as shown in Figure 19.Figure 19
Amplitude spectrum of UCNs and background. The temperature of the source was 1.35 K.This count rate is one order of magnitude lower than expected in this configuration and two orders of magnitude lower than targeted. Potential explanations include (i) a defect in the extraction which, because of the modifications of the configuration, did not exist in the previous measurements, (ii) defects in the DLC-coated storage volume, which has not been characterized on its own due to its particular geometry, and (iii) the collection of impurities on a cold spot on the extraction window. This opens the way to further improvements and tests in the next cycles, which could not be done immediately because of the long reactor shutdown right after this measurement. In addition, the background, estimated to14×10-3 events/s, is high compared to the signal because of background fast neutrons in the reactor building. Extra-shielding will be added to the installation in order to reduce this background.


## 4.1. Description

UCNs are transported through the extraction guide from the He-II source to the spectrometer as shown in Figure6. The extraction guides and the extraction window are the same as the previous configuration. The guides connect the source to an intermediate storage volume depicted in Figure 15. Optical elements, at the heart of the spectrometer, are connected at the exit of the storage volume. They are sets of silica pieces with different coatings and roughness states, providing different conditions to the transport of UCNs. They are described in Figures 16 and 17. All these parts have been provided by the SESO Company in Aix en Provence (France), and Diamond-Like Carbon (DLC) coatings were applied by Advanced Material Laboratory (IN2P3, France).Figure 15
The DLC-coated storage volume. The inner sizes are40×40×340 mm3.Figure 16
Schematics and pictures of the optical elements. The extraction mirror and the scatterer have the same sizes (300×100×20 mm3) and each consists of 3 pieces of 100×100×20 mm3 glued together. The mean roughness of the extraction mirror (floor) is very low (0.5 nm) to allow for specular bouncing, whereas the mean roughness of the scatterer (ceiling) is intentionally high (5.6 μm) to allow UCNs with too high vertical velocity to be diffused back. The transport mirror (300×250×70 mm3) has a surface mean roughness of 0.5 nm and a planarity of 80 nm. The absorber mirror is 300×90×30 mm3.Figure 17
The optical elements on the granit table. The extraction mirrors assembly and the transport mirror are placed on two separate adjustable supports. Their adjustment could be done with3+3 micrometric screws. To adjust the height and the orientation of the surface of the transport mirror with a great accuracy, we use 3 piezo-electric elements. The distance between the absorber and the transport mirror is adjustable as well using 3 piezo-electric elements. The piezos are driven from the control computer with a Labview application.The intermediate storage volume, made of several aluminum parts coated with DLC, is designed to randomize the UCNs trajectories. In order to close the storage volume, a nickel-coated butterfly valve is used, preventing the UCNs from returning into the source. The exit of the intermediate storage volume is closed by the first optical elements of the spectrometer, forming a slit and allowing only UCNs with negligible vertical velocity to be transmitted. The slit is composed of an extraction mirror and a scatterer. The scatterer is placed above the extraction mirror at a height between 50 and 200μm. Both are coated with DLC in order to provide reflection of UCNs with a broad velocity range from the surfaces and thus to assure proper operation of the so-called semidiffusive slit [18, 19].The second part of the spectrometer is a transport mirror, aligned with the slit (or positionned a few micrometers above/below the slit, allowing the selection of quantum states). The magnetic excitation array [16], described in Figure 18, can be placed above it. To induce resonant transitions between quantum states, a periodic magnetic field gradient will be produced with an array of wires located above the transport mirror. The third part of the spectrometer is an absorber, placed above and at the end of the transport mirror at an adjustable height to filter the quantum states.Figure 18
The magnetic excitation array is made of 4 modules of 32 wires each (1 mm2) constituting a 128-wire array. The wires are spaced by 0.25 mm.


## 4.2. Cooling of the Magnetic Excitation Array

The magnetic excitation array was built and tested before set-up. The current needed (5 A in average) to generate a sufficient magnetic field gradient produces 50 W of power that has to be dissipated without affecting too much the transport mirror. The wire system must be cooled down. Circulation of cooled gaseous N2 will be used for that purpose. To test the cooling in conditions as close as in the experimental set-up, the magnetic excitation array was placed in a vacuum chamber and connected to a power supply (5 A, 8 V). The heat exchanger was placed in a LN2 tank at a distance of 9 m from the magnetic excitation array, and a circulating N2 circuit linked the heat exchanger to the magnetic excitation array. Results of that test were very good, and we could easily stabilize the temperature of the wires at 15°C during 12 h without human intervention. The magnetic excitation array equipped with the cooling system will be installed during autumn 2014.


## 4.3. Preliminary Results

In summer 2013, we conducted the first tests of the full extraction chain. During this cycle, the source temperature could not be cooled lower than 1.35 K. Removing the transport mirror and the absorber, and placing the detector at the exit of the extraction slit of height 127μm, we measured the first UCNs in the GRANIT spectrometer. The 3He detector was equipped with a titanium window. In this configuration, we measured (10.9±1.5)×10-3 UCN/s as shown in Figure 19.Figure 19
Amplitude spectrum of UCNs and background. The temperature of the source was 1.35 K.This count rate is one order of magnitude lower than expected in this configuration and two orders of magnitude lower than targeted. Potential explanations include (i) a defect in the extraction which, because of the modifications of the configuration, did not exist in the previous measurements, (ii) defects in the DLC-coated storage volume, which has not been characterized on its own due to its particular geometry, and (iii) the collection of impurities on a cold spot on the extraction window. This opens the way to further improvements and tests in the next cycles, which could not be done immediately because of the long reactor shutdown right after this measurement. In addition, the background, estimated to14×10-3 events/s, is high compared to the signal because of background fast neutrons in the reactor building. Extra-shielding will be added to the installation in order to reduce this background.


## 5. Simulation

STARucn (Simulation of Transmission, Absorption and Reflection of ultracold neutrons) is a public Monte-Carlo software designed to simulate experimental set-ups and guides for UCNs, developped at LPSC Grenoble [20]. It relies heavily on CERN’s ROOT packages [21]. Its main features are modularity, easy configuration of geometry and simulation, propagation of UCNs, with or without gravity, interaction in volumes through effective lifetime, and interaction at surface (quantum reflection or transmission, surface losses due to absorbtion upon reflection, and specular and diffuse reflection). These features had been successfully benchmarked against analytical calculations and existing codes such as MCUCN of PSI and GEANT4UCN [22].In STARucn, a material is described with its Fermi potentialVF, its fraction of diffuse reflections d, its loss factor η defined as the ratio of imaginary to real part of the Fermi potential, and its effective lifetime τ in the material.Several measurements of the source characteristics can be used to estimate the unknown values of the parameters of the various materials. Our simulation (with geometries such as in Figure20) shows that the build-up time within the source is only driven by the losses on the Be and BeO surfaces of the source. Using results from 2010 measurement [12, 13], we set ηBe=2.7×10-4. In a similar way, the emptying time depends mostly on the reflective properties of the stainless steel extraction. Assuming this property is uniform in all the extraction, the best fit is found at dsteelextraction=3%. The parameter dBe was arbitrarily set to 0.1 because the simulation showed that it had no significant influence on any measurement.Figure 20
Geometry described and used in the simulation.Assuming a production of 38000 UCN/s isotropically distributed in phase space in the source (corresponding to the reactor power of 56 MW), the expected number of UCNs detected in the simulation is around 5 times higher than what is measured. No satisfactory explanation was found: for instance, a 12 cm2 hole in the extraction combined with a high ηsteelextraction would be needed. The most relevant figures are in Table 3.Table 3
Comparison of previous and recent measurements with simulations. The number of extractible UCNs correponds to the number of UCNs available in the source after the UCN density has been saturated according to the simulation. The number of extracted UCNs corresponds to the number of UCNs detected during the experiments. The measured (resp., simulated) extraction efficiency is the ratio of the number of UCNs detected during the measurement (resp., the simulation) to the simulated number of extractible UCNs. The missing factor is calculated according to the discrepancy between these two efficiencies.
Configuration
2010 [12, 13]
2013
Storage time (meas.)
67 s
30 s
Extractible UCNs (sim.)
2.2
×
1
0
6
1.0
×
1
0
6
Extracted UCNs (meas.)
274000
21000
Extraction efficiency (meas.)
13%
2.1%
Extraction efficiency (sim.)
65–50%
15–5%
Missing factor
4-5
2.5–7Possible candidates to explain these losses are a larger angular divergency of the incoming beam of cold neutrons in the source than estimated, an eventual misalignement of the monochromator reducing the production rate, eventual defects or losses in the extraction, and/or degradation of the BeO source volume. These options will be investigated during the next ILL reactor cycles.Moreover, simulating the experiment of Section4.3, with the corresponding count rate as measured in Section  2.1 and the spectrum as measured in Section  3.5, the number of neutrons measured is 10 times lower than simulated. However, no other measurement during this cycle can rule out a problem with the source or the extraction.


## 6. Improvements

These tests of all elements of the GRANIT experiment prompt us to modify some parts of the set-up. Some of improvements are quite easy to perform. The cleanliness of the extraction guides could be easily improved, starting with a rigorous and defined production process, then having good procedures for washing, stoving, and packaging each part just before installation.Moreover, due to the narrow UCN velocity range produced in the UCN source, DLC coating in extraction, mirrors, and intermediate storage volume is no longer crucial. For the next steps, we will replace the stainless steel foils of the extraction (that were dued to be DLC-coated) by electropolished tubes. This will decrease the UCN losses due to the properties of the surface (mainly roughness and hydrogen adsorption) of these foils and will considerably help for the assembly of the guide, as well as for connecting source to spectrometer, thus diminishing risks of damage.For the same reasons the intermediate storage volume will be replaced with a new one made of oxygen free high conductivity copper and in a cylindrical shape thus increasing the life time of UCNs in the volume and decreasing UCN losses. Another effort must be done to adjust the optical elements with more reliability without the contact of a probe. A solution using laser sensors is under study.


## 7. Conclusion

We have overviewed the current status of the GRANIT facility. First complete test of the GRANIT UCN source and spectrometer was performed during the last reactor cycle in 2013. Further improvements are identified based on the measured results, and they are being implemented.



---

*Source: 730437-2015-03-09.xml*