JHEPGC  Vol.7 No.2 , April 2021
High Index SMES Device for Gravitomagnetic Field Generation
Abstract: A method is described for creating a measurable unbalanced gravitational acceleration using a gravitomagnetic field surrounding a superconducting toroid as described by Forward (1962). An experimental Superconducting Magnetic Energy Storage (SMES) toroid configuration of wound superconducting nanowire is proposed to create a measurable acceleration field along the axis of symmetry, providing experimental confirmation of the additive nature of a Lense-Thirring derived gravitomagnetic field. In the present paper, gravitational coupling enhancement of this effect is explored using a high index or high permittivity material, as predicted by Sarfatti (2020) using his modification to Einstein’s General Relativity Field Equations for gravitational coupling in matter.

1. Introduction

When Forward [1] first proposed a gravitomagnetic toroid for unbalanced gravitational force production in 1962 that any experimental realization was quite impractical. However, recent advances in high temperature superconducting (HTSC) nanorod wire (nanowire) technology, described recently by Stephenson et al. [2] and Rieken et al. [3], have enabled a new class of superconducting magnetic energy storage (SMES) devices operating at current densities sufficient to develop measurable gravitomagnetic fields.

In the present paper, an experimental SMES toroid configuration is proposed that uses a super dielectric material (e.g. high permittivity super capacitor) to substantially improve gravitational coupling of mass flow to the curvature of spacetime as predicted by Sarfatti [4]. Depending on the nature of gravitational coupling, it is predicted that a set of standard accelerometers could measure acceleration fields along the axis of symmetry of the toroidal coil, thus providing experimental confirmation of the additive nature of the gravitomagnetic fields, as well as the production of a linear component of the overall acceleration field, see Figure 1 for details.

In the instantiation of Forward’s gravitational generation coil described in this paper, superconducting electron flow provides the change in mass current in the toroid. A high permittivity material such as a super dielectric material (SDM) [5] is added to the space in the center of the toroid to improve gravitational coupling. Alternatively or additionally high permittivity insulators may be used as an exterior jacket to the conducting portion of the nanowire so as to further boost gravitational coupling, as gravitomagnetic effects are largest adjacent to the mass flow.

2. Gravitomagnetic Force Equation for Toroid Mass Flow

As first developed in Forward 1962 Ref. [1] the linear force Gf developed by gravitomagnetic force in the mass flow toroid of Figure 1 is given by Equation (1):

G f = ( η 4 π ) ( N T ˙ r 2 R 2 ) (1)

Figure 1. Gravitational force generation coil from Forward 1962 [1] with an inspiralling mass current, with a vector potential P, creating gravitomagnetic field G, with high index material in the center to improve coupling.

where η is gravitomagnetic permeability, η = η o η r as described in reference [2], N is number of turns in the toroid coil winding, T ˙ represents the mass flow, r is the cross sectional radius of the toroid, and R is the overall radius of the toroid.

Mass flow can be generated by means of generating an electrical current by virtue of putting the mass of electrons in motion.

Electron Mass Flow

Single electron mass flow shown in Figure 2 is given by mass momentum in Equation (2):

T e = p e = ( Ω × r ) m e (2)

where Ω = angular rate, angular velocity is v = Ω × r in the classical case [6].

Change in mass flow for the single electron flow shown in Figure 2 is given by Equation (3):

T ˙ e = p ˙ e = a m e = ( Ω × v ) m e (3)

This is equivalent to centripetal force shown in Equation (4):

T ˙ e = F e = m e v 2 r = m e a e = m e ( ω 2 r ) (4)

where ω is the angular rate, given in Equation (5):

ω = 2 π t p = d θ d t (5)

3. Current in Idealized SMES

We now attempt to estimate the possible inspiralling current flux enabled by the emerging SMES technology as it relates to the core geometry constraints described in Figure 3, a toroid with torus geometry. We start with the assumptions needed to calculate the number of turns N.

From Equation (1) the torus assumptions made in Figure 3 can be factored into Equation (6) as follows:

