activation time sensitivity to hyperpolarization
(
8
). The S2-S3 linker is also directed toward the
S4/C-linker interface and may play a yet-unknown
role in channel gating.
Summary
The structure of rEag1 reveals a non
–
domain
swapped architecture of the S1 to S6 that is due
to a short five-residue S4-S5 linker. This repre-
sents a divergence from the domain-swapped
architecture of previous voltage-gated ion chan-
nel structures (Fig. 4) (
1
–
4
) and suggests a new
paradigm for voltage-dependent gating for the
Eag family of K
v
channels. On the basis of the
structure of Eag1, we propose a gating mecha-
nism in which S4 enters the cytoplasm in a
down or hyperpolarized state to interact with
and induce a rotation of the C-linker and S6
in a direction that tightens the helical bundle
to close the channel (Fig. 6B). In the up or de-
polarized state of the VS, S4 moves into the
membrane, which allows the C-linker and S6
to rotate in a direction that loosens the helical
bundle and thus relieves the high-energy bend
in S6 to open the channel (Fig. 6B).
Two important consequences result from a
gating mechanism in which the VS interacts
with the cytoplasmic domains to gate the chan-
nel. First, this allows the cytoplasmic domains
to close the channel independent of VS con-
formation. This is observed in the structure of
Eag1, as binding of CaM to the cytoplasmic
domains closes the pore, but the VS is in the
up or depolarized conformation. Second, this
provides an added level of regulation through
the interaction of intracellular domains with
the voltage-dependent gating machinery. In
Eag1, the N terminus of the PAS domain, which
confers sensitivity to hyperpolarization (
8
), is
poised to interact with the S4 and S4-S5 linker
in a closed conformation (Fig. 5D).
REFERENCES AND NOTES
1. S. B. Long, E. B. Campbell, R. Mackinnon,
Science
309
,
897
–
903 (2005).
2. S. B. Long, X. Tao, E. B. Campbell, R. MacKinnon,
Nature
450
,
376
–
382 (2007).
3. J. Payandeh, T. Scheuer, N. Zheng, W. A. Catterall,
Nature
475
,
353
–
358 (2011).
4. J. Wu
et al
.,
Science
350
, aad2395 (2015).
5. S. B. Long, E. B. Campbell, R. Mackinnon,
Science
309
,
903
–
908 (2005).
6. H. R. Guy, S. R. Durell, J. Warmke, R. Drysdale, B. Ganetzky,
Science
254
, 730 (1991).
7. M. Ju, D. Wray,
Biochem. Biophys. Res. Commun.
342
,
1088
–
1097 (2006).
8. H. Terlau, S. H. Heinemann, W. Stühmer, O. Pongs, J. Ludwig,
J. Physiol.
502
, 537
–
543 (1997).
9. J. H. Morais Cabral
et al
.,
Cell
95
, 649
–
655 (1998).
10. F. W. Muskett
et al
.,
J. Biol. Chem.
286
, 6184
–
6191 (2011).
11. J. Wang, M. C. Trudeau, A. M. Zappia, G. A. Robertson,
J. Gen.
Physiol.
112
, 637
–
647 (1998).
12. W. N. Zagotta
et al
.,
Nature
425
, 200
–
205 (2003).
13. T. I. Brelidze, A. E. Carlson, W. N. Zagotta,
J. Biol. Chem.
284
,
27989
–
27997 (2009).
14. T. I. Brelidze, A. E. Carlson, B. Sankaran, W. N. Zagotta,
Nature
481
, 530
–
533 (2012).
15. M. J. Marques-Carvalho
et al
.,
J. Mol. Biol.
423
,34
–
46
(2012).
16. Y. Haitin, A. E. Carlson, W. N. Zagotta,
Nature
501
, 444
–
448
(2013).
17. J. Ludwig
et al
.,
EMBO J.
13
, 4451
–
4458 (1994).
18. R. Schönherr, K. Löber, S. H. Heinemann,
EMBO J.
19
,
3263
–
3271 (2000).
19. U. Ziechner
et al
.,
FEBS J.
273
, 1074
–
1086 (2006).
20. L. S. Mortensen
et al
.,
J. Physiol.
593
, 181
–
196 (2014).
21. P. Bijlenga
et al
.,
J. Physiol.
512
, 317
–
323 (1998).
22. L. A. Pardo, A. Brüggemann, J. Camacho, W. Stühmer,
J. Cell
Biol.
143
, 767
–
775 (1998).
23. L. A. Pardo
et al
.,
EMBO J.
18
, 5540
–
5547 (1999).
24. B. Hemmerlein
et al
.,
Mol. Cancer
5
, 41 (2006).
25. J. R. Agarwal, F. Griesinger, W. Stühmer, L. A. Pardo,
Mol.
Cancer
9
, 18 (2010).
26. F. Mello de Queiroz, G. Suarez-Kurtz, W. Stühmer, L. A. Pardo,
Mol. Cancer
5
, 42 (2006).
27. B. R. Downie
et al
.,
J. Biol. Chem.
283
,36234
–
36240
(2008).
28. D. Gómez-Varela
et al
.,
Cancer Res.
67
, 7343
–
7349 (2007).
29. J. García-Quiroz
et al
.,
PLOS ONE
7
, e45063 (2012).
30.J.Ludwig,D.Owen,O.Pongs,
EMBO J.
16
, 6337
–
6345
(1997).
31. M. Ju, D. Wray,
FEBS Lett.
524
, 204
–
210 (2002).
32. Y. Jiang
et al
.,
Nature
423
,33
–
41 (2003).
33. Single-letter abbreviations for the amino acid residues are as
follows: A, Ala; C, Cys; D, Asp; E, Glu; F, Phe; G, Gly; H, His;
I, Ile; K, Lys; L, Leu; M, Met; N, Asn; P, Pro; Q, Gln; R, Arg;
S, Ser; T, Thr; V, Val; W, Trp; and Y, Tyr.
34. P. H. Hsu
et al
.,
PLOS ONE
7
, e41203 (2012).
35. X. Tao, A. Lee, W. Limapichat, D. A. Dougherty, R. MacKinnon,
Science
328
,67
–
73 (2010).
36. M. Zhang, J. Liu, G. N. Tseng,
J. Gen. Physiol.
124
, 703
–
718
(2004).
37. S. K. Aggarwal, R. MacKinnon,
Neuron
16
,1169
–
1177
(1996).
38. S. A. Seoh, D. Sigg, D. M. Papazian, F. Bezanilla,
Neuron
16
,
1159
–
1167 (1996).
39. Z. Lu, A. M. Klem, Y. Ramu,
Nature
413
, 809
–
813 (2001).
40. D. del Camino, M. Holmgren, Y. Liu, G. Yellen,
Nature
403
,
321
–
325 (2000).
41. É. Lörinczi
et al
.,
Nat Commun
6
, 6672 (2015).
42. R. M. Hardman, P. J. Stansfeld, S. Dalibalta, M. J. Sutcliffe,
J. S. Mitcheson,
J. Biol. Chem.
282
, 31972
–
31981 (2007).
43. G. Barbato, M. Ikura, L. E. Kay, R. W. Pastor, A. Bax,
Biochemistry
31
, 5269
–
5278 (1992).
ACKNOWLEDGMENTS
We thank M. Ebrahim at the Rockefeller University cryo-EM
resource center for help with data collection and J. Chen and
members of the MacKinnon laboratory for helpful discussions. This
work was supported in part by NIH grant GM43949. J.R.W. is a
Damon Runyon Fellow supported by the Damon Runyon Cancer
Research Foundation (DRG-2212-15) and R.M. is an investigator in
the Howard Hughes Medical Institute. The low-pass filtered and
amplitude-modified three-dimensional cryo-EM density maps for
rEag1 have been deposited in the EM Data Bank with accession code
EMD-8215. Atomic coordinates for rEag1 have been deposited in the
Protein Data Bank under accession code 5K7L. J.R.W. performed
functional experiments and expressed, purified, and determined
single-particle cryo-EM structure of rEag1. J.R.W. and R.M. designed
experiments, analyzed and interpreted results, and wrote the
manuscript.
SUPPLEMENTARY MATERIALS
www.sciencemag.org/content/353/6300/664/suppl/DC1
Materials and Methods
Figs. S1 to S9
References (
44
–
65
)
1 April 2016; accepted 22 June 2016
10.1126/science.aaf8070
PHYSICS
Laser spectroscopy of
muonic deuterium
Randolf Pohl,
1,2
*
François Nez,
3
Luis M. P. Fernandes,
4
Fernando D. Amaro,
4
François Biraben,
3
João M. R. Cardoso,
4
Daniel S. Covita,
5
Andreas Dax,
6
Satish Dhawan,
6
Marc Diepold,
1
Adolf Giesen,
7,8
†
Andrea L. Gouvea,
4
Thomas Graf,
7
Theodor W. Hänsch,
1,9
Paul Indelicato,
3
Lucile Julien,
3
Paul Knowles,
10
‡
Franz Kottmann,
11
Eric-Olivier Le Bigot,
3
Yi-Wei Liu,
12
José A. M. Lopes,
4,13
Livia Ludhova,
10
§
Cristina M. B. Monteiro,
4
Françoise Mulhauser,
10,1
|| Tobias Nebel,
1
¶
Paul Rabinowitz,
14
Joaquim M. F. dos Santos,
4
Lukas A. Schaller,
10
Karsten Schuhmann,
11,8,15
Catherine Schwob,
3
David Taqqu,
15
João F. C. A. Veloso,
5
Aldo Antognini,
1,11,15
The CREMA Collaboration
The deuteron is the simplest compound nucleus, composed of one proton and one
neutron. Deuteron properties such as
theroot-mean-squarechargeradius
r
d
and the
polarizability serve as important benchmarks for understanding the nuclear forces and
structure. Muonic deuterium
m
d is the exotic atom formed by a deuteron and a negative
muon
m
–
. We measured three 2S-2P transitions in
m
d and obtain
r
d
=2
:
12562
ð
78
Þ
fm,
which is 2.7 times more accurate but 7.5
s
smaller than the CODATA-2010 value
r
d
=2
:
1424
ð
21
Þ
fm. The
m
dvalueisalso3.5
s
smaller than the
r
d
value from electronic
deuterium spectroscopy. The smaller
r
d
, when combined with the electronic isotope
shift, yields a
“
small
”
proton radius
r
p
, similar to the one from muonic hydrogen,
amplifying the proton radius puzzle.
P
recision spectroscopy of atomic energy
levels can be used to determine properties
of the nucleus (
1
). Deuterium (D), for exam-
ple, is a heavier isotope of hydrogen (H),
with a nucleus, the deuteron (d), composed
of one proton and one neutron (
2
). D was dis-
covered through a tiny shift in the Balmer
spectral lines of D-enriched hydrogen (
3
). This
shift is caused mainly by the mass difference
between the proton and the deuteron. Today, the
nuclear masses are accurately known from cyclo-
tron frequency measurements in a Penning trap
SCIENCE
sciencemag.org
12 AUGUST 2016
•
VOL 353 ISSUE 6300
669
RESEARCH
|
RESEARCH ARTICLES
on August 26, 2016
http://science.sciencemag.org/
Downloaded from
(
1
), and the measured isotope shift of the 1S-2S
transition in H and D (
4
) determines the (squared)
deuteron-proton charge radius difference (
5
)
d
ð
2
Þ
ð
H
;
D
Þ
≡
r
2
d
−
r
2
p
¼
3
:
82007
ð
65
Þ
fm
2
ð
1
Þ
This is because the wave function of atomic S
states is maximal at the origin, where the nu-
cleus resides, and the wave function overlap
with the extended nuclear charge distribution
reduces the atomic binding energy. Equation 1
links measurements of transition frequencies
in H and D. These, together with elastic electron
scattering on protons (
6
) and deuterons (
7
), de-
termine the Rydberg constant
R
∞
,
r
p
and
r
d
in the
CODATA adjustment of the fundamental physical
constants (
1
).
Muonic atoms are a special class of
“
exotic
”
atoms that offer access to nuclear properties with
much higher accuracy. In a muonic atom, the nu-
cleusisorbitedbyonenegativemuon
m
–
,instead
of the usual electrons
e
–
.Themuon
’
slargermass
m
m
¼
207
m
e
results in a muonic Bohr radius that
is smaller than the corresponding electronic Bohr
radius by the ratio of reduced masses
m
red
¼
m
ℓ
m
nuc
=
ð
m
ℓ
þ
m
nuc
Þ
.Here
m
ℓ
is the mass of the
lepton (muon
m
–
or electron
e
−
), and
m
nuc
is the
mass of the nucleus. As the Bohr radius shrinks
proportionally to 1
=
m
red
, the overlap of the muon
’
s
wave function with the nuclear charge distribution
increases as
m
3
red
.For
m
d,
m
red
¼
196
m
e
,andthe
wave function overlap is
ð
m
red
=
m
e
Þ
3
≈
10
7
larger
in
m
dthaninD.AmeasurementoftheLambshift
(2P-2S energy difference) in
m
disthereforeex-
tremely sensitive to the deuteron charge radius
r
d
.
Our recent measurements of the Lamb shift in
muonic hydrogen
m
phaveresultedinavalueof
the proton charge radius
r
p
= 0.84087(39) fm,
which is 10 times more accurate, but 4%, or 7
s
,
smaller (
8
,
9
) than the CODATA-2010 value (
1
),
which is the most recently published CODATA
compilation. This so-called
“
proton radius puzzle
”
has questioned the correctness of various experi-
ments or quantum electrodynamics (QED) calcu-
lations, the value of the Rydberg constant, our
understanding of the proton structure, or the
standard model of particle physics (
10
,
11
).
Here we present measurements of the three
2S-2P transitions in
m
d highlighted in Fig. 1, yield-
ing a precise value of
r
d
. The principle of the ex-
periment is to form
m
d atoms in the metastable 2S
state (
12
) and to measure the 2S-2P transitions by
pulsed laser spectroscopy. Comparison with theory
(
13
)reveals
r
d
. The muonic deuterium data pre-
sented here were acquired in the same measure-
ment period as the muonic hydrogen data in (
8
,
9
).
Independent and reliable calculations of QED
(
14
–
17
) and nuclear structure effects (
18
–
22
)in
m
d, which are required to interpret the experi-
ment, have recently become available and are
summarized in (
13
).
Measurement of the spectral lines of
muonic deuterium
The experiment has been described before (
8
,
9
).
In brief, a 5 × 12 mm
2
beam of low-energy neg-
ative muons
m
−
(3-keV kinetic energy, average
rate 600/s) is stopped in a 20-cm-long target
filled with 1 hPa of D
2
gas at 20°C. A pulsed laser
system (
23
,
24
) is triggered on the detection of a
single arriving muon and provides pulses with
an energy of ~0.25 mJ, tunable around a wave-
length of 6
m
m,andcalibratedagainstwatervapor
absorption lines known within a few megahertz
(
25
). A multipass mirror cavity (
26
)ensuresgood
laser illumination of the muon stop volume. Large-
area avalanche photo diodes (
27
,
28
) detect the
2-keV
K
a
x-rays from the radiative 2P
→
1S tran-
sition that follows the laser-induced 2S
→
2P exci-
tation of
m
d. The laser frequency is changed
every few hours, and the resonances displayed
in Fig. 2 are obtained by plotting the number of
2-keV x-rays (normalized to the number of stopped
muons) detected in time coincidence with the
laser pulse, as a function of laser frequency. On
the peak of the resonance, we recorded up to
10 laser-induced x-rays (
“
events
”
)perhourwith
all data reduction cuts (
9
) applied. The back-
ground level of about 2 events per hour origi-
nates mainly from misidentified muon decay
electrons. About a third of the recorded events
are without laser light, providing the expected
background level shown as horizontal bands in
Fig. 2. The resonances are fitted with a flat back-
ground plus a Lorentzian line shape model that
takes into account varying laser pulse energies
and saturation effects.
The three resonances shown in Fig. 2 are the
m
d
transitions 2
S
F
¼
3
=
2
1
=
2
→
2
P
F
¼
5
=
2
3
=
2
,2
S
F
¼
1
=
2
1
=
2
→
2
P
F
¼
3
=
2
3
=
2
,
and 2
S
F
¼
1
=
2
1
=
2
→
2
P
F
¼
1
=
2
3
=
2
, abbreviated as #1, #2, and
#3, respectively. Their positions and uncertainties
are
n
1
¼
50816
:
27
T
0
:
84
ð
stat
Þ
T
0
:
35
ð
syst
Þ
GHz
ð
2
Þ
n
2
¼
52061
:
2
T
2
:
0
ð
stat
Þ
T
0
:
35
ð
syst
Þ
GHz
ð
3
Þ
n
3
¼
52154
:
1
T
2
:
2
ð
stat
Þ
T
0
:
35
ð
syst
Þ
GHz
ð
4
Þ
The systematic uncertainties of 0.35 GHz arise
from laser frequency fluctuations (
8
)andZeeman
shifts from a conceivable small admixture of circular
polarizedlightandthe5Tmagneticfieldofthe
muon beam line. Line-pulling effects from off-resonant
excitation of neighboring levels are negligible (
29
).
Deuteron charge radius
For the fit of line #1, the Lorentzian width was
fixed to the natural radiative line width of
G
¼
19
:
5 GHz (
8
,
9
), as the freely fitted value
G
¼
13
:
1 GHz is 2
:
6
s
too small. Both fits agreed on
the line center within 0.33 GHz, and the uncer-
tainty quoted in Eq. 2 is the larger one from the
fit with fixed width. The difference
n
3
−
n
2
¼
92
:
9
T
3
:
0 GHz from the fit is in good agreement
(1
:
5
s
) with the theoretical value of 88.045 GHz
(
13
). The amplitude of line #3 is larger than zero
only with a significance of 4
:
5
s
, but it serves to
identify line #2 unambiguously. The alternative
—
namely, that the left peak in Fig. 2 (bottom) is in
fact line #3
—
is disfavored with 6
:
9
s
significance
thanks to the absence of a peak with twice the
amplitude ~90 GHz left of line #2.
Combining the three measured frequencies
and using the theoretical 2P fine structure and
2P
3/2
hyperfine splittings (
13
), we determine the
2P-2S Lamb shift (LS) and 2S hyperfine splitting
(HFS) in
m
d
D
E
exp
LS
¼
202
:
8785
ð
31
Þ
stat
ð
14
Þ
syst
meV
ð
5
Þ
D
E
exp
HFS
¼
6
:
2747
ð
70
Þ
stat
ð
20
Þ
syst
meV
ð
6
Þ
with total experimental uncertainties of 0
:
0034
and 0
:
0073 meV, respectively. The measured 2S
670
12 AUGUST 2016
•
VOL 353 ISSUE 6300
sciencemag.org
SCIENCE
1
Max-Planck-Institut für Quantenoptik, 85748 Garching,
Germany.
2
Johannes Gutenberg-Universität Mainz,
QUANTUM, Institut für Physik & Exzellenzcluster PRISMA,
55099 Mainz, Germany.
3
Laboratoire Kastler Brossel, UPMC-
Sorbonne Universités, CNRS, Ecole Normale Supérieure
–
PSL
Research University, Collège de France, 75005 Paris, France.
4
LIBPhys, Department of Physics, University of Coimbra, 3004-
516 Coimbra, Portugal.
5
I3N, Departamento de Física,
Universidade de Aveiro, 3810-193 Aveiro, Portugal.
6
Physics
Department, Yale University, New Haven, CT 06520-8121, USA.
7
Institut für Strahlwerkzeuge, Universität Stuttgart, 70569
Stuttgart, Germany.
8
Dausinger + Giesen GmbH, Rotebühlstrasse
87, 70178 Stuttgart, Germany.
9
Ludwig-Maximilians-Universität,
Munich, Germany.
10
Département de Physique, Université de
Fribourg, 1700 Fribourg, Switzerland.
11
Institute for Particle
Physics, ETH Zurich, 8093 Zurich, Switzerland.
12
Physics
Department, National Tsing Hua University, Hsinchu 300, Taiwan.
13
Instituto Politécnico de Coimbra, ISEC, 3030
–
199, Portugal.
14
Department of Chemistry, Princeton University, Princeton, NJ
08544-1009, USA.
15
Paul Scherrer Institute, 5232 Villigen
–
PSI,
Switzerland.
*Corresponding author. Email: pohl@uni-mainz.de
†
Present
address: Deutsches Zentrum für Luft- und Raumfahrt e.V. in der
Helmholtz-Gemeinschaft, 70569 Stuttgart, Germany.
‡
Present
address: LogrusData, Vienna, Austria. §Present address:
Forschungzentrum Jülich IKP-2 and RWTH Aachen University,
Germany. ||Present address: International Atomic Energy Agency,
Vienna, Austria. ¶Present address: Honeywell Process Solutions Inc,
500 Brooksbank Avenue, North Vancouver BC V7J 3S4, Canada.
Fig. 1.
n
= 2 levels in muonic deuterium.
The order
of the 2P
3/2
sublevels is changed by the nuclear
quadrupole moment (
13
). The three measured
transitions are indicated.
RESEARCH
|
RESEARCH ARTICLES
on August 26, 2016
http://science.sciencemag.org/
Downloaded from
HFS is in excellent agreement with the theoret-
ical value,
D
E
theo
HFS
¼
6
:
2791
ð
50
Þ
meV (
13
).
The Lamb shift in
m
d is extraordinarily sen-
sitive (
13
) to the root mean square (RMS) deu-
teron charge radius
D
E
theo
LS
¼
228
:
7766
ð
10
Þ
meV
þ
D
E
TPE
LS
ð
7
Þ
−
6
:
1103
ð
3
Þ
r
2
d
meV
=
fm
2
where
D
E
TPE
LS
ð
theo
Þ¼
1
:
7096
ð
200
Þ
meV
ð
8
Þ
is the deuteron polarizability contribution (
13
)from
two-photon exchange (TPE), recently calculated
with good accuracy (
18
–
22
). The charge radius
effectinEq.7contributesasmuchas14%tothe
2P-2S Lamb shift, which explains the excellent
sensitivity of our measurement to
r
d
.Weobtain
r
d
from equating Eqs. 5 and 7, and using Eq. 8,
which yields
r
d
ð
m
d
Þ¼
2
:
12562
ð
13
Þ
exp
ð
77
Þ
theo
fm
ð
9
Þ
where the theory uncertainty is almost ex-
clusively from
D
E
TPE
LS
(Eq. 8). This radius is in
7
:
5
s
disagreement with the CODATA value (
1
),
which is the best estimate of the deuteron ra-
dius obtained from precision spectroscopy of H
and D and electron scattering on protons and
deuterons,
r
d
ð
CODATA
Þ¼
2
:
1424
ð
21
Þ
fm
ð
10
Þ
(seeFig.3).Wearehencefacedwiththefactthat
precision determinations of the Lamb shift in
both
m
p and
m
d, from a total of five measured
resonances, each show a
≥
7
s
discrepancy to the
predictions based on fundamental physical con-
stants from the self-consistent CODATA world
average (
1
), carefully checked QED calculations
(
13
,
30
), and physics within the standard model.
The CODATA deuteron radius
r
d
is tightly
linked to the CODATA proton radius
r
p
,byvirtue
ofEq.1.However,asdetailedin(
31
), we have
deduced a deuteron charge radius considering
spectroscopy data in regular deuterium alone
—
i.e.,
without relying on the value of the proton radius.
This yields a value of
r
d
ð
Dspectroscopy
Þ¼
2
:
1415
ð
45
Þ
fm
ð
11
Þ
in excellent agreement with the CODATA value,
but 3
:
5
s
larger than the value obtained here
from muonic deuterium (see Fig. 3, blue point,
“
D
spectroscopy
”
).
This distinct 3
:
5
s
discrepancy between the
atomic physics determinations of
r
d
from D and
m
d is almost as severe as the 4
:
0
s
atomic physics
discrepancy between the
r
p
values from H spec-
troscopy [see (
1
), table XXXVIII, adjustment 8]
and
m
p(
9
) (see Fig. 4). These two discrepancies
are independent, as explained in (
31
).
The difference between the deuteron radii from
the spectroscopy of electronic and muonic deute-
rium is only 0.017 fm, or 0.8%. Thus, even though the
deuteron charge radius
r
d
ð
e
−
dscatt
:
Þ¼
2
:
130
ð
10
Þ
fm, extracted from elastic electron-deuteron scatter-
ing (
7
), is accurate to 0.5%, it is unfortunately not
accurate enough to distinguish between the values
from
m
dandCODATA.
Proton and deuteron radius puzzle
Many attempts to explain the proton radius dis-
crepancy exist (
10
,
11
). Our muonic deuterium
result provides fresh insight, as the so-called
“
proton radius puzzle
”
is in fact not limited to
the proton; there is a distinct deuteron radius
puzzle. Using
r
d
(CODATA) in Eq. 7 yields a Lamb
shift that is
e
L
S
ð
m
d
Þ¼
0
:
438
ð
59
Þ
meV smaller than
the measured value, Eq. 5, and hence resonance
frequencies that are
∼
104 GHz smaller than ob-
served (Fig. 2). The
e
LS
ð
m
d
Þ
is even somewhat larger
than the proton radius discrepancy
e
LS
ð
m
p
Þ¼
0
:
329
ð
47
Þ
meV between the LS we observed in
m
p and the one calculated with the CODATA value
of
r
p
(
9
).
The ratio of discrepancies in
m
dand
m
p,
e
LS
ð
m
d
Þ
=
e
LS
ð
m
p
Þ¼
1
:
3
ð
2
Þ
is in agreement with the ratio of
the wave-function overlap from the reduced mass
ratio,
½
m
red
ð
m
d
Þ
=
m
red
ð
m
p
Þ
3
¼
1
:
17. Such a scaling
is expected for several beyond
–
standard model
(BSM) physics scenarios (
10
,
11
,
32
–
34
), where
a new force between muons and protons is
SCIENCE
sciencemag.org
12 AUGUST 2016
•
VOL 353 ISSUE 6300
671
Fig. 2. Three measured resonances in muonic deuterium.
The resonances are labeled #1 (
A
), and #2
and #3 (
B
). The signal (
y
axis) is
“
normalized number of events
”
as described in (
8
). Predicted resonance
positions are shown based on Eqs. 7 and 8: The CODATA-2010 deuteron radius (pink, Eq. 10) would
correspond to ~104 GHz lower resonance positions, which is a difference of 7.5
s
.The
“
expected
”
deuteron
radius Eq. 13, (
“
m
p + iso,
”
brown) obtained by combining the proton radius from muonic hydrogen (
9
)
and the electronic isotope shift (
“
iso
”
), Eq. 1, is consistent with the observed resonance positions
within ~2.6
s
.The top and bottom panel
’
s data were recorded in 1 week and 2 days, respectively. As an
example, the three highest points around the peak of resonance #1 contain a total of 260 events,
recorded in 21 hours.
RESEARCH
|
RESEARCH ARTICLES
on August 26, 2016
http://science.sciencemag.org/
Downloaded from