In general relativity the energy of a test-body moving in a spherically symmetric gravitational field can be written as:
$$E=mc^2\left(\frac{\sqrt{1-\frac{2GM}{rc^2}}}{\sqrt{1-\frac{v^2}{c^2\left((1-\frac{2GM}{rc^2})^2(\hat{r}\cdot\hat{v})^2+(1-\frac{2GM}{rc^2})|\hat{r}\times\hat{v}|^2\right)}}}\right).$$
The strange part is just because the speed of light is not the same in the radial direction as it is transverse to the radial direction. To be precise I mean velocity in coordinate time and treat the r-parameter of the Schwarzschild solution as a real distance. For a pure non-radial motion the expression above becomes a little bit simpler:
$$E=mc^2\left(\frac{\sqrt{1-\frac{2GM}{rc^2}}}{\sqrt{1-\frac{v^2}{c^2\left(1-\frac{2GM}{rc^2}\right)}}}\right). $$
This can be rewritten as:
$$E=mc^2\left(\frac{{1-\frac{2GM}{rc^2}}}{\sqrt{1-\frac{2GM}{rc^2}-\frac{v^2}{c^2}}}\right).$$
For a circular motion you have $v=\sqrt{GM/r}$ in general relativity as well as classically, and by replacing $v^2$ with $GM/r$ above you get an expression where you find that you need infinite energy to stay in a circular orbit at the photon sphere radius and by taking the derivative of that expression with respect to "r" you find the innermost stable radius.
If you define $c(r)=c\sqrt{1-\frac{2GM}{rc^2}}$ and $m(r,v)=\frac{m}{\sqrt{1-\frac{2GM}{rc^2}-\frac{v^2}{c^2\left(1-\frac{2GM}{rc^2})(\hat{r}\cdot\hat{v})^2)+|\hat{r}\times\hat{v}|^2\right)}}}$
which for pure non-radial motion reduces to: $$m(r,v)=\frac{m}{\sqrt{1-\frac{2GM}{rc^2}-\frac{v^2}{c^2}}}.$$
you can write: $E=m(r,v)(c(r))^2$
Question: Can you write, is it correct to write, the energy of a test-particle moving under Schwarzschild conditions in the form $E=m(r,v)(c(r))^2$ like I have outlined above?
There is a specific reason why I am asking and that is that I want to generalize Newtons gravitational law by writing:
$$\frac{d(m\gamma\bar{v})}{dt}=-\frac{GMm\gamma}{r^2}\hat{r}$$ and when I use $\gamma=\frac{1}{\sqrt{1-\frac{v^2}{c^2}}}$ I can only explain one third of the "anomalous perihelion shift". In order to explain all of it I can use $\gamma=\frac{1}{\sqrt{1-\frac{2GM}{rc^2}-\frac{v^2}{c^2}}}{\frac{1}{\sqrt{1-\frac{2GM}{rc^2}}}}$. If you do this you get an expression that looks like:
$$\frac{d\bar{v}}{dt}=-\frac{GM}{r^2}\left(\hat{r}-3\frac{v^2(\hat{r}\cdot\hat{v})\hat{v}}{c^2(1-\frac{2GM}{rc^2})} +\frac{v^4(\hat{r}\cdot\hat{v})\hat{v}}{c^4(1-\frac{2GM}{rc^2})^2}\right)$$.
(Note: here $\frac{\bar{r}}{r}=\hat{r}$ and $\frac{\bar{v}}{v}=\hat{v}$)
(Edit: Maybe it possible to find the acceleration by the condition that the time derivative of the energy expression must be zero, given that the expression for the energy is correct.)
The resulting orbits from using this expression is close to what is expected from GR but not quite right, I think it could be improved. In the plot of numerical integrations using the last expression below the green circle represents the Schwarzschild radius and the red circle represents the innermost stable circular orbit.