G f = ( η o η r 4 π ) S ( N T ˙ r 2 R 2 ) (6)

Figure 2. Electron orbit around one loop of SMES toroid.

Figure 3. Toroid with a torus shaped core geometry, including high index scaling matter in center.


Gf = gravitomagnetic force;

ηo = absolute gravitomagnetic permeability;

ηr = relative gravitomagnetic permeability;

N = number of turns in the coil of the torus;

S = Sarfatti scaling factor;

T ˙ = change in mass, or mass flow;

r = cross section radius of torus;

R = centerline radius of torus.

Note that Equation (6) differs from Forward’s original Equation (1) formulation in that it includes the scaling factor S which accounts for the scaling of gravitational coupling due to high index matter or metamaterials, which was first expressed by Sarfatti as an additional zero rank tensor in Einstein’s Field Equations [4] [7]:

R u v 1 2 R g u v = ( 8 π G c 4 ) S T u v (7)


Ruv = Ricci curvature tensor;

Rguv = Ricci scalar curvature;

G = gravitational constant;

c = speed of light;

S = Sarfatti scaling tensor;

Tuv = Stress Energy Tensor.

For the purposes of describing an idealized case with a realistic geometry we develop a description of a device bounded by a 10 cm toroid centerline diameter, shown in Figure 3 and Figure 4, with a cross-sectional diameter of 1 cm. Furthermore 16 sectors are defined as shown in Figure 4.

We further add additional assumptions regarding to what extent conductors are wrapped around the toroid shaped device to determine constraints on the number of conductive loops that can be accommodated using the technology described in [2]. As shown in Figure 5 via cross section we assume here a conductor winding depth of 0.1 cm.

Figure 4. Torus section definition.

Figure 5. Conductor cross section in each sector.

The segmented share of the inside of cross section Cs is 1/16th of the overall inner circumference as given in Equation (8) as follows:

C s = 2 π r i 16 = 1.77 cm (8)

with a depth D the minimum inner loop cross sectional area can be described in Equation (9) as follows:

A s e c = D C s = 0.1 cm × 1.77 cm = 0.177 cm 2 (9)

This area is shown packed with conductors in Figure 6. Windings are depicted as packed in depth and along sector circumference.

Assuming each nanowire conductor has a diameter dc of 100 μm, then the cross-sectional area of each conductor will be given by Equation (10) as follows:

Figure 6. Conductor packing in each sector.

A c = π r c 2 = 7 . 854 × 10 9 m 2 (10)

For packing nanowire conductors in a cross-sectional area described in Figure 6, assume as a worst case rectangular area described by the shortest edges such that a number of conductors in depth, Nd, may be packed in one dimension, with the number of conductors, Ncs, packed in the other dimension. These packing counts may be calculated in Equation (11a) and Equation (11b) as follows:

N d = D c / d c = 0.1 cm 100 μ m = 10 (11a)

N c s = C s / d c = 1.77 cm 100 μ m = 177 (11b)

The total number of windings by sector will therefore be the product of Nd and Ncs:

N s e c = N d N c s = 10 × 177 = 1770 (12)

And with 16 sectors the total number of windings for the entire toroid will be:

N = 16 N c s = 28320 (13)

what is T ˙ with the forgoing assumptions?

In this idealized case electrons circulate about a coil of circumference cr, or slightly larger, as described by Equation (14):

c r = 2 π r = 3.14 cm (14)

Electron motion

Assume further a supply voltage of 16 KV, resulting in 16 KeV of kinetic energy for each electron, which corresponds with the upper limit of a non-relativistic case, where v = 0.25c, so thatγ = 1.06 - 1.0.

Then from Equation (4) for non-relativistic circular motion, the vector change in DC current flow is:

T ˙ e = m e v 2 r (15)

which for a single electron has the following values:

me = mass of the electron = 9.11 × 10−31 Kg;

v = velocity of the electron = 0.25c = 0.75 × 108 m/s;

r = 0.5 cm for the assumed geometry.

And where the angular acceleration of the electron is:

