Nonvanishing magnetic field in pure electrostatic systems at rest in the curved space-time of Earth

TL;DR

This paper demonstrates that electrostatic systems at rest on Earth's curved spacetime can generate a nonvanishing magnetic field of purely geometric origin, impacting precision measurements like EDM searches.

Contribution

It reveals the existence of a geometric magnetic field in electrostatic systems on Earth, which cannot be shielded and affects sensitive EDM experiments.

Findings

Electrostatic systems exhibit a nonvanishing magnetic field due to spacetime curvature.

This magnetic field cannot be eliminated by magnetic shielding.

Implications for EDM experiments in electric storage rings and neutron EDM tests.

Abstract

Solutions of the Maxwell equations for electrostatic systems at rest on the rotating Earth's surface are shown to exhibit a nonvanishing magnetic field despite zero electric currents in the system. Such a field is of pure geometric origin, and in contrast to the conventional magnetic field of the Earth it can not be screened away by magnetic shielding. A practical significance of this magnetic field is that its EDM-like background signal would persist in the much discussed searches for the EDM of charged particles of magic energy trapped in all electric storage rings, as well as in the neutron EDM experiments.

Equations

$$F_{\mu \nu }={\partial }_{\mu }{A}_{\nu }-{\partial }_{\nu }{A}_{\mu }.$$

$${\epsilon }_{\mu \nu \lambda \rho }{\partial }_{\nu }{F}_{\lambda \rho }=0$$

$$\frac{1}{\sqrt{-g}}{\partial }_{\nu }\left(\sqrt{-g}{F}^{\nu \mu }\right)=4\pi J^{\mu },$$

$$g\equiv \det g_{\mu \nu },F^{\mu \nu }=g^{\mu \lambda }{g}^{\nu \rho }{F}_{\lambda \rho }.$$

$${\mathbf{E}}^{\alpha }=F^{\alpha 0}=e_{\mu }^{\alpha }{e}_{\nu }^0{F}^{\mu \nu }=e_i^{\alpha }{e}_0^0{F}^{i0}+e_i^{\alpha }{e}_j^0{F}^{ij},$$

$${\epsilon }_{\alpha \beta \gamma }{\mathbf{H}}^{\gamma }=-F^{\alpha \beta }=-e_{\mu }^{\alpha }{e}_{\nu }^{\beta }{F}^{\mu \nu },{\epsilon }_{123}=1.$$

$$\frac{D{X}^a}{ds}=\frac{d{X}^a}{ds}+\left({\gamma }_{b\mu }^a\frac{d{x}^{\mu }}{ds}\right)X^b=\frac{d{X}^a}{ds}+\left({\gamma }_{bc}^a{u}^c\right)X^b.$$

$$mc\frac{D{u}^{\mu }}{ds}=\frac{q}{c}{F}^{\mu \nu }{u}_{\nu }.$$

$$\frac{d\gamma }{dt}=\frac{q}{mc}\mathbf{E}\mathbf{\beta },$$

$$\frac{d(\gamma {\mathbf{\beta }}^{\alpha })}{dt}=\frac{q}{mc}{\left(\mathbf{E}+[\mathbf{\beta },\mathbf{H}]\right)}^{\alpha },$$

$$u^a=\gamma (1,\mathbf{\beta }),\mathbf{\beta }\equiv \mathbf{v}/c,ds=\frac{c}{\gamma }dt.$$

$$u^a=e_{\mu }^a\frac{d{x}^{\mu }}{ds}.$$

$$\displaystyle\frac{\mathop{\rm d}\nolimits{\bf S}}{\mathop{\rm d}\nolimits t}=\big{[}{\bm{\Omega}},\,{\bf S}\big{]},$$

$$\displaystyle{\bm{\Omega}}=\frac{q}{mc}\Bigg{\{}-\left(G+\frac{1}{\gamma}\right){\bf H}+\frac{\gamma G({\bf H}{\bm{\beta}})}{\gamma+1}{\bm{\beta}}+\left(G+\frac{1}{\gamma+1}\right)[{\bm{\beta}},\,{\bf E}]+$$

$$\displaystyle+\eta_{EDM}\left(-{\bf E}+\frac{\gamma}{\gamma+1}({\bf E}{\bm{\beta}}){\bm{\beta}}+[{\bf H},\,{\bm{\beta}}]\right)\Bigg{\}}.$$

$$\mu =(G+1)\frac{q\hslash }{2mc},d={\eta }_{EDM}\frac{q\hslash }{2mc}$$

$$\frac{d\gamma }{dt}=\frac{q}{mc}\mathbf{E}\mathbf{\beta }+c({\gamma }_{\alpha 0c}{u}^c){\mathbf{\beta }}^{\alpha },$$

$$\displaystyle\frac{\mathop{\rm d}\nolimits(\gamma{\bm{\beta}}^{\alpha})}{\mathop{\rm d}\nolimits t}=\frac{q}{mc}\left({\bf E}^{\alpha}+[{\bm{\beta}},\,{\bf H}]^{\alpha}\right)+c\Big{(}(\gamma_{\alpha 0c}u^{c})+(\gamma_{\alpha\beta c}u^{c}){\bm{\beta}}^{\beta}\Big{)}.$$

$$\displaystyle{\bm{\Omega}}^{\alpha}=-\frac{q}{mc}\bigg{\{}\big{(}G+1\big{)}{\bf H}+\eta_{EDM}{\bm{E}}\bigg{\}}^{\alpha}-\frac{c}{2}\varepsilon_{\alpha\beta\rho}\gamma_{\beta\rho 0}.$$

$$g_{0i}\ne 0.$$

$$(\partial /\partial x^0)J^{\mu }(x)=0,J^0\ne 0,J^i=0.$$

