New Experimental Constraint on an Exotic Spin-dependent Interaction
Young Jin Kim,* Ping-Han Chu,† and Igor Savukov
Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
We present experimental results to probe an exotic spin-dependent interaction for polarized electrons. The
experiment employs a spin-exchange relaxation-free (SERF) magnetometer, as both a source of 100 % polarized
electrons and a detector. In this experimental approach, an interaction between the SERF polarized electrons and
unpolarized atoms of a solid-state mass was studied. We report a first experimental constraint on the interaction in
the interaction range below 1 cm.
The extreme small value of the neutron electric dipole
moment suggested the existence of new fundamental
bosons [1-3] such as the axion [1], axion-like-particles
(ALPs), dark photons [4,5], which are prominent darkmatter candidates [6]. The exotic particles are very light
and can mediate new macroscopic interactions between
fermions, beyond electromagnetism and gravity, with an
interaction range from micrometers to centimeters. There
are fifteen possible such interactions that contain static
spin-dependent operators or both spin- and velocitydependent operators [7,8]. Some of the interactions are
not invariant under parity (P) or time-reversal (T)
symmetries [9]; therefore their discovery would provide
new sources for P and T symmetry violations, which are
essential for the matter-antimatter asymmetry of the
Universe that cannot be explained by the Standard Model
of particle physics.
To explore all the fifteen exotic spin-dependent
interactions, we proposed a novel experimental approach
based on a spin-exchange relaxation-free (SERF)
magnetometer [9]. The SERF magnetometer is the most
sensitive cryogenic-free magnetic field sensor reaching
femotesla sensitivity [10] and contains a large number of
almost 100% optically polarized electrons of Rb atoms
(1014 cm-3 at about 160°C cell temperature) located in a
vapor cell. The interactions between the polarized Rb
electrons and atoms of a solid-state mass can tilt the Rb
electron spins, similar as an external magnetic field. By
detecting the tilt of the polarized Rb spins, the exotic
interactions can be precisely determined or constrained.
Of the fifteen exotic spin-dependent interactions, there
are three possible ones between polarized electrons and
unpolarized atoms [9] that can be written as (in SI units,
adopting the numbering schemes in [8,11]):
V4+5 = − f4+5
$1 1'
!
!2
[σˆ i ⋅ (υ × r̂)] &% + 2 )( e−r/λ ,
8π me c
λr r
(1)
"1 1%
!2 ˆ
(σ i ⋅ r̂) $ + 2 ' e−r/λ ,
# λr r &
8π me
(2)
V9+10 = f9+10
*
†
Email address: [email protected]
Email address: [email protected]
V12+13 = f12+13
! "1%
!
(σˆ i ⋅ υ ) $ ' e−r/λ ,
#r&
8π
(3)
where ! is Planck’s constant, me is the mass of the
polarized electrons, c is the speed of light in vacuum, σˆ i
is the spin vector of the ith polarized electron while
!
!
σ i = !σˆ i / 2 , r̂ = r / r is a unit vector in the direction
!
between the polarized electrons and unpolarized atoms, υ
is their relative velocity vector, and λ = ! / mALP c is the
interaction range, mALP is the boson mass. Here fi is the
coupling strength constant for the interaction Vi . While
the static spin-dependent interaction V9+10 for polarized
electrons has been experimentally investigated and its
experimental limits have been set [12-15], the other
interactions V4+5 and V12+13 for polarized electrons,
dependent on both spins and their relative velocity, have
never been constrained by experiments. For that reason,
we first focus on the interaction V4+5 using our proposed
experiments [9], which will be a subject in the remainder
of this paper.
The experimental setup to investigate the V4+5
interaction is shown in Fig. 1. For an unpolarized mass,
we used a 2 × 2 × 2 cm 3 cube-shaped nonmagnetic
bismuth germinate insulator (BGO, Bi4Ge3O12) which is
advantageous for the high number density of nucleons
( 4.3×10 24 cm-3). We employed a cm-scale SERF
magnetometer produced by QuSpin [16] which contains a
3× 3× 3 mm 3 87Rb vapor cell and a single laser for both
optical pumping and probing [17]. This magnetometer is a
compact, self-contained unit with all the necessary optical
components [17] and can be readily operated using a
control software. The magnetometer is surrounded with a
cylindrical ferrite shield with end-caps (18 cm diameter
and 38 cm height), which is inserted into a two-layer open
mu-metal concentric cylindrical shield (26 cm inner
diameter, 29 cm outer diameter, and 69 cm height) to
suppress the ambient magnetic noise. The BGO mass,
69 cm
1.4 cm
SERF
OPM
BGO
G10 rod
motor
Wood support
y
x
z
mu-metal cylindrical shield
FIG. 1. Side view of a schematic of the experimental setup (not
scaled) to probe the exotic spin-dependent interaction V4+5. An
unpolarized BGO mass is placed next to a Rb vapor cell,
located inside a SERF magnetometer head, where the Rb
electron spins are aligned along the x -axis. The mass is rotated
clockwise and counterclockwise around the z -axis.
attached to a rigid plastic rod connected to a stepper motor
(DMX-J-SA-17 provided by Arcus Technology) fixed on
a three-axis translation stage (Thorlabs PT3), was moved
closely to the magnetometer head by the translation stage.
The gap between the nearest vapor cell wall and the
nearest part of the sample was set to ~ 5 mm, which is
smaller than what was considered in Ref. [9] (this is one
of advantages of employing the QuSpin magnetometer).
The centers of the Rb cell and the mass were aligned in
the z -axis. The motor and the translation stage containing
magnetic parts were placed outside the mu-metal shield to
prevent them from negatively affecting the sensitive
magnetic measurements, including systematic signal shifts.
The Rb electron spins inside the vapor cell aligned
along the x -axis in Fig. 1 and the magnetometer are
sensitive to magnetic fields in both directions,
perpendicular to the spin alignment, i.e., y - and z -axis.
The intrinsic field noise levels of the magnetometer in
these directions were measured to be 15 fT/Hz1/2. To
create the relative velocity term, the BGO mass was
rotated around the z -axis by the motor. Specially, to
reduce systematic errors mainly due to trace magnetic
contamination of the BGO mass and the DC offset of the
magnetometer we compared the magnetometer signals
between clockwise and counterclockwise mass rotations,
because the systematic effects will be the same for
opposite rotations while the sign of the interactioninduced signals will be reversed (see Eq. (1)). In this
!
configuration, only the component of (υ × r̂)z remains,
which tilts the Rb electron spins by a small angle similarly
as external magnetic fields [18]. The tilt can be observed
with the laser beam (probing) with very high sensitivity.
For that reason, we focused on the measurement of the
only z direction out of the two available sensitive
directions of the magnetometer because unlike the
y component, the z field measurement is immune to the
rotation of tiny magnetic impurities such as iron particles
possibly buried in the BGO mass.
+ω
-ω
-ω
+ω
160
Magnetic field (pT)
29 cm
26 cm
18 cm
38 cm
Rb cell
162
3-axis
Ferrite shield
translation stage
modulation
158
156
154
152
150
3.45
1
50%
50%
50%
50%
+1
-3
+3
-1
3.46
1.5
3.47
2
3.48
3.49
3.5
3.51
3
4
2.5
3.5
Sample number
3.52
4.5
3.53
5
Time (s)
3.54
5.5
4
x 10
FIG. 2. Time traces of the SERF magnetometer readout showing
two adjacent cycles of the mass rotation reversal. The angular
velocity ω was chosen 1 Hz.
The BGO mass was rotated clockwise for 1 s and then
rotated counterclockwise for 1 s at an angular velocity ω
of 1 Hz. This motion was continuously repeated and the
SERF magnetometer signals were collected for 2 hours.
Figure 2 indicates typical time traces of the magnetometer
readout presenting two adjacent cycles of BGO mass
rotation reversal. During the rotation reversals, the
magnetometer detected spikes with amplitudes on the
order of pT associated with quickly reversing currents in
the motor, which were chosen to be reference points to
chop each cycles. The V4+5 interaction can be extracted by
taking the difference of the magnetic field signals at +ω
and −ω within one cycle. However, there still exist
inevitable DC drifts on the magnetometer readout, mainly
due to the magnetometer electronics. To remove this DC
drifts from the data, we used the “drift-correction
algorithm” which takes the average of the data for each
half cycle in two adjacent cycles of rotation reversal and
applies a [+1 -3 +3 -1] weighting to the averaged data [19].
This algorithm is able to expand the effect of the DC drift
into a polynomial function of the time and cancel it up to
second order. In the assumption that the data average for
each half cycle is B(t) = at + bt 2 ± V4+5 where t is the time
and the first two terms are the DC drifts up to second
order while the last term is the interaction of interest (the
sign will be determined by rotation directions), then the
weighting calculation is given by
ΔB = B(τ / 2) − 3B(τ ) + 3B(3τ / 2) − B(2τ )
= 8V4+5
(4)
where τ is the period of the rotation reversal cycle; as a
result, the DC drifts up to second order is effectively
removed. Note that to obtain the V4+5 , the resulting value
ΔB should be divided by a factor of 8. In addition, in
order to sufficiently eliminate the systematic effects from
the spikes and the following transient regions, only the
last 50 % of the data in each half cycle were used for the
data average in this algorithm.
180
m
ALP
Data
Gaussian fit
160
140
−1
0
10
(eV)
−3
10
−4
10
−5
10
10
3.8 fT sensitivity
0.35 fT sensitivity
χ 2/dof=1.895
120
−2
−5
Log(f4+5)
100
80
60
−10
−15
40
−20
20
0
-1
-0.5
0
V
4+5
0.5
(T)
1
-12
−25 −6
10
−5
10
−4
FIG. 3. Distribution of the V4+5 interaction processed with the
drift-correction algorithm. The data were collected for 2 hours.
Figure 3 shows a histogram of the V4+5 interaction
obtained by the drift-correction algorithm, for data
collected for 2 hours. This distribution was fit to the
Gaussian distribution, giving V4+5 = (2.56 ± 3.81) ×10 −15 T .
This value corresponds to the energy shift of
(4.65 ± 6.92) ×10 −19 eV for Rb atoms, with the Rb
gyromagnetic ratio of 2π × 7.0 ×10 9 Hz/T .
To constrain the coupling strength f4+5 of the
interaction, we performed a Monte Carlo simulation
which averages the interaction potential, operators to the
right of f4+5 in Eq. (1), over the volumes of the BGO
mass and the Rb cell containing the polarized electron
spins [9]:
1) we generated random points inside both the volumes
of the BGO mass and the Rb cell;
2) an interaction range between 10-1 to 10-6 m was
defined to calculate the potential between two randomly
generated points inside both the volumes;
3) the potential was integrated over all the volumes and
normalized to give the average potential for the nucleon
density of BGO mass;
4) the experimental sensitivity to the coupling strength
was derived by dividing the experimental limit of the
energy shift for Rb atoms in the above by the calculated
average potential.
The solid curve in Fig. 4 shows the experimentally set
limit on the coupling strength of the V4+5 interaction
between unpolarized nucleons of the BGO mass and
polarized Rb electron spins in the vapor cell, with the
statistical uncertainty of 3.81×10 −15 T . Our experimental
−3
10
10
−2
10
−1
10
λ (m)
x 10
FIG. 4. Constraints of the coupling strength of the V4+5
interaction between unpolarized nucleons in the BGO mass and
polarized Rb electron spins in the SERF magnetometer as a
function of the interaction range (bottom axis) and the ALP
mass (top axis). The solid curve indicates the experimentally set
limit based on the data while the dashed curve shows the
estimated limit on the assumption of 10 days data collection
time.
method is sensitive to the interaction in the interaction
range below 10-2 m. The statistical uncertainty can be
enhanced by longer data collection: for 10 days collection
time, the statistical uncertainty will be improved to
0.35 ×10 −15 T . The dashed curve in Fig. 4 indicates this
case. Note that the experimental constraints on the
interaction are set for the first time.
In conclusion, we investigated the exotic spin-dependent
V4+5 interaction for polarized electrons using a novel
experimental method based on a SERF magnetometer. We
report the first experimental constraint on the V4+5
interaction of electrons in the interaction range below 1
cm. In near future, we plan to probe the other interactions,
V9+10 and V12+13 , by properly moving the BGO mass next
to the Rb vapor cell. In the case of V12+13 , because it has
never been investigated experimentally, our experimental
results will be quite valuable.
ACKNOWLEDGMENTS
The authors thank Dr. Pulak Nath for help in experiment.
The authors gratefully acknowledge this work was
supported by the U.S. DOE through the LANL/LDRD
program.
__________________________________________________
[1] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38,
1440 (1977).
[2] S. Weinberg, Phys. Rev. Lett. 40, 223 (1978).
[3] F. Wilczek, Phys. Rev. Lett. 40, 279 (1978).
[4] J. Jaeckel and A. Ringwald, Annu. Rev. Nucl. Part.
Sci. 60, 405 (2010).
[5] K. A. Olive et al. (Particle Data Group), Chin. Phys. C
38, 090001 (2014).
[6] G. Bertone, D. Hooper, and J. Silk, Phys. Rep. 405,
279 (2005).
[7] J. E. Moody and F. Wilczek, Phys. Rev. D 30, 130
(1984).
[8] B. A. Dobrescu and I. Mocioiu, J. High Energy Phys.
11, (2006) 005.
[9] P.-H. Chu, Y. J. Kim, and I. Savukov, Phys. Rev. D
94, 036002 (2016).
[10] I. K. Kominis, T. W. Kornack, J. C. Allred, and M.
V. Romalis, Nature 422, 596 (2003).
[11] T. M. Leslie, E. Weisman, R. Khatiwada, and J. C.
Long, Phys. Rev. D 89, 114022 (2014).
[12] G. D. Hammond, C. C. Speake, C. Trenkel, and A. P.
Patón, Phys. Rev. Lett. 98, 081101 (2007).
[13] S. A. Hoedel, F. Fleischer, E. G. Adelberger, and B.
R. Heckel, Phys. Rev. Lett. 106, 041801 (2011).
[14] W. A. Terrano, E. G. Adelberger, J. G. Lee, and B. R.
Heckel, Phys. Rev. Lett. 115, 201801 (2015).
[15] W.-T. Ni, S.-S. Pan, H.-C. Yeh, L.-S. Hou, and J.
Wan, Phys. Rev. Lett. 82, 2439 (1999).
[16] QuSpin Inc. Available at: http://www.quspin.com/.
[17] I. M. Savukov, Y. J. Kim, V. Shah, and M. Boshier,
Meas. Sci. Technol. 28, 035104 (2017).
[18] T. Karaulanov, I. Savukov, and Y. J. Kim, Meas. Sci.
Technol. 27, 055002 (2016).
[19] Y. J. Kim, C.-Y. Liu, S. K. Lamoreaux, G. Visser, B.
Kunkler, A. N. Matlashov, J. C. Long, and T. G.
Reddy, Phys. Rev. D 91, 102004 (2015).