Every time I am faced with an analysis of the spectral function it looks like a "new" unintuitive set of jugglery with expectation values and I am unable to see a general picture of what this construction means. I am not sure that I can frame a coherent single question and hence I shall try to put down a set of questions that I have about the idea of spectral functions in QFT.
I guess in a Dirac field theory one defines the spectral function as follows, $$ \rho_{ab}(x-y) = \frac{1}{2} \langle0|\{\psi_a(x),\bar{\psi}_b(y)\}|0\rangle. $$ Also I see this other definition in the momentum space as, $$ S_{Fab}\,(p) = \int _0 ^\infty d\mu^2 \frac{\rho_{ab}(\mu^2)}{p^2 - \mu ^2 + i\epsilon} $$ Are these two the same things conceptually? I tried but couldn't prove an equivalence. (I define the Feynman propagator as $S_{Fab} = \langle0|T\, \psi_a(x) \, \bar{\psi}_b(y)|0\rangle$.)
Much of the algebraic complication I see is in being able to
handle the quirky minus sign in the time-ordering of the fermionic
fields which is not there in the definition of the Feynman propagator of the Klein-Gordon field (..which apparently is seen by all
theories!..) and to see how the expectation values that one gets like $\langle0|\psi_a(0)|n\rangle \langle n|\bar{\psi}_b(0)|0\rangle$ and $\langle0|\bar{\psi}_b(0)|n\rangle \langle n|\psi_a(0)|0\rangle$ and how these are in anyway related to the Dirac operator $(i\gamma^\mu p_\mu +m)_{ab}$ that will come-up in the far more easily doable
calculation of the spectral function for the free Dirac theory.Is it true that for any QFT given its Feynman propagator $S_F(p)$ there will have to exist a positive definite function $\rho(p^2)$ such that the relation, $$ S_F(p) = \int _0^\infty d\mu ^2 \frac{\rho(\mu^2)} {p^2 - \mu^2 +i\epsilon} $$ is satisfied?
So no matter how complicatedly interacting a theory for whatever spin it is, its Feynman propagator will always "see" the Feynman propagator for the Klein-Gordon field at some level? (..all the interaction and spin intricacy being seen by the spectral function weighting it?..)
One seems to say that it is always possible to split the above integral into two parts heuristically as,
$$\begin{eqnarray}S_F(p) &=& \sum (\text{free propagators for the bound states})\\ &&+ \int_\text{states} \big( (\text{Feynman propagator of the Klein-Gordon field of a certain mass})\\ && ~~~~~~~~~~~~~~\times(\text{a spectral function at that mass})\big)\end{eqnarray}$$
Is this splitting guaranteed irrespective of whether one makes the usual assumption of "adiabatic continuity" as in the LSZ formalism or in scattering theory that there is a bijection between the asymptotic states and the states of the interacting theory - which naively would have seemed to ruled out all bound states?
To put it another way - does the spectral function see the bound states irrespective of or despite the assumption of adiabatic continuity?