Precision Measurement of the Position-space Wave Functions of Gravitationally Bound Ultracold Neutrons

Gravity is the most familiar force at our natural length scale. However, it is still exotic from the view point of particle physics. The first experimental study of quantum effects under gravity was performed using a cold neutron beam in 1975. Following this, an investigation of gravitationally bound quantum states using ultracold neutrons was started in 2002. This quantum bound system is now well understood, and one can use it as a tunable tool to probe gravity. In this paper, we review a recent measurement of position-space wave functions of such gravitationally bound states, and discuss issues related to this analysis, such as neutron loss models in a thin neutron guide, the formulation of phase space quantum mechanics, and UCN position sensitive detectors. The quantum modulation of neutron bound states measured in this experiment shows good agreement with the prediction from quantum mechanics.


Introduction
Phenomena due to the gravitational field have been well understood at the macroscopic scales. However, there are only a few cases of experiments at microscopic scales, due to gravity's extreme weakness when compared to the other forces, such as electromagnetic and nuclear forces.
Even though quantum mechanics was established in the early 1900s, the first experiment to investigate a quantum effect under gravity was reported in 1975 by the group of R. Colella [1], in where a neutron interference pattern induced by a gravitational potential was observed. The major systematic uncertainty was attributed to a bending effect of an interferometer when one rotates the system to manipulate a relative quantum phase difference between two neutron paths. A subsequent series of experiments were carried out with careful studies of systematic effects by measuring the bending effect using X-rays and considering the Sagnac Effect [2] due to the rotation of the Earth [3,4]. The latest attempt was performed using normal-symmetric and skew-symmetric interferometers and a two-wavelength difference measurement method. The agreement with theoretical calculation is at the 1% level [5].
Another approach to observe quantum effects under gravity was carried out by the group of V. V. Nesvizhevsky [6,7]. In this experiment, one measured a transmission of ultracold neutrons (UCNs) through a vertically thin neutron guide in the terrestrial gravitational field as a function of the slit thickness, and showed the evidence of quantum bound states in the gravitational potential by analyzing the minimum height of a guide through which UCN can be transmitted. In a later detailed analysis, the characteristic sizes of position-space wave functions for the first and second quantum states were evaluated with an uncertainty of around 10% [8]. The major systematic uncertainty in this evaluation came from the difficulty in modeling the rough surface of a scatterer used on the ceiling of the guide, which causes some deformation of the wave functions. To minimize the effects of the deformation, an "differential method" was designed, in which one directly measures the wave functions using a position sensitive detector. In this method, a neutron guide of moderate height can be used and the deforming effect was kept rather small compared to the previous "integrated method". The first attempt of the differential method was reported in [8], where a plastic nuclear track detector (CR39) with an uranium coating [9] was used as the detector.
These quantum systems, in which the UCN bounces on a smooth floor, are expressed by the time-independent Schrödinger equation, with the linear potential where ψ n and E n are the eigenfunctions and eigenenergies,h is the reduced Planck constant, m is the neutron mass, and g is the standard gravitational acceleration. Equation (1) can be rewritten in dimensionless form as where ξ n ≡ z/z 0 − E n /E 0 and the characteristic length and energy are given by This is the Airy equation and solutions are described by Airy special functions, Ai(ξ n ) and Bi(ξ n ). Calculated probability distributions for the first five states are illustrated in Fig. 1. Note that ξ n = 0, where z n ≡ z 0 E n /E 0 , are the classical turning points of a classical bouncing ball, which are denoted in the cross points of each eigenenergy and the potential line of mgz in Fig. 1. This system is a suitable device with which to test the inverse square law of standard Newtonian gravity and to search for new gravity-like shortrange interactions around these scales. The first limit to a non-Newtonian force with a Yukawa-type interaction potential using this microscopic system was reported in [11]. The limits for a CP-violating Yukawa-type potential were shown in [12,13]. Several experimental schemes have been proposed to improve the sensitivity to such hypothetical new physics. One idea to achieve a better resolution on the measurement of the characteristic length is to utilize a convex reflection mirror to magnify the neutron distributions [14,15]. In this review, we discuss an experiment using this scheme to precisely measure the UCN position-space wave functions, performed by the group of S. Komamiya of the University of Tokyo [16]. Experiments exploiting other ideas for measuring the energy scale, i.e. the energy differences between quantum states, are proposed by observing resonance transitions induced by a magnetic field [17][18][19][20] and mechanical vibrations [21]. These projects are called GRANIT and qBounce, respectively. The first measurement of the resonance transition from the ground to the third state was reported by the group of H. Abele in [22] 1 . For details of those projects, please see review papers [23] and [24] in this special issue.

