Quantized Ultracold Neutrons in Rough Waveguides: GRANIT Experiments and Beyond

We apply our general theory of transport in systems with random rough boundaries to gravitationally quantized ultracold neutrons in rough waveguides as in GRANIT experiments (ILL, Grenoble). We consider waveguides with roughness in both two and one dimensions(2Dand1D).Inthebiaseddiffusionapproximationthedepletiontimesforthegravitationalquantumstatescanbeeasilyexpressedviaeachotherirrespectiveofthesystemparameters.Thecalculationoftheexitneutroncountreducestoevaluationofasingleconstantwhichcontainsacomplicatedintegralofthecorrelationfunctionofsurfaceroughness.Inthecaseof1Droughness(randomgrating)thisconstantiscalculatedanalyticallyforcommontypesofthecorrelationfunctions.Theresultsobeysimplescalingrelationswhichareslightlydifferentin1Dand2D.WepredicttheexitneutroncountforthenewGRANITcell.


Introduction
One of the most interesting recent achievements in neutron physics is a series of GRANIT [1][2][3][4] and GRANIT-inspired [5][6][7] experiments aimed at detection of quantization of neutron motion by the Earth gravitational field.The experiment hinges on sending a beam of ultracold neutrons between rough and flat mirrors.The rough mirror scatters away and, eventually, absorbs the neutrons in higher gravitational quantum states leaving only the particles in the lowest states which reach the neutron counter.The observation of these quantized ultralow energy levels in the peV range opens a door to using neutrons for probing weak fundamental forces.Success of these experiments stimulates similar proposals for studying atomic and antimatter beams [8][9][10][11][12][13].
Another potential application area is the use of such experiment as a test for a quantum transport theory in systems with random rough boundaries (see, for example, [14][15][16][17][18][19][20][21][22][23]; a brief review of relevant theoretical methods can be found in the beginning of [24]).The neutrons are sensitive only to geometrical and statistical properties of random surface inhomogeneities and their behavior in the rough waveguide can serve as a perfect application of our transport theory including a model-free description of the GRANIT experiments which is the main goal of this paper.Some of our earlier results for a neutron waveguide with 1D roughness (a mirror with random grating) can be found in [25,26].

General Equations
Recently we developed a consistent perturbative approach to quantum transport along rough surfaces [24].Within our approach, the roughness-driven transition probabilities    (q − q  ) between the states (, q) → (  , q  ) decouple into a product of the Fourier image of the correlation function of surface roughness (q − q  ) and the boundary values of the wave functions in the absence of roughness (q is the particle momentum along the surface; ,   are the quantum numbers describing the spatial quantization for confined motion perpendicular to the walls).In essence, the correlation function of surface roughness plays the role similar to the impurity cross-section for transport in systems with bulk impurities.

