Abstract
Experimental techniques to manipulate cold molecules have seen great development in recent years. The precision measurements of cold molecules are expected to give insights into fundamental physics. Here we use a rovibrationally pure sample of ultracold KRb molecules to improve the measurement on the stability of electrontoproton mass ratio \(\left( {\mu = \frac{{m_{\mathrm{e}}}}{{M_{\mathrm{p}}}}} \right)\). The measurement is based upon a large sensitivity coefficient of the molecular spectroscopy, which utilizes a transition between a nearly degenerate pair of vibrational levels each associated with a different electronic potential. Observed limit on temporal variation of μ is \(\frac{1}{\mu }\frac{{d\mu }}{{dt}} = (0.30 \pm 1.0) \times 10^{  14} \, {\mathrm{year}}^{  1}\), which is better by a factor of five compared with the most stringent laboratory molecular limits to date. Further improvements should be straightforward, because our measurement was only limited by statistical errors.
Introduction
Cold molecules are becoming a popular tool for precision measurements. The search for the electron electric dipole moment (EDM) is a good example^{1,2,3}. Historically, hot atoms in a vapor cell or in an atomic beam were gradually replaced with cold molecular sources and trapped molecular ions. Cold molecules have a number of advantages for increasing the precision, like longer interaction time, smaller motional decoherence, and it is easier to make use of its rich internal degrees of freedom^{4}. In the case of EDM experiments, molecules (and molecular ions) were placed in Ω doublet levels, where large internal electric field and systematic error rejection are available^{5}. Similarly, one can expect to improve the measurement of the variation of fundamental constants by selecting molecules with a level structure that enhances the sensitivity.
Thus far, various physical systems and experimental techniques have been used to measure the stability of the electrontoproton mass ratio μ = m_{e}/M_{p}: each of them provide a crucial independent check^{6, 7}. Molecular spectroscopy is often called modelfree, because it directly reflects the stability of the inertial mass of the nucleus. Shelkovnikov et al. studied the absorption lines of SF_{6} and obtained a limit on the current variation of μ as \(\frac{1}{\mu }\frac{{d\mu }}{{dt}} = (  3.8 \pm 5.6) \times 10^{  14}\, {\mathrm{year}}^{  1}\)^{8}.
In the mean time, impressive results of \(\frac{1}{\mu }\frac{{d\mu }}{{dt}} = 0.2(1.1) \times 10^{  16}\,{\mathrm{year}}^{  1}\)^{9} and \(\frac{1}{\mu }\frac{{d\mu }}{{dt}} =  0.5(1.6) \times 10^{  16}\,{\mathrm{year}}^{  1}\)^{10} were obtained by using atomic clocks. These measurements are essentially measuring the variation of the electrontoproton magnetic moment ratio.
Here we report a measurement on the variation of μ using ultracold KRb molecules. Diatomic alkali molecules like KRb have two lowlying, overlapping potentials: the deep X^{1}Σ^{+} ground state and the shallower a^{3}Σ^{+} state. It has been pointed out that using a transition between a nearly degenerate pair of vibrational levels each associated with X^{1}Σ^{+} and a^{3}Σ^{+} potentials, one can realize an enhanced sensitivity to the variation of μ^{11}. We realize such an enhancement (on the order of ~10^{4}) by using ultracold KRb molecules produced by photoassociation in a magnetooptical trap followed by stimulated Raman adiabatic passage^{12}. We obtain a better result than previous measurements using molecules^{8}. That is solely limited by counting statistics and demonstrates the power of the method.
Results
Sensitivity to the variation of μ
To quantify the sensitivity of the transition to the variation of fundamental constants, we define the sensitivity W_{μ} and the sensitivity coefficient K_{μ} as:
W_{μ} and K_{μ} determine the absolute and relative accuracy needed for frequency standards used in the experiment. In general, K_{μ}s for vibrational and rotational transitions of diatomic molecules are 1/2 and 1, respectively. Large sensitivity coefficients can be realized using transitions between nearly degenerate states with different symmetries. Such transitions were successfully used in previously mentioned astronomical observations, where mixed torsion–rotation transitions with K_{μ} = −32.8 were investigated^{13}. A similar idea was used to improve the limit on the variation of the finestructure constant^{14}.
Target transition
In this study, we used such transitions in ultracold KRb molecules. As pointed out in ref. ^{11}, closely lying states in X^{1}Σ^{+} and a^{3}Σ^{+} potentials give rise to a large K_{μ}. Figure 1 shows the actual energy curves calculated for ^{41}K^{87}Rb^{15}. From these potential energy curves we can calculate the sensitivity for each vibrational state with the accuracy only being limited by the uncertainty of the potential energy curve (“Methods” section). Although the vibrational states with the binding energies about 1000 cm^{−1} in X^{1}Σ^{+} have the highest sensitivities of W_{μ} = 35 THz, their energies are far from that of the vibrational states in a^{3}Σ^{+}, whose largest binding energy is about 240 cm^{−1}.
We found that one of the best combinations is that of the v = 86 state in the X^{1}Σ^{+} potential and the v = 16 state in the a^{3}Σ^{+} potential because the transition frequency is in the microwave range, the sensitivity coefficient is large (~10^{4}), and the strength of the transition is not too small owing to the finite singlettriplet mixing (~2%). Though its sensitivity is not the highest, it is as high as about one fourth of that. Moreover, it is another advantage that the experiments with microwaves are easier than the experiments with frequencies in the optical domain.
Understanding the details of the closely lying hyperfine manifolds is mandatory for estimating the strength of each microwave transition. In the following discussion, we express the hyperfine states as S, F_{1}, F, m_{F}〉 in short, where the detail of labeling is in the “Method” section. Note that S = 0 (S = 1) signifies the v = 86 state in X^{1}Σ^{+} (v = 16 state in a^{3}Σ^{+}).
We obtained the hyperfine spectra by sweeping the frequency of the stimulated Raman adiabatic passage (STIRAP) laser (Fig. 2). The results are summarized in the Supplementary Note 1 and Supplementary Table 1. We identified the most suitable transition to test the stability of μ as i〉 − f〉 (Fig. 1) where
The i〉 – f〉 transition fulfills the required conditions: the sensitivity to the variation of μ is large, the magnetic sublevels are resolved by a small magnetic field, and there is no firstorder Zeeman shift. The transition frequency is ν_{i−f} = 634.96 MHz, and the sensitivity of the transition to the variation of μ is W_{i−f} = ∂ν_{i−f}/∂ (ln μ) = −9.45(4) THz (“Methods” section). The resulting sensitivity coefficient is:
Spectroscopy
Figure 3 shows a typical spectrum of the i〉 − f〉 transition obtained by 6 h of data collection. A Gaussian fit to the data provided the fullwidthofhalfmaximum (FWHM) of ~50 Hz, which is consistent with the ideal spectrum obtained by the microwave πpulse of 16 ms duration. Longer pulse durations degraded the signal because ballistically expanding clouds of molecules with mean velocity ~130 mm s^{−1} started to leave the area of the ionization beam whose radius was ~2 mm. To access only the m_{F} = 0 state in the S = 1, F_{1} = 1/2, F = 1〉 manifold a magnetic field of ~1.5 μT was applied. The measurements were repeated every 120 ms. The statistical uncertainty of the central frequency was ~93 mHz.
Statistical and systematic errors
We paid special attention to the magnetic field during the measurement as it drifts because of the magnetization of the metal chamber and fluctuates due to other environmental factors. The i〉 − f〉 transition is immune to the firstorder Zeeman shift; however, it is affected by the secondorder Zeeman shift (~0.82 Hz μT^{−2}). To evaluate this shift, we simultaneously measured the transition frequency to the hyperfine state f_{1}〉 ≡ S = 1, F_{1} = 1/2, F = 1, m_{F} = 1〉 whose energy is sensitive to the firstorder Zeeman shift \((\sim \! 3.3{\kern 1pt} {\mathrm{kHz}}\,\mu {\mathrm{T}}^{  1})\). The observations for i〉 − f〉 and i〉 − f_{1}〉 alternated every 2.4 s, and the frequency of the i〉 − f_{1}〉 transition was used to subtract the secondorder Zeeman shift from the measured i〉 − f〉 transition frequencies. The statistical and systematic errors caused by the correction step were evaluated as 15 and 1 mHz, respectively, where the systematic error is caused by the error of the Zeeman coefficients.
The blackbody radiation (BBR) shift can fluctuate owing to the temperature fluctuation of the vacuum chamber (26.8 ± 2 °C). The coefficient was calculated using ab initio calculations for the transition dipole moments^{16}. The BBR shifts for i〉 and f〉 states are on the order of 1 Hz, but the effect on the transition frequency is far smaller since two shifts almost cancel out. Because of the difficulty of the calculation, the uncertainty of the BBR shift on the transition frequency is larger than the central value of the estimation, where the BBR shifts induced by the coupling with other vibrational states in X^{1}Σ^{+} and a^{3}Σ^{+} are much smaller than that induced by the coupling with other electronic states. From both of the temperature fluctuation and the calculated BBR coefficients with the uncertainty (0.08(20) Hz at 300 K), we estimated the fluctuation of the blackbody shift during the measurements as <10 mHz.
There are several systematic effects with fluctuations <1 mHz; hence, we can neglect them in the uncertainty budget table (Table 1). The list includes the Stark shift, the density shift, and the reference clock. The Stark shift is caused by a static electric field in the vacuum chamber, which effectively shields environmental electric field fluctuations. The main source of the electric field in the detection region is the microchannel plate (MCP). We canceled this field during the measurements by applying a compensation field from electrodes placed inside the chamber, where we measured the electric field by observing the microwave transition between the rotational states of X^{1}Σ^{+}, v = 0, N = 0〉 and N = 1〉, whose sensitivity to the electric field (~12.4 Hz cm^{2} V^{−2}) is one order higher than that of the i〉 − f〉 transition (~1.08 Hz cm^{2} V^{−2}). However, a finite residual electric field of 0.19(6) V cm^{−1} persists and it causes a dcStark shift of ~40 mHz, where the stark shift coefficients are estimated by using the ab initio calculation for the dipole moment^{17}. This issue, however, is not important because the fluctuation of the Stark shift, which comes from the fluctuation of the MCP voltage, should be <1 mHz. The fluctuation of the frequency clock is also negligible because we use a commercial Rbfrequency standard (RSR FS725) locked to the Global Positioning System signal \((\delta f/f\sim 10^{  12})\). The density shift is also <1 mHz because we blast off the atomic clouds by shining the resonant lasers before the microwave application. Note that though the photoassociation recoil induces a finite mean velocity to the molecular cloud and that causes a frequency shift of about 10 mHz at most by the Doppler effect for microwave, the shift does not affect the measurement of the temporal variation of μ, because all measurements were performed in the identical configuration and the shift is also identical for all measurements.
Stability of μ
We accumulated data intermittently for about 16 months and present the result in Fig. 4. From the linear fit to the data, we obtain the temporal variation of f as:
and the temporal variation of μ as:
This is the most accurate molecular test for the temporal variation of μ to date. In a precise discussion, we measured the frequency ratio of f/f_{Cs}, where the f_{Cs} is the clock transition frequency of Cs atom. However, since the precision of our measurement is about 10^{−10}, which is much lower than that of the other stability tests using f_{Cs}, we ignored the shift of f_{Cs} for simplicity.
Discussion
We have derived a molecular limit on the temporal variation of μ by observing the microwave transition of photoassociated KRb molecules. The result is noteworthy for the following reasons. First, because of the large sensitivity coefficient of about 10^{4} the variation of μ was evaluated with four orders of magnitude higher accuracy than that of the frequency in our measurement. This can be a great benefit to the improvement of the precision spectroscopy because we can greatly reduce the necessity of the precision of the reference clock. Second, the sensitivities to the variation of μ of the molecular levels can be calculated from the experimentally determined molecular energy curves, thus they can be determined almost model independently. Improving the measurements using the molecular levels is important because it is independent of the measurement related to the atomic nuclear magnetic moments. Third, since ultracold molecules synthesized from ultracold atoms are used, the application to the next generation of spectroscopic methods is expected. We started to construct a machine to incorporate the idea in ref. ^{18} to cool alkali–metal dimers to submicrokelvin temperatures by narrowline laser cooling. We can make the best use of submicrokelvin molecules by building a molecular fountain setup^{11}, where narrowing the spectral linewidth to about 1 Hz is expected. Another solution to achieve a narrower linewidth is holding the molecules in an optical lattice with a magic wavelength^{19}.
By means of these techniques, spectral linewidths of about 1 Hz will be achieved. In addition, by making the efforts for a larger number of molecules and more data accumulations, the improvement of the precision by twothree orders of magnitude is expected in future experiments. Though the uncertainty of the calculation of BBR shift gives the biggest systematic error in our results, the BBR shift coefficients can be experimentally determined with good accuracy by the measurements in the longwavelength infrared fields like CO_{2} lasers, as performed for the ion clock transition^{20}.
Methods
Ultracold molecular spectroscopy
Our experimental procedure for producing KRb molecules is described in detail in^{12}(Fig. 2). Briefly, KRb molecules in the X^{1}Σ^{+}, v = 91, N = 0〉 state were produced by photoassociation in a magnetooptical trap (MOT) of ^{87}Rb and ^{41}K atoms. The temperature of the molecules was about 140 μK, which is close to the temperature of atoms. Then the MOT field is turned off and both the Rb and K atoms were blasted off by shining the resonant lasers in order to suppress the energy shift caused by the molecule–atom collisions. After that, molecules were transferred to the initial state (i.e., v = 86 state in the X^{1}Σ^{+} potential) by STIRAP. We used an intermediate state for the Raman transition of (3)^{1}Σ^{+}, v = 41, J = 1〉 with Raman laser wavelengths of 875 and 872 nm. After irradiating the molecules with a microwave pulse, molecules in the initial state were optically pumped to different levels by applying the 872 nm laser. Then the molecules in the final state were selectively observed by resonance enhanced multiphoton ionization by applying a pulsed laser whose wavenumber was 16218.1 cm^{−1}.
Labeling of the hyperfine states
We express the hyperfine states based on Hund’s case b_{βS}^{21}. The hyperfine states were labeled S, I_{Rb}, F_{1}, I_{K}, F_{2}, N, F, m_{F}〉, where S is the total electron spin, I_{Rb}(I_{K}) is the nuclear spin of rubidium (potassium), N is the rotational angular momentum, and F(m_{F}) is (the z component of) the total angular momentum of the molecule. F_{1} and F_{2} are defined as F_{1} ≡ S + I_{Rb} and F_{2} ≡ F_{1} + I_{K}, and F satisfies the relation F ≡ F_{2} + N. Note that I_{Rb} = I_{K} = 3/2. We only discuss N = 0 states in this letter; thus, F = F_{2}. In the “Results” section, we express the hyperfine states as S, F_{1}, F, m_{F}〉 in short.
Uncertainty of the sensitivity
We calculated the sensitivities of i〉 and f〉 by using the experimentally determined potential energy curves reported in^{15}. Their uncertainties were estimated from the uncertainties of the vibrational level intervals^{11}, which were estimated from the difference between the vibrational energies predicted by^{15} and our experimental observations.
Data availability
The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.
References
 1.
