Notes on the Motion of Celestial Bodies

A novel method for the computation of the motion of multi-body systems is proposed against the traditional one, based on the dynamic exchange of attraction forces or using complex field equations, that hardly face two-body problems. The Newton gravitational model is interpreted as the emission of neutrino/gravitons from celestial bodies that combine to yield a cumulative flux that interacts with single bodies through a momentum balance. The neutrino was first found by Fermi to justify the energy conservation in β decay and, using his model; we found that the emission of neutrino from matter is almost constant independently from the nuclides involved. This flux can be correlated to Gauss constant G, allowing the rebuilding of Newton law on the basis of nuclear data, the neutrino weight and the speed of light. Similarly to nature, we can therefore separate in the calculations the neutrino flux, that represents the gravitational field, is dependent on masses and is not bound to the number of bodies involved, from the motion of each body that, given the field, is independent of the mass of bodies themselves. The conflict between exchanges of forces is avoided, the mathematics is simplified, the computational time is reduced to seconds and the stability of result is guaranteed. The example of computation of the solar system including the Sun and eight planets over a period of one to one hundred years is reported, together with the evolution of the shape of the orbits.


Introduction
The problem computing the interactive forces of a group of celestial objects and consequently of predicting their orbital motion for all future times has been, after Newton, the headache of scientists and mathematicians till nowadays. restricted three-body problem in which the third mass is negligible.
The general n-body problem fell in the hand of mathematicians: the elegant Hamilton's approach transforms the original Newtonian system of 3 n second order differential equations of motion in a system of 6n first order differential equations with 3n position coordinates and 3n momentum values; however no analytical or numerical solution has been found.
After the extensive work of Poincare [1] [2], many proposals and models of Celestial Mechanics, both numerical and analytical or a mixture of the two, are present in the literature, and also full books have been written.
Unfortunately all works are based on the idea of interacting forces, that, after the failure of eminent mathematicians and on the basis of our modest experience, we consider an incorrect representation of gravity.
We did not find in the literature an alternative view as that presented in this paper and therefore, instead of boring the reader with incorrect mathematics and unsuccessful results, we want only to cite the recent opinion of Celletti [3], an expert researcher working in the field, that sounds like a voice outside the crowd.
"An accurate prediction of the dynamics of the objects of the solar system often requires very long computations; even the easiest problem provided by the two-body model deserves computational skill in solving Kepler's equation, which allows deriving the Keplerian elements of the orbit as a function of time.
In passing from the two to the three-body problem, an increasing complexity, due to nonintegrability and chaos, is a trademark of Celestial Mechanics; the three body problem motivated the development of perturbation theories aimed to find approximate solutions of the equations of motion." Therefore the fault is not of the Newton model, nor of mathematics, but of the way mathematic is applied.
For those that are interested in the relation between mathematics and physics we suggest the last recent book of Penrose [4] The Road to Reality. The book is just over 1100 pages, of which the first 383 are dedicated to mathematics: Penrose's goal is to acquaint inquisitive readers with the mathematical tools needed to understand nature.
Only on the last page of the conclusion, he wisely invites the reader to open a different window for understanding the unknowns of the Universe. Journal of Applied Mathematics and Physics In the case under study, of celestial mechanics, the escape from the exchange of forces model was found in the so called field equations models.
Every engineer expert in fluid dynamics had the opportunity to deal with the Navier-Stokes field equations, that describe the motion of viscous fluids through mass, momentum and energy balance, and is aware of the problems encountered in their use to solve real problems.
Remember that these equations were written applying balances to a cell of fluid and then transforming finite differences to derivates by the limit to zero procedure.
The analytical solution of these equations is, except for simpler cases, impossible, due to the large nonlinearities, the nature and properties of fluids (forgetting the turbulence tensor) and the boundary conditions that are the cause of instabilities.
The numerical solution involves the opposite procedure of transforming differential equations into finite difference equations, that are usually different from the original elementary cell balances, are numerically less stable and need a convergence with boundary conditions. The use of parametric functions or series of functions followed by parameters optimization, partially solves the instability problems but is quite approximate. Therefore, in solving scientific or engineering problems, we decided to maintain the original cell balances, selecting the cell size, the cell equations according to the nature of phenomena involved and including boundary conditions: the resulting system of nonlinear equations can be easily solved by available computer programs [5], almost independently from the problem complexity and size.
The use of computers, not available at the times of Newton and Poncarè, helps in the simplification of mathematical algorithms and computation procedures and, contrarily to the common opinion, to bring researchers nearer to physical phenomena, if supported by experimental data.
It would have been impossible for Maxwell to write his electromagnetic field equations, without the anticipation, the experiments and the investigations of Faraday, a genius with poor format education and limited knowledge of mathematics.
Maxwell's equations are again partial differential equations that relate the electric and magnetic fields to each other and, in spite of the elegant and exhaustive description of the phenomena involved, they are very difficult to solve due to nonlinearities and boundary conditions. As far as multiple attractions or repulsion are concerned, the two bodies problem is again a limit for analytical calculations and numerical calculations do no go far ahead.
When Einstein was in trouble with the problem gravitation, the example of Maxwell, whose equations are in line with Lorenz transformation and hence with special relativity, paved the way. Journal of Applied Mathematics and Physics Einstein field equation theory, as explained in [6], are hard to swallow but, after the theoretical presentation of tensor properties, the elegance of mathematics is impressive and the capability of including both fluid dynamics (Navier-Stokes) and electrodynamics (Maxwell) is surprising.
For this reason Einstein field equation was considered, during the last century, a mathematical monument comprehensive of all the laws of the Universe.
However we feel uncomfortable with its use for solving real problems, even with the aid of supercomputers and there is no knowledge of its use for systems exceeding two bodies.
We have also to recognize that the theory is far to be simple, in contradiction with Einstein main thought: "most of the fundamental ideas of science are essentially simple and may, as a rule, be expressed in a language comprehensible to everyone".
The concept of exchange of forces has been overcome by the idea of the curvature of space-time, a sort of mobile mathematical basin built around the bodies that constrains all physical bodies and make them move along specific lines.
In the absence of the knowledge of the nature of gravity we hold this as a sophisticated representation of a phenomenon that Galileo tried to explain to Simplicio making him rotate a physical bucket full of water, in the Dialogue Concerning the Chief World Systems (Dialogo dei Massimi Sistemi) [7].
Einstein thought dominated the scene of the last century generating a school of scientists that under his name are proposing new investigations and advance new theories.
My feeling is that we are shortening times and enlarging figures beyond what Nature is telling us.
What follows is an example of joining experimental facts taken from mechanics, nuclear data and quantum mechanics theories, to change our view of the motion of bodies in the sky.