Advances in High Energy Physics
In the case of neutron beams propagating between one rough mirror and one flat mirror, the transition probabilities have the form where Ψ  () is the value of the wave function in quantum state  on the mirror in the absence of roughness and   is the neutron absorption barrier for the mirror material.If the absorption barrier is very high,   → ∞ (1) becomes The wave functions Ψ  and the gravitational quantum states  q of neutrons between two horizontal mirrors are described in detail in [25,27].The transition probabilities    (q, q  ) determine the collision operator in the transport equation which in this case is a set of coupled Boltzmann-like equations for the distribution functions  q : The contribution from transitions to the states within continuous spectrum above the absorption threshold   , for which there are no reverse processes, is negligible [25].Integration with the energy -function in (3) reduces these equations in the relaxation time approximation  q (q) = ( −   )  to where  2  /2 =  −   and  is the angle between q  and q   .We use the dimensionless variables, which are common to the field (for details see [25]).All distances are measured in units of  0 , where  0 = ℏ 2/3 (2 2 ) −1/3 ∼ 5.871 m is the size of the lowest quantum state for neutrons in the infinite gravitational trap (open geometry without an upper mirror).The dimensionless distance between the mirrors ℎ = / 0 in experiment typically does not go down below 2. The average amplitude and the correlation radius of surface roughness  = ℓ/ 0 and  = / 0 are usually within the 0.1 ÷ 1 range.The energies are scaled by  0 =  0 ∼ 0.602 peV∼ 9.6366 × 10 −32 J which is the energy of a neutron in the lowest gravitational quantum state.The quantized levels in the gravity field   (ℎ) =   ()/ 0 start from about 2. In these units, the typical kinetic energy of particles in the beam  = / 0 and the absorption threshold   =   / 0 are very large, ,   ∼ 10 5 , and low gravitational states   ≪  are not sensitive to the behavior of the potential near the absorption threshold   .Velocities (momenta) along the wall V  in the beam direction  are measured in units of The characteristic times can be measured in units of  0 , which provides the scale for the oscillation frequency of neutrons in the gravitational well.In the original GRANIT cell [1][2][3]   / ∼ 0.16 and the time of flight of neutrons through the cell   is   / 0 ≃ 23.In the new cell [28]   ∼  and   / 0 ≃ 26.Diffusion of neutrons between discrete states  has a strong directional bias upward, towards higher states [25].The bias is explained by the rapid growth of the boundary values of the product of wave functions in (1) and ( 2) with increasing  and   (roughly, as  2  2 ).This increase in the rate of jumps  →   is checked by the decay of the correlation function (|q  − q   |) at large |q  − q   | which is determined by the value of the correlation radius , |q  − q   | ≲ 1/.In the end, the transition rates  →   represent, as a function of   , a relatively narrow peak around some  1 ≫ .
This strong upward bias has two consequences.First, almost all the time   necessary for a neutron, which is initially in a low gravitational state , to go up in states and disappear over the absorption barrier   , is spent on the first transition upwards.Further transitions to higher and higher states are going faster and faster.And second, this bias allows one to neglect the return of particles back to the lowest states.This results in decoupling of (4): Since the high narrow peak for transitions rates  →   is centered around some  1 ≫ , the absorption times   for neutrons that initially occupy low gravitational states  differ from each other, according to (1), only by the values of the wave functions on the rough mirror Ψ 2  (), where  1 is the depletion time for the first-the lowestgravitational state, and the coefficient 10 5 is inserted purely for the computational convenience.The ratios   (ℎ)/ 1 (ℎ) (7) are plotted in Figure 1 for  = 2; 3; 4; 5 as a function of the distance between the mirrors ℎ.Note that these ratios   / 1 in the biased diffusion approximation are the same for 1D and 2D roughness and do not depend on the roughness parameters at all.Though all the depletion times   rapidly decrease with decreasing ℎ (see below), the ratios   (ℎ)/ 1 (ℎ) are increasing.
In [29,30] the linewidths 1/  were considered as independent fitting parameters.Equation (7) shows that these   are not independent and are trivially related to each other.In the next two sections we calculate these   and express them via the parameters of the correlation function of surface roughness.
The fact that the depletion times   can be relatively close to each other illustrates the difficulty in observing the stepwise dependence of the exit neutron count   on ℎ, where   (0) is the number of neutrons in state  entering the waveguide of length .For the lowest gravitational states the velocities V  are more or less the same, V  ≈ V 0 √, and the exit neutron count is If all the states in front of the waveguide are equally populated,   (0) =  0 (9) becomes and all parameters of the experimental setup collapse into the single constant Φ (or  1 ).What makes Φ the most important parameter in the problem, technically more important than even  1 , is that it does not depend on the waveguide width ℎ while the relaxation time  1 is inversely proportional to  1 (ℎ).Constant Φ is obviously different for waveguides with 1D and 2D roughness, Φ 1 and Φ 2 , and these two situations should be considered separately.The value of the dimensionless constant Φ depends on the properties of the waveguide, but the moment we calculate Φ we know the dependence of the exit neutron count on the waveguide width ℎ, Figure 2, irrespective of the origin of these particular values of Φ.In Figure 2 taken from [26] we plot the exit neutron count   (ℎ)/ 0 (11) for several values of Φ.The noticeable quantum steps on the curves start appearing for Φ > 40.These steps correspond to distinct consecutive depletions of the gravitational quantum states; at lower Φ the depletion processes overlap.The experimental data from [1][2][3][4] are fitted to the curve with Φ = 23.5 which assumes 1D Gaussian roughness with  = 0.19 and  = 1.19 (see below).The fitting parameters are the number of neutrons entering the waveguide  0 and, to a lesser degree, the width of the waveguide which cannot be measured precisely.This is much smaller than the set of fitting parameters used in [29,30].The quality of the fit is reasonably good taking into account that in earlier experiments [1][2][3][4] the absolute value of the width of the rough waveguide ℎ could not be measured with an accuracy better than 10% and that the correlation function of surface roughness and the distribution of neutrons entering the waveguide   (0) have not been measured at all.
The only remaining task is to calculate  1 , and, by extension, the dimensionless parameter Φ (10).

