Comments on gravitoelectromagnetism of Ummarino and Gallerati in “Superconductor in a weak static gravitational field” vs other versions

Recently reported [Eur. Phys. J. C., 77, 549 (2017). https://doi.org/10.1140/epjc/s10052-017-5116-y] gravitoelectromagnetic equations of Ummarino and Gallerati (UG) in their linearized version of general relativity (GR) are shown to match with (a) our previously reported special relativistic Maxwellian Gravity equations in the non-relativistic limit and with (b) the non-relativistic equations derived here, when the speed of gravity cg\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c_g$$\end{document} (an undetermined parameter of the theory here) is set equal to c (the speed of light in vacuum). Seen in the light of our new results, the UG equations satisfy the Correspondence Principle (cp), while many other versions of linearized GR equations that are being (or may be) used to interpret the experimental data defy the cp. Such new findings assume significance and relevance in the contexts of recent detection of gravitational waves and the gravitomagnetic field of the spinning earth and their interpretations. Being well-founded and self-consistent, the equations may be of interest and useful to researchers exploring the phenomenology of gravitomagnetism, gravitational waves and the novel interplay of gravity with different states of matter in flat space-time like UG’s interesting work on superconductors in weak gravitational fields.


Introduction
In a recent interesting theoretical study on the interplay of superconductivity and weak static gravitational field, Ummarino and Gallerati [1] concluded that the reduction of the gravitational field in a superconductor, if it exists, is a transient phenomenon and depends strongly on the parameters that characterize the superconductor. The gravitational equations used by the authors in their study are represented by the following Gravito-Maxwell Equations: a e-mail: behera.hh@gmail.com and the equation of motion of a particle moving with nonrelativistic velocity, v, is given by Gravito-Lorentz force law: In the above equations ρ g = ρ 0 is the (rest) mass density, j g = ρ g v mass current density, the speed of gravitational waves in vacuum c g = c (the speed of light in vacuum) in [1], E g is the usual Newtonian gravitational field (called gravitoelectric field) and B g is the gravitational analogue of magnetic induction field (called gravito-magnetic field). The E g and B g fields are related to gravitoelectric scalar potential φ g and gravitomagnetic vector potential A g as Ummarino and Gallerati [1] derived these equations from Einstein's GR by linearization procedure in the weak field and slow motion approximations. We name the Eqs. (1)(2)(3)(4)(5)(6)(7) to represent the General Relativistic Maxwellian Gravity of Ummarino-Gallerati as (GRMG-UG) to differentiate it from other existing (or future) formulations that (may) result from different methods of study. For instance, in this communication we show how one can derive these equations in a nonrelativistic approach without invoking the space-time curvature and the linearization schemes of GR. By adopting, in essence, the non-relativistic approach of Schwinger et al. [2] for their derivation of Maxwell's equations and the Lorentz force law in the electromagnetic case, we derive the fundamental equations of Time-Dependent Galileo-Newtonian Gravitodynamics of moving bodies (called here as Non-Relativistic Maxwellian Gravity or NRMG in short) within the Galileo-Newtonian relativity physics by combining the following ingredients: (a) The validity of Newton's laws of gravitostatics, (b) The validity of the equation of continuity that expresses the law of conservation of mass, 1 and (c) Postulating the existence of gravitational waves traveling in free space with a finite velocity c g (which is an undetermined parameter of the theory but whose value may be fixed in measurement of physical quantities involving c g or by comparing the field equations with those obtainable from more advanced theories).
As will be shown here, our derived equations of NRMG match with the Eqs. (1-7) of GRMG-UG, if c g = c. Further, if c g = c, the equations of NRMG also match with those of Special Relativistic Maxwellian Gravity (SRMG) [3] 2 in flat space-time, where c g = c is a natural outcome of the theory there (here discussed in Sect. 3). From McDonald's [4] report of little known Heaviside's Gravity (HG) [4,5] of 1893, 3 we find that the equations of NRMG also match with those of HG in which Heaviside thought c g might be equal to c. These findings assume significance and relevance in the contexts of recent experimental detection of gravitational waves [6][7][8][9], gravito-magnetic field of the spinning Earth using the orbital data of two laser-ranged satellites (LAGEOS and LAGEOS II) [10][11][12][13][14] and the Gavity Probe B (GP-B) experimental results [15][16][17] vindicating Einstein's GR, because all of these results are being interpreted in the 1 Schwinger et al. [2] derived (b) using the Galileo-Newton principle of relativity (masses at rest and masses with a common velocity viewed by a co-moving observer are physically indistinguishable). Here, we will use this relativity principle for deriving the gravitational analogue of the Lorentz force law. 2 We named our relativistic gravity as Maxwellian Gravity since J. C. Maxwell [J.C Maxwell, Phil. Trans. Roy. Soc. London 155, 459-512 (1865), sec. 82: Note on the Attraction of Gravitation] first tried to construct a field theory of gravity analogous to Classical Electromagnetic theory and did not work on the matter further as "He was dissatisfied with his results because the potential energy of a static configuration is always negative but he felt this should be re-expressible as an integral over field energy density which, being the square of the gravitational field, is positive [3,4]. Maxwell concluded his note on the attraction of gravitation with the statement,"As I am unable to understand in what way a medium can possess such properties, I cannot go any further in this direction in searching for the cause of gravitation". 3  literature as new crucial tests of general relativity having no Galileo-Newtonian counterpart. For instance, Ciufolini et al. [10], in their report of measurement of the Lense-Thirring effect 4 stated: Newton's law of gravitation has a formal counterpart in Coulomb's law of electrostatics; however, Newton's theory has no phenomenon formally analogous to magnetism. On the other hand, Einstein's theory of gravitation predicts that the force generated by a current of electrical charge, described by Ampère's law, should also have a formal counterpart "force" generated by a current of mass.

Non-relativistic Maxwellian gravity (NRMG)
In Galileo-Newtonian physics, gravitational mass, 5 m g , is the source of Newtonian gravitational field, E g , which obeys the following two equations, viz., where ρ g = ρ 0 is the (positive) rest mass density (to be shown shortly in this section and also in the next section in the relativistic case) and G is Newton's universal gravitational constant having the value 6.674 × 10 −11 N · m 2 /kg 2 .
In an inertial frame, the gravitational force on a point particle having gravitational mass m g in a gravitational field E g is expressed as where φ g is gravitational potential at the position of m g . Further, the mass that appears in Newton's second law of motion, which Newton actually wrote in The Principia, is of the following form (valid only in inertial frames): is called the inertial mass, m, which is a measure of inertia of a body's resistance to any change in its state of motion represented by its velocity v or linear momentum p = mv. If m does not change with time, then the most powerful and profound Eq. (11) of Newton in classical physics, whose form survived special relativistic revolution, takes the familiar non-relativistic (v << c) form, where a is the acceleration of m as measured in an inertial frame. In Galileo-Newtonian physics, inertia is of two types, viz., inertia of rest and inertia of motion. Accordingly, inertial mass is of two types: rest mass (m 0 = a measure of inertia of a body at rest) and inertial 6 mass (m = a measure of inertia of a body in motion). Thus a qualitative distinction between two types of inertial masses exists in Newtonian physics, but there is no quantitative distinction, i.e., m = m 0 in Newtonian world of physics. 7 With this understanding, if the only force acting on the particle is the gravitational force, then the equation of motion of a particle of rest mass m 0 and gravitational mass m g can be obtained from Eqs. (10) and (12) as Equation (13) shows us that if m g = m 0 (or ρ g = ρ 0 for continuous mass distribution), then any particle of whatever rest mass will fall with the same acceleration in a given gravitational field. This is the law of universality of free fallknown to be true since the time of Galileo. In view of these observations, in this section term 'mass' is to be understood as the rest mass which represents the gravitational charge (or mass) of a particle. In Newtonian physics, mass obeys a local conservation law expressed by the equation of continuity: where j 0 = ρ 0 (r, t)v is the rest mass current density. With ρ g = ρ 0 , the three laws expressed in Eqs. (8), (9) and (14) have their individual indisputable validity in Newtonian Physics. But they are simultaneously valid only for the systems or situations where remain valid (here in (15) the time dependence of the physical quantities is made explicit). The Eqs. (8)(9)(10)(11)(12)(13)(14)(15) describe the familiar Galileo-Newtonian gravitodynamics. Newton's gravitational force law (10) implies "action-at-a-distance"the gravitational force seems to act instantaneously at a distance -which Newton considered as an absurd element of his theory. Despite his great efforts, Newton could not offer a plausible mechanism for resolving the problem of action-ata-distance and left the problem for others to resolve as evident from his 3rd letter to Bentley, dated February 5, 1692/3, where he wrote [18]: It is inconceivable, that inanimate brute Matter should, without the Mediation of something else, which is not material, operate upon, and affect other Matter without mutual Contact, as it must be, if Gravitation in the Sense of Epicurus, be essential and inherent in it. And this is one Reason why I desired you would not ascribe innate Gravity to me. That Gravity should be innate, inherent and essential to Matter, so that one Body may act upon another at a Distance thro' a Vacuum, without the Mediation of any thing else, by and through which their Action and Force may be conveyed from one to another, is to me so great an Absurdity, that I believe no Man who has in philosophical Matters a competent Faculty of thinking, can ever fall into it. Gravity must be caused by an Agent acting constantly according to certain Laws; but whether this Agent be material or immaterial, I have left to the Consideration of my Readers.
In our attempt to resolve Newton's action-at-a-distance problem within his domain of physics, we first wish to explore a system where the Gauss's law of gravitostatics (8) and the equation of continuity (14) work peacefully or co-exist simultaneously but now with the restrictions in (15) removed by the following conditions: For this purpose, we introduce the time dependence in (8) by taking its time derivative and then write the result as Now using Eq. (17) in Eq. (14) we obtain The quantity inside the parenthesis of (18) is a vector whose divergence is zero. Since ∇ · (∇ × X) = 0 for any vector X, the vector inside the parenthesis of (18) can be expressed as the curl of some other vector, say h. Mathematically speaking, the Eq. (18) admits of two independent solutions: Except for the sign ambiguity in (19) (to be clear up soon), we now realize that we have just arrived at (in some form) the gravitational analogue of the Ampère-Maxwell law in classical electrodynamics: where H is called the magnetic field, B = μ 0 H is called the magnetic induction field in vacuum, E represents the electric field, j e = ρ e v = the electric charge current density (ρ e being electric charge density), 0 and μ 0 respectively represents the electric permitivity and the magnetic permeability of free space or vacuum and they are related to the speed of electromagnetic wave in vacuum, c , by the relation and c is a universal constant of nature. In fact, there are four field equations in electrodynamics which are collectively known as Maxwell's equations: Maxwell's Eqs. (22a-22d) form the basis of all classical electromagnetic phenomena including the production and transmission of electromagnetic waves carrying energy and momentum even in vacuum. When combined with the Lorentz force equation, and Newton's 2nd law of motion (11), these equations provide a complete description of classical dynamics of interacting charged particles and electromagnetic fields. Do analogous phenomena occur in gravitational physics? For instance, Maxwell's equations predict the existence of electromagnetic waves that travel through vacuum (where ρ e = 0 and j e = 0) at a universal speed c = 3 × 10 8 m/s and the wave equations for the fields E and B in vacuum are obtainable from (22a-22d) as Before we explore analogous wave equations for the E g and h fields, we note that one of the solutions in (19) would correspond to the reality and the other might be a mathematical possibility having no or new physical significance in which we are not presently interested in. We should note that Nature does not behave in two different ways at the same place and time -A physical quantity has to take a unique value at a particular place and time to be real. For this reason and for our present purpose, we have to choose one of the two solutions that may correspond to the reality. Which one to choose? We can answer this question by following the rule of study by analogy. A close look at Maxwell's in-homogeneous Eqs.
(22a, 22b) suggests us: (a) that the source terms (ρ e ) and (j e ) are of the same sign (b) that the analogue of 0 in gravity is 0g : which we may call the gravito-electric permitivity of vacuum or whatever better name one may assign to it, and (c) to introduce a new constant μ 0g as which is the gravitational analogue of μ 0 and c g is a new universal constant for vacuum having the dimension of velocity (which will turn out as the speed of gravitational waves in vacuum, if they exist) such that we get a relation an analogous to (21) as Accepting the suggestion (a), we are naturally led to choose the solution (19) having a negative sign before j 0 as the real solution in view of Eq. (8), where the source term ρ 0 has a negative sign before it. Now, we multiply 8 μ 0g with our chosen Eq. (19) to get where we considered (μ 0g = 4π G/c 2 g ) and defined B g = μ 0g h as the gravitomagnetic induction field. Now taking the curl of (28) and using the vector identity ∇ × (∇ × B) = ∇(∇ · B) − ∇ 2 B, we get In vacuum (where j 0 = 0) Eq. (29) will reduce to a wave equation for the B g field, viz., provided the following two conditions: are fulfilled. Now the taking curl of the Eq. (32) and using (28) we get the following equations for vacuum (where j 0 = 0) Since ρ 0 = 0 = ∇ · E g in vacuum, this equation yields a wave equation for the E g field: Thus, we notice that we have arrived at gravitational wave producing Gravito-Maxwell equations representing what call here as

Non-relativistic Maxwellian gravity (NRMG):
To complete the dynamical picture, we explore the gravitational the analogue of Lorentz force law (23) below adopting the logic and methods of Schwinger et al. [2] in their derivation of the Lorentz force law. Consider two inertial frames S and S having a relative velocity v between them. Let all the masses are in static arrangement in one of these frames, say S , which is moving with velocity v. This way we introduce the time dependence of ρ and E g in the simplest way by assuming all masses are in uniform motion with a common velocity v with respect to Sframe. Here we use the Galileo-Newton principle of relativity (masses at rest and masses with a common velocity viewed by a co-moving observer are physically indistinguishable) and insist that physical laws are the same in the two inertial frames. Further, we assume that |v| << c g = c. To catch up with the moving masses one would have to move with velocity v. Accordingly, the time derivative in the co-moving system, in which the masses are at rest, is the sum of explicit time dependent and co-ordinate dependent contributions, so, in going from static system to uniformly moving system, we make the replacement The equation for constancy of gravitational field in Eq. (15c) becomes, in the moving system The vector identity for constant v is Now using Gauss's law of gravitostatics (35a) in the above identity we get Again, from Eqs. (38) and (40), we get Multiplication of equation (41) by c −2 g gives us Equation (42) will agree with the Eq. (35b) of NRMG, only when Now consider the vector identity for constant v where we have used ∇ ·B g = 0. Moreover, in the co-moving system where the masses are at rest -static -the B g field should also not change with time: From (44)and (45) we get Again, the vector identity gives us the left hand side of the wave Eq. (34) for E g in vacuum (i.e. outside the mass distribution where ∇ ·E g = 0) as By means of the Eq. (35b) of NRMG and Eq. (46), the right side of the wave Eq. (34) becomes (j 0 = 0 outside the charge distribution) Equations (48,49) show that the wave equation for E g (34) will hold if But this cannot be completely correct since as v → 0 ⇒ E g → 0. No gravitostatics ! However, all that is necessary is that the curl of this E g in (50) should be valid: or, if we use Eq. (46), This is consistent with gravitostatics since it generalizes ∇ × E g = 0 to time-dependent situation. Now we are ready to address the question: What replaces the equation F = m 0 E g to describe the force on a point particle of mass m 0 , when that particle moves with some non-relativistic velocity v in given E g and B g field? For this purpose, consider two coordinate systems, one in which the particle is at rest (co-moving coordinate system) and one in which it moves at velocity v. Suppose, in the later coordinate system, the gravitational (or gravito-electric) and gravitomagnetic fields are given by E g and B g , respectively. In the co-moving frame, the force on the particle is where E eff g is the gravitational (or gravito-electric) field in this frame. In transforming to the co-moving frame, all the other masses -those responsible for E g and B g -have been given an additional counter velocity −v. From Eq. (50), we then infer that (+v × B g ) has the character of an additional gravito-electric field in the co-moving frame. Hence, the sug- leading to the gravitational analogue of the Lorentz force law that we call Gravito-Lorentz force law: This gravito-Lorentz force law when used in Newton's 2nd law of motion, we get the following equation of motion of NRMG Since ∇ · B g = 0, B g can be defined as the curl of some vector function, say A g : where A g is the vector potential for NRMG. Using Eq. (58) in the Gravito-Faraday law of NRMG (35d), we get the following expression for the E g in terms of φ g and A g as We now recognize that the equations of NRMG correspond to the Eqs. (1-7) of GRMG-UG provided we set c g = c. In this derivation, c g is an undetermined parameter of the theory, whose value may be obtained from some more advanced theory or from experiments.
Substituting the expression for E g given by Eq. (59) and the expression for B g defined by (58) in the in-homogeneous field Eqs. (35a, 35b) of NRMG, we get the following expressions for the in-homogeneous Eqs. (35a, 35b) in terms of potentials (φ g , A g ) as if the following gravitational Lorenz gauge condition, is imposed. These will determine the generation of gravitational waves by prescribed gravitational mass and mass current distributions. Particular solutions of (60) and (61) in vacuum are where t = t −|r−r |/c g is the retarded time and dv is an elementary volume element at r . Thus, we saw that retardation in gravity is possible in Newtonian space and time in the same procedure as we adopt in electrodynamics. Hence, we have reasons to strongly disagree with Rohrlich's conclusion [19]: "Because the Newtonian theory is entirely static, retardation is not possible until the correction due to deviations from Minkowski space is considered". According to the present field theoretical view, gravitation, like electromagnetism and all other fundamental interactions, acts locally through fields: A mass at one point produces a field, and this field acts on whatever masses with which it comes into contact [20]. Because of the finite propagation speed of the fields, gravitational effects/information propagate at finite speed, that is the cause behind retarded interaction. However, in case of static or quasi-static mass distributions, retardation effects are negligible and hence no distinction can be made between local interaction and action-at-a-distance. By introducing the concept of physical fields carring energy and momentum, 9 one can address Newton's "action-at-a-distance" problem within Newton's world of physics by extending his field equations to time-dependent fields, sources and searching for the conditions for the existence of gravitational waves in free space traveling at a finite speed. This is what Heaviside had done in 1893 [4,5]. May be Einstein had not seen Heaviside's field equations when he was working on his relativistic theory of gravity. Had Einstein seen Heaviside's field equations, his remark on Newton's theory of gravity would have been different than what he made before the 1913 congress of natural scientists in Vienna [21],viz., After the un-tenability of the theory of action at distance had thus been proved in the domain of electrodynamics, confidence in the correctness of Newton's action-at-adistance theory of gravitation was shaken. One had to believe that Newton's law of gravity could not embrace the phenomena of gravity in their entirety, any more than Coulomb's law of electrostatics embraced the theory of electromagnetic processes.

Special relativistic Maxwellian gravity (SRMG)
With the establishment of special relativity (SR) theory and the equivalence of mass and energy, the meaning of the inertial mass and gravitational mass became ambiguous, because SR suggests two inertial mass-energy concepts: (1) the Lorentz invariant rest-mass m 0 = E 0 /c 2 (E 0 = restenergy, which is the sum total of all forms of energy in the rest frame of a body or particle) and (2) the mass attributed to the relativistic energy m = E/c 2 (E = sum of all forms of energy at rest and motion) which is not Lorentz-invariant. The qualitative distinction that existed between two inertial mass concepts in Newtonian mechanics became quantitatively distinct and clear in SR. Now, one fundamental question arises, "What form of mass (or energy) should be the source of gravity (m g ) in a relativistic version of Newtonian gravity?" In any construction of a field theory of gravity compatible with SR and the correspondence principle by which a relativistic theory gravity is reducible to Newtonian gravity, a decision on which form of "mass" (or energy) is the source of gravity has to be taken. Such a decision, as Price [22] has rightly pointed out, will be crucial not only to the resolution of the ambiguity mentioned above but also to the issue of the nonlinear nature of gravity. One of the Eddington's [23] four reasons to feel dissatisfied with Newton's Law of gravitation is appropriate here to quote: The most serious objection against the Newtonian Law as an exact law was that it had become ambiguous. The law refers to the product of the masses of the two bodies; but the mass depends on the velocity-a fact unknown in Newton's days. Are we to take the variable mass, or the mass reduced to rest? Perhaps a learned judge, interpreting Newton's statement like a last will and testament, could give a decision; but that is scarcely the way to settle an important point in scientific theory.
In his construction of a relativistic theory of gravity popularly known as General Relativity (GR), Einstein has taken a decision in favor of the equality of m with m g . For a theoretical justification of this decision, Einstein by writing Newton's equation of motion in a gravitational field (in our present mathematical notation) as (wrongly!) inferred from it [24]: It is only when there is numerical equality between the inertial and gravitational mass that the acceleration is independent of the nature of the body.
This inference is often expressed in one of the two ways: (S1) that the particle's motion is mass independent, or (S2) that the particle's inertial mass m = its gravitational mass m g .
The two statements (S1) and (S2) are sometimes used interchangeably as the weak equivalence principle (WEP) in the literature [25][26][27]. This use of terminology is rather confusing, as the two statements are logically independent [28]. To explore this possibility, to get new insights for making Newtonian gravity compatible with the SR, to regard old problems from a new angle, we re-examined [3] an often cited [29,30] Salisbury-Menzel's [31] 10 thought experiment (SMTE) from a new perspective as discussed in the following subsection. Before that the author would like to remark that perhaps Einstein, himself, was not satisfied with his above inference of m g = m, as we can sense from his another statement on the equality of m g with m [32,33]: The proportionality between the inertial and gravitational masses holds for all bodies without exception, with the (experimental) accuracy achieved thus far, so 10 The correct special relativistic result is obtained by P. Lorrain that we may assume its general validity until proved otherwise.
The last three words of Einstein's above statement,'until proved otherwise', show that he was very cautious and not very confident of what he was stating. Based on the experimental results available up to 1993, Mashhoon [30] noted that the observational evidence for the principle of equivalence of gravitational and inertial masses is not yet precise enough to reflect the wave nature of matter and radiation in their interactions with gravity (see other references on equivalence principle in [3,30]).
3.1 Re-examination of SMTE to show m 0 = m g Consider a system of two non-spinning point-like charged particles with charges q 1 and q 2 with respective rest masses m 01 (= E 01 /c 2 ) and m 02 (= E 02 /c 2 ) such that they are at rest in an inertial frame S under equilibrium condition due to a mutual balance of the force of Coulombic repulsion (F C ) and the Newtonian gravitostatic attraction (F N ) between them. Our aim is to investigate the condition of equilibrium of this two-particle system (realizable in a Laboratory by taking two perfectly identical spherical metallic spheres having requisite masses and charges so that they are in equilibrium) and in different inertial frames in relative motion. For our re-examination purpose, suppose that the particles are positively charged and they are in empty space. Let the particle No. 2 be positioned at the origin of S -frame and r be the position vector of the particle No. 1 with respect to the particle No. 2. In this S -frame the condition of equilibrium is fulfilled by where r = |r | and other symbols have their usual meanings. From (66) we get Equation (67) represents the condition of equilibrium, in terms of the charges and rest masses (or rest energies) of the particles, under which an equilibrium can be ensured in the S -frame. For example, if each metallic sphere is given a charge of 1 × 10 −6 Coulomb, then the rest mass of each sphere should be 1.162 × 10 4 kg, to fulfill the equilibrium condition (67) in a laboratory experiment. Now, let us investigate the problem of equilibrium of the said particle system from the point of view of an observer in another inertial frame S, in uniform relative motion with respect to the S -frame. To simplify the investigation, let the relative velocity v of S and S -frame be along a common X/ X -axis with corresponding planes parallel as usual. Since the particles are at rest in S -frame, both of them have the same uniform velocity v relative to the S-frame. Let the position vector of the particle No. 1 with respect to the particle No. 2 as observed in the S-frame be r and the angle between v and r be θ .
For an observer in the S-frame, the force of electric origin on either particle (say on particle No. 1 due to particle No. 2) is no more simply a Coulomb force, but a Lorentz force, viz., where What about the force of gravitational interaction as observed in the S-frame? It can not simply be a Newtonian force but something else, otherwise the particle system will not remain in equilibrium in the S-frame. Such a situation will amount to a violation of the principle of relativity in special relativity. A null force should remain null in all inertial frames. Therefore, a new force law of gravity has to be invoked so that the equilibrium is maintained in accordance with the principle of relativity (Lorentz invariance of physical laws). Let this new unknown force be represented by F gL such that the equilibrium condition in S-frame is satisfied as: Taking into account the Eqs, (68)-(71), F gL in Eq. (72) can be expressed as: Now, using Eq. (67), we can eliminate q 1 q 2 from Eq. (73) to get the expression for F gL in terms of m 01 , m 02 and G as: By comparing the quantities in Eq. (75) with those in Eqs. (25)(26)(27), we immediately find that Now, (74) may be rearranged to the following form to represent the Gravito-Lorentz force law of SRMG: where Equations (77-79) are in complete formal analogy with the Eqs. (68-70) of classical electromagnetism in its relativistic version. Thus, from the requirement of the frameindependence of the equilibrium condition, we not only obtained a gravitational analogue of the Lorentz-force law expressed by Eq. (77) but also unexpectedly found the Lorentz-invariant rest mass as the gravitational analogue of the electric charge by electromagnetic analogy. From this analysis, the gravitational charge (or rest mass) invariance may be interpreted as a consequence of the Lorentzinvariance of the physical laws. These findings are in conformity with Poincare's [34] remark that if equilibrium is to be a frame-independent condition, it is necessary for all forces of non-electromagnetic origin to have precisely the same transformation law as that of the Lorentz-force. Now, following Rosser's [35] approach to classical electromagnetism via relativity one can obtain the field equations of SRMG as represented in Eqs. (35a-35d) with c g = c, from the Eqs. (77)(78)(79). Alternatively, after recognizing our new findings from the above thought experiment, especially m g = m 0 and c g = c, one may follow the procedure we followed in NRMG to arrive at desired gravitational wave producing field equations. It is to be noted that the Lorrain's (See Lorrain in footnote 10 of this paper) exact special relativistic derivation of gravitational analogue of the magnetic force from SMTE matches with our SRMG results.
The relativistic equation of motion (83) is independent of the mass of the particle moving in an external gravitoelectromagnetic (GEM) field f αβ . Thus we saw that the motion of a particle in an external GEM field can be independent of its mass without any postulation on the equality of gravitational mass with frame-dependent inertial mass. Equation (83) is the relativistic generalization of Galileo's law of Universality of Free Fall (UFF) expressed through the non-relativistic equations of motion (13) and (57) and known to be true both theoretically and experimentally since Galileo's time. Now, if we introduce the energy momentum four vector: where p 0 = E/c and U α is the 4-velocity, then with this p α we can re-write equation (83) as Thus, in SRMG the fields f αβ couple to the energymomentum 4-vector of all particles of whatever rest masses they have, provided m g = m 0 holds exactly. It is to be noted that the equation of motion (83 or 87) holds only in an inertial frame. Appropriate modifications are necessary for its application in non-inertial frames, as is done in non-relativistic physics by introducing pseudo-forces.

Original analysis of SMTE with assumption of
In the original analysis of SMTE [31] Salisbury and Menzel (SM) axiomatically used flat space-time and assumed m g = m = m 0 / 1 − v 2 /c 2 for their thought experimental demonstration of gravito-magnetic field (they called it Gyron field) and the gravitational analogue of Lorentz force law. From the analysis of their results, one can find that in the slow motion approximation, if the gravito-Lorentz force law is written in the following form which yields On the other hand if one considers c SM g = c, then equation (88) has to be written in the following form: We designate this type of gravity as linearized version of nonlinear SRMG (SRGM-N) in flat space-time. The origins of the non-linearity of (SRGM-N), the appearance of the spurious value of c g = c/ √ 2 or a factor of "2" in the gravitomagnetic force term (due to a supposed value of c g = c) are all now traced to the adoption of Einstein's doubtful postulate on the equality of gravitational mass with the velocity dependent inertial mass. As such, SRMG-N does not correspond either to the NRMG when c g = c or to the non-relativistic limit of SRMG.

Maxwellian Gravity of others without GR
Sciama [36], in 1953, hypothetically adopted SRMG (assuming m g = m 0 ) to explain the origin of inertia, calling it a toy model theory of gravity which differs from GR principally in three respects: (a) It enables the amount of matter in the universe to be estimated from a knowledge of the gravitational constant, (b) The principle of equivalence is a consequence of the theory, not an initial axiom and (c) It implies that gravitation must be attractive. However, he concluded his paper mentioning three limitations of such a theory: (i) It is incomplete because the relativistic form of Newton's law must be derived from a tensor potential, 11 not from a vector potential, (ii) It is difficult to give a consistent relativistic discussion of the structure of the universe as a whole 12 and (iii) It is also difficult to describe the motion of light in a gravitational field. 13 Carstoiu [37][38][39], in 1969, rediscovered Heaviside's gravitational Eqs. (35a-35d) (in our present notation as per the report of Brilloiun [39]) assuming the existence of a second gravitational field called gravitational vortex (here called gravito-magnetic field) and assumed c g = c by electromagnetic analogy [39]. In 1980, Cattani [40] considered linear equations for the gravitational field by introducing a new field by calling it the Heavisidian field which depends on the velocities of gravitational charges in the same way as a magnetic field depends on the velocities of electric charges and shown that a gravitational field may be written with linear co-variant equations in the same way as for the electromagnetic field. Cattani's equations differ from some important formulae of general relativity such as the gravitational radiation, Coriolis force by a factor of 4. In 1982, Singh [41] considered a vector gravitational theory having formal symmetry with the electromagnetic theory and explained the (a) precession of the perihelion of a planet (b) bending of light in the gravitational field and (c) gravitational red-shift by postulating the self-interaction between a particle velocity and its vector potential. In 2004, Flanders and Japaridze [42] axiomatically used the field equations of SRMG (albeit without reference to [3]) and special relativity to explain the photon deflection and perihelion advance of Mercury in the gravitational field of the Sun. Borodikhin [43] explained the perihelion advance of Mercury, gravitational deflection of light as well as Shapiro time delay by postulating a vector theory of gravity in flat space-time that is nothing but SRMG. Borodikhin also showed that in a vector theory of gravity, there exists a model for an expanding Universe. Jefimenko [44,45] also deduced the equations of NRMG by extending Newton's gravitational theory to timedependent sources and fields and using the causality principle. Jefimenko assumed c g = c and postulated a gravito-Lorentz force. Recently, Heras [46], by recognizing the general validity of the axiomatic approach to Maxwell's equations of electromagnetic theory, used those axioms to derive only the field equations (leaving out gravito-Lorentz force law) of SRMG, where the in-variance of gravitational charge is considered. Other recent derivations of SRMG equations from different approaches include the works of Nyambuya [47], Sattinger [48], Vieira and Brentan [49]. The historical objections of several researchers, starting from Maxwell (see footnote 2 of this paper) upto Misner, Thorne and Wheeler (MTW, Sect. 7.2) [54], concerning negative energy density of gravitational field ('Maxwell's Enigma' as Sattinger puts it) in a linear Lorentz invariant field theory of gravity are also refuted by Sattinger [48], who considered negative field energy density for SRMG in agreement with our result [3]. In the discussion on the Dark Matter problem, Sattinger further noted: "The Maxwell-Heaviside equations of gravitation constitute a linear, relativistic correction to Newton's equations of motion; they interpolate between Newton's and Einstein's theories of gravitation, and are therefore a natural mathematical model on which to build a dynamical theory of galactic structures". Our unique, important and new findings reported and discussed here, again confirmed by our very recent work on "Attractive Heaviside-Maxwellian (vector) Gravity from Quantum Field Theory" [50] (where gravitational energy density for free fields is fixed positive by choice to address the objection of MTW (Sect. 7.2) [54] without any inconsistency with the field Eqs. (81-83) of SRMG), corroborate all the above suggested or derived linear vector gravitational equations in flat space-time which are seen to satisfy the correspondence principle (cp) in its correct sense: Newtonian Gravity ⇔ NRMG ⇔ SRMG. This means that the field equations of SRMG have room for both positive and negative energy solutions -a discussion of which is left out here.

Maxwellian gravity From GR (GRMG)
Different Maxwell-Lorentz-type equations for gravity obtained from GR by different researchers following different linearization procedures are not isomorphic [51] as seen below: some contain non-linear terms and do not satisfy the cp from the perspective of our present new NRGM results substantiated by the cp-respecting SRMG results. Some samples of this type of General Relativistic Maxwellian Gravity (GRMG) are listed below for discussion. the the non-relativistic Gravito-Lorentz force of NRMG or of the SRMG. Further equation of continuity does not follow from GRMG-BF field equations because of the existence of some non-linear terms in Eq. (92a).

GRMG of Harris (GRMG-H)
Instead of using the PPN formalism of GRMG-BF, Harris [55] derived a set of gravitational equations for slowly moving particles in weak gravitational fields starting from the equations of GR. The resulting equations have some resemblance to those in electromagnetism: The Gravito-Lorentz force equation of Harris is of the following form The field E g is related to the scalar potential (φ g ) and vector potential (A g ) as In empty space (where ρ 0 = 0), the field equations (100a-100d) give us the following wave equations for the fields (E g , B g ): where c g = c/ √ 2. So GRMG-H is also defective like GRMG-BF.
Since the divergence of curl of any vector is identically zero, the divergence of Eq. (105d) gives us ∇ · j 0 = 0. Thus Eq. (105d) has the same limitation as that of the Ampère's law of electromagnetism. Therefore, it needs a correction like Maxwell's correction to the Amp'ere's law. While the conditions ∇ · j 0 = 0 and ∂ρ 0 /∂t = 0 are valid for steadystate problems, the general situation, where ∇ · j 0 = 0 and ∂ρ 0 /∂t = 0, is given by the continuity equation for mass and mass current or momentum density: The Eqs. (105a) and (105d) will be consistent with the continuity Eq. (106), if we make the following Maxwell-like correction to the gravito-Ampère law (105d): Further, without the above correction to the Eq. (105d), there can not be gravitational waves. Now, the corrected selfconsistent field Eqs. (105a-105c, 107) yield transverse gravitational waves; the wave equations for the g and b fields of GRMG-OR, in vacuum, take the following forms: where c g = c/ √ 2. This means that the transverse gravitational waves originating from slowly varying fields and weak sources travel through vacuum not at the light speed but at a reduced speed c g = c/ √ 2. Thus GRMG-OR will have exact correspondence with NRMG provided c g = c/ √ 2. However, GRMG-OR will not correspond to the slow motion approximation of SRMG where c g = c in vacuum even at relativistic motion of the fields and sources. Now, if we define a new filed for GRMG-OR as then the Gravito-Lorentz force law of GRMG-OR will take the form: and the field Eqs. (105a-105c, 107) will take the form Again these equations yield wave equations for the g andb fields for which c g = c/ √ 2. Further, if we define B g = b/4, the redefined field equations (not written here) yield wave equations for the g and B g fields for which c g = c/ √ 2 in vacuum, while the gravito-Lorentz force law of GRMG-OR takes the following form The defects of GRMG-OR are apparent from the perspectives of cp violation and the suspected value of c g in Einstein's linearized equations.

