To study stable bound orbits, we regard the geodesic motion of massive and massless particles as a two-dimensional potential problem (see ref.[1]} about the detail). The Hamiltonian of a free particle with mass \(m\)  is written as
\(\mathcal {H}=g^{\mu \nu }p_\mu p_\nu +m^2, \) 
