I thought that the source of this 2 is the square-root appearing in the proper time formula:
[tex]dT = \sqrt{g_{\mu,\nu} dx^{\mu} dx^{\nu}}[/tex]
In weak, time-independent gravitational fields, [itex]g_{\mu,\nu}[/itex] takes the form
[tex]g_{0,0} = 1 + h[/tex]
[tex]g_{1,1} = g_{2,2} = g_{3,3} = -1[/tex]
[tex]g_{\mu,\nu} = 0 \quad \quad (\mu \ne \nu)[/tex]
So using
[tex]\left(\frac{v}{c}\right)^2 = \left(\frac{dx^1}{dx^0}\right)^2 +\left(\frac{dx^2}{dx^0}\right)^2 +\left(\frac{dx^3}{dx^0}\right)^2[/tex]
the formula for proper time becomes
[tex]dT = \sqrt{1+h - \frac{v}{c}^2} dt[/tex]
For [itex]A[/itex] and [itex]v[/itex] small, we can expand the square-root to get
[tex]dT = \left( 1 + \frac{1}{2} h - \frac{1}{2} \left(\frac{v}{c}\right)^2 \right) dt[/tex]
The equations of motion for a test particle are obtained by maximizing the integral of [itex]\left(1 + \frac{1}{2} h - \frac{1}{2} \left(\frac{v}{c}\right)^2 \right) dt[/itex]. Since the constant 1 is irrelevant for the equations of motion, this is the same as maximizing the integral of [itex]D dt[/itex], where
[tex]D = \frac{1}{2} h - \frac{1}{2} \left(\frac{v}{c}\right)^2[/tex]
But, in this limit (small [itex]A[/itex] and small [itex]v[/itex]), Newtonian mechanics can be used, so we should also be able to obtain the equations of motion by minimizing the integral of [itex]L dt[/itex], where
[tex]L = \frac{1}{2} mv^2 - m \Phi[/tex]
(where [itex]\Phi[/itex] is the Newtonian gravitational potential). This suggests that [itex]D[/itex] and [itex]L[/itex] are linearly related. The only way for this to be true is if
[tex]h = 2 \Phi/c^2[/tex]