3.1. FSU measurement of H2+/D+ by alternating between large and small cyclotron orbits [39]
With the aim of obtaining an improved result for
md/
mp, the CFR of H
2+ to D
+ was measured at FSU by simultaneously trapping a D
+ and H
2+ and alternating them between the trap center and a 2 mm radius parking orbit [
39]. The D
+ was produced by injecting CD
4 while the H
2+ was produced by simply operating the FEP for a few seconds, which presumably desorbed H
2 from either the FEP itself or the holes in the endcaps. Because H
2 and H
2+ have different internuclear separations, H
2+ produced by ionization of H
2 can be produced in any of the bound vibrational levels up to
v = 19 [
40]. The vibrationally excited levels are all highly metastable, with lifetimes against spontaneous decay, which occurs primarily by electric quadrupole transitions, between 7 days (for more excited levels) to 22 days (for
v = 1), [
41] (see Table I of [
39]). The extra mass-energy due to the rovibrational energy is significant. For instance, the energy difference between
v = 0 and
v = 1 increases the H
2+ mass by approximately 1.4 x 10
-10. In a run of 7 hours, 15 alternate measurements of
fc for each ion were obtained, resulting in a statistical precision per run for the H
2+(
v,
N)/D
+ CFR as low as 4 x 10
-11. Hence different vibrational levels of H
2+ were partially resolved by their difference in
fc. This was the first mass spectroscopy of molecular vibrational energy.
Since the CFR resolution was not sufficient to determine the H
2+ vibrational state in all runs, Stark-quenching was used to increase the rate of rovibrational decay rate to the ground state [
42]. In the large cyclotron orbit, the H
2+ ion experiences a
v x
B motional electric field. This electric field mixes the ground and excited electronic states. This results in a small electric dipole moment which increases the rate of rovibrational decay. For ρ
c = 2 mm and
B = 8.5 T the lifetime of
v = 1 is reduced to 2.13 days, while the lifetimes of higher excited levels are reduced to a few hours. In this way, simply by placing the H
2+ in a 2-mm-radius cyclotron orbit for ~1 week, it was possible to measure the CFR with 7 H
2+ ions that were almost certainly in the vibrational groundstate. However, since the spacing between rotational energy levels was less than the CFR resolution,
e.g., the spacing between
N = 0 and
N = 2 changes the CFR by only 11.5 x 10
-12, and because, even with Stark quenching, the mean lifetimes of the rotational levels are months or years, it was not possible to directly determine the H
2+ rotational state. Hence, estimates of the shift and uncertainty in the final CFR due to the rotational energy of the 7 H
2+ ions were made by assuming a Boltzmann rotational distribution for the parent H
2, and then modeling the rovibrational cascade to the ground vibrational level. The resulting shift agreed with an estimate based on the spread of the measured CFRs. Including a contribution to allow for the possibility that collisions with neutrals might also change the rotational level, the overall correction applied to the H
2+/D
+ CFR due to H
2+ rotational energy was 16(16) x 10
-12.
The largest systematic correction and second largest uncertainty to the CFR was from imbalance in the cyclotron radii used in the PnP measurements between H
2+ and D
+, coupled with special relativity (SR). To obtain this shift and its uncertainty, ρ
c of both ions were systematically varied by varying the length of the cyclotron drive pulse
Td at constant amplitude, and then extrapolating the plot of CFR against
Td2 to zero. Except for possible imbalance in the initial cyclotron energy, this extrapolation gives the CFR corrected for SR. The correction determined with this procedure was 41(7) x10
-12. A third significant systematic resulted from the fact that the
fc measurements of the H
2+ and D
+ were carried out with the trap voltages set so that the H
2+ and D
+ axial frequencies were respectively 80 Hz below and above the detector resonant frequency. This was done so that the change in trap voltage between the PNP measurements on the ions was reduced, hence reducing the shift in the CFR due to the change in ion equilibrium position coupled with magnetic field gradient,
i.e., the “
B1Δ
V” shift. Since the ions were on different sides of the detector resonance, their measured axial frequencies were “pushed” in opposite directions due to the ion-detector interaction, which shifts the
fc’s obtained using the invariance theorem. The required “coil-pushing” correction was 8.2(1.0) x 10
-11. The remaining systematic shifts were the residual
B1Δ
V shift, needing a correction -0.6(0.6) x 10
-12, and that due to the polarizability of the H
2+ [
43,
44], needing a correction 1.1(0.3) x 10
-12. The resulting total systematic correction was 65(18) x 10
-12. With a statistical uncertainty of 6.3 x 10
-12, the final result was
M[D
+]/
M[H
2+(0,0] = 0.999 231 660 004(19).
3.2. MPIK measurement of the atomic mass of the deuteron and HD+ [45]
Using the LIONTRAP apparatus previously used to measure the proton against
12C
6+ [
32], the MPIK collaboration measured the CFR of the deuteron against
12C
6+ and of HD
+ against
12C
4+ [
45]. Compared to the proton measurements, the quadratic magnetic field inhomogeneity was reduced from
B2/
B0 = −7.2(4) x 10
-8 to 6.5(6.5) x 10
-10 mm
-2. The stability of the magnetic field was also improved by stabilization of the pressure of the liquid nitrogen and liquid helium reservoirs and by improved trap alignment. A single axial detector with resonance frequency near 461 kHz was used.
In order to load deuterons and HD+ ions, a surface layer of a deuterated organic compound was printed onto the surface of the carbon-fiber-loaded PEEK target. Unlike the proton measurement, in both cases the ion pairs form a near m/q doublet, so the axial motion was detected using a single tuned circuit. As for the proton measurement, the ions were shuttled into the measurement trap from the adjacent storage traps and measurements of fc for each ion were obtained using the PnA method. Unlike the FSU procedure, where measurements were alternated between the ions, and a CFR measurement derived from a polynomial fit to both sets of fc data for the entire run, in the LIONTRAP procedure a CFR measurement was considered to be the result of single measurements of fc on each ion in the pair. The first ion was chosen at random, so successive measurements could be on the same ion. Each run typically produced 27 CFRs, which were then averaged to give a CFR for the whole run.
Over the D+/12C6+ measurement campaign 41 runs were obtained using 4 ion pairs, each trapped for 1 to 4 months. Analogous to the FSU H2+/D+ measurements, to allow for amplitude dependent shifts due to SR, the cyclotron drive amplitudes Ai of both ions were varied, and an extrapolation made to zero Ai2. Due to the lower magnetic field and the smaller minimum ρc of 10 microns in the PnA, the relativistic shifts were an order of magnitude smaller than for the FSU H2+/D+ measurements. Feedback cooling was also used to reduce Tz to 1.2(5) K, which reduced the initial thermal ρc in the PnA. From the fit to the CFR data with different driven ρc, a D+/12C6+ CFR with statistical uncertainty of 5.4 x10-12 was obtained.
Because the ions in each pair had different mass, the initial thermal cyclotron energy did not cancel in the CFR, even if the ions had the same cyclotron temperature. This required a SR correction of −2.9(1.2) x 10-12 to the D+/12C6+ CFR. Overall, the largest systematic correction was due to the unequal image charges (again resulting from the ions’ different mass), 82.1(4.1) x10-12. However, the largest systematic uncertainty in the CFR overall, 4.7 x10-12, was in the determination of the axial frequency, which was done using the dip technique. Since, on resonance with the detector, the FWHM of the dip due to the 12C6+ was 3 Hz, determination of the ion’s fz to sufficient accuracy required a subdivision of the linewidth by a factor of 500. Due to the ion-detector pushing effect, the measured fz was also sensitive to uncertainty in the detector resonance frequency. Because of the small B2/B0, the correction to the CFR for magnetic field imperfections was only 0.3(0.6) x10-12. Between measurements of fc on each ion, the detector resonance frequency was shifted using a varactor, so the measurements were carried out at the same trap voltage, hence eliminating any B1ΔV shift. The combined systematic uncertainty was 6.5 x 10-12, and the final value for the mass ratio 6M[D+]/M[12C6+] was 1.007 052 737 911 7(85). This is the most precise result for a CFR directly relating to 12C to date.
Similar techniques were used for the HD+/12C4+ measurement. As in work at FSU, the HD+ was assumed to be in its rovibrational groundstate and a correction was made for its polarizability. From one ion pair trapped for 7 weeks 4M[HD+]/M[12C4+] = 1.007 310 263 905(19)(8)(20) (stat)(sys)(total) was obtained.
3.3. FSU measurement of H2+/D+ using simultaneous measurement of cyclotron frequencies in coupled magnetron orbits [46]
In order to eliminate the uncertainty in a CFR measurement due to variation in the magnetic field, in the 1990’s the MIT mass spectrometry group developed a technique in which the modified cyclotron frequencies of a pair of ions were measured simultaneously [
47,
48]. In this method, which is applicable to ion pairs with fractional mass difference in the range 10
-4 < Δ
m/
m < 10
-3, the ions are placed in coupled magnetron orbits, such that the ions orbit the center of the trap, 180° apart, with nearly equal radii of ~0.5 mm. In this configuration, due to the Coulomb interaction between the ions, the magnetron modes of the ions are strongly coupled, while the axial and modified cyclotron modes, though perturbed, remain largely independent. Simultaneous PnP measurements can then be performed on the two ions. In 2002-2003 this technique was applied at MIT to ions with
m/
q near 30, producing 4 CFRs with world record uncertainties of 7 x 10
-12. After a 20-year hiatus, the method was re-developed at FSU and applied to a second measurement of the H
2+/D
+ CFR, the first application to light ions.
More formally, the normal modes of the coupled magnetron motion are a “common-like mode”, which approximates the motion of the center-of-charge of the ions, and a “separation-like mode”, which approximates the vector difference between the ions. The ideal configuration corresponds to minimizing the amplitude of the common-like mode, while setting the amplitude of the separation-like mode,
i.e., the ion-ion separation ρ
s, to its optimal value. As shown in [
31], the CFR can then be derived from a precise value for the difference in the modified cyclotron frequencies of the two ions, Δ
fct =
fct1 –
fct2, combined with less precise values of
fct1 and
fz1 (or
fct2 and
fz2). The fractional uncertainties for
fct1 and
fz1 can be larger than the fractional uncertainty in the CFR by factors of
m/Δ
m and (
fct/
fz)
2 (
m/Δ
m), respectively. (A precise measurement of Δ
fz is not required since the ions follow similar paths in the magnetic and electrostatic fields. Hence, effectively, Δ
fz is determined by Δ
fct.) Applying the PNP technique simultaneously to both ions, the CFR measurement is essentially reduced to a precise measurement of the phase difference Δφ = 2πΔ
fctTevol, as determined from the phases of the simultaneous axial ring-down signals. Importantly, the sensitivity to shifts to
fz that would otherwise affect the CFR is greatly relaxed. Implementing this method required re-developing the important tool of “phase-locked driven axial motion” [
31]. This allowed the continuous measurement of an ion’s
fz in real time, and was essential for monitoring the amplitudes of the common and separation modes of the ion pair, and for cooling the common-mode motion.
A run began with a (typically) 15-minute period of “phase-lock” cooling of any common-mode motion that had been produced in the previous run. The actual CFR measurement then consisted of cycles of simultaneous PnPs on the two ions, with a longest Tevol of 10.1 s, interleaved with PnPs with Tevol of 0.1, 0.3, 1.1 and 3.3 s, which were needed for phase unwrapping the individual fct. Throughout the run, phase coherence was maintained between all synthesizers used for the PnPs. Hence, the phases for different Tevol could be averaged over the whole run, and phase unwrapping applied to the averaged phases. After trials with different ion-ion separations it was found that ρs = 0.8 mm was optimum. This gave the best compromise between stability of the coupled magnetron motion, which improved with reduced ρs due to increased ion-ion coupling, and the need to minimize ion-ion induced axial anharmonicity, which could only be partially compensated by applying C4.
The improvement in precision using the simultaneous technique was less dramatic than at MIT with
m/
q = 30. This was partly because the ambient magnetic field at FSU was more stable than at MIT, but also because, at low
m/
q, noise on
fct due to fluctuations in ρ
c combined with SR was comparable in magnitude to noise due to magnetic field fluctuations. This SR noise on
fct is given by σ(
fct)/
fct = (2π
fct/
c)
2σ(ρ
c2)/2, where σ(ρ
c2) is the rms fluctuation in ρ
c2 from PnP to PnP. σ(ρ
c2) originates from the cyclotron motion at the start of the PnPs, which varies randomly from PnP to PnP, and which combines by phasor addition with the driven cyclotron motion, with the result σ(ρ
c2) = 2
1/2ρ
cthρ
cdrive, where ρ
cth is the rms value of the initial cyclotron radius, and ρ
cdrive is the radius produced by the drive. ρ
cth is given by ρ
cth = (2
kBTc/
m)
1/2/(2π
fct), where
Tc is the ion’s effective cyclotron temperature resulting from cyclotron-to-axial coupling. In the ideal case,
Tc = (
fc/
fz)
Tz, where
Tz is the ion’s axial temperature. Hence, ρ
cth, and the minimum ρ
cdrive (~5ρ
cth) for adequate phase initialization in the PnP, are essentially independent of the ion’s mass. So, overall, this relativistic noise varies as
fct2 and so is a more serious issue for light ions. In order to reduce this relativistic noise,
Tz was reduced by a factor of 2 by applying electronic feedback to the axial motion of each ion, using the scheme described in [
49]. This was done with
fz shifted to resonance with the detector by changing the trap voltage. However, even with feedback the overall gain in statistical precision in a 6-hour run was only about a factor of two compared to a run with the alternating technique.
With ρ
s = 0.8 mm both ions were outside the axial line of sight to room temperature. Further, the ions were not in large cyclotron orbits during the measurement. Hence, the average ion lifetime against collision with neutrals was considerably longer than with the alternating technique and excited vibrational levels did not undergo Stark quenching. Combined with the factor-of-two improved resolution, this enabled the tracking of the rovibrational decay of 3 different H
2+ ions to the vibrational ground state. The rovibrational decays manifested as discrete jumps in the H
2+/D
+ CFR between plateaus corresponding to a given rovibrational state. In one case an H
2+ was tracked from
v = 9 to
v = 0 over a period of more than two months. Taking account of the electric-quadrupole selection rule for H
2+ rovibrational decay, Δ
N = 0, ±2, it was possible to fit the plateaus in CFR to calculated shifts using the theoretical rovibrational energies [
50], and so assign certain plateaus to unique rotational levels on a probabilistic basis. Hence, to the extent that the assignment was correct, the uncertainty due to rotational energy was eliminated. Moreover, because the fit averaged over more than 300 runs, a very small statistical uncertainty of 2 x 10
-12 was obtained.
As with the alternating technique, in order to correct for the systematic shift due to SR and imbalance in ρc, CFR measurements were made with a range of ρc in the PnPs. This resulted in a correction of 29.5(1.4) x 10-12. In the simultaneous method the trap voltage was set so that fz of the H2+ and D+ ions were symmetrically below and above the detector resonance frequency. To cool the axial motion before the PnP, each ion was shifted to resonance by changing the trap voltage. This process is necessarily asymmetric, and, because of possible noise spikes or other asymmetries in the detector noise, there was concern that Tz and, so Tc, at the start of the PnPs could be different between the ions, leading to a systematic SR shift to the CFR. In order to estimate a possible difference in Tc, use was made of the fact that this would also result in a difference in the rms fluctuations in the individual cyclotron frequencies fct1 and fct2, as discussed above. These frequency fluctuations were determined from the Allan deviation of the long-Tevol phases of the individual ions. From this a correction of 2.9(2.9) x 10-12 to the CFR was derived, which was the largest source of systematic uncertainty.
Since fz did not need to be known precisely, the detector-pushing effect, and in fact all effects that shift the individual fz’s, including ion-ion interaction, had negligible effect on the CFR. Because of the symmetry between the ions, the ion-ion interaction effects on Δfct and hence the CFR were <10-13 and so negligible. There was no B1ΔV shift. However, because the ions did not have identical mass, the magnetron radii of the two ions in the coupled magnetron motion were not identical. The resulting correction for trap imperfections and rms magnetron radius difference was −1.1(0.2) x 10-12. Finally, allowing for a possible difference in the rms axial amplitudes of the ions due to a difference in Tz during the cyclotron phase evolution, which produces a shift by interacting with B2, there was a correction of 0.5(0.5) x 10-12. The final result for the mass ratio M[D+]/M[H2+(0,0)] was 0.999 231 660 003 0(21)(37)(43), (stat)(sys)(total). This result is in excellent agreement with the alternating method. It is also the most precise mass ratio to date. A caveat is that the rotational state identification was probabilistic. If one of the two possible but less probable assignments is chosen, the mass ratio shifts down by 2.7 or 3.6 sigma.
3.4. LIONTRAP mass of 4He [51]
Following the measurement of the atomic mass of the deuteron, the LIONTRAP apparatus was used by the MPIK collaboration to measure the atomic mass of
4He [
51]. Since in LIONTRAP the trap is completely enclosed, a He source was developed that loads gas from a reservoir inside the trap chamber into the creation trap in front of the FEP. Although it was initially planned to measure
3He to help resolve the light-ion puzzle, due to a technical issue only
4He could be loaded, after which it was decided to measure the CFR
4He
2+/
12C
6+. This was serendipitous. Their result, using methods that have by now been well validated, was in more than 6-sigma disagreement with the previously accepted result, published by the UW group nearly 20 years earlier.
As in the measurement of the D+/12C6+ ratio, the ions were trapped in different traps in the electrode stack, each shuttled into the measurement trap for the fc measurements. The PnA method was used and a CFR measurement consisted of an fc measurement on each ion, with the first being chosen randomly, the complete CFR measurement taking 3800 s. A single axial detector with resonant frequency near 468 kHz was used, again shifted in frequency using a varactor to match the respective fzi of the two ions at the same trap voltage, eliminating the B1ΔV shift. The correction for amplitude-dependent shifts due to SR and trap imperfections was obtained using ρc from 10 to 80 μm and fitting the CFR versus the squares of the drive strengths of the respective ions. With a total data set of 482 cycles this gave a CFR with statistical uncertainty of 9 x10-12. By using feedback to reduce the ion’s axial temperature to 1.7(3) K, the correction and uncertainty due to the ions’ cyclotron energy before the cyclotron drive pulse was only -1.8(0.3) x 10-12. Again similar to the D+/12C6+ measurement, the largest systematic correction, at -65.8(3.3) x 10-12, was due to image charge effects, while the largest contribution to the systematic uncertainty, 7.1 x 10-12, was from the determination of fz, by fitting the dip in the detector noise signal. Additional corrections due to magnetic field inhomogeneity and electrostatic anharmonicity were essentially negligible. The final result for 3M[4He2+]/M[12C6+] was 1.000 650 921 192 8(90)(78)(119) (stat, sys, total).
3.5. FSU measurement of HD+/3He+, HD+/T+ and T+/3He+ for the beta-decay Q-value of tritium and improved masses of T and 3He [52]
Previously, in 2014-15, the FSU group measured the HD
+/
3He
+ and HD
+/T
+ CFRs and from the double ratio obtained a
Q-value for tritium beta-decay, with uncertainty of 0.07 eV [
53]. This was the first measurement on light ions by the FSU group and also the start of the so-called light ion (or
3He) mass puzzle. This was the 4-sigma discrepancy between
M[HD
+]/
M[
3He
+] derived from the atomic masses of
p,
d and
h individually referenced to
12C, and the same mass ratio as measured by FSU. Expressed as a mass difference,
M[
p] +
M[
d] –
M[
h] obtained using
M[
d] and
M[
h] from the University of Washington (UW) group [
54], and
M[
p] from CODATA-2010 [
55] (itself mainly derived from earlier measurements by UW), was greater than that obtained from the FSU HD
+/
3He
+ mass ratio [
53] by 0.79(18) nu. This was the first indication that some previously accepted values of light ion masses, obtained with single ion Penning trap techniques, might have significantly underestimated uncertainties.
Two years later, using a rebuilt set-up with an improved detector and a more homogeneous magnetic field, and an outer ion radius increased from 1.07 to 2 mm, the FSU group re-measured the HD
+/
3He
+ ratio, both directly [
2] and also using H
3+ as an intermediary [
56]. This confirmed the original HD
+/
3He
+ CFR of [
53] and reduced its uncertainty. (The measurements against H
3+ were complicated by the mass shift due to highly-excited, metastable rotational states of H
3+, and so only produced a lower limit for 2
M[
p]
– M[
d].) The discrepancy in
M[
p] +
M[
d] –
M[
h] was also partly resolved by the MPIK collaboration’s measurements of
M[
p] [
32] and
M[
d] [
45] (see 3.1 above). If these replaced the CODATA-2010 [
55] and UW [
54] values,
M[
p] +
M[
d] –
M[
h] differed from the value from the HD
+/
3He
+ ratio of [
53] by 0.35(15) nu, and from that of [
2] by 0.26(9) nu. Nevertheless, given the remaining discrepancies, and the importance of the tritium
Q-value, the FSU group decided to repeat the measurements with tritium using the improved apparatus.
Although the simultaneous method was considered, the measurements used the alternating technique. In the case of HD
+/
3He
+ and HD
+/T
+ the ions in the pairs are separated in mass by a fraction of 2 x 10
-3, which resulted in an axial frequency difference of 670 Hz. Consequently, if the trap voltage was set so the ions were positioned symmetrically above and below the detector resonance as required for the simultaneous method, the ions would be each separated by 16 FWHM from the center of the coil resonance, significantly reducing the signal-to-noise for detection of the axial motion. Neither was the simultaneous method applicable to directly measuring T
+/
3He
+ since the fractional mass difference is only 6.6 x10
-6. At the optimum ion-ion separation of 0.8 mm, this would have caused the axial motions of the two ions to be strongly coupled, preventing application of the PnP method. However, with the alternating method, and with the outer ion in a 2 mm radius cyclotron orbit, the separation in
fz between the inner and outer ion was increased to close to 20 Hz due to the residual
C6 and
B2. This enabled PnPs with negligible interference from ion-ion coupling. Compared to using HD
+ as an intermediary, the direct measurement of the T
+/
3He
+ CFR reduced the time required to achieve a given statistical uncertainty by a factor of 4. The improved detector compared to [
53] enabled the use of a smaller ρ
c, and in combination with a x30 reduction in
B2, to −3.7(7) x 10
-9 mm
-2, allowed ρ
c to be varied to quantify the systematic due to special relativity and cyclotron radius imbalance. Additionally, with a parking radius of 2 mm the effects of ion-ion interaction on the CFR were negligible.
Similar to the alternating D+/H2+ measurement, a run typically consisted of 7 hours of data taking with 15 interchanges, and yielded a statistical uncertainty of 4 x 10-11 for the best runs. However, this statistical uncertainty was degraded for approximately 50% of the runs due to rapid changes in the ambient magnetic field due to the operation of a magnetic spectrograph in a nearby laboratory, and also due to electromagnetic interference on the detector signal. The final results were based on 84 runs of HD+/3He+, 74 of HD+/T+ and 79 of T+/3He+, with additional runs for calibrating the cyclotron drives and investigating systematic errors. From independent fits to the HD+/3He+, HD+/T+ and T+/3He+ CFRs vs Td2, non-correlated statistical uncertainties of 11.4, 13.2 and 8.6 x10-12, respectively, were obtained.
In contrast to the above H
2+/D
+ measurements, the PnPs were done at the same
fz. Hence, it could be assumed that the thermal cyclotron energies were balanced, eliminating any residual relativistic shift after the extrapolation to zero
Td2. The detector-pushing effect on
fz between the ions was also balanced, and so had negligible effect on the CFR. To calibrate the
B1Δ
V shift, measurements were carried out with a T
+/H
2+ pair, with the H
2+ having been previously stored in a 2 mm cyclotron radius orbit for more than 3 days, so that it could be assumed to be in the
v = 0 or
v = 1 vibrational state. Making use of an adequately precise prediction for the T
+/H
2+ CFR, a systematic correction of −1.5(4) x10
-12 to be applied to the HD
+/
3He
+ and HD
+/T
+ CFRs was determined. A correction of 94.3(1) x 10
-12 was also applied to these two CFRs to allow for the polarizability of HD
+ [
43,
44]. All other systematics, including those due to ion-ion interaction were at the level of 10
-13 or less. After applying the systematic corrections and uncertainties, a least-squares adjustment (LSA) to the three ratios resulted in
M[
3He
+]/
M[HD
+] = 0.998 048 085 131 8(92),
M[T
+]/
M[HD
+] = 0.998 054 687 290 2(97), and
M[
3He
+]/
M[T
+] = 0.999 993 384 973 2(77), with correlation coefficients (labeling the three ratios as 1,2,3)
r12 = 0.67,
r13 = 0.36, and
r23 = −0.46.