Precision measurement and issues related to this analysis
The precision measurement with a convex magnification mirror [16] was performed using a UCN source provided at ILL(Institut Laue-Langevin) [25]. The velocity distribution was measured by a TOF method and is well characterized by a Gaussian distribution with mean of 9.4 m/s and standard deviation of 2.8 m/s. Figure 2(a) shows a schematic drawing of the experiment. The entire system was mounted on an anti-vibration table and magnetically shielded by a Permalloy sheet. The main components can be separated into three parts: a vertically thin collimating guide, a magnification mirror(rod), and a pixelated position sensitive detector. The layout of these components is shown in Fig. 2

Collimating guide
The guide settles gravitationally bound quantum states. To clearly distinguish each quantum state, which have energy differences of order of ∆E ∼ 1 peV, a resolving time of ∆t ∼h/∆E ∼ 1 msec is required. Considering that the horizontal velocities of UCNs which we use are less than 15 m/s, the collimating guide was designed to have a sufficient length (192 mm in our setup) to form each quantum states. The bottom mirror is made of polished glass with a roughness of R a (arithmetic average) = 0.03 µm. The ceiling scatterer is a Gd-Ti-Zr alloy (54/35/11) deposited on glass, with a Fermi potential tuned to be nearly zero, and a roughness of R a = 0.4 µm [30], which scatters out neutrons in the higher states. The guide selects lower states with appropriate populations to improve the contrast of the quantum spatial modulation 2 . The neutron loss models for the scatterer or rough surface have been discussed in detail in [31][32][33], and they are still interesting issues not only for UCN guiding applications but also for storage experiments such as neutron EDM measurements. In the experiment reviewed in this paper, empirical models of the loss rates are adopted [16]. The loss rate by the scatterer, Γ n , is assumed to be proportional to the probability of finding neutron in the roughness region, and is given by where γ is a scaling constant, estimated from data to be 9.5 +0.7 −0.9 × 10 4 s −1 , h is the height of the guide (100 µm), δ is the roughness of the scatterer (0.4 µm), andψ n are deformed wave functions in the guide. Neutron losses at the bottom mirror due to absorption, non-specular reflection, up scattering and other processes are modeled empirically as in which the loss rate is assumed to be proportional to the classical bouncing number per unit time, g/2ṽ z,n,max , where β is a scaling constant (estimated to be 0.38 +0.04 −0.03 ),ṽ z,n,max ≡ 2Ẽ n /m is the maximum vertical velocity of a neutron in the nth state, andẼ n is the eigenenergy of the deformed wave functions. The transmissivity of the guide for each state can be written as where < ... > vx indicate the average over the neutron horizontal velocities. By applying a diabatic transition from region II to III, the population distribution (a probability in the paper [16]) of neutrons for each state at the end of region III is estimated as in Fig. 3(a).  [16].

