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

Authors: Y. Kamiya; G. Ichikawa; S. Komamiya
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

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.


---

## Body

## 1. 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 Colella [1], 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 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, a “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 a 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(1){-ħ22md2dz2+V(z)}ψn(z)=Enψn(z)
with the linear potential
(2)V(z)={mgz,z≥0,∞,z≤0,
where ψn and En are the eigenfunctions and eigenenergies, ħ 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
(3)(d2dξn2-ξn)ψn(ξn)=0,
where ξn≡z/z0-En/E0 and the characteristic length and energy are given by
(4)z0=(ħ22m2g)1/3~6μm,E0=(mg2ħ22)1/3~0.6peV.
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 Figure 1. Note that ξn=0, where zn≡z0En/E0, 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 Figure 1.Figure 1
Probability distributions of gravitationally bound states of neutrons for the first five states [35]. Horizontal lines indicate eigenenergies En and vertical lines show the classical turning points zn for each state.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 short-range 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 [10]. The limits for a CP-violating Yukawa-type potential were shown in [11, 12]. 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 [13, 14]. In this review, we discuss an experiment using this scheme to precisely measure the UCN position-space wave functions, performed by the group of Komamiya of the University of Tokyo [15]. Experiments exploiting other ideas for measuring the energy scale, that is, the energy differences between quantum states, are proposed by observing resonance transitions induced by a magnetic field [16–19] and mechanical vibrations [20]. These projects are called GRANIT [21] and qBounce, respectively. The first measurement of the resonance transition from the ground to the third state was reported by the group of Abele in [22]. Recently, a new limit for the CP-violating Yukawa-type potential using the resonance method was reported in [23]. It also shows a limit for the chameleon field [24–26], a dark energy candidate.


## 2. Precision Measurement and Issues Related to This Analysis

The precision measurement with a convex magnification mirror [15] was performed using a UCN source provided at ILL (Institut Laue-Langevin) [27]. 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 antivibration 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 Figure 2(b).Schematic drawing of the precision measurement experiment using a convex magnification mirror [15]. (a) is an overview of the system and (b) shows the geometry of the main components.
(a)
(b)


### 2.1. Collimating Guide

The guide settles gravitationally bound quantum states. To clearly distinguish each quantum state, which has energy differences of order ofΔE~1peV, a resolving time of Δt~ħ/ΔE~1msec 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 state. The bottom mirror is made of polished glass with roughness of Ra (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 roughness of Ra=0.4μm [28], 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. (The eigenfunctions modulate somewhat coherently, especially around the lower region, except for the ground state, as seen in Figure 1. To improve the contrast of the quantum modulation, a method using a negative step of several tens of μm, which transfers neutrons from the ground state to higher states, was considered in [9]. Results were reported in [29–31].) The neutron loss models for the scatterer or rough surface have been discussed in detail in [32–34], 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 [15]. 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
(5)Γn=γ∫h-2δh‍|ψn~|2dz,
where γ is a scaling constant, estimated from data to be 9.5-0.9+0.7×104s-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, nonspecular reflection, upscattering, and other processes are modeled empirically as
(6)Bn=βg22mEn~,
in which the loss rate is assumed to be proportional to the classical bouncing number per unit time, g/2v~z,n,max⁡, where β is a scaling constant (estimated to be 0.38-0.03+0.04), v~z,n,max⁡≡2En~/m is the maximum vertical velocity of a neutron in the nth state, and En~ is the eigenenergy of the deformed wave functions. The transmissivity of the guide for each state can be written as
(7)pn~∝〈exp⁡⁡[-lvx(Γn+Bn)]〉vx,
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 [15]) of neutrons for each state at the end of region III is estimated as in Figure 3(a).(a) is an estimated population distribution (probability distribution of the guide transmission) at the end of region III and (b) shows the corresponding Wigner phase space distribution [15].
(a)
(b)


### 2.2. 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 [35]. 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 Ra=1.9nm, two orders of magnitude smaller than the neutron wavelength.Figure 4
Magnification powerM of the cylindrical rod as a function of the height z [35]. It is about 20 times magnification around the lower region.To calculate the detected position distribution on the detector surface, the wave function at the end of the bottom floor is reexpressed by the Wigner phase space distribution [36, 37]:
(8)W(z,pz)≡12πħ∫-∞∞‍δηe-(i/ħ)pzη〈z+12η|ρ^|z-12η〉,
where pz 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 [38]. 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 [39]. Figure 3(b) shows the Wigner distribution constructed from the estimated populations in our experiment [15]. The time evolution of the Wigner distribution is calculated by the evolution of the density operator described by the Liouville-von Neumann equation
(9)∂ρ^∂t=-iħ[H^,ρ^],
where
(10)H^=pz2^2m+V(z^).
By evaluating the kinetic part as
(11)T=-pzm∂∂zW(z,pz)
and the potential part as
(12)∑l=0∞‍Ul=∑l=0∞‍(-1)l(ħ/2)2l(2l+1)!d2l+1V(z)dz2l+1∂2l+1∂pz2l+1W(z,pz),
one can obtainthe quantum Liouville equation for the Wigner distribution [40]
(13)(∂∂t+pzm∂∂z-dV(z)dz∂∂pz)W(z,pz)=∑l=1∞‍Ul.
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 [15].Figure 5
Black crosses show measured neutron position distribution [15]. 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 [15].


### 2.3. Position Sensitive Detector for UCNs

A back-thinned CCD (HAMAMATSU S7030-1008) with thin Ti-10B-Ti layers is used for the position sensitive detector in the experiment [15]. Its pixel size is 24μm×24μm and the thickness of the active volume is about 20μm. The 10B layer converts neutrons into charged particles by the nuclear reaction 10B(n,α)7Li. 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 thicknesses of the layers are 20 nm and 200 nm for Ti and 10B, 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 Figure 6) using very cold neutron beams at ILL [35]. For the details of the evaluation scheme, see [14, 41]. A neutron converter of 6Li using 6Li(n,α)3H reaction was also investigated. It is concluded that the use of 10B gives better spatial resolution [41].Figure 6
Fitting results of a Gd shadow pattern [35]. Spatial resolution is measured to be 3.35±0.09μm by evaluating these line spread functions. For details of evaluation method, see [14, 41].Other detectors using a silicon pixel device, Timepix [42], with 6LiF and 10B converter were studied in [43]. Its pixel pitch is 55μm×55μm and the thickness of its silicon layer is about 300μm. The spatial resolution for the 6LiF 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 10B 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, 29–31, 44]. Usually two fission fragments are emitted from a thin 235U 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 [44].


