Which classes of (scalar or systems of) linear first or second-order order hyperbolic PDEs $Lf=g$ in $n$ variables and spaces $S$ of functions have the property that there is a retarded Green's function $G:S\to S$ with $LG=1$. Retarded means that for every $g\in S$ whose support is disjoint from the past causal cone of $x$, the support of the retarded solution $g=Gf$ of $Lf=g$ is disjoint from the past causal cone of $x$.
When can one choose $S$ as a space of smooth (i.e., $C^\infty$) functions?
Where can I read about the mathematical tools for studying this and similar questions?
A particular case I am interested in is hyperbolic $L$ of the form $$L=\pmatrix{\alpha & (a-d)^* \cr b-d & B},$$ where $d$ is the exterior derivative, ${}^*$ is the Minkowski adjoint, $\alpha$ is a scalar field on Minkowski space, $a,b$ are covector fields, and $B$ is a matrix field mapping covectors to covectors. Hyperbolicity is defined as in Wikipedia.