zotero/storage/6IJRQCJB/.zotero-ft-cache

829 lines
21 KiB
Plaintext
Raw Permalink Blame History

This file contains ambiguous Unicode characters

This file contains Unicode characters that might be confused with other characters. If you think that this is intentional, you can safely ignore this warning. Use the Escape button to reveal them.

arXiv:1301.6045v2 [physics.class-ph] 6 Feb 2013
Newton-Machian analysis of Neo-tychonian model of planetary motions
Luka Popov
University of Zagreb, Department of Physics, Bijeniˇcka cesta 32, Zagreb, Croatia E-mail: lpopov@phy.hr
Abstract. The calculation of the trajectories in the Sun-Earth-Mars system will be performed in two different models, both in the framework of Newtonian mechanics. First model is well-known Copernican system, which assumes the Sun is at rest and all the planets orbit around it. Second one is less-known model developed by Tycho Brahe (1546-1601), according to which the Earth stands still, the Sun orbits around the Earth, and other planets orbit around the Sun. The term “Neo-tychonian system” refers to the assumption that orbits of distant masses around the Earth are synchronized with the Suns orbit. It is the aim of this paper to show the kinematical and dynamical equivalence of these systems, under the assumption of Machs principle.
PACS numbers: 45.50.Pk, 96.15.De, 04.20.-q, 45.20.D-, 01.65.+g
Submitted to: Eur. J. Phys.
Accepted for publication: 16 Jan 2013
Newton-Machian analysis of Neo-tychonian model of planetary motions
2
1. Introduction
The discussion of motion of celestial bodies is one of the most interesting episodes in the history of science. There are two diametrically opposite schools of thought: one that assumes that the Sun stands still, and Earth and other planets orbit around it; and another that assumes that the Earth stands still, and Sun and other planets in some manner orbit around the Earth. The first school of thought comes from Aristarchus (310230 bc) and is generally addressed as heliocentrism, another from Ptolemy (90-168 bc) and is generally known as geocentrism. Since Aristotle, the ultimate authority in science for more than two millennia, accepted the geocentric assumption, it became dominant viewpoint among scientists of the time. The turnover came with Copernicus (so-called “Copernican revolution”) who in his work De Revolutionibus proposed a hypothesis that the Sun stands in the middle of the known Universe, and that Earth orbits around it, together with other planets. Copernicus system was merely better than Ptolemys, because Copernicus assumed the trajectories of the planets are perfect circles, and required the same number of epicycles (sometimes even more) as Ptolemys model [1]. The accuracy of Ptolomys model is still a subject of vivid debates among historians of science [2].
The next episode in this controversy is Keplers system with elliptical orbits of planets around the Sun. That system did not require epicycles, it was precise and elegant. It is therefore general view that Keplers work finally settled the question whether it is the Sun or the Earth that moves. But what is less known is that Tycho Brahe, Keplers tutor, developed a geostatic system that was just as accurate and elegant as Keplers: the Sun orbits around the Earth, and all the other planets orbit around the Sun. The trajectories are ellipses, and all the Keplers laws are satisfied. In that moment of history, the Keplers and Brahes models were completely equivalent and equally elegant, since neither of them could explain the mechanism and reason why the orbits are the way they are. It had to wait for Newton.
Sir Isaac Newton, as it is generally considered, gave ultimate explanation of planetary motions that was in accord with Keplers model, and excluded Brahes one. The laws of motions and the inverse square law of gravity could reproduce all the observed data only with the assumption that the Sun (i.e. the center of mass of the system, which can be very well approximated by the center of the Sun) stands still, and all planets move around it. According to Newtons laws, it is impossible for small Earth to keep the big Sun in its orbit: the gravitational pull is just too weak. This argument is very strong, and it seemed to settle the question for good.
But in the end of 19th century, the famous physicist and philosopher Ernst Mach (1839-1916) came with the principle which states the equivalence of non-inertial frames. Using the famous “Newtons bucket” argument, Mach argues that all so-called pseudoforces (forces which result from accelerated motion of the reference frame) are in fact real forces originating form the accelerated motion of distant masses in the Universe, as observed by the observer in the non-inertial frame. Some go even further, stating that
Newton-Machian analysis of Neo-tychonian model of planetary motions
3
“every single physical property and behavioral aspect of isolated systems is determined by the whole Universe” [3]. According to Machs principle, the Earth could be considered as the “pivot point” of the Universe: the fact that the Universe is orbiting around the Earth will create the exact same forces that we usually ascribe to the motion of the Earth.
Machs principle played a major role in the development of the Einsteins General Theory of Relavity [4], as well as other developments in gravitation theory, and has inspired some interesting experiments [5]. This principle still serves as a guideline for some physicists who attempt to reformulate (“Machianize”) Newtonian dynamics [6, 7], or try to construct new theories of mechanics [8]. Some arguments and critiques against Machs principle have also been raised [9]. Since the time of its original appearance [10, 11, 12], Machs principle has been reformulated in numbers of different ways [13, 14]. For the purpose of this paper, we will only focus on the one of the consequences of Machs principle: that the inertial forces can be seen as resulting from real interactions with distant matter in the Universe, as was for example shown by A. Zylbersztajn [15].
The only question remains: are these forces by themselves enough to explain all translational motions that we observe from Earth, and can they reproduce the Tycho Brahes model? The discussion in this paper will show that the answer to this question is positive. In order to demonstrate it, we will consider the Sun-Earth-Mars system.
The paper is organized as follows. In section 2 an overview of two-body problem in the central potential and of Keplers problem is given. In section 3 the calculations of Earths and Mars trajectories are performed in the heliocentric system, both analytically (by applying the results from previous section) and numerically. In section 4 the calculations of Suns and Mars trajectories are performed in geocentric system, due to the presence of pseudo-potential originating from the fact of accelerated motion of the Universe. Finally, the conclusion of the analysis is given.
2. Two-body problem in the central potential
2.1. General overview
We start with the overview of two body problem in Newtonian mechanics. Although
there are alternative and simpler ways to solve this problem [16, 17], we will follow the
usual textbook approach [18, 19]. The Lagrangian of the system reads
L
=
1 2
m1r˙ 21
+
1 2
m2r˙ 22
U (|r1
r2|) ,
(2.1)
where U is potential energy that depends only on the magnitude of the difference of
radii vectors (so-called central potential ). We can easily rewrite this equation in terms
of relative position vector r ≡ r1 r2, and let the origin be at the centre of mass, i.e.
m1r1 + m2r2 ≡ 0. Solution of these equations are
r1
=
m2 m1 + m2
r
,
r2
=
m1 m1 + m2
r
.
(2.2)
Newton-Machian analysis of Neo-tychonian model of planetary motions
4
The Lagrangian (2.1) so becomes
L
=
1 2
µr˙ 2
U(r) ,
where r ≡ |r| and µ is the reduced mass,
(2.3)
1 µ
=
1 m1
+
1 m2
(2.4)
In that manner, the two-body problem is reduced to one-body problem of particle with
coordinate r and mass µ in the potential U(r).
Using polar coordinates, the Lagrangian (2.3) can be written as:
L = 1 µ(r˙2 + r2φ˙2) U (r) 2
(2.5)
One can immediately notice that variable φ is cyclic (it does not appear in the
Lagrangian explicitly). Consequence of that fact is momentum conservation law, since (∂/∂t)(∂L/∂φ˙) = ∂L/∂φ = 0. Therefore,
∂L ∂φ˙
=
µr2φ˙
=
const.
(2.6)
is the integral of motion. In order to find a solution for the trajectory of a particle, it is not necessary to
explicitly write down the Euler-Lagrange equations. Instead, one can use the energy conservation law,
E
=
1 2
µ(r˙2
+
r2φ˙ 2)
+
U (r)
=
1 2
µr˙2
+
2 2µr2
+
U (r)
(2.7)
Straightforward integration of (2.7) gives equation for the trajectory,
φ(r) =
dr/r2 2µ [E U (r)] 2/r2
(2.8)
2.2. Keplers problem
Let us now consider the particle in the potential
U
(r)
=
k r
,
(2.9)
generally known as Keplers problem. Since our primary interest is in the planetary
motions under the influence of gravity, we will take k > 0. Integration of eq. (2.8) for
that potential gives:
p r
=
1
+
e
cos
φ
,
(2.10)
where 2p is called lactus rectum of the orbit, and e is eccentricity. These quantities are
given by
p
=
2 µk
,
e=
1
+
2E2 µk2
(2.11)
Newton-Machian analysis of Neo-tychonian model of planetary motions
5
Expression (2.10) is the equation of a conic section with one focus in the origin. For
E < 0 and e < 1 the orbit is an ellipse.
One can also determine minimal and maximal distances from the source of the
potential, called perihelion and aphelion, respectively:
rmin
=
1
p +
e
,
rmax
=
1
p
e
.
(2.12)
These parameters can be directly observed, and often are used to test a model or a
theory regarding planetary motions.
3. Earth and Mars in heliocentric perspective
According to Newtons law of gravity, the force between two massive objects reads:
F
=
Gm1m2 |r1 r2|3
(r1
r2)
.
(3.1)
Which leads to a potential (F = ∇U)
U
(|r1
r2|)
=
Gm1m2 |r1 r2|
.
(3.2)
This is obviously Keplers potential (2.9) with k = Gm1m2, where G is Newtons gravitational constant.
Since the Sun is more than 5 orders of magnitude more massive than Earth and Mars, we will in all future analysis use the approximation
µ ≈ mi ,
(3.3)
where mi is mass of the observed planet. For the same reason, gravitational interaction between Earth and Mars can be neglected, since it is negligible compared with the
interaction between Earth/Mars and the Sun.
Using these assumptions, we can write down corresponding Lagrangians,
LES
=
1 2
mE
r˙ 2E
S
+
GmE MS rES
,
LM S
=
1 2
mM
r˙ 2M
S
+
GmM MS rM S
,
(3.4)
where mE and mM are masses of Earth and Mars, respectively. Subscripts ES (MS) correspond to the motion of Earth (Mars) with respect to the Sun. These trajectories
can be calculated using the exact solution (2.10) with appropriate strength constants
k and initial conditions which determine E and . Another way is to solve the Euler-
Lagrange equations numerically, using astronomical parameters [20] (e.g. aphelion and
perihelion of Earth/Mars) to choose the inital conditions that fit the observed data.
The former has been done using Wolfram Mathematica package. The result is shown
on Fig. 1.
For latter comparison, one could write out the expressions for the e and p parameters
for the Earth. Putting the expressions for energy (2.7) and momentum (2.6) into
Newton-Machian analysis of Neo-tychonian model of planetary motions
6
2.0
1.5
1.0
0.5
0.0
0.5
1.0
1.5
2.0
2.0
1.5
1.0
0.5
0.0
0.5
1.0
1.5
2.0
Astron om ical u n its
Figure 1. Trajectories of Earth and Mars in heliocentric system over the period of 2 years. Blue and red lines represent Earths and Mars orbits, respectively (color online).
Equations (2.11) it is straightforward to obtain
p
=
φ˙ 2 r4 GMS
,
e=
1
2GMSφ˙2r3 r˙2φ˙2r4 G2MS2
φ˙ 4 r6
,
(3.5)
where φ˙, r˙ and r are angular velocity, radial velocity and distance respectively, taken in
the same moment of time (e.g. in t = 0).
Fig. 2 displays motion of the Mars as viewed from the Earth, gained by trivial
coordinate transformation
rME(t) = rES(t) + rMS(t) ,
(3.6)
where rES(t) and rMS(t) are solutions of Euler-Lagrange equations for the Lagrangians (3.4). Equation (3.6) is just the mathematical expression of the Tycho Brahes claim.
Retrograde motion of the Mars can be useful in the attempt to understand and determine
orbital parameters, as was qualitative and quantitative shown by B. Thompson [21].
The acceleration that Earth experiences due to the gravitational force of the Sun
is usually referred as centripetal acceleration and is given by
acp
=
Fcp mE
=
GMS rE2 S
ˆrES
,
(3.7)
Newton-Machian analysis of Neo-tychonian model of planetary motions
7
2
1
0
1
2
2
1
0
1
2
Astron om ical u n its
Figure 2. Trajectory of the Mars as seen from the Earth over the period of 7 years. Calculation of this trajectory is done numerically in the heliocentric system.
where ˆr is the unit vector in the direction of vector r, rES(t) is radius vector describing motion of Earth with respect to the Sun, and Fcp is centripetal force, i.e. the force that causes the motion.
4. Sun and Mars in geocentric perspective
4.1. The pseudo-potential
From the heliocentric perspective, the fact that the Earth moves around the Sun results
with centrifugal pseudo-force, observed only by the observer on the Earth. But if we
apply Machs principle to the geocentric viewpoint, one is obliged to speak about the real
forces resulting from the fact that the Universe as a whole moves around the observer
sitting on the stationary Earth. Although these forces will further be considered as the
real forces, we well keep the usual terminology and call them pseudo-forces, for the sake
of convenience. Our focus here will be on the annual orbits, not on diurnal rotation
which requires some additional physical assumptions [8] [22] that are beyond the scope
of this paper.
The Universe is regarded as an (N + 1)-particle system (N celestial bodies plus
planet Earth). From the point of a stationary Earth, one can write down the Lagrangian
that describes the motions of celestial bodies:
L
=
1 2
N
mir˙ 2i
i=1
1 2
N i=1
Gmimj rij
N
i=1
GmE mi ri
Ups ,
(4.1)
Newton-Machian analysis of Neo-tychonian model of planetary motions
8
where rij ≡ |ri rj|, Ups stands for pseudo-potential, satisfying Fps = ∇Ups. Fps is the pseudo-force given by
N
Fps = m acp,i ,
i=1
(4.2)
where acp,i is centripetal acceleration for given celestial body (with respect to the Earth) and m is a mass of the object that is subjected to this force. Its easy to notice that the
dominant contribution in these sums comes from the Sun. The close objects (planets,
moons, etc) are much less massive than the Sun, and massive object are much further
away. The same approximation is implicitly used in section 3.
In the Machian picture, the centripetal acceleration is a mere relative quantity,
describing the rate of change of relative velocity. Therefore, centripetal acceleration of
the Sun with respect to Earth is given by Equation (3.7), with rES = rSE. All that considered, Equation (4.2) becomes
Fps
=
GmMS rS2 E
ˆrSE
,
(4.3)
where rSE(t) describes the motion of the Sun around the Earth, and m is the mass of the body under consideration.
We can now finally write down the pseudo-potential which influences every body
observed by still observer on Earth:
Ups(r)
=
GmMS rS2 E
ˆrSE
·
r,
(4.4)
where r(t) describes motion of particle of mass m with respect to the Earth. Notice
that this is not a central potential.
4.2. Sun in Earths pseudo-potential
In order to determine Suns orbit in Earths pseudo-potential, one needs to take dominant contributions of the Lagrangian (4.1), as was explained earlier. Taking into account the expression for pseudo-potential given in Equation (4.4), one ends up with
LSE
=
1 2
MS
r˙ 2SE
GMS2 rSE
.
(4.5)
This Lagrangian has the exact same form as the reduced Lagrangian (2.3). That
means that we can immediately determine the orbit by means of Equations (2.11) by
substituting µ = MS and k = GMS2. This leads to the following result (subscript SE will be omitted):
p
=
φ˙ 2 r4 GMS
,
e=
1
2GMSφ˙2r3 r˙2φ˙2r4 G2MS2
φ˙ 4 r6
,
(4.6)
Newton-Machian analysis of Neo-tychonian model of planetary motions
9
which is the exact equivalent as the previous result given in Equations (3.5), since φ˙, r˙ and r are relative quantities, by definition equivalent in both models. We can therefore conclude that the Suns orbit in the Earths pseudo-potential is equivalent to that observed from the Earth in the heliocentric system.
It remains to show the same thing for Mars orbit.
4.3. Mars in Earths pseudo-potential
In the similar way as before, we take the dominant contributions of Lagrangian (4.1)
together with Equation (4.4) and form the following Lagrangian
LM E
=
1 2
mM
r˙ 2M
E
+
GmM MS |rME rSE|
GmM MS rS2 E
ˆrSE
·
rM E
,
(4.7)
where subscript ME refers to the motion of Mars with respect to the Earth, and rSE(t) is the solution of Euler-Lagrange equations for the Lagrangian (4.5).
2
1
0
1
2
2
1
0
1
2
Astron om ical u n its
Figure 3. Trajectories of the Sun (yellow) and the Mars (red) moving in Earths pseudo-potential over the period of 7 years (color online). Calculation of this trajectory is performed numerically in the geocentric system.
The Euler-Lagrange equations for rME(t) using Lagrangian (4.7) are too complicated to be solved analytically, but they can easily be solved numerically. The numerical solutions for the equations of motion for both the Sun and Mars are displayed in Fig. 3. The equivalence of trajectories gained in two different ways is obvious, justifying the model proposed by Tycho Brahe.
Newton-Machian analysis of Neo-tychonian model of planetary motions
10
5. Conclusion
The analysis of planetary motions has been performed in the Newtonian framework with the assumption of Machs principle. The kinematical equivalence of the Copernican (heliocentric) and the Neo-tychonian (geocentric) systems is shown to be a consequence of the presence of pseudo-potential (4.4) in the geocentric system, which, according to Mach, must be regarded as the real potential originating from the fact of the simultaneous acceleration of the Universe. This analysis can be done on any other celestial body observed from the Earth. Since Sun and Mars are chosen arbitrarily, and there is nothing special about Mars, one can expect to come up with the same general conclusion.
There is another interesting remark that follows from this analysis. If one could put the whole Universe in accelerated motion around the Earth, the pseudo-potential corresponding to pseudo-force (4.2) will immediately be generated. That same pseudopotential then causes the Universe to stay in that very state of motion, without any need of exterior forces acting on it.
Acknowledgments
This work is supported by Ministry of Science, Sports and Technology under contract number 119-0982930-1016.
References
[1] Koestler A 1959 The Sleepwalkers: A History of Mans Changing Vision of the Universe (Hutchinson: MacMillan), pp. 194-195
[2] Rawlins D 1987 “Ancient heliocentrists, Ptolemy, and the equant” Am. J. Phys. 55 235239 [3] Rosen J. 1981 “Extended Mach principle” Am. J. Phys. 49 258-264 [4] Newburgh R 2007 “Inertial forces, absolute space, and Machs principle: The genesis of relativity”
Am. J. Phys. 75 427430 [5] Lichtenegger H and Mashhoon B 2004 “Machs principle” arXiv:physics/0407078 [physics.hist-ph] [6] Hood C G 1970 “A Reformulation of Newtonian Dynamics” Am. J. Phys. 38 438442 [7] Barbour J 1974 “Relative-distance Machian theories” Nature 249 328329 [8] Assis A K T 1999 Relational Mechanics (Montreal: Aperion) [9] Hartman H I and Nissim-Sabat C 2003 “On Machs critique of Newton and Copernicus” Am. J.
Phys. 71 11631168 [10] Mach E 1872 Die Geschichte und die Wurzel des Satzes von der Erhaltung der Arbeit, (Prague:
Calve) [11] Mach E 1883 Die Mechanik in ihrer Entwickelung historisch-kritisch dargestellt, (Leipzig:
Brockhaus) [12] Mach E 1911 History and Root of the Principle of the Conservation of Energy (Chicago: Open
Court) [13] Rovelli C 2004 Quantum Gravity (Cambridge University Press) p 75 [14] Barbour J 2010 “The definition of Machs principle” arXiv:1007.3368 [gr-qc] [15] Zylbersztajn A 1994 “Newtons absolute space, Machs principle and the possible reality of fictitious
forces”, Eur. J. Phys. 15, 18 [16] Hauser W 1985 “On planetary motion” Am. J. Phys. 53 905907
Newton-Machian analysis of Neo-tychonian model of planetary motions
11
[17] Gauthier N 1986 “Planetary orbits” Am. J. Phys. 54 203 [18] Landau L D and Lifshiz E M 1976 Mechanics 3rd ed (Oxford: Butterworth-Heinemann) pp. 2540 [19] Goldstein H 1980 Classical Mechanics 2nd ed (Addison-Wesley) pp. 70102 [20] Handbook of Chemistry and Physics (49th ed) ed Weast R C 1968 (Cleveland Ohio: The Chemical
Rubber Company), pp. F 145146 [21] Thompson B G 2005 “Using retrograde motion to understand and determine orbital parameters”
Am. J. Phys. 73 10231029 [22] Vet˝o B 2011, “Gravitomagnetic field of the universe and Coriolis force on the rotating Earth” Eur.
J. Phys. 32 13231329