The Gravity Field Equation Simple
The number of important laws of physics can be counted on the fingers of one hand and may change with time; the experiments, if properly performed, will step in science forever, ready to be explained and used in different periods of time, One of these important experiments is that performed by Galileo with the falling of objects of different shapes, weight and size at the same time, due to Earth gravity.
Galileo could not appreciate the secret hidden under his experiment and he spent his life, even when arrested in home, to look at the sky with his telescope; the sky was also the preferred reference of Copernicus, Kepler and Newton for gravitational phenomena.
Till nowadays gravitation has been segregated in the sky, but the Galileo experiment suggests another direction that even Einstein could not imagine: if the action of gravity is the same on a ball of iron and on a piece of wood it should Journal of Applied Mathematics and Physics not depend on the chemical and type of nuclide involved but has to relight to the elementary particles, neutrons and protons that are common to all matter. Therefore, if we want to understand the nature of gravity, our investigation should be moved away from the fantastic objects appearing in the sky and fix our attention on the nuclear structure of matter.
We had to wait four centuries, after Galileo, before Fermi [8] developed a theory of β decay and justified the existence of the neutrino ν, an elusive particle that can cross matter without interaction.
He assumed the following nuclear dynamic transformations of protons p, neutron n and electrons/positrons β − /β + with the neutrino ν compensating for the lack of energy in the electrons emission spectrum: This scheme also justifies the number of nuclear bonds and the nuclear binding energy of nuclides [9] and, if we argue that this reaction scheme is valid not only for radioactive nuclides but can be extended to all nuclides at all temperatures, we have a model and a new law for the emission of radiation from the bulk of matter.
From Equation (1) we can compute the neutron and proton dynamics in all nuclides [10]: We assume k 1 = 0.0009625, experimentally known, being related to the 12 minutes half life of neutron, and derive the other constants by fitting the neutron/ proton distribution of the 1806 known nuclides including stable nuclei and γ emitters.
The least square fit of the proton and neutron distribution gives the following values with an astonishing determination index equal to 99.7%: If we use these constants to compute the rate of neutrino emission, we discover that it is almost constant for all nuclides with a mean value of F o = 6.668E+23 neutrino per kg and per second and this value does not significantly change from light to heavy nuclides as first reported in [10]. We remember that, following Plank, matter emits radiation from the surface independently from its chemical or physical composition (another experimental fact), in agreement with our result based on fitting data of all known nuclides.
The radiation is generated and transferred at various temperatures but we know that the nuclear clock is not affected by temperature and therefore the constancy of atomic matter emissivity is an important new law of physics.
The radiation is generated by nuclei in the bulk of a body, only part of it reaches a small layer near the surface where is emitted in the form of light, with Journal of Applied Mathematics and Physics a Black Body spectrum, the remaining part, the most important one is filtered and degraded by matter and escapes undetected from the body in the form of the elusive neutrino and is the one justifying Galileo experiment and the one responsible for "moving the Sun and the other Stars".
However only a small portion of the neutrino flux interacts with the physical bodies, the majority being dispersed in the space: the neutrino is therefore responsible for the increase of the entropy of the Universe, with the stars loosing mass, the attraction between celestial bodies decreasing and the expansion of the Universe continuing slow in our time frame.
Anyhow, once known the position and the mass M i of celestial bodies, the vector flux of neutrino F s per unit area in a point of the space at a distance R i and unit vector r i can be easily computed by vector summation of bodies flux contributions: The neutrino flux is spread on the surface of the sphere of radius R i but we are interested in the surface available that is the cross section σ n of nucleon (proton-neutron) having mass m n (kg) and radius r n (m): ( ) that is equivalent to Newton attraction force.
The Gauss constant G therefore can be easily related to nuclear parameters:  Here we use, as example, a simplified solar system with the Sun and only the eight planets, with the data, taken from literature, reported in Table 1.
To make it simpler we start with the origin in the Sun and with the planets positioned along the x axis at the given distance from the Sun.       Everybody can appreciate how easily and comfortably this calculation can be done without the anxiety of mutual interactive competing forces.
The neutrino flux is therefore the real field equation for proceeding to celestial mechanics calculations: the cause of motion is a physical reality, separated from the effect, the motion itself.
Therefore Newton law should not be interpreted as an exchange of forces but as the emission by bodies of neutrino in the space that causes a neutrino unbalance on the receiving bodies.
Newton law is still valid, the mathematics is formally the same but the use of equations is different: cause and effect are separated.

