Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Math light bending by gravity

  1. Mar 3, 2007 #1
    Hey all.
    I was wondering what is the math behind calculating how much light is bent by gravity/curvature of space. Can someone teach me??

  2. jcsd
  3. Mar 3, 2007 #2
    Are you interested in the the bending of starlight from conservation of energy and angular momentum, or do you want to see the relativistic approach?
  4. Mar 3, 2007 #3
    If your question pertains to the amount of space curvature for a given amount of matter, the excess radius is MG/3c^2 where M is the mass of the object (say earth), G is the gravitational constant and c^2 is the square of the velocity of light. If you put in the numbers, it is seen that the entire mass of the earth only modifies the effective radius by about 1.5 millimeters
  5. Mar 4, 2007 #4
    I'd prefer the relativistic approach.
  6. Mar 4, 2007 #5
    I kind of don't know what your saying. What does that equation do:confused:? I'm curious now.

    Anyway my real question is how can you calculate the bending of light passing a massive object, for example the sun??
  7. Mar 4, 2007 #6
    Sure! I'd like to know how that gives you the result, could you at least outline it?
  8. Mar 4, 2007 #7


    User Avatar
    Staff Emeritus
    Science Advisor

    You'll need at least basic differential equations. There's a series of formulas using different variables in for instance, MTW "gravitation" pg 673 for the orbit of light in the Schwarzschild geometry (that is, near a single large mass).

    The orbit is determined by a single constant, the "impact parameter" b.

    It has various forms, it will look something like the following differential equation which expresses r as a function of [itex]\phi[/itex] in the [itex]\theta=0[/itex] plane. r, [itex]\theta[/itex], and [itex]\phi[/itex] are all Schwarzschild coordinates, which are derived from the familar "spherical coordinates" by the same names.

    \left( \frac{1}{r^2} \frac{dr}{d\phi} \right)^2 + \frac{1-2M/r}{r^2} = \frac{1}{b^2}
    Last edited: Mar 4, 2007
  9. Mar 5, 2007 #8
    Hey, cesiumfrog,
    Thanks for asking. Here's how I got the experimentally determined value of bending using angular momentum and energy conservation http://alzasol.spaces.live.com/blog/?&_c02_owner=1 [Broken]

    There's a drawing at the bottom of the text. You'll have to klick on it.
    Last edited by a moderator: May 2, 2017
  10. Mar 5, 2007 #9
    There is a lot of math that goes into the derivation of the deviation of a light beam passing by the surface of a massive object - during the development of GR, Einstein had first reached a wrong conclusion about the deviation - later he realized his error and published the correct deviation which has been verified by subsequent experiments - the correct value is twice that which would be calculated using Newtonian approximations - in a recent thread, I think it was Space Tiger that posted a clever Newtonian derivation using the notion of an impulse function which is a short cut way of arriving at would happen if one calculates the force at each instant using Newtonian gravity. But even these methods require a bit of mathematical maturity to be able to appreciate what is being said -

    The bottom line is that the spatial deviations are very small even for massive bodies - for the earth the change in the spatial curvature to which I referred is very small - and in the case of the Sun the deviation of light grazing its surface is was barely measurable in 1919 when the theory of GR was first put to the test. But according to GR, both space and time are affected by matter - and it is not at all difficult to set up experiments that measure the temporal changes - for example GPS satellites are corrected for altitude before launch - the change in clock rates for even small differences in height are easily measured and routinely correct in order to synchronize the satellite clocks with GPS receivers on the surface of the earth
  11. Mar 11, 2007 #10
    I saw that document. There is something I would like to ask you: is it really possible to write c^2 = c_x^2 + c_y^2? If a body has velocity components v_x = 0.9c and v_y = 0.9c, then its speed should be c*Sqrt(2*0.9^2) = 1.27c > c So, it's impossible for a body to have such velocity components?
    Last edited by a moderator: May 2, 2017
  12. Mar 11, 2007 #11


    User Avatar
    Staff Emeritus
    Science Advisor

    Unfortunately, the whole derivation in that document is a mess as far as getting any accurate answers go.

    If you consder the energy-momentum 4-vector, it does turn out that there is a conserved energy [itex]p_0[/itex] and, assuming that the orbit is in the [itex]\theta=0[/itex] plane, a conserved anguler momentum [itex]p_\phi[/itex], but the correct relativistic forms of these conserved quantities aren't given by the URL above, which is someones ambitious but unfortunately not technically correct attempt at solving the problem.

    One good online source that gives the right answers for the orbits of material bodies http://www.fourmilab.ch/gravitation/orbits/, though it doesn't really go into the correct derviation of the formulas for energy and angular momentum.

    Unfortunately it doesn't directly handle the OP's question of the orbit of light, though one can get that orbit by taking the orbit of a material body in the limit as it approaches 'c'.

    As far as lightarrow's question about velocities goes, in flat space-time it is true that v^2 = (dx/dt)^2 + (dy/dt)^2, however the space-time for the Schwarzschild geometry isn't flat, and this formula doesn't work right there.
  13. Mar 11, 2007 #12


    User Avatar
    Staff Emeritus
    Science Advisor

    http://arxiv.org/abs/gr-qc/0506075" [Broken] is another reference for a correct way of deriving the equations of motion of a particle in the Schwarzschild geometry near a large mass.

    Basically, the easiest way to get the right answer is to use a Lagrangian approach, and to find the path that extremizes "proper time" - this has also been called the "principle of maximal aging".

    http://www.eftaylor.com/pub/GRtoPLA.pdf is also useful in describing the idea, but rather than find an exact solution it demonstrates that the principle of least action (maximal aging) yields the Newtonian answer in the non-relativistic limit.

    There are some more references at http://www.eftaylor.com/leastaction.html#actionsummary

    about the general idea of the "least action" approach to physics.
    Last edited by a moderator: May 2, 2017
  14. Mar 11, 2007 #13
    Yeah, you are right.
  15. Mar 11, 2007 #14
    One interesting thing here to ponder on is, since the proper time of a light beam is always zero how would we find a single path that extremizes proper time. :wink:

    Would you agree that that approach has some issues?
  16. Mar 13, 2007 #15

    George Jones

    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    A geodesic can be parametrized by something other than proper time.
  17. Mar 13, 2007 #16
    Which, including their physical interpretations, do you consider good parameterization options for null curves?
  18. Mar 13, 2007 #17
    The mistake you're making is an implicit assumption that the action principle for the equations of motion of a particle are unique. This is obviously untrue. Besides the obvious choices for massive particles where one supposes that the value of the action is proportional to the proper length, there are myriad other possibilities. One can, if one wishes to be fancy about it, construct an action principle by demanding that the action be stationary with respect to harmonic functions of compact support defined along the world line. Alternatively, you could just go the old route of introducing an einbein along the world line and demanding that the action be stationary with respect to variations in it; this is covered, for example, in the first string theory lecture every student receives as a motivating example for the Brink-Di Vecchia-Howe action for a bosonic string.
  19. Mar 13, 2007 #18


    User Avatar
    Staff Emeritus
    Science Advisor

    Another approach I've seen is to simply take the limit for the path of a massive particle as m->0. For instance, a neutrino has mass according to recent observations, but the path it follows is essentially the same path that light follows.

    Using an einbein (sidenote: we've disucssed einbein's in https://www.physicsforums.com/showthread.php?p=1054196 ) has a lot of advantages, too. Even in the case of a massive particle, taking the direct route of minimizing the action (which can be expressed as the square root of a quadratic form, because proper time squared is a quadratic form) gets very messy due to the square root.

    However, while it's convenient computationally, einbein's may not be as familiar to people as the approach using the calculus of variations, so it might not be the first choice if one is trying to present the material at an undergraduate college level. Unfortunately, I don't think there is any way to approach the topic using only high-school level math - one really needs calculus!
  20. Mar 13, 2007 #19
    Pervect, I particularly have trouble with the word "essentially" in this context. To me it is far from obvious that we can say that. I would be more inclined to say that a lightlike path is essentialy different from a timelike path.

    Now if someone would refer me to some paper that proves that a lightlike path is simply the limit of a timelike path on a manifold with a Minkowski like metric, then I would be happy. But honestly, I am not holding my breath. :smile:

    Actually, while we are at it, I would not be surprised that someone conjured up a proof that the contrary is true!
    Last edited: Mar 13, 2007
  21. Mar 14, 2007 #20

    George Jones

    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    Not all good parameterizations have obvious physical interpretations, but for null curves, one that that does is the parameterization that makes the tangent 4-vector to the curve equal to the four-momentum.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook