Asymptotic behavior of geodesics in conformally compact manifolds

We study the asymptotic behavior of geodesics near the boundary of a conformally compact Riemannian manifold $(X,g)$. In the case where the sectional curvature at infinity is constant (the asymptotically hyperbolic case) it is known that non-trapped geodesics extend to the conformal infinity as smoothly immersed curves in $\overline{X}$ and that the conformal infinity can be locally smoothly parametrized by the initial directions of such non-trapped geodesics starting at a given point. We show that in the general case non-trapped geodesics typically only extend to the conformal infinity with $C^{1,α}$ regularity, due to the presence of a logarithmic singularity sourced by the gradient of the limiting sectional curvature. In spite of this singular behavior, we show that the endpoints of such geodesics depend smoothly on the initial conditions, so that the conformal infinity can still be locally smoothly parametrized using the initial velocities of geodesics. In particular, the smooth structure of the conformal infinity is canonically determined by $(X,g)$. We also show how to specify `initial data' for geodesics at the conformal infinity and obtain an asymptotic expansion that parametrizes the geodesics tending toward a given boundary point.