In Minkowski, the metric covariant derivative is just $\nabla = \partial$, and therefore geodesics (solutions $\gamma(t)$ to the geodesic equation $\nabla_{\dot\gamma}\dot\gamma = 0$) are just straight lines. However, when considering a background electromagnetic field, the covariant derivative is given by the connection over the U(1) bundle so that $D_\mu = \partial_\mu + ieA_\mu$, with $A_\mu$ the electromagnetic vector potential.
I am having trouble in seeing the behavior of geodesics in this space. Are they still straight lines? Can one use the geodesic equation $$D_{\dot \gamma} \dot \gamma =0 ?$$
Having an $i$ in the $D_\mu$ makes everything even weirder, as it feels like this 'geodesic' will take complex values.
My guess from my very basic understanding of bundles is that this geodesic should be studied not in my flat spacetime $\mathbb{R}^4$, but rather on the total space $U(1)\times \mathbb{R}^4$, and somehow what I should be interested in is the projection of this geodesic onto my manifold, or something like that.
Context of the question:
In QFT in curved spacetimes (QFTCS) the vacuum expectation value of observables such as $\phi(x)^2$ is defined through a point splitting prescription, i.e. $$ \left<\phi(x)^2\right> := \lim_{x'\to x} \langle0\rvert\phi(x)\phi(x')\lvert 0\rangle. $$ This prescription can be extended to the stress-energy tensor by parallel transporting the $T_{ab}$ along the (unique) geodesic that joins x' and x.
In my case, I am studying the charge current density $\left<j_\mu(x)\right> $ on flat spacetime with a background EM following the same procedure. I want to know if this geodesic changes due to having to do this parallel transport, and along which path it is done.
I also don't really understand what exactly parallel transport in that limit means.