Waveguides with 1D Roughness
We start from waveguides with 1D roughness for which most of the calculations can be carried out analytically.Towards the end of the paper we will mention why 1D roughness is important though the existing rough mirrors exhibit 2D roughness [31].
It is more convenient to start not from ( 12) but from (6). 1D roughness is a random grating perpendicular to the direction of the beam and the 2D momentum both before and after scattering has only one component  along the beam, q = (, 0).The transition probabilities    and the distribution function have the form where   is the number of particles in state  per unit length of the beam and   = √ −   , and  is the overall energy of the neutrons.Then (6) reduce to In our dimensionless variables the scattering probabilities (1) and ( 13) in ( 15) obtain the form where  1 stands for the dimensionless part of the Fourier image of 1D correlation function of surface roughness: the lower index 1 indicates that this equation describes the case of 1D roughness.Since the transition rate  →   rapidly increases with increasing   , only the terms with large   make a noticeable contribution to the sum (16) which can be replaced by the integration: In these notations our main parameter Φ 1 obtains a very simple form: and the problem reduces to the integrations ( 21) and ( 23) of the correlation functions of various functional forms.If, as it is often assumed, the correlation function is Gaussian,  1 (/) =  2 0  2 exp(− 2 /2 2 ), its Fourier image is also Gaussian,  1 () = exp(− 2 /2) and ) .
This equation has been used in [26] when fitting the experimental data from [1][2][3][4] (see Figure 2).The only information we have about the mirror roughness in those earlier experiments is that the average lateral size of inhomogeneities is of the order 1.19 and the average amplitude is approximately 0.119.Assuming that these numbers give the values of  and  and that we are dealing with the 1D Gaussian roughness, ( 25) yields Φ 1 ∼ 23.5.This value is used in the fit in Figure 1.Since these assumptions are somewhat arbitrary, the quality of the fit is actually much better than one might expect with all the uncertainties.The Fourier image of a power law correlation function, behaves like an exponential function: In the opposite case, when the power spectrum of roughness  1 () is given by a power law function, the correlation function itself behaves exponentially: and the constant Φ 1 has the form The purely exponential correlation functions in configuration or momentum spaces emerge from (30) or (27) when  = 1/2 or  = 1/2 : The parameter Φ 1 and the relaxation times   exhibit a universal dependence on the amplitude  and the correlation radius  of the surface roughness and on the overall kinetic energy of particles in the beam : This scaling is precise only at   → ∞.For finite   the powers of  and  are slightly different (e.g., at   ≈ 10 5 the scaling for the exponential correlator is  −0.493 ).These equations allow one to find the depletion times   for each gravitational state, (7) and (10).Figure 3 presents inverse depletion times 1/  (ℎ) in kHz for  = 1; 2; 3; 4; 5 in a waveguide with exponential roughness  =  = 1, (31) and (32).The values of the neutron velocity and time of flight entering  1 (24) correspond to the original GRANIT cell [1][2][3][4]:   / 0 ≈ 23,   / ≈ 0.16.The scaling (33) allows one to get 1/  for other values of , .The depletion becomes faster or slower depending on the value of  −2  1/2 .
The depletion times   determine both the disappearance of neutrons and broadening of the gravitational energy states   + ℏ/  (ℎ) with decreasing width of the waveguide ℎ.The broadening of the state  increases dramatically and rapidly reaches the separation between the states after the width ℎ becomes smaller than some critical value ℎ  .Below ℎ  the gravitational state   looses its discrete quantum nature; this happens simultaneously with depletion of the neutron population of this state.

