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.
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)ħ22md2dz2ψn+mgzψn=Enψ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 En. 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 z0=(ħ2/2m2g)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 f0=mgz0/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 En to the state of energy Em is given by fnm=(En-Em)/2πħ. Solving the problem (1), one can show that
(2)fnm=f0(ϵ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 Section 2, 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 Figure 1. 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 vx velocity along the beam to be distributed with a mean value of v0=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 vx 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.
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).
UCNs are first prepared in an excited state by going down step (1) of height 15μm. The populations pn of the quantum states after the step are expected to be about p1=0.02,p2=p3=p4=0.3. Thus, the population of the ground quantum state is suppressed as compared to the populations of excited states.
Next, transitions between quantum states are induced with a periodic magnetic field gradient. The length of the transition region is L=16 cm, corresponding to an average passage time t0=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 f21=254 Hz, corresponding to the resonant horizontal velocity of v21=df21=2.54 m/s; for the 3→1 case, we expect f31=462 Hz and v31=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.
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.
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 field B→ 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→|=βzcos(ωt). The quantum mechanical excitation potential reads
(3)V^=sμβz^cos(ωt).
When the excitation frequency is close to a resonance, ω≈2πfnm, a Rabi oscillation between states |n〉 and |m〉 will take place at the angular frequency
(4)Ωnm=μħ〈n|z^|m〉β,
where the matrix elements of z^ can be expressed as (see, e.g., [14])
(5)〈n|z^|m〉=2z0(ϵn-ϵm)2(n≠m),〈n|z^|n〉=23z0ϵn.
To maximize the transition probability at resonance, the excitation strength β should verify Ωt0=π, where t0 is the excitation time. This condition can be expressed using (4) and (5) as
(6)βneeded=π2ħμz0(fnmf0)21t0.
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 Figure 1. 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 I1,I2,I3,I4,-I1,-I2,-I3,-I4,I1,…. 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 components Bx(x),Bz(x) and gradients ∂zBx(x),∂zBz(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|=Bx∂zBx+Bz∂zBz|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 I1=I4=1.4 A and I2=I3=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.
Magnetic field gradient ∂z|B| produced at the surface of the main mirror, without any external magnetic field.
Note that the result shown in Figure 2 assumes that no external field is applied. By applying an external field B→0, with, for example, B0,x=B0,z=1.5 mT and B0,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.
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 the x,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 B0,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 frequency f. 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 B0,y is applied in the transverse y direction. The purpose of this section is to calculate the minimum B0,y field to apply in order to guarantee the adiabaticity of neutron spin transport when passing by the magnetic excitation.
Magnetic field generated at the surface of the main mirror assuming the benchmark configuration I1=I4=1.4 A and I2=I3=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 velocity v along the x direction sees a time dependent magnetic field in its rest frame given by
(8)Bx(t)=B1cos(2πft+ϕ)sin(2πvtd),By(t)=B0,y,Bz(t)=-B1cos(2πft+ϕ)cos(2πvtd).
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 B1=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Π→dt=γΠ→×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 pmax, 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 B0,y,f,v, and ϕ, the criterion pmax 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 B0,y.
Numerically calculated neutron spin-flip probability during the passage below the wire array in the benchmark configuration I1=I4=1.4 A and I2=I3=3.5 A, as a function of the driving frequency, for different values of the external B0,y field.
As a conclusion of the numerical study, the value of the external field B0,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 B0,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 frequency f and fixed amplitudes I1=I4=1.4 A and I2=I3=3.5 A) and the external field (B0,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→|=β^B1cos2(2πft+ϕ)B12cos2(2πft+ϕ)+B0,y2,
where f is the driving frequency of the current in the wire array, B1=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).
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 Section 2, 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)〉=∑nan(t)|n〉,
which is the solution of the time dependent Schrödinger equation
(12)idandt=Enħan+∑ms2γβ(t)〈n|z^|m〉am.
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 |a1(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.
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 the 2→1 transition, with frequencies
(13)f+=141.5Hz,f-=113.5Hz,
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+β1cos(4πft+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)mg±=mg±μβ0.
Thus, we expect a spin-dependent resonance frequency of the n→m transition given by
(16)fnm±=(mg±)2/32π(2mħ)1/3(ϵn-ϵm)=fnm(1±μβ0mg)2/3,
where fnm 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 f21, namely,
(17)f21=((2f+)3/2+(2f-)3/2)2)2/3.
Remember that f+ and f- refer to driving frequencies, that correspond to excitation frequencies of 2f+ and 2f-. When applied to the maxima of the curve shown in Figure 7, one extracts f21=255.8 Hz which differs from the true resonance frequency f21,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 the x 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.
AppendixMagnetic Field of a Square Wire
Here, we provide formulas for the magnetic field generated by an infinitely long square wire with current I flowing uniformly in the wire in the y direction. The magnetic field lies in the (x,z) plane. We assume that (x=0,y=0) corresponds to the center of the wire. The length of the square is denoted by c (c=1 mm for the purpose of this paper). Formulas are valid outside the wire. We define the following quantities:
(A.1)xm=x-c2,zm=z-c2,xp=x+c2,zp=z+c2,L1=ln((xm2+zm2)(xp2+zp2)(xp2+zm2)(xm2+zp2)),L2=ln((xm2+zm2)(xp2+zm2)(xm2+zp2)(xp2+zp2)),L3=ln((xm2+zm2)(xm2+zp2)(xp2+zm2)(xp2+zp2)),Amm=arctan(xmzm),Apm=arctan(xpzm),Amp=arctan(xmzp),App=arctan(xpzp).
By integrating the Biot-Savart law over the volume of the wire, we find the following expressions for the field components:
(A.2)Bx(x,z)=-μ0I4πc2[xL1-c2L2+2zm(Amm-Apm)ggggggggg+2zp(App-Amp)c2],Bz(x,z)=μ0I4πc2[zL1-c2L3+2xm(Amp-Amm)gggggggg+2xp(Apm-App)c2].
By taking the derivative, we find the following expressions for the gradients:
(A.3)∂zBx=μ0I2πc2[Amp-Amm+Apm-App],∂zBz=μ0I4πc2L1.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
NesvizhevskyV. V.Near-surface quantum states of neutrons in the gravitational and centrifugal potentials201053764567510.3367/UFNe.0180.201007a.06732-s2.0-77958613661NesvizhevskyV. V.BörnerH. G.PetukhovA. K.AbeleH.BaeßlerS.RueßF. J.StöferleT.WestphalA.GagarskiA. M.GuennadyP.StrelkovA.Quantum states of neutrons in the Earth's gravitational field2002415686929729910.1038/415297a2-s2.0-18244395624AntoniadisI.BaesslerS.BuchnerM.Short-range fundamental forces20111275577810.1016/j.crhy.2011.05.004BraxP.PignolG.Strongly coupled chameleons and the neutronic quantum bouncer201110711130110.1103/PhysRevLett.107.111301BraxP.PignolG.RoulierD.Probing strongly coupled chameleons with slow neutrons20138808300410.1103/PhysRevD.88.083004JenkeT.CronenbergG.BurgdörferJ.Gravity resonance spectroscopy constrains dark energy and dark matter scenarios20141121515110510.1103/PhysRevLett.112.151105JenkeT.GeltenbortP.LemmelH.AbeleH.Realization of a gravity-resonance-spectroscopy technique20117646847210.1038/nphys19702-s2.0-79958014220NesvizhevskyV. V.ProtasovK. V.MooreD. C.Quantum states of neutrons in the Earth’s gravitational field: state of the art, applications, perspectives2006New York, NY, USANova Science65107BaesslerS.BeauM.KreuzM.The GRANIT spectrometer20111270772810.1016/j.crhy.2011.04.010KreuzM.NesvizhevskyV. V.Schmidt-WellenburgP.A method to measure the resonance transitions between the gravitationally bound quantum states of neutrons in the GRANIT spectrometer20096112-332633010.1016/j.nima.2009.07.059Schmidt-WellenburgP.Ultracold-neutron infrastructure for the gravitational spectrometer GRANIT200961126727110.1016/j.nima.2009.07.096ZimmerO.PiegsaF. M.IvanovS. N.Superthermal source of ultracold neutrons for fundamental physics experiments201110713480110.1103/PhysRevLett.107.134801RoulierD.VezzuF.BaesslerS.Status of the GRANIT facility2014201473043710.1155/2014/730437CaudauC.NesvizhevskyV. V.FertlM.PignolG.ProtasovK. V.Transitions between levels of a quantum bouncer induced by a noise-like perturbation2012677101310.1016/j.nima.2012.02.037