$$\sqrt{-g}{F}^{ij}={\epsilon }_{ijk}{\partial }_k\psi +\left({\partial }_i{\mathcal{A}}_j-{\partial }_j{\mathcal{A}}_i\right).$$

$$F^{ij}={\tilde{e}}_a^i{\tilde{e}}_b^j{F}^{ab}=-{\epsilon }_{\alpha \beta \gamma }{\tilde{e}}_{\alpha }^i{\tilde{e}}_{\beta }^j{\mathbf{H}}^{\gamma }=-\frac{e_0^0}{\sqrt{-g}}{\epsilon }_{ijk}{e}_k^{\alpha }{\mathbf{H}}^{\alpha }.$$

$${\mathbf{H}}^{\alpha }=-(e_0^0{)}^{-1}{\tilde{e}}_{\alpha }^i{\partial }_i\psi .$$

$$F_{i0}=e_i^a{e}_0^b{F}_{ab}=e_0^0{e}_i^{\alpha }{F}_{\alpha 0}=-e_0^0{e}_i^{\alpha }{\mathbf{E}}^{\alpha }={\partial }_i{A}_0,\text{mbox}or{\mathbf{E}}^{\alpha }=-\frac{{\tilde{e}}_{\alpha }^i}{e_0^0}{F}_{i0},$$

$$F_{ij}=e_i^a{e}_j^b{F}_{ab}=-{\epsilon }_{\alpha \beta \gamma }{e}_i^{\alpha }{e}_j^{\beta }{\mathbf{H}}^{\gamma }-\left(e_j^0{e}_i^{\alpha }-e_i^0{e}_j^{\alpha }\right){\mathbf{E}}^{\alpha }.$$

$$\displaystyle\partial_{i}\bigg{(}\sqrt{-g}\left(e^{0}_{0}\right)^{-2}g^{ij}\partial_{j}\psi\bigg{)}=-\varepsilon_{ijk}e^{0}_{0}e^{\alpha}_{i}{\bf E}^{\alpha}\partial_{k}\left(\frac{e^{0}_{j}}{e^{0}_{0}}\right).$$

$$F^{i0}={\tilde{e}}_a^i{\tilde{e}}_b^0{F}^{ab}={\tilde{e}}_0^0{\tilde{e}}_{\alpha }^i{\mathbf{E}}^{\alpha }-{\epsilon }_{\alpha \beta \gamma }{\tilde{e}}_{\alpha }^i{\tilde{e}}_{\beta }^0{\mathbf{H}}^{\gamma }=\frac{1}{e_0^0}{\tilde{e}}_{\alpha }^i{\mathbf{E}}^{\alpha }-\frac{1}{\sqrt{-g}{e}_0^0}{\epsilon }_{ijk}{e}_j^0{\partial }_k\psi .$$

$${\partial }_i\left(\frac{\sqrt{-g}}{e_0^0}{\tilde{e}}_{\alpha }^i{\mathbf{E}}^{\alpha }\right)-{\epsilon }_{ijk}{\partial }_i\left(\frac{e_j^0}{e_0^0}\right)\cdot {\partial }_k\psi =4\pi \sqrt{-g}{J}^0.$$

$$\frac{\sqrt{-g}{\tilde{e}}_{\alpha }^i}{e_0^0}=\left(1+\frac{r_g}{R}\right){\delta }_{\alpha }^i+\left(\frac{[\mathbf{\omega },\mathbf{R}{]}^2}{2c^2}{\delta }_{\alpha }^i-\frac{[\mathbf{\omega },\mathbf{R}{]}^i[\mathbf{\omega },\mathbf{R}{]}^{\alpha }}{2c^2}\right),$$

$$\frac{{\tilde{e}}_{\alpha }^i}{e_0^0}={\delta }_{\alpha }^i+\left(\frac{[\mathbf{\omega },\mathbf{R}{]}^2}{2c^2}{\delta }_{\alpha }^i-\frac{[\mathbf{\omega },\mathbf{R}{]}^i[\mathbf{\omega },\mathbf{R}{]}^{\alpha }}{2c^2}\right).$$

$$\frac{\sqrt{-g}}{{\left(e_0^0\right)}^2}{g}^{ij}=-\left(1+\frac{r_g}{R}\right){\delta }^{ij},e_0^0{e}_i^{\alpha }={\delta }_i^{\alpha }+O(\mathbf{\omega }\otimes \mathbf{\omega }),$$

$$\displaystyle\frac{e^{0}_{j}}{e^{0}_{0}}=-\Bigg{\{}1+\frac{r_{g}}{R}\left(2-\frac{IR^{2}_{\oplus}}{R^{2}}\right)\Bigg{\}}\frac{[{\bm{\omega}},\,{\bf R}]^{j}}{c},$$

$$\displaystyle\partial_{i}\bigg{\{}\left(1+\frac{r_{g}}{R}\right)\partial_{i}\psi\bigg{\}}$$

Full text

I Introduction

In this work we raise one generic issue with Maxwell equations in curved space-time. Specifically, we focus on the case of pure electrostatic systems in the noninertial reference frame with a globally time-independent metric tensor with $g_{0i}\neq 0$. Evidently, for an observer in a flat space-time and inertial frame the charges in motion do generate both electric and magnetic fields. The issue is what are the fields seen by an observer in the noninertial frame comoving with the electric charges, so that to this observer the electric current is a vanishing one. Will he see pure electric field or shall there be some trace of the noninertial motion in the form of a residual magnetic field? From the point of view of Einstein’s general relativity theory such a magnetic field if exists would be of geometric origin and will have a peculiar property of not to be subjected to screening by the conventional magnetic shields.