a e = v 2 r = 1.125 × 10 18 m / s 2 (16)

Thus, change in mass flow represents centripetal acceleration in the case of circular motion:

T ˙ e = m e a e = 10.25 × 10 13 N (17)

Equation (17) corresponds to the change in mass flow for one electron in one loop of coil. Total mass flow change is therefore the mass flow change per electron times the number of electrons:

T ˙ = T ˙ e N e (18)

Electron current

What is the number of electrons Ne in one loop in motion (part of the mass flow) at a given time for an assumed velocity of v = 0.25c? Ne in one loop can be described by the current I times the period of a single loop circulation ∆t:

N e = I Δ t (19)

where the period of an orbit can be described by:

Δ t = c r v = 2 π r v = 0.4189 ns (20)

what is the possible current inside the idealized device for the case where the entire winding is in series? We make the assumption about max current to stay below critical current density of 250 MA/m2as described in [2] [3].

Current is limited by the maximum permissible current density and the cross section of the conductort:

I = J A c (21)

where J is material dependent. For the nanowire assumed in Ref. [2], J = 250 MA/m2. Cross sectional area Ac = 7.854 × 10−9 m2 as given in Equation (10). Therefore, maximum current for this conductor diameter isI = 1.96 Amps.

Expanding on Equation (21) the number of electrons Ne in circulation in one loop may be calculated by noting that there are 6.2415 × 1018 electrons per Coulomb:

N e = ( electrons Coulomb ) I ( Coulombs sec ) Δ t = 5.12 × 10 9 electrons (22)

4. Forces in Idealized SMES

Expressing Equation (18) as force per electron times the number of electrons in motion in one loop:

