<?xml version="1.0" encoding="UTF-8"?>
<article xmlns="http://specifications.silverchair.com/xsd/article/1/0/SCJATS-journalpublishing1-0.xsd" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">ptep</journal-id>
<journal-title-group>
<journal-title>Progress of Theoretical and Experimental Physics</journal-title>
</journal-title-group>
<issn pub-type="epub">2050-3911</issn>
<publisher>
<publisher-name>Oxford University Press</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.1093/ptep/pty066</article-id>
<article-id pub-id-type="publisher-id">pty066</article-id>
<article-id pub-id-type="arxiv">arXiv:1801.10244</article-id>
<article-categories>
<subj-group subj-group-type="category-toc-heading">
<subject>Papers</subject>
<subj-group subj-group-type="category-toc-heading">
<subject>Theoretical Particle Physics</subject>
</subj-group>
</subj-group>
<subj-group subj-group-type="category-journal-collection">
<subject>PTEP/B30</subject>
<subject>PTEP/B31</subject>
<subject>PTEP/E00</subject>
<subject>PTEP/E01</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>The general relativistic effects on the magnetic moment in Earth&#x2019;s gravity</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Morishima</surname><given-names>Takahiro</given-names></name>
<xref ref-type="aff" rid="AFF1"/>
<xref ref-type="corresp" rid="COR1"/>
<email xlink:type="simple">moris@nagoya-u.jp</email>
</contrib>
<contrib contrib-type="author">
<name><surname>Futamase</surname><given-names>Toshifumi</given-names></name>
<xref ref-type="fn" rid="FN100"/>
<xref ref-type="aff" rid="AFF2"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Shimizu</surname><given-names>Hirohiko M</given-names></name>
<xref ref-type="aff" rid="AFF1"/>
<xref ref-type="fn" rid="FN100"/>
</contrib>
</contrib-group>
<aff id="AFF1">Laboratory for Particle Properties, Department of Physics, Nagoya University, Furocho, Chikusa, Nagoya 464-8602, Japan</aff>
<aff id="AFF2">Department of Astrophysics and Atmospheric Science, Kyoto Sangyo University, Motoyama, Kamigamo, Kita, Kyoto 603-8555, Japan</aff>
<author-notes>
<corresp id="COR1">E-mail: <email>moris@nagoya-u.jp</email></corresp>
<fn id="FN100"><p>Toshifumi Futamase, <email>tof@cc.kyoto-su.ac.jp</email>; Hirohiko M Shimizu, <email>hirohiko.shimizu@nagoya-u.jp</email></p></fn>
</author-notes>
<pub-date pub-type="cover">
<month>06</month>
<year>2018</year>
</pub-date>
<pub-date pub-type="collection">
<day>01</day>
<month>06</month>
<year>2018</year>
</pub-date>
<pub-date pub-type="epub" iso-8601-date="2018-06-28">
<day>28</day>
<month>06</month>
<year>2018</year>
</pub-date>
<volume>2018</volume>
<issue>6</issue>
<elocation-id>063B07</elocation-id>
<history>
<date date-type="received">
<day>31</day>
<month>01</month>
<year>2018</year>
</date>
<date date-type="rev-recd">
<day>27</day>
<month>04</month>
<year>2018</year>
</date>
<date date-type="accepted">
<day>11</day>
<month>05</month>
<year>2018</year>
</date>
</history>
<permissions>
<copyright-statement>&#x00A9; The Author(s) 2018. Published by Oxford University Press on behalf of the Physical Society of Japan.</copyright-statement>
<copyright-year>2018</copyright-year>
<license license-type="cc-by" xlink:href="http://creativecommons.org/licenses/by/4.0/">
<license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
<license-p>Funded by SCOAP<sup>3</sup></license-p>
</license>
</permissions>
<self-uri xlink:href="pty066.pdf"/>
<abstract abstract-type="abstract"><title>Abstract</title>
<p>The magnetic moment of free fermions in the Earth&#x2019;s gravitational field has been studied on the basis of general relativity. Adopting the Schwarzschild metric for the background spacetime, the dipole coupling between the magnetic moment and the magnetic field has been found to be dependent on the gravity in the calculation up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation1"><![CDATA[$O(1/c^2)$]]></tex-math></inline-formula>. The gravity dependence can be formulated by employing the effective value of the magnetic moment as a gravity-dependent quantity <inline-formula><tex-math notation="LaTeX" id="ImEquation2"><![CDATA[$\boldsymbol{\mu}_{\rm m}^{\rm eff}= (1+3\phi/c^2) \,\,\boldsymbol{\mu}_{\rm m}$]]></tex-math></inline-formula> for the cases of minimal coupling, non-minimal coupling, and a mixture of the two. The gravitationally induced anomaly is found to be canceled in the experimental values of the anomalous magnetic moment measured in the Penning trap and storage ring methods.</p>
</abstract>
<kwd-group kwd-group-type="jel">
<kwd>B30</kwd>
<kwd>B31</kwd>
<kwd>E00</kwd>
<kwd>E01</kwd>
</kwd-group>
<counts>
<page-count count="25"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="SEC1"><title>1. Introduction</title>
<p>The magnetic moment is regarded as one of the fundamental properties of elementary particles. The magnetic moment of fermions with a rest mass of <inline-formula><tex-math notation="LaTeX" id="ImEquation3"><![CDATA[$m$]]></tex-math></inline-formula>, a charge of <inline-formula><tex-math notation="LaTeX" id="ImEquation4"><![CDATA[$e$]]></tex-math></inline-formula>, and a spin of <inline-formula><tex-math notation="LaTeX" id="ImEquation5"><![CDATA[$\boldsymbol{S} = \boldsymbol{\sigma}/2$]]></tex-math></inline-formula> is defined as
<disp-formula id="pty066-M-1"><label>(1)</label><tex-math notation="LaTeX" id="Equation1"><![CDATA[
\begin{eqnarray}
\boldsymbol{\mu}_{\rm m} \equiv \frac{{\rm g}}{2}\, \frac{e}{m} \boldsymbol{S}
,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation6"><![CDATA[${\rm g}$]]></tex-math></inline-formula> is the constant factor specific to each elementary particle, referred to as the gyromagnetic ratio or Land&#x00E9; g-factor<xref ref-type="fn" rid="FN1"><sup>1</sup></xref>. Here we consider particles with a spin of 1/2 such as electrons, protons, and atoms placed in the Earth&#x2019;s gravitational field and consider if their magnetic moments, defined by Eq. (<xref ref-type="disp-formula" rid="pty066-M-1">1</xref>), are independent of gravity in the same way as the mass or the electric charge. We focus on the influence of the spacetime distortion on the basis of general relativity, which is not a quantum mechanical effect mediated by the graviton at the Planck scale but a classical mechanical effect remaining even in low-energy regions.</p>
<p>Below, we consider observables measured using instruments fixed on the Earth&#x2019;s surface. General relativity requires that the translational motion of a charged particle with an electric charge of <inline-formula><tex-math notation="LaTeX" id="ImEquation7"><![CDATA[$e$]]></tex-math></inline-formula> obeys the geodesic equation
<disp-formula id="pty066-M-2"><label>(2)</label><tex-math notation="LaTeX" id="Equation2"><![CDATA[
\begin{eqnarray}
\frac{D u^\mu}{d \tau}
&\equiv&
\frac{d u^\mu}{d \tau} + \Gamma^\mu_{\nu\lambda}u^\nu u^\lambda
= \frac{e}{m} F^{\mu\nu}u_\nu
,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation8"><![CDATA[$\tau$]]></tex-math></inline-formula> is the proper time, <inline-formula><tex-math notation="LaTeX" id="ImEquation9"><![CDATA[$u^\mu=(u^0, \boldsymbol{u})$]]></tex-math></inline-formula> is the four velocity vector, <inline-formula><tex-math notation="LaTeX" id="ImEquation10"><![CDATA[$F^{\mu\nu}$]]></tex-math></inline-formula> is the electromagnetic tensor, and <inline-formula><tex-math notation="LaTeX" id="ImEquation11"><![CDATA[$\Gamma^\mu_{\nu\lambda}$]]></tex-math></inline-formula> is the Christoffel symbol. The Earth&#x2019;s gravitational field can be described by the Schwarzschild metric, approximated up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation12"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> in the expansion series of <inline-formula><tex-math notation="LaTeX" id="ImEquation13"><![CDATA[$\epsilon = 1/c$]]></tex-math></inline-formula>, as
<disp-formula id="pty066-M-3"><label>(3)</label><tex-math notation="LaTeX" id="Equation3"><![CDATA[
\begin{eqnarray}
ds^2 &=& g_{\mu\nu}\, dx^{\mu}\, dx^{\nu} \nonumber\\
&=& \epsilon^{-2}\bigl(1 + \epsilon^2 2\phi + \epsilon^4 2\phi^2\bigr) dt^2
- (1 - \epsilon^2 2\phi)\,\, (dx^2 + dy^2 + dz^2) + O(\epsilon^4),
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation14"><![CDATA[$\phi = -GM/r$]]></tex-math></inline-formula> is the Earth&#x2019;s gravitational potential at a distance of <inline-formula><tex-math notation="LaTeX" id="ImEquation15"><![CDATA[$r = \sqrt{x^2 + y^2 + z^2}$]]></tex-math></inline-formula> from the center of a uniformly distributed spherical mass <inline-formula><tex-math notation="LaTeX" id="ImEquation16"><![CDATA[$M$]]></tex-math></inline-formula> with a radius of <inline-formula><tex-math notation="LaTeX" id="ImEquation17"><![CDATA[$R$]]></tex-math></inline-formula> (see <xref ref-type="sec" rid="APP2">Appendix B</xref>). Equation (<xref ref-type="disp-formula" rid="pty066-M-2">2</xref>) becomes
<disp-formula id="pty066-M-4"><label>(4)</label><tex-math notation="LaTeX" id="Equation4"><![CDATA[
\begin{eqnarray}
\frac{d \boldsymbol{\beta}}{d t}
&=&
\big(1 + (2\gamma^2 + 1)\epsilon^2\phi\big) \,
\frac{e}{\gamma\,m}\,
\big(
\boldsymbol{E}
+ \boldsymbol{\beta} \times \boldsymbol{B}
- (1 - 4\gamma^2\epsilon^2\phi)\,\boldsymbol{\beta}\,(\boldsymbol{\beta} \cdot \boldsymbol{E})
\big)
\nonumber\\
& & \hspace{10mm}
- \epsilon (1 + \beta^2)\nabla\phi
+ \epsilon (4\boldsymbol{\beta} \cdot \nabla\phi)\boldsymbol{\beta}
\,\,\,\, + O(\epsilon^4)
,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation18"><![CDATA[$\boldsymbol{\beta} = \boldsymbol{v}/c$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation19"><![CDATA[$\gamma = 1/\sqrt{1 - \beta^2}$]]></tex-math></inline-formula>, and
<disp-formula id="pty066-M-5"><label>(5)</label><tex-math notation="LaTeX" id="Equation5"><![CDATA[
\begin{eqnarray}
u^0 &=& \frac{dt}{d\tau}
= 1/\sqrt{(1+2\epsilon^2\phi+2\epsilon^4\phi^2)-(1-2\epsilon^2\phi)\beta^2}
\nonumber\\
&=& \gamma \big( 1-(2\gamma^2-1)\epsilon^2\phi \big) \,\,\, +O(\epsilon^4).
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Here we suggest that the <inline-formula><tex-math notation="LaTeX" id="ImEquation20"><![CDATA[${\bf E}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation21"><![CDATA[${\bf B}$]]></tex-math></inline-formula> represent the electric and magnetic fields to be measured by the observer on the Earth&#x2019;s surface, and correspond to the field vectors in the flat spacetime.</p>
<p>At the limit of <inline-formula><tex-math notation="LaTeX" id="ImEquation22"><![CDATA[$\phi \rightarrow 0$]]></tex-math></inline-formula>, Eq. (<xref ref-type="disp-formula" rid="pty066-M-2">2</xref>) results in
<disp-formula id="pty066-M-6"><label>(6)</label><tex-math notation="LaTeX" id="Equation6"><![CDATA[
\begin{eqnarray}
\frac{d \boldsymbol{\beta}}{d t}
&=&
\frac{e}{\gamma\, m}\,
\bigl(
\boldsymbol{E}
+ \boldsymbol{\beta} \times \boldsymbol{B}
- (\boldsymbol{E} \cdot \boldsymbol{\beta}) \boldsymbol{\beta}
\bigr)
,
\end{eqnarray}]]></tex-math></disp-formula>
which is the well-known equation of motion of charged particles in the flat spacetime (Minkowski spacetime) consistently with special relativity. Additionally, in the region of <inline-formula><tex-math notation="LaTeX" id="ImEquation23"><![CDATA[$\beta = v/c \ll 1$]]></tex-math></inline-formula>, it results in the Lorentz equation of motion
<disp-formula id="pty066-M-7"><label>(7)</label><tex-math notation="LaTeX" id="Equation7"><![CDATA[
\begin{eqnarray}
\frac{d (m\boldsymbol{v})}{d t} = e ( \boldsymbol{E} + \boldsymbol{v} \times \boldsymbol{B})
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>For simplicity, we consider the case where the electric field is zero and the motion is limited on the horizontal plane and obtain
<disp-formula id="pty066-M-8"><label>(8)</label><tex-math notation="LaTeX" id="Equation8"><![CDATA[
\begin{eqnarray}
\frac{d \boldsymbol{\beta}}{d t}
&\simeq&
(1 + 3\epsilon^2\phi) \, \frac{e}{m} \boldsymbol{\beta} \times \boldsymbol{B}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>According to the correspondence in classical mechanics, the magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation24"><![CDATA[$\boldsymbol{\mu}_{\rm m}$]]></tex-math></inline-formula> can be related to the angular momentum of a charged particle moving on a circle under a constant magnetic field, which leads to <inline-formula><tex-math notation="LaTeX" id="ImEquation25"><![CDATA[$\left\vert\boldsymbol{\mu}_{\rm m}\right\vert=e/(2m)$]]></tex-math></inline-formula>. Equation (<xref ref-type="disp-formula" rid="pty066-M-8">8</xref>) can be written as
<disp-formula id="pty066-M-9"><label>(9)</label><tex-math notation="LaTeX" id="Equation9"><![CDATA[
\begin{eqnarray}
\frac{d \boldsymbol{\beta}}{d t}
\simeq
(1 + 3\epsilon^2\phi) \,\, 2\left\vert\boldsymbol{\mu}_{\rm m}\right\vert \boldsymbol{\beta} \times \boldsymbol{B}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The factor <inline-formula><tex-math notation="LaTeX" id="ImEquation26"><![CDATA[$(1+3\epsilon^2\phi)$]]></tex-math></inline-formula> represents the effect of the spacetime curvature due to the Earth&#x2019;s gravitational field. Here we note that <inline-formula><tex-math notation="LaTeX" id="ImEquation27"><![CDATA[$\boldsymbol{\mu}_{\rm m}$]]></tex-math></inline-formula> is defined in the flat spacetime, while the magnetic field <inline-formula><tex-math notation="LaTeX" id="ImEquation28"><![CDATA[$\boldsymbol{B}$]]></tex-math></inline-formula> is measured by the observer on the Earth&#x2019;s surface. Therefore, in this paper, we define the effective magnetic moment as
<disp-formula id="pty066-M-10"><label>(10)</label><tex-math notation="LaTeX" id="Equation10"><![CDATA[
\begin{eqnarray}
\boldsymbol{\mu}_{\rm m}^{\rm eff} \simeq (1 + 3\epsilon^2\phi) \,\, \boldsymbol{\mu}_{\rm m} ,
\end{eqnarray}]]></tex-math></disp-formula>
to include the gravitational effect. The above discussion remains within the classical picture since the lepton is not treated as a fermion with quantized spin and it is not trivial to extend the influence of gravity for quantum mechanical particles. Therefore the effective magnetic moment should be examined based on the formalism compatible with both quantum mechanics and general relativity.</p>
<p>In quantum field theory, the magnetic moment of a fermion is related to the form factors <inline-formula><tex-math notation="LaTeX" id="ImEquation29"><![CDATA[$F_1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation30"><![CDATA[$F_2$]]></tex-math></inline-formula> defined in the vertex function
(see Refs. [<xref ref-type="bibr" rid="B1">1</xref>,<xref ref-type="bibr" rid="B2">2</xref>])
<disp-formula id="pty066-M-11"><label>(11)</label><tex-math notation="LaTeX" id="Equation11"><![CDATA[
\begin{eqnarray}
{\cal M}^\mu (p',p) =
\gamma^\mu F_1(q^2)
+
\frac{i}{2m} \sigma^{\mu\nu}q_\nu F_2(q^2)
,
\end{eqnarray}]]></tex-math></disp-formula>
and the g-factor in the flat spacetime is defined as
<disp-formula id="pty066-M-12"><label>(12)</label><tex-math notation="LaTeX" id="Equation12"><![CDATA[
\begin{eqnarray}
{\rm g}
= 2\,(F_1(0)+F_2(0))
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The origin of the magnetic moment varies according to the type of fermions as <inline-formula><tex-math notation="LaTeX" id="ImEquation31"><![CDATA[$F_1(0) = 1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation32"><![CDATA[$F_2(0) = 0$]]></tex-math></inline-formula> for Dirac particles, i.e., fermions with <inline-formula><tex-math notation="LaTeX" id="ImEquation33"><![CDATA[$\mathrm{g}=2$]]></tex-math></inline-formula>, the internal structure and radiative corrections of which can be ignored; <inline-formula><tex-math notation="LaTeX" id="ImEquation34"><![CDATA[$F_1(0) = 0$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation35"><![CDATA[$F_2(0) \ne 0$]]></tex-math></inline-formula> for neutrons and neutral atoms, i.e., fermions with <inline-formula><tex-math notation="LaTeX" id="ImEquation36"><![CDATA[$\mathrm{g}\ne2$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation37"><![CDATA[$e=0$]]></tex-math></inline-formula>; and <inline-formula><tex-math notation="LaTeX" id="ImEquation38"><![CDATA[$F_1(0) = 1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation39"><![CDATA[$F_2(0) \ne 0$]]></tex-math></inline-formula> for protons and ions, i.e., fermions with <inline-formula><tex-math notation="LaTeX" id="ImEquation40"><![CDATA[$\mathrm{g}\ne2$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation41"><![CDATA[$e\ne0$]]></tex-math></inline-formula>. We applied quantum mechanics in the curved spacetime in general relativity for the three types:
<list list-type="simple">
<list-item><p><bold>Type A:</bold> <inline-formula><tex-math notation="LaTeX" id="ImEquation42"><![CDATA[$\quad$]]></tex-math></inline-formula> Dirac particles with <inline-formula><tex-math notation="LaTeX" id="ImEquation43"><![CDATA[$\mathrm{g}=2$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation44"><![CDATA[$e\ne0$]]></tex-math></inline-formula> (<inline-formula><tex-math notation="LaTeX" id="ImEquation45"><![CDATA[$F_1(0)=1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation46"><![CDATA[$F_2(0)=0 $]]></tex-math></inline-formula>),</p></list-item>
<list-item><p><bold>Type B:</bold> <inline-formula><tex-math notation="LaTeX" id="ImEquation47"><![CDATA[$\quad$]]></tex-math></inline-formula> fermions with <inline-formula><tex-math notation="LaTeX" id="ImEquation48"><![CDATA[$\mathrm{g}\ne2$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation49"><![CDATA[$e=0$]]></tex-math></inline-formula> <inline-formula><tex-math notation="LaTeX" id="ImEquation50"><![CDATA[$\qquad$]]></tex-math></inline-formula> (<inline-formula><tex-math notation="LaTeX" id="ImEquation51"><![CDATA[$F_1(0)=0$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation52"><![CDATA[$F_2(0)\ne$]]></tex-math></inline-formula> 0),</p></list-item>
<list-item><p><bold>Type C:</bold> <inline-formula><tex-math notation="LaTeX" id="ImEquation53"><![CDATA[$\quad$]]></tex-math></inline-formula> fermions with <inline-formula><tex-math notation="LaTeX" id="ImEquation54"><![CDATA[$\mathrm{g}\ne2$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation55"><![CDATA[$e\ne0$]]></tex-math></inline-formula> <inline-formula><tex-math notation="LaTeX" id="ImEquation56"><![CDATA[$\qquad$]]></tex-math></inline-formula> (<inline-formula><tex-math notation="LaTeX" id="ImEquation57"><![CDATA[$F_1(0)=1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation58"><![CDATA[$F_2(0)\ne$]]></tex-math></inline-formula> 0).</p></list-item>
</list></p>
<p>The coupling between the fermions and the electromagnetic field corresponds to minimal coupling for Type A, non-minimal coupling for Type B, and a mixture of minimal and non-minimal coupling for Type C.</p>
<p>In this paper, we generalize the Dirac equation in curved spacetime by employing the covariant formulation. For simplicity, we employ the Schwarzschild metric as the background spacetime, derive the Hamiltonian of free fermions to be observed using instruments fixed on the Earth&#x2019;s surface up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation59"><![CDATA[$O(1/c^2)$]]></tex-math></inline-formula> for the above three types and evaluate the effective value of the magnetic moment. We introduce the <italic>effective</italic> magnetic moment of fermions in the Earth&#x2019;s gravitational field and discuss the influence on precision measurements of the anomalous magnetic moment of electrons and muons.</p>
</sec>
<sec id="SEC2"><title>2. Effective magnetic moment in curved spacetime</title>
<p>The influence of general relativity in the Hamiltonian of fermions has been considered and reported in several papers [<xref ref-type="bibr" rid="B3">3</xref>&#x2013;<xref ref-type="bibr" rid="B5">5</xref>]. Wajima et al. [<xref ref-type="bibr" rid="B5">5</xref>] employed the covariant Dirac equation
<disp-formula id="pty066-M-13"><label>(13)</label><tex-math notation="LaTeX" id="Equation13"><![CDATA[
\begin{equation}
\bigl( \gamma^\mu (i\nabla_\mu ) - \epsilon^{-1} m \bigr) \Psi=0
\end{equation}]]></tex-math></disp-formula>
and regarded its time component
<disp-formula id="pty066-M-14"><label>(14)</label><tex-math notation="LaTeX" id="Equation14"><![CDATA[
\begin{eqnarray}
{\cal H}_{}
&\equiv& i\partial_0 \nonumber\\
&=& - i \Gamma_0 + \frac{1}{g^{00}}\gamma^0\gamma^i (-i\partial_i - i \Gamma_i)
+ \epsilon^{-1}\frac{1}{g^{00}}\gamma^0 m
\,\,\,\,\,\,\,\,
\end{eqnarray}]]></tex-math></disp-formula>
as the Hamiltonian and calculated its expectation value up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation60"><![CDATA[$O(1/c^2)$]]></tex-math></inline-formula>, where <inline-formula><tex-math notation="LaTeX" id="ImEquation61"><![CDATA[$\Psi$]]></tex-math></inline-formula> is the spinor function, <inline-formula><tex-math notation="LaTeX" id="ImEquation62"><![CDATA[$\nabla_\mu = \partial_\mu + \Gamma_\mu$]]></tex-math></inline-formula> the covariant derivative, <inline-formula><tex-math notation="LaTeX" id="ImEquation63"><![CDATA[$\Gamma_\mu$]]></tex-math></inline-formula> the spin connection, <inline-formula><tex-math notation="LaTeX" id="ImEquation64"><![CDATA[$\gamma^\mu$]]></tex-math></inline-formula> the <inline-formula><tex-math notation="LaTeX" id="ImEquation65"><![CDATA[$\gamma$]]></tex-math></inline-formula>-matrices in the curved spacetime (<xref ref-type="sec" rid="APP1">Appendix A</xref>), and <inline-formula><tex-math notation="LaTeX" id="ImEquation66"><![CDATA[$\epsilon = 1/c$]]></tex-math></inline-formula> a small quantity to be used as the variable for the power series expansion to identify the post-Newtonian approximation. Since the electric charge and the magnetic moment were not involved in their method, we develop and generalize their method to derive the Hamiltonian containing gravitational effects in general relativity on the basis of the covariant Dirac equation with the electric charge and the magnetic moment. Here we adopt the Heisenberg equation for the magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation67"><![CDATA[$\boldsymbol{\mu}_{\rm m}$]]></tex-math></inline-formula>
<disp-formula id="pty066-M-15"><label>(15)</label><tex-math notation="LaTeX" id="Equation15"><![CDATA[
\begin{eqnarray}
\frac{d\boldsymbol{S}}{dt}
= \frac{1}{i} \Bigl[ \boldsymbol{S}, {\cal H} \Bigr]
= \boldsymbol{\mu}_{\rm m} \times \boldsymbol{B}
\end{eqnarray}]]></tex-math></disp-formula>
as the definition of the magnetic moment in the flat spacetime. We generalize Eq. (<xref ref-type="disp-formula" rid="pty066-M-15">15</xref>) using the general relativistic Hamiltonian to evaluate the effective magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation68"><![CDATA[$\boldsymbol{\mu}_{\rm m}^{\rm eff}$]]></tex-math></inline-formula>.</p>
<sec id="SEC2.1"><title>2.1. Dirac particles with <inline-formula><tex-math notation="LaTeX" id="ImEquation69"><![CDATA[$\mathrm{g}=2$]]></tex-math></inline-formula> (Type A)</title> 
<p>In this section, we consider Type A, for which the form factors are <inline-formula><tex-math notation="LaTeX" id="ImEquation70"><![CDATA[$F_1(0) = 1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation71"><![CDATA[$F_2(0) = 0$]]></tex-math></inline-formula>, corresponding to point-like fermions, i.e., electrons and muons without radiative corrections. The covariant Dirac equation for this case is given with minimal coupling with the four potential <inline-formula><tex-math notation="LaTeX" id="ImEquation72"><![CDATA[$A_{\mu}$]]></tex-math></inline-formula> as
<disp-formula id="pty066-M-16"><label>(16)</label><tex-math notation="LaTeX" id="Equation16"><![CDATA[
\begin{eqnarray}
\bigl( \,\gamma^\mu\, i\,(\nabla_\mu + i\, e A_\mu ) - \epsilon^{-1} m \, \bigr) \,\,\Psi=0
,
\end{eqnarray}]]></tex-math></disp-formula>
which leads to the Hamiltonian
<disp-formula id="pty066-M-17"><label>(17)</label><tex-math notation="LaTeX" id="Equation17"><![CDATA[
\begin{eqnarray}
{\cal H}^{^{\rm A}}
&=&
- i\, \Gamma_0 + e A_0
+ \frac{1}{g^{00}}\gamma^0\gamma^i (-i\,\partial_i - i\,\Gamma_i + e A_i)
+ \epsilon^{-1} \frac{1}{g^{00}}\gamma^0 m
,
\end{eqnarray}]]></tex-math></disp-formula>
which is the general form of the Hamiltonian of Dirac particles with <inline-formula><tex-math notation="LaTeX" id="ImEquation73"><![CDATA[$\mathrm{g}=2$]]></tex-math></inline-formula> with the electric charge and the magnetic moment. It is uniquely determined for a given metric of the spacetime. Adopting the Schwarzschild metric (<xref ref-type="disp-formula" rid="pty066-M-3">3</xref>), we can rewrite Eq. (<xref ref-type="disp-formula" rid="pty066-M-17">17</xref>), after some algebra to rearrange it (see <xref ref-type="sec" rid="APP4">Appendix D</xref>), as
<disp-formula id="pty066-M-18"><label>(18)</label><tex-math notation="LaTeX" id="Equation18"><![CDATA[
\begin{eqnarray}
{\cal H}^{^{\rm A}}
&=& \epsilon^{-2} \rho_3 m
+ \epsilon^{-1} \rho_1 \boldsymbol{\sigma} \cdot ( \boldsymbol{p} - e \boldsymbol{A})
+ e A_0
+ \rho_3 m\phi
+ \epsilon\, \rho_1
\Bigl\{
\boldsymbol{\sigma} \cdot ( \boldsymbol{p} - e \boldsymbol{A}),\, \phi
\Bigr\}
+ \epsilon^2
\rho_3 \frac{m\phi^2}{2}
\nonumber\\
& &
\,\,\,\,\,+ O(\epsilon^4)
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>This Hamiltonian is a <inline-formula><tex-math notation="LaTeX" id="ImEquation74"><![CDATA[$4\times4$]]></tex-math></inline-formula> matrix. We can obtain a <inline-formula><tex-math notation="LaTeX" id="ImEquation75"><![CDATA[$2\times2$]]></tex-math></inline-formula> even matrix by applying the Foldy&#x2013;Wouthuysen&#x2013;Tani transformation [<xref ref-type="bibr" rid="B6">6</xref>&#x2013;<xref ref-type="bibr" rid="B8">8</xref>]. Its left-up <inline-formula><tex-math notation="LaTeX" id="ImEquation76"><![CDATA[$2\times2$]]></tex-math></inline-formula> matrix (the large component) <inline-formula><tex-math notation="LaTeX" id="ImEquation77"><![CDATA[${\cal H}_{+}^{^{\rm A}}$]]></tex-math></inline-formula> given as
<disp-formula id="pty066-M-19"><label>(19)</label><tex-math notation="LaTeX" id="Equation19"><![CDATA[
\begin{eqnarray}
{\cal H}_{+}^{^{\rm A}}
&=& \epsilon^{-2} m
\nonumber\\
& & + \frac{(\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}))^2}{2m} + m\phi
+ e A_0
\nonumber\\
& &
+ \epsilon^2
\Bigg(
\frac{m\phi^2}{2}
- \frac{(\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}))^4 }{8m^3}
+ \frac{3}{8m}
\Big\{ \boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}),\,
\big\{\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}),\,\phi \big\}\Big\}
\nonumber\\
& & \hspace{0.7cm}
- \frac{1}{8m^2}
\Big[ \boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}), \big[\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}), e A_0 \big] \Big]
\Bigg)+ O(\epsilon^4)
\end{eqnarray}]]></tex-math></disp-formula>
corresponds to the Hamiltonian of the Dirac particle in curved spacetime. Substituting <inline-formula><tex-math notation="LaTeX" id="ImEquation78"><![CDATA[$\boldsymbol{\mu}_{\rm m}=e\boldsymbol{S}/m=e\boldsymbol{\sigma}/2m$]]></tex-math></inline-formula>, Eq. (<xref ref-type="disp-formula" rid="pty066-M-19">19</xref>) leads to
<disp-formula id="pty066-M-20"><label>(20)</label><tex-math notation="LaTeX" id="Equation20"><![CDATA[
\begin{eqnarray}
{\cal H}_{+}^{^{\rm A}}
&=& \epsilon^{-2}\,m \nonumber\\
& & +
\frac{(\boldsymbol{p} - e \boldsymbol{A})^2}{2m}
- \boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B}
+ m\phi
+ e A_0
\nonumber\\
& &
+ \epsilon^2
\left(
\frac{m\phi^2}{2}
- \frac{(\boldsymbol{p} - e \boldsymbol{A})^4}{8m^3}
- \frac{(\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B})^2}{2m}
- \frac{1}{4m^2} \nabla^2(\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B})
- 3\phi\,\, \boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B}
\right.\nonumber\\
& & \hspace{0.7cm}
+ \frac{1}{2m^2}
(\boldsymbol{p} - e \boldsymbol{A}) \cdot \bigl(\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B}\, (\boldsymbol{p} - e \boldsymbol{A})\bigr)
+ \frac{3}{2m}
(\boldsymbol{p} - e\boldsymbol{A}) \cdot \bigl(\phi(\boldsymbol{p} - e\boldsymbol{A}) \bigr)
\nonumber\\
& & \hspace{0.7cm}
- \frac{3\phi}{2m r^2}\boldsymbol{S} \cdot \big(\boldsymbol{r} \times (\boldsymbol{p} - e\boldsymbol{A})\big)
\nonumber\\
& & \hspace{0.7cm}\left.
- \frac{{\mu}_{\rm m}}{4m}\,\, \nabla\cdot\left(\boldsymbol{E}+\frac{\partial\boldsymbol{A}}{\partial t} \right)
- \frac{1}{2m}\boldsymbol{\mu}_{\rm m} \cdot \left((\boldsymbol{E}+\frac{\partial \boldsymbol{A}}{\partial t}) \times (\boldsymbol{p} - e\boldsymbol{A})\right)\right) + O(\epsilon^4)
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Terms of the order of <inline-formula><tex-math notation="LaTeX" id="ImEquation79"><![CDATA[$\epsilon^{-2}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation80"><![CDATA[$\epsilon^{0}$]]></tex-math></inline-formula> are well known, while some of terms of the order of <inline-formula><tex-math notation="LaTeX" id="ImEquation81"><![CDATA[$\epsilon^{2}$]]></tex-math></inline-formula> are newly introduced coupling terms between the gravitational potential <inline-formula><tex-math notation="LaTeX" id="ImEquation82"><![CDATA[$\phi$]]></tex-math></inline-formula>, the magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation83"><![CDATA[$\boldsymbol{\mu}_{\rm m}$]]></tex-math></inline-formula>, and the magnetic field <inline-formula><tex-math notation="LaTeX" id="ImEquation84"><![CDATA[$\boldsymbol{B}$]]></tex-math></inline-formula>. The effective magnetic moment is given as
<disp-formula id="pty066-M-21"><label>(21)</label><tex-math notation="LaTeX" id="Equation21"><![CDATA[
\begin{eqnarray}
&& \boldsymbol{\mu}_{\rm m}^{\rm eff} \simeq \, (1+ 3\epsilon^2\phi) \,\,\boldsymbol{\mu}_{\rm m}
\hspace{1cm}
(\mbox{for Type A})
\end{eqnarray}]]></tex-math></disp-formula>
for sufficiently slow particles moving in a static spacetime. The result is consistent with the <italic>gravitationally modified magnetic moment</italic> obtained in the <inline-formula><tex-math notation="LaTeX" id="ImEquation85"><![CDATA[$TH\epsilon\mu$]]></tex-math></inline-formula> formalism [<xref ref-type="bibr" rid="B9">9</xref>,<xref ref-type="bibr" rid="B10">10</xref>]<xref ref-type="fn" rid="FN2"><sup>2</sup></xref>. The above discussions are restricted to Dirac fermions minimally coupled to electromagnetic fields. Non-minimal coupling cases such as neutrons and spatially spread cases such as fermions dressed with renormalization and quantum radiative corrections are not involved.</p>
</sec>
<sec id="SEC2.2"><title>2.2. Fermions with <inline-formula><tex-math notation="LaTeX" id="ImEquation86"><![CDATA[$\mathrm{g}\ne2, e=0$]]></tex-math></inline-formula> (Type B)</title> 
<p>Neutrons and neutral atoms correspond to the case of form factors of <inline-formula><tex-math notation="LaTeX" id="ImEquation87"><![CDATA[$F_1(0) = 0$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation88"><![CDATA[$F_2(0) \ne 0$]]></tex-math></inline-formula> since <inline-formula><tex-math notation="LaTeX" id="ImEquation89"><![CDATA[$e=0$]]></tex-math></inline-formula>. Since they cannot minimally couple to the electromagnetic field, their covariant Dirac equation contains the following non-minimal coupling between the magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation90"><![CDATA[$\mu_{\rm m}$]]></tex-math></inline-formula> and the electromagnetic tensor <inline-formula><tex-math notation="LaTeX" id="ImEquation91"><![CDATA[$F_{\mu\nu}$]]></tex-math></inline-formula>:
<disp-formula id="pty066-M-22"><label>(22)</label><tex-math notation="LaTeX" id="Equation22"><![CDATA[
\begin{eqnarray}
\left(\gamma^\mu \,i\, \nabla_\mu - \epsilon^{-1} m
- \epsilon\,\frac{1}{2}{\mu}_{\rm m}\,\sigma^{\mu\nu} F_{\mu\nu} \right) \Psi=0
,
\end{eqnarray}]]></tex-math></disp-formula>
  where <inline-formula><tex-math notation="LaTeX" id="ImEquation92"><![CDATA[$\sigma^{\mu\nu}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation93"><![CDATA[$F_{\mu\nu}$]]></tex-math></inline-formula> are second-rank covariant tensors (see <xref ref-type="sec" rid="APP1">Appendix A</xref>). The Hamiltonian is calculated as
<disp-formula id="pty066-M-23"><label>(23)</label><tex-math notation="LaTeX" id="Equation23"><![CDATA[
\begin{eqnarray}
{\cal H}^{^{\rm B}}
&=& - i\, \Gamma_0 + \frac{1}{g^{00}}\gamma^0\gamma^i (-i\,\partial_i - i\,\Gamma_i)
+ \epsilon^{-1}\frac{1}{g^{00}}\gamma^0 m \left(1+\epsilon^2\,\frac{\mu_{\rm m}}{2m}\,\sigma^{\mu\nu} F_{\mu\nu}\right)\!.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>In the same manner as described in the previous section, the large component in the Schwarzschild spacetime <inline-formula><tex-math notation="LaTeX" id="ImEquation94"><![CDATA[${\cal H}_{+}^{^{\rm B}}$]]></tex-math></inline-formula> is obtained as
<disp-formula id="pty066-M-24"><label>(24)</label><tex-math notation="LaTeX" id="Equation24"><![CDATA[
\begin{eqnarray}
{\cal H}_{+}^{^{\rm B}}
&=& \epsilon^{-2} m \nonumber\\
& &
+ \frac{(\boldsymbol{\sigma} \cdot \boldsymbol{p})^2}{2m}
+ m\phi
- \boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B}
\nonumber\\
& & + \epsilon^2 \Bigg(
\frac{m\phi^2}{2}
- \frac{(\boldsymbol{\sigma} \cdot \boldsymbol{p})^4}{8m^3}
- 3\,\phi\,\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B}
+ \frac{3}{8m}\big\{\boldsymbol{\sigma} \cdot \boldsymbol{p}, \, \{\boldsymbol{\sigma} \cdot \boldsymbol{p},\,\phi\} \big\}
\nonumber\\
& & \hspace{10mm}
+ \frac{1}{8m^2}\Big\{\boldsymbol{\sigma} \cdot \boldsymbol{p}, \,\{\boldsymbol{\sigma} \cdot \boldsymbol{p},\,\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B} \} \Big\}
\Bigg)
\nonumber\\
& &
- \epsilon^2 \frac{i}{2m} \big[ \boldsymbol{\sigma} \cdot \boldsymbol{p},\,\, \boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{E}\big]
+ O(\epsilon^4)
\nonumber\\
&=&
\epsilon^{-2} m
\nonumber\\
& &
+ \frac{\boldsymbol{p}^2 }{2 m} + m\phi - \boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B}
\nonumber\\
& &
+ \epsilon^2
\Bigg(
\frac{m\phi^2}{2}
- \frac{\boldsymbol{p}^4 }{8 m^3}
- 3\phi \,\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B}
+ \frac{3}{2m}\, \boldsymbol{p} \cdot (\phi\,\boldsymbol{p}) - \frac{3\phi}{2m r^2} \boldsymbol{S} \cdot (\boldsymbol{r} \times \boldsymbol{p})
\nonumber\\
& &
\hspace{10mm}
- \frac{1}{8m^2} \nabla^2 (\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{B})
+ \frac{1}{4m^2}(\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{p})\, (\boldsymbol{B} \cdot \boldsymbol{p} )
\nonumber\\
& &
\hspace{10mm}
+ \frac{1}{4m^2} (\boldsymbol{B} \cdot \boldsymbol{p}) \, (\boldsymbol{\mu}_{\rm m} \cdot \boldsymbol{p})
+ \frac{\mu_{\rm m}}{4m^2} (\nabla \times \boldsymbol{B}) \times \boldsymbol{p}
\nonumber\\
& &
\hspace{10mm}
- \frac{1}{2m} {\mu}_{\rm m}\, \nabla \cdot \boldsymbol{E}
- \frac{1}{m} \boldsymbol{\mu}_{\rm m} \cdot (\boldsymbol{E} \times \boldsymbol{p})
- \frac{i}{2m}\, \boldsymbol{\mu}_{\rm m} \cdot (\nabla \times \boldsymbol{E})
\Bigg)
+ O(\epsilon^4)
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Accordingly, the effective magnetic moment is obtained as
<disp-formula id="pty066-M-25"><label>(25)</label><tex-math notation="LaTeX" id="Equation25"><![CDATA[
\begin{eqnarray}
\boldsymbol{\mu}_{\rm m}^{\rm eff} \simeq (1+ 3\,\epsilon^2\,\phi) \,\,\boldsymbol{\mu}_{\rm m}
\hspace{1cm}
(\mbox{for Type B}).
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>This result shows that the effective magnetic moment with the general relativistic influence up to the order of <inline-formula><tex-math notation="LaTeX" id="ImEquation95"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> has a common form for both Type A and Type B, in which the minimal and non-minimal coupling terms are employed in the Dirac equations (<xref ref-type="disp-formula" rid="pty066-M-16">16</xref>) and (<xref ref-type="disp-formula" rid="pty066-M-22">22</xref>), respectively.</p>
</sec>
<sec id="SEC2.3"><title>2.3. Fermions with <inline-formula><tex-math notation="LaTeX" id="ImEquation96"><![CDATA[$g\ne2$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation97"><![CDATA[$e\ne0$]]></tex-math></inline-formula> (Type C)</title> 
<p>Spread charged fermions have the form factors of <inline-formula><tex-math notation="LaTeX" id="ImEquation98"><![CDATA[$F_1(0) = 1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation99"><![CDATA[$F_2(0) \ne 0$]]></tex-math></inline-formula>. This case corresponds to electrons and muons dressed with quantum radiative corrections and also composite particles such as protons and atoms containing internal structure. The corresponding covariant Dirac equation contains both minimal coupling and non-minimal coupling terms as
<disp-formula id="pty066-M-26"><label>(26)</label><tex-math notation="LaTeX" id="Equation26"><![CDATA[
\begin{eqnarray}
\left(
\gamma^\mu \,i\,(\nabla_\mu + i\,e A_\mu) - \epsilon^{-1} m
- \epsilon\,\,\frac{{\rm g} - 2}{2}\frac{e}{4m}\sigma^{\mu\nu}F_{\mu\nu}
\right) \Psi=0
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The influence of radiative corrections and internal structure is contained in the last term in Eq. (<xref ref-type="disp-formula" rid="pty066-M-26">26</xref>), which can be confirmed to observe that it returns to Eq. (<xref ref-type="disp-formula" rid="pty066-M-16">16</xref>) to eliminate the last term corresponding to the spread of fermions. The Hamiltonian in <inline-formula><tex-math notation="LaTeX" id="ImEquation100"><![CDATA[$4\times4$]]></tex-math></inline-formula>-matrix format is calculated as
<disp-formula id="pty066-M-27"><label>(27)</label><tex-math notation="LaTeX" id="Equation27"><![CDATA[
\begin{eqnarray}
{\cal H}^{\rm C}
&=& - i\, \Gamma_0 + e A_0
+ \frac{1}{g^{00}}\gamma^0\gamma^i (-i\,\partial_i - i\,\Gamma_i + e A_i)
+ \epsilon^{-1} \frac{1}{g^{00}}\gamma^0 m
{\left(1 + \epsilon^2\,\frac{{\rm g} - 2}{2}\frac{e}{4m^2}\sigma^{\mu\nu}F_{\mu\nu}\right)}
,
\nonumber\\
\end{eqnarray}]]></tex-math></disp-formula>
and the large component in the <inline-formula><tex-math notation="LaTeX" id="ImEquation101"><![CDATA[$2\times2$]]></tex-math></inline-formula>-matrix format Hamiltonian <inline-formula><tex-math notation="LaTeX" id="ImEquation102"><![CDATA[${\cal H}_{+}^{^{\rm C}}$]]></tex-math></inline-formula> as
<disp-formula id="pty066-M-28"><label>(28)</label><tex-math notation="LaTeX" id="Equation28"><![CDATA[
\begin{eqnarray}
{\cal H}_{+}^{^{\rm C}}
&=& \epsilon^{-2} m
\nonumber\\
& & + \frac{(\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}))^2}{2m} + m\phi
+ e A_0
- {\it \Delta}\boldsymbol{\mu}_0 \cdot \boldsymbol{B}
\nonumber\\
& &
+ \epsilon^2
\Bigg(
\frac{m\phi^2}{2}
- \frac{(\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}))^4 }{8m^3}
- 3\,\phi\,{\it \Delta}\boldsymbol{\mu}_0 \cdot \boldsymbol{B}
\nonumber\\
& & \hspace{1cm}
+ \frac{3}{8m}
\Big\{ \boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}),\,
\big\{\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}),\,\phi \big\}\Big\}
- \frac{1}{8m^2}
\Big[ \boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}), \big[\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}), e A_0 \big] \Big]
\nonumber\\
& & \hspace{1cm}
+ \frac{1}{8m^2}\Big\{\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}), \,\{\boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}),\,{\it \Delta}\boldsymbol{\mu}_0 \cdot \boldsymbol{B} \} \Big\}
\Bigg)
\nonumber\\
& &
- \epsilon^2 \frac{i}{2m} \big[ \boldsymbol{\sigma} \cdot (\boldsymbol{p} - e\boldsymbol{A}), {\it \Delta}\boldsymbol{\mu}_0 \cdot \boldsymbol{E}\big]
+ O(\epsilon^4)
\nonumber\\
&=&
\epsilon^{-2} m
\nonumber\\
& &
+ \frac{ (\boldsymbol{p} - e \boldsymbol{A})^2 }{2m}
+ m\phi
+ e A_0
- ( \boldsymbol{\mu}_0 + \Delta\boldsymbol{\mu}_0 ) \cdot \boldsymbol{B}
\nonumber\\
& &
+ \epsilon^2
\Bigg(
\frac{m\phi^2}{2}
- \frac{(\boldsymbol{p} - e \boldsymbol{A})^4}{8 m^3}
- 3 \phi\,\, ( \boldsymbol{\mu}_0 + \Delta\boldsymbol{\mu}_0 ) \cdot \boldsymbol{B}
\nonumber\\
& & \hspace{1cm}
- \frac{(\boldsymbol{\mu}_0 \cdot \boldsymbol{B})^2}{2 m}
- \frac{1}{4 m^2} \nabla^2 (\boldsymbol{\mu}_0 \cdot \boldsymbol{B})
+ \frac{1}{2 m^2} (\boldsymbol{p} - e \boldsymbol{A}) \cdot \big( \boldsymbol{\mu}_0 \cdot \boldsymbol{B}\,\,(\boldsymbol{p} - e \boldsymbol{A}) \big)
\nonumber\\
& & \hspace{1cm}
+ {\frac{3}{2m}} (\boldsymbol{p} - e\boldsymbol{A}) \cdot \big( \phi(\boldsymbol{p} - e\boldsymbol{A}) \big)
- {\frac{3\phi}{2m r^2}} \boldsymbol{S} \cdot \big(\boldsymbol{r} \times (\boldsymbol{p} - e\boldsymbol{A})\big)
\nonumber\\
& &
\hspace{1cm}
- \frac{1}{8m^2} \nabla^2 (\Delta\boldsymbol{\mu}_0 \cdot \boldsymbol{B})
+ \frac{1}{4m^2} \big(\Delta\boldsymbol{\mu}_0 \cdot (\boldsymbol{p} - e\boldsymbol{A})\big)\,\, \big(\boldsymbol{B} \cdot (\boldsymbol{p} - e\boldsymbol{A})\big)
\nonumber\\
& &
\hspace{1cm}
+ \frac{1}{4m^2} \big(\boldsymbol{B} \cdot (\boldsymbol{p} - e\boldsymbol{A})\big) \,\, \big(\Delta\boldsymbol{\mu}_0 \cdot (\boldsymbol{p} - e\boldsymbol{A})\big)
+ \frac{\Delta\mu_0}{4m^2}\, (\nabla \times \boldsymbol{B}) \times \boldsymbol{p}
\nonumber\\
& &
\hspace{1cm}
- \frac{\mu_0}{4m} \,\nabla\cdot (\boldsymbol{E} + \frac{\partial \boldsymbol{A}}{\partial t} )
- \frac{1}{2m} \boldsymbol{\mu}_0 \cdot \big((\boldsymbol{E} + \frac{\partial \boldsymbol{A}}{\partial t}) \times (\boldsymbol{p} - e \boldsymbol{A})\big)
\nonumber\\
& &
\hspace{1cm}
- \frac{\Delta \mu_0}{2m} \nabla \cdot \boldsymbol{E}
- \frac{1}{m} \Delta\boldsymbol{\mu}_0 \cdot (\boldsymbol{E} \times \boldsymbol{p})
- \frac{i}{2m} \Delta\boldsymbol{\mu}_0 \cdot (\nabla \times \boldsymbol{E})
\Bigg)+ O(\epsilon^4)
,
\end{eqnarray}]]></tex-math></disp-formula>
where
<disp-formula id="pty066-M-29"><label>(29)</label><tex-math notation="LaTeX" id="Equation29"><![CDATA[
\begin{eqnarray}
\boldsymbol{\mu}_0 &\equiv& \frac{e}{m}\boldsymbol{S} = \frac{e}{2m}\boldsymbol{\sigma} \hspace{5mm}\mbox{(Bohr magneton)},\nonumber\\
\Delta\boldsymbol{\mu}_0 &\equiv& \boldsymbol{\mu}_m - \boldsymbol{\mu}_0 = \frac{{\rm g}-2}{2}\frac{e}{2m}\boldsymbol{\sigma}.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Accordingly, the effective magnetic moment is
<disp-formula id="pty066-M-30"><label>(30)</label><tex-math notation="LaTeX" id="Equation30"><![CDATA[
\begin{eqnarray}
\boldsymbol{\mu}_{\rm m}^{\rm eff}
&\simeq&
( 1 + 3\epsilon^2\phi) \,\, \boldsymbol{\mu}_{\rm m}
\hspace{1cm}
(\mbox{for Type C}).
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Equation (<xref ref-type="disp-formula" rid="pty066-M-30">30</xref>) can be understood as the general expression of the effective magnetic moment of fermions including leptons and hadrons. Although the g-factor of proton or atoms largely deviates from <inline-formula><tex-math notation="LaTeX" id="ImEquation103"><![CDATA[$\mathrm{g}=2$]]></tex-math></inline-formula>, the effective magnetic moment is given in the same form as Type A: Dirac particles with minimal coupling (<xref ref-type="disp-formula" rid="pty066-M-16">16</xref>) and Type B: neutral fermions with non-minimal coupling (<xref ref-type="disp-formula" rid="pty066-M-22">22</xref>).</p>
</sec>
<sec id="SEC2.4"><title>2.4. Effective anomalous magnetic moment in curved spacetime</title>
<p>The effective magnetic moment to be observed with the instrument fixed on the Earth&#x2019;s surface has been obtained as Eq. (<xref ref-type="disp-formula" rid="pty066-M-30">30</xref>) up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation104"><![CDATA[$O(1/c^2)$]]></tex-math></inline-formula> adopting the Schwarzschild metric as the background spacetime, common to leptons and hadrons of Type A, Type B, and Type C. The effective values of the gyromagnetic ratio (g-factor) and the anomalous magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation105"><![CDATA[${a\equiv{\rm g}/2 - 1={\mu}_{\rm m}/{\mu}_0-1}$]]></tex-math></inline-formula>, to be measured in ground-based experiments, can be defined as
<disp-formula id="pty066-M-31"><label>(31)</label><tex-math notation="LaTeX" id="Equation31"><![CDATA[
\begin{eqnarray}
{\rm g}^{\rm eff}
&\simeq& (1 + 3\,\epsilon^2\,\phi) \,\, {\rm g}
,
\nonumber\\
{a}^{\rm eff}
&\simeq& a \, + 3\, (1+a)\epsilon^2\phi
,
\end{eqnarray}]]></tex-math></disp-formula>
which implies that the <inline-formula><tex-math notation="LaTeX" id="ImEquation106"><![CDATA[${\rm g}$]]></tex-math></inline-formula>-factor and the anomalous magnetic moment depend on the gravitational field as a result of the influence of the curved spacetime in general relativity. For an example, the anomalous magnetic moment of fermions with <inline-formula><tex-math notation="LaTeX" id="ImEquation107"><![CDATA[$\mathrm{g}\simeq2$]]></tex-math></inline-formula> is evaluated as <inline-formula><tex-math notation="LaTeX" id="ImEquation108"><![CDATA[$|a^{\rm eff}| \simeq |3\epsilon^2\phi| \simeq 2.1\times 10^{-9}$]]></tex-math></inline-formula> on the Earth&#x2019;s surface.</p>
<p>This implies that the Earth&#x2019;s gravitational field induces an anomalous magnetic moment of <inline-formula><tex-math notation="LaTeX" id="ImEquation109"><![CDATA[$2.1 \times 10^{-9}$]]></tex-math></inline-formula> for electrons, muons, etc. regardless of the type of fermions, in addition to quantum radiative corrections.</p>
</sec>
<sec id="SEC2.5"><title>2.5. Parametric post-Newtonian (PPN) approximation</title>
<p>Here we consider the case in which the parametric post-Newtonian (PPN) approximation is adopted as the background spacetime:
<disp-formula id="pty066-M-32"><label>(32)</label><tex-math notation="LaTeX" id="Equation32"><![CDATA[
\begin{eqnarray}
ds^2 &=& \epsilon^{-2} \big( 1 +\epsilon^2 2\phi +\epsilon^4 2 \beta_* \phi^2 \big) dt^2
- (1-\epsilon^2 2\gamma_* \phi)\,\, (dx^2 + dy^2 + dz^2)
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>This is a primitive version of one of the possible extensions of general relativity [<xref ref-type="bibr" rid="B11">11</xref>], in which the parameter choice of <inline-formula><tex-math notation="LaTeX" id="ImEquation110"><![CDATA[$\beta_* \ne 1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation111"><![CDATA[$\gamma_* \ne 1$]]></tex-math></inline-formula> corresponds to the scalar&#x2013;tensor gravity model and the choice of <inline-formula><tex-math notation="LaTeX" id="ImEquation112"><![CDATA[$\beta_* = 1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation113"><![CDATA[$\gamma_* = 1$]]></tex-math></inline-formula> corresponds to standard general relativity. In the same manner as in the case of general relativity, the large-component Hamiltonian can be obtained via the covariant Dirac equation (<xref ref-type="disp-formula" rid="pty066-M-26">26</xref>) using the PPN metric as
<disp-formula id="pty066-M-33"><label>(33)</label><tex-math notation="LaTeX" id="Equation33"><![CDATA[
\begin{align}
{\cal H}^{\rm PPN}_{+}
&=
\epsilon^{-2} m
\nonumber\\
& \quad{} + e A_0
+ m\phi
+ \frac{ (\boldsymbol{p} - e \boldsymbol{A})^2 }{2m}
- ( \boldsymbol{\mu}_0 + \Delta\boldsymbol{\mu}_0 ) \cdot \boldsymbol{B}
\nonumber\\
& \quad{}
+ \epsilon^2
\Bigg(
(\beta_* - \frac{1}{2}) m\phi^2
- (1 + 2\gamma_*) \phi\,\, ( \boldsymbol{\mu}_0 + \Delta\boldsymbol{\mu}_0 ) \cdot \boldsymbol{B}
\nonumber\\
& \quad{} \hspace{1cm}
- \frac{(\boldsymbol{p} - e \boldsymbol{A})^4}{8 m^3}
- \frac{(\boldsymbol{\mu}_0 \cdot \boldsymbol{B})^2}{2 m}
- \frac{1}{4 m^2} \nabla^2 (\boldsymbol{\mu}_0 \cdot \boldsymbol{B})
\nonumber\\
& \quad{} \hspace{1cm}
+ \frac{1}{2 m^2} (\boldsymbol{p} - e \boldsymbol{A}) \cdot \big( \boldsymbol{\mu}_0 \cdot \boldsymbol{B}\,\,(\boldsymbol{p} - e \boldsymbol{A}) \big)
\nonumber\\
& \quad{} \hspace{1cm}
+ \frac{(1 + 2\gamma_*)}{2m} (\boldsymbol{p}-e\boldsymbol{A}) \cdot \big( \phi(\boldsymbol{p}-e\boldsymbol{A}) \big)
- \frac{(1 + 2\gamma_*)\phi}{2m r^2} \boldsymbol{S} \cdot \big( \boldsymbol{r} \times (\boldsymbol{p}-e\boldsymbol{A})\big)
\nonumber\\
& \quad{}
\hspace{1cm}
- \frac{1}{8m^2} \nabla^2 (\Delta\boldsymbol{\mu}_0 \cdot \boldsymbol{B})
+ \frac{\Delta\mu_0}{4m^2}\, (\nabla \times \boldsymbol{B}) \times \boldsymbol{p}
\nonumber\\
& \quad{}
\hspace{1cm}
+ \frac{1}{4m^2} \big(\Delta\boldsymbol{\mu}_0 \cdot (\boldsymbol{p} - e\boldsymbol{A})\big)\, \big(\boldsymbol{B} \cdot (\boldsymbol{p} - e\boldsymbol{A})\big)
+ \frac{1}{4m^2} \big(\boldsymbol{B} \cdot (\boldsymbol{p} - e\boldsymbol{A})\big) \, \big(\Delta\boldsymbol{\mu}_0 \cdot (\boldsymbol{p} - e\boldsymbol{A})\big)
\nonumber\\
& \quad{}
\hspace{1cm}
- \frac{\mu_0}{4m} \,\nabla\cdot (\boldsymbol{E} + \frac{\partial \boldsymbol{A}}{\partial t} )
- \frac{1}{2m} \boldsymbol{\mu}_0 \cdot \Big((\boldsymbol{E} + \frac{\partial \boldsymbol{A}}{\partial t}) \times (\boldsymbol{p} - e \boldsymbol{A})\Big)
\nonumber\\
& \quad{}
\hspace{1cm}
- \frac{\Delta \mu_0}{2m} \nabla \cdot \boldsymbol{E}
- \frac{1}{m} \Delta\boldsymbol{\mu}_0 \cdot (\boldsymbol{E} \times \boldsymbol{p})
- \frac{i}{2m} \Delta\boldsymbol{\mu}_0 \cdot (\nabla \times \boldsymbol{E})
\Bigg)
+ O(\epsilon^4)
.
\end{align}]]></tex-math></disp-formula></p>
<p>Accordingly, the effective magnetic moment can be derived, in the same manner, as
<disp-formula id="pty066-M-34"><label>(34)</label><tex-math notation="LaTeX" id="Equation34"><![CDATA[
\begin{eqnarray}
\boldsymbol{\mu}_{\rm m}^{\rm eff}
&\simeq&
\big( 1 + (1+2\gamma_*)\,\epsilon^2\phi \big) \,\, \boldsymbol{\mu}_{\rm m}
\hspace{1cm}
(\mbox{for PPN}).
 
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>We note that the parameter of the theoretical model of gravity <inline-formula><tex-math notation="LaTeX" id="ImEquation114"><![CDATA[$\gamma_*$]]></tex-math></inline-formula> explicitly appears in the gravitationally induced correction term. This implies that the precise determination of the magnetic moment in a specific local circumstance may determine the model parameter <inline-formula><tex-math notation="LaTeX" id="ImEquation115"><![CDATA[$\gamma_*$]]></tex-math></inline-formula>.</p>
</sec>
</sec>
<sec id="SEC3"><title>3. Penning trap experiment to measure electron <inline-formula><tex-math notation="LaTeX" id="ImEquation116"><![CDATA[$\mathrm{g}-2$]]></tex-math></inline-formula> in the Earth&#x2019;s gravity</title>
<p>The g-factor of leptons exactly equals 2 according to the Dirac theory (<inline-formula><tex-math notation="LaTeX" id="ImEquation117"><![CDATA[${\rm g}_l$]]></tex-math></inline-formula>=2, for <inline-formula><tex-math notation="LaTeX" id="ImEquation118"><![CDATA[$l = {e},\mu,\tau$]]></tex-math></inline-formula>), which is one of the important consequences of the Dirac equation for free fermions with an electric charge of <inline-formula><tex-math notation="LaTeX" id="ImEquation119"><![CDATA[$e$]]></tex-math></inline-formula> and a spin of <inline-formula><tex-math notation="LaTeX" id="ImEquation120"><![CDATA[$1/2$]]></tex-math></inline-formula> minimally coupling to an electromagnetic field. However, the real value of <inline-formula><tex-math notation="LaTeX" id="ImEquation121"><![CDATA[${\rm g}_l$]]></tex-math></inline-formula> deviates from 2 by a fraction of about <inline-formula><tex-math notation="LaTeX" id="ImEquation122"><![CDATA[$0.002$]]></tex-math></inline-formula> as a result of the quantum radiative corrections. This deviation is referred to as the anomalous magnetic moment, defined as
<disp-formula id="pty066-M-35"><label>(35)</label><tex-math notation="LaTeX" id="Equation35"><![CDATA[
\begin{equation}
a_{l}
\equiv \frac{{\rm g}_{l}}{2} -1
= \frac{\mu_{\rm m}^l}{\mu_{\rm B}^l} -1 \hspace{3mm} ( = 0.001\cdots )
,
\end{equation}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation123"><![CDATA[$\mu_{\rm m}^l$]]></tex-math></inline-formula> is the magnetic moment and <inline-formula><tex-math notation="LaTeX" id="ImEquation124"><![CDATA[${\mu_{\rm B}^l}$]]></tex-math></inline-formula> is the Bohr magneton. The Bohr magneton is defined as <inline-formula><tex-math notation="LaTeX" id="ImEquation125"><![CDATA[${\mu_{\rm B}^l}$]]></tex-math></inline-formula>=<inline-formula><tex-math notation="LaTeX" id="ImEquation126"><![CDATA[$(e/m_l)\,S$]]></tex-math></inline-formula> for the lepton mass <inline-formula><tex-math notation="LaTeX" id="ImEquation127"><![CDATA[$m_l$]]></tex-math></inline-formula> and spin <inline-formula><tex-math notation="LaTeX" id="ImEquation128"><![CDATA[$S$]]></tex-math></inline-formula>. For electrons, a comparison between the theoretical value, which was obtained with higher orders of quantum radiative corrections based on the standard model of elementary particles (SM) [<xref ref-type="bibr" rid="B12">12</xref>,<xref ref-type="bibr" rid="B15">15</xref>], and the experimental value, which was obtained in the precision measurement [<xref ref-type="bibr" rid="B13">13</xref>&#x2013;<xref ref-type="bibr" rid="B16">16</xref>], shows the extremely precise agreement up to the 12th digit:
<disp-formula id="pty066-M-36"><label>(36)</label><tex-math notation="LaTeX" id="Equation36"><![CDATA[
\begin{eqnarray}
a_{e{\rm (EXP)}} - a_{e{\rm (SM)}} = -0.91 \,\, (0.82) \times 10^{-12},
\end{eqnarray}]]></tex-math></disp-formula>
as shown in <xref ref-type="table" rid="T1">Table 1</xref>, which clearly validates the accuracy of the quantum electrodynamics (QED) employed in the standard model of elementary particles [<xref ref-type="bibr" rid="B12">12</xref>,<xref ref-type="bibr" rid="B15">15</xref>,<xref ref-type="bibr" rid="B16">16</xref>].</p>
<p><table-wrap id="T1" orientation="portrait" position="float"><label>Table 1.</label><caption><p>Comparison of the theoretical value and the experimental value of the anomalous magnetic moment of the electron.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left">&#x00A0;</th>
<th align="left">Anomalous magnetic moment (<inline-formula><tex-math notation="LaTeX" id="ImEquation129"><![CDATA[$a_e\equiv {\rm g}_e/2-1$]]></tex-math></inline-formula>)</th>
<th align="center">Ref.</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left">Experiment</td>
<td align="left"><inline-formula><tex-math notation="LaTeX" id="ImEquation130"><![CDATA[$a_{e{\rm (EXP)}} = 1159\,652\,180.73 \, (0.28) \times10^{-12}$]]></tex-math></inline-formula></td>
<td align="center">[<xref ref-type="bibr" rid="B13">13</xref>&#x2013;<xref ref-type="bibr" rid="B16">16</xref>]</td>
</tr>
<tr>
<td align="left">Theory</td>
<td align="left"><inline-formula><tex-math notation="LaTeX" id="ImEquation131"><![CDATA[$a_{e{\rm (SM)}} = 1159\,652\,181.643 (0.77) \times10^{-12}$]]></tex-math></inline-formula></td>
<td align="center">[<xref ref-type="bibr" rid="B12">12</xref>,<xref ref-type="bibr" rid="B15">15</xref>]</td>
</tr>
<tr>
<td align="left">Difference</td>
<td align="left"><inline-formula><tex-math notation="LaTeX" id="ImEquation132"><![CDATA[$a_{e{\rm (EXP)}} - a_{e{\rm (SM)}} = -0.91 (0.82) \times10^{-12}$]]></tex-math></inline-formula></td>
<td align="left">&#x00A0;</td>
</tr>
<tr>
<td align="left">(Deviation)</td>
<td align="left">(1.1 <inline-formula><tex-math notation="LaTeX" id="ImEquation133"><![CDATA[$\sigma$]]></tex-math></inline-formula>)</td>
<td align="left">&#x00A0;</td>
</tr>
</tbody>
</table>
</table-wrap></p>
<sec id="SEC3.1"><title>3.1. Interpretation of the experimental result of electron <inline-formula><tex-math notation="LaTeX" id="ImEquation134"><![CDATA[$\mathrm{g}-2$]]></tex-math></inline-formula></title>
<p>Below, we examine the interpretation of the experimental result of the Penning trap experiment [<xref ref-type="bibr" rid="B13">13</xref>&#x2013;<xref ref-type="bibr" rid="B15">15</xref>].</p>
<p>The cyclotron frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation135"><![CDATA[$\Omega_{\rm c}$]]></tex-math></inline-formula> and the spin precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation136"><![CDATA[$\Omega_{\rm s}$]]></tex-math></inline-formula> of an electron moving in a uniform magnetic field are measured in the Penning trap experiment [<xref ref-type="bibr" rid="B13">13</xref>&#x2013;<xref ref-type="bibr" rid="B15">15</xref>]. The anomalous magnetic moment of electrons <inline-formula><tex-math notation="LaTeX" id="ImEquation137"><![CDATA[$a_{e}$]]></tex-math></inline-formula> is regarded as being equal to the ratio of the anomalous spin precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation138"><![CDATA[$\Omega_{\rm a} \equiv \Omega_{\rm s} - \Omega_{\rm c}$]]></tex-math></inline-formula> to the cyclotron frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation139"><![CDATA[$\Omega_{\rm c}$]]></tex-math></inline-formula> as
<disp-formula id="pty066-M-37"><label>(37)</label><tex-math notation="LaTeX" id="Equation37"><![CDATA[
\begin{eqnarray}
a_{e}^{\rm conv} {}_{({\rm EXP})} = \frac{\Omega_{\rm a}}{\Omega_{\rm c}}= \frac{\Omega_{\rm s}}{\Omega_{\rm c}}-1
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The <inline-formula><tex-math notation="LaTeX" id="ImEquation140"><![CDATA[${a_{e}^{\rm conv}}_{\rm (EXP)}$]]></tex-math></inline-formula> has been treated as the electron anomalous magnetic moment in the conventional analysis without consideration of the gravitational effect. However, this experiment was carried out in the Earth&#x2019;s gravitational field. Thus, the experimentally measured frequencies in Eq. (<xref ref-type="disp-formula" rid="pty066-M-37">37</xref>) are the effective frequencies in a curved spacetime. Hereafter, we derive the effective values of the cyclotron frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation141"><![CDATA[$\Omega_{\rm c}^{\rm eff}$]]></tex-math></inline-formula> and the spin precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation142"><![CDATA[$\Omega_{\rm s}^{\rm eff}$]]></tex-math></inline-formula> in a curved spacetime.</p>
<p>General relativity requires that the translational motion of a free particle with an electric charge of <inline-formula><tex-math notation="LaTeX" id="ImEquation143"><![CDATA[$e$]]></tex-math></inline-formula> is described by the covariant equation (<xref ref-type="disp-formula" rid="pty066-M-2">2</xref>). Assuming that the Earth&#x2019;s gravitational field can be described with the Schwarzschild metric shown in Eq. (<xref ref-type="disp-formula" rid="pty066-M-3">3</xref>), the equation of translational motion of a free electron can be written as Eq. (<xref ref-type="disp-formula" rid="pty066-M-4">4</xref>). In general, the solution of this equation gives a spiral motion that is a combination of translational acceleration and rotational motion. The cyclotron frequency is a measure of the rotational motion. In other words, the cyclotron frequency corresponds to the frequency of the rotational motion of the constant-norm vector parallel to the velocity. We define the unit vector parallel to the 3D velocity <inline-formula><tex-math notation="LaTeX" id="ImEquation144"><![CDATA[$\boldsymbol{\beta}$]]></tex-math></inline-formula> as
<disp-formula id="pty066-M-38"><label>(38)</label><tex-math notation="LaTeX" id="Equation38"><![CDATA[
\begin{eqnarray}
\boldsymbol{\beta} \equiv \beta\,\,\boldsymbol{{\hat \beta}}, \,\,\,\,\,\,\,\,\,
\boldsymbol{{\hat \beta}} \cdot \boldsymbol{{\hat \beta}} = \frac{\boldsymbol{\beta}\cdot\boldsymbol{\beta}}{\beta^2} = 1
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Here we use <inline-formula><tex-math notation="LaTeX" id="ImEquation145"><![CDATA[$\boldsymbol{\beta} \cdot \nabla\phi=0$]]></tex-math></inline-formula> and assume that the effects of the gravitational gradient are negligibly small; the rotational component of the translational motion is given as
<disp-formula id="pty066-M-39"><label>(39)</label><tex-math notation="LaTeX" id="Equation39"><![CDATA[
\begin{eqnarray}
\frac{d \boldsymbol{{\hat \beta}} }{dt}
&=&
- \big( 1 + (2\gamma^2 + 1)\epsilon^2\phi \big)\,\,
\frac{e}{m} \left(\frac{\boldsymbol{B}}{\gamma} - \frac{\gamma}{\gamma^2 - 1} (\boldsymbol{\beta} \times \boldsymbol{E}) \right)
\times \boldsymbol{{\hat \beta}}
\nonumber\\
&=&
\boldsymbol{\Omega_{\rm c}}^{\rm eff} \,\times \, \boldsymbol{{\hat \beta}}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Using the cyclotron frequency in the flat spacetime <inline-formula><tex-math notation="LaTeX" id="ImEquation146"><![CDATA[$ \boldsymbol{\Omega}_{\rm c} = -e/m \big( \boldsymbol{B}/\gamma - \gamma/(\gamma^2 - 1)\, (\boldsymbol{\beta} \times \boldsymbol{E}) \big) , $]]></tex-math></inline-formula> the effective cyclotron frequency can be written as
<disp-formula id="pty066-M-40"><label>(40)</label><tex-math notation="LaTeX" id="Equation40"><![CDATA[
\begin{eqnarray}
\boldsymbol{\Omega}_{\rm c}^{\rm eff}
&=&
\,\,\,\, \big( 1 + (2\gamma^2 + 1)\epsilon^2\phi \big)\,\, \boldsymbol{\Omega}_{\rm c}
+ O(\epsilon^4)
,
\end{eqnarray}]]></tex-math></disp-formula>
which represents the cyclotron frequency in the Earth&#x2019;s gravitational field based on general relativity.</p>
<p>In the measurement of the electron <inline-formula><tex-math notation="LaTeX" id="ImEquation147"><![CDATA[${\rm g}_e - 2$]]></tex-math></inline-formula> experiment (i.e., the Penning trap experiment), the electron velocity is sufficiently small as <inline-formula><tex-math notation="LaTeX" id="ImEquation148"><![CDATA[$\beta\ll1$]]></tex-math></inline-formula> and the effective value of the electron&#x2019;s cyclotron frequency is given as
<disp-formula id="pty066-M-41"><label>(41)</label><tex-math notation="LaTeX" id="Equation41"><![CDATA[
\begin{eqnarray}
\boldsymbol{\Omega}_{\rm c}^{\rm eff} = (1 + 3\epsilon^2\phi) \,\, \boldsymbol{\Omega}_{\rm c}.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>We can interpret that the <inline-formula><tex-math notation="LaTeX" id="ImEquation149"><![CDATA[$\boldsymbol{\Omega}_{\rm c}^{\rm eff}$]]></tex-math></inline-formula> is used as the experimentally measured cyclotron frequency to be observed with the ground-based instrument.</p>
<p>Here we analyze the behavior of an electron spin in the Earth&#x2019;s gravitational field. The time evolution of an electron spin in the flat spacetime is described by the equation
<disp-formula id="pty066-M-42"><label>(42)</label><tex-math notation="LaTeX" id="Equation42"><![CDATA[
\begin{eqnarray}
\frac{d \boldsymbol{S}}{d t}
&=& \boldsymbol{\mu}_{\rm m} \times \boldsymbol{B}
= \boldsymbol{\Omega}_{\rm s} \times \boldsymbol{S}
,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation150"><![CDATA[$\boldsymbol{\Omega}_{\rm s}=-({\rm g}_{e}/2)(e/m)\,\boldsymbol{B}$]]></tex-math></inline-formula>. Rewriting this equation with the four vector, we obtain
<disp-formula id="pty066-M-43"><label>(43)</label><tex-math notation="LaTeX" id="Equation43"><![CDATA[
\begin{eqnarray}
\frac{d S^\mu}{d \tau}
&=&
\frac{{\rm g}_e}{2}\frac{e}{m}
\Bigg(
F^{\mu\nu} S_\nu
+ \frac{1}{c^2} u^\mu (S_\lambda F^{\lambda \nu} u_\nu )
\Bigg)
- \frac{1}{c^2} u^\mu (S_\lambda \frac{d u^\lambda}{d \tau} )
,
\end{eqnarray}]]></tex-math></disp-formula>
which is known as the Bargmann&#x2013;Michel&#x2013;Telegdi (BMT) equation [<xref ref-type="bibr" rid="B22">22</xref>]. In the same manner as in Eq. (<xref ref-type="disp-formula" rid="pty066-M-2">2</xref>), employing the covariant derivative in the BMT equation, we obtain the covariant equation of spin motion in general relativity as
<disp-formula id="pty066-M-44"><label>(44)</label><tex-math notation="LaTeX" id="Equation44"><![CDATA[
\begin{eqnarray}
\frac{D S^\mu}{d \tau}
&=&
\frac{{\rm g}_e}{2}\frac{e}{m}
\Bigg(
F^{\mu\nu} S_\nu
+ \frac{1}{c^2} u^\mu (S_\lambda F^{\lambda \nu} u_\nu )
\Bigg)
- \frac{1}{c^2} u^\mu (S_\lambda \frac{D u^\lambda}{d \tau} )
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Here we substitute the Schwarzschild metric defined by Eq. (<xref ref-type="disp-formula" rid="pty066-M-3">3</xref>) into the covariant BMT equation. In the same way as we have derived the cyclotron frequency based on the equation of translational motion of electrons, we consider the case where the gradient of the gravitational potential <inline-formula><tex-math notation="LaTeX" id="ImEquation151"><![CDATA[$\nabla\phi$]]></tex-math></inline-formula> is negligibly small, such as the case where the motion is limited on a horizontal plane. We define the four spin vector as <inline-formula><tex-math notation="LaTeX" id="ImEquation152"><![CDATA[$s^{\mu}=(S^0,\boldsymbol{S})$]]></tex-math></inline-formula>; the spatial component of the generalized BMT equation (<xref ref-type="disp-formula" rid="pty066-M-44">44</xref>) up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation153"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> can be written as
<disp-formula id="pty066-M-45"><label>(45)</label><tex-math notation="LaTeX" id="Equation45"><![CDATA[
\begin{eqnarray}
\frac{d \boldsymbol{S}}{d t}
=
\boldsymbol{S} \times \frac{\big(1+\epsilon^2\phi(2\gamma^2+1)\big)}{\gamma}\frac{{\rm g}_e}{2m}\boldsymbol{B}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>We can take <inline-formula><tex-math notation="LaTeX" id="ImEquation154"><![CDATA[$\beta\ll1$]]></tex-math></inline-formula> in the measurement of electron <inline-formula><tex-math notation="LaTeX" id="ImEquation155"><![CDATA[${\rm g}_e - 2$]]></tex-math></inline-formula> (i.e., the Penning trap experiment), which implies that the effective value of the electron spin precession frequency to be observed with the ground-based instrument is given as
<disp-formula id="pty066-M-46"><label>(46)</label><tex-math notation="LaTeX" id="Equation46"><![CDATA[
\begin{eqnarray}
\boldsymbol{\Omega}_{\rm s}^{\rm eff}
&=&
( 1+3\epsilon^2\phi )\,\, \boldsymbol{\Omega}_{\rm s}.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The experimental value of the anomalous magnetic moment obtained in the Penning trap experiment, denoted as <inline-formula><tex-math notation="LaTeX" id="ImEquation156"><![CDATA[$a_{e{\rm (EXP)}}^{\rm eff}$]]></tex-math></inline-formula> below, should be compared with the ratio of the kinematically calculated effective values <inline-formula><tex-math notation="LaTeX" id="ImEquation157"><![CDATA[$\Omega_{\rm c}^{\rm eff}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation158"><![CDATA[$\Omega_{\rm s}^{\rm eff}$]]></tex-math></inline-formula>. Substituting the effective values <inline-formula><tex-math notation="LaTeX" id="ImEquation159"><![CDATA[$\Omega_{\rm c}^{\rm eff}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation160"><![CDATA[$\Omega_{\rm s}^{\rm eff}$]]></tex-math></inline-formula> into Eq. (<xref ref-type="disp-formula" rid="pty066-M-37">37</xref>), we obtain
<disp-formula id="pty066-M-47"><label>(47)</label><tex-math notation="LaTeX" id="Equation47"><![CDATA[
\begin{eqnarray}
a_e^{\rm eff} {}_{({\rm EXP})}
&=&
\frac{\Omega_{\rm s}^{\rm eff}}{\Omega_{\rm c}^{\rm eff}}-1
= \frac{(1 + 3\epsilon^2\phi)\,\,{\Omega}_{\rm s}}{(1 + 3\epsilon^2\phi)\,\,{\Omega}_{\rm c}}-1
\nonumber\\
&=& \frac{\Omega_{\rm s}}{\Omega_{\rm c}}-1
\nonumber\\
&=& a_e^{\rm conv} {}_{({\rm EXP})}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Although both effective values <inline-formula><tex-math notation="LaTeX" id="ImEquation161"><![CDATA[$\Omega_{\rm c}^{\rm eff}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation162"><![CDATA[$\Omega_{\rm s}^{\rm eff}$]]></tex-math></inline-formula> are dependent on the Earth&#x2019;s gravitational field, the gravitational effects are canceled in their ratio and <inline-formula><tex-math notation="LaTeX" id="ImEquation163"><![CDATA[$a_{e{\rm (EXP)}}^{\rm eff}$]]></tex-math></inline-formula> coincides with <inline-formula><tex-math notation="LaTeX" id="ImEquation164"><![CDATA[$a_{e{\rm (EXP)}}^{\rm conv}$]]></tex-math></inline-formula>. It is reasonable to regard the high-precision agreement between the experimental and theoretical values as confirming the cancellation of gravitational effects. It also suggests that, if we take the ratio of the conventionally defined <inline-formula><tex-math notation="LaTeX" id="ImEquation165"><![CDATA[$a_e$]]></tex-math></inline-formula> in Eq. (<xref ref-type="disp-formula" rid="pty066-M-35">35</xref>) to the Bohr magneton, the &#x201C;true&#x201D; anomalous magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation166"><![CDATA[$a_e^{\rm eff} {}_{({\rm EXP})}$]]></tex-math></inline-formula> deviates from the &#x201C;conventional&#x201D; value <inline-formula><tex-math notation="LaTeX" id="ImEquation167"><![CDATA[$a_e^{\rm conv} {}_{({\rm EXP})}$]]></tex-math></inline-formula> as
<disp-formula id="pty066-M-48"><label>(48)</label><tex-math notation="LaTeX" id="Equation48"><![CDATA[
\begin{eqnarray}
a_e^{\rm eff} {}_{({\rm EXP})}
&\equiv& \frac{{\rm g}_e^{\rm eff}}{2} - 1
= \frac{\mu_{\rm m}^{\rm eff}}{\mu_{\rm B}}-1
= \frac{\Omega_{\rm s}^{\rm eff}}{\Omega_{\rm c}}-1
\nonumber\\
&=& (1 + 3\epsilon^2\phi)\,\frac{\Omega_{\rm s}}{\Omega_{\rm c}}-1
\nonumber\\
&\ne& \frac{\Omega_{\rm s}}{\Omega_{\rm c}}-1
= a_e^{\rm conv} {}_{({\rm EXP})}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Since the Earth&#x2019;s gravity induces an anomaly even for <inline-formula><tex-math notation="LaTeX" id="ImEquation168"><![CDATA[${\rm g}_e=2$]]></tex-math></inline-formula> additionally to radiative corrections, the comparison between the experimental value measured in the curved spacetime and the theoretical value calculated in the flat spacetime results in a difference of
<disp-formula id="pty066-M-49"><label>(49)</label><tex-math notation="LaTeX" id="Equation49"><![CDATA[
\begin{eqnarray}
\left\vert a_e^{\rm eff} {}_{({\rm EXP})} - a_e{}_{({\rm SM})} \right\vert
&\simeq& 3\epsilon^2|\phi|
\sim 2.1 \times 10^{-9}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Here we introduce the redefinition of the effective Bohr magneton including the general relativistic effects of the Earth&#x2019;s gravitational field as
<disp-formula id="pty066-M-50"><label>(50)</label><tex-math notation="LaTeX" id="Equation50"><![CDATA[
\begin{eqnarray}
\mu_{\rm B}^{\rm eff}
&\equiv&
\mu_{\rm B}^{\rm eff}(\phi)
\simeq
(1 + 3\epsilon^2\phi)\,\mu_{\rm B}.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>In this case, we can treat <inline-formula><tex-math notation="LaTeX" id="ImEquation169"><![CDATA[$a_{e{\rm(EXP)}}^{\rm eff}$]]></tex-math></inline-formula> as an invariant quantity independent of the gravitational field if we cancel the general relativistic effects by redefining the anomalous magnetic moment using the ratio of the effective magnetic moment to the effective Bohr magneton as
<disp-formula id="pty066-M-51"><label>(51)</label><tex-math notation="LaTeX" id="Equation51"><![CDATA[
\begin{eqnarray}
a_e^{\rm eff}{}_{({\rm EXP})}
&\equiv&
\frac{\mu_{\rm m}^{\rm eff}}{\mu_{\rm B}^{\rm eff}}-1
= \frac{\Omega_{\rm s}^{\rm eff}}{\Omega_{\rm c}^{\rm eff}}-1
\nonumber\\
&=& a_e^{\rm conv}{}_{({\rm EXP})},
\end{eqnarray}]]></tex-math></disp-formula>
instead of its ratio to the Bohr magneton in the flat spacetime. This redefinition restores the validity of the test of the standard model of elementary particles via the comparison of the theoretical value <inline-formula><tex-math notation="LaTeX" id="ImEquation170"><![CDATA[$a_{e{\rm (SM)}}$]]></tex-math></inline-formula> and the experimental value <inline-formula><tex-math notation="LaTeX" id="ImEquation171"><![CDATA[$a_e^{\rm eff} {}_{({\rm EXP})}$]]></tex-math></inline-formula>.</p>
</sec>
</sec>
<sec id="SEC4"><title>4. Storage ring experiment to measure muon <inline-formula><tex-math notation="LaTeX" id="ImEquation172"><![CDATA[$\mathrm{g}-2$]]></tex-math></inline-formula> in the Earth&#x2019;s gravity</title>
<p>In the storage ring experiments for <inline-formula><tex-math notation="LaTeX" id="ImEquation173"><![CDATA[${\rm g}_\mu - 2$]]></tex-math></inline-formula> such as CERN [<xref ref-type="bibr" rid="B17">17</xref>,<xref ref-type="bibr" rid="B18">18</xref>] and BNL E821 [<xref ref-type="bibr" rid="B19">19</xref>&#x2013;<xref ref-type="bibr" rid="B21">21</xref>], the muon anomalous magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation174"><![CDATA[$a_{\mu ({\rm EXP})}$]]></tex-math></inline-formula> was determined from the relation
<disp-formula id="pty066-M-52"><label>(52)</label><tex-math notation="LaTeX" id="Equation52"><![CDATA[
\begin{eqnarray}
a_{\mu ({\rm EXP})}
&\equiv& \frac{R}{\lambda - R}
,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation175"><![CDATA[$R \equiv {\Omega_{\rm a}}/{\Omega_p}$]]></tex-math></inline-formula> is the ratio of the anomalous spin precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation176"><![CDATA[$\Omega_{\rm a}$]]></tex-math></inline-formula> of muons to the proton spin precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation177"><![CDATA[$\Omega_p$]]></tex-math></inline-formula> measured using the nuclear magnetic resonance, and <inline-formula><tex-math notation="LaTeX" id="ImEquation178"><![CDATA[$\lambda\equiv {\mu_\mu }/{\mu_p}$]]></tex-math></inline-formula> is the ratio of the magnetic moment of muons to that of protons. This can be understood as <inline-formula><tex-math notation="LaTeX" id="ImEquation179"><![CDATA[$a_{\mu ({\rm EXP})}$]]></tex-math></inline-formula> being determined from the relation
<disp-formula id="pty066-M-53"><label>(53)</label><tex-math notation="LaTeX" id="Equation53"><![CDATA[
\begin{eqnarray}
a_{\mu ({\rm EXP})}
&=& \frac{\Omega_{\rm a}}{\Omega_{\rm L} - \Omega_{\rm a}}
,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation180"><![CDATA[$\Omega_{\rm L}$]]></tex-math></inline-formula> is the Larmor precession frequency of muons. This interpretation can be obtained using the classical kinematical framework and is the basis of the design of the experiment to measure the muon <inline-formula><tex-math notation="LaTeX" id="ImEquation181"><![CDATA[${\rm g}_\mu - 2$]]></tex-math></inline-formula> [<xref ref-type="bibr" rid="B15">15</xref>&#x2013;<xref ref-type="bibr" rid="B21">21</xref>]. We reconsider this interpretation including the spacetime curvature in the Earth&#x2019;s gravitational field on the basis of general relativity.</p>
<sec id="SEC4.1"><title>4.1. Spin motion in flat spacetime</title>
<p>Here we give an overview of the procedure to evaluate the experimental value of the anomalous magnetic moment using the equation of motion of spin in the flat spacetime (see, e.g., Chap. 11 in Ref. [<xref ref-type="bibr" rid="B23">23</xref>]). The BMT equation (<xref ref-type="disp-formula" rid="pty066-M-43">43</xref>) described in the laboratory frame (inertial frame) can be decomposed into a set of equations for time and spatial components as
<disp-formula id="pty066-M-54"><label>(54)</label><tex-math notation="LaTeX" id="Equation54"><![CDATA[
\begin{eqnarray}
\frac{d S^0}{d \tau}
&=&
F_0
+
\,\,\gamma^2 \, \left(\boldsymbol{S} \cdot \frac{d\boldsymbol{\beta}}{d\tau}\right)
\nonumber\\
\frac{d \boldsymbol{S}}{d \tau}
&=&
\boldsymbol{F}
+
\,\,\gamma^2 \boldsymbol{\beta}\, \left(\boldsymbol{S} \cdot \frac{d\boldsymbol{\beta}}{d\tau}\right)\!.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The 3D spin vector <inline-formula><tex-math notation="LaTeX" id="ImEquation182"><![CDATA[$\boldsymbol{S}$]]></tex-math></inline-formula> in the inertial frame and the 3D spin vector <inline-formula><tex-math notation="LaTeX" id="ImEquation183"><![CDATA[$\boldsymbol{s}$]]></tex-math></inline-formula> are related through the Lorentz transformation
<disp-formula id="pty066-M-55"><label>(55)</label><tex-math notation="LaTeX" id="Equation55"><![CDATA[
\begin{eqnarray}
\boldsymbol{s} &=& \boldsymbol{S} - \frac{\gamma}{\gamma+1}(\boldsymbol{\beta}\cdot\boldsymbol{S}) \boldsymbol{\beta} \nonumber\\
&=& \boldsymbol{S} - \frac{\gamma}{\gamma+1} S^0 \boldsymbol{\beta}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The equation of motion of the spin in the rest frame <inline-formula><tex-math notation="LaTeX" id="ImEquation184"><![CDATA[$\boldsymbol{s}$]]></tex-math></inline-formula> can be derived as
<disp-formula id="pty066-M-56"><label>(56)</label><tex-math notation="LaTeX" id="Equation56"><![CDATA[
\begin{eqnarray}
\frac{d\boldsymbol{s}}{d\tau}
&=& \left(
\boldsymbol{F}
- \frac{\gamma}{\gamma+1} F_0 \, \boldsymbol{\beta}
\right)
+ \frac{\gamma^2}{\gamma+1}\,\,\boldsymbol{s}
\times
\left(\boldsymbol{\beta} \times {\frac{d\boldsymbol{\beta}}{d\tau}}\right) \nonumber\\
&=& \boldsymbol{F}'
+ \frac{\gamma^2}{\gamma+1}\,\,\boldsymbol{s}
\times
\left(\boldsymbol{\beta} \times {\frac{d\boldsymbol{\beta}}{d\tau}}\right)\!.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Rewriting the proper time <inline-formula><tex-math notation="LaTeX" id="ImEquation185"><![CDATA[$\tau$]]></tex-math></inline-formula> with the time in the inertial frame <inline-formula><tex-math notation="LaTeX" id="ImEquation186"><![CDATA[$t$]]></tex-math></inline-formula>, we can reproduce the familiar equation
<disp-formula id="pty066-M-57"><label>(57)</label><tex-math notation="LaTeX" id="Equation57"><![CDATA[
\begin{eqnarray}
\frac{d\boldsymbol{s}}{dt}
&=&
\frac{\boldsymbol{F'}}{\gamma} + \frac{\gamma^2}{\gamma+1}\,\,\boldsymbol{s}
\times
\left(\boldsymbol{\beta} \times {\frac{d\boldsymbol{\beta}}{dt}}\right) \nonumber\\
&=&
\frac{{\rm g}_\mu}{2} \frac{e}{\gamma\,m}\,\, \boldsymbol{s} \times \boldsymbol{B'}
+ \frac{\gamma^2}{\gamma+1}\,\,\boldsymbol{s}
\times
\left(\boldsymbol{\beta} \times {\frac{d\boldsymbol{\beta}}{dt}}\right)\!.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>We note that the last term in the equation of motion of spin in the rest frame (<xref ref-type="disp-formula" rid="pty066-M-57">57</xref>) represents the relativistic term known as Thomas precession, which does not appear in the non-relativistic equation of motion of spin (<xref ref-type="disp-formula" rid="pty066-M-42">42</xref>). The electric and magnetic fields in the inertial frame <inline-formula><tex-math notation="LaTeX" id="ImEquation187"><![CDATA[$\boldsymbol{E}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation188"><![CDATA[$\boldsymbol{B}$]]></tex-math></inline-formula> are related to those in the rest frame <inline-formula><tex-math notation="LaTeX" id="ImEquation189"><![CDATA[$\boldsymbol{E}'$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation190"><![CDATA[$\boldsymbol{B}'$]]></tex-math></inline-formula> through the Lorentz transformation as
<disp-formula id="pty066-M-58"><label>(58)</label><tex-math notation="LaTeX" id="Equation58"><![CDATA[
\begin{eqnarray}
\boldsymbol{E'} &=& \gamma (\boldsymbol{E}+\boldsymbol{\beta} \times \boldsymbol{B})
- \frac{\gamma^2}{\gamma+1}\,\,(\boldsymbol{\beta} \cdot \boldsymbol{E})\,\boldsymbol{\beta}\,
\nonumber\\
\boldsymbol{B'} &=& \gamma (\boldsymbol{B}-\boldsymbol{\beta} \times \boldsymbol{E})
- \frac{\gamma^2}{\gamma+1}\,\,(\boldsymbol{\beta} \cdot \boldsymbol{B})\,\boldsymbol{\beta}.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The equation of motion of spin in the rest frame can be rewritten with the quantities in the inertial frame by substituting Eqs. (<xref ref-type="disp-formula" rid="pty066-M-6">6</xref>) and (<xref ref-type="disp-formula" rid="pty066-M-58">58</xref>), which results in
<disp-formula id="pty066-M-59"><label>(59)</label><tex-math notation="LaTeX" id="Equation59"><![CDATA[
\begin{eqnarray}
\frac{d\boldsymbol{s}}{dt}
&=& \frac{\rm g}{2} \frac{e}{m}\,\, \boldsymbol{s}
\times
\left(
\boldsymbol{B}-\boldsymbol{\beta} \times \boldsymbol{E}
- \frac{\gamma}{\gamma+1}\,\,(\boldsymbol{\beta} \cdot \boldsymbol{B})\,\boldsymbol{\beta}
\right)
\nonumber\\
& &
\hspace{1cm}
+ \frac{\gamma^2}{\gamma + 1}\,\,\boldsymbol{s}
\times
\left(
\boldsymbol{\beta}
\times
\frac{e}{m\gamma}\,\,
\big(
\boldsymbol{E}
+ \boldsymbol{\beta} \times \boldsymbol{B}
- (\boldsymbol{\beta} \cdot \boldsymbol{E})\,\boldsymbol{\beta}
\big)
\right)
\nonumber\\
&=& \boldsymbol{s}
\times
\frac{e}{m}
\left(
(\frac{{\rm g}}{2} - 1 + \frac{1}{\gamma})\boldsymbol{B}
- (\frac{{\rm g}}{2} - \frac{\gamma}{\gamma+1})\boldsymbol{\beta} \times \boldsymbol{E}
- (\frac{{\rm g}}{2} - 1)\frac{\gamma}{\gamma+1}\,(\boldsymbol{\beta} \cdot \boldsymbol{B}) \,\boldsymbol{\beta}\right)\nonumber\\
&=& \boldsymbol{\Omega}_{\rm s} \times \boldsymbol{s}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>By using the spin precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation191"><![CDATA[$\boldsymbol{\Omega}_{\rm s}$]]></tex-math></inline-formula>, the anomalous spin precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation192"><![CDATA[$\boldsymbol{\Omega}_{\rm a}$]]></tex-math></inline-formula> can be written as
<disp-formula id="pty066-M-60"><label>(60)</label><tex-math notation="LaTeX" id="Equation60"><![CDATA[
\begin{eqnarray}
\boldsymbol{\Omega}_{\rm a}
&=&
\boldsymbol{\Omega}_{\rm s}
- \boldsymbol{\Omega}_{\rm c}
\nonumber\\
&=&
- \frac{e}{m}
\left(
a_\mu\,\boldsymbol{B}
- (a_\mu - \frac{1}{\gamma^2 - 1})\,\boldsymbol{\beta} \times \boldsymbol{E}
- a_\mu \frac{\gamma}{\gamma + 1} \,(\boldsymbol{\beta} \cdot \boldsymbol{B}) \,\boldsymbol{\beta}
\right)\!,
\end{eqnarray}]]></tex-math></disp-formula>
with <inline-formula><tex-math notation="LaTeX" id="ImEquation193"><![CDATA[$a_\mu\equiv {\rm g}_\mu/2 - 1$]]></tex-math></inline-formula>. We assume that the motion is confined in the horizontal plane and the magnetic field is applied vertically. In addition, we assume that the electric term (the second term) can be canceled as <inline-formula><tex-math notation="LaTeX" id="ImEquation194"><![CDATA[$\big( a_\mu - 1/(\gamma^2 - 1) \big) = 0$]]></tex-math></inline-formula> by selecting a muon momentum of <inline-formula><tex-math notation="LaTeX" id="ImEquation195"><![CDATA[$\gamma\simeq$]]></tex-math></inline-formula> 29.3, which is referred to as the magic momentum. The anomalous magnetic moment of the muon in the flat spacetime, which is defined as <inline-formula><tex-math notation="LaTeX" id="ImEquation196"><![CDATA[$a_{\mu {\rm (EXP)}} \equiv -{ \Omega_{\rm a} }/{ (e\,B/m) } $]]></tex-math></inline-formula>, can be written as
<disp-formula id="pty066-M-61"><label>(61)</label><tex-math notation="LaTeX" id="Equation61"><![CDATA[
\begin{eqnarray}
a_{\mu ({\rm EXP})}
&\equiv& - \frac{ { \Omega}_{\rm a} }{ (\frac{e { B}}{m}) }
= \frac{ { \Omega}_{\rm a} }{ { \Omega}_{\rm L} - { \Omega}_{\rm a} }
= \frac{ { \Omega_{\rm a}}/{\Omega_p} }{ { \Omega_{\rm L}}/{\Omega_p} - {\Omega_{\rm a}}/{\Omega_p} }
\nonumber\\
&=& \frac{R}{\lambda - R}
,
\end{eqnarray}]]></tex-math></disp-formula>
by rewriting the magnetic field <inline-formula><tex-math notation="LaTeX" id="ImEquation197"><![CDATA[$\boldsymbol{B}$]]></tex-math></inline-formula> using the Larmor precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation198"><![CDATA[$\Omega_{\rm L}=-({\rm g}_\mu/2)(e/m)\,{B}=-(1+a_\mu)(e/m){B}$]]></tex-math></inline-formula>. Equation (<xref ref-type="disp-formula" rid="pty066-M-61">61</xref>) represents the experimental value in the conventional interpretation. The experimental values of the anomalous magnetic moment <inline-formula><tex-math notation="LaTeX" id="ImEquation199"><![CDATA[$a_{\mu ({\rm EXP})}$]]></tex-math></inline-formula> using Eqs. (<xref ref-type="disp-formula" rid="pty066-M-60">60</xref>) and (<xref ref-type="disp-formula" rid="pty066-M-61">61</xref>) are evaluated within the classical kinematical framework. The relations were applied to achieve high precision by employing the magic momentum to suppress the contribution of electric fields.</p>
</sec>
<sec id="SEC4.2"><title>4.2. Spin motion in curved spacetime</title>
<p>In this section, we apply the procedure of the evaluation of the experimental value of the anomalous magnetic moment in the flat spacetime to that in the curved spacetime.</p>
<p>The BMT equation in the curved spacetime in the inertial frame (<xref ref-type="disp-formula" rid="pty066-M-44">44</xref>) can be decomposed into the time and spatial components as
<disp-formula id="pty066-M-62"><label>(62)</label><tex-math notation="LaTeX" id="Equation62"><![CDATA[
\begin{align}
& \hspace{0mm} \frac{d S^0}{d \tau}
+ \epsilon^2 ( S^0 \boldsymbol{u} + u^0 \boldsymbol{S} ) \cdot \nabla\phi
\nonumber\\
& \hspace{1cm} =
F_0
+ ( 1 - 4\gamma^2\epsilon^2\phi )
\,\,\gamma^2 \, \left(\boldsymbol{S} \cdot \frac{d\boldsymbol{\beta}}{d\tau}\right)
+ O(\epsilon^4)
\nonumber\\
& \hspace{0mm} \frac{d \boldsymbol{S}}{d \tau}
+ \epsilon \,\, u^0 S^0 \, \nabla\phi
+ \epsilon^2 \big(
(\boldsymbol{u} \cdot \boldsymbol{S})\nabla\phi - (\boldsymbol{S} \cdot \nabla\phi)\boldsymbol{u} - (\boldsymbol{u} \cdot \nabla\phi)\boldsymbol{S}
\big) \nonumber\\
& \hspace{1cm} =
\boldsymbol{F}
+ ( 1 - 4\gamma^2\epsilon^2\phi )
\,\,\gamma^2 \boldsymbol{\beta}\, \left(\boldsymbol{S} \cdot \frac{d\boldsymbol{\beta}}
{d\tau}\right)
+ O(\epsilon^4).
\end{align}]]></tex-math></disp-formula></p>
<p>Considering that the contribution arising from the gradient of the gravitational potential <inline-formula><tex-math notation="LaTeX" id="ImEquation200"><![CDATA[$\nabla\phi$]]></tex-math></inline-formula> is sufficiently small to be ignored in the post-Newtonian order (see <xref ref-type="sec" rid="APP5">Appendix E</xref>), Eq. (<xref ref-type="disp-formula" rid="pty066-M-62">62</xref>) can be simplified as
<disp-formula id="pty066-M-63"><label>(63)</label><tex-math notation="LaTeX" id="Equation63"><![CDATA[
\begin{eqnarray}
\frac{d S^0}{d \tau}
&=&
F_0
+ ( 1 - 4\gamma^2\epsilon^2\phi)
\,\,\gamma^2 \, \left(\boldsymbol{S} \cdot \frac{d\boldsymbol{\beta}}{d\tau}\right)
+ O(\epsilon^4)
\nonumber\\
\frac{d \boldsymbol{S}}{d \tau}
&=&
\boldsymbol{F}
+ ( 1 - 4\gamma^2\epsilon^2\phi )
\,\,\gamma^2 \boldsymbol{\beta}\, \left(\boldsymbol{S} \cdot \frac{d\boldsymbol{\beta}}{d\tau}\right) + O(\epsilon^4)
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Following the procedure in the flat spacetime, we derive the relation between the spins in the inertial frame and the rest frame using the Lorentz transformation. The Lorentz transformation can be naturally defined in the flat spacetime since the Lorentz invariance is satisfied at any point. However, the Lorentz transformation cannot be defined in the general coordinate system (curved spacetime), which disturbs the applicability of the same procedure. We avoid the difficulty by employing the local inertial frame to satisfy the local Lorentz invariance (see <xref ref-type="sec" rid="APP3">Appendix C</xref>). The four vector in the general coordinate system <inline-formula><tex-math notation="LaTeX" id="ImEquation201"><![CDATA[$x^\mu = (x^0, \boldsymbol{x})$]]></tex-math></inline-formula> can be related to the four vector in the local inertial frame <inline-formula><tex-math notation="LaTeX" id="ImEquation202"><![CDATA[$x^{(a)} = (\tilde{x}^0, \tilde{\boldsymbol{x}})$]]></tex-math></inline-formula>, using the tetrad <inline-formula><tex-math notation="LaTeX" id="ImEquation203"><![CDATA[$e^{(a)}_\mu$]]></tex-math></inline-formula>, as
<disp-formula id="pty066-M-64"><label>(64)</label><tex-math notation="LaTeX" id="Equation64"><![CDATA[
\begin{eqnarray}
x^{(a)} &\equiv& e^{(a)}_{\mu} \,x^{\mu} = (\tilde{x}^{0}, \tilde{\boldsymbol{x}})
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Using the explicit expression of the tetrad up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation204"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> in the Schwarzschild metric (shown in Eq. (<xref ref-type="disp-formula" rid="pty066-MB-7">B.7</xref>)), the transformation from the general coordinate system into the local inertial frame can be simplified as
<disp-formula id="pty066-M-65"><label>(65)</label><tex-math notation="LaTeX" id="Equation65"><![CDATA[
\begin{eqnarray}
\tilde{S}^0 &=& (1+\epsilon^2\phi) \, {S}^0 \nonumber\\
\tilde{\boldsymbol{S}} &=& (1-\epsilon^2\phi) \, \boldsymbol{S} \nonumber\\
\tilde{\boldsymbol{\beta}} &=& (1-2\epsilon^2\phi) \,{\boldsymbol{\beta }} \nonumber\\
\tilde{ t} &=& (1+\epsilon^2\phi) \,\,t \nonumber\\
\tilde{ \gamma } &=& \big( 1-2\epsilon^2\phi\, (\gamma^2-1) \big) \,{{ \gamma}} \nonumber\\
\tilde{\boldsymbol{E}} &=& \boldsymbol{E} \nonumber\\
\tilde{\boldsymbol{B}} &=& (1+2\epsilon^2\phi) \boldsymbol{B} \nonumber\\
& & \cdots \,\,\,\, {\rm etc.}
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>We rewrite the BMT equation in the curved spacetime (<xref ref-type="disp-formula" rid="pty066-M-63">63</xref>) using the local inertial frame and obtain
<disp-formula id="pty066-M-66"><label>(66)</label><tex-math notation="LaTeX" id="Equation66"><![CDATA[
\begin{eqnarray}
&&
\frac{d \tilde{S}^0}{d \tau}
=
\tilde{F}_0
+
\,\tilde{\gamma}^2 \,\left(\tilde{\boldsymbol{S}} \cdot \frac{d \tilde{\boldsymbol{\beta}}}{d\tau}\right)
+ O(\epsilon^4)
\nonumber\\
&&
\frac{d \tilde{\boldsymbol{S}}}{d \tau}
=
\tilde{\boldsymbol{F}}
+ \tilde{\gamma}^2
\tilde{\boldsymbol{\beta}}\, \left( \tilde{\boldsymbol{S}} \cdot \frac{d \tilde{\boldsymbol{\beta}}}{d\tau}\right)
+ O(\epsilon^4)
,
\end{eqnarray}]]></tex-math></disp-formula>
which is an analogue to the relation in the flat spacetime. The local Lorentz transformation is obtained as
<disp-formula id="pty066-M-67"><label>(67)</label><tex-math notation="LaTeX" id="Equation67"><![CDATA[
\begin{eqnarray}
\tilde{\boldsymbol{s}}
&=& \tilde{\boldsymbol{S}} - \frac{\tilde{\gamma}}{\tilde{\gamma}+1} (\tilde{\boldsymbol{\beta}}\cdot \tilde{\boldsymbol{S}}) \tilde{\boldsymbol{\beta}} \nonumber\\
&=& \tilde{\boldsymbol{S}} - \frac{\tilde{\gamma}}{\tilde{\gamma}+1} \tilde{S}^0 \tilde{\boldsymbol{\beta}}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>In the same manner as the flat spacetime, using the BMT equation in the curved spacetime (<xref ref-type="disp-formula" rid="pty066-M-66">66</xref>) and the local Lorentz transformation (<xref ref-type="disp-formula" rid="pty066-M-67">67</xref>), the equation of motion of spin in the rest frame can be written as
<disp-formula id="pty066-M-68"><label>(68)</label><tex-math notation="LaTeX" id="Equation68"><![CDATA[
\begin{eqnarray}
\frac{d\tilde{\boldsymbol{s}}}{d\tau}
&=& \left(
\tilde{\boldsymbol{F}}
- \frac{\tilde{\gamma}}{\tilde{\gamma}+1} \tilde{F}_0 \, \tilde{\boldsymbol{\beta}}
\right)
+ \frac{\tilde{\gamma}^2}{\tilde{\gamma}+1}\,\,\tilde{\boldsymbol{s}}
\times
\left(\tilde{\boldsymbol{\beta}} \times {\frac{d \tilde{\boldsymbol{\beta}}}{d\tau}}\right) \nonumber\\
&=& \tilde{\boldsymbol{F}}'
+ \frac{\tilde{\gamma}^2}{\tilde{\gamma}+1}\,\,\tilde{\boldsymbol{s}}
\times
\left(\tilde{\boldsymbol{\beta}} \times {\frac{d \tilde{\boldsymbol{\beta}}}{d\tau}}\right)\!,
\end{eqnarray}]]></tex-math></disp-formula>
which is similar to the case of the flat spacetime. This relation gives a relation similar to the case of the flat spacetime:
<disp-formula id="pty066-M-69"><label>(69)</label><tex-math notation="LaTeX" id="Equation69"><![CDATA[
\begin{eqnarray}
\frac{d \tilde{\boldsymbol{s}}}{d \tilde{t}}
&=&
\frac{\tilde{\boldsymbol{F}}'}{\tilde{\gamma}}
+ \frac{\tilde{\gamma}^2}{\tilde{\gamma}+1}\,\,\tilde{\boldsymbol{s}}
\times
\left(\tilde{\boldsymbol{\beta}} \times {\frac{d \tilde{\boldsymbol{\beta}}}{d \tilde{t}}}\right) \nonumber\\
&=&
\frac{{\rm g}_\mu}{2} \frac{e}{\tilde{\gamma}\,m}\,\, \tilde{\boldsymbol{s}} \times \tilde{\boldsymbol{B'}}
+ \frac{\tilde{\gamma}^2}{\tilde{\gamma}+1}\,\,\tilde{\boldsymbol{s}}
\times
\left(\tilde{\boldsymbol{\beta}} \times {\frac{d \tilde{\boldsymbol{\beta}}}{d \tilde{t}}}\right)\!,
\end{eqnarray}]]></tex-math></disp-formula>
by replacing the proper time <inline-formula><tex-math notation="LaTeX" id="ImEquation205"><![CDATA[$\tau$]]></tex-math></inline-formula> with the time in the inertial frame <inline-formula><tex-math notation="LaTeX" id="ImEquation206"><![CDATA[$\tilde{t}$]]></tex-math></inline-formula>. Here we notice that the electric and magnetic fields in the rest frame <inline-formula><tex-math notation="LaTeX" id="ImEquation207"><![CDATA[$\tilde{\boldsymbol{E'}}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation208"><![CDATA[$\tilde{\boldsymbol{B'}}$]]></tex-math></inline-formula> are given through a local Lorentz transformation, using the electric ands magnetic fields in the inertial frame <inline-formula><tex-math notation="LaTeX" id="ImEquation209"><![CDATA[$\tilde{\boldsymbol{E}}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation210"><![CDATA[$\tilde{\boldsymbol{B}}$]]></tex-math></inline-formula>, as
<disp-formula id="pty066-M-70"><label>(70)</label><tex-math notation="LaTeX" id="Equation70"><![CDATA[
\begin{eqnarray}
\tilde{\boldsymbol{E'}} &=& \tilde{\gamma} (\tilde{\boldsymbol{E}}+\tilde{\boldsymbol{\beta}} \times \tilde{\boldsymbol{B}})
- \frac{\tilde{\gamma}^2}{\tilde{\gamma}+1}\,\,(\tilde{\boldsymbol{\beta}} \cdot \tilde{\boldsymbol{E}})\,\tilde{\boldsymbol{\beta}}
\nonumber\\
\tilde{\boldsymbol{B'}} &=& \tilde{\gamma} (\tilde{\boldsymbol{B}}-\tilde{\boldsymbol{\beta}} \times \tilde{\boldsymbol{E}})
- \frac{\tilde{\gamma}^2}{\tilde{\gamma}+1}\,\,(\tilde{\boldsymbol{\beta}} \cdot \tilde{\boldsymbol{B}})\,\tilde{\boldsymbol{\beta}}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Thus, the equation of motion of spin in the rest frame <inline-formula><tex-math notation="LaTeX" id="ImEquation211"><![CDATA[$\tilde{\boldsymbol{s}}$]]></tex-math></inline-formula> can be written using the electric and magnetic fields in the inertial frame <inline-formula><tex-math notation="LaTeX" id="ImEquation212"><![CDATA[$\tilde{\boldsymbol{E}}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation213"><![CDATA[$\tilde{\boldsymbol{B}}$]]></tex-math></inline-formula> as
<disp-formula id="pty066-M-71"><label>(71)</label><tex-math notation="LaTeX" id="Equation71"><![CDATA[
\begin{eqnarray}
\frac{d \tilde{\boldsymbol{s}}}{d \tilde{t}}
&=& \frac{\rm g_\mu}{2} \frac{e}{m}\,\, \tilde{\boldsymbol{s}}
\times
\left(
\tilde{\boldsymbol{B}}-\tilde{\boldsymbol{\beta}} \times \tilde{\boldsymbol{E}}
- \frac{\tilde{\gamma}}{\tilde{\gamma}+1}\,\,(\tilde{\boldsymbol{\beta}} \cdot \tilde{\boldsymbol{B}})\, \tilde{\boldsymbol{\beta}}
\right)
\nonumber\\
& &
\hspace{1cm}
+ \frac{\tilde{\gamma}^2}{\tilde{\gamma} + 1}\,\,\tilde{\boldsymbol{s}}
\times
\left(
\tilde{\boldsymbol{\beta}}
\times
\frac{e}{m\tilde{\gamma}}\,\,
\left(
\tilde{\boldsymbol{E}}
+ \tilde{\boldsymbol{\beta}} \times \tilde{\boldsymbol{B}}
- (\tilde{\boldsymbol{\beta}} \cdot \tilde{\boldsymbol{E}})\,\tilde{\boldsymbol{\beta}}
\right)
\right)
\nonumber\\
&=& \tilde{\boldsymbol{s}}
\times
\frac{e}{m}
\left(
(\frac{{\rm g_\mu}}{2} - 1 + \frac{1}{\tilde{\gamma}})\tilde{\boldsymbol{B}}
- (\frac{{\rm g_\mu}}{2} - \frac{\tilde{\gamma}}{\tilde{\gamma}+1})\tilde{\boldsymbol{\beta}} \times \tilde{\boldsymbol{E}}
- (\frac{{\rm g_\mu}}{2} - 1)\frac{\tilde{\gamma}}{\tilde{\gamma}+1} (\tilde{\boldsymbol{\beta}} \cdot \tilde{\boldsymbol{B}})\,\tilde{\boldsymbol{\beta}}
\right)\!.\quad
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Equation (<xref ref-type="disp-formula" rid="pty066-M-71">71</xref>) has the same form as the equation of spin motion (<xref ref-type="disp-formula" rid="pty066-M-59">59</xref>). We convert the electric and magnetic fields <inline-formula><tex-math notation="LaTeX" id="ImEquation214"><![CDATA[$\tilde{\boldsymbol{E}}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation215"><![CDATA[$\tilde{\boldsymbol{B}}$]]></tex-math></inline-formula> in the local inertial frame to those in the general coordinate system to consider the experimental values with the ground-based instrument.</p>
<p>The equation of motion of spin <inline-formula><tex-math notation="LaTeX" id="ImEquation216"><![CDATA[$\boldsymbol{s}$]]></tex-math></inline-formula> in the rest frame represented in the general coordinate system is obtained as
<disp-formula id="pty066-M-72"><label>(72)</label><tex-math notation="LaTeX" id="Equation72"><![CDATA[
\begin{eqnarray}
\frac{d \boldsymbol{s}}{d t}
&=& \boldsymbol{s}
\times
\frac{e}{m}
\Bigg(
(\frac{{\rm g_\mu}}{2} - 1 + \frac{1}{\tilde{\gamma}})
(1 + 3\epsilon^2\phi) \boldsymbol{B}
\nonumber\\
& & \hspace{12mm}
- (\frac{{\rm g_\mu}}{2} - \frac{\tilde{\gamma}}{\tilde{\gamma}+1})
(1 - \epsilon^2\phi) \boldsymbol{\beta} \times \boldsymbol{E}
\nonumber\\
& & \hspace{12mm}
- (\frac{{\rm g_\mu}}{2} - 1)\frac{\tilde{\gamma}}{\tilde{\gamma}+1}
(1 - \epsilon^2\phi) (\boldsymbol{\beta} \cdot \boldsymbol{B})\,\boldsymbol{\beta}
\Bigg)
\nonumber\\
&=& \boldsymbol{s}
\times
\frac{e}{m}
\Bigg(
( \frac{{\rm g_\mu}}{2} - 1 + \frac{1}{\gamma}[1+2\epsilon^2\phi(\gamma^2 - 1)] )
(1 + 3\epsilon^2\phi) \boldsymbol{B}
\nonumber\\
& & \hspace{12mm}
- ( \frac{{\rm g_\mu}}{2} - \frac{\gamma}{\gamma + 1}[1 - 2\epsilon^2\phi(\gamma - 1)] )
(1 - \epsilon^2\phi) \boldsymbol{\beta} \times \boldsymbol{E}
\nonumber\\
& & \hspace{12mm}
- ( \frac{{\rm g_\mu}}{2} - 1)\frac{\gamma}{\gamma + 1}
[1 - \epsilon^2\phi(2\gamma - 1)] (\boldsymbol{\beta} \cdot \boldsymbol{B})\,\boldsymbol{\beta}
\Bigg)
\nonumber\\
&=& \boldsymbol{\Omega}_{\rm s}^{\rm eff} \times \boldsymbol{s}.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Putting <inline-formula><tex-math notation="LaTeX" id="ImEquation217"><![CDATA[$a_\mu\equiv {\rm g}_\mu/2 - 1$]]></tex-math></inline-formula> and using the effective value of the spin precession frequency in the curved spacetime (<xref ref-type="disp-formula" rid="pty066-M-72">72</xref>) and the effective value of the cyclotron frequency in the curved spacetime (<xref ref-type="disp-formula" rid="pty066-M-40">40</xref>), the effective value of the anomalous spin precession frequency <inline-formula><tex-math notation="LaTeX" id="ImEquation218"><![CDATA[$\boldsymbol{\Omega}_{\rm a}^{\rm eff}$]]></tex-math></inline-formula> is obtained as
<disp-formula id="pty066-M-73"><label>(73)</label><tex-math notation="LaTeX" id="Equation73"><![CDATA[
\begin{eqnarray}
\boldsymbol{\Omega}_{\rm a}^{\rm eff}
&=&
\boldsymbol{\Omega}_{\rm s}^{\rm eff}
-\boldsymbol{\Omega}_{\rm c}^{\rm eff}
\nonumber\\
&=&
- \frac{e}{m}
\Bigg(
(1 + 3\epsilon^2\phi)\, a_\mu {\bf B}
\nonumber\\
& & \hspace{10mm}
- \Big[
a_\mu - \frac{1}{\gamma^2 - 1}
- \epsilon^2\phi \Big(
4 + a_\mu
+ \frac{3}{\gamma^2 - 1}
\Big)
\Big] \,\,
\boldsymbol{\beta} \times {\bf E} \,\,\,\,
\nonumber\\
& & \hspace{10mm}
- a_\mu\, \frac{\gamma}{\gamma + 1}
\big( 1 - \epsilon^2\phi(2\gamma - 1) \big)
(\boldsymbol{\beta} \cdot {\bf B})\,\boldsymbol{\beta}
\Bigg)
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>When the motion is confined on the horizontal plane, a magnetic field is applied vertically, and the momentum is chosen to cancel the contribution of the electric field, the effective value of the anomalous magnetic moment of the muon is interpreted as
<disp-formula id="pty066-M-74"><label>(74)</label><tex-math notation="LaTeX" id="Equation74"><![CDATA[
\begin{eqnarray}
a_{\mu ({\rm EXP})}^{\rm eff}
&=& \frac{ { \Omega}_{\rm a}^{\rm eff} }{ { \Omega}_{\rm L}^{\rm eff} - { \Omega}_{\rm a}^{\rm eff} }
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The Larmor precession frequency in the curved spacetime can be interpreted as a special solution of the BMT equation with <inline-formula><tex-math notation="LaTeX" id="ImEquation219"><![CDATA[$u^\mu$]]></tex-math></inline-formula>=<inline-formula><tex-math notation="LaTeX" id="ImEquation220"><![CDATA[$(u^0,\boldsymbol{0})$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation221"><![CDATA[$s^\mu$]]></tex-math></inline-formula>=<inline-formula><tex-math notation="LaTeX" id="ImEquation222"><![CDATA[$(0, \boldsymbol{s})$]]></tex-math></inline-formula>, which corresponds to the null translational motion, as
<disp-formula id="pty066-M-75"><label>(75)</label><tex-math notation="LaTeX" id="Equation75"><![CDATA[
\begin{eqnarray}
\boldsymbol{\Omega}_{\rm L}^{\rm eff}
= - (1+3\epsilon^2\phi)\,\, \frac{\rm g_\mu}{2} \frac{e}{m} \boldsymbol{B}
= (1+3\epsilon^2\phi)\,\, \boldsymbol{\Omega}_{\rm L}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Therefore, we obtain
<disp-formula id="pty066-M-76"><label>(76)</label><tex-math notation="LaTeX" id="Equation76"><![CDATA[
\begin{eqnarray}
a_{\mu ({\rm EXP})}^{\rm eff}
&=& \frac{ { \Omega}_{\rm a}^{\rm eff} }{ { \Omega}_{\rm L}^{\rm eff} - { \Omega}_{\rm a}^{\rm eff} } \nonumber\\
&=& \frac{ (1+3\epsilon^2\phi)\,\, {\Omega}_{\rm a} }
{ (1+3\epsilon^2\phi)\,\, {\Omega}_{\rm L} - (1+3\epsilon^2\phi)\,\, {\Omega}_{\rm a} } \nonumber\\
&=& \frac{ {\Omega}_{\rm a} }{ {\Omega}_{\rm L} - {\Omega}_{\rm a} } \nonumber\\
&=& a_{\mu ({\rm EXP})}
,
\end{eqnarray}]]></tex-math></disp-formula>
which show that the effective value of the anomalous magnetic moment in the curved spacetime is equal to that in the flat spacetime up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation223"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> of general relativity. Thus the effective values of cyclotron frequency, spin precession frequency, and Larmor precession frequency have been derived based on the primitive consideration based on a kinematical framework. However those effective values in curved spacetime are, respectively, different from those values in the flat spacetime, the gravitational contribution is canceled in the ratio (<xref ref-type="disp-formula" rid="pty066-M-76">76</xref>) and the anomalous magnetic moment in the curved spacetime coincides with the case in the flat spacetime.</p>
</sec>
<sec id="SEC4.3"><title>4.3. Gravitational influence on the storage ring experiment to measure muon <inline-formula><tex-math notation="LaTeX" id="ImEquation224"><![CDATA[$\mathrm{g}-2$]]></tex-math></inline-formula></title>
<p>We have found that the gravitational influence on the magnetic moment of muons in the storage ring experiment is canceled in the ratio based on the kinematical consideration up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation225"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> using the Schwarzschild metric if the muon momentum is chosen to satisfy
<disp-formula id="pty066-M-77"><label>(77)</label><tex-math notation="LaTeX" id="Equation77"><![CDATA[
\begin{eqnarray}
a_\mu - \frac{1}{\gamma^2 - 1} =0 ,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation226"><![CDATA[$a_\mu$]]></tex-math></inline-formula> is the experimental value of the anomalous magnetic moment. Consequently, the comparison between the experimental and theoretical values remains valid as before. However, Eq. (<xref ref-type="disp-formula" rid="pty066-M-73">73</xref>) implies that general relativity modifies the contribution of the <inline-formula><tex-math notation="LaTeX" id="ImEquation227"><![CDATA[$\boldsymbol{\beta} \times \boldsymbol{E}$]]></tex-math></inline-formula> term. The perfect cancellation of the <inline-formula><tex-math notation="LaTeX" id="ImEquation228"><![CDATA[$\boldsymbol{\beta} \times \boldsymbol{E}$]]></tex-math></inline-formula> term corresponds to
<disp-formula id="pty066-M-78"><label>(78)</label><tex-math notation="LaTeX" id="Equation78"><![CDATA[
\begin{eqnarray}
a_\mu - \frac{1}{\gamma^2 - 1} - \epsilon^2\phi \left( 4 + a_\mu + \frac{3}{\gamma^2 - 1} \right)
\,\,\,\,\sim\,\,\,\,
(a_\mu + 2.8\times 10^{-9}) - \frac{1}{\gamma^2 - 1}
= 0
,
\end{eqnarray}]]></tex-math></disp-formula>
which introduces an offset to the experimental value of the anomalous magnetic moment as large as <inline-formula><tex-math notation="LaTeX" id="ImEquation229"><![CDATA[$2.8\times 10^{-9}$]]></tex-math></inline-formula>. Therefore, if the <inline-formula><tex-math notation="LaTeX" id="ImEquation230"><![CDATA[$\boldsymbol{\beta} \times \boldsymbol{E}$]]></tex-math></inline-formula> term is perfectly canceled in the experiment, <inline-formula><tex-math notation="LaTeX" id="ImEquation231"><![CDATA[$a_\mu + 2.8\times 10^{-9}$]]></tex-math></inline-formula> is measured instead of <inline-formula><tex-math notation="LaTeX" id="ImEquation232"><![CDATA[$a_\mu$]]></tex-math></inline-formula> and the offset may be misidentified as a discrepancy between the theory and the experiment. The general relativistic offset may be clarified in future experiments with improved accuracy.</p>
</sec>
</sec>
<sec id="SEC5"><title>5. Conclusion</title>
<p>The magnetic moment of free fermions to be observed using instruments fixed on the Earth&#x2019;s surface has been examined by evaluating the effective magnetic moment including the influence of the gravitational field up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation233"><![CDATA[$O(1/c^2)$]]></tex-math></inline-formula>, adopting the Schwarzschild metric as the background spacetime. We found that the magnetic moments of fermions are influenced by the spacetime curvature as <inline-formula><tex-math notation="LaTeX" id="ImEquation234"><![CDATA[$\boldsymbol{\mu}_{\rm m}^{\rm eff}= (1+3\phi/c^2) \,\,\boldsymbol{\mu}_{\rm m} $]]></tex-math></inline-formula> for the cases of minimal coupling, non-minimal coupling, and a mixture of the two. In the same manner, the effective values of the gyromagnetic ratio and the anomalous magnetic moment depend on the Earth&#x2019;s gravitational field as <inline-formula><tex-math notation="LaTeX" id="ImEquation235"><![CDATA[${\rm g}^{\rm eff} \simeq (1 + 3\phi/c^2) \,\, {\rm g} $]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation236"><![CDATA[${a}^{\rm eff} \simeq a + 3 (1 + a) \,\,\phi/c^2$]]></tex-math></inline-formula>, which seems inconsistent with the precise agreement between the experimental values and the theoretical values calculated in the flat spacetime.</p>
<p>However, experimental values of the anomalous magnetic moment were obtained as a ratio of the spin precession frequency and the cyclotron frequency in the Penning trap and the storage ring experiments, and the gravitational influence is canceled to recover the consistency between the experiment and the theory. Proper treatment of the experimental offset of the anomalous magnetic moment measured in the storage ring method would be necessary to compare the theoretical and experimental values, taking into consideration the change of the cancellation condition of the electric field contribution.</p>
</sec>
</body>
<back>
<ack>
<title>Acknowledgements</title>
<p>This work is supported in part by a Grant-in-Aid for Science Research from Japan Society for the Promotion of Science 501100001691 (JSPS) (No. 17K05453 to T.F).</p>
</ack>
<sec><title>Funding</title>
<p>Open Access funding: SCOAP<inline-formula><tex-math notation="LaTeX" id="ImEquation237"><![CDATA[$^3$]]></tex-math></inline-formula>.</p>
</sec>
<app-group><app id="APP1"><title>Appendix A. <inline-formula><tex-math notation="LaTeX" id="ImEquation238"><![CDATA[$\boldsymbol{\gamma}$]]></tex-math></inline-formula>-matrix</title>
<sec><title/>
<p>The <inline-formula><tex-math notation="LaTeX" id="ImEquation239"><![CDATA[$\gamma$]]></tex-math></inline-formula>-matrices in a curved spacetime satisfy
<disp-formula id="pty066-MA-1"><label>(A.1)</label><tex-math notation="LaTeX" id="Equation79"><![CDATA[
\begin{eqnarray}
\{\gamma^\mu, \gamma^\nu\} &=& \gamma^\mu\gamma^\nu+\gamma^\nu\gamma^\mu=2 g^{\mu\nu} I_4
,
\nonumber\\
\{\gamma_\mu, \gamma_\nu\} &=& \gamma_\mu\gamma_\nu+\gamma_\nu\gamma_\mu=2 g_{\mu\nu} I_4
,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation240"><![CDATA[$I_4$]]></tex-math></inline-formula> is the <inline-formula><tex-math notation="LaTeX" id="ImEquation241"><![CDATA[$4\times4$]]></tex-math></inline-formula> unit matrix, <inline-formula><tex-math notation="LaTeX" id="ImEquation242"><![CDATA[$g_{\mu\nu}$]]></tex-math></inline-formula> the metric tensor of the curved spacetime. In the flat spacetime, the above relation leads to
<disp-formula id="pty066-MA-2"><label>(A.2)</label><tex-math notation="LaTeX" id="Equation80"><![CDATA[
\begin{eqnarray}
\{\gamma^{(\alpha)}, \gamma^{(\beta)}\} &=& \gamma^{(\alpha)}\gamma^{(\beta)}+\gamma^{(\beta)}\gamma^{(\alpha)}=2 \eta^{\alpha\beta} I_4, \nonumber\\
\{\gamma_{(\alpha)}, \gamma_{(\beta)}\} &=& \gamma_{(\alpha)}\gamma_{(\beta)}+\gamma_{(\beta)}\gamma_{(\alpha)}=2 \eta_{\alpha\beta} I_4,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation243"><![CDATA[$\eta_{\alpha\beta}$]]></tex-math></inline-formula> is the metric tensor of the flat spacetime (Minkowski metric). We define the tetrad <inline-formula><tex-math notation="LaTeX" id="ImEquation244"><![CDATA[$e^{\mu}_{(\alpha)}$]]></tex-math></inline-formula> satisfying the relation
<disp-formula id="pty066-MA-3"><label>(A.3)</label><tex-math notation="LaTeX" id="Equation81"><![CDATA[
\begin{equation}
g_{\mu\nu} \,\,e^\mu_{(\alpha)} e^\nu_{(\beta)} = \eta_{\alpha\beta} =
\left(
\begin{array}{cccc}
\epsilon^{-2} & 0 & 0 & 0 \\
0 & -1 & 0 & 0 \\
0 & 0 & -1 & 0 \\
0 & 0 & 0 & -1
\end{array}
\right)\!.\,\,\,\,\,\,
\end{equation}]]></tex-math></disp-formula></p>
<p>The relation between the <inline-formula><tex-math notation="LaTeX" id="ImEquation245"><![CDATA[$\gamma$]]></tex-math></inline-formula>-matrices in the curved spacetime and the flat spacetime can be written as
<disp-formula id="pty066-MA-4"><label>(A.4)</label><tex-math notation="LaTeX" id="Equation82"><![CDATA[
\begin{equation}
\gamma^{\mu}=\gamma^{(\alpha)} \,\,e^\mu_{(\alpha)} ,\,\,\,\,\,
\gamma_{\mu}=\gamma_{(\alpha)} \,\,e_\mu^{(\alpha)}.
 
\end{equation}]]></tex-math></disp-formula></p>
<p>Using the tetrad <inline-formula><tex-math notation="LaTeX" id="ImEquation246"><![CDATA[$e^{\mu}_{(\alpha)}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation247"><![CDATA[$\gamma$]]></tex-math></inline-formula>-matrices, the spin connection <inline-formula><tex-math notation="LaTeX" id="ImEquation248"><![CDATA[$\Gamma_\mu$]]></tex-math></inline-formula> can be written as
<disp-formula id="pty066-MA-5"><label>(A.5)</label><tex-math notation="LaTeX" id="Equation83"><![CDATA[
\begin{eqnarray}
{\Gamma}_\mu &=& \frac{1}{4}\, g^{\lambda\nu} \, e_\lambda^{(\alpha)} \, \nabla_\mu \, e_\nu^{(\beta)}
\,\,\frac{1}{2}\, \Big[ \gamma_{(\alpha)}, \gamma_{(\beta)} \Big]
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Here we show the explicit representation of <inline-formula><tex-math notation="LaTeX" id="ImEquation249"><![CDATA[$\gamma$]]></tex-math></inline-formula>-matrices in the flat spacetime using the standard Dirac representation as
<disp-formula id="pty066-MA-6"><label>(A.6)</label><tex-math notation="LaTeX" id="Equation84"><![CDATA[
\begin{eqnarray}
&&
\gamma^{(0)} = {\epsilon} \beta =
\epsilon
\left(
\begin{array}{cc}
I & 0 \\
0 & -I
\end{array}
\right)\!,\,\,\,
\gamma^{(i)} =
\left(
\begin{array}{cc}
0 & \sigma_i \\
- \sigma_i & 0
\end{array}
\right)\!, \nonumber\\
&&
\gamma_{(0)} = \frac{1}{\epsilon}
\left(
\begin{array}{cc}
I & 0 \\
0 & -I
\end{array}
\right)\!,\,\,\,\,\,
\gamma_{(i)} =
\left(
\begin{array}{cc}
0 & -\sigma_i \\
\sigma_i & 0
\end{array}
\right)\!,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation250"><![CDATA[$\sigma_i$]]></tex-math></inline-formula> is the Pauli matrix satisfying
<disp-formula id="pty066-MA-7"><label>(A.7)</label><tex-math notation="LaTeX" id="Equation85"><![CDATA[
\begin{equation}
\sigma_i \, \sigma_j = \delta_{ij} + i\, \epsilon_{ijk} \,\sigma_k
\end{equation}]]></tex-math></disp-formula>
and their explicit representation is
<disp-formula id="pty066-MA-8"><label>(A.8)</label><tex-math notation="LaTeX" id="Equation86"><![CDATA[
\begin{equation}
\sigma_1 =
\left(
\begin{array}{cc}
0 & 1 \\
1 & 0
\end{array}
\right)\!,\,\,\,\,
\sigma_2 =
\left(
\begin{array}{cc}
0 & -i \\
i & 0
\end{array}
\right)\!,\,\,\,\,
\sigma_3 =
\left(
\begin{array}{cc}
1 & 0 \\
0 & -1
\end{array}
\right)\!.\,\,\,\,\,\,\,
\end{equation}]]></tex-math></disp-formula></p>
<p>The <inline-formula><tex-math notation="LaTeX" id="ImEquation251"><![CDATA[$4\times4$]]></tex-math></inline-formula> <inline-formula><tex-math notation="LaTeX" id="ImEquation252"><![CDATA[$\rho_i$]]></tex-math></inline-formula> matrices are Dirac matrices satisfying
<disp-formula id="pty066-MA-9"><label>(A.9)</label><tex-math notation="LaTeX" id="Equation87"><![CDATA[
\begin{equation}
\rho_i \, \rho_j = \delta_{ij} + i\, \epsilon_{ijk} \,\rho_k
\end{equation}]]></tex-math></disp-formula>
and their explicit representation is
<disp-formula id="pty066-MA-10"><label>(A.10)</label><tex-math notation="LaTeX" id="Equation88"><![CDATA[
\begin{equation}
\rho_1 =
\left(
\begin{array}{cc}
0 & I \\
I & 0
\end{array}
\right)\!,\,\,\,\,
\rho_2 =
\left(
\begin{array}{cc}
0 & -iI \\
iI & 0
\end{array}
\right)\!,\,\,\,\,
\rho_3 =
\left(
\begin{array}{cc}
I & 0 \\
0 & -I
\end{array}
\right)\!,\,\,\,\,\,\,\,
\end{equation}]]></tex-math></disp-formula>
which leads to
<disp-formula id="pty066-MA-11"><label>(A.11)</label><tex-math notation="LaTeX" id="Equation89"><![CDATA[
\begin{equation}
{\bf \alpha}_i = \rho_1 \sigma_i =
\left(
\begin{array}{cc}
0 & \sigma_i \\
\sigma_i & 0
\end{array}
\right)\!, \,\,\,\,\,\,\,\,
{\bf \beta} = \rho_3 =
\left(
\begin{array}{cc}
I & 0 \\
0 & -I
\end{array}
\right)\!.
\end{equation}]]></tex-math></disp-formula></p>
<p>Here we define the second-rank tensor <inline-formula><tex-math notation="LaTeX" id="ImEquation253"><![CDATA[$\sigma^{(\alpha)(\beta)}$]]></tex-math></inline-formula> as
<disp-formula id="pty066-MA-12"><label>(A.12)</label><tex-math notation="LaTeX" id="Equation90"><![CDATA[
\begin{gather}
\sigma^{(\alpha)(\beta)} = \frac{i}{2}\,\big[\gamma^{(\alpha)}, \gamma^{(\beta)}\big] ,\,\,\,\,\,\,\,\,
\gamma^{(\alpha)} \gamma^{(\beta)} = \eta^{\alpha\beta}- i \,\sigma^{(\alpha)(\beta)} \\
\end{gather}]]></tex-math></disp-formula>
<disp-formula id="pty066-MA-13"><label>(A.13)</label><tex-math notation="LaTeX" id="Equation91"><![CDATA[
\begin{gather}
\sigma^{(0)(k)} = \frac{i}{2}\,\big[\gamma^{(0)}, \gamma^{(k)}\big] = i\epsilon\rho_1\sigma_k = i\epsilon\alpha_k \nonumber\\
\sigma^{(i)(j)} = \frac{i}{2}\,\big[\gamma^{(i)}, \gamma^{(j)}\big] = \epsilon^{ijk} \sigma_k\\
\end{gather}]]></tex-math></disp-formula>
<disp-formula id="pty066-MA-14"><label>(A.14)</label><tex-math notation="LaTeX" id="Equation92"><![CDATA[
\begin{gather}
\sigma_{(0)(k)} = \frac{i}{2}\,\big[\gamma_{(0)}, \gamma_{(k)}\big] = -i\epsilon^{-1}\rho_1\sigma_k = -i\epsilon^{-1}\alpha_k \nonumber\\
\sigma_{(i)(j)} = \frac{i}{2}\,\big[\gamma_{(i)}, \gamma_{(j)}\big] = \epsilon_{ijk} \sigma_k
.
\end{gather}]]></tex-math></disp-formula></p>
<p>Using the tetrad <inline-formula><tex-math notation="LaTeX" id="ImEquation254"><![CDATA[$e^{(\alpha)}_{\mu}$]]></tex-math></inline-formula>, the second-rank tensors <inline-formula><tex-math notation="LaTeX" id="ImEquation255"><![CDATA[$\sigma^{\mu\nu}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation256"><![CDATA[$F_{\mu\nu}$]]></tex-math></inline-formula> can be written as
<disp-formula id="pty066-MA-15"><label>(A.15)</label><tex-math notation="LaTeX" id="Equation93"><![CDATA[
\begin{eqnarray}
\sigma^{\mu\nu} &=& e^{\mu}_{(\alpha)} \, e^{\nu}_{(\beta)} \sigma^{(\alpha)(\beta)} \nonumber\\
F_{\mu\nu} &=& e_{\mu}^{(\alpha)} \, e_{\nu}^{(\beta)} F_{(\alpha)(\beta)} 
\end{eqnarray}]]></tex-math></disp-formula>
and
<disp-formula id="pty066-MA-16"><label>(A.16)</label><tex-math notation="LaTeX" id="Equation94"><![CDATA[
\begin{eqnarray}
\sigma^{\mu\nu} F_{\mu\nu} &=& \sigma^{(\alpha)(\beta)}\,F_{(\alpha)(\beta)}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
</sec></app>
<app id="APP2"><title>Appendix B. Geometrical values in the Schwarzschild metric</title>
<sec><title/>
<p>The Schwarzschild metric can be written as
<disp-formula id="pty066-MB-1"><label>(B.1)</label><tex-math notation="LaTeX" id="Equation95"><![CDATA[
\begin{eqnarray}
ds^2 &=& \left(1 - \epsilon^2\frac{2GM}{r}\right) c^2 dt^2
- \frac{dr^2 }{(1 - \epsilon^2\frac{2GM}{r})}
- r^2(d\theta^2+\sin^2 \theta\, d{\varphi}^2)
\end{eqnarray}]]></tex-math></disp-formula>
in spherical coordinates. This can be rewritten, using isotropic coordinates and expanding up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation257"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula>, as
<disp-formula id="pty066-MB-2"><label>(B.2)</label><tex-math notation="LaTeX" id="Equation96"><![CDATA[
\begin{eqnarray}
ds^2
&=& \epsilon^{-2} (1 + \epsilon^2 2\phi+\epsilon^4 2\phi^2) \, dt^2
- (1 - \epsilon^2 2\phi) (dx^2 + dy^2 + dz^2)
\nonumber\\
& & \hspace{1cm} + O(\epsilon^4)
,
\end{eqnarray}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation258"><![CDATA[$\phi = -GM/r$]]></tex-math></inline-formula> is the Earth&#x2019;s gravitational potential. When the metric tensor <inline-formula><tex-math notation="LaTeX" id="ImEquation259"><![CDATA[$g_{\mu\nu}$]]></tex-math></inline-formula> is uniquely defined, the Christoffel symbol <inline-formula><tex-math notation="LaTeX" id="ImEquation260"><![CDATA[$\Gamma^\lambda_{\mu\nu}$]]></tex-math></inline-formula> and the tetrad <inline-formula><tex-math notation="LaTeX" id="ImEquation261"><![CDATA[$e^{\mu}_{(\alpha)}$]]></tex-math></inline-formula> are also uniquely determined in the following.</p>
<sec id="SECB.1"><title>B.1. Metric tensor <inline-formula><tex-math notation="LaTeX" id="ImEquation262"><![CDATA[$g_{\mu\nu}$]]></tex-math></inline-formula></title>
<p>Using the definition and Eq. (<xref ref-type="disp-formula" rid="pty066-MB-2">B.2</xref>), the explicit expression of each component of the four metric tensor <inline-formula><tex-math notation="LaTeX" id="ImEquation263"><![CDATA[$g_{\mu\nu}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation264"><![CDATA[$g^{\mu\nu}$]]></tex-math></inline-formula> is given up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation265"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> as
<disp-formula id="pty066-MB-3"><label>(B.3)</label><tex-math notation="LaTeX" id="Equation97"><![CDATA[
\begin{eqnarray}
g_{00} &=& \epsilon^{-2}\bigl( 1 + \epsilon^2 2\phi + \epsilon^4 2\phi^2 \bigr) \nonumber\\
g_{01} &=& g_{02} = g_{03} = 0 \nonumber\\
g_{ij} &=& -(1-2\epsilon^2\phi)\,\, \delta_{ij} \nonumber \\
g^{00} &=& \epsilon^{2}\big( 1 - 2\epsilon^2\phi + \epsilon^4 2\phi^2 \big) \nonumber\\
g^{01} &=& g^{02} = g^{03} = 0 \nonumber\\
g^{ij} &=& -(1+2\epsilon^2\phi)\,\, \delta^{ij} .
\end{eqnarray}]]></tex-math></disp-formula></p>
</sec>
<sec id="SECB.2"><title>B.2. Christoffel symbol <inline-formula><tex-math notation="LaTeX" id="ImEquation266"><![CDATA[$\Gamma^\lambda_{\mu\nu}$]]></tex-math></inline-formula></title>
<p>The Christoffel symbols are defined as
<disp-formula id="pty066-MB-4"><label>(B.4)</label><tex-math notation="LaTeX" id="Equation98"><![CDATA[
\begin{equation}
\Gamma^\lambda_{\mu\nu} = \frac{1}{2}g^{\lambda\kappa}(g_{\kappa\mu;\nu}+g_{\kappa\nu;\mu}-g_{\mu\nu;\kappa})
.
\end{equation}]]></tex-math></disp-formula></p>
<p>Substituting Eq. (<xref ref-type="disp-formula" rid="pty066-MB-3">B.3</xref>), we obtain explicit expressions for the Christoffel symbols up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation267"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> as
<disp-formula id="pty066-MB-5"><label>(B.5)</label><tex-math notation="LaTeX" id="Equation99"><![CDATA[
\begin{align}
& \Gamma^0_{00}=0 \nonumber\\
& \Gamma^1_{00}=\phi_1\hspace{1.3cm}\Gamma^2_{00}=\phi_2\hspace{1.3cm}\Gamma^3_{00}=\phi_3 \nonumber\\
& \Gamma^0_{01}=\epsilon^2\phi_1\hspace{1.cm}\Gamma^0_{02}=\epsilon^2\phi_2\hspace{1.cm}\Gamma^0_{03}=\epsilon^2\phi_3 \nonumber\\
& \Gamma^i_{0j}=0 \nonumber\\
& \Gamma^0_{ij}=0 \nonumber\\
& \Gamma^1_{11}=-\epsilon^2\phi_1\hspace{0.7cm}\Gamma^2_{11}=\epsilon^2\phi_2\hspace{1.0cm}\Gamma^3_{11}=\epsilon^2\phi_3\ \nonumber\\
& \Gamma^1_{12}=-\epsilon^2\phi_2\hspace{0.7cm}\Gamma^2_{12}=-\epsilon^2\phi_1\hspace{0.7cm}\Gamma^3_{12}=0\nonumber\\
& \Gamma^1_{13}=-\epsilon^2\phi_3\hspace{0.7cm}\Gamma^2_{13}=0\hspace{1.5cm}\Gamma^3_{13}=-\epsilon^2\phi_1\nonumber\\
& \Gamma^1_{22}= \epsilon^2\phi_1\hspace{1.0cm}\Gamma^2_{22}=-\epsilon^2\phi_2\hspace{0.7cm}\Gamma^3_{22}=\epsilon^2\phi_3\nonumber\\
& \Gamma^1_{23}=0\hspace{1.5cm}\Gamma^2_{23}=-\epsilon^2\phi_3\hspace{0.7cm}\Gamma^3_{23}=-\epsilon^2\phi_2\nonumber\\
& \Gamma^1_{33}=\epsilon^2\phi_1\hspace{1.0cm}\Gamma^2_{33}=\epsilon^2\phi_2\hspace{1.0cm}\Gamma^3_{33}=-\epsilon^2\phi_3
,
\end{align}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation268"><![CDATA[$\phi_i=\partial_i\phi=-\frac{\phi}{r^2}x_i$]]></tex-math></inline-formula>.</p>
</sec>
<sec id="SECB.3"><title>B.3. Tetrad <inline-formula><tex-math notation="LaTeX" id="ImEquation269"><![CDATA[$e^\mu_{(\alpha)}$]]></tex-math></inline-formula></title>
<p>Using Eqs. (<xref ref-type="disp-formula" rid="pty066-MB-3">B.3</xref>) and (<xref ref-type="disp-formula" rid="pty066-MB-5">B.5</xref>), the tetrad <inline-formula><tex-math notation="LaTeX" id="ImEquation270"><![CDATA[$e^\mu_{(a)}$]]></tex-math></inline-formula> satisfying
<disp-formula id="pty066-MB-6"><label>(B.6)</label><tex-math notation="LaTeX" id="Equation100"><![CDATA[
\begin{equation}
g_{\mu\nu} \,\,e^\mu_{(a)} e^\nu_{(b)} = \eta_{(a)(b)} =
\left(
\begin{array}{cccc}
\epsilon^{-2} & 0 & 0 & 0 \\
0 & -1 & 0 & 0 \\
0 & 0 & -1 & 0 \\
0 & 0 & 0 & -1
\end{array}
\right)\,\,\,\,\,\,
\end{equation}]]></tex-math></disp-formula>
can be obtained up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation271"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> as
<disp-formula id="pty066-MB-7"><label>(B.7)</label><tex-math notation="LaTeX" id="Equation101"><![CDATA[
\begin{eqnarray}
e^0_{(0)} &=& 1 - \epsilon^2 \phi+\epsilon^4\frac{1}{2}\phi^2 \nonumber\\
e^i_{(0)} &=& e^0_{(i)} = 0 \nonumber\\
e^i_{(j)} &=& (1 + \epsilon^2 \phi) \delta^i_j \nonumber\\
e^{(0)}_0 &=& 1 + \epsilon^2 \phi +\epsilon^4\frac{1}{2}\phi^2 \nonumber\\
e^{(i)}_0 &=& e^{(0)}_i = 0 \nonumber\\
e^{(i)}_j &=& (1 - \epsilon^2 \phi) \delta^i_j \nonumber\\
e^{(0)0} &=& \epsilon^2 \left(1-\epsilon^2 \phi +\epsilon^4\frac{1}{2}\phi^2\right) \nonumber\\
e^{(0)i} &=& e^{(i)0} = 0 \nonumber\\
e^{(i)j} &=& - (1 + \epsilon^2 \phi) \delta^i_j.
\end{eqnarray}]]></tex-math></disp-formula></p>
</sec>
</sec>
</app>
<app id="APP3"><title>Appendix C. Local inertial frame</title>
<sec><title/>
<p>The tetrad enables us to define the local inertial frame in which the local Lorentz invariance (LLI) is satisfied. The four vectors in the local inertial frame <inline-formula><tex-math notation="LaTeX" id="ImEquation272"><![CDATA[$x^{(\mu)}$]]></tex-math></inline-formula> and in the general coordinate system <inline-formula><tex-math notation="LaTeX" id="ImEquation273"><![CDATA[$x^{\mu}$]]></tex-math></inline-formula> are related to each other using the tetrad as
<disp-formula id="pty066-MC-1"><label>(C.1)</label><tex-math notation="LaTeX" id="Equation102"><![CDATA[
\begin{eqnarray}
x^{(a)} &\equiv& e^{(a)}_{\mu} x^{\mu} = (\tilde{x}^{0}, \tilde{\boldsymbol{x}})
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Using the representation of the tetrad up to the post-Newton order <inline-formula><tex-math notation="LaTeX" id="ImEquation274"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula> shown in Eq. (<xref ref-type="disp-formula" rid="pty066-MB-7">B.7</xref>), the physical quantities relevant in this paper can be expressed as
<disp-formula id="pty066-MC-2"><label>(C.2)</label><tex-math notation="LaTeX" id="Equation103"><![CDATA[
\begin{eqnarray}
\tilde{S}^0 &=& (1+\epsilon^2\phi) \, {S}^0 \nonumber\\
\tilde{\boldsymbol{S}} &=& (1-\epsilon^2\phi) \, \boldsymbol{S} \nonumber\\
\tilde{\boldsymbol{\beta}} &=& \frac{ \tilde{\boldsymbol{u}} }{ \tilde{u}^0 }
= (1-2\epsilon^2\phi) \,{\boldsymbol{\beta }} \nonumber\\
\tilde{ t} &=& (1+\epsilon^2\phi) \,\,t \nonumber\\
\tilde{\boldsymbol{x}} &=& (1-\epsilon^2\phi) \,\,\boldsymbol{x} \nonumber\\
\tilde{ \gamma } &=& \big( 1-2\epsilon^2\phi\, (\gamma^2-1) \big) \,{{ \gamma}} . 
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>The conversion of a tensor of rank 2 (electromagnetic tensor) is given as
<disp-formula id="pty066-MC-3"><label>(C.3)</label><tex-math notation="LaTeX" id="Equation104"><![CDATA[
\begin{eqnarray}
F_{(a)(b)} &=& e_{(a)}^{\mu} e_{(a)}^{\nu} F_{\mu\nu} 
,
\end{eqnarray}]]></tex-math></disp-formula>
which results in
<disp-formula id="pty066-MC-4"><label>(C.4)</label><tex-math notation="LaTeX" id="Equation105"><![CDATA[
\begin{eqnarray}
\tilde{\boldsymbol{E}} &=& F_{(0)(\alpha)}
= e_{(0)}^{0} e_{(\alpha)}^{\alpha} F_{0\alpha}
= \boldsymbol{E}
\nonumber\\
\tilde{\boldsymbol{B}} &=& F_{(\alpha)(\beta)}
= e_{(\alpha)}^{\alpha} e_{(\beta)}^{\beta} F_{\mu\nu}
= (1+2\epsilon^2\phi) \boldsymbol{B}
.
\end{eqnarray}]]></tex-math></disp-formula></p>
</sec>
</app>
<app id="APP4"><title>Appendix D. Normalization of the Hamiltonian</title>
<sec><title/>
<p>The scalar product of a spinor function <inline-formula><tex-math notation="LaTeX" id="ImEquation275"><![CDATA[$\Psi$]]></tex-math></inline-formula> is defined as
<disp-formula id="pty066-MD-1"><label>(D.1)</label><tex-math notation="LaTeX" id="Equation106"><![CDATA[
\begin{equation}
\left\langle \Psi \,{|}\, \Psi \right\rangle = \int \Psi^{\dagger} \Psi\, \sqrt{h}\, d^3x
,
\end{equation}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation276"><![CDATA[$h_{ij}$]]></tex-math></inline-formula> are the spatial components of the 4D tensor <inline-formula><tex-math notation="LaTeX" id="ImEquation277"><![CDATA[$g_{\mu\nu}$]]></tex-math></inline-formula> : <inline-formula><tex-math notation="LaTeX" id="ImEquation278"><![CDATA[$h_{ij}=-g_{ij}$]]></tex-math></inline-formula>. The expectation value of the Hamiltonian is given as
<disp-formula id="pty066-MD-2"><label>(D.2)</label><tex-math notation="LaTeX" id="Equation107"><![CDATA[
\begin{equation}
\left\langle {\cal H} \right\rangle
=
\left\langle \Psi \,{|}\, {\cal H} \,{|}\, \Psi \right\rangle
=
\int \Psi^{\dagger}\, {\cal H} \,\Psi\, \sqrt{h}\, d^3x
.
\end{equation}]]></tex-math></disp-formula></p>
<p>Here we put
<disp-formula id="pty066-MD-3"><label>(D.3)</label><tex-math notation="LaTeX" id="Equation108"><![CDATA[
\begin{eqnarray}
\Psi' = h^{\frac{1}{4}}\, \Psi
,\,\,\,\,\, \,\,\,\,\,
{\cal H}' = h^{\frac{1}{4}}\, {\cal H} \, h^{-\frac{1}{4}}\,
\end{eqnarray}]]></tex-math></disp-formula>
to obtain
<disp-formula id="pty066-MD-4"><label>(D.4)</label><tex-math notation="LaTeX" id="Equation109"><![CDATA[
\begin{eqnarray}
\left\langle \Psi \,{|}\, {\cal H} \,{|}\, \Psi \right\rangle
&=& \int \Psi^{\dagger}\, {\cal H} \,\Psi\, \sqrt{h}\, d^3x \nonumber\\
&=& \int {\Psi'}^{\dagger}\, {\cal H}' \,{\Psi'}\, d^3x \nonumber\\
&=& \left\langle \Psi' \,{|}\, {\cal H}' \,{|}\, \Psi' \right\rangle
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>This corresponds to the normalization of the expectation value of the Hamiltonian into the expectation value in the local flat spacetime as <inline-formula><tex-math notation="LaTeX" id="ImEquation279"><![CDATA[${\cal H}'$]]></tex-math></inline-formula>, which justifies the description of <inline-formula><tex-math notation="LaTeX" id="ImEquation280"><![CDATA[$\epsilon$]]></tex-math></inline-formula>-dependent terms such as the post-Newtonian terms as a perturbation to the terms in the flat spacetime [<xref ref-type="bibr" rid="B5">5</xref>].</p>
</sec>
</app>
<app id="APP5"><title>Appendix E. Typical scale of the gradient terms of the gravitational field</title>
<sec><title/>
<p>In general the equation of translational motion of a charged particle in curved spacetime involving general relativistic effects can be written as in Eq. (<xref ref-type="disp-formula" rid="pty066-M-2">2</xref>). Using the Schwarzschild metric and considering up to the post-Newtonian order <inline-formula><tex-math notation="LaTeX" id="ImEquation281"><![CDATA[$O(\epsilon^2)$]]></tex-math></inline-formula>, the equation can be written as in Eq. (<xref ref-type="disp-formula" rid="pty066-M-4">4</xref>). Here we estimate the magnitude of the contribution of the terms containing the gradient of the gravitational potential <inline-formula><tex-math notation="LaTeX" id="ImEquation282"><![CDATA[$\nabla\phi$]]></tex-math></inline-formula> (last two terms in Eq. (<xref ref-type="disp-formula" rid="pty066-M-4">4</xref>)) in the precise measurement of <inline-formula><tex-math notation="LaTeX" id="ImEquation283"><![CDATA[${\rm g}_{\mu} - 2$]]></tex-math></inline-formula> in the storage ring apparatus.</p>
<p>The BNL E821 experiment [<xref ref-type="bibr" rid="B19">19</xref>&#x2013;<xref ref-type="bibr" rid="B21">21</xref>] has apparatus parameters of <inline-formula><tex-math notation="LaTeX" id="ImEquation284"><![CDATA[$\gamma$]]></tex-math></inline-formula> = 29.3, <inline-formula><tex-math notation="LaTeX" id="ImEquation285"><![CDATA[$B$]]></tex-math></inline-formula> = 1.45 [T], <inline-formula><tex-math notation="LaTeX" id="ImEquation286"><![CDATA[$r$]]></tex-math></inline-formula> = 7112 [mm] and <inline-formula><tex-math notation="LaTeX" id="ImEquation287"><![CDATA[$\Omega_{\rm c} = 42 $]]></tex-math></inline-formula> [MHz], and the magnitudes of the electromagnetic interaction terms are estimated as
<disp-formula id="pty066-ME-1"><label>(E.1)</label><tex-math notation="LaTeX" id="Equation110"><![CDATA[
\begin{eqnarray}
f_{\rm EM}
&=& \Bigg|\, \big( 1 + (2\gamma^2 + 1) \epsilon^2\phi \big) \,\frac{e}{\gamma\,m} \,\boldsymbol{\beta} \times \boldsymbol{B} \,\Bigg| \nonumber\\
&=& \Big|\, \big( 1 + (2\gamma^2 + 1) \epsilon^2\phi \big) \,\, \Omega_{\rm c} \times \beta \,\Big| \nonumber\\
&\simeq& 4.2 \times 10^{7} \,\,\,{\rm [s^{-1}] }
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>On the other hand, the magnitudes of terms containing the gradient of the gravitational potential are
<disp-formula id="pty066-ME-2"><label>(E.2)</label><tex-math notation="LaTeX" id="Equation111"><![CDATA[
\begin{eqnarray}
f_{\nabla\phi}
&=& \Big|\, \epsilon (1+\beta^2)\nabla\phi \Big|\, \nonumber\\
&\simeq& 2 \times \frac{GM/R^2 }{c} \nonumber\\
&\simeq& 6.5 \times 10^{-8} \,\,{\rm [s^{-1}] } 
.
\end{eqnarray}]]></tex-math></disp-formula></p>
<p>Therefore, the relative magnitudes of the contributions of the gravity gradient and the electromagnetic interaction are
<disp-formula id="pty066-ME-3"><label>(E.3)</label><tex-math notation="LaTeX" id="Equation112"><![CDATA[
\begin{eqnarray}
\frac{f_{\nabla\phi}}{f_{\rm EM}} &\simeq& 1.5\times 10^{-15}
,
\end{eqnarray}]]></tex-math></disp-formula>
which shows that the gravity gradient contribution is <inline-formula><tex-math notation="LaTeX" id="ImEquation288"><![CDATA[$10^{-15}$]]></tex-math></inline-formula> times smaller than the electromagnetic interaction and is even <inline-formula><tex-math notation="LaTeX" id="ImEquation289"><![CDATA[$10^{-5}$]]></tex-math></inline-formula> times smaller than the post-Newtonian effects. Consequently, the contribution of the gravity gradient is negligible as long as the magnitude of the relevant contribution is <inline-formula><tex-math notation="LaTeX" id="ImEquation290"><![CDATA[$10^{-10}$]]></tex-math></inline-formula> times that of the main contribution.</p>
</sec></app></app-group>
<fn-group>
<title>Footnotes</title>
<fn id="FN1"><p><sup>1</sup> In this paper we use units <inline-formula><tex-math notation="LaTeX" id="ImEquation291"><![CDATA[$\hbar=c=1$]]></tex-math></inline-formula>.</p></fn>
<fn id="FN2"><p><sup>2</sup> The <inline-formula><tex-math notation="LaTeX" id="ImEquation292"><![CDATA[$TH\epsilon\mu$]]></tex-math></inline-formula> formalism with the Schwarzschild metric corresponds to <inline-formula><tex-math notation="LaTeX" id="ImEquation293"><![CDATA[${T = 1 + \epsilon^2 2\phi + \epsilon^4 2\phi^2}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation294"><![CDATA[${H = 1 - \epsilon^2 2\phi}$]]></tex-math></inline-formula>, which leads to the gravitational dependence of <inline-formula><tex-math notation="LaTeX" id="ImEquation295"><![CDATA[${T^{1/2}/H = 1 + 3\epsilon^2\phi}$]]></tex-math></inline-formula> for the magnetic moment [<xref ref-type="bibr" rid="B10">10</xref>].</p></fn>
</fn-group>
<ref-list>
<title>References</title>
<ref id="B1"><label>[1]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name name-style="western"><surname>Itzykson</surname><given-names>C.</given-names></string-name> and <string-name name-style="western"><surname>Zuber</surname><given-names>J. B.</given-names></string-name></person-group>, <source>Quantum Field Theory</source> (<publisher-name>McGraw-Hill</publisher-name>, <publisher-loc>New York</publisher-loc>, <year>1985</year>).</mixed-citation></ref>
<ref id="B2"><label>[2]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name name-style="western"><surname>Peskin</surname><given-names>M. E.</given-names></string-name> and <string-name name-style="western"><surname>Schroeder</surname><given-names>D. V.</given-names></string-name></person-group>, <source>An Introduction to Quantum Field Theory</source> (<publisher-name>Addison-Wesley</publisher-name>, <publisher-loc>New York</publisher-loc>, <year>1997</year>).</mixed-citation></ref>
<ref id="B3"><label>[3]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Anandan</surname><given-names>J.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>30</volume>, <fpage>1615</fpage> (<year>1984</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRevD.30.1615">http://dx.doi.org/10.1103/PhysRevD.30.1615</ext-link></comment>)</mixed-citation></ref>
<ref id="B4"><label>[4]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Parker</surname><given-names>L.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>22</volume>, <fpage>1922</fpage> (<year>1980</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRevD.22.1922">http://dx.doi.org/10.1103/PhysRevD.22.1922</ext-link></comment>)</mixed-citation></ref>
<ref id="B5"><label>[5]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Wajima</surname><given-names>S.</given-names></string-name>, <string-name name-style="western"><surname>Kasai</surname><given-names>M.</given-names></string-name>, and <string-name name-style="western"><surname>Futamase</surname><given-names>T.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>55</volume>, <fpage>1964</fpage> (<year>1997</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRevD.55.1964">http://dx.doi.org/10.1103/PhysRevD.55.1964</ext-link></comment>)</mixed-citation></ref>
<ref id="B6"><label>[6]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Foldy</surname><given-names>L. L.</given-names></string-name> and <string-name name-style="western"><surname>Wouthuysen</surname><given-names>S. A.</given-names></string-name></person-group>, <source>Phys. Rev.</source> <volume>78</volume>, <fpage>29</fpage> (<year>1950</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRev.78.29">http://dx.doi.org/10.1103/PhysRev.78.29</ext-link></comment>)</mixed-citation></ref>
<ref id="B7"><label>[7]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Tani</surname><given-names>S.</given-names></string-name></person-group>, <source>Prog. Theor. Phys.</source> <volume>6</volume>, <fpage>267</fpage> (<year>1951</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1143/ptp/6.3.267">http://dx.doi.org/10.1143/ptp/6.3.267</ext-link></comment>)</mixed-citation></ref>
<ref id="B8"><label>[8]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Case</surname><given-names>K. M.</given-names></string-name></person-group>, <source>Phys. Rev.</source> <volume>95</volume>, <fpage>1323</fpage> (<year>1954</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRev.95.1323">http://dx.doi.org/10.1103/PhysRev.95.1323</ext-link></comment>)</mixed-citation></ref>
<ref id="B9"><label>[9]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Lightman</surname><given-names>A. P.</given-names></string-name> and <string-name name-style="western"><surname>Lee</surname><given-names>D. L.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>8</volume>, <fpage>364</fpage> (<year>1973</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRevD.8.364">http://dx.doi.org/10.1103/PhysRevD.8.364</ext-link></comment>)</mixed-citation></ref>
<ref id="B10"><label>[10]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Will</surname><given-names>C. M.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>10</volume>, <fpage>2330</fpage> (<year>1974</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRevD.10.2330">http://dx.doi.org/10.1103/PhysRevD.10.2330</ext-link></comment>)</mixed-citation></ref>
<ref id="B11"><label>[11]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name name-style="western"><surname>Will</surname><given-names>C. M.</given-names></string-name></person-group>, <source>Theory and Experiment in Gravitational Physics</source> (<publisher-name>Cambridge University Press</publisher-name>, <publisher-loc>Cambridge, UK</publisher-loc>, <year>1993</year>).</mixed-citation></ref>
<ref id="B12"><label>[12]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Aoyama</surname><given-names>T.</given-names></string-name>, <string-name name-style="western"><surname>Hayakawa</surname><given-names>M.</given-names></string-name>, <string-name name-style="western"><surname>Kinoshita</surname><given-names>T.</given-names></string-name>, and <string-name name-style="western"><surname>Nio</surname><given-names>M.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>91</volume>, <fpage>033006</fpage> (<year>2015</year>). (<comment><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1103/PhysRevD.91.033006">https://doi.org/10.1103/PhysRevD.91.033006</ext-link></comment>)</mixed-citation></ref>
<ref id="B13"><label>[13]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Hanneke</surname><given-names>D.</given-names></string-name>, <string-name name-style="western"><surname>Fogwell</surname><given-names>S.</given-names></string-name>, and <string-name name-style="western"><surname>Gabrielse</surname><given-names>G.</given-names></string-name></person-group>, <source>Phys. Rev. Lett.</source> <volume>100</volume>, <fpage>120801</fpage> (<year>2008</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRevLett.100.120801">http://dx.doi.org/10.1103/PhysRevLett.100.120801</ext-link></comment>)</mixed-citation></ref>
<ref id="B14"><label>[14]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Hanneke</surname><given-names>D.</given-names></string-name>, <string-name name-style="western"><surname>Fogwell Hoogerheide</surname><given-names>S.</given-names></string-name>, and <string-name name-style="western"><surname>Gabrielse</surname><given-names>G.</given-names></string-name></person-group>, <source>Phys. Rev. A</source> <volume>83</volume>, <fpage>052122</fpage> (<year>2011</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRevA.83.052122">http://dx.doi.org/10.1103/PhysRevA.83.052122</ext-link></comment>)</mixed-citation></ref>
<ref id="B15"><label>[15]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Mohr</surname><given-names>P. J.</given-names></string-name>, <string-name name-style="western"><surname>Newell</surname><given-names>D. B.</given-names></string-name>, and <string-name name-style="western"><surname>Taylor</surname><given-names>B. N.</given-names></string-name></person-group>, <source>Rev. Mod. Phys.</source> <volume>88</volume>, <fpage>035009</fpage> (<year>2016</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/RevModPhys.88.035009">http://dx.doi.org/10.1103/RevModPhys.88.035009</ext-link></comment>)</mixed-citation></ref>
<ref id="B16"><label>[16]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Patrignani</surname><given-names>C.</given-names></string-name></person-group> <etal>et al</etal> [<collab>Particle Data Group</collab>], <source>Chin. Phys. C</source> <volume>40</volume>, <fpage>100001</fpage> (<year>2016</year>). (<comment><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1088/1674-1137/40/10/100001">https://doi.org/10.1088/1674-1137/40/10/100001</ext-link></comment>)</mixed-citation></ref>
<ref id="B17"><label>[17]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Bailey</surname><given-names>J.</given-names></string-name></person-group> <etal>et al</etal>, <source>Nucl. Phys. B</source> <volume>150</volume>, <fpage>1</fpage> (<year>1979</year>). (<comment><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/0550-3213(79)90292-X">https://doi.org/10.1016/0550-3213(79)90292-X</ext-link></comment>)</mixed-citation></ref>
<ref id="B18"><label>[18]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Bailey</surname><given-names>J.</given-names></string-name></person-group> <etal>et al</etal>, <source>Phys. Lett. B</source> <volume>68</volume>, <fpage>191</fpage> (<year>1977</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/0370-2693(77)90199-X">http://dx.doi.org/10.1016/0370-2693(77)90199-X</ext-link></comment>)</mixed-citation></ref>
<ref id="B19"><label>[19]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Bennett</surname><given-names>G. W.</given-names></string-name></person-group> <etal>et al</etal>, <source>Phys. Rev. D</source> <volume>73</volume>, <fpage>072003</fpage> (<year>2006</year>). (<comment><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1103/PhysRevD.73.072003">https://doi.org/10.1103/PhysRevD.73.072003</ext-link></comment>)</mixed-citation></ref>
<ref id="B20"><label>[20]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Miller</surname><given-names>J. P.</given-names></string-name>, <string-name name-style="western"><surname>de Rafael</surname><given-names>E.</given-names></string-name> and <string-name name-style="western"><surname>Roberts</surname><given-names>B. L.</given-names></string-name></person-group>, <source>Rep. Prog. Phys.</source> <volume>70</volume>, <fpage>795</fpage> (<year>2007</year>). (<comment><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1088/0034-4885/70/5/R03">https://doi.org/10.1088/0034-4885/70/5/R03</ext-link></comment>)</mixed-citation></ref>
<ref id="B21"><label>[21]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Paley</surname><given-names>J. M.</given-names></string-name></person-group>, Measurement of the anomalous magnetic moment of the negative muon to 0.7 parts per million, <source>Ph.D. Thesis</source>, Boston University (2004).</mixed-citation></ref>
<ref id="B22"><label>[22]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Bargmann</surname><given-names>V.</given-names></string-name>, <string-name name-style="western"><surname>Michel</surname><given-names>L.</given-names></string-name>, and <string-name name-style="western"><surname>Telegdi</surname><given-names>V. L.</given-names></string-name></person-group>, <source>Phys. Rev. Lett.</source> <volume>2</volume>, <fpage>435</fpage> (<year>1959</year>). (<comment><ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1103/PhysRevLett.2.435">http://dx.doi.org/10.1103/PhysRevLett.2.435</ext-link></comment>)</mixed-citation></ref>
<ref id="B23"><label>[23]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name name-style="western"><surname>Jackson</surname><given-names>J. D.</given-names></string-name></person-group>, <source>Classical Electrodynamics</source> (<publisher-loc>Wiley, New York</publisher-loc>, <year>1998</year>), <edition>3rd ed</edition>.</mixed-citation></ref>
</ref-list>
</back>
</article>