A focus of this paper will be on the special case of noninertial motion of practical interest: magnetic fields in electrostatic systems residing at rest in the curved space-time of the rotating body. Our principal statement is that the full system of Maxwell equations entails a non-vanishing magnetic field despite zero currents. In terms of the metric properties of the relevant space-time, the effect is due to the nonvanishing off-diagonal components of the metric tensor $g_{0i},\ i=1,2,3$. In cases of practical interest, this geometric magnetic field is driven by the Earth rotation.

Our primary motivation was the recent interest in the impact of the gravity induced spin rotations in the ultimate sensitivity searches for the electric dipole moments (EDM) of neutrons and of charged particles in the storage ring experiments Silenko and Teryaev (2016)-Abusaif et al. (2019). The effect of the Earth rotation in the neutron EDM experiments has already received much attention Baker et al. (2007), Lamoreaux and Golub (2007), Serebrov, Kolomenskiy, Pirozhkov, Krasnoschekova, Vassiljev, Polyushkin, Lasakov, Murashkina, Solovey, Fomin, Shoka and Zherebtsov (2015). The impact of the geometric magnetic field has not been discussed before. As we shall see, this field has a strong dependence on the configuration of the static electric field. In the geometry of the neutron EDM experiments, the geometric magnetic field is a nonuniform one with the nonvanishing gradient. Such a gradient of the magnetic fields is of particular concern in the comagnetometry which is crucial in the neutron EDM experiments (Abel, Ayres, Baker, Ban, Bison, Bodek, Bondar, Crawford, Chiu, Chanel, Chowdhuri, Daum, Dechenaux, Emmenegger, Ferraris-Bouches, Flaux, Geltenbort, Green, Griffith, van der Grinten, Harris, Henneck, Hild, Iaydjiev, Ivanov, Kasprzak, Kermaidic, Kirch, Koch, Komposch, Koss, Kozela, Krempel, Lauss, Lefort, Lemiere, Leredde, Mohanmurthy, Pais, Piegsa, Pignol, Quemener, Rawlik, Rebreyend, Ries, Roccia, Rozpedzic, Schmidt-Wellenburg, Schnabel, Severijns, Virot, Weis, Wursten, Wyszynski, Zejma and Zsigmond (2019) and references therein) . To this end, the salient feature of the geometric magnetic field is that it changes a sign when the external electric field is flipped. Consequently, interaction of the magnetic moment of the particle with the geometric magnetic field will generate a false signal of EDM. Numerically, the false EDM is still below the sensitivity of the current neutron EDM experiments, but will be in the ballpark of the ultimate sensitivity neutron EDM experiments aiming atvimproving the existing upper bound by one order in magnitude and beyond Chupp, Fierlinger, Ramsey-Musolf, Singh (2019).

The further presentation is organized as follows. A treatment of physics in the noninertial frames demands for a consistent use of Einsten’s general relativity (GR) formalism. In Section II we invoke GR to derive the results for the geometry of the curved spfce-time and dynamics of spinning charged particles in such spaces. In particular, since the space-time is curved, one needs to elucidate the meaning of the magnetic field. Specifically, we check that the magnetic field excerts on a charged particle a force which rotates its velocity ${\bf v}$ (defined in a certain orthonormal basis) preserving its magnitude, and induces the conventional precession of its magnetic moment. Trivial though may they sound, these constraints are important, when the nonvanishing geometrical magnetic field is present.

A subject of Section III is a derivation of the geometric magnetic field in the pure electrostatic systems at rest in the noninertial reference frames. We demonstrate that the geometric magnetic field is a salient feature of Maxwell equations in the stationary metrics with $g_{0i}\neq 0$. In the simplest case of the laboratory residing on the rotating Earth, a connection between the geometric magnetic field, which is an axial vector, and the polar external electric field, contains the angular velocity of Earth’s rotation, so that the parity conservation is warranted.

In Section IV we comment on three examples of the geometric magnetic field. The most interesting case at present is a search for the EDM of neutrons via precession of the neutron spin in the external electric field. Here the geometric magnetic field happens to have a finite gradient along the electric field and gives rise to the false EDM signal. Although tiny, this false EDM is in the ballpark of future ultimate precision neutron EDM experiments. The second example is the geometric magnetic field in all electric storage rings. Such a ring is discussed as a dedicated proton EDM machine Abusaif et al. (2019); Anastassopoulos et al. (2016). In this case, in contrast to the constant residual magnetic field of the Earth, the geometric magnetic moment has a quadrupole-like behaviour along the circumference of the storage ring. Finally, for the sake of academic completeness, we report an exact solution for the uniformly charged sphere at rest on Earth’s surface.

The major points of this study are summarized in the Conclusions. The Appendix provides a brief survey of the GR formalism our derivations are based upon.

II Maxwell equations and particle and spin dynamics

II.1 Maxwell equations

The local coordinates are denoted as $x^{\mu},\,\mu=0,1,2,3=(0,i)$, $i=1,2,3$. Let the electromagnetic field (2-form) be expressed through 4-potential $A_{\mu}$ (1-form) in local coordinates as:

[TABLE]

Then the homogeneous Maxwell equations

[TABLE]

are satisfied automatically. The field (1) is a holonomic one.

We have also the inhomogeneous Maxwell equations in the local coordinates:

[TABLE]

The electric and magnetic fields in ONB are defined by the usual rules (see Appendix D):

[TABLE]

[TABLE]

II.2 The particle and spin dynamics

Let $\mathop{\rm d}\nolimits s\equiv\sqrt{g_{\mu\nu}\mathop{\rm d}\nolimits x^{\mu}\mathop{\rm d}\nolimits x^{\nu}}$ be infinitely small interval when particle is moving along its world line. We use the standard notations $DX^{\mu}/\mathop{\rm d}\nolimits s\equiv u^{\nu}\nabla_{\nu}X^{\mu}$, where $u^{\nu}\equiv\mathop{\rm d}\nolimits x^{\nu}/\mathop{\rm d}\nolimits s$, so that $g_{\mu\nu}u^{\mu}u^{\nu}=1$. According to the definition and Eqs. (71), (72) we have