T ˙ = T ˙ e ( Newtons electron ) N e ( #electrons ) = 5.25 mN (23)

Thus, each loop experiences about 5.25 mN of integrated centripetal force ( T ˙ ) due to the electrons in circulation within.

We now describe the scale factor to couple this force to the gravitomagnetic effect.

Revisiting Equation (6) which describes the overall linear force developed at the center of the toroidal coil, total gravitomagnetically developed force will be:

G f = ( η o η r ) ( S ) ( N T ˙ r 2 4 π R 2 ) = ( η o η r ) ( S ) ( 0.118 N ) (24)

where known variables have been grouped on the right and unknown variables have beencollected on the left. This raises the question what are the correct values for ηo and ηr?

Predicted force using gravitational potential scaling

If ηo goes as G/2c as does gravitomagnetic potential (Ref. [8], Equation (1.5)), then:

η o = G 2 c = 1.11 × 10 19 (25)

In this case G f = ( η r ) 1.31 × 10 20 . Values of ηr are experimentally unknown at this time. However if values of ηr track values of ε r 2 then values as high as ηr = 1016 may be possible, yielding Gf = 0.131 mN. For a test mass of 1 Kg this is equivalent to an easily measurable 13.4 micro-g’s. This is for a scaling factor of G/2c. If the scaling factor is G/c2 then the effect drops to an unmeasurably low level of 4.47 × 10−14 g’s, even for ηr = 1016.

Predicted force using gravitational field scaling

If the scaling factor is G/c3 or G/c4 then the effect is negligible for ηr = 1016, and a much higher index material will need to be found to improve coupling, or other means must be found for improving either the S tensor and/or gravitomagnetic permeability ηr.

Consider the case where ηo goes as 8πG/c4 as does the coupling constant in Einstein’s Field Equations, then:

η o = 8 π G c 4 = 2.07 × 10 43 (26)

which gives the following expression for total gravitational force developed:

G f = ( η o η r ) ( S ) ( 0.118 N ) = ( η r ) ( S ) 2.44 × 10 44 N (27)

A value of ηr=1016 gives G f = ( S ) 4.43 × 10 23 Newtons, which requires a scalar value of S of at least an additional 1016 to obtain a measurable value for Gf. Per [4] for an isotropic material S can be written as:

S = 1 2 ( ϵ 2 + 1 μ 2 ) = 1 2 ( n 4 + 1 ) / μ 2 (28)

Therefore either a permittivity of ε = 108 or an index of n = 104 with μ = 1 would be sufficient to provide the necessary 16 orders of magnitude in coupling improvement. However care should be taken that effects captured by the S factor are not duplicative of those accounted for in the value of ηr.

Per [5] super dielectrics can provide up to ε = 109 which, if we simultaneously had ηr = 1016 via some unrelated phenomena, would improve the overall generated gravitational force to a marginally measurable level:

G f = ( η r ) ( S ) 2.44 × 10 44 = 0.244 nN (29)

Alternatively metamaterials could also be investigated to establish whether non-isotropic materials may create a more advantageous gravitational coupling improvement.

Lenz’s Law Implications

Is there an equivalent of Lenz’s Law for gravitomagnetics? Consider a hypothetical case where a Forward toroid that is capable of supporting sufficient current and coupling to develop 1 G of acceleration. As a thought experiment imagine disconnecting the power supply and shorting the toroidal coil to itself, such that it is one continuous closed loop of nanowire, and equip a flying craft with such a coil.

At the surface of the Earth such a properly oriented toroidal coil would be immersed in a 1g gravitational field. Should not such a field generate a counter-current in a toroidal coil capable of supporting a current of that level, essentially providing a self-powered counter field? If so, a craft so equipped would be able to hover using a counter-field equal but opposite to the Earth’s gravitational field without need of additional power sources.

Similarly, additional closed loop toroidal coils, properly oriented, could also be used to “current charge,” collecting energy from a static gravitational field, potentially providing power that could be directed to other propulsive or non-propulsive functions of such a hypothetical craft.

5. Conclusion

An argument is made for using high permittivity materials in the core or donut hole and/or nanowire insulators of an SMES toroid to improve gravitomagnetic coupling for the creation of an acceleration field, possibly of measurable amplitude. Improved coupling will be beneficial for shrinking device scale and complexity and overcoming the very weak coupling between mass flow and gravitomagnetic spacetime curvature.


The authors wish to acknowledge Jack Sarfatti for helpful correspondence on his insights into altering Einstein’s Field Equations in the presence of matter and/or metamaterials. The financial support of Seculine Consulting is gratefully acknowledged.


HTSC: High temperature superconductor

SC: Superconductor

SDM: Super dielectric material

SMES: Superconducting magnetic energy storage

Cite this paper: Stephenson, G. (2021) High Index SMES Device for Gravitomagnetic Field Generation. Journal of High Energy Physics, Gravitation and Cosmology, 7, 367-376. doi: 10.4236/jhepgc.2021.72020.

[1]   Forward, R. (1962) Guidelines to Antigravity. Proceedings of the Gravity Research Foundation, New Boston.

[2]   Stephenson, G., Rieken, W. and Bhargava, A. (2019) Extended Cases of Laboratory Generated Gravitomagnetic Field Measurement Devices. Journal of High Energy Physics, Gravitation and Cosmology, 5, 375-394.

[3]   Rieken, W., Bhargava, A., et al. (2018) YBa2Cu3Ox Superconducting NanoRods. Japanese Journal of Applied Physics, 57, 023101.

[4]   Stephenson, G. (2020) 2000-2020 Summary of Gravitational Work. APEC Proceedings, Blaine, 12 December 2020.

[5]   Fromille, S. and Phillips, J. (2014) Super Dielectric Materials. Materials, 7, 8197-8212.

[6]   Knudsen, J.M. and Hjorth, P.G. (2000) Elements of Newtonian Mechanics: Including Nonlinear Dynamics. 3rd Edition, Springer, Berlin, 96.

[7]   Sarfatti, J. (2020) Explaining US Navy Close Encounters with Tic Tac UAV Metric Engineering. Proceedings of the Estes Park Advanced Propulsion Workshop, Space Studies Institute, North Hollywood, October 2020.

[8]   Mashhoon, B. (2008) Gravitoelectromagnetism: A Brief Review.