I am interested in the effects of the discrete symmetry operation that reverses the time-orientation of a manifold, i.e., flipping past by future. To give some context, I'm thinking about classical field theories in curved spacetimes. Let me start by stating the things I know:
Time-reversed worldlines
A worldline $z^\alpha(\tau)$ is a map that takes a real number (the proper time $\tau$) and returns a spacetime point in the manifold $\mathcal{M}$
$$ \tag{1} z: \mathbb{R}\rightarrow \mathcal{M} $$
The full path that $z$ "draws" on the manifold is invariant under time reversal, but the direction in which we traverse it changes. I'll use a bar to denote "time-reversed" quantities, then
$$ \tag{2} \bar{z}^\alpha (\tau)=z^\alpha(-\tau) $$
Similarly, the 4-velocity satisfies
\begin{equation} \begin{aligned} \bar{u}^\alpha (\tau)\equiv& \frac{d\bar{z}(\tau)}{d\tau}\\ =&-u^\alpha(-\tau) \end{aligned} \end{equation}
Time-reversed Green functions
Consider a Green function satisfying some simply equation
$$ \tag{4} D_x G(x,y)=\delta(x,y)$$
where $D_x$ is some nice differential operator like $\nabla_\mu \nabla^\mu$ or such. Let's focus on the retarded and advanced solutions, which satisfy
- $G_{\text{ret}}(x,y)=0$ when $y$ is not in the past of $x$.
- $G_{\text{ret}}(x,y)=G_{\text{adv}}(y,x)$.
The retarded and advanced solutions clearly depend on the time-orientation of the manifold, since they distinguish past and future. It follows from their causality that they are time-reversed of each other
$$ \tag{5} \bar{G}_{\text{ret}}(x,y)=G_{\text{adv}}(x,y) $$
Confusion about Green functions evaluated on worldlines
Consider now a retarded Green function evaluated along a wordline $G_{\text{ret}}\big[z(\tau_1),z(\tau_2)\big]$. I would think that time-reversing the worldlines would be equivalent to time-reversing the Green function itself
$$ \tag{6} G_\text{ret}\big[z(-\tau_1),z(-\tau_2)\big]\stackrel{?}{=}G_\text{adv}\big[z(\tau_1),z(\tau_2)\big]$$
but I don't think that is true, at least not in general spacetimes that are not homogeneous.
I believe it would be true only if we flip time orientation specifically along the mid point $\tau_0=\frac{\tau_1+\tau_2}{2}$ between $\tau_1$ and $\tau_2$
$$\tag{7} G_\text{ret}\big[z(\tau_0+\Delta \tau),z(\tau_0-\Delta \tau)\big]=G_\text{adv}\big[z(\tau_0-\Delta \tau),z(\tau_0 +\Delta \tau)\big]$$
which essentially corresponds to flipping the two arguments, just as in the previous section.
Deeper confusion about time-reversed higher $n$-point functions
It is common in perturbation theory to find higher-order $n$-point functions that are sourced by nonlinearities, such as
$$ \tag{8} D_x J(x,y,z)=G(x,y)G(x,z)$$
which can admits retarded and advanced solutions
$$ \tag{9} J_{\text{ret/adv}}(x,y,z)=\int dV' G_{\text{ret/adv}}(x,x')G_{\text{ret/adv}}(x',y)G_{\text{ret/adv}}(x',z)$$
In this case, I am not even sure a property such as $G_{\text{ret}}(x,y)=G_{\text{adv}}(y,x)$ exists...
It is clear that flipping time orientation and exchanging arguments produce the same result for 2-point functions. The question is, how does the relation between these two operations generalize to higher order n-point functions?