[TABLE]

Then the equation of motion of a charged particle in the local coordinates has the form

[TABLE]

Evidently, inside the ”freely falling elevator” in the Riemann normal coordinates (84), the dynamic equations are of the same form as in a flat space-time with Cartesian coordinates:

[TABLE]

Here $\mathop{\rm d}\nolimits t$ is the proper time interval in the laboratory frame of reference in which the particle moves with the velocity ${\bf v}$. In the general case

[TABLE]

A complete description of spin precession of relativistic particle with MDM and EDM in the gravity field is found in works Orlov, Flanagan and Semertzidis (2012)-Vergeles and Nikolaev (2019). Let ${\bf S}$ be the polarization vector in the particle rest frame. In the Riemann coordinates, the precession angular velocity ${\bf\Omega}$ equals Frenkel (1926)-Fukuyama and Silenko (2013)

[TABLE]

Here

[TABLE]

are MDM and EDM, correspondingly.

A formulation of the spin precession problem demands the orthonormal basis (ONB). It is important that the fields ${\bf E}$ and ${\bf H}$ in different ONB differ only by the Lorentz transformation. However, when describing dynamics in arbitrary ONB, a connection appears in the equations. As a result, Eqs. (8) take the following form:

[TABLE]

Since the polarisation vector ${\bf S}$ is defined in the particle rest frame, the explicit expression for a contribution of the connection to the precession frequency (10) for a relativistic particle is too lengthy to be reproduced here Vergeles and Nikolaev (2019), we only cite a result for the nonrelativistic case:

[TABLE]

The general conclusion from this brief discussion is that the fields ${\bf E}$ and ${\bf H}$ in ONB defined according to (4) and (5) possess all the dynamic properties of the electric and magnetic fields, correspondingly. We use hereafter the ONB (67) which is a physically sensible one.

III Magnetic field in the pure electrostatic system in the nonertial frame

Let’s consider a noninertial space-time with the stationary metric in the framework of Einstein’s general theory of relativity. The stationarity means that all components of the metric tensor $g_{\mu\nu}(x)$ are time-independent, i.e. $(\partial/\partial x^{0})g_{\mu\nu}=0$. A generic feature of the correspontent metric is nonvanishig off-diagonal elements,

[TABLE]

The outlined coordinate frame is denoted as $K$. In the cases of practical interest treated in Section IV the reference frame K is the laboratory frame used in the description of terrestrial experiments. We use also the inertial reference frame of distant stars $K^{\prime}$. The gravitating body and the corresponding frame $K$ rotate w.r.t. to the frame $K^{\prime}$ with the angular velocity ${\bm{\omega}}$. Our electrostatic system resides at rest on the rotating body and the corresponding electric 4-current, as defined in the frame $K$, is the stationary one:

[TABLE]

In this Section we derive the main result of the paper: For the stationary metric (13) in the case of stationary electric 4-current (14) the Maxwell equations can not be satisfied with the vanishing magnetic field (see below Eq. (20)).

A proof of this property proceeds as follows. Note that any antisymmetric field in the three-dimensional space $\sqrt{-g}F^{ij}=-\sqrt{-g}F^{ji}$ can be represented as

[TABLE]

The decomposition (15) would be unique in three-dimensional Euclidean space for fields decreasing at infinity. As well known, in the rotating frames one faces the formal issue of the horizon. However, in all the cases of practical interest, the charge and current distributions are localized way inside the horizon radius. Consequently, the decomposition (15) is well defined and unique, and $J_{i}=0$ entails ${\cal A}_{i}=0$.

According to the rules (70), the field definition (5) and the fact that $\tilde{e}^{i}_{0}=0$, we obtain:

[TABLE]

Here we used one of the relations (83). Making use of (15) (with ${\cal A}_{i}=0$) in (16) yields

[TABLE]

Now we turn to the homogeneous Maxwell equations (2). We express the holonomic field $F_{\mu\nu}$ in terms of the electric and magnetic fields:

[TABLE]

[TABLE]

Eq. (2) with $\mu=0$ imply the identity $\varepsilon_{ijk}\partial_{k}F_{ij}=0$. Therefore, applying the operator $\varepsilon_{ijk}\partial_{k}$ to the Eq. (19), using (17) and the identity $\varepsilon_{ijk}\partial_{k}F_{i0}=-\varepsilon_{ijk}\partial_{k}\left(e^{0}_{0}e^{\alpha}_{i}{\bf E}^{\alpha}\right)=0$ (see (18)), we obtain the equation

[TABLE]

Recall that $e^{0}_{j}=g_{0j}/\sqrt{g_{00}}\neq 0$ according to (81) and (13). Therefore if ${\bf E}\neq 0$, then (20) entails $\partial_{j}\psi\neq 0$ and, according to Eq. (17), ${\bf H}\neq 0$ as well. We shall refer to the field (17) as the geometric magnetic field ${\bf H}_{\bm{\omega}}$. As we shall see, the geometric magnetic field has its origin in the rotation of the frame $K$, hence the subscript ${\bm{\omega}}$.

Next we look into the equation for the electric field. Let’s express $F^{i0}$ in terms of physical fields:

[TABLE]

The substitution of the right-hand side of (21) into Eq. (3) with $\mu=0$ leads to

[TABLE]

Here we have used representation (17) and one of the relations (83).

The system of equations (17), (20) and (22) is complete and exact. The electric field ${\bf E}$ is expressed in terms of the potential $A_{0}$ with the help of relation (18).

