love_42 said:
The equation
ai = -Γijk vj vk
is perfectly capable of explaining the deflection of light by the sun.
Agreed. But calculating the deflection requires, as others have said, is a more involved process than the shortcuts you are apparently taking.
One can calculate the deflection in Schwarzschild coordiantes, but it doesn't give as much insight into the cause as using PPN coordinates.
A bit of background. The PPN paramter ##\gamma## is described as:
\gamma How much space curvature ##g_{ij}## is produced by unit rest mass
There is another nonzero PPN paramter ##\beta## in GR, but it doesn't affect the deflection of light.
The details of the calculation can be found in MTW's (Misner, Thonre, Wheeler) "Gravitation", pg 1102, though you'll need to go to the earlier pages to get the PPN metric. Which is a bit involved, as PPN as described in MTW is used to accommodate any metric theory of gravity to allow a comparsion framework for multiple theories. However, GR has only two nonzero PPN paramters, ##\beta## and ##\gamma##, and I'll give the metric coefficients for those at the end of the post.
Using the geodesic equation, one can calculate the trajectory of the light path in the ##r, \phi## plane that light takes:
$$\frac{b}{r} = \sin \, \phi + \frac{\left( 1+\gamma \right) M}{b} \left( 1 - \cos \phi \right)$$
That's just from the textbook, if you don't trust the textbook, you'll probably have to wade through it on your own, but it should be possible by solving the geodesic equations.
b can be regarded as the "impact parameter". It is approximately the distance of closest approach, which is
$$r_{min} = b\left[ 1 - \frac{ (1 + \gamma) M}{b}\right]$$
This can be expanded in a series form around around ##\phi=0##, which will be regard as the ingoing ray, and around ##\phi=\pi##, which will be regard as the outgoing ray:
$$\frac{b}{r} = \phi + O(\theta^2) \quad \frac{b}{r} = \frac{2 \left(1+ \gamma \right) M}{b} + \left( \pi - \phi \right) + O((\pi - \phi)^2)$$
The end result is that the spatial curvature term ##\gamma=1## doubles the deflection of the light for small angles ##\phi## compared to a hypothetical deflection for a metric theory with no spatial curvature. An example of such a theory would be the Newtonian approximation to GR which has no spatial curvature. This approximation is also discussed in MTW.
The metric elements for the Newtonian approximation in geometric units where G=c=1 are:
$$g_{00} = -1 + 2U \quad g_{ij} = \delta_{ij}$$
The metric elements for the PPN approximation in geometric units are
$$g_{00} = -1 + 2U - 2\beta U^2 \quad g_{ij} = \delta_{ij}(1+2\gamma U)$$
where U = M/r, and which can be regarded as the Newtonian potential energy. Also ##\delta_{ij} =## 1 if i=j and 0 otherwise
This is of course a lot of work. You'll most likely need to consult a textbook (such as MTW) along with a lot of time to go through it all.
But it is a textbook description of why we say that "spatial curvature" causes the doubling of the deflection of light. Implicit in talking about "spatial curvature" is a specific time slicing. The time slicing used in PPN and in the Schwarzschild metric are the same, they're the natural choice for the split which make the metric coordinate expressions time independent.