Fast and Stable Planetary Orbit Calculation
Now that we have a precise and simple representation of the gravitational field, an object, immersed in this field, is willing to move and we know that our planets move following a quasi-circular orbit or better an elliptic orbit around the Journal of Applied Mathematics and Physics

Sun.
Due to planets motion, the Sun also moves together with them around a barycenter whose coordinates are: The velocity and acceleration of all orbiting bodies are better defined using polar coordinates (r, ϑ ) where r and ϑ are the radial and tangential versors. The coefficient of radial acceleration, following Newton, is equal to F g and that of the tangential acceleration is zero and is related to the momentum conservation: For a circular orbit de second derivative of r, r′′ is equal zero, Equation (14) can be used for a first local approximation but if we insist for a Keplerian orbit we can use the Vis-viva equation [11], derived from the conservation of energy and momentum, that we can rewrite in the following form: from Equations (13) and (15) we can compute the second derivative of r, that is related to a property of the ellipse, the semi-major axis a: The gravitational field F g can be rigorously computed, once known the position of masses M i while Equations (14) and (15) contain hypothesis about the orbit, are useful locally for the calculation, but the final orbit shape is mainly influenced by the change of F g from time to time.
The only additional variable we need is the radial velocity v r in Equation (12), that is zero for circular orbits and is proportional to v ϑ for the elliptical ones, depending on the eccentricity e, that, like the semimajor axis, should be a result of the calculation.
Therefore we assume and find the value of η for the planets imposing the closure of the orbit after the first annual cycles.
The last line of Table 1 shows the values of η used for the calculations. The vector v r is directed along the radius r toward the barycenter b and using Journal of Applied Mathematics and Physics The vector v ϑ is perpendicular to F g and thus using Equations (14) or (15) and the velocity of each celestial bodies is simply defined by: The problem of the orbit calculation is now reduced to elementary algebra through the following simple steps: 1) Given planets positions, compute the barycenter using (10).
2) Compute the fluxes F s (F sx , F sy ) and F g (F gx , F gy ) from (9)  Some trial runs can be used to adjust η in Equation (17). Note that steps 1 and 2 are rigorous and there is no conflict between bodies.
Step 3 and 4 are independent of the mass of the bodies that has been accounted in step 1 and 2 and therefore no conflict is present here too.
Consequently the computation is very accurate and stable with a computer time that can be measured in seconds for processing the elementary algebra (21), most of the time being dedicated to printing results. The calculation has been made using locally the circular orbit approximation that is maintaining the second radial derivative equal zero and computing v ϑ with Equation (14), but introducing the correction factors η of v r of Table 1.
The orbits are very clean: almost one revolution and one line for the majority of planets and shape apparently indistinguishable from a circle.
If we introduce the second derivative in Equation (15) and maintain the same correction factor for v r the results appear very similar, the controlling factor being the field equation F s (4) representing the flux or momentum F g (7) of neutrino.
We report for comparison in Figure 5 the terrestrial planets orbits using Equation (15) for the computation of v ϑ .  The orbits of Figure 5 are very similar to those of Figure 4(c) as if the gravitational field is prevailing over the tricks of the calculation.
Changing the integration step from one to two days has no effect and the computation time is almost instantaneous.
The computation over an extended time of one hundred Earth years is equivalently stable allowing many rounds of the terrestrial planets and some complete run of the gaseous ones; the computational time again is less than ten seconds.
The calculation has been made with a time step of 1.2 days exceedingly small for the phenomena involved and in agreement with the graphical ability of Excel. for the terrestrial and gaseous planets using Equation (14) for v ϑ and the η correction for v r of Table 1.
In Figure 6(a) one can distinguish the effect, on the motion of barycenter, of the rapid rotation of the terrestrial planets and the slow motion of the gaseous ones, In Figure 6(b) one can appreciate the slow motion of the pachyderm Sun to compensate the planets merry go-round.
In Figure 6(c) the orbits of the terrestrial planets appear to shift revolution by revolution, as if each planet is trying to find the right position compatible with the position of the other planets.
We have selected some orbits within the bunch and they appear identical at the beginning, in the center and at the end of calculations and they appear identical but shifted in space (not reported here).
The orbits of the gaseous planets appear regular, even if they had not enough time to stabilize themselves and to balance the entire solar system. There is no way, even superimposing figures, to distinguish about the results and again it is clear that the dominating force is the neutrino flux F g , in spite of the tricks used to force the Keplerian elliptic orbits.
To better understand the evolution of the orbits we computed, as probably Kepler did, the distance of the planets from the Sun in function of time.
We discover that the time necessary for a planet to settle in stable orbits and find its place in the sky, depends on their frequency and is less for Mercury than for Mars, and that the orbits, appearing circles, are really elliptic. Figure 8 shows that Mecury rapidly reaches a steady state distance and Figure   9 gives a gross variation of the distance in a large period of time.
Going to details Figure 10 shows the variation of distance from the Sun within each Mercury year in a sinusoidal behavior typical of elliptic orbits.         additional time to reach a stable position and still shows an elliptic shape for its orbit.
For the gaseous planets the time needed for stabilization will be longer, but this should not be a problem for our calculation method.
The method presented virtually has no limits in stability and computational time and can be used for studying different starting points, adding additional celestial bodies within or outside the solar system and simulating the evolution over millennia of known clusters of rotating bodies.
We have given only a flavor of what can be done with this tool for analyzing the celestial motion in detail and of the quantity of data available under the simple graphical representation.
We have also shown the merits of this new approach to gravity for multiple celestial bodies calculations.
The comparison with experimental data will be of help both for refining the method and for acquiring a quantitative picture of small and large pieces of the Universe.

Conclusions
Gravity has always been seen as an exchange of forces, between massive bodies, that decreases with the square of the distance coherently with Newton Universal Gravitational Law.
The transmission of this force is considered by Newton instantaneous but, due to the large astronomical distance and assuming that the action of this forces propagates with the velocity of light, the interval between the start and the arriv- As a general rule, derived by the experience on many other applications, we can say that, if a model or a set of mathematical expressions do not fit with physical phenomena, either mathematics are wrong or incorrectly applied.
Therefore we feel comfortable with Newton Law, living on our Earth, and should be also comfortable, travelling in the near and far space, provided we understand and use gravity the right way.

Conflicts of Interest
The author declares no conflicts of interest regarding the publication of this paper.