Reference: Moore, Thomas A., A General Relativity Workbook, University Science Books (2013) – Chapter 13; Problem P13.7.
In addition to being deflected when it passes a mass, photons also experience a delay (relative to the travel time in the absence of the deflecting mass) known as the Shapiro delay, after Irwin Shapiro who first observed the effect in the 1960s. Crudely speaking, the delay is due to the curvature of space-time resulting in the photons having to travel a longer path, although of course the actual reason is a bit more subtle, being related to the revised notions of space and time in the Schwarzschild metric.
Our starting point for calculating the delay in the travel time is the photon radial equation of motion:
The closest approach of the photon to the mass occurs when is a minimum at , and thus when . Thus
We can express the original differential equation in terms of by eliminating .
Trying to integrate this is impossible as it looks (except numerically), so we need to simplify it by making a few approximations. As usual, we can assume that (that is, the photon never gets close to the Schwarzschild radius ). As a first attempt, we might just ignore these terms (that is, take ). If we do that, we get
where is the time at which the photon is at the closest approach . Thus is the time to travel to or from the point of closest approach to some point in flat space (since we’ve taken ), so this obviously can’t include the delay factor. We need a better approximation, and one that is hopefully still integrable.
We’ll introduce the variable and rewrite the equation above:
Since we’re assuming and , as well, so we can expand the operand of the square root in a Taylor series and keep only the first order term. We define to simplify the notation.
Thus we’re now faced with
(Note: the minus sign in the last square root is correct; the plus sign given by Moore in equation 13.27 is wrong. See the list of errata in the book.)
Needless to say, this integral still can’t be done, but we can approximate yet again by noting that
Thus to first order, we now have
These integrals can be done (using software or tables), and we get (substituting back for ):
The first term in this answer is the flat space term, so the latter two terms represent the delay. Both these terms are non-negative (since the argument of the logarithm is so it is non-negative).
As an example, we can work out the delay for a photon travelling from Venus to Earth when Venus is on the opposite side of the sun from Earth. The radius of Venus’s orbit is 0.723 AU (astronomical units) and the radius of the Earth’s orbit is 1 AU. If the photon just grazes the surface of the sun, is the sun’s radius which is . for the sun is 1.477 km.
The delays in travelling from Venus to the Earth are then (divide by the speed of light to convert from km to sec):