Let $\hat{U}$ be a unitary transformation that acts on fermionic fields (defined on a lattice) as $$\hat\psi_j \rightarrow \hat{U}\hat\psi_j\hat{U}^{-1},$$ $$\hat\psi_j^\dagger \rightarrow \hat{U}\hat\psi_j^\dagger\hat{U}^{-1}.$$
I consider it to be anti-unitary, just because it makes more explicit my argument, so $$i\rightarrow \hat{U}\hat i\hat{U}^{-1}=-i,$$ where $i$ is the imaginary unit.
Now I want to understand how this transformation modify correlation functions of the type $\left<\hat\psi_i\hat\psi_j^\dagger\right>$. I can think about this in two ways, but I don't really know which one is the correct approach:
Since $\left<\hat\psi_i\hat\psi_j^\dagger\right>$ is just a scalar, one way of proceeding is to only use the anti-unitary property of the operator, i.e. $$\left<\hat\psi_i\hat\psi_j^\dagger\right>\rightarrow \hat{U}\left<\hat\psi_i\hat\psi_j^\dagger\right>\hat{U}^{-1}=\left<\hat\psi_i\hat\psi_j^\dagger\right>^{*},$$ so the transformation only complex-conjugate the correlation function. This does not look good to me, because e.g. for a Gaussian states different transformations which only share the anti-unitary property would modify the state, determined thy the correlation matrix, in the same way, and there would be no difference between e.g. time-reversal and sublattice symmetry.
The other, more sensible, approach, would be to get the operators $\hat{U}$ into the braket, so $$\left<\hat\psi_i\hat\psi_j^\dagger\right>\rightarrow \hat{U}\left<\hat\psi_i\hat\psi_j^\dagger\right>\hat{U}^{-1}=\left<\hat{U}\hat\psi_i\hat\psi_j^\dagger\hat{U}^{-1}\right>=\left<\hat{U}\hat\psi_i\hat{U}^{-1}\hat{U}\hat\psi_j^\dagger\hat{U}^{-1}\right>.$$ This looks better, but still does not satisfy me, because with these approach the state wouldn't be affected by whether the operator $\hat{U}$ is unitary or anti-unitary, in the sense that the anti-unitarity in this calculation doesn't complex conjugate anything.
So in conclusion, neither of this approaches satisfy me, as the first one is insensitive to the specific properties of the transformation, while the second one is insensitive to the unitary or anti-unitary character of it. Is any of this approaches correct, or maybe some intermediate approach?