The ACME Collaboration. Improved limit on the electric dipole moment of the electron. Nature 562, 355–360 (2018).
 2.
Cairncross, W. B. et al. Precision measurement of the electron’s electric dipole moment using trapped molecular ions. Phys. Rev. Lett. 119, 153001 (2017).
 3.
Hudson, J. J. et al. Improved measurement of the shape of the electron. Nature 473, 493 (2011).
 4.
Carr, L. D., DeMille, D., Krems, R. V. & Ye, J. Cold and ultracold molecules: science, technology and applications. New J. Phys. 11, 055049 (2009).
 5.
Eckel, S., Hamilton, P., Kirilov, E., Smith, H. W. & DeMille, D. Search for the electron electric dipole moment using Ωdoublet levels in PbO. Phys. Rev. A 87, 052130 (2013).
 6.
DeMille, D., Doyle, J. M. & Sushkov, A. O. Probing the frontiers of particle physics with tabletopscale experiments. Science 357, 990–994 (2017).
 7.
Uzan, J.P. The fundamental constants and their variation: observational and theoretical status. Rev. Mod. Phys. 75, 403–455 (2003).
 8.
Shelkovnikov, A., Butcher, R. J., Chardonnet, C. & AmyKlein, A. Stability of the protontoelectron mass ratio. Phys. Rev. Lett. 100, 150801 (2008).
 9.