Magnification mirror
A Ni coated cylindrical rod is used as a magnifying convex mirror. Its radius is 3 mm and it is placed to have a grazing angle of 20 deg. for horizontally moving neutrons at the bottom floor level (z = 0). Figure 4 shows the magnification power as a function of the height [10]. It gives about 20 times magnification around the lower region of z ∼ 20 µm. Before depositing the Ni coating, the cylinder was polished at the Research Center for Ultra-Precision Science and Technology, Osaka University. The roughness of the rod after the depositing was measured to be R a = 1.9 nm, two orders of magnitude smaller than the neutron wavelength.
To calculate the detected position distribution on the detector surface, the wave function at the end of the bottom floor is re-expressed by the Wigner phase space distribution [37,38], where p z is the momentum andρ is a density operator. The Wigner distribution is known as a phase space formulation of quantum mechanics, and is widely used for quantum optics, for example in the study of decoherence [39].
As an application to massive particles, one can find a paper which shows a phase-space tomography of the Wigner distribution for a coherent atomic beam in a double-slit experiment [40]. Figure 3(b) shows the Wigner distribution constructed from the estimated populations in our experiment [16]. The time evolution of the Wigner distribution is calculated by the evolution of the density operator described by the Liouville -von Neumann equation whereĤ By evaluating the kinetic part as and the potential part as one can obtain the Quantum Liouville Equation for the Wigner distribution In the case of V (z) = mgz, the right hand side of the equation vanishes and it becomes the classical Liouville equation. Therefore, in region IV, one can treat the evolution of each phase point of the Wigner distribution as a classical path under gravity. Figure 5 shows the measured data and the best theoretical estimation using this model. The corresponding p-value is 0.715, and the experimental data support the quantum features described by the phase space formulation using the Wigner distribution [16]. Position at the detector surface (mm) Figure 5: Black crosses show measured neutron position distribution [16]. Modulated distribution was clearly measured. The result shows good agreement with the quantum expectation calculated using phase space formalization by Wigner distribution (p-value is 0.715). The best estimated line is shown by a red solid curve and the other lines indicate contributions of each state to the distribution. For definitions of fitting parameters, best fit values, systematic uncertainties, and the other details, see [16].

Position sensitive detector for UCNs
A back-thinned CCD (HAMAMATSU S7030-1008) with thin Ti-10 B-Ti layers is used for the position sensitive detector in the experiment [16]. Its pixel size is 24 µm × 24 µm and the thickness of the active volume is about 20 µm. The 10 B layer converts neutrons into charged particles by the nuclear reaction 10 B(n,α) 7 Li. The secondary particles are emitted in a nearly back-to-back configuration. One of them deposits its kinetic energy in the active area and creates a charge cluster, which typically spreads into nine pixels. The weighted center of the charge cluster is a good estimation of the incident neutron position. The thickness of the layers are 20 nm and 200 nm for Ti and 10 B, respectively, and they are formed by evaporating directly on the CCD surface. The spatial resolution is measured to be 3.35 ± 0.09 µm by evaluating the line spread function (LSF) (see Fig. 6) using very cold neutron beams at ILL [10]. For the detail of the evaluation scheme, see [15,42]. A neutron converter of 6 Li using 6 Li(n,α) 3 H reaction was also investigated. It is concluded that the use of 10 B gives better spatial resolution [42].  Figure 6: Fitting results of a Gd shadow pattern [10]. Spatial resolution is measured to be 3.35 ± 0.09 µm by evaluating these line spread functions. For details of evaluation method, see [15,42].
Other detectors using a silicon pixel device, Timepix [43], with 6 LiF and 10 B converter were studied in [44]. Its pixel pitch is 55 µm × 55 µm and the thickness of its silicon layer is about 300 µm. The spatial resolution for the 6 LiF converter was evaluated by LSF to be 2.3 µm, corresponding to about 5.3 µm in FWHM of the point spread function (PSF). The performance of the 10 B converter is estimated by Monte-Carlo simulations to be better than 3 µm in FWHM of the PSF.
Another detector concept of Uranium coated plastic nuclear tracker (CR39) [9] was used in several experiments [8,[34][35][36]45]. Usually two fission fragments are emitted from a thin 235 U coating, and one of the daughter nuclei makes a track of defects in the CR39. By chemical etching, the diameter of track points is increased up to 1 µm, allowing us to scan these vertexes using an optical microscope. Position resolution is around 1 µm [9]. By carefully analyzing the vertex shape, the spatial resolution can be improved to 0.7 µm [45].

Summary
The quantum system of a gravitationally bound neutron is one of the most suitable tools to investigate gravity or gravity-like hypothetical interactions around the scale of 10 µm in length or 1 peV in energy. After establishing this research field by the pioneering work in observation of the quantum state [6], experimental schemes for precision measurements of these characteristic scales have developed rapidly, and nowadays, one can establish limits for parameter spaces of new physics experimentally [11][12][13]26]. Furthermore one can start discussing a phase space formulation of quantum physics for the gravitationally bound quantum state [16]. An experiment for a possible test of phase space formalization using spatial interference is under preparation.