To be specific, we describe the curved space-time of the rotating body by the Kerr metric described in Appendix B. We report here the iterative solutions for the electromagnetic potentials $A_{0}$ and $\psi$ and the corresponding fields, treating the angular velocity of the rotating body ${\bm{\omega}}$ and its gravitational radius $r_{g}$ as small parameters. Making use of (79)-(80), we rewrite Eqs. (17), (20) and (22) keeping the terms of the order of ${\bm{\omega}}$, $r_{g}$, $r_{g}{\bm{\omega}}$, ${\bm{\omega}}\otimes{\bm{\omega}}$.

There emerges a simple hierarchy of corrections in powers of ${\bm{\omega}}$. We are interested in the electrostatic system with the initial potential $A_{0}^{(0)}$ and the corresponding Minkowski space defined electric field ${\bm{E}}^{(0)}=-\nabla A_{0}^{(0)}$. Concerning the geometric magnetic field, the crucial points are Eq. (B3) for the off diagonal $g_{0i}$, Eq. (B5) for $e^{0}_{i}$ and Eq. (B6) for $\tilde{e}^{0}_{\alpha}$. The three related quantities are all proportional to the angular velocity of the rotating body. Consequently $e^{0}_{j}/e^{0}_{0}=\mathop{\rm O}\nolimits({\bm{\omega}})$ and, according to (20), the expansion for the magnetic potential $\psi$ starts with the linear term $\psi^{(1)}=\mathop{\rm O}\nolimits({\bm{\omega}})$.

Now we proceed to the electric field. According to the expansion of the Kerr metric in Appendix B, we have

[TABLE]

which does not contain the linear term. Then, equations (22) and (23) tell that the electric field acquires the first correction only to the second order in ${\bm{\omega}}$, i.e., ${\bf E}^{(1)}=0,\,{\bf E}^{(2)}\neq 0$. In the due turn, Eq. (20) guarantees that the quadratic correction to the magnetic potential vanishes: $\psi^{(2)}=0$.

Going back to the magnetic potential $\psi$, we make use of

[TABLE]

and invoke the relations

[TABLE]

Then Eq. (20) takes the final form

[TABLE]

A departure of the metric of the curved space-time from the Minkowski one changes the relationship between the potentail $A_{0}$ and the charge distribution and the relationship between the electric field ${\bf E}^{(2)}$ and the gradient of $A_{0}^{(2)}$. Although these second order corrections are unlikely to be of any practical significance in the terrestrial laboratories, we cite them for the sake of completeness. To the zeroth order in ${\bf\omega}$, but keeping the terms linear in $r_{g}$, we have

[TABLE]

After some algebra, Eq. (22) with the account for Eqs. (23), (25) yields the equation for the second order correction to the scalar potential,

[TABLE]

In this derivation, to the desired accuracy, we made use of $e^{0}_{j}/e^{0}_{0}=-[{\bm{\omega}},\,{\bf R}]/c$. Note how the second order correction to the electric potential couplles to the first order potential for the geometric magnetic field.

Finally, upon using Eq. (24) and gathering together all second order corrections, the Eq. (18) yields the second order correction to the electric field

[TABLE]

Apart from the gradient of the second order scalar potential, here emerge corrections quadratic in the angular velocity ${\bf\omega}$. The system of equations (26)-(29) is complete to the desired accuracy.

The rotating body of the practical interest is the Earth. It is the case of weak gravity. On the terrestrial surface $r_{g}/R_{\oplus}\sim 10^{-9}$ and $\omega R_{\oplus}\sim 1.5\cdot 10^{-6}$. Hence we keep the terms $\mathop{\rm O}\nolimits(|{\bm{\omega}}|)$ and neglect the former corrections. In this approximation equations (79)-(80) simplify to

[TABLE]

The second order corrections to the electric potential and electric field can be neglected and we have the familiar ${\bf E}=-\nabla A_{0}$ and the Poisson equation $\mathop{\rm div}\nolimits{\bf E}=-\Delta A_{0}=4\pi J^{0}$, while the Poisson equation for the potential of the geometric magnetic field takes a simple form,

[TABLE]

to be used in the subsequent analysis of terrestrial experiments.

IV Manifestations of the geometric magnetic field

IV.1 False EDM signal in the neutron EDM experiments

Here we comment on the possible implications of the geometrical magnetic field for the neutron EDM experiments. The fundamental observable is the change of the Larmor precession frequency

[TABLE]

from the parallel to antiparallel uniform fields. The EDM is extracted from the frequency shift

[TABLE]

The implicit assumption is that flipping the electric field does not change the magnetic one. Our point is that this is not the case with the geometric magnetic field.

In practice the electric field is generated in the plane capacitor with the gap much smaller than the plate size. The relevant solution of the one-dimensional problem for the geometric magnetic field proceeds as follows. In the gap in between the plates we have

[TABLE]

Now we solve the Poisson equation (31) for the magnetic potential, representing the electric field through $A_{0}(z)$ and integrating by parts:

[TABLE]

where $\Delta^{-1}(z)=|z|/2$ is the inverse to the Laplace operator. One more differentiation yields

[TABLE]

The geometric field ${\bf H}_{\bm{\omega}}$ is parallel to the external electric field ${\bf E}_{0}$. Its salient feature is the nonvanishing constant gradient

[TABLE]

At the mid plane the geometric field vanishes: ${\bf H}_{\bm{\omega}}(0)=0$.

The crucial feature of the neutron EDM experiments is the comagnetometry: one measures the neutron spin precession frequency with respect to that of the mercury comagnetometer. The mercury atoms are evenly distributed in the volume of the neutron storage cell and the average geometric magnetic field acting on the mercury comagnetometer vanishes: $\langle{\bf H}_{\bm{\omega}}^{(Hg)}\rangle={\bf H}_{\bm{\omega}}(0)=0$. The centre of mass of neutrons differs from that of the mercury by the offset $\langle z\rangle$, what entails the nonvanishing average geometric magnetic field acting on the magnetic moment of neutrons