## 2.1. Collimating Guide

The guide settles gravitationally bound quantum states. To clearly distinguish each quantum state, which has energy differences of order ofΔE~1peV, a resolving time of Δt~ħ/ΔE~1msec 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 state. The bottom mirror is made of polished glass with roughness of Ra (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 roughness of Ra=0.4μm [28], 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. (The eigenfunctions modulate somewhat coherently, especially around the lower region, except for the ground state, as seen in Figure 1. To improve the contrast of the quantum modulation, a method using a negative step of several tens of μm, which transfers neutrons from the ground state to higher states, was considered in [9]. Results were reported in [29–31].) The neutron loss models for the scatterer or rough surface have been discussed in detail in [32–34], 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 [15]. 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
(5)Γn=γ∫h-2δh‍|ψn~|2dz,
where γ is a scaling constant, estimated from data to be 9.5-0.9+0.7×104s-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, nonspecular reflection, upscattering, and other processes are modeled empirically as
(6)Bn=βg22mEn~,
in which the loss rate is assumed to be proportional to the classical bouncing number per unit time, g/2v~z,n,max⁡, where β is a scaling constant (estimated to be 0.38-0.03+0.04), v~z,n,max⁡≡2En~/m is the maximum vertical velocity of a neutron in the nth state, and En~ is the eigenenergy of the deformed wave functions. The transmissivity of the guide for each state can be written as
(7)pn~∝〈exp⁡⁡[-lvx(Γn+Bn)]〉vx,
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 [15]) of neutrons for each state at the end of region III is estimated as in Figure 3(a).(a) is an estimated population distribution (probability distribution of the guide transmission) at the end of region III and (b) shows the corresponding Wigner phase space distribution [15].
(a)
(b)


## 2.2. 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 [35]. 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 Ra=1.9nm, two orders of magnitude smaller than the neutron wavelength.Figure 4
Magnification powerM of the cylindrical rod as a function of the height z [35]. It is about 20 times magnification around the lower region.To calculate the detected position distribution on the detector surface, the wave function at the end of the bottom floor is reexpressed by the Wigner phase space distribution [36, 37]:
(8)W(z,pz)≡12πħ∫-∞∞‍δηe-(i/ħ)pzη〈z+12η|ρ^|z-12η〉,
where pz 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 [38]. 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 [39]. Figure 3(b) shows the Wigner distribution constructed from the estimated populations in our experiment [15]. The time evolution of the Wigner distribution is calculated by the evolution of the density operator described by the Liouville-von Neumann equation
(9)∂ρ^∂t=-iħ[H^,ρ^],
where
(10)H^=pz2^2m+V(z^).
By evaluating the kinetic part as
(11)T=-pzm∂∂zW(z,pz)
and the potential part as
(12)∑l=0∞‍Ul=∑l=0∞‍(-1)l(ħ/2)2l(2l+1)!d2l+1V(z)dz2l+1∂2l+1∂pz2l+1W(z,pz),
one can obtainthe quantum Liouville equation for the Wigner distribution [40]
(13)(∂∂t+pzm∂∂z-dV(z)dz∂∂pz)W(z,pz)=∑l=1∞‍Ul.
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 [15].Figure 5
Black crosses show measured neutron position distribution [15]. 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 [15].


## 2.3. Position Sensitive Detector for UCNs

A back-thinned CCD (HAMAMATSU S7030-1008) with thin Ti-10B-Ti layers is used for the position sensitive detector in the experiment [15]. Its pixel size is 24μm×24μm and the thickness of the active volume is about 20μm. The 10B layer converts neutrons into charged particles by the nuclear reaction 10B(n,α)7Li. 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 thicknesses of the layers are 20 nm and 200 nm for Ti and 10B, 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 Figure 6) using very cold neutron beams at ILL [35]. For the details of the evaluation scheme, see [14, 41]. A neutron converter of 6Li using 6Li(n,α)3H reaction was also investigated. It is concluded that the use of 10B gives better spatial resolution [41].Figure 6
Fitting results of a Gd shadow pattern [35]. Spatial resolution is measured to be 3.35±0.09μm by evaluating these line spread functions. For details of evaluation method, see [14, 41].Other detectors using a silicon pixel device, Timepix [42], with 6LiF and 10B converter were studied in [43]. Its pixel pitch is 55μm×55μm and the thickness of its silicon layer is about 300μm. The spatial resolution for the 6LiF 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 10B 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, 29–31, 44]. Usually two fission fragments are emitted from a thin 235U 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 [44].


## 3. 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 of10μ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 [10–12, 23]. Furthermore one can start discussing a phase space formulation of quantum physics for the gravitationally bound quantum state [15]. An experiment for a possible test of phase space formalization using spatial interference is under preparation.



---

*Source: 859241-2014-09-10.xml*