Gravitational field is one of the most important fields created by Nature which has (so far) pure attractive effect among anything with mass or energy   , and its associated force is radial and inverse proportional to the square of the separation of the objects. For an arbitrary mass density distribution which is not depending explicitly on time, it is known that the gravitational field in the space can be found by solving Poisson equation with proper boundary conditions  . However, Einstein  saw that gravity can be seen as the component of a pseudo-metric defined in a space-time manifold, where the shape of this manifold is due to the mass and energy contained in it  . This was called the General Relativity Theory and had enormous success explaining and predicting many gravitational phenomena in Nature    . However, this classical theory has had serious difficulties to construct a more fundamental theory (as many thinks) which is related with the quantization of gravity (the union of General Relativity with Quantum Mechanics  ), where a graviton appears as a quantum particle associated to the interaction with gravity which has spun two   . This direction of research has its own value, and there are many scientists working on this topic.
In addition, a big effort has been made to unify gravity with electrodynamics  -  , but still a consistent mathematical and geometrical theory is missing, alternatively, using perturbation theory on GR   , or the decomposition of the Weyl and Maxwell tensors in electric and magnetic parts  -  . The so called Gravito-Electromagnetism (GEM) emerges as a result of these approaches, which is the expression of the equations to describe the gravitational field as a Maxwell-like type of equations. The similarity of gravity and electrodynamics equations appeared long time ago when people realized that, in the static case, both fields satisfy Poisson’s equations. In fact, Heaviside   was one of the first persons who saw the similarities between Electromagnetism and Gravity through Poisson’s equation, and in addition  , for a current of mass motion as continuous fluid ( ), he postulated a gravity field h having a curl analog to magnetic field in Maxwell’s equations for Electrodynamics, in order to get a Maxwell-like equations for gravity, which, of course lead us to gravitational waves. Since his approach lacks of mathematical justification, it cannot be considered as a deduction, and it is done for a very particular case (fluid), the idea settled there was that Maxwell-like equations for gravity could be arisen without GR theory. In a similar axiomatic form, Yaroslav  has recently show the deduction of Maxwell-like equations for gravity without GR, where perihelion of Mercury and bending of light by massive object are presented. It also worths to mention works based on Hodge theory and differential forms  , where a similarity with the electromagnetic permittivity and permeability was found. The idea of having a magnetic-like gravitational field in GR was reinforced due to Lense-Thirring effect   , where the angular momentum of a rotating body can be interpreted as a gravitational magnetic field, at large distances. Then, one could have the conclusion that GR is needed in order to see a gravitational field like GEM.
On this paper, one will be focused on other classical aspects of the gravitational field, to show that a full Maxwell-like equations for gravitational field arises as a consequence of the continuity equation for the mass density and the usual Gauss’ theorem.
2. Maxwell-Like Equation for Gravity
For a time explicitly depending density of mass, , where the following continuity equation is satisfied
with J being the mass current density, which for a fluid-like system is given by , being v the fluid velocity. This continuity equation brings about the total mass conservation of the system, . Now, it is well known  that if the law of forces between to mass elements is of the form , where and are the vector position of the two elements, and is their distance separation, then it follows that the Gauss’ law for the gravitational field R in a compact region of the space is
where G is the constant of gravitation ( ). Differentiating this expression with respect the time, and using (1) in the region W, one has
This equation implies that there must exist a vector field W such that
Of course, there is a well known indetermination of this expression that can be absorbed in the definition of R or W. It is necessary to mention that up to know non experiment has been made to try to find the W gravitational field (something similar when one has a current in a wire and see what happen to the motion of a charge nearby), and this is the gravitational field so called h by Heaviside  , defined on purpose to obtain Maxwell-like equations within his approach. So far, one knows and of the vector fields R and W. To know and , let us take the curl of (4) and use the known expression for any vector field a,
To decouple equations for R and W, one can choose the following relations
Denoting the constant of proportionality of (6a) as , one has for (5) the following expression
Now, taking the curl of (6b), it follows that
which can be written as
Equations (27) and (28) represent non homogeneous wave equations with known sources. The wave propagates with speed “l”, and by selecting on (6a), the equations which determine the gravitational fields R and W have a Maxwell-like form
Thus, one may say in summary that if the Gauss’ theorem is valid for the gravitational field R, and the continuity equation is satisfied for r and J, then, there must exist a gravitational vector field W such that W and R satisfy inhomogeneous wave equations which propagate with an unknown speed “l”. One must mention that there has not been a direct measurement of the speed of gravitational waves, but it is assume that gravity also travel to the speed of light  . In addition, none experiment has been so far trying to find the gravitational field W, and general controversy about Maxwell equations for gravitational field can be found in reference  .
Now, in terms of the vector and scalar potentials, and , which is well known where they come from, the gravity fields can be written as
Of course, and identical relation with the electrodynamics Maxwell’s equation can be gotten by choosing (11) and (13) as
If one applies the same initial hypothesis to the electric field in the case of electrodynamics, one see that existence of the magnetic field and the displacement current must appear as a direct consequence of the continuity equation for the density of charge and current. The decouple of the resulting equations would bring about the usual Maxwell’s equations and the non homogeneous wave equations for the electric and magnetic fields, with the waves traveling at the speed of light.
3. Lorentz-Like Force
Firstly, one needs to say that in the MKS system of units, one has the following units of our quantities and vectors
Now, given the gravitational field defined in the space-time, the force felt by and object of mass M located at the point x at the time “t” due to the field R is given by
which is well known. However, given the field in the space-time, the only experimental evidence that this field could be associated to a magnet-like force is the Lense-Thirring effect  . In this way, the force associated to the gravitational field W would be
If one considers that these are true, one gets a Lorentz-like force expression for an object of mass M under the gravitational interaction with the fields R and W,
where v is the velocity of the object.
4. Energy Density, Poynting Vector, Stress Matrix
Let us write the Maxwell-like equations for the gravitational fields R and W (MKS units) as
where the resulting decoupled equations for R and W are the inhomogeneous wave equations
and write the Lorentz-like law of force, acting over an object of mass m, as
where v is the velocity of the object. The units of the vector fields R and W are the same, m/s2. As pointed out before, the field W, the speed l, and the force still require direct experimental verification.
Now, let us proceed as one does in Electrodynamics Theory  . From (22), one knows that there exists a vector field A (vector potential) such that
Substituting this in (21), it follows that there exists a function F (scalar potential) such that
Using (27) and (28) in the expressions (20) and (23), and choosing the Lorentz-like gauge
one gets decoupled equations for the potentials F and A which also satisfy inhomogeneous wave equations
The power density is related with the mechanical energy density as
which can be written using (23) and (21) as
bringing about the relation
where S and are the Poynting vector and the gravitational energy density,
As one can see, the concept of gravitational energy is well defined, and the Poynting vector S has the usual meaning of flux of density energy per unit area per unit time. For pure gravitational field ( ), one obtains the continuity equation
implying that the gravitational energy velocity is
In vacuum ( and ), a solution of (24) and (25) are plane gravitational waves, , where k and w are the vector wave number and frequency of gravitational wave, and one has the properties: , , and . From (37) it follows that the velocity of the gravitational plane wave is .
Now, assuming the current density of mass as a fluid in motion, where every element of the fluid is moving with velocity v, the current density is determined by , and from the expression (26), the density of force at any point on the fluid is
This expression can be written using (20) and (23) as
where is the density of mechanical linear momentum of the system. After adding the zero parts and and doing some manipulations, this expression can be written as
where and T are the gravitational linear momentum and tensor stress,
being Ä the tensorial product, and I being the -identity matrix. So, the gravitational field has a well defined linear momentum and stress tensor.
5. Formulation in Space-Time Â4
Consider the following vectors in Â4 with components
Consider also the following vector fields defined on Â4 (space-time) with components
and the anti-symmetric matrix defined by the following matrix elements
or, from the expressions (27) and (28), this matrix is explicitly given by
Using the operators
and the operator (assuming Einstein summation convention)
one can have the following identifications
In fact, the last expression follows necessarily due to anti-symmetric definition of F. Wave Equations (24), (25), (30a) and (30b) can be written as
This formulation is convenient to see the change of form of these expression under a change of inertial reference system.
6. Change of Inertial Reference System
Let S and S’ be two inertial reference systems which at have a common origin, and S’ is moving with respect to S with a constant velocity V, the reference system S is fixed with respect the far away Galaxies in our Universe. Let us defined as
the normalized velocity of S’. Then, there must be a nonsingular matrix such that the space-time of the system S can be transformed on the space-time of the system S’,
In other words, the set has the structure of group. So, the transformation of the fields defined in that one would expect are of the form
It is important to point out that this transformation is not necessarily the know Poincaré-Lorentz transformation since the mass is not a Poincaré-Lorentz invariant. Of course, (48)-(54), and (55) are invariant of form (they transform covariantly) under a change of an inertial reference system, although this transformation has not given explicitly. In addition, one has that
that is, is an invariant under these transformations. Therefore, one can define the Lagrangian density for the gravitational field as
and it is known from Electrodynamics  that Euler-Lagrange equations applied to this Lagrangian leads us to Maxwell-like Equations (51) and (53).
This group of transformations also will leave some metric (defined in ) invariant. If this metric is given in the inertial reference system S by
where , is the transpose vector, and g is a real constant nonsingular matrix. In the inertial reference system S’, this metric is given by
So, one will have invariance of the metric ( ) if it follows that
Then, one chooses as the reference systems S’ that one which goes with the object ( ), and the time defined on this system (even it is not really an inertial system) is called “proper time” t. Thus, one gets
Factorizing dt from the right hand side, one can obtain a relation between the proper time and the time measured in the system S of the form
where is the resulting function of this factorization, and is related with the velocity of the object (seen from the system S), . The equations
of motion invariant under the set of transformations must be of similar form as in special theory of relativity    ,
This formulation will be useful later on when radiation reaction force for gravitational field be studied. For our purposes, it is not necessary to know the set of transformations neither the function , if , however, one assumes that Poincaré-Lorentz-like transformations can be used here, the velocity of the inertial system S’ is (with ), and the normalized velocity of the object is (with ), one would have
7. Gravitational Waves Emission
It is known that the particular solution of the Equations (30a) and (30b) are given by the convolution of the inhomogeneity with the fundamental solution of the wave equations  ,
where is the fundamental solution,
This brings about the known retarded potentials
where is the domain where r and J are defined. For a point object of mass m which is moving arbitrarily, having the position and velocity , with and , a Liénard-Wiechert potentials are gotten, and the resulting gravitational fields are of the form
where is the vector going from the object position x’, to the observer position x, is the unitary vector, and the gravitational fields and are
being the normalized velocity of the object, and is the retarded time. The gravitational power emitted by the accelerated object per solid angle is
and using (81), (82), and (83), one gets
Thus, any object of mass m which is accelerated will emit gravitational radiation. Of course, the object must have a huge mass in order for this radiation to be observed. This result contrast a lot with the given by General Relativity since in this theory the object must have a quadrupole configuration in order to emits gravitational energy    , and the gravitational waves emitted are of quadrupole type waves.
If one has a charged particle of charge q and mass m which is accelerated, this particle will emits electromagnetic  and gravitational energy such that the ratio of electromagnetic to gravitational energy emitted is
where is the constant dielectric of the vacuum. For the electron, , this ratio is of the order of 1039, whatever be its acceleration, that is, the electromagnetic radiation totally dominate the gravitational radiation. In order for a object of mass m and charge q to emits the same gravitational energy as the electromagnetic energy emission, it would require that its mass would be . For and object having an electron charge, this mass would be of the order of 21 kg.
8. Radiation Reaction Force for Gravitational Field
In this section, one will follow the approach given in reference   for the radiation reaction force in terms of the external force. An accelerated object of mass m dissipate energy due to gravitational waves emission that the object does during its accelerated motion. From the expressions (69) and (70), the modified equation of motion due to to the transformation that leave invariant of form the Maxwell-like equations is
where F is the external force, and . Although the function is really unknown, it is possible to assume the differentiation and to write this expression as the Newtonian equation of motion
where is some matrix function depending on . Using this expression in (79), it follows that
where and are the direction and the angle resulting from the expression , and are the direction and the angle resulting from the expression , and q is the angle between and . Making the integration of this expression with respect the solid angle and over the time intervale within the external force is acting, , one gets the energy dissipated by the accelerated object as
Now, one assumes that this energy is due to the work done by a nonconservative force which moves the object from the position at to the position x at a time . Because , one obtains
Equaling (84) and (85), knowing that the angle between and v must be p (dissipative direction), and having with , one gets finally and expression for the radiation reaction force of gravitational field as
where v is the magnitude of the velocity of the object. Therefore, the dynamical equations of motion which take into account the energy lost by gravitational radiation due to the accelerated object is
with given by (86). Note that whenever , one has that too. If the external force is zero, there is not acceleration, implying that there is not radiation of gravitational waves, one expects that this happen experimentally.
Although there have been several indirect facts    that indicates that the gravitational waves speed could be the same speed of light
(  ), it is necessary a direct measurement of this important parameter for gravitational waves. The gravitational field W appears naturally from mathematical analysis, but requires experimental verification together with the law of force that it may define, . It looks like that the transformations between the inertial reference systems S
and S’ must be of the type Poincaré-Lorentz transformation, but one might bring the objection that the mass is not a Lorentz invariant as the charge is. However, this does not represent an obstacle to propose a relation of the type (40) since one must have a linear relation between these vector fields defined in S and S’. If Poincaré-Lorentz transformation is assumed valid here, then, one could use the function as the same function used on special relativity theory
9. Conclusion and Comments
It has been shown that Maxwell-like equations for gravitational field can appear without having any relation at all with General Relativity. In fact, these types of equation could have been appeared much earlier than the formulation of the electrodynamics. Thus, it is absolutely astonishing the existing closely relation between post Newtonian gravity theory and electrodynamics theory. Several implications of Maxwell-like gravitational field formulation has been shown based on what we already know from Electrodynamics formulation. There are several results that one needs to point out: First, the energy associated to gravitational field or gravitational waves is a well defined concept. Second, linear gravitational waves appear from this extended Newtonian gravity theory, and General Relativity is not needed at all to explain their existence and their properties. Third, gravitational waves radiation appears for any accelerated object of mass m, and dipole type of radiation is always expected. Using this extended Newtonian approach for gravity, it is possible to include in the body equations of motion the dissipative effect caused by the emission of gravitational energy due to the acceleration of the body. This has been done through the gravitational radiation reaction force depending of the external force. Finally, due to this vector formulation of the linear gravity, one could say that to unify gravity with quantum mechanics, gravity could be considered as a vector field quantity with spin one boson as interacting quantum particle, instead of tensor quantity with spin two bosons as interacting quantum particle (General Relativity), and one could do Quantum Gravitodynamics (QGD) under this approach through (64), although a concern could arise for not having a small parameter on this expression. With respect to the resent claimed gravitational waves detections    , maybe it is possible to see a component of the Maxwell-like gravitational waves here, or these same waves be of Maxwell-type gravitational waves. However, these statements would require a deep careful analysis which could be made later on.
 Vollick, D.N. (2004) Palatini Approach to Born-Infeld-Einstein Theory and a Geometric Description of Electrodynamics. Physical Review D, 69, Article ID: 064030.
 Kopeikin, S. and Bashhoon, B. (2002) Gravitomagnetic Effects in the Propagation of Electromagnetic Waves in Variable Gravitational Fields of Arbitrary-Moving and Spinning Bodies. Physical Review D, 65, Article ID: 064025.
 Iorio, L. and Luchesi, D.M. (2003) LAGEOS-Type Satellites in Critical Supplementary Orbit Configuration and the Lense-Thirring Effect Detection. Classical and Quantum Gravity, 20, 2477.
 Thirring. H. (1918) über die formale Analogie zwischen den elektromagnetischen Grundgleichungen und den Einsteinschen Gravitationsgleichungen erster Naherung. Physikalische Zeitschrift, 19, 204-205.
 Lense. J. and Thirring. H. (1918) Uber den EinfluB der Eigenrotation der Zentralkorper auf die Bewegung der Planeten und Monde nach der Einsteinschen Gravitationstheorie. Physikalische Zeitschrift, 19, 156-163.
 Taylor. J.H. and Weisberg. J.M. (1982) A New Test of General Relativity—Gravitational Radiation and the Binary Pulsar PSR 1913+16. Astrophysical Journal, 253, 908-920.
 Li, L.-S. (2014) Gravitational Radiation Damping and Evolution of the Orbit of Compact Binary Stars (Solution by the Second Perturbation Method). Journal of Astrophysics and Astronomy, 35, 189-200.