GRMG of Pascual-Sànchez, Moore, (GRMG-PS-M):
Following Huei [56], Wald [57] and Ohanian and Ruffini [58], Pascual-Sànchez [59], by using some approximation of GR, obtained the following set of Lorentz-Maxwell-like gravitomagnetic equations that match with Moore's [60] equations from GR in the the weak field and slow motion approximation: Although the field Eqs. (114a-114d) match with those of GRMG-UG [1], SRMG where c g = c exactly and with NRMG equations conditionally when c g = c as shwon here, the Gravito-Lorentz force, containing a factor of 4 in the gravitomagnetic interaction, does not satisfy the cp in the sense discussed here. It is to be noted that Ciubotariu [61] considered Peng's [62] version of GRMG equations 16 in the prediction of absorption of gravitational waves, while Minter et al. [63] considered GRMG-PS-M version in their investigation on the question of the existence of mirrors for gravitational waves.

Maxwellian Gravity of Mashhoon (MG-Mashhoon)
Mashhoon [64,65], in his general linear solution of the gravitational field equations of Einstein, obtained following gravitational analogues of Maxwell's equations: By defining E g and B g fields in terms of scalar and vector potentials (φ g , A g ) as he wrote the Lagrangian, L, for the motion of a test particle of rest mass m 0 (to linear order in φ g and A g ) as where γ = (1 − v 2 /c 2 ) −1/2 is the Lorentz factor. The equation of motion, dp/dt = F, where p = γ mv is the kinetic momentum, is expressed as if ∂A g /∂t = 0 and F is expressed to lowest order in v/c, φ g and A g . In empty space, the field Eqs. (115a-115d) give us wave equations for E g and B g fields with c g = c. But because of the factor of 2 appearing in the relativistic Gravito-Lorentz force law (118), it does not correspond to the cp respecting special relativistic Gravito-Lorentz force law of SRMG expressed in equations (77,83).
written here in Mashhoon's convention for the gravitostatics Gauss's law (115a) and Gravito-Lorentz force Law (118). Further, in the non-relativistic case also, the Eq. (118) does not match with our non-relativistic result here, if c g = c.
In our above discussions, we found that the speed of gravitational waves in different linearized versions of Einstein's equations is not unique; it depends on the thought of a scientist concerning the mangagement of the spurious factor of 4 by splitting it 4 = 2 × 2 and moving a factor of 2 to some other place of the equations for consistency. 17 Eddington [66] in his textbook "The Mathematical Theory of Relativity, which Einstein suggested was "the finest presentation of the subject in any language", rightly found and said that weak-field solutions of the wave equation obtained from Einstein's field equations were just coordinate changes which we can "propagate" with the speed of thought. Apart from this defect of yielding a non-unique value of c g , the linearized theory of gravitational waves has its limits because the linear approximation is not valid for sources where gravitational self-energy is not negligible [67], as in the case of merging of highly compact objects like the Neutron stars and the 'so called' Black Holes. It is important to note that the current experimental data on observation of gravitational waves [6][7][8][9] and grvitomagnetic phenomena are being (or may be) explained by using any one of these cp-defying linearized versions of GR. Further, these generally perceived general relativistic phenomena are being interpreted as having no counterpart in the Newtonian world, which we found not to be true and satisfactory. The author wishes to stress that just by connecting the Gauss Law of gravitostatics and equation of continuity in a consistent way, one gets gravitational analogue of Ampère-Maxwell law, which is being hailed as one of the most important predictions of GR (where the mathematical trees seem to obscure the physical forest). Any talk of gravitomagnetism has long been the prerogative of general relativists. With our present report, even undergraduates, untrained in the mathematical gymnastics of general relativity, can now talk and think of the generation, transmission and detection of gravitational waves [68][69][70] 18 akin to electromagnetic waves as matter of scientific curiosity and 17 Such a scheme my be made general by thinking 4 = 2x × 2 x with x being any non-zero real number to get other sorts of spurious results in a self-consistent way. 18 Which, in the framework of gravito-electromagnetic theory, cannot be just any other wave, but the waves whose nature is dictated by equations of gravitodynamics: Gravito-Maxwell Equations. By contrast, the investigate the role of gravitoelectromagnetism in different fields of study, ranging from classical physics to quantum physics [1,3] and to cosmology [36,43] even quantum cosmology to test the validity (or the domain of validity) of the proposed gravitoelectromagnetic theory from experimental point of view.
Basically, we offer a new pack of beautiful cards (or a beautifully simple, self-consistent toy model theory of gravity) to play with, if one likes. However, it remains to be seen how the recently observed phenomena concerning gravitomagnetism and gravitational waves may be interpreted with our new findings with c g = c, since these phenomena could well be explained by NRMG just by adjusting the undetermined parameter to c g = c/2 or c/ √ 2 or by cp-violating SRMG without invoking the curved space-time concept of Einstein and keeping other experimental parameters intact or by SRMG (c g = c with normal gravito-Lorentz force) and varying other parameters of the experiment in a self consistent way, or re-analyzing the sources of possible theoretical and experimental errors in the interpretation of the experimental data. The author aims to address these questions in future, if his odd situation permits and wishes the young minds do this, if they can.

