Are there any results about regularity of the generalized eigenfunctions used in spectral representations of unbounded differential operators?

I am using the "Direct integral" version of the spectral theorem, as given e.g. in Hall. It says that a diagonal representation for an unbounded operator on the Hilbert space $L^2(mathbb{R}^n)$ can be constructed as a "direct integral", or

$$
mathcal{H}_{hat{A}} = int^oplus_{sigma(hat{A})}dmu(a) mathcal{H}_a.
$$

States in the Hilbert space are formed as square-integrable "sections", or as functions of $a$ with values in the Hilbert spaces $mathcal{H}_a$. Since the $mathcal{H}_a$ spaces are not true subspaces, vectors in them are not actually in $mathcal{H}$, so are a "generalised eigenvectors". They can be defined as a functional on a dense subset of the Hilbert space by means of

$$
v in mathcal{H}_a: ; psi(a^prime) ; mapsto (psi(a),v) in mathbb{C},
$$

which can then of course be extended back to functionals on $L^2(mathbb{R}^n)$. I’m fairly sure that this chapter’s results can be used to demonstrate that if $hat{A}$ is in the algebra generated by $hat{x}$ and $hat{p}$ then these functionals will necessarily be tempered distributions.

Now, the specific $hat{A}$ I’m interested in is a differential operator on $L^2(mathbb{R}^{3oplus1})$. (Specifically, it’s the "Klein-Gordon-like" operator used for spin-half particles in an applied electromagnetic field, formed by squaring the Dirac operator.) What I would really like to be able to do is, given that $hat{A}$ is a differential operator, use some results about PDEs on these generalised eigenvectors.

This, though, requires some notions about the regularity of the generalized eigenfunctions, which I have just not been able to find. My issue is that the operator $hat{A}$ I’m using is second-order Hyperbolic, and so there aren’t many good results about regularity in the PDE maths literature for me. This seems to be because the mere fact that $psi$ obeys the differential equation $(hat{A}-a)psi=0$ isn’t enough to enforce any sort of regularity, for a similar reason that $delta(x-ct)$ is a solution to the wave equation. I think that I therefore need regularity enforced by the use of the solution to the differential equation as a generalized eigenfunction, or as a basis state in the direct integral spectral representation.

I realise this might be a bit of a long shot, but I’ve struggled for days to find such a result, and it would be extremely useful for me. Does anyone know of any literature which deals with this question of the regularity of the basis of spectral representations of differential operators?