Besides translation, spin around its axis and rotation around center of the Milky Way, the Sun performs relative motion in the solar system Laplacian plane, also. This motion was anticipated by Newton himself, in his Principia.
The form of the Sun’s orbit is substantially different from the other solar system bodies’ orbits. Namely, the Sun moves along the path composed of the chain of large and small loops [1, 2, 6, 9]. This chain is situated within the circular outline with the diameter approximately twice as large as the Sun’s is. Under supposition that the solar system is stable, the Sun is going to move along it, in the same region, for eternity, never reitereiting the same path.
It was also shown in this work that velocity and acceleration of the Sun’s center of mass are completely defined by the relative velocities and accelerations of the planets with respect to the Sun.
As said, Isaac Newton was the first who pointed to the fact that the Sun moves around the center of mass of the solar system (Philosophiae Naturalis Principia Mathematica): ” …since that centre of gravity (the solar system mass centre) is constantly at rest, the sun, according to the various positions of the planets, must continually move every way, but will never recede far from that centre”.
If the solar system is treated as a stable, isolated system of the point mass particles moving under mutual gravitational interactions, two dynamic conservation principles may be used for the study of its motion: conservation of the momentum and of the angular momentum of the system.
The consequence of the first rule is uniform motion of the system’s mass center C, while the consequence of the second is motion of the system in one, Laplacian, or invariant plane .
This plane is within 0, 5 of the Jupiter’s orbital plane and may be regarded as the weighted average of all planetary orbital planes. The point mass particles solar system model, involves necessity of neglecting differences between the orbital planes mainly originated in the transfer of the (small, but changeable) Sun’s and planet’s spin angular momenta to its total angular momentum.
Neglecting rotation of the solar system, as a whole, around the center of our galaxy, this plane moves translatory, together with the mass center C through the space.
Existence of the invariant plane permits introduction of an “inertial” reference frame xCy lying in it. Another Cartesian coordinate system x’0y’ (0x’ and 0y’ parallel with Cx and Cy) was adopted as the relative, that is, the heliocentric frame of reference (Fig. 1).
Fig. 1 – Inertial and Relative Frame of Reference
Position, Velocity and Acceleration Vectors of the Sun
It seems that, using astronometrical data positioning the outer planets from 1653. to 2060., P. D. Jose was the first one to determine the Sun’s orbit, in his work , 1965. This paper was abundantly cited later in the works treating the sunspot activities.
In his work  R. Bitsch has determined the Sun’s orbit integrating differential equations of motions of this celestial body exposed to the resultant of the planets’ attracting gravitational forces. It was supposed that these heavenly bodies are in circular, uniform motions around the Sun. The initial planet’s configuration was adopted arbitrary, because the author assumed that, in the long term, it does not affect the shape of the orbit.
One would be tempted to use term Kepler’s for the model employed in the work  and Copernicus’ for the model in , but it would not be correct, since the Sun moves in both models, of course.
In fact, the choice of either model in determining the trajectory, velocity or acceleration of the Sun is completely irrelevant, since the orbit, velocity and acceleration of the Sun’s center of mass are very small compared with the correspondent kinematic parameters characterizing motions of the planets. Correctly defining the initial conditions is all that matters.
For that reason “Copernican” model is the adopted here: the planets move at an average distance, with average angular velocity around the Sun.
Taking the planets’ configuration on 21 March 1978  as initial conditions and the model in which all planets, excepting Mercury (its orbit is not stabilized in the solar system’s invariant plane yet) move uniformly around the Sun, M. Marjanov  obtained, practically, the same form of the Sun’s path as Jose had. This trajectory covers the time interval of fifty years: 21. 03. 1978. + 50.
Since, in the inertial plane of reference must be
it follows out that
Thus, position, velocity and acceleration of the Sun’s center of mass are completely defined by the relative positions, velocities and accelerations of the planets with respect to the Sun.
Fig.2 shows the orbit of the Sun from 21. 03.2000. to 21. 03.2040. It, as mentioned, corresponds to the orbit that got Jose /1/, but is substantially different from the one given by Bish /2/, as the latter had wrongly assumed that initial conditions do not affect the form of a path.
The contour of the Sun is given as referential and the dot marks denote years.
Fig. 2 – Sun’s Path from 21.03.2000. – 21.03.2040.
Concerning the influence of the inner planets on the form of the path, Jose was right: it is quite negligible. Their influence would become visible only if a part of seemingly smooth path was magnified thousand to ten thousand times. The influence of the inner planets actions is far more evident when speed and acceleration of the Sun are considered.
The maximal distance from Sun to the center of mass of C ~ 1,5 10 km is obtainable when all planets are lined up on the same side of the star and the minimal ~21000 km, when Jupiter is on one side and all the other planets on the other side of the Sun, in the same direction.
Of course, chances for exactly such alignments of the celestial bodies are reduced to zero /11/.
The anticipated path of the Sun over the next 2000 years is represented in the Figure 3. It is situated within the circular outline with the diameter at least twice size of the Sun’s. Again, the Sun’s disk is given for comparison.
Provided that the solar system is stable, Sun is going to move in this region for eternity, never reitereiting the same path.
Fig. 3 – Sun’s Orbit Over 2000 Years
Sun’s velocity from 21.03.2000. to 21.03.2040. is represented in the Figure 4. Here we can see again that the outer, giant planets, especially Jupiter and Saturn have dominant influence on this kinematical quantity. When these two planets are in conjunction (~ every 20 years), the speed of the Sun is maximal ~ 16 (15,7), and when in opposition, the velocity is minimal, about 9 (8,8).
Fig. 4 – Sun’s Velocity from 21.03.2000. to 21.03.2040.
The acceleration diagram from 21.03.2000. to 21.03.2040. is given in the Figure 5. It shows much more irregularities and roughness than the previous diagram and represents, of course, a measure of the resulting planets’ gravitational attractions. No need to say that the peaks in this diagram correspond to the different conjunctional combinations of two or more planets.
Fig. 5 – Sun’s Acceleration from 21.03.2000. to 21.03.2040.
If all the planets fell in the same direction, on the same side of the Sun, therefore, when they were all in conjunction, the Sun would have a maximal acceleration ~0,27 /.
When Jupiter is on one side and all the other planets on the other side of the Sun, in the same direction, acceleration of the Sun is minimal ~0,14 /.
Solar System Bodies’ Gravitational Interactions
Previous diagram will be more understandable if one looks at the following Table of the of the mean Sun – planets and the maximal planet- planets interactions.
Table: Average Sun – Planets & Maximal Planet- Planets Interactions
The exposed table contains two kinds of data. The first row and the first column represent , = 1,2, … ,8 ,that is, the average Sun – planets interactions. These forces are average because of the adopted “Copernicus’” and not of the “Kepler’s” model.
Of course, the Sun-Jupiter interaction is the greatest one and it was taken to be the referential: its value is 100. An interesting fact is that the second one is the Sun – Venus interaction and that the Sun-Earth is of the same order of magnitude as the Sun-Saturn interaction, although the Saturn’s mass is 95 times greater than the Earth’s.
All the other values in the Table 1 represent , , = 1,2, … ,8 : maximal, that is, conjunctional interactions of the planets. In the range of the here adopted five digit precision, the only existing is the one between Jupiter and Saturn.
It is evident that motions of the inner planets, especially that of Venus and Earth are the main cause of the fluctuations around a smooth acceleration curve representing the influence of the outer planets, and above all, the influence of Jupiter Saturn.
The Sun and the planets in the paper were simulated by the model of stable, isolated system of the point mass particles, moving under mutual gravitational interactions in the averaged Laplacian plane.
Shape of the solar path is entirely different from the orbits of the other bodies moving around it. Namely, the Sun moves along the path in the form of chain consisting of large and small loops. This chain is situated within the circular outline with the diameter approximately twice as large as the Sun’s is. Under supposition that the solar system is stable, the Sun is going to move along it, in the same region, for eternity, never reitereiting the same path.
Then the velocity and acceleration diagrams of center of mass of the Sun were given.
At last, for better understanding of the obtained results, Table of the of the mean Sun – planets and the maximal planet- planets interactions was contributed.
 P. D. Jose: Sun’s Motion and Sunspots, The Astronomical Journal, Vol. 70, Nr. 3, 1965.
 R. Bitsch: The Sun Moves! email@example.com, 2003.
 M. Milanković: Osnovi nebeske mehanike, Naučna knjiga, Beograd, 1988.
 T. Anđelić: Uvod u astrodinamiku, Matematički institut, Beograd, 1983.
 1978 Yearbook of Astronomy, Ed. P. Moore, W.W.Northon&Company Inc, New York.
 K. Croswell: Planet Quest, Oxford University Press, Oxford, 1999.
 Heliocentric Ephemeris – 2000 – Achenar – UK Mirror.
 Marjanov M.: Gravitational Resonances, Bull. Astron., Belgrade, 156, 9 – 19, 1997.
 Marjanov M.: The Sun’s Orbit and Centers of Gravitation of the Solar System, Lecture on the Mathematical Institute SANU, Department of Mechanics, March 2005.
 Marjanov M.: On the Cause of Resonant Motions of Celestial Bodies, Facta Universitatis, Series: Mechanics, Automatic Control and Robotics, vol. 4, No 17, 2005.
 Marjanov M.: Pythagoras Numbers and Quasiperiodicity of the Solar System, Zbornik radova Gradjevinskog fakulteta u Subotici, 2006.
 Marjanov M.: One Solution of the Newton’s Three Body Problem, Lecture on the Mathematical Institute SANU, Department of Mechanics, Dec. 2008.