[TABLE]

The most important point is that this geometric field changes the sign when the electric field is flipped. The net effect is that the apparent EDM of neutrons, $d_{n}^{obs}$, as given by the procedure (33), will acquire the false component, $d_{n}^{obs}=d_{n}+d_{false}$, where

[TABLE]

In the experiment Pendlebury, Afach, Ayres, Baker, Ban, Bison, Bodek, Burghoff, Geltenbort, Green, Griffith, van der Grinten, Grujic, Harris, Helaine, Iaydjiev, Ivanov, Kasprzak, Kermaidic, Kirch, Koch, Komposch, Kozela, Krempel, Lauss, Lefort, Lemiere, May, Musgrave, Naviliat-Cuncic, Piegsa, Pignol, Prashanth, Quemener, Rawlik, Rebreyend, Richardson, Ries, Roccia, Rozpedzic, Schnabel, Schmidt-Wellenburg, Severijns, Shiers, Thome, Weis, Winston, Wursten, Zejma and Zsigmond (2015) the neutron center of mass offset was $\langle z\rangle\simeq 2.8$mm, the more recent experiment reports $\langle z\rangle\simeq 3.9$mm. Taking the former, we find $d_{false}\approx 2.5\times 10^{-28}$ e$\cdot$cm. It is still way below the recently reported best upper bound on the neutron EDM, $d_{n}=(0.0\pm 1.1_{stat}\pm 0.2_{sys})\times 10^{-26}$e$\cdot$cm Abel, Afach, Ayres, Baker, Ban, Bison, Bodek, Bondar, Burghoff, Chanel, Chowdhuri, Chiu, Clement, Crawford, Daum, Emmenegger, Ferraris-Bouchez, Fertl, Flaux, Franke, Fratangelo, Geltenbort, Green, Griffith, van der Grinten, Grujić, Harris, Hayen, Heil, Henneck, Hélaine, Hild, Hodge, Horras, Iaydjiev, Ivanov, Kasprzak, Kermaidic, Kirch, Knecht, Knowles, Koch, Koss, Komposch, Kozela, Kraft, Krempel, Kuźniak, Lauss, Lefort, Lemière, Leredde, Mohanmurthy, Mtchedlishvili, Musgrave, Naviliat-Cuncic, Pais, Piegsa, Pierre, Pignol, Plonka-Spehr, Prashanth, Quéméner, Rawlik, Rebreyend, Rienäcker, Ries, Roccia, Rogel, Rozpedzik, Schnabel, Schmidt-Wellenburg, Severijns, Shiers, Tavakoli Dinani, Thorne, Virot, Voigt, Weis, Wursten, Wyszynski, Zejma, Zenner and Zsigmond (2020), but can become sizeable in the next generation of the neutron EDM experiments aiming at $d_{n}<10^{27}$ e$\cdot$cm Chupp, Fierlinger, Ramsey-Musolf, Singh (2019). With the neutron storage cell of height 12 cm, the geometric magnetic field induced spread of the false EDM within the ensemble of stored neutrons can be as large as

[TABLE]

IV.2 Geometric magnetic field as a background in all electric proton EDM storage rings

The principal idea of searches for the proton EDM in the all electric ring, run at the so-called magic energy, is to eliminate the magnetic field acting on the proton magnetic moment. Then the sole rotation of the proton spin will be due to interaction of its EDM with the radial electric field that confines protons in the storage ring, and very ambitious sensitivity to the proton EDM,

[TABLE]

is in sight Anastassopoulos et al. (2016); Abusaif et al. (2019). Here we comment on implications of the geometric magnetic field for such proton EDM experiments.

The storage ring is a cylinder capacitor with the gap $d$ which is much smaller compared to the height of cylinders $h$, which in its turn is much smaller than radii of cylinders $r_{1,2}=\rho\mp d/2$. In view of $d\ll h\ll\rho$ we neglect the dependence on the vertical coordinate and have the two-dimensional geometry. The beam trajectory is in the midplane of the storage ring at the orbit radius $|{\bf r}|=\rho$. The electric field in the gap is given by

[TABLE]

It is instructive to start with the storage ring located on the North pole. From the viewpoint of distant observer in the reference frame K’, the ring rotates with the Earth’s rotation angular velocity ${\bf\omega}$. The static charges on the two rotating cylinders produce the opposite currents and generate in the gap the magnetic fields of the same sign and magnitude. The net result is the magnetic field

[TABLE]

At first sight, it introduces an asymmetry between the clockwise and anticlockwise beams in the all electric storage ring, parasitic from the viewpoint of searches for the proton EDM. However, it is basically the motional magnetic field and, to the experimenter in the terrestrial laboratory K, it vanishes entirely. Such an exact cancellation only holds at the North and South poles, and at an arbitrary latitude it does not work. In the generic case, the result (43) for the magnetic field suggests the small parameter $\eta_{\bm{\omega}}={|{\bm{\omega}}|\rho}/{c}$, similar to that appearing in Eq. (39). For the storage ring of radius $\rho\sim 50$m we have

[TABLE]

which is some four orders in magnitude larger than the target value $\eta_{EDM}\sim 10^{-15}$.

Now we solve for the geometric magnetic field following the formalism of Section III. The electric field (42) suggests for the magnetic potential $\psi$ the Ansatz

[TABLE]

where ${\bm{\omega}}_{t}$ is a projection of the Earth’s angular velocity onto the ring plane. A generic solution to Eq. (31) is

[TABLE]

and

[TABLE]

where ${\bf n}={\bf r}/r$ .

The constant $\zeta$ is fixed by the boundary condition that the electric potential $A_{0}$ vanishes rapidly beyond the capacitor, so that in the integral representation for $\psi$ one can perform the integration by parts:

[TABLE]

where $\Delta^{-1}({\bf r})=(1/2\pi)\ln|{\bf r}|$ is the inverse to the Laplace operator. The resulting equation for the symmetric matrix $A_{ij}({\bf r})$

[TABLE]

entails

[TABLE]

Hence the expansion of $A_{ij}$ into irreducible tensor structures is of the form

[TABLE]

and a comparison to (47) gives immediately

[TABLE]

Our final result for the geometric magnetic field in the gap of the storage ring is

[TABLE]

where in the last step we neglected $|\log(r/\rho)|<d/(2\rho)\ll 1$.

The background magnetic fields are of prime concern to the planned searches for the proton EDM in the all electric magic storage rings, for a detailed discussion see the recent monographic document by the CPEDM (Charged Particles EDM) collaboration Abusaif et al. (2019). Important virtue of the all electric rings is a cancellation of many systematic effects when one compares spin rotations of simultaneously stored clockwise (CW) and anticlockwise (ACW) rotating protons. The magnetic Lorentz forces split the orbits of the CW and ACW beams. As discussed extensively in Abusaif et al. (2019), the modern techniques allow a very strong, but as yet incomplete, screening of the Earth’s magnetic field. Despite much work, reported in Abusaif et al. (2019), the analysis of the magnetic imperfection effects is still in the formative stage.

The Earth’s magnetic field ${\bf H}_{\oplus}$ and the geometric magnetic field ${\bf H}_{\bf\omega}$ do differ markedly. In contrast to the Earth’s magnetic field, the geometric one is not subject to screening by magnetic shields. On the scale of the storage ring, the Earth’s magnetic field can be regarded as a uniform one and has the constant projection onto the ring plane. It is pointing along the (magnetic) meridian, which we take for the y-axis : ${\bf H}_{\oplus}^{t}=(0,H_{\oplus}^{t})$. In contrast to that, the geometric magnetic field has the quadrupole-like behaviour along the particle orbit, ${\bf H}_{\bf\omega}=H_{\bf\omega}(\sin 2\theta,\cos 2\theta)$. The angular position of the particle in a ring, $\theta$, is defined by ${\bf n}=(\cos\theta,-\sin\theta)$.

In the all electric magnetic rings the most dangerous ones are the radial magnetic fields. In the above two cases they are equal to $H_{\oplus}^{(r)}=({\bf n}\cdot{\bf H}_{\oplus})=-H_{\oplus}\sin\theta$ and $H_{\bf\omega}^{(r)}=({\bf n}\cdot{\bf H}_{\bf\omega})=H_{\bf\omega}\sin\theta$  According to Abusaif et al. (2019), to the first approximation the rotation of the proton spin is proportional to the one-turn integral $\oint d\theta H_{\oplus}^{(r)}$. .Obviously, both the Earth’s and geometric magnetic fields share the property

[TABLE]

The argument of Ref. Abusaif et al. (2019) about vanishing false EDM signal from is then applicable to the geometric magnetic field as well.

The above consideration is somewhat naive and must be complemented by a consistent treatment of the spin-orbit coupling, though. Namely, this cancellation of the false EDM effect might become incomplete because of the orbit distortions which are very much distinct in the two cases. To be on the safe side, one needs a dedicated analysis of the false spin rotations with simultaneous allowance for the orbit distortions. Furthremore, one needs to pay an attention to a possible cross talk between the impact of the geometric magnetic field and the residual Earth’s magnetic field. It is an important complex issue on its own to be addressed to in the future, it goes beyond the scope of the present communication.

IV.3 Geometric magnetic field of the conducting charged sphere

For the sake of completeness, we comment on the charged sphere at rest in the rotating system $K$. Inside the sphere we have ${\bf E}=0$ and $A_{0}=\mathop{\rm const}\nolimits$ thereof:

[TABLE]

Here $Q$ is the charge of the sphere, and the radius-vector ${\bf r}=0$ at the centre of the sphere.

According to (31)

[TABLE]

and

[TABLE]

Here ${\bf n}={\bf r}/{r}$ and

[TABLE]

is the geometric magnetic moment of the conducting charged sphere, induced by the Earth’s rotation. The second term in (63) is the familiar motional magnetic field in the rotating frame $K$ which is entailed by the electric field in inertial frame $K^{\prime}$.

V Conclusions

We have shown that in the pure electrostatic systems at rest on the rotating bodies there can exist a geometric magnetic field. From the general relativity point of view, it originates from the nonvanishing off-diagonal elements $g_{0i}$ of the metric tensor which are proportional to the angular velocity of rotation of the gravitating body. We presented several specific examples of the geometric field for different configurations of the static electric field on the rotating body described by the Kerr metric. In the configuration of experimental setups used in the terrestrial searches for the EDM of neutrons, the geometric magnetic field changes the sign when the electric field is flipped. Consequently, its interaction with the magnetic moment of the neutron can imitate the neutron EDM and that can become a sizable background in the next generation of the neutron EDM experiments. We found a fairly large background geometric magnetic field in all electric magic storage rings considered a dedicated machine for searches of the proton EDM. The symmetry properties of the geometric magnetic field suggest strong cancellations of its contribution to the proton spin rotations. Stil, its impact on the signal of EDM remains an open issue - here one needs a dedicated analysis with full allowance for the spin-orbit dynamics in the storage ring.

Acknowledgements.

We are grateful to A.Ya. Maltsev, A.A. Starobinsky and S.S. Vergeles for valuable comments and discussions. This work was carried out as a part of the State Program 0033-2019-0005.

Appendix A Geometry

We consider stationary metric in the reference frame K. Following the textbook Landau and Lifshitz (1971), let’s diagonalize this quadratic form:

[TABLE]

where the field $e^{a}_{\mu}(x)$ is called tetrad, $a,b,\ldots=0,1,2,3,\,a=(0,\alpha),\,\alpha=1,2,3,\,\eta_{ab}=\mathop{\rm diag}\nolimits(1,-1,-1,-1))$. Two infinitesimally close events are simultaneous if 1-form

[TABLE]

According to (64) and (65) the squared interval between simultaneous events is

[TABLE]

The local orthonormal basis (ONB) $\tilde{e}^{\mu}_{a}(x)$ is defined by the equations

[TABLE]

Since according to (66)

[TABLE]

one readily finds:

[TABLE]

The rules of the tensor component transitions from coordinate basis to ONB and vice versa are standard. For example

[TABLE]

In ONB the tensor indices are lowered and raised with the help of metric tensors $\eta_{ab}$ and $\eta^{ab}$. With the above chosen tetrad there is a complete equivaence between $J^{i}=0$ and $J^{\alpha}=0$:

[TABLE]

The covariant derivatives $\nabla_{\mu}$ in the coordinate basis and in ONB are related as

[TABLE]

where $\gamma^{a}_{bc}\equiv\tilde{e}^{\mu}_{c}\gamma^{a}_{b\mu}$ is the totality of the connection coefficients, and

[TABLE]

The fact that the connection is free of torsion is fixed by the equations

[TABLE]

The connection coefficients are determined uniquely by Eqs. (72) and (73):

[TABLE]

Under the local Lorentz transformation

[TABLE]

the connection 1-form $\gamma^{a}_{b\mu}\mathop{\rm d}\nolimits x^{\mu}$ transforms as follows:

[TABLE]

Here $\Lambda^{a}_{b}(x)$ is a local Lorentz transformation matrix.

Appendix B Metric and tetrad

To proceed further, we need to choose the appropriate metric. Let the rotating reference frame K be defined for the Earth rotating with constant angular velocity ${\bm{\omega}}$ in the inertial frame of distant stars $K^{\prime}$ . The local coordinates, vectors etc. in $K^{\prime}$ are denoted as $x^{\prime\mu}$, ${\bf R}^{\prime}$ and so on. In the frame $K^{\prime}$ we have the Kerr metric of the rotating Earth. For the purposes of our analysis it is sufficient yo use a limit of weak gravity and nonrelativistic rotation velocity and we expand the Kerr metric retaining the terms linear in ${\bm{\omega}}$ and $r_{g}$, bilinear in $r_{g}$ and ${\bm{\omega}}$ and quadratic in ${\bm{\omega}}$:

[TABLE]

where $R^{\prime}=|{\bf R}^{\prime}|$, $k=6,674\cdot 10^{-8}\mbox{cm}^{3}\cdot\mbox{g}^{-1}\cdot\mbox{sec}^{-2}$ is the gravitation constant, $M_{\oplus}$ and $C=IM_{\oplus}R_{\oplus}^{2}$ are the Earth mass and moment of inertia relative to polar axis, $I=0.3307$, $M_{\oplus}=5,972\cdot 10^{27}\mbox{g}$, $R_{\oplus}=6,378\cdot 10^{8}\mbox{cm}$, $r_{g}=2kM_{\oplus}/c^{2}=0,887\mbox{cm}$. Next we transform the metric (77) into the metric in the frame $K$ rotating with angular velocity ${\bm{\omega}}$ relative to the frame $K^{\prime}$. We take the coordinates in $K$ and $K^{\prime}$ having the same origin at the centre of Earth. The local coordinates in the frame $K$ are denoted as $x^{\mu}$ and, by definition, they are connected with coordinates $x^{\prime\mu}$ as follows,

[TABLE]

and the metric $g_{\mu\nu}$ in the frame $K$ equals

[TABLE]

The inverse metric tensor $g^{\mu\nu}=\eta^{ab}\tilde{e}^{\mu}_{a}\tilde{e}^{\nu}_{b}$ is

[TABLE]

To the same approximation the tetrad equals

[TABLE]

and the ONB vector fields are

[TABLE]

Appendix C Useful relations

The following relations are used in the main body of the text:

[TABLE]

which imply that

[TABLE]

and so forth.

Appendix D Riemann normal coordinates

For the correct interpretation of the electromagnetic fields in a curved space-time it is useful to keep in mind the form of dynamic equations in the Riemann normal coordinates. The Riemann normal coordinates $y^{\mu}$ can be introduced in the vicinity of any point $p$, so that $y^{\mu}(p)=0$ and

[TABLE]

where $\mathop{{\mathfrak{R}}}\nolimits_{ab\,\nu\mu}$ is the Riemann curvature tensor. Since near the Earth surface $|\mathop{{\mathfrak{R}}}\nolimits_{ab\,\nu\mu}|\sim r_{g}/R_{\oplus}^{3}\sim 0,5\cdot 10^{-26}\mbox{cm}^{-2}$, one can ignore the space-time curvature in the small vicinity of point $p$. This vicinity with Riemann normal coordinates $\{y^{\mu}\}$ is the mathematical model of a ”freely falling elevator”. Therefore, when writing differential equations in the center of normal coordinates, the curvature can be neglected. It follows from here that differential dynamic equations in normal Riemann coordinates $\{y^{\mu}\}$ inside the ”elevator” have the same form as in the Cartesian coordinates in the Minkowski space. Obviously, in the Riemann normal coordinates inside the ”elevator”, all field components in the ONB $F^{ab}=e^{a}_{\mu}e^{b}_{\nu}F^{\mu\nu}$ coincide with of the same named field components in the Riemann coordinates. Upon transition to arbitrary ONB, all tensors are transformed in accordance with the usual Lorentz rules, and a connection determined by ONB appears in the equations of motion.