Godun, R. M. et al. Frequency ratio of two optical clock transitions in ^{171}Yb ^{+} and constraints on the Time variation of fundamental constants. Phys. Rev. Lett. 113, 210801 (2014).
 10.
Huntemann, N. et al. Improved limit on a temporal variation of m _{p}/m _{e} from comparisons of Yb ^{+} and Cs atomic clocks. Phys. Rev. Lett. 113, 210802 (2014).
 11.
DeMille, D. et al. Enhanced sensitivity to variation of m _{e}/m _{p} in molecular spectra. Phys. Rev. Lett. 100, 043202 (2008).
 12.
Aikawa, K. et al. Coherent transfer of photoassociated molecules into the rovibrational ground state. Phys. Rev. Lett. 105, 203001 (2010).
 13.
Bagdonaite, J. et al. A stringent limit on a drifting protontoelectron mass ratio from alcohol in the early universe. Science 339, 46–48 (2013).
 14.
Leefer, N., Weber, C. T. M., Cingöz, A., Torgerson, J. R. & Budker, D. New limits on variation of the finestructure constant using atomic dysprosium. Phys. Rev. Lett. 111, 060801 (2013).
 15.
Pashov, A. et al. Coupling of the X ^{1}Σ^{+} and a ^{3}Σ^{+} states of KRb. Phys. Rev. A 76, 022511 (2007).
 16.
