# Gravitational Resonance Spectroscopy with an Oscillating Magnetic Field Gradient in the GRANIT Flow through Arrangement

Authors: G. Pignol; S. Baeßler; V. V. Nesvizhevsky; K. Protasov; D. Rebreyend; A. Voronin
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

Gravitational resonance spectroscopy consists in measuring the energy spectrum of bouncing ultracold neutrons above a mirror by inducing resonant transitions between different discrete quantum
levels. We discuss how to induce the resonances with a flow through arrangement in the GRANIT spectrometer, excited by an oscillating magnetic field gradient. The spectroscopy could be realized in two distinct modes  (so called DC and AC) using the same device to produce the magnetic excitation. We present calculations demonstrating the feasibility of the newly proposed AC mode.


---

## Body

## 1. Introduction

Ultracold neutrons bouncing over a horizontal mirror are used to probe quantum effects of a particle in the gravitational field [1]. The vertical motion of such neutrons bouncing at submillimeter distance from the mirror has discrete energy spectrum that can be calculated from the stationary Schrödinger equation:
(1)
ħ
2
2
m
d
2
d
z
2
ψ
n
+
m
g
z
ψ
n
=
E
n
ψ
n
,
where m is the neutron mass, g
=
9.81 m/s2 is the local gravitational acceleration, and ψ
n
(
z
)
n
=
1,2
,
… are the stationary wavefunctions with associated energy E
n. The existence of the quantization of the vertical motion was demonstrated a decade ago [2], profiting from the relatively large spatial extension of the ground state wavefunction characterized by z
0
=
(
ħ
2
/
2
m
2
g
)
1
/
3
≈
5.87 μm.Precision study of the quantum states is motivated by their sensitivity to extra short range interactions (see [3] and references therein), in particular, those induced by Chameleon Dark Energy [4–6]. In addition, the neutron quantum states provide a unique test of the weak equivalence principle in a quantum regime, since the inertial and gravitational masses in (1) do not cancel.High precision measurements can be achieved with the gravitational resonance spectroscopy technique, where transitions between quantum states are induced by a periodic excitation [7]. The characteristic frequencies of the transitions, as low as f
0
=
m
g
z
0
/
2
π
ħ
≈
145 Hz, are accessible by electrical as well as mechanical oscillators. More precisely, the quantum frequencies for a transition n
→
m from the state of energy E
n to the state of energy E
m is given by f
n
m
=
(
E
n
-
E
m
)
/
2
π
ħ. Solving the problem (1), one can show that
(2)
f
n
m
=
f
0
(
ϵ
n
-
ϵ
m
)
,
where ϵ
n
=
[
2.338,4.088,5.521,6.787
,
…
] is the series of the negative zeros of the Airy function. To perform the spectroscopy of the bouncing neutron, the interaction that couples different quantum states can be a vibration of the bottom mirror or an oscillating magnetic field gradient [8]. The former has been used by the QBounce collaboration [7]; the latter will be used in the GRANIT spectrometer [9].As a first step of GRANIT, a flow through measurement of the resonant transitions magnetically excited between the first three quantum states will be realized as first proposed in [10]. According to this proposal, a space periodic (but static) magnetic field gradient will be generated at the surface of the bottom mirror. The frequency of the excitation seen by a neutron will thus vary according to its horizontal velocity. A detailed analysis of this scheme is provided in [9]. Another mode of operation (the AC mode) can be implemented with the same setup, consisting in generating a homogeneous gradient but oscillating in time, which should allow a more direct probe of the resonances.The paper is organized as follows: in Section2, we present the flow through setup, in Section 3, we describe in some details the magnetic excitation, in Section 4, we calculate the conditions for the adiabaticity of spin transport, and, in Section 5, we present a theoretical description of the magnetically induced transitions in the AC mode.


## 2. The GRANIT Flow through Setup