Spin-1 vector gravity vs spin-2 tensor gravity
Many quantum field theorists, like Gupta [71], Feynman [72], 19 Zee [73] and Gasperini [74], 20 to name a few, have rejected spin-1 vector theory of gravity on the ground that Footnote 18 continued general description of gravitational waves in GR is different [67]. However, in spite of the interesting work by Mead [68], attempting to provide a non-general relativistic explanation of the generation and detection of gravitational waves, which Isi et al. [69] referred to, a proper description of gravitational waves within a purely gravitoelectromagnetic formalism is still far from being established. 19 On page 30 of ref. [72], one finds: A spin-1 theory would be essentially the same as electrodynamics. There is nothing to forbid the existence of two spin-1 fields, but gravity can't be one of them, because one consequence of the spin 1 is that likes repel, and un-likes attract. This is in fact a property of all odd-spin theories; conversely, it is also found that even spins lead to attractive forces, so that we need to consider only spins 0 and 2, and perhaps 4 if 2 fails; there is no need to work out the more complicated theories until the simpler ones are found inadequate. 20 On page 27, Gasperini noted: A correct description of gravity in the relativistic regime thus requires an appropriate generalization of Newtons theory. Which kind of generalization? A natural answer seems to be suggested by the close formal analogy existing between the Newton force among static masses and the Coulomb electrostatic force among electric charges. In the same way as the Coulomb potential corresponds to the fourth component of the electromagnetic vector potential, the Newton potential might correspond to the component of a four-vector, and the relativistic gravitational interaction might be represented by an appropriate vector field, in close analogy with the electromagnetic theory. Such an attractive speculation, however, has to be immediately discarded for a very simple reason: vector-like interactions produce repulsive static interactions between sources of the same sign, while -if gravitation is described by a vector field theory like Maxwell's electromagnetic theory as discussed here, then vector-like interactions will produce repulsive static interactions between sources of the same sign, while -according to Newton's gravitational theory -the static gravitational interaction between masses of the same sign is attractive. However, in one of our recent work [50] on quantum field theoretical rediscovery of Heaviside-Maxwellian gravity (in flat space-time) (same as SRMG), we have shown that this not true because of the existence of a fundamental difference in the sign before the source term in the in-homogeneous field equations of SRMG and relativistic Maxwell's electromagnetism (RMEM) as seen below in SI units: where = ∂ μ ∂ μ = 1 c 2 ∂ 2 ∂t 2 − ∇ 2 , A μ g and j μ 0 respectively represents the 4-potential and 4-mass current density of SRMG while A μ and j μ e respectively represents the 4potential and 4-charge current density of RMEM. This fact that a vector gravitational theory, as proposed first by Heaviside and later rediscovered by many others following different appoaches, implies attractive interaction between static masses was transparently clear to Sciama [36]. In spin-1 SRMG, contrary to the electromagnetic cases, like masses (or gravitational charges) should attract and unlike masses (if they exist) should repel each other under static conditions, while like (parallel) mass currents repel and un-like mass currents attract each other [3]. Similarly, there should be attraction between like gravitomagnetic poles and repulsion between un-like gravitomagnetic poles: opposite to the case in electromagnetism where like magnetic poles repel and un-like magnetic poles attract each other [3]. Following Dirac's scheme of predicting the spin magnetic moment of an electron, in our previous work [3], we have shown that the gravitomagnetic moment of a Dirac (spin 1/2) fermion is exactly equal to its spin angular momentum: μ sg =h 2 , which can just be inferred from the magnetic moment of a Dirac (spin 1/2) fermion, μ se = qh 2m 0 , by replacing the electric charge q with the gravitational charge m 0 of the Dirac fermion as per SRMG. However, in GR, μ sg value is not unique as found by different authors and noted in [3].
Regarding the idea of spin-2 graviton, Wald [57] noted that the linearized Einstein's equations in vacuum are precisely the equations written down by Fierz and Pauli [75], in 1939, to describe a massless spin-2 field propagating in flat space-time. Thus, in the linear approximation, general relativity reduces to the theory of a massless spin-2 field which undergoes a non-linear self-interaction. It should be noted, Footnote 20 continued as is well known -the static gravitational interaction between masses of the same sign is attractive. however, that the notion of the mass and spin of a field require the presence of a flat back ground metric η ab which one has in the linear approximation but not in the full theory, so the statement that, in general relativity, gravity is treated as a mass-less spin-2 field is not one that can be given precise meaning outside the context of the linear approximation [57]. Even in the context linear approximation, the original idea of spin-2 graviton gets obscured due to the several faces of Gravito-Maxwell equations seen here. This may be seen as another limitation of GR for not making a unique and unambiguous prediction on the spin of graviton.

Little known Heaviside's work on gravity
Heaviside's work on gravity, which McDonald [4] called a low velocity, weak-field approximation to general relativity, is little known and has not received as much attention as it deserves. This is because, in many leading papers and books exploring gravitomagnetic phenomena and gravitational waves, one finds rare or no mention of Heaviside's name, although Heaviside [5] predicted gravitomagnetic effects and considered the necessity of possible existence of gravitational waves (for which there must be some wave equations first written down by him) almost 20 years before Einstein's prediction of gravitational waves [76,77]. Further, it is more surprising not to find Heaviside's name in the Scientific Background on the Nobel Prize in Physics 2017 [78] when the 2017 Physics Nobel Prize was declared to be awarded to Rainer Weiss, Barry C. Barish and Kip S. Thorne for (their) decisive contributions to the LIGO detector and the observation of gravitational waves. Brillouin [39], in his final remark on Carstoiu's [37,38] suggestions for gravity waves aptly stated, "It is very strange that such an important paper had been practically ignored for so many years, but the reader may remember that Heaviside was the forgotten genius of physics, abandoned by everybody except a few faithful friends." version and named it Special Relativistic Maxwellian Gravity