Beuc, R. et al. Predictions for the observation of KRb spectra under cold conditions. J. Phys. B 39, S1191 (2006).
 17.
Kotochigova, S., Julienne, P. S. & Tiesinga, E. Ab initio calculation of KRb dipole moments. Phys. Rev. A 68, 022501 (2003).
 18.
Kobayashi, J., Aikawa, K., Oasa, K. & Inouye, S. Prospects for narrowline cooling of KRb molecules in the rovibrational ground state. Phys. Rev. A 89, 021401 (2014).
 19.
Zelevinsky, T., Kotochigova, S. & Ye, J. Precision test of massratio variations with latticeconfined ultracold molecules. Phys. Rev. Lett. 100, 043201 (2008).
 20.
Arnold, K. J., Kaewuam, R., Roy, A., Tan, T. R. & Barrett, M. D. Blackbody radiation shift assessment for lutetium in clock. Nat. Comm. 9, 1650 (2018).
 21.
Townes, C. H. & Schawlow, A. L. Microwave Spectroscopy. (Dover Publications, New York, 1975).
Acknowledgements
This work was supported by a GrantinAid for Young Scientists (A) of JSPS (No. 23684034), a GrantinAid for Scientific Research on Innovative Areas of JSPS (No. 24104702), JST PRESTO Grant Number JPMJPR17P5, Japan and the Matsuo Foundation.
Author information
Affiliations
Contributions
J.K. and A.O. performed the experiment and analyzed the data. J.K. and S.I. supervised the project and prepared the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information: Nature Communications would like to thank Colin Parker, Marianna Safronova, and other anonymous reviewers for their contributions to the peer review of this work.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kobayashi, J., Ogino, A. & Inouye, S. Measurement of the variation of electrontoproton mass ratio using ultracold molecules produced from lasercooled atoms. Nat Commun 10, 3771 (2019). https://doi.org/10.1038/s41467019117611
Received:
Accepted:
Published:
Further reading

Simultaneous Zeeman deceleration of polyatomic free radical with lithium atoms
Frontiers of Physics (2021)

Molecular collisions: From nearcold to ultracold
Frontiers of Physics (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.