Waveguides with 2D Roughness
In contrast to systems with 1D roughness, most of the calculations in 2D cases can be done only numerically.We start from (12) which in our notations acquires the form similar to (20): where After replacing the summation by the integration, we get equations similar to ( 23) and ( 24): The calculation of the zeroth angular harmonic of the 2D correlation function in momentum space  (0) 2 can be done analytically [32] for the Gaussian and exponential correlation functions.Even in these cases all further calculations should be done numerically.
In contrast to the earlier experiments, the roughness correlation function for the ongoing GRANIT experiments in a new cell has been measured [31].This correlation function is very close to the isotropic 2D exponential function with the correlation radius  = 0.65 and amplitude  = 1.02 while the time of flight   / 0 ≃ 26.This roughness yields Φ 2 = 5.22 × 10 3 .The numerical examples below assume this setup.Our results for other correlation functions can be found in [31] in which we discuss the identification of the correlation function for the new rough mirror.
As in the case of 1D roughness, both Φ and 1/  remain proportional to  2 .This scaling is obvious.The dependence on the correlation radius is more elusive and we cannot get an analytical expression similar to (33).The reason is the presence of  1/2 in the argument of  (0) 2 (40) in the integrand in (36).Our numerical data show that the scaling for  > 0.3 remains similar to (33): At   → ∞ the index,  = −1/2, the same as in the case of 1D roughness.The deviation from  = −1/2 at finite   is more pronounced than in 1D case: at   = 10 5 , as in experiment, the index  becomes −0.465 instead of −0.493 as in 1D.The scaling index for energy in the 2D case is obviously different from ( 33) because of different dependencies on velocities in the integrand.Instead of  = −3/4 as in (33), we find  ≃ −1.165 though, in contrast to , the quality of the fit is not very good.The critical values of ℎ  , below which the depletion and the broadening of the gravitational state   (ℎ) explode, are slightly larger than those for the 1D roughness.This is explained by the role of the sideway scattering which is absent in the 1D case.Otherwise, Figures 3 and 4 are similar.
Our prediction for the neutron count   (ℎ)/ 0 for the new GRANIT cell with the 2D exponential roughness with  = 1.03 and  = 0.65 (Φ 2 = 5.22 × 10 3 ) is given in Figure 5.The curve exhibits much more pronounced quantum steps than those in Figure 1.It is worth mentioning that the neutron count   is much less sensitive to the value of Φ when Φ is in (2 ÷ 8) × 10 3 range than when it is in the range 5 ÷ 40 as in Figure 1.
The presence of well-developed steps on the curve, which correspond to consecutive depletion of the gravitational quantum states, is explained mostly by a relatively large amplitude of roughness , several times bigger than for the older cell.This large value of the amplitude of roughness  presents challenges for both theory and experiment and can potentially degrade the accuracy of information extracted from the experimental data.On a theoretical side, the main assumption of the theory,  ≪ ,ℎ, is violated.For experiment, the large amplitude of roughness means that it becomes virtually impossible to accurately measure the distance between the mirrors ℎ.It also means that the energy levels become broad and not very well defined even if one disregards the interlevel transitions.The most likely consequence of all these factors would be some smearing of the steps in comparison to those in Figure 5.

Conclusions
In summary, we developed a quantitative theory of propagation of ultracold neutrons through a rough waveguide.The immediate applications are the ongoing GRANIT experiments at ILL (Grenoble) aimed at analysis of quantization of neutrons by the gravity field.There are also other experimental groups exploring similar setups.If successful, these experiments will produce neutrons in well defined ultralow energy states in the peV range which can be used for precise measurements of fundamental forces.We analyzed waveguides with 1D and 2D roughness.The ratios of the depletion times (line broadenings) in the biased diffusion approximation were the universal functions of the waveguide width and did not depend on the waveguide parameters.All relevant waveguide and roughness parameters collapsed into a single constant (essentially, a linewidth of the lowest quantum state), which was responsible for the exit neutron count.This constant strongly depended on the functional form of the roughness correlation function.We calculated this constant for various waveguides.In waveguides with 1D roughness the calculations could be carried out analytically for the most common types of the correlation functions; the 2D calculations were mostly numerical.
Our results were in good agreement with earlier experimental data despite the lack of experimental information about many important parameters.The predicted neutron count for the new experimental setup, for which the roughness profile was accurately measured, exhibited well-developed quantum step corresponding to consecutive depletion of the lower and lower gravitational states.Large amplitude of roughness in this setup could degrade the usability of the results.One of the possible ways to circumvent these difficulties and produce a much more controllable environment would be the use of a radically new design for a rough mirror which we called an Ising mirror [31,33].

1 Figure 1 :
Figure1: The ratios   (ℎ)/ 1 (ℎ)(7) for the lowest quantum levels  = 2; 3; 4; 5 as a function of the slit width ℎ.The curves are marked by the values of .These functions are the same for 1D and 2D roughness and in the biased diffusion approximation do not depend on the roughness parameters.

Figure 2 :
Figure 2: The dependence of the exit neutron count   / 0 (11) on the waveguide width ℎ for different values of Φ.The curves are marked by the values of Φ.The experimental data are fitted to the curve with Φ = 23.5 which corresponds to 1D Gaussian roughness with  = 0.19 and  = 1.19.

Figure 3 :
Figure 3: The inverse depletion times 1/  for the lowest gravitational states  = 1; 2; 3; 4; 5 for a waveguide with 1D exponential roughness with  =  = 1.The curves are marked by the values of .The vertical scale is in kHz.Equation (33) extends the results to other values of  and .

Figure 4 :
Figure 4: The same as Figure 3 for 2D exponential roughness in the new cell.

Figure 5 :
Figure 5: The prediction for the neutron count for the new waveguide with 2D exponential roughness with parameters  = 0.65,  = 1.02.