The sketch of the flow through setup is shown in Figure1. Ultracold neutrons are produced from a dedicated superthermal source installed at a cold beamline of the high flux reactor of the Institut Laue-Langevin. The source relies on downscattering of neutrons with a wavelength of 0.89 nm in a superfluid helium bath cooled down to 0.8 K [11, 12]. According to a preliminary measurement of the UCN velocity spectrum [13], we expect the v
x velocity along the beam to be distributed with a mean value of v
0
=
4 m/s and a standard deviation of 1.5 m/s. In the rest of this paper, we will assume a Gaussian profile for the v
x distribution with those parameters. UCNs are extracted using a narrow slit to accept only those with practically no vertical velocity. The spectroscopy is performed with four steps: (1) state preparation, (2) resonant transition, (3) state analysis, and (4) detection of transmitted flux.Figure 1
Sketch of the flow through setup. Ultracold neutrons enter from the left; they go through the step to depopulate the ground quantum state (1), 16 cm long transition region (2), 9 cm long analyzer (3), and detector (4).(i)
UCNs are first prepared in an excited state by going down step (1) of height15 μm. The populations p
n of the quantum states after the step are expected to be about p
1
=
0.02
,
p
2
=
p
3
=
p
4
=
0.3. Thus, the population of the ground quantum state is suppressed as compared to the populations of excited states.
(ii)
Next, transitions between quantum states are induced with a periodic magnetic field gradient. The length of the transition region isL
=
16 cm, corresponding to an average passage time t
0
=
40 ms. Two different schemes could be implemented in principle: the AC excitation and the DC excitation. In the DC mode, the field gradient is static and spatially oscillating in the x direction with a period of d
=
1 cm. In this case, only neutrons with specific horizontal velocities meet the resonance condition. The deexcitation 2
→
1 is expected to be induced by an excitation frequency f
21
=
254 Hz, corresponding to the resonant horizontal velocity of v
21
=
d
f
21
=
2.54 m/s; for the 3
→
1 case, we expect f
31
=
462 Hz and v
31
=
4.62 m/s. In the AC mode, the field gradient is spatially uniform, oscillating in time. One would then find the resonances by directly scanning the excitation frequency.
(iii)
A second horizontal mirror above the main mirror serves as a state analyzer. For a slit opening of about 25μm, only ground state neutrons are accepted; higher quantum states are rejected. The length of the analyzer in the x direction is 9 cm.
(iv)
Finally, neutrons are detected at the exit of the analyzer. In AC mode, the flux of transmitted neutrons should display a resonance pattern as a function of the excitation frequency. In DC mode, one has to measure the horizontal velocity of the transmitted neutrons to deduce the resonant frequency. This is achieved by measuring the height of the neutrons with a position sensitive detector after a free fall distance of 30 cm.Before developing the details, let us estimate the strength of the needed oscillating magnetic field gradient. The interaction of a neutron with a magnetic fieldB
→ is described by the operator -
μ
→
^
·
B
→, where μ
→
^ is the neutron magnetic moment observable acting on the internal spin degree of freedom. When the variation of the magnetic field is slow compared to the Larmor frequency, the spin will follow the direction of the magnetic field. This adiabaticity condition will be addressed in Section 4. In this case, the neutron trajectory and the spin dynamics are effectively decoupled. Then, the motion of the neutron is determined by the potential s
μ
|
B
→
| where μ
=
60.3 neV/T is the magnetic moment of the neutron and s
=
1 for “spin up” neutrons and s
=
-
1 for “spin down” neutrons. Classically, a vertical force is applied on the neutron by the field gradient ∂
z
|
B
→
|. Let us now assume a magnetic excitation of the form |
B
→
|
=
β
z
cos
⁡
(
ω
t
). The quantum mechanical excitation potential reads
(3)
V
^
=
s
μ
β
z
^
cos
⁡
(
ω
t
)
.
When the excitation frequency is close to a resonance, ω
≈
2
π
f
n
m, a Rabi oscillation between states |
n
〉 and |
m
〉 will take place at the angular frequency
(4)
Ω
n
m
=
μ
ħ
〈
n
|
z
^
|
m
〉
β
,
where the matrix elements of z
^ can be expressed as (see, e.g., [14])
(5)
〈
n
|
z
^
|
m
〉
=
2
z
0
(
ϵ
n
-
ϵ
m
)
2
(
n
≠
m
)
,
〈
n
|
z
^
|
n
〉
=
2
3
z
0
ϵ
n
.To maximize the transition probability at resonance, the excitation strengthβ should verify Ω
t
0
=
π, where t
0 is the excitation time. This condition can be expressed using (4) and (5) as
(6)
β
needed
=
π
2
ħ
μ
z
0
(
f
n
m
f
0
)
2
1
t
0
.
One finds a needed field gradient of β
=
0.22 T/m to induce the 2
→
1 transition and β
=
0.74 T/m to induce the 3
→
1 transition.


## 3. The Magnetic Excitation

