An exact solution to the gravitational bending of starlight is derived without resorting to general relativity. The accurate solution can be reduced to the result as obtained by general relativity under special conditions. For objects with a larger mass and a smaller radius, the approximated result as obtained by Schwartzchild’s metric based on general relativity will be invalid. A century-long misunderstanding about the deflection of light by gravity based on the classical theory will also be clarified.
Gravitational deflection of starlight by a massive object such as the Sun, causing the light trajectory to bend from a straight line, has been an interesting topic for over three centuries [
It has been a main-stream understanding that, Newton’s classical theory of gravity yields a deflection angle of about 0.87" for light beam from a remote star passing near the surface of the Sun, and the deflection angle as calculated by Schwartzchild’s metric based on general relativity is 1.75", twice the value as obtained by the classical theory [
In order to demonstrate the theoretical predictions, an expedition led by Eddington during the total solar eclipse in 1919, reported the observation of light deflection with an angle consistent with general relativity [
Along with the successful explanation for the anomalous precession of the perihelion of Mercury, direct observations of light bending by gravity in the trajectory so as to compare the observed deflection angle with the theoretical results will be a strong supplement to fully appreciate the acclaimed theory of general relativity. Due to the profound nature of space-time curvature in the theory, calculations for the deflection angle by Schwartzchild’s metric based on general relativity require complicated mathematics [
We shall present detail formulations with an exact solution for the trajectory of light from a remote background star grazing a massive object, without resorting to general relativity. We shall demonstrate that, the approximated result for the deflection angle obtained by Schwartzchild’s metric based on general relativity is valid only when the mass to radius ratio ( M / R ) is small, and the approximated result will deviate largely from the accurate value when the mass to radius ratio of the massive object is close to or larger than one.
We start from Newton’s law of gravity. The gravitational force on a mass m0 from a central mass M is F = G m 0 M / r 2 , with G the gravitational constant, and r the distance between the two objects. Based on the mass-energy relation ε = m 0 c 2 and Plank’s principle of radiation quantization, light (or photons) of energy h ν will exhibit an effective mass of m 0 = h v / c 2 , where c is the velocity of light, h is Plank’s constant. By solving the related orbital equation under proper boundary conditions, an exact deflection angle of α = 2 arcsin [ 2 k 0 / ( 1 + 2 k 0 ) ] can be obtained, with k0 a constant related to the speed of light and the solar mass and the solar radius. Following is the detail derivation.
Assuming the origin O of polar coordinates is located at the centre of the massive object such as the Sun as illustrated in
where H is the angular momentum per unit mass, H = r 2 ( d φ / d t ) , M is the mass of the Sun, and the photon mass m0 is not involved in the equation.
The general solution to the above equation is
with u 1 = G M / ( 2 H 2 ) = u 0 / 2 = 1 / p 1 , H 2 = p 1 ( G M / 2 ) . Based on the boundary condition and asymptotic behavior at specific angles and positions, the constants A and B can be determined.
At the closest spot S from the light beam to the surface of the Sun, φ = π/2, r = R, with R the radius of the Sun, or the closest distance between the light beam and the center of the Sun, the orbital equation yields, u = A cos φ + B sin φ + u 1 = B + u 1 = 1 / R , and
For the light beam traveling from a remote star towards the Sun, it has already been deflected before reaching the Sun. As a result, under the condition φ = π, r is not infinite ( r ≠ ∞ ). Only when φ is larger than π by a deflection angle α1, which is the angle between the line OQ and the X-axis, so that the asymptotic line PE tangential to the hyperbola of the light beam at the negative infinite distance (from the remote star) tends to be parallel to OQ, with φ = π + α 1 , then r = ∞ , and u = 0, the orbital equation yields
A cos α 1 + B sin α 1 − u 1 = 0 . (4)
The same consideration also applies for the light beam travelling away from the Sun to the other side of infinite distance (such as the Earth). Under the condition φ = 0, r is not infinite ( r ≠ ∞ ), rather a finite value, only when φ is smaller than 0 by a deflection angle α1 with φ = 0 − α 1 = − α 1 , then r = ∞ , and u = 0, the orbital equation yields
A cos α 1 − B sin α 1 + u 1 = 0 . (5)
Combining the above two conditions of Equation (4) and Equation (5), B sin α 1 − u 1 = 0 , and
sin α 1 = u 1 / B = u 1 / ( 1 R − u 1 ) . (6)
The eccentricity of the hyperbola orbit is ε = [ 1 + ( E / m 0 ) H 2 / ( G M ) 2 ] 1 / 2 , with E the total energy, H 2 = p 1 ( G M / 2 ) , and u 1 = G M / ( 2 H 2 ) = 1 / p 1 . At the remote infinite distance, the gravitational potential energy is zero, and E = m0c2. At the closest spot S from the light beam to the Sun, φ = π / 2 , R = p 1 / ( 1 + ε ) . So,
ε = [ 1 + ( m 0 c 2 / m 0 ) H 2 / ( G M ) 2 ] 1 / 2 = [ 1 + c 2 H 2 / ( G M ) 2 ] 1 / 2 = [ 1 + 1 2 c 2 p 1 G M / ( G M ) 2 ] 1 / 2 = { 1 + 1 2 R c 2 ( 1 + ε ) / ( G M ) } 1 / 2 = [ 1 + 1 2 ( 1 + ε ) / k 0 ] 1 / 2
where k 0 = G M / ( R c 2 ) , which is a constant related to the speed of light, and the solar mass and the solar radius.
Solving the above quadratic equation we get, ε 2 − ( 1 2 k 0 ) ε − ( 1 + 1 2 k 0 ) = 0 , and
ε = { 1 2 k 0 + [ ( 1 2 k 0 ) 2 + 4 ( 1 + 1 2 k 0 ) ] 1 / 2 } / 2 = 1 2 k 0 + 1 . (7)
The deflection angle α1 at one side of the Sun is
sin α 1 = 1 / ε = 2 k 0 / ( 1 + 2 k 0 ) . (8)
The total deflection angle at both sides is
α = 2 α 1 = 2 arcsin [ 2 k 0 / ( 1 + 2 k 0 ) ] . (9)
This solution is analytically exact without simplification, with related assumptions for the photon energy of the form m0c2, and the gravitational force of the form Gm0M/r2.
Under the small angle approximation, with a small value of k0 for the Sun of about 10−6, which is much smaller than one, then sin α 1 ~ α 1 , 1 + 2 k 0 ~ 1 , sin ( α 1 ) ~ α 1 = 1 / ε = 2 k 0 , the total deflection angle reduces to
α = 2 α 1 = 2 arcsin [ 2 k 0 / ( 1 + 2 k 0 ) ] ~ 2 arcsin ( 2 k 0 ) = 4 k 0 . (10)
The approximated deflection angle of 4k0 is the same as given by general relativity based on Schwartzchild’s metric, and equal to 1.751" for the Sun. Most of the reported values for the deflection angle by the Sun are about 1.7" [
The above derivation procedure for the light trajectory is the same as solving Kepler’s orbital problem in the fundamental mechanics [
For nearly a century, it has been a widely accepted opinion that, the deflection angle based on the classical Newtonian theory is half the value as given by general relativity [
One of the errors during the conventional calculation based on the classical Newton’s theory is mistakenly ignoring one side of the deflection for incoming light from the remote star, resulting in half of the deflection angle, with another mistake of a wrong asymptotic behavior in the orbit, as shown in
based on the classical theory.
A simple way to modify the mistakes in the conventional calculation based on the classical theory, is to multiply the deflection angle by a factor of two, resulting in a deflection angle of 4k0, which is the same as that obtained by general relativity. In fact, in the calculations based on general relativity, Schwartzchild’s metric also considered a multiplying factor of two.
The reason for the good agreement between the observations and the theoretical value of 4k0 as given by general relativity for light passing near the Sun, if the observations are reliable and other deflection mechanisms are properly excluded, is due to the specific parameters of the Sun. Based on the mass and the radius of the Sun, k 0 = G M / ( R c 2 ) ~ 10 − 6 ≪ 1 , then 2 k 0 / ( 1 + 2 k 0 ) ~ 2 k 0 , the deflection angle of 4k0 as given by general relativity will be highly accurate.
For other stars or galaxies of much larger masses and smaller radii, with k0 close to one or even larger than one, the deflection angle of 4k0 as given by general relativity will deviate largely from the accurate value as given above. Observations of starlight deflected by a massive objects with a high density such as neutron stars, quasars, or black holes, can be used to test the above differences.
Since distortions of radio signals caused by the plasma fluctuations in the solar corona has been observed [
In summary, the deflection angle of light from a remote star grazing a massive object has been derived with an exact solution by solving Binet’s orbital equation, based on the classical Newton’s law of gravity F = Gm0M/r2, with a photon energy in the form of m0c2. The present solution is more accurate than the result as given by general relativity based on Schwartzchild’s metric. As a key candidate for testing the effectiveness of general relativity, direct observations of light bending in the trajectory by the Sun with an accurate deflection angle, or preferably by other massive objects with a known mass and a larger mass to radius ratio, will continue to be a challenge for experimental studies.
Huang, F.Y. (2017) Accurate Solution to the Gravitational Bending of Starlight by a Massive Object. Journal of Modern Physics, 8, 1894-1900. https://doi.org/10.4236/jmp.2017.811112