{'fraction_non_alphanumeric': 0.0452707017668346, 'fraction_numerical': 0.022990136947420706, 'mean_word_length': 4.100977198697068, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Limits on temporal variation of fine structure constant, quark masses and strong interaction from atomic clock experiments

We perform calculations of the dependence of nuclear magnetic moments on quark masses and obtain limits on the variation of the ratio of quark mass and strong interaction scale (m_q/Lambda_{QCD}) from recent atomic clock experiments with hyperfine transitions in H, Rb, Cs, Yb+, Hg+ and optical transition in Hg+. Experiments with Cd+, deuterium/hydrogen, molecule SF_6, Zeeman transitions in He_3/Xe are also discussed.


I. INTRODUCTION
Interest in the temporal and spatial variation of major constants of physics has been recently revived by astronomical data which seem to suggest a variation of the electromagnetic constant α = e 2 /hc at the 10 −5 level for the time scale 10 billion years, see [1] (a discussion of other limits can be found in the review [2] and references therein). However, an independent experimental confirmation is needed.
The hypothetical unification of all interactions implies that variation of the electromagnetic interaction constant α should be accompanied by the variation of masses and the strong interaction constant. Specific predictions need a model. For example, the grand unification model discussed in [3] predicts that the quantum chromodynamic (QCD) scale Λ QCD (defined as the position of the Landau pole in the logarithm for the running strong coupling constant) is modified as follows: δΛ QCD /Λ QCD ≈ 34δα/α. The variation of quark and electron masses in this model is given by δm/m ∼ 70δα/α. This gives an estimate for the variation of the dimensionless ratio This result is strongly model-dependent. However, the large coefficients in these expressions are generic for grand unification models, in which modifications come from high energy scales: they appear because the running strong coupling constant and Higgs constants (related to mass) run faster than α. This means that if these models are correct the variation of masses and strong interaction may be easier to detect than the variation of α.
One can measure only variation of the dimensionless quantities, therefore we want to extract from the measurements variation of the dimensionless ratio m q /Λ QCD where m q is the quark mass (with the dependence on the normalization point removed). A number of limits on variation of m q /Λ QCD have been obtained recently from consideration of Big Bang Nucleosynthesis, quasar absorption spectra and Oklo natural nuclear reactor which was active about 1.8 billion years ago [4,5,6,7] (see also [8,9,10,11,12]). Below we consider the limits which follow from laboratory atomic clock comparison. Laboratory limits with a time base about a year are especially sensitive to oscillatory variation of fundamental constants. A number of relevant measurements have been performed already and even larger number have been started or planned. The increase in precision is very fast.
It has been pointed out by Karshenboim [13] that measurements of ratio of hyperfine structure intervals in different atoms are sensitive to variation of nuclear magnetic moments. First rough estimates of the dependence of nuclear magnetic moments on m q /Λ QCD and limits on time variation of this ratio have been obtained in our paper [4]. Using H, Cs and Hg + measurements [14,15], we obtained the limit on variation of m q /Λ QCD about 5 · 10 −13 per year. Below we calculate the dependence of nuclear magnetic moments on m q /Λ QCD and obtain the limits from recent atomic clock experiments with hyperfine transitions in H, Rb, Cs,Yb + ,Hg + and optical transition in Hg + . It is convenient to assume that the strong interaction scale Λ QCD does not vary, so we will speak about variation of masses. We shall restore Λ QCD in final answers.
The hyperfine structure constant can be presented in the following form The factor in the first bracket is an atomic unit of energy. The second "electromagnetic" bracket determines the dependence on α. An approximate expression for the relativistic correction factor (Casimir factor) for s-wave electron is the following where γ = 1 − (Zα) 2 , Z is the nuclear charge. Variation of α leads to the following variation of F rel [14]: More accurate numerical many-body calculations [16] of the dependence of the hyperfine structure on α have shown that the coefficient K is slightly larger than that given by this formula. For Cs (Z=55) K= 0.83 (instead of 0.74), for Rb K=0.34 (instead of 0.29), for Hg + K=2.28 (instead of 2.18). The last bracket in eq. (2) contains the dimensionless nuclear magnetic moment µ in nuclear magnetons ( the nuclear magnetic moment M = µ eh 2mpc ), electron mass m e and proton mass m p . We may also include a small correction due to the finite nuclear size. However, its contribution is insignificant.
Recent experiments measured time dependence of the ratios of hyperfine structure intervals of 199 Hg + and H [14], 133 Cs and 87 Rb [17] and ratio of optical frequency in Hg + and 133 Cs hyperfine frequency [18]. In the ratio of two hyperfine structure constants for different atoms time dependence may appear from the ratio of the factors F rel (depending on α) and ratio of nuclear magnetic moments (depending on m q /Λ QCD ). Magnetic moments in a single-particle approximation (one unpaired nucleon) are: for j = l − 1/2. Here the orbital g-factors are g l = 1 for valence proton and g l = 0 for valence neutron. The present values of spin g-factors g s are g p = 5.586 for proton and g n = −3.826 for neutron. They depend on m q /Λ QCD . The light quark masses are only about 1% of the nucleon mass (m q = (m u + m d )/2 ≈ 5 MeV). The nucleon magnetic moment remains finite in the chiral limit of m u = m d = 0. Therefore, one may think that the corrections to g s due to the finite quark masses are very small. However, there is a mechanism which enhances quark mass contribution: π-meson loop corrections to the nucleon magnetic moments which are proportional to π-meson mass m π ∼ m q Λ QCD ; m π =140 MeV is not so small. According to calculation in Ref. [19] dependence of the nucleon g-factors on π-meson mass m π can be approximated by the following equation where a= 1.37/GeV , b= 0.452/GeV 2 for proton and a= 1.85/GeV , b= 0.271/GeV 2 for neutron. This formula is a fit of the numerical calculations performed using the Cloudy Bag Model (pion field coupled to nucleon quarks). This model reproduces leading terms of chiral perturbation theory for nucleon magnetic moments [20] (small m π limit). The results also agree with lattice calculations [19] (performed at m π ∼ 500 MeV). Eq. (8) leads to the following estimates: δg n g n = −0.213 δm π m π = −0.107 δm q m q Eqs. (6,7,9,10) give variation of nuclear magnetic moments as functions of m q /Λ QCD . The measured time variations of ratios of hyperfine frequencies depend on two parameters: the ratio of proton spin magnetic moment and proton orbital magnetic moment (this ratio is proportional to g p ) and the ratio of proton and neutron spin magnetic moments (M p /M n = g p /g n ). According to eqs. (9, 10) the ratio g p /g n practically does not depend on m q and seems to be not sensitive to variation of quark masses and strong interaction. However, this conclusion may be misleading. Magnetic moments depend also on the strange quark mass m s . In the minimal order of the chiral perturbation theory this dependence is dominated by mixing p − K + Λ. Similar process involving K 0 meson loop and Λ does not contribute to the neutron magnetic moment since K 0 does not carry electric charge. As a result lowest order corrections for proton and neutron magnetic moments are very different: for proton a(K)/a(π) = 0.4, for neutron a(K)/a(π) = −0.03 [20,21] (see eq. (8) for definition of the coefficient a). The proton (or neutron) g-factor is a ratio of the proton (or neutron) spin magnetic moment to the nuclear magneton eh/2m p c (quantum for proton orbital magnetic moment). Therefore, dependence of g p ∼ M p m p and g n ∼ M n m p on the strange quark mass may also appear from the relatively large (∼ 20%) contribution of m s to the proton mass m p [4,6]. In the case of g p this contribution is of opposite sign to that of the lowest order in the chiral perturbation theory (a(K)m K ) and has comparable magnitude. This cancellation means that we can not reliably determine dependence of g p on the strange quark mass m s . In this case it would be safer to neglect contribution of the strange quark mass to g p . However, we want to keep dependence on the strange quark mass in the ratio of neutron and proton magnetic moments since in a number of cases this is the only effect which provides dependence on fundamental masses and strong interaction. Using the lowest order chiral perturbation theory results [20,21] presented above we obtain the dependence on the strange quark mass δ(gn/gp) gn/gp ≈ 0.1 δms ms . This result should possibly be treated as an order of magnitude estimate since K-meson mass is not as small as π-meson mass and accuracy of the chiral theory may be low in this case. Thus, we arrive to the following equations: Note that transferring the entire m s -dependence into the neutron g-factor is mostly a matter of convenience, all that we can say is that g n /g p depends on m s /Λ QCD . Now we can find variation of nuclear magnetic moments using eqs. (6,7). For all even Z nuclei with valence neutron ( 199 Hg, 171 Yb, 111 Cd, etc) we obtain δµ µ = δgn gn . For 133 Cs we have valence proton with j=7/2, l=4 and δµ µ = 0.22 δm π m π = 0.11 δm q m q (13) For 87 Rb we have valence proton with j=3/2, l=1 and As an intermediate result it is convenient to present dependence of the ratio of hyperfine constant A to atomic unit of energy E = mee 4 h 2 (or energy of 1s-2s transition in hydrogen equal to 3/8 E) on variation of the fundamental constants. We introduce a parameter V defined by the realation We start from the hyperfine structure of 133 Cs which is used as a frequency standard. Using eqs.(2,13) we obtain Here we have taken into account that the proton mass m p ∝ Λ QCD . For hyperfine transition frequencies in other atoms we obtain Now we can use these results to find limits on variation of the fundamental constants from the measurements of the time dependence of hyperfine structure intervals. The dependence of ratio of frequencies A( 133 Cs)/A( 87 Rb) can be presented in the following form Therefore, the result of the measurement [17] may be presented as a limit on variation of the parameter X: (24) Note that if the relation (1) is correct, variation of X(CsRb) would be dominated by variation of [m q /Λ QCD ]. The relation (1) would give X(CsRb) ∝ α 7 . For A( 133 Cs)/A(H) we have Therefore, the result of the measurements [15] may be presented as For A( 199 Hg)/A(H) we have The result of measurement [14] may be presented as In Ref. [13] the limit on variation of the ratio of hyperfine transition frequencies 171 Yb + / 133 Cs has been obtained (this limit is based on measurements [22]). Using eqs. (16,21) we can present the result as a limit on X(Y bCs) = α 0.7 [m s /Λ QCD ] 0.11 : The optical clock transition energy E(Hg) (λ=282 nm) in Hg + ion can be presented in the following form: Numerical calculation of the relative variation of E(Hg) has given [16]: This corresponds to V (HgOpt) = α −3.2 . Variation of the ratio of the Cs hyperfine splitting A(Cs) to this optical transition energy is described by (32) The work [18] gives the limit on variation of this parameter: Molecular rotational transitions frequencies are proportional to m e /m p , therefore one should assume V = m e /Λ QCD . For vibrational molecular transtions V = (m e /Λ QCD ) 1/2 .
Comparison of Cs hyperfine standard with SF 6 molecular vibration frequencies was discussed in Ref. [23]. In this case X(CsV ibrations) = The measurements of hyperfine constant ratio in different isotopes of the same atom depends on the ratio of magnetic moments and is sensitive to m q /Λ QCD . For example, it would be interesting to measure the rate of change for hydrogen/deuterium ratio where X(HD) = [m q /Λ QCD ] −0.05 [m s /Λ QCD ] 0.23 .
R. Walsworth suggested to measure the ratio of the Zeeman transition frequencies in noble gases which gives us time dependence of ratio of nuclear magnetic moments. Consider, for example 129 Xe/ 3 He. For 3 He the magnetic moment is very close to that of neutron. For other noble gases nuclear magnetic moment is also given by valence neutron, however, there are significant manybody corrections. For 129 Xe the valence neutron is in s 1/2 state which corresponds to the single-particle value of nuclear magnetic moment µ = µ n = −1.913. The measured value is µ = −0.778. The magnetic moment of the nucleus changes most efficiently due to the spin-spin interaction because valence neutron transfers a part of its spin < s z > to core protons and proton magnetic moment is large and has opposite sign. In this approximation µ = (1 − b)µ n + bµ p . This gives b=0.24 and the ratio of magnetic moments Y ≡ µ( 129 Xe)/µ( 3 He)≈ 0.76 + 0.24g p /g n . Using eqs. (11,12) we obtain an estimate for the variation of this ratio δY /Y ≈ 0.1δm s /m s . Therefore, in this case one can measure variation of µ( 129 Xe)/µ( 3 He) corresponding to variation of X = [m s /Λ QCD ] 0.1 .
Note that an accuracy of the results presented in this paper depends strongly on fundamental constant we are studying. The accuracy for the dependence on α is few percent. The accuracy for m q /Λ QCD is about 30%, it is limited by the accuracy of a single-particle approximation for nuclear magnetic moments. The results for m s /Λ QCD should possibly be treated as order of magnitude estimates. However, the accuracy here may be improved using proper QCD calculation [24]. Relation (1) between variation of α and m/Λ QCD has been used as an illustration only.
The author is grateful to C. Chardonnet, S. Karshenboim, D.B. Leinweber, A.W. Thomas and R. Walsworth for valuable discussions. This work is supported by the Australian Research Council.