The magnetic field excitation will be generated by an array of 128 copper wires with square section arranged as shown in Figure1. In practice, the system is constituted of four modules, each one holding 32 adjacent wires. A wire has a section of 1 mm2 and a length of 30 cm in the y direction. Adjacent wires are separated by a gap of 0.25 mm. Electrical connectors are arranged so that the following 8-periodic pattern current could be applied I
1
,
I
2
,
I
3
,
I
4
,
-
I
1
,
-
I
2
,
-
I
3
,
-
I
4
,
I
1
,
…. Thus, the magnetic field will be 1 cm periodic. The system will be placed above the horizontal mirror in the transition region as shown in Figure 1 at a distance of 0.8 mm from the mirror.The magnetic field generated by a single infinitely long square wire can be calculated analytically; the corresponding formulas are reported in the appendix. The magnetic field componentsB
x
(
x
)
,
B
z
(
x
) and gradients ∂
z
B
x
(
x
)
,
∂
z
B
z
(
x
) at the surface of the mirror are obtained by summing the corresponding quantities for each of the 128 wires weighted by the individual currents. Then, the field gradient is calculated according to
(7)
∂
z
|
B
|
=
B
x
∂
z
B
x
+
B
z
∂
z
B
z
|
B
|
.
It is possible to tune the currents to obtain a homogeneous gradient at the surface of the mirror. We show, in Figure 2, the result for I
1
=
I
4
=
1.4 A and I
2
=
I
3
=
3.5 A where a field gradient of 0.52 T/m is obtained, as needed to induce resonant transitions between quantum states. This will be the benchmark configuration for the rest of the paper. The residual “noise” seen in Figure 2 has an amplitude of 0.02 T/m; this noise would increase for a wire array closer to the mirror. The frequency of this noise seen by a neutron passing at 4 m/s is about 2 kHz, much higher than the frequencies of interest for resonant transitions between low lying quantum states.Figure 2
Magnetic field gradient∂
z
|
B
| produced at the surface of the main mirror, without any external magnetic field.Note that the result shown in Figure2 assumes that no external field is applied. By applying an external field B
→
0, with, for example, B
0
,
x
=
B
0
,
z
=
1.5 mT and B
0
,
y
=
0, the situation changes dramatically as shown in Figure 3. With a strong external field applied, a gradient oscillating in the x direction with a period of 1 cm is generated.Figure 3
Magnetic field gradient∂
z
|
B
| produced at the surface of the main mirror with external magnetic field applied in the x
,
z plane.As a result, the wire array is a versatile device to generate the field gradient that can be used for the AC excitation mode as well as for the DC excitation mode. In the DC mode, we apply DC current in the wire array and apply a strong external field in thex
,
z direction. The vertical force exerted by the field gradient on the neutron will oscillate in space. In the AC mode, we apply AC current in the wire array and a small external field B
0
,
y to satisfy the adiabaticity of spin transport, as detailed in the next section.


## 4. The Adiabaticity of Spin Transport

The magnetic field gradient produced by the array of wires described in the previous section will exert a force on the passing neutrons. The sign of the force depends on the relative orientation between the neutron spin and the magnetic field. To induce resonant transitions between quantum states, a neutron should feel a well-defined oscillating vertical force. Thus, one must make sure that the magnetic field is strong enough to hold the neutron spins parallel or antiparallel to the magnetic field at any time. If the adiabaticity condition for spin transport along the wire array is fulfilled, then the spin dynamics and the neutron trajectory are decoupled.In the AC mode of excitation, the input current in the wire array is oscillating with driving frequencyf. The magnetic field amplitude created by the wire array at the surface of the mirror is given by the pattern shown in Figure 4, with the whole pattern oscillating in time at the driving frequency f. Thus, the magnetic field generated by the wire array will not be sufficient to hold the neutron spin, since the magnitude of the field crosses zero at a frequency f. To maintain a nonzero value of the field magnitude at any time, a static, homogeneous external field B
0
,
y is applied in the transverse y direction. The purpose of this section is to calculate the minimum B
0
,
y field to apply in order to guarantee the adiabaticity of neutron spin transport when passing by the magnetic excitation.Figure 4
Magnetic field generated at the surface of the main mirror assuming the benchmark configurationI
1
=
I
4
=
1.4 A and I
2
=
I
3
=
3.5 A.Here, we calculate the spin dynamics only, assuming that a neutron pass below the wire array in a straight horizontal line trajectory, at the surface of the bottom mirror. A neutron with velocityv along the x direction sees a time dependent magnetic field in its rest frame given by
(8)
B
x
(
t
)
=
B
1
cos
⁡
(
2
π
f
t
+
ϕ
)
sin
(
2
π
v
t
d
)
,
B
y
(
t
)
=
B
0
,
y
,
B
z
(
t
)
=
-
B
1
cos
⁡
(
2
π
f
t
+
ϕ
)
cos
⁡
(
2
π
v
t
d
)
.
It results from the combination of the oscillation of the field in space with period d and the oscillation in time with frequency f. We will set B
1
=
0.8 mT for the benchmark wire currents described in the previous section.The spin dynamics is given by the Bloch equation for the polarization vectorΠ
→:
(9)
d
Π
→
d
t
=
γ
Π
→
×
B
→
(
t
)
,
where γ
=
2
μ
/
ħ
=
183 kHz/mT is the neutron gyromagnetic ratio.We have solved numerically the Bloch equation using a Runge Kutta solver. The initial condition for the polarization vectorΠ
→
(
0
) was set to the unit vector aligned with B
→
(
0
), describing a “spin up” neutron. We define the spin-flip probability at time t as p
(
t
)
=
(
1
-
Π
→
·
B
→
/
|
B
→
|
)
/
2. With this definition, p
(
t
)
=
0 if the spin is aligned with the magnetic field at time t and p
(
t
)
=
1 if the spin has reversed its direction relative to the magnetic field at time t. As an adiabaticity criterion, we take p
max
⁡, the maximum spin-flip probability during the passage of a neutron below the wire array of duration L
/
v. For a given set of parameters B
0
,
y
,
f
,
v, and ϕ, the criterion p
max
⁡ was numerically calculated. The result was then averaged over the phase ϕ and the velocity spectrum v. The final result is presented in Figure 5 as a function of the driving frequency f, for different values of the external field B
0
,
y.Figure 5
Numerically calculated neutron spin-flip probability during the passage below the wire array in the benchmark configurationI
1
=
I
4
=
1.4 A and I
2
=
I
3
=
3.5 A, as a function of the driving frequency, for different values of the external B
0
,
y field.As a conclusion of the numerical study, the value of the external fieldB
0
,
y
=
0.3mT is sufficient to hold the neutron spin with an accuracy better than one percent, in the frequency range of interest between 0 and 300 Hz. In the following calculations, the external holding transverse field will thus be set to B
0
,
y
=
0.3 mT.


## 5. The Resonant Transitions in the AC Mode

We have now defined a magnetic configuration for the AC mode with the wire array (oscillating currents at variable frequencyf and fixed amplitudes I
1
=
I
4
=
1.4 A and I
2
=
I
3
=
3.5 A) and the external field (B
0
,
y
=
0.3 mT) that (i) holds the neutron spin and (ii) generates an oscillating gradient with the required amplitude to induce resonant transitions between quantum states. Note that the time dependent gradient β
(
t
) seen by the neutrons is not perfectly harmonic. Following (7), the expression of the time dependent gradient is
(10)
β
(
t
)
=
∂
z
|
B
→
|
=
β
^
B
1
cos
⁡
2
(
2
π
f
t
+
ϕ
)
B
1
2
cos
⁡
2
(
2
π
f
t
+
ϕ
)
+
B
0
,
y
2
,
where f is the driving frequency of the current in the wire array, B
1
=
0.8 mT, and β
^
=
0.52 T/m. We plot β
(
t
) in Figure 6, where it is apparent that the excitation frequency (the frequency of the β
(
t
) excitation) is twice the driving frequency f (the frequency of the oscillating currents in the wire array).Figure 6
Time dependent gradientβ
(
t
) for a driving frequency of f
=
100 Hz, corresponding to an excitation frequency of 200 Hz. Solid black line: (10), dashed red line: first order Fourier expansion equation (14).We will now simulate the transition probabilities that could be observed in the GRANIT flow through setup as described in Section2, assuming the benchmark magnetic configuration resulting in the excitation (10). The problem consists in calculating the time evolution of a neutron quantum state
(11)
|
ψ
(
t
)
〉
=
∑
n
‍
a
n
(
t
)
|
n
〉
,
which is the solution of the time dependent Schrödinger equation
(12)
i
d
a
n
d
t
=
E
n
ħ
a
n
+
∑
m
‍
s
2
γ
β
(
t
)
〈
n
|
z
^
|
m
〉
a
m
.
Here, we assume that the adiabaticity condition is fulfilled. The initial condition of the state is chosen to be |
ψ
(
0
)
〉
=
|
2
〉 immediately after the preparation step. Then, we solved (12) with a Runge-Kutta algorithm for a given set of parameters f
,
v
,
s, and ϕ, where the sum is restricted to the first four quantum states. For a given horizontal velocity v, the probability |
a
1
(
L
/
v
)
|
2 for the neutron to be detected in the ground state at the exit of the magnetic excitation is calculated. The result is then averaged over the excitation phase ϕ, the spin state s
=
±
1, and the velocity spectrum. We plot in Figure 7 the transition probability as a function of the driving frequency f.Figure 7
Numerical solution of the Schrödinger equation (12). The probability of the 2
→
1 transition is plotted as a function of the driving frequency.According to this numerical calculation, we expect to see two resonances in the transmitted UCN flux associated with the2
→
1 transition, with frequencies
(13)
f
+
=
141.5
H
z
,
f
-
=
113.5
H
z
,
where f
+ and f
- are the maxima of the resonance curve corresponding to spin “up” neutrons (s
=
1) and spin “down” neutrons (s
=
-
1), respectively.The splitting of the resonances depending on the spin state could be interpreted as a “Stern-Gerlach” split due to a constant magnetic field gradient. To see this, it is useful to perform the Fourier expansion of the gradient excitationβ
(
t
) given by (10)
(14)
β
(
t
)
=
β
0
+
β
1
cos
⁡
(
4
π
f
t
+
2
ϕ
)
+
⋯
,
where β
0
=
0.289 T/m and β
1
=
0.228 T/m. The constant term β
0 of the excitation should be thought of as an spin-dependent effective modification of g. This is done by identification of the total constant vertical force (gravity plus constant gradient) to an effective gravitational force:
(15)
m
g
±
=
m
g
±
μ
β
0
.
Thus, we expect a spin-dependent resonance frequency of the n
→
m transition given by
(16)
f
n
m
±
=
(
m
g
±
)
2
/
3
2
π
(
2
m
ħ
)
1
/
3
(
ϵ
n
-
ϵ
m
)
=
f
n
m
(
1
±
μ
β
0
m
g
)
2
/
3
,
where f
n
m is the unperturbed transition frequency. The splitting of the resonances obtained by the numerical resolution of the full Schrödinger equation (12) is in quantitative agreement with (16).In addition, (16) motivates a combination of the two resonant frequencies f
+ and f
- to extract the unperturbed transition frequency f
21, namely,
(17)
f
21
=
(
(
2
f
+
)
3
/
2
+
(
2
f
-
)
3
/
2
)
2
)
2
/
3
.
Remember that f
+ and f
- refer to driving frequencies, that correspond to excitation frequencies of 2
f
+ and 2
f
-. When applied to the maxima of the curve shown in Figure 7, one extracts f
21
=
255.8 Hz which differs from the true resonance frequency f
21
,
true
=
253.8 Hz given by (2) by 2 Hz. In fact, there are several features of the full problem (given by (12)) that the simple estimate (16) does not catch. The simplified formula (16) can be obtained by assuming a two-level system (states |
1
〉 and |
2
〉) excited by a harmonic force given by (14) that couple states |
1
〉 and |
2
〉, neglecting the self-couplings of the types 〈
1
|
z
^
|
1
〉 and 〈
2
|
z
^
|
2
〉. The full calculation (12) takes into account the nonharmonic excitation given by (10) that couples all states |
1
〉, |
2
〉, |
3
〉, and |
4
〉 including self-couplings. All these complications are potential sources of frequency shifts of the resonance line. We then conclude that these shifts are below the percent level.


## 6. Conclusion

In the GRANIT flow through arrangement, two possible modes to induce resonant transitions between the quantum states could be used. In the DC mode, where the magnetic field gradient oscillates spatially along thex direction, the excitation frequency is controlled by the horizontal neutron velocity. In the AC mode, where the gradient is homogeneous in space and oscillates in time, the excitation frequency is selected directly by the frequency of the current driving the magnetic excitation. We have shown that the condition of adiabaticity of spin transport can also be fulfilled in the AC mode using a moderate horizontal magnetic field normal to the neutron propagation axis. Finally, a calculation of the expected resonance line for the 2
→
1 transition indicates that a measurement of the transition frequency at a precision better than a percent is possible. A detailed comparison describing the relative merits of the two methods with the associated systematic effects is left for a future work.



---

*Source: 628125-2014-08-24.xml*