{\rm par}{{\rm par}}
\def{\rm sch}{{\rm sch}}
\def{\rm Im}{{\rm Im}}
\def{\rm Re}{{\rm Re}}
\def{\rm s-}\lim{{\rm s-}\lim}
\def{\rm w-}\lim{{\rm w-}\lim}
\def{\rm s-}C_\infty{\rm -}\lim{{\rm s-}C_\infty{\rm -}\lim}
\def{\rm cl}{{\rm cl}}
\def{\rm Bog}{{\rm Bog}}
\def{\rm H}{{\rm H}}
\def\scriptstyle{/}{\scriptstyle{/}}
\def{\rm cpl}{{\rm cpl}}
\def{\rm Cl}{{\rm Cl}}
\def\overline{\overline}
\defC_\infty{\rm -}\lim{C_\infty{\rm -}\lim}
\def{\rm nt}{{\rm nt}}
\def{\rm reg}{{\rm reg}}
\def{\rm sg}{{\rm sg}}
\def\frac{d}{dt}{\frac{d}{dt}}
\def{\rm pp}{{\rm pp}}
\def{\rm sep}{{\rm sep}}
\def{\rm gen}{{\rm gen}}
\def{\rm Pf}{{\rm Pf}}
\def{\rm B}{}
\def\bar{z}{\overline{z}}
\defl{l}
\def{\rm disc}{{\rm disc}}
\def{\rm alg}{{\rm alg}}
\def{\rm r}{{\rm r}}
\defC^\infty{C^\infty}
\defC_0^\infty{C_0^\infty}
\def{\rm q}{{\rm q}}
\defs{s}
\def{\rm s/a}{{\rm s/a}}
\def{\rm graph}{{\rm graph}}
\def{\rm ren}{{\rm ren}}
\def{\rm rank}{{\rm rank}}
\def{\rm sp}{{\rm sp}}
\def{\rm span}{{\rm span}}
\def\proof{
\noindent{\bf Proof.}\ \ }
\def{\rm imp}{{\rm imp}}
\def\skew4\hat f{\skew4\hat f}
\def{\rm al}{{\rm al}}
\def\skew4\tilde f{\skew4\tilde f}
\def{\rm h}{{\rm h}}
\def{\mathfrak h}{{\mathfrak h}}
\def{\rm sh}{{\rm sh}}
\def{\rm th}{{\rm th}}
\def{\rm tg}{{\rm tg}}
\def{\rm ess\:inf}{{\rm ess\:inf}}
\def{\rm ess\:sup}{{\rm ess\:sup}}
\def{\rm Pol}{{\rm Pol}}
\def{\rm pol}{{\rm pol}}
\def{\rm FBI}{{\rm FBI}}
\def{\odot}{{\odot}}
\def{\rm an}{{\rm an}}
\def{a\in{\mathcal A}}{{a\in{\mathcal A}}}
\def{\mathcal Z}{{\mathcal Z}}
\def{\mathcal Y}{{\mathcal Y}}
\def{\mathcal V}{{\mathcal V}}
\def{\mathcal L}{{\mathcal L}}
\def{\mathcal S}{{\mathcal S}}
\def{\mathcal Q}{{\mathcal Q}}
\def{\mathcal D}{{\mathcal D}}
\def{\mathcal U}{{\mathcal U}}
\def{\mathcal M}{{\mathcal M}}
\def{\mathcal N}{{\mathcal N}}
\def{\mathcal O}{{\mathcal O}}
\def{\mathcal C}{{\mathcal C}}
\def{\mathcal W}{{\mathcal W}}
\def{\mathfrak F}{{\mathfrak F}}
\def\fM{{\mathfrak M}}
\def{\rm cw}{{\rm cw}}
\def{\rm fin}{{\rm fin}}
\def{\mathfrak f}{{\mathfrak f}}
\def{a_{\rm min}}{{a_{\rm min}}}
\def{a_{\rm max}}{{a_{\rm max}}}
\def{\rm ham}{{\rm ham}}
\def\frac{x}{c}{\frac{x}{c}}
\def{\rm CCR}{{\rm CCR}}
\def{\rm CAR}{{\rm CAR}}
\def{\rm Cliff}{{\rm Cliff}}
\def{\rm Pair}{{\rm Pair}}
\def{\rm a}{{\rm a}}
\def{\rm cm}{{\rm cm}}
\def{\rm red}{{\rm red}}
\def{\rm inert}{{\rm inert}}
\def{\rm free}{{\rm free}}
\def{\rm mod}{{\rm mod}}
\def{\rm vec}{{\rm vec}}
\def{\rm i}{{\rm i}}
\def{\rm ad}{{\rm ad}}
\def{\rm Ad}{{\rm Ad}}
\def{\rm aut}{{\rm aut}}
\def{\rm der}{{\rm der}}
\def{\rm Span}{{\rm Span}}
\def{\rm Dom}{{\rm Dom}}
\def{\mathop{Mp}\limits^\circ}{{\mathop{Mp}\limits^\circ}}
\def{\mathop{Gl}\limits^\circ}{{\mathop{Gl}\limits^\circ}}
\def{\mathop{gl}\limits^\circ}{{\mathop{gl}\limits^\circ}}
\def{\rm Aut}{{\rm Aut}}
\def{\mathop{Sl}\limits^\circ}{{\mathop{Sl}\limits^\circ}}
\def{\mathop{sl}\limits^\circ}{{\mathop{sl}\limits^\circ}}
\def{\mathop{L}\limits^\circ}{{\mathop{L}\limits^\circ}}
\def{\mathop{O}\limits^\circ}{{\mathop{O}\limits^\circ}}
\def{\mathop{SO}\limits^\circ}{{\mathop{SO}\limits^\circ}}
\def{\mathop{o}\limits^\circ}{{\mathop{o}\limits^\circ}}
\def{\rm KN}{{\rm KN}}
\def{\rm aKN}{{\rm aKN}}
\def{|\!\times}{{|\!\times}}
\def{\rm KN}{{\rm KN}}
\def{\rm aKN}{{\rm aKN}}
\def{\rm cv}{{\rm cv}}
\def{\rm ct}{{\rm ct}}
\def{\rm Fock}{{\rm Fock}}
\def\Omega{\Omega}
\def{\rm loc}{{\rm loc}}
\def\mathop{\to}\limits{\mathop{\to}\limits}
\def\mathop{\otimes}\limits{\mathop{\otimes}\limits}
\def\mathop{\cup}\limits{\mathop{\cup}\limits}
\def\mathop{\cap}\limits{\mathop{\cap}\limits}
\def\mathop{\times}\limits{\mathop{\times}\limits}
\def\mathop{\oplus}\limits{\mathop{\oplus}\limits}
\def\mathop{\bigsqcup}\limits{\mathop{\bigsqcup}\limits}
\def\mathop{\Pi}\limits{\mathop{\Pi}\limits}
\def{\rm sgn}{{\rm sgn}}
\def{\rm id}{{\rm id}}
\def{\rm Tr}{{\rm Tr}}
\def{\rm T}{{\rm T}}
\def{\rm Or}{{\rm Or}}
\def{\rm Ker}{{\rm Ker}}
\def{\rm Gr}{{\rm Gras}}
\def{\rm c}{{\rm c}}

\def$\Box$\medskip{$\Box$\medskip}
\def{\rm A}{{\rm A}}

\def\tilde{\eta}{\tilde{\eta}}
\def\tilde{y}{\tilde{y}}
\def\tilde{z}{\tilde{z}}
\def\tilde{U}{\ttilde{U}}
\def\ttilde{P}{\ttilde{P}}
\def\tilde{J}{\tilde{J}}
\def{\tilde{D}}{{\tilde{D}}}
\def\ttilde\Xi{\ttilde\Xi}
\def\rm Int{\rm Int}
\def  \partial{ \partial}
\def{\mathop{sp}\limits^\circ}{{\mathop{sp}\limits^\circ}}
\def{\mathop{Sp}\limits^\circ}{{\mathop{Sp}\limits^\circ}}
\def{\mathop{U}\limits^\circ}{{\mathop{U}\limits^\circ}}
\def{\mathop{u}\limits^\circ}{{\mathop{u}\limits^\circ}}
\def\langle u \rangle{\langle u \rangle}
\def{\mathcal Sp}{{\mathcal Sp}}
\def\frac{1}{2}{\frac{1}{2}}
\def\frac{1}{4}{\frac{1}{4}}
\def\langle x \rangle{\langle x \rangle}
\def\frac{x}{t}{\frac{x}{t}}
\def{\rm supp}{{\rm supp}}
\def{\bf D}{{\bf D}}
\def{\rm hc}{{\rm hc}}
\def{\rm ad}{{\rm ad}}
\def{\rm e}{{\rm e}}
\def\mathcal D{\mathcal D}
\def{\tilde Y}{{\tilde Y}}
\def{\rm d}{{\rm d}}
\def{\rm Ran}{{\rm Ran}}
\def\tilde{S}{\tilde{S}}
\def{\rm ex}{{\rm ex}}
\def{\mathchoice {\rm 1\mskip-4mu l} {\rm 1\mskip-4mu l{{\mathchoice {\rm 1\mskip-4mu l} {\rm 1\mskip-4mu l}
{\rm 1\mskip-4.5mu l} {\rm 1\mskip-5mu l}}}
\def\bbbone{{\mathchoice {\rm 1\mskip-4mu l} {\rm 1\mskip-4mu l}
\def{\mathcal H}{{\mathcal H}}
\def{\rm as}{{\rm as}}
\def{\rm tq}{{\rm tq}}
\def{\rm can}{{\rm can}}
\def{\rm j}{{\rm j}}
\def {\rm j}_1{ {\rm j}_1}
\defl{l}
\def{\rm r}{{\rm r}}
\def{\rm w}{{\rm w}}
\def{\widetilde{Sp}}{{\tilde{Sp}}}
\def{\tilde X}{{\tilde X}}
\def{\tilde\Lambda}{{\tilde\Lambda}}
\def{\rm cm}{{\rm cm}}
\def{\rm comp}{{\rm comp}}
\def{\rm dual}{{\rm dual}}
\def\mathop{\otimes}\limits{\mathop{\otimes}\limits}
\def\mathop{\times}\limits{\mathop{\times}\limits}
\def\mathop{\oplus}\limits{\mathop{\oplus}\limits}
\def\mathop{\bigsqcup}\limits{\mathop{\bigsqcup}\limits}
\def\mathop{\Pi}\limits{\mathop{\Pi}\limits}
\def{\rm sgn}{{\rm sgn}}
\def{\rm id}{{\rm id}}
\def{\rm Tr}{{\rm Tr}}
\def{\rm Or}{{\rm Or}}
\def{\rm Ker}{{\rm Ker}}
\def{\rm Gr}{{\rm Gr}}
\def{\rm c}{{\rm c}}
\def{\rm nondeg}{{\rm nondeg}}
\def{\rm Mas}{{\rm Mas}}
\defC_0^\infty{C_0^\infty}
\def{\mathcal P}{{\mathcal P}}
\def{\mathcal J}}\def\C{{\mathcal C}{{\mathcal J}}\def\C{{\mathcal C}}
\def{\mathcal T}{{\mathcal T}}
\def{\mathcal F}{{\mathcal F}}
\def{\mathcal G}{{\mathcal G}}
\def{\mathcal X}{{\mathcal X}}
\def{\mathcal I}{{\mathcal I}}
\def{\mathcal N}{{\mathcal N}}
\def{\mathcal S}{{\mathcal S}}
\def{\mathcal Q}{{\mathcal Q}}
\def{\mathcal X}{{\mathcal X}}
\def{\mathcal K}{{\mathcal K}}
\def{\langle R\rangle}{{\langle R\rangle}}
\def{\rm rw}{{\rm rw}}
\def{\rm cw}{{\rm cw}}
\def{\rm rw1}{{\rm rw1}}
\def{\rm cw1}{{\rm cw1}}

\def\tilde{y}{\tilde{y}}
\def\tilde{z}{\tilde{z}}
\def\tilde{U}{\ttilde{U}}
\def\ttilde{P}{\ttilde{P}}
\def\tilde{J}{\tilde{J}}
\def{\tilde{D}}{{\tilde{D}}}
\def\ttilde\Gamma{\ttilde\Gamma}
\def\ttilde\gamma{\ttilde\gamma}
\def\ttilde\Xi{\ttilde\Xi}
\def\ttilde{\dd}{\ttilde{{\bf D}}}
\def\rm Int{\rm Int}
\def  \partial{ \partial}
\def\langle u \rangle{\langle u \rangle}
\def{\mathcal X}{{\mathcal X}}
\def\frac{1}{2}{\frac{1}{2}}
\def\langle x \rangle{\langle x \rangle}
\def\frac{x}{t}{\frac{x}{t}}
\def{\rm supp}{{\rm supp}}
\def{\bf D}{{\bf D}}
\def{\rm hc}{{\rm hc}}
\def{\rm ad}{{\rm ad}}
\def{\rm e}{{\rm e}}
\def{\tilde Y}{{\tilde Y}}
\def{\rm d}{{\rm d}}
\def{\rm Ran}{{\rm Ran}}
\def\tilde{S}{\tilde{S}}
\def{\rm ex}{{\rm ex}}

\def{\mathcal H}{{\mathcal H}}
\def{\rm as}{{\rm as}}
\def{\rm tq}{{\rm tq}}
\def{\rm tame}{{\rm tame}}
\def{\rm ie}{{\rm ie}}

\def{\rm Diff}{{\rm Diff}}
\def{\mathcal W}{{\mathcal W}}
\def\begin{proposition}{\begin{proposition}}
\def\end{proposition}{\end{proposition}}
\def{\rm Op}^{a^*\!,a}{{\rm Op}^{a^*\!,a}}
\def{\rm Op}^{a,a^*}{{\rm Op}^{a,a^*}}
\def{\rm KN}{{\rm KN}}
\def{\rm aKN}{{\rm aKN}}
\def{\rm b}{{\rm b}}
\def{\rm rev}{{\rm rev}}
\def{\widetilde{Sp}}{{\widetilde{Sp}}}
\def{\rm w}{{\rm w}}
\def{\widetilde{Sp}}{{\widetilde{Sp}}}
\def\tilde{U}{\tilde{U}}
\def{\rm res}{{\rm res}}

\def{\rm Op}^{\rm w}{{\rm Op}^{\rm w}}
\def{\rm Op}^{1, 0}{{\rm Op}^{1, 0}}
\def{\rm Op}{{\rm Op}}
\defs{s}
\def{\scriptscriptstyle\#}{{\scriptscriptstyle\#}}
\newcommand{\mat}[4]{\left(\begin{array}{cc}#1 &#2  \\ #3 &#4 \end{array}\right)}
\newcommand{\lin}[2]{\left(\begin{array}{c}#1 \\#2\end{array}\right)}
\def{\rm ext}{{\rm ext}}

\def\CARal{{\rm C\hskip 0.25 em \hbox{\raise 1.72 ex
\hbox{$\scriptscriptstyle\rm al$}\kern -0.57 em A}R}}

\def{\scriptscriptstyle\#}{{\scriptscriptstyle\#}}
\def\otimesal{\mathop{\hbox{\raise 1.5 ex
  \hbox{$\scriptscriptstyle\rm al$}
\kern -0.92 em \hbox{$\otimes$}}}}
\def\oplusal{\mathop{\hbox{\raise 1.5 ex
  \hbox{$\scriptscriptstyle\rm al$}
\kern -0.92 em \hbox{$\oplus$}}}}
\def\Gammal{\hbox{\raise 1.68 ex
\hbox{$\scriptscriptstyle\rm al$}\kern -0.50 em $\Gamma$}}
\def\Bal{\hbox{\raise 1.68 ex
\hbox{$\scriptscriptstyle\rm  al$}\kern -0.50 em $B$}}
\def\CARal{{\rm C\hskip 0.25 em \hbox{\raise 1.72 ex
\hbox{$\scriptscriptstyle\rm al$}\kern -0.57 em A}R}}
\def{\scriptscriptstyle\#}{{\scriptscriptstyle\#}}

\def{\mathcal E}{{\mathcal E}}

\def{\rm x}{{\rm x}}
\def\rA{{\rm A}}\def
\frac{1}{4}{\frac{1}{4}}
\newcommand\Texp{{\rm  Texp}}
\def{\rm WF}{{\rm WF}}
\newcommand{\co}[1]{T^{*}\! #1}
\newcommand{\coo}[1]{T^{*}\! #1\backslash Z}
\newcommand{\cooo}[1]{T^{*} #1\backslash \{0\}}

\newcommand{{\textstyle \frac{1}{2}}}{{\textstyle \frac{1}{2}}}
\newcommand{{\frac{1}{2}}}{{\frac{1}{2}}}
\newcommand{{\rm (KG)}}{{\rm (KG)}}
\newcommand{( \mu {\rm sc})}{( \mu {\rm sc})}
\newcommand{{\rm mod} \ C^{\infty}}{{\rm mod} \ C^{\infty}}
\newcommand{\bra}{\langle}
\newcommand{\rangle}{\rangle}
\newenvironment{commw}{\hspace{5pt}\itshape !` }{ --- M.W. ! \hspace{5pt}}
\def(\rr^{d})\otimes \cc^{2}{({\mathbb R}^{d})\otimes {\mathbb C}^{2}}

\begin{document}
\title[Construction of Hadamard states]{Construction of Hadamard states by pseudo-differential calculus}
\author{C. G\'erard}
\address{D\'epartement de Math\'ematiques, Universit\'e de Paris XI, 91405 Orsay Cedex France}
\email{christian.gerard@math.u-psud.fr}
\author{M. Wrochna}
\address{Mathematisches Institut, Universit\"at G\"ottingen, Bunsenstr. 3-5, D - 37073 G\"ottingen, Germany}
\email{wrochna@uni-math.gwdg.de}
\keywords{Hadamard states, microlocal spectrum condition, pseudo-differential calculus}
\subjclass[2010]{81T20, 35S05, 35S35}
\begin{abstract}
 We give a new construction based on pseudo-differential calculus of quasi-free Hadamard states for Klein-Gordon equations on a class of space-times whose metric is well-behaved at spatial infinity. In particular on this class of space-times,  we construct all pure Hadamard states whose two-point function (expressed in terms of Cauchy data on a Cauchy surface) is a matrix of pseudo-differential operators.  We also study their covariance under symplectic transformations.

 As an aside, we give a new construction of Hadamard states on arbitrary globally hyperbolic space-times which is an alternative to the classical construction by Fulling, Narcowich and Wald.
\end{abstract}
\maketitle

\section{Introduction}\label{sec:introduction}
\subsection{Hadamard states}
Hadamard states are nowadays widely accepted as possible physical states of the non-interacting quantum field theory on a curved space-time. One of the main reason is their applicability to renormalization of the stress-energy tensor, a necessary step in the formulation of semi-classical Einstein equations. Moreover, the Hadamard condition plays an essential role in the perturbative construction of interacting quantum field theory \cite{BF00}. Other related concepts making use of Hadamard states include local thermal equilibrium \cite{schlemmer} and quantum energy inequalities \cite{fewster}.

Since the work of Radzikowski \cite{radzikowski}, the Hadamard condition (renamed \emph{microlocal spectrum condition}), is formulated as a requirement for the wave front set of the associated two-point function $\Lambda$, which is necessarily a bi-solution of the free equations of motion. It is therefore natural to try to construct such states using the standard apparatus of microlocal analysis, based on pseudo-differential calculus. Although a construction is already known for space-times with compact Cauchy surface \cite{junker}, it does certainly not cover many cases of physical interest and lacks the capability to produce many states on a fixed space-time with distinct properties.

 In this paper we address these questions for  a class of space-times whose metric components are suitably well-behaved at spatial infinity, allowing also for external potentials. In this case  it is possible to obtain rather complete and transparent results.
 Namely we construct a large class of quasi-free Hadamard states, whose two-point functions, expressed in terms of Cauchy data on a fixed Cauchy surface, are matrices of pseudo-differential operators.
 In particular we can construct \emph{all } such \emph{pure} Hadamard states and study their covariance under symplectic transformations.

 As an additional result we give a new construction of Hadamard states on arbitrary globally hyperbolic space-times which is an alternative to the classical construction by Fulling, Narcowich and Wald \cite{FNW}. Our method turns out not to produce pure states in general, it allows however to keep track of their local properties.

\subsection{Methods}
Our analysis is set on three levels:
\begin{enumerate}\setlength{\leftmargin}{5ex}\setlength{\itemindent}{0cm}
\item\label{it:i1} Our starting point are normally hyperbolic operators on ${\mathbb R}\times{\mathbb R}^d$ of the form
\begin{equation}\label{asto}
 \partial_{t}^{2}+ a(t, {\rm x}, D_{{\rm x}})= \partial_t^2 - \partial_{{\rm x}^{j}}a^{jk} \partial_{{\rm x}^{k}}+  b^{j} \partial_{{\rm x}^{j}}-  \partial_{{\rm x}^{j}}\overline{b}^{j}+ m,
\end{equation}
where
\begin{equation}\label{asto-hyp}
\begin{array}{l}
a^{jk}, b^{j}, m \in C^{\infty}({\mathbb R}, C^{\infty}_{\rm bd}({\mathbb R}^{d})), \   m(x)\in {\mathbb R},\\[2mm]
[a^{jk}](x)\geq  c(t)\bbbone\hbox{ \  uniformly on }{\mathbb R}^{1+d}, \ c(t)>0.
\end{array}.
\end{equation}
We refer to this case as the \emph{model Klein--Gordon equation} and give a construction of the associated parametrix  for  the Cauchy problem, in such way that the propagation of positive-frequency and negative-frequency singularities is under control. This allows us to reformulate the microlocal spectrum condition in terms of Cauchy data. We show how to construct many non-necessarily pure Hadamard states and then characterize pure ones. We also describe classes of symplectic transformations which preserve the microlocal spectrum condition.

\item\label{it:i2} The above results are easily extended to operators of the form $f ( \partial_{t}^{2}+ a(t) )g$, where $f$ and $g$ are smooth densities. This way, we show that the problem of constructing Hadamard states is reduced to the model case above if $M={\mathbb R}\times{\mathbb R}^d$, the metric is given by
\begin{equation}\label{asti}
g= -c (x)dt^{2}+ h_{jk}(x)d{\rm x}^{j}d{\rm x}^{k},
\end{equation}
and the Klein-Gordon operator is of the form
\[\begin{array}{rl}
P(x, D_{x})=&\ c^{-\frac{1}{2}}|h|^{-\frac{1}{2}}( \partial_{t}+ {\rm i} V) c^{-\frac{1}{2}}|h|^{\frac{1}{2}}( \partial_{t}+{\rm i} V)\\[2mm]
&- c^{-\frac{1}{2}}|h|^{-\frac{1}{2}}( \partial_{j}+{\rm i} A_{j})c^{\frac{1}{2}}| h|^{\frac{1}{2}}h^{jk}( \partial_{k}+{\rm i} A_{k})+ \rho,
\end{array}
\]
where
$A_{\mu}(x)= (V(x), \rA_{j}(x))$, $|h|= \det[h_{jk}]$, $ [h^{jk}]= [h_{jk}]^{-1}$ and the following hypotheses are assumed:
\begin{equation}\label{asti-hyp}
 \begin{array}{rl}
&\forall \ I\subset {\mathbb R}\hbox{ compact interval }\exists\  C>0\hbox{ such that }\\[2mm]
&C\leq c(x), \ C\bbbone\leq [h_{jk}(x)],\hbox{uniformly for }x\in I\times {\mathbb R}^{d},\\[2mm]
&  h_{jk}(x), \ c(x), \ \rho(x), A_{\mu}(x)\in C^{\infty}({\mathbb R}, C^{\infty}_{\rm bd}({\mathbb R}^{d})).
\end{array}
\end{equation}
\item\label{it:i3} For arbitrary space-times (and external potentials), using a suitable partition of unity, we explain how to glue together two-point functions of Hadamard states on smaller regions of the space-time into a globally-defined one. Using the results obtained for the special case above, this yields a new construction of Hadamard states on arbitrary globally-hyperbolic space-times.
\end{enumerate}

Let us mention that beside the construction for space-times with compact Cauchy surface due to Junker, a general existence result for Hadamard states is known \cite{FNW}, as well as a collection of various examples on specific classes of space-times. However, the existence argument of Fulling, Narcowich and Wald has the disadvantage of being highly non-explicit and requires non-local information on the space-time as an input. Those drawbacks are to a large extent avoided in our approach, as explained in Section \ref{sec7}. As for the known examples of Hadamard states for the Klein-Gordon equation, these include:
\begin{romanenumerate}
\item passive states for stationary space-times (this includes ground- and KMS states) \cite{passivity},
\item states constructed in \cite{DMPas} for a subclass of asymptotically flat vacuum space-times at null infinity (see \cite{moretti} for the proof),
\item states constructed in \cite{DMPcosm} for a class of cosmological space-times (this includes the Bunch-Davies state on de-Sitter space-time),
\item so-called states of low energy for FLRW space-times \cite{olbermann},
\item the so-called Unruh state \cite{unruh},
\item ground states and over-critical states for static potentials on Minkowski space-time considered in \cite{wrochna}.
\end{romanenumerate}
A short inquiry shows that the sets of assumptions (\ref{asto-hyp}) and (\ref{asti-hyp}) studied by us in greater detail are only partially covered by the examples above.

\subsection{Plan of the paper}

The paper is organized as follows.

In Section \ref{sec0} we recall basic facts on bosonic quasi-free states. A special emphasis is put on explaining the relation between the neutral and the (less often discussed) {\em charged} case. In our problem, the use of microlocal analysis makes much more natural to work with complex quantities. In order to cover both cases it is sufficient to consider gauge-invariant charged quasi-free states.

Section \ref{sec22} contains basic definitions and facts on the Klein-Gordon equation, wave front sets and Hadamard states.

 In Section \ref{sec1} we recall mostly well-known results on pseudo-differential calculus needed later on. This includes theorems on the pseudo-differential property of functions of pseudo-differential  operators and several results related to Egorov's theorem. In Section \ref{sec2} we specify our assumptions for the space-times (\ref{asti}) and explain the reduction to the model Klein-Gordon equation (\ref{asto}).

Section \ref{sec3} contains the key technical ingredient used to construct Hadamard states in Section \ref{sec4},  namely a sufficiently explicit construction of a {\em parametrix} for the Cauchy problem.  Let us explain in  more details  this construction and its use.

As is well-known, the appropriate phase-space for the quantization of the Klein-Gordon equation is the space of its smooth {\em space-compact solutions}, which can be identified with  the image of ${\mathcal D}(M)$  under the Pauli-Jordan commutator function $E$.
Alternatively one can fix  a Cauchy surface $\Sigma\subset M$ and identify a space-compact solution $\phi$ with its Cauchy data $\rho\phi= (\rho_{0}\phi, \rho_{1}\phi)\in {\mathcal D}(\Sigma)\oplus {\mathcal D}(\Sigma)$.

The Hadamard condition on a state $\omega$ is formulated in terms of its  two-point function $\Lambda$, which is an element of ${\mathcal D}'(M\times M)$.  For practical purposes it is more convenient to fix a Cauchy surface $\Sigma\subset M$ and to consider instead  the two-point function $\lambda$ of $\omega$ in terms of Cauchy data. Both two-point functions are related by:
\begin{equation}\label{covcov}
 \Lambda= (\rho\circ E)^{*} \circ \lambda \circ (\rho\circ E).
\end{equation}
Since the Hadamard condition singles out the {\em positive energy} component  ${\mathcal N}_{+}$ of the characteristic manifold ${\mathcal N}= \{(x, \xi)\in T^{*}M: \xi_{\nu}g^{\mu\nu}(x)\xi_{\mu}=0\}$, we see from (\ref{covcov}) that the essential step is to formulate conditions on the Cauchy data $\rho \phi$ which ensure that a (distributional or finite energy) solution $\phi$ will have its wave front set contained in ${\mathcal N}_{+}$.
The natural way to do this is to construct a {\em parametrix} for the Cauchy problem on $\Sigma$, i.e. a solution  modulo smooth errors of the Cauchy problem with data $(f_{0}, f_{1})$ compactly supported distributions.

Such a construction is well-known in microlocal analysis and is usually done with  {\em Fourier integral operators}.   To obtain more transparent results later on,  we use instead  a more abstract construction, relying on Hilbert space methods, which leads to a smaller number of arbitrary choices. In our case the arbitrariness is described by the choice of a pseudo-differential operator $r$ acting on $\Sigma$ (see Def. \ref{defo}).  A concrete consequence of our parametrix construction is as follows:  let us consider for definiteness the space ${\rm Sol}_{E}(P)$ of {\em finite energy solutions} of the Klein-Gordon equation(see Subsect. \ref{sec3.3}). Then for each choice of $r$ this space splits into the  direct sum of two spaces ${\rm Sol}_{E}^{\pm}(P, r)$, whose elements have wavefront sets in ${\mathcal N}_{\pm}$ and which are {\em symplectically orthogonal}.  We will use this fact through an appropriate reparametrization of Cauchy data in which the symplectic form still has a very simple form (see Lemma \ref{4.1}).

Section \ref{sec4} contains the key results of the paper. Theorem \ref{4.2} characterizes Hadamard states for the model Klein--Gordon equation in terms of their symplectically smeared two-point function $\lambda$. This allows us to construct a large class of Hadamard states in Subsect. \ref{4.3}. In Subsect. \ref{sec4.4} we describe classes of symplectic transformations which preserve the microlocal spectrum condition. Stronger results are derived for pure quasi-free states in Subsect. \ref{sec4.5}, in particular {\em all} pure Hadamard states with pseudo-differential two-point functions are characterized.

Among these pure states, there exists  a `canonical' Hadamard state  with two-point function $\lambda(r)$, which is distinguished modulo the choice of the pseudo-differential operator $r$ appearing in the construction of the parametrix. We explicitely find the symplectic transformation relating two canonical two-point functions  $\lambda(r_{1})$, $\lambda(r_{2})$ for two different choices of $r$. In this way the remaining ambiguities in the construction of $\lambda(r)$ are completely understood.

In Subsect. \ref{for-record} we briefly discuss the static case and show how the ground state and KMS states fit in our construction.

In Section \ref{sec7} we present our  alternative construction of Hadamard states on an arbitrary globally hyperbolic space-time. Instead of the classical deformation argument  of \cite{FNW}, we fix a Cauchy surface $\Sigma$ and use the fact that the symplectic form expressed in terms of Cauchy data is a local operator.

We include some remarks on the case of a compact Cauchy surface and compare our results with the construction from \cite{junker} in Subsect. \ref{ssec:junker}. Various proofs are collected in Appendix \ref{sec-app1}.

\section{Bosonic quasi-free states}\label{sec0}\setcounter{equation}{0}

In this section we recall well-known facts about bosonic quasi-free states, following \cite {DG}.  We choose to  work with complex symplectic spaces and gauge invariant states  (or equivalently charged fields), possibly after complexification of a  real situation. In this respect the structures obtained by complexifying a complex symplectic space deserves some attention, see  Subsect. \ref{complix}.

\subsection{Notation}\label{notation}
\def\!\cdot\!{\!\cdot\!}
If ${\mathcal X}$ is a real or complex vector space we denote by ${\mathcal X}^{{\scriptscriptstyle\#}}$ its dual. Bilinear forms on ${\mathcal X}$ are identified with elements of $L({\mathcal X}, {\mathcal X}^{{\scriptscriptstyle\#}})$, which leads to the notation $x_{1}\!\cdot \!bx_{2}$ for $b\in L({\mathcal X}, {\mathcal X}^{{\scriptscriptstyle\#}})$, $x_{1}, x_{2}\in {\mathcal X}$. The space of symmetric (resp. anti-symmetric) bilinear forms on ${\mathcal X}$ is denoted by $L_{\rm s}({\mathcal X}, {\mathcal X}^{{\scriptscriptstyle\#}})$ (resp. $L_{\rm a}({\mathcal X}, {\mathcal X}^{{\scriptscriptstyle\#}})$).

If $\sigma\in L_{\rm s}({\mathcal X}, {\mathcal X}^{{\scriptscriptstyle\#}})$, we denote by $O({\mathcal X}, \sigma)$ the linear (pseudo-)orthogonal group on ${\mathcal X}$. Similarly if
 $\sigma\in L_{\rm a}({\mathcal X},{\mathcal X}^{{\scriptscriptstyle\#}})$ is non-degenerate, i.e. $({\mathcal X}, \sigma)$ is a symplectic space, we denote by $Sp({\mathcal X}, \sigma)$ the linear symplectic group on ${\mathcal X}$.

If ${\mathcal X}$ is a complex vector space,  we denote by ${\mathcal X}_{{\mathbb R}}$ its {\em real form}, i.e. ${\mathcal X}$ considered as a real vector space.
We denote by $\overline{{\mathcal X}}$ a {\em conjugate vector space} to ${\mathcal X}$, i.e. a complex vector space $\overline{{\mathcal X}}$ with an anti-linear isomorphism
${\mathcal X}\ni x\mapsto \overline{x}\in \overline{{\mathcal X}}$. The {\em canonical conjugate vector space} to ${\mathcal X}$ is simply  the real vector space ${\mathcal X}_{{\mathbb R}}$ equipped with the complex structure $-{\rm i}$, if ${\rm i}$ is the complex structure of ${\mathcal X}$. In this case the map $x\to \overline{x}$ is the identity.
If $a\in L({\mathcal X}_{1}, {\mathcal X}_{2})$, we denote by $\overline{a}\in L(\overline{{\mathcal X}}_{1}, \overline{{\mathcal X}}_{2})$ the linear map defined by:
\begin{equation}\label{toti}
\overline{a}\overline{x}_{1}:= \overline{ax_{1}}, \ \overline{x}_{1}\in \overline{{\mathcal X}}_{1}.
\end{equation}

We denote by ${\mathcal X}^{*}$ the {\em anti-dual} of ${\mathcal X}$, i.e. the space of anti-linear forms on ${\mathcal X}$.
Clearly ${\mathcal X}^{*}$ can be identified with $\overline{{\mathcal X}^{{\scriptscriptstyle\#}}}\sim \overline{{\mathcal X}}^{{\scriptscriptstyle\#}}$.

Sesquilinear forms on ${\mathcal X}$ are identified  with elements of $L({\mathcal X}, {\mathcal X}^{*})$, and we use the notation $(x_{1}|b x_{2})$ or sometimes $\overline{x_{1}} \!\cdot\! b x_{2}$ for $b\in L({\mathcal X}, {\mathcal X}^{*})$, $x_{1}, x_{2}\in {\mathcal X}$.

The space of hermitian (resp. anti-hermitian) sesquilinear forms on ${\mathcal X}$ is denoted by $L_{\rm s}({\mathcal X}, {\mathcal X}^{*})$ (resp. $L_{\rm a}({\mathcal X}, {\mathcal X}^{*})$).

If $q\in L_{\rm s}({\mathcal X}, {\mathcal X}^{*})$ is non-degenerate, i.e. $({\mathcal X}, q)$ is a pseudo-unitary space, we denote by $U({\mathcal X}, q)$ the linear pseudo-unitary group on ${\mathcal X}$.

If $b$ is a bilinear form on  the real vector space ${\mathcal X}$, its canonical sesquilinear extension to ${\mathbb C}{\mathcal X}$ is by definition the sesquilinear form $b_{{\mathbb C}}$ on ${\mathbb C}{\mathcal X}$ given by
\[
(w_1|b_{{\mathbb C}}w_2):=x_1 \!\cdot\! b x_2+y_1  \!\cdot\! b y_2+{\rm i} x_1 \!\cdot\!  b y_2 - {\rm i} y_1  \!\cdot\! b x_2, \quad w_i=x_i+{\rm i} y_i
\]
for $x_i,y_i\in{\mathcal X}$, $i=1,2$. This extension maps (anti-)symmetric forms on ${\mathcal X}$ onto (anti-)hermitian forms on ${\mathbb C}{\mathcal X}$.

Conversely if ${\mathcal X}$ is a complex vector space and ${\mathcal X}_{{\mathbb R}}$ is its real form, i.e. ${\mathcal X}$ considered as a real vector space,  then for $b\in L_{\rm s/a}({\mathcal X}, {\mathcal X}^{*})$ the form ${\rm Re }b$ belongs to $L_{\rm s/a}({\mathcal X}_{{\mathbb R}}, {\mathcal X}_{{\mathbb R}}^{{\scriptscriptstyle\#}})$.

\subsection{Bosonic quasi-free states:  neutral case}\label{sec0.1}
Let $({\mathcal X}, \sigma)$ be a real symplectic space, i.e. a pair consisting of a real vector space ${\mathcal X}$ and a non-degenerate anti-symmetric form $\sigma\in L_{\rm a}({\mathcal X}, {\mathcal X}^{{\scriptscriptstyle\#}})$.

We denote $\mathcal{A}({\mathcal X},\sigma)$ the Weyl CCR $C^*$-algebra of $({\mathcal X}, \sigma)$, formally generated by elements of the form $W(y)$ for $y\in{\mathcal X}$,  with:
\[
W(y)^* = W(-y), \ \ \  W(x)W(y) = {\rm e}^{-{\rm i} ( x \!\cdot\!  \sigma y)/2} W(x+y), \ \ \ x,y\in{\mathcal X}.
\]

\begin{definition}A state $\omega$ on $\mathcal{A}({\mathcal X},\sigma)$ is called a (bosonic, neutral) quasi-free state if there is a symmetric form $\eta$ (called  the {\em covariance} of $\omega$) on ${\mathcal X}$ such that
\[
\omega(W(x))={\rm e}^{-{\frac{1}{2}} x \!\cdot\! \eta x}, \ \ \ \ \ x\in{\mathcal X}.
\]
\end{definition}
A quasi-free state $\omega$ on $\mathcal{A}({\mathcal X}, \sigma)$ is regular, i.e. the field operators $\phi(x)$ are well-defined as selfadjoint operators in the GNS representation of $\omega$ with:
\[
[\phi(x_{1}), \phi(x_{2})]= {\rm i} x_{1} \!\cdot\!  \sigma x_{2}\bbbone,\hbox{ as quadratic forms on }{\rm Dom} \phi(x_{1})\cap {\rm Dom} \phi(x_{2}),
\]
and:
\begin{equation}
\label{e0.000}
\omega(\phi(x_{1})\phi(x_{2}))= x_{1}  \!\cdot\! \eta x_{2}+ \frac{{\rm i}}{2} x_{1} \!\cdot\!  \sigma x_{2}, \quad x_{1}, x_{2}\in {\mathcal X}.
\end{equation}
It is convenient to introduce the sesquilinear hermitian form
\[
q:= {\rm i} \sigma_{{\mathbb C}},
\]
usually called the {\em charge} and
\[
\lambda_{\pm}:= \eta_{{\mathbb C}}\pm  \frac{1}{2} q \in L_{\rm s}({\mathbb C}{\mathcal X}, {\mathbb C}{\mathcal X}^{*}).
\]
The following results are well-known (see e.g. \cite[Chap. 17]{DG}).
\begin{proposition}\label{prop:defomega} Let $\eta\in L_{\rm s}({\mathcal X}, {\mathcal X}^{{\scriptscriptstyle\#}})$. Then  the following are equivalent:
\ben
\item $\eta$ is the covariance of a quasi-free state on $\mathcal{A}({\mathcal X},\sigma)$,
\item
$x   \!\cdot\! \eta x\geq 0$, $|x_{1}  \!\cdot\! \sigma x_{2}|\leq 2 (x_{1} \!\cdot\! \eta x_{1})^{\frac{1}{2}}(x_{2} \eta x_{2})^{\frac{1}{2}}$, $x_{1},x_{2}\in{\mathcal X}$,
\item $\lambda_{\pm}\geq 0$ on ${\mathbb C}{\mathcal X}$.
\een
\end{proposition}

\begin{proposition}\label{ista}
Let $\eta\in L_{\rm s}({\mathcal X}, {\mathcal X}^{{\scriptscriptstyle\#}})$.  Then  the following are equivalent:
\ben
\item $\eta$ is the covariance of a {\em pure} quasi-free state on $\mathcal{A}({\mathcal X},\sigma)$,
\item   $(2\eta, \sigma)$ is {\em K\"{a}hler}, i.e.  there exists an anti-involution $ {\rm j}_1\in Sp({\mathcal X}, \sigma)$ such that  $2\eta= \sigma  {\rm j}_1$.
\een
\end{proposition}
\begin{proposition}\label{istaa}
 Let $\eta_{1}$, $\eta_{2}$ be covariances of two pure quasi-free states on $\mathcal{A}({\mathcal X},\sigma)$. Then there exists $r\in Sp({\mathcal X}, \sigma)$ such that $\eta_{2}= r^{{\scriptscriptstyle\#}}\eta_{1}r$.
\end{proposition}

\subsection{ Bosonic quasi-free states: charged case}\label{sec0.2}

Let us now consider the case of a complex symplectic space $({\mathcal Y}, \sigma)$, i.e.  a pair consisting of a complex vector space ${\mathcal Y}$ and a non-degenerate anti-hermitian form $\sigma\in L_{\rm a}({\mathcal Y}, {\mathcal Y}^{*})$.
The complex structure on $ {\mathcal Y}$ will be denoted by $ {\rm j}$, to distinguish it from the complex number $ {\rm i}\in {\mathbb C}$.

 As before we introduce  the charge $q:= {\rm i} \sigma$ which is  hermitian.

Note that  $({\mathcal Y}_{{\mathbb R}}, {\rm Re}\sigma)$ is a real symplectic space called the {\em real form } of $({\mathcal Y}, \sigma)$, with ${\rm j}\in Sp({\mathcal Y}_{{\mathbb R}}, {\rm Re}\sigma)$.

Conversely if $({\mathcal X}, \sigma)$ is a real symplectic space equipped with an anti-involution  ${\rm j}\in Sp({\mathcal X}, \sigma)$, then denoting by ${\mathcal Y}$ the space ${\mathcal X}$ equipped with the complex structure ${\rm j}$ and setting $(x_{1}| \hat{\sigma}x_{2}):= x_{1} \!\cdot\! \sigma x_{2}- {\rm i}  x_{1}  \!\cdot\! \sigma {\rm j} x_{2}$, the space $({\mathcal Y}, \hat{\sigma})$ is a complex symplectic space whose real form is $({\mathcal X}, \sigma)$.

For coherence of notation we will denote the Weyl CCR algebra $\mathcal{A}({\mathcal Y}_{{\mathbb R}}, {\rm Re} \sigma)$ by $\mathcal{A}({\mathcal Y}, \sigma)$.

Let us now consider a quasi-free state $\omega$ on $\mathcal{A}({\mathcal Y}, q)$, as in Subsect. \ref{sec0.1}. The state $\omega$ is called \emph{gauge-invariant} if
\[
\omega(W(y))=\omega(W({\rm e}^{{\rm j}\theta}y)), \ \ 0\leq\theta<2\pi, \ \ y\in{\mathcal Y}.
\]

If the state $\omega$ is not gauge-invariant, the complex structure ${\rm j}$ plays no role and one can forget it. One is then reduced to the situation of Subsect. \ref{sec0.1}.

If $\eta$ is the covariance of $\omega$ then $\omega$ is gauge-invariant iff ${\rm j}\in U({\mathcal Y}_{{\mathbb R}}, \eta)$. One can then uniquely associate to $\eta$ a ${\rm j}-$sesquilinear hermitian form $\hat{\eta}$ defined by
\begin{equation}
\label{tati}
(y_{1}| \hat{\eta}y_{2}):= y_{1} \!\cdot\! \eta y_{2}-{\rm i} y_{1} \!\cdot\!  \eta{\rm j} y_{2}, \ y_{1}, y_{2}\in {\mathcal Y}.
\end{equation}

Let now $\phi(y)$ for $y\in {\mathcal Y}$ be  the  selfadjoint fields in the GNS representation of $\omega$. One can  introduce  the {\em charged fields}:
\[
\psi(y):=\frac{1}{\sqrt{2}}(\phi(y)+ {\rm i} \phi({\rm j} y)), \ \psi^{*}(y):=\frac{1}{\sqrt{2}}(\phi(y)- {\rm i} \phi({\rm j} y)), \ y\in {\mathcal Y}.
\]
The map ${\mathcal Y}\ni y\mapsto \psi^{*}(y)$ (resp. ${\mathcal Y}\ni y\mapsto \psi(y)$) is ${\mathbb C}-$linear (resp. ${\mathbb C}-$anti-linear). The commutation relations take the form:
 \[
[\psi(y_{1}), \psi(y_{2})]= [\psi^{*}(y_{1}), \psi^{*}(y_{2})]=0,  \ [\psi(y_{1}), \psi^{*}(y_{2})]=  (y_{1}| qy_{2})\bbbone, \ y_{1}, y_{2}\in {\mathcal Y}.
\]
If $\omega$ is a gauge-invariant quasi-free state on $\mathcal{A}({\mathcal Y}, q)$, then:
\[
\omega(\psi(y_{1})\psi(y_{2}))= \omega(\psi^{*}(y_{1})\psi^{*}(y_{2}))=0, \ y_{1}, y_{2}\in {\mathcal Y},
\]
and we set:
\begin{equation}\label{e0.tutu}
\begin{array}{rl}
\omega(\psi(y_{1})\psi^{*}(y_{2}))=: &(y_{1}| \lambda_{+} y_{2}),\\[2mm]
\omega(\psi^{*}(y_{2})\psi(y_{1}))=: &(y_{1}| \lambda_{-} y_{2}), \ y_{1}, y_{2}\in {\mathcal Y}.
\end{array}
\end{equation}
Clearly $\lambda_{+}-\lambda_{-}=q$.

We will call $\lambda_{\pm}\in  L_{\rm s}({\mathcal Y}, {\mathcal Y}^{*})$ the
{\em complex covariances } of the gauge invariant quasi-free state $\omega$.

 Introducing the selfadjoint fields $\phi(y)$ we obtain that
\[
\omega(\phi(y_{1})\phi(y_{2}))= {\rm Re}(y_{1}| (\lambda- \frac{1}{2} q)y_{2})+ \frac{{\rm i}}{2}{\rm Re}(y_{1}| \sigma y_{2}).
\]
Therefore we have \begin{equation}\label{ido}
\eta= {\rm Re}( \lambda_{\pm}\mp\frac{1}{2} q), \ \hat{\eta}= \lambda_{\pm}\mp \frac{1}{2} q.
\end{equation}
In this situation we will call $\eta$ the {\em real covariance} of the state $\omega$, to distinguish it from the complex covariances $\lambda_{\pm}$.

 \begin{remark} the state $\omega$ is of course  uniquely determined by either $\lambda_{+}$ or $\lambda_{-}$, but later  conditions on a state $\omega$ look nicer  when formulated in terms of the pair of  covariances $\lambda_{\pm}$.
Note that  $\lambda_{-}$ is usually called the {\em  charge density} associated to  $\omega$.
\end{remark}

The following  propositions are the analogues of Props. \ref{prop:defomega}, \ref{ista}, \ref{istaa}. We sketch their proofs for the reader's convenience.
\begin{proposition}\label{isto}
 Let $\lambda_{\pm}\in L_{\rm s}({\mathcal Y}, {\mathcal Y}^{*})$. Then the following are equivalent:
 \ben
 \item $\lambda_{\pm}$ are  the covariances of a gauge-invariant quasi-free state on $\mathcal{A}({\mathcal Y}, q)$,
 \item $\lambda_{\pm}\geq 0$ and $\lambda_{+}- \lambda_{-}=q$.
 \een
\end{proposition}
\proof
Since $\omega$ is gauge-invariant we have
\[
{\rm j}\in O({\mathcal Y}_{{\mathbb R}}, \eta)\cap Sp({\mathcal Y}_{{\mathbb R}}, {\rm Re}\sigma)= O({\mathcal Y}_{{\mathbb R}}, \eta)\cap O({\mathcal Y}_{{\mathbb R}}, {\rm Re}q).
\]
From this fact  and (\ref{ido}) we deduce that $\eta\geq 0\Leftrightarrow\lambda_{+}\geq \frac{1}{2} q$, and that the  second condition in Prop. \ref{prop:defomega} (with $\sigma$ replaced by ${\rm Re}\sigma$) is equivalent to
\[
\pm q\leq 2\lambda_{+} -q\ \Leftrightarrow \ \lambda_{\pm}\geq 0.
\]
This completes the proof of the proposition. $\Box$\medskip
\begin{proposition}\label{istu}
Let $\lambda_{\pm}\in L_{\rm s}({\mathcal Y}, {\mathcal Y}^{*})$. Then the following are equivalent:
 \ben
 \item $\lambda_{\pm}$ are  the covariances of a {\em pure} gauge-invariant quasi-free state on $\mathcal{A}({\mathcal Y}, q)$,
 \item there exists  an involution $\kappa\in U({\mathcal Y}, q)$  such that $q\kappa\geq 0$ and $\lambda_{\pm}= \frac{1}{2} q(\kappa\pm \bbbone)$.
 \item $\lambda_{\pm}\geq \pm\frac{1}{2} q$, $ \lambda_{\pm} q^{-1}\lambda_{\pm}= \pm \lambda_{\pm}$, $\lambda_{+}- \lambda_{-}=q$.
 \een
\end{proposition}
\proof   By Prop. \ref{ista} the state $\omega$ is pure iff there exists  an anti-involution $ {\rm j}_1\in Sp({\mathcal Y}_{{\mathbb R}}, {\rm Re}\sigma)$ such that \begin{equation}\label{toto}
2\eta= ({\rm Re}\sigma)  {\rm j}_1.
\end{equation}  Since  ${\rm j}\in O({\mathcal Y}_{{\mathbb R}}, \eta)\cap Sp({\mathcal Y}_{{\mathbb R}}, {\rm Re}\sigma)$ we obtain  that $ {\rm j}_1\in U({\mathcal Y}, q)$, i.e. $ {\rm j}_1$ is ${\mathbb C}-$linear and pseudo-unitary for $q$. From (\ref{toto}) we then get that $2\lambda_{+}-q= \sigma {\rm j}_1$. Setting $\kappa= -{\rm j} {\rm j}_1$ we see that $\kappa\in U({\mathcal Y}, q)$ and $\lambda_{+}= \frac{1}{2} q(\bbbone + \kappa)$.
Therefore (1) is equivalent to
\[
(4)\ \lambda_{+}\geq 0, \ \lambda_{+}\geq q, \ \lambda_{+}= \frac{1}{2} q(\bbbone + \kappa), \ \kappa^{2}= \bbbone,  \kappa\in U({\mathcal Y}, q).
\]
(4) clearly implies (2). Let us prove the converse implication.  Set $P_{\pm}:= \frac{1}{2}(\bbbone \pm \kappa)$. Clearly $P_{\pm}$ are projections with $P_{\pm}^{*}q= qP_{\pm}$,  $\kappa P_{\pm}= \pm P_{\pm}$, and
\[
\lambda_{+}\geq 0, \ \lambda_{+}\geq q \Leftrightarrow \pm qP_{\pm}\geq 0.
 \]
 Now we have
 \[
 qP_{\pm}= qP_{\pm}^{2}= P_{\pm}^{*}qP_{\pm}= \pm P_{\pm}^{*}q\kappa P_{\pm},
\]
which completes the proof since $q\kappa\geq 0$.  The fact that $(2)$ and $(3)$ are equivalent  is an easy computation. $\Box$\medskip

\begin{proposition}
 Let $\lambda_{\pm}$, $\tilde{\lambda}_{\pm}$ be the covariances of two pure, gauge-invariant quasi-free states on $\mathcal{A}({\mathcal Y}, q)$. Then there exists $r\in U({\mathcal Y}, q)$ such that $\tilde{\lambda}_{\pm}= r^{*}\lambda_{\pm}r$.
\end{proposition}
\proof We introduce the real covariances $\eta$, $\tilde\eta$. By Prop. \ref{istaa} there exists $r\in Sp({\mathcal Y}_{{\mathbb R}}, {\rm Re}\sigma)$ with $\tilde{\eta}= r^{{\scriptscriptstyle\#}}\eta r$. Using the gauge-invariance of the two states  we obtain that $r{\rm j}= {\rm j} r$, hence $r\in U({\mathcal Y}, q)$.  From this we easily obtain the proposition.  $\Box$\medskip
\subsection{Charge reversal}
\begin{definition}
 Let $({\mathcal Y}, \sigma)$ a complex symplectic space. \ben
 \item a map $\chi\in L({\mathcal Y}_{{\mathbb R}})$ is called a {\em charge reversal} if $\chi^{2}= \bbbone$ or $\chi^{2}= -\bbbone$ and
 \[
(\chi y_{1}| \sigma \chi y_{2})= (y_{2}| \sigma y_{1}), \ y_{1}, y_{2}\in{\mathcal Y}.
\]
\item
a gauge invariant quasi-free state $\omega$ on $\mathcal{A}({\mathcal Y}, q)$ is {\em invariant }under  the charge reversal $\chi$ if
\[
\omega(\phi(\chi y_{1})\phi(\chi y_{2}))= \omega(\phi(y_{1})\phi(y_{2})), \ y_{1}, y_{2}\in {\mathcal Y},
\]
or equivalently:
\[
\omega(\psi(\chi y_{1})\psi^{*}(\chi y_{2}))= \omega(\psi^{*}(y_{1})\psi(y_{2})),  \ y_{1}, y_{2}\in {\mathcal Y}.
\]
\een
\end{definition}
Note that  a charge reversal is automatically anti-linear. The last condition above can be rephrased as
\begin{equation}\label{reversal}
 \chi^{*}\lambda_{\pm} \chi= \lambda_{\mp}.
\end{equation}
\subsection{Complexification of   bosonic quasi-free states}\label{complix}
\subsubsection{Neutral case}
Let now $({\mathcal X}, \sigma)$ be a real symplectic space.  We   equip ${\mathbb C}{\mathcal X}$ with $\sigma_{{\mathbb C}}$, obtaining a complex symplectic space.  We set as in Subsect. \ref{sec0.2} $q= {\rm i} \sigma_{{\mathbb C}}$. The canonical complex conjugation on ${\mathbb C}{\mathcal X}$ is a charge reversal on $({\mathbb C}{\mathcal X}, \sigma_{{\mathbb C}})$.

Clearly $(({\mathbb C}{\mathcal X})_{{\mathbb R}}, {\rm Re} \sigma_{{\mathbb C}})$ is isomorphic to $({\mathcal X}\oplus {\mathcal X}, \sigma\oplus \sigma)$ as  real symplectic spaces.
If $\omega$ is a quasi-free state on $({\mathcal X}, \sigma)$ with covariance $\eta$, then we can consider the quasi-free state $\tilde{\omega}$ on $\mathcal{A}(({\mathbb C}{\mathcal X})_{{\mathbb R}}, {\rm Re} \sigma_{{\mathbb C}})$ with covariance ${\rm Re} \eta_{{\mathbb C}}$.

It is easy to see that $\tilde{\omega}$ is gauge-invariant with  covariances $\lambda_{\pm}$  equal to
\[
\lambda_{\pm}= \eta_{{\mathbb C}}\pm \frac{1}{2} q.
\]
Moreover $\tilde{\omega}$ is invariant under charge reversal.

Therefore by complexifying a quasi-free state $\omega$ on a real symplectic space $({\mathcal X}, \sigma)$, we obtain a gauge-invariant quasi-free state $\tilde{\omega}$ on $\mathcal{A}({\mathbb C}{\mathcal X}, \sigma_{{\mathbb C}})$. It follows that, possibly after complexifying the real symplectic space $({\mathcal X}, \sigma)$, one can always restrict the discussion to gauge-invariant  quasi-free states.

\subsubsection{Charged case}\label{chargi}
Assume now that $({\mathcal X}, \sigma)$ is a complex symplectic space, with complex structure ${\rm j}$ and  $\omega$ be a gauge-invariant quasi-free state on $\mathcal{A}({\mathcal X}, \sigma)$, with real covariance $\eta$.

We can apply the above procedure to  the real symplectic space $({\mathcal X}_{{\mathbb R}}, {\rm Re} \sigma)$.  Then ${\mathbb C} {\mathcal X}$ has two complex structures, the canonical one ${\rm i}$ and ${\rm j}_{{\mathbb C}}$,  (the complexification of ${\rm j}$). As is well known, (see e.g. \cite[Sect. 1.3.6]{DG}) ${\mathbb C}{\mathcal X}$ splits as the direct sum
\[
{\mathbb C}{\mathcal X}= {\mathcal Z}\oplus \overline{{\mathcal Z}},
\]
where
\[
 {\mathcal Z}:= {\rm Ker}({\rm j}_{{\mathbb C}} -{\rm i} ), \ \overline{{\mathcal Z}}:= {\rm Ker}({\rm j}_{{\mathbb C}}+{\rm i})
\]
are called the {\em holomorphic} resp. {\em anti-holomorphic subspaces} of ${\mathbb C}{\mathcal X}$.   Note that the natural conjugation on ${\mathbb C}{\mathcal X}$ maps bijectively ${\mathcal Z}$ onto $\overline{{\mathcal Z}}$. Then (see \cite[Sect. 17.1.2]{DG}) the sequilinear extensions $\eta_{{\mathbb C}}$ and $({\rm Re} \sigma)_{{\mathbb C}}$ are reduced w.r.t. the direct sum ${\mathcal Z}\oplus \overline{{\mathcal Z}}$, i.e.:
\[
({\rm Re} \sigma)_{{\mathbb C}}:= \mat{\sigma_{{\mathcal Z}}}{0}{0}{\overline{\sigma}_{{\mathcal Z}}}, \ \eta_{{\mathbb C}}:= \mat{\eta_{{\mathcal Z}}}{0}{0}{\overline{\eta}_{{\mathcal Z}}},
\]
where $\sigma_{{\mathcal Z}}\in L_{\rm a}({\mathcal Z}, {\mathcal Z}^{*})$, $\eta_{{\mathcal Z}}\in L_{\rm s}({\mathcal Z}, {\mathcal Z}^{*})$ and if $a\in L({\mathcal Z}, {\mathcal Z}^{*})$ we define $\overline{a}\in L(\overline{{\mathcal Z}}, \overline{{\mathcal Z}}^{*})$ by $(\overline{z}_{1}| \overline{a}\overline{z}_{1}):= \overline{(z_{1}| a z_{2})}$.

Let us  remark that one can identify ${\mathcal X}$ with ${\mathcal Z}$ by the (${\mathbb C}-$linear) map
\[
T: {\mathcal X}\ni x\mapsto Tx= \frac{1}{\sqrt{2}}(x- {\rm i} {\rm j} x)\in {\mathcal Z},
\]
and also $\overline{{\mathcal X}}$ with $\overline{{\mathcal Z}}$ by $\overline{T}$.
Under these identifications,  $\sigma_{{\mathcal Z}}$  is identified with $\sigma$,  $\eta_{{\mathcal Z}}$ with $\hat{\eta}$, where $\hat{\eta}$ is the ${\rm j}-$sesquilinear extension of $\eta$ defined in (\ref{tati}). It follows then from (\ref{ido}) that   the sesquilinear form $\eta_{{\mathbb C}}+ \frac{{\rm i}}{2}({\rm Re} \sigma)_{{\mathbb C}}$ acting on ${\mathbb C}{\mathcal X}$ is identified with
\[
\mat{\lambda_{+}}{0}{0}{\overline{\lambda_{-}}}, \hbox{ acting on }{\mathcal X}\oplus \overline{{\mathcal X}}.
\]

\section{Hadamard  states}\label{sec22}
\subsection{Klein-Gordon equations on a globally hyperbolic space time}\label{sec22.0}
\subsubsection{Notation}
 Let $M$ be a smooth manifold.  As usual  ${\mathcal E}(M)$ is the space of smooth, (complex valued) functions on $M$,  ${\mathcal D}(M)\subset {\mathcal E}(M)$ the space of smooth  compactly supported functions on $M$, ${\mathcal D}'(M)$ the space of distributions on $M$ and ${\mathcal E}'(M)$ the space of compactly supported distributions.   We denote by $\langle u|v\rangle$, for $u\in {\mathcal D}(M)$ (resp. ${\mathcal E}(M)$), $v\in {\mathcal D}'(M)$ (resp. ${\mathcal E}'(M)$) the (bilinear) duality bracket.
 \subsubsection{Klein-Gordon equations}
 Assume $M$ is equipped with a Lorentzian metric $g=g_{\mu\nu}dx^{\mu}dx^{\nu}$ such that $(M, g)$  is a globally hyperbolic space-time.  We use the convention $(-, +,\cdots ,+)$ for the signature.

We use the notations \[
|g|:= {\rm det}[g_{\mu\nu}], \quad [g^{\mu\nu}]:= [g_{\mu\nu}]^{-1}, \quad dv:= |g|^{\frac{1}{2}}dx.
\]

If $S$ is a Cauchy hypersurface, we denote by $n^{\nu}$ the unit  future directed normal vector field to $S$ (after choosing  a time orientation), and by $ds$ the surface measure on $S$ obtained from $dv$.

We fix a smooth vector potential $A_{\mu}(x) dx^{\mu}$ and a smooth function $\rho: M\to {\mathbb R}$.
The associated Klein-Gordon operator is:
\begin{equation}\label{e0.00}
P(x, D_{x})= |g|^{-\frac{1}{2}}( \partial_{\mu}+ {\rm i} A_{\mu})|g|^{\frac{1}{2}}g^{\mu\nu} ( \partial_{\nu}+ {\rm i} A_{\nu})+\rho.
\end{equation}
We equip ${\mathcal D}(M)$ with the scalar product
\[
(u_{1}|u_{2})= \int_{M} \overline{u}_{1}u_{2} dv,
\]
so that $P(x, D_{x})$ is formally selfadjoint. We denote by
$E_{\pm}$ the {\em retarded/advanced fundamental solutions} of $P(x, D_{x})$,
 and by $E= E_{+}- E_{-}$ the {\em Pauli-Jordan commutator function} Recall that $ E_{\pm}^{*}= E_{\mp}$ hence $E= - E^{*}$.

A function $u$ on $M$ is called {\em space-compact} if the intersection of ${\rm supp}\, u$ with any Cauchy hypersurface of $M$ is compact. The space of smooth space-compact functions will be denoted by $C^\infty_{\rm sc}(M)$.

We denote by ${\rm Sol}_{\rm sc}(P)\subset C^\infty(M)$ the space of   smooth space-compact solutions of
\[
{\rm (KG)} \quad P(x, D_{x})\phi=0.
\]
One has (see e.g. \cite{BGP}):
\begin{equation}\label{e0.1}
 E {\mathcal D}(M)={\rm Sol}_{\rm sc}(P), \ {\rm Ker} E= P{\mathcal D}(M).
\end{equation}
Moreover if  we  fix a Cauchy hypersurface $S$ and set
\[
\begin{array}{rcl}
\rho:& {\rm Sol}_{\rm sc}(P)&\to {\mathcal D}(S)\oplus {\mathcal D}(S)\\[2mm]
&\phi&\mapsto (\phi_{\mid S}, {\rm i}^{-1}n^{\mu}(\nabla_{\mu}+ {\rm i} A_{\mu})\phi_{\mid S})=: (\rho_{0}\phi, \rho_{1}\phi),
\end{array}
\]
then $\rho: {\rm Sol}_{\rm sc}(P)\to {\mathcal D}(S)\oplus {\mathcal D}(S)$ is bijective (see e.g. \cite{BGP}).

Let $\varsigma$ be the sesquilinear form on ${\mathcal D}(M)/P{\mathcal D}(M)$ defined by
\[
([u]|\varsigma[v]):=\bra \overline{u}, E v\rangle =\bra E , \overline{u}\otimes v\rangle, \quad u,v\in{\mathcal D}(M).
\]
By construction $({\mathcal D}(M)/P{\mathcal D}(M),\varsigma)$ is a complex symplectic space.  Setting also
\[
(f|\sigma g):= -{\rm i}\int_{S}( \overline{f_{0}}g_{1}+ \overline{f_{1}}g_{0})ds, \  f,g\in{\mathcal D}(S)\oplus {\mathcal D}(S),
\]
we have
\begin{equation}\label{e0.2}
([u]|\varsigma [v])= (\rho \circ E u|\sigma\rho \circ E v),
\end{equation}
i.e.
\[
\rho\circ E: ({\mathcal D}(M)/P{\mathcal D}(M), \varsigma)\to ({\mathcal D}(S)\oplus {\mathcal D}(S), \sigma)
\]
is a symplectomorphism.

\subsection{The wave front set}\label{sec-1}
Let $M$ be a smooth manifold,  $\co{M}$ its cotangent bundle of $M$. The zero section  of $\co{M}$  will be denoted by $Z$.

\subsubsection{Operations on conic sets}\label{sec-1.1}
A set $\Gamma\subset \coo{M}$ is {\em conic} if
\[
(x, \xi)\in \Gamma\Rightarrow (x, t \xi)\in \Gamma, \ \forall \ t>0.
\]
If  $\Gamma\subset \coo{M}$ is conic, we set:
\[
\overline{\Gamma}:=\{(x, -\xi): \  (x, \xi)\in \Gamma\}.
\]
Let $M_{i}$, $i=1,2$ be two manifolds, $Z_{i}$ the zero section of $\co{M_{i}}$ and $\Gamma\subset \coo{(M_{1}\times M_{2})}$ be a conic set.   The elements of $\coo{(M_{1}\times M_{2})}$
 will be denoted by $(x_{1}, \xi_{1}, x_{2}, \xi_{2})$, which allows to
 consider $\Gamma$ as a relation  between
 ${\rm T}^{*}M_{2}$ and $T^{*}M_{1}$, still denoted by $\Gamma$. Clearly $\Gamma$ maps conic sets into conic sets.
 We set:
 \[
\begin{array}{l}
\Gamma':= \{(x_{1}, \xi_{1}, x_{2}, - \xi_{2}) :\ (x_{1}, \xi_{1}, x_{2}, \xi_{2})\in \Gamma\}\subset \coo{(M_{1}\times M_{2})},\\[2mm]
{\rm Exch}(\Gamma):= \{(x_{2}, \xi_{2}, x_{1}, \xi_{1}) : \ (x_{1}, \xi_{1}, x_{2}, \xi_{2})\in \Gamma\}\subset \coo{(M_{2}\times M_{1})},\\[2mm]
_{M_{1}}\!\Gamma:=\{(x_{1}, \xi_{1}) : \ \exists \ x_{2}\hbox{ such that } (x_{1}, \xi_{1}, x_{2},0)\in \Gamma\}= \Gamma(Z_{2})\subset \coo{M_{1}}_{1},\\[2mm]
\Gamma\!_{M_{2}}:=\{(x_{2}, \xi_{2}) : \ \exists \ x_{1}\hbox{ such that } (x_{1}, 0, x_{2},\xi_{2})\in \Gamma\}= \Gamma^{-1}(Z_{1})\subset\coo{M_{2}}_{2}.
\end{array}
\]

\subsubsection{Properties of the wave front set}\label{sec-1.2}
Recall that if $u\in {\mathcal D}'(M)$ then the wave front set ${\rm WF}(u)$ is a conic subset of $\coo{M}$. We refer to \cite{hoermander1} for the exact definition and the proof of the following basic properties:

\medskip

\ben
\item {\bf Complex conjugation}: if $u\in {\mathcal D}'(M)$ then ${\rm WF}(\overline{u})= \overline{{\rm WF}(u)}$.

\item {\bf Restriction to a sub-manifold}: let $S\subset M$ a sub-manifold. The {\em co-normal bundle} to $S$ in $M$ is:
\[
T_{S}^{*}M:=\{(x, \xi)\in \coo{M} : \ x\in S, \ \xi\cdot v=0\ \forall v\in T_{x}S\}.
\]
If $u\in{\mathcal D}'(M)$, the restriction $u_{|S}$ of $u$ to $S$ is well defined if
${\rm WF}(u)\cap T_{S}^{*}M=\emptyset$.  One has
\[
{\rm WF}(u_{|S})\subset \{(x, \xi_{|T_{x}S}) : \ x\in S, \ (x, \xi)\in {\rm WF}(u)\}.
\]
\item {\bf Kernels}: let $K: {\mathcal D}(M_{2})\to {\mathcal D}'(M_{1})$ be linear continuous and denote  by $K(x_{1}, x_{2})\in {\mathcal D}'(M_{1}\times M_{2})$ its distributional kernel. Then $Ku$ is well defined for $u\in {\mathcal E}'(M_{2})$ if ${\rm WF}(u)\cap {\rm WF}(K)'\!_{M_{2}}=\emptyset$ and in such case
\[
{\rm WF}(Ku)\subset \, _{M_{1}}\!{\rm WF}(K)\cup  {\rm WF}(K)'\circ {\rm WF}(u).
\]
\item {\bf Composition}: let  $K_{1}\in {\mathcal D}'(M_{1}\times M_{2})$, $K_{2}\in {\mathcal D}'(M_{2}\times M_{3})$, where $K_{2}$ is properly supported, i.e. the projection:
${\rm supp}\,K_{2}\to M_{2}$ is proper. Then $K_{1}\circ K_{2}$ is well defined if
\[
{\rm WF}(K_{1})'\!_{M_{2}}\, \cap \, _{M_{2}}\!{\rm WF}(K_{2})'= \emptyset,
\]
and in such case
\[
{\rm WF}(K_{1}\circ K_{2})'\subset {\rm WF}(K_{1})'\circ {\rm WF}(K_{2})'\ \cup \ _{M_{1}}\!{\rm WF}(K_{1})'\times Z_{3}\ \cup \ Z_{1}\times {\rm WF}(K_{2})'\!_{M_{3}}.
\]
\item {\bf Adjoint:} let us denote by $K^{*}$ the adjoint of $K$ with respect to any smooth non-vanishing density $dx$ on $M$. Then
\[
{\rm WF}(K^{*})'= {\rm Exch}({\rm WF} (K)').
\]
\een

\subsection{Distinguished parametrices and microlocal spectrum condition}
\subsubsection{Distinguished parametrices}
Let us recall basic elements of the theory of distinguished parametrices of Duistermaat and H\"ormander  \cite{duistermaat-hormander}  for the case of the Klein-Gordon operator $P(x,D)$. The {\em characteristic manifold} of $P(x,D)$ is
\[
{\mathcal N}:=\{ (x,\xi)\in \coo{M} : \ p(x,\xi)=0 \},
\]
where $p(x,\xi)=g^{\mu\nu}(x)\xi_{\mu}\xi_{\nu}$ is the principal symbol of $P(x,D)$.

We use the notation $X=(x,\xi)$ for points in $\coo{M}$. We write $X_1\sim X_2$ if $X_1=(x_1,\xi_1)$ and $X_2=(x_2,\xi_2)$ are in ${\mathcal N}$ and $X_1$ and $X_2$ are on the same Hamiltonian curve for $p$.

Let us denote  by $V_{x\pm}\subset T_x M$ for $x\in M$, the open future/past light cones and $V_{x\pm}^*$ the dual cones
\[
V^{*}_{x\pm}:=\{ \xi\in T^*_{x}M : \ \xi\cdot v >0, \ \forall v\in V_{x\pm}, \ v\neq 0 \}.
\]
The set ${\mathcal N}$ has two connected components invariant under the Hamiltonian flow of $p$, namely:
\[
{\mathcal N}_{\pm}:=\{ X\in {\mathcal N} : \ \xi\in V^*_{x\pm}\}.
\]

Recall that $E_{\pm}$ denote respectively the retarded and advanced fundamental solution. We denote respectively $E_{\rm  F}$, $E_{\overline{\rm  F}}$ the Feynman and anti-Feynman parametrix. The theory of Duistermaat and H\"ormander provides among others a description of the wave front sets of the parametrices $E_{\pm}$, $E_{\rm  F}$, $E_{\overline{\rm  F}}$ and establishes their uniqueness up to smooth functions.
The proof of the next lemma can be found for instance in \cite[Thm. 2.29]{junker}.

\begin{lemma}\label{lem:raff}We have:
\begin{eqnarray*}
(1) & {\rm WF}(E)'=\{(X_{1}, X_{2})\in {\mathcal N}\times {\mathcal N} : \  \ X_{1}\sim X_{2}\},\\
(2) & {\rm WF}(E_{\rm +}-E_{\rm F})'=\{(X_{1}, X_{2})\in {\mathcal N}_{-}\times {\mathcal N}_{-} : \  \ X_{1}\sim X_{2}\},\\
(3) & {\rm WF}(E_{\rm -}-E_{\rm F})'=\{(X_{1}, X_{2})\in {\mathcal N}_{+}\times {\mathcal N}_{+} : \  \ X_{1}\sim X_{2}\}.
\end{eqnarray*}
\end{lemma}
\subsubsection{Microlocal spectrum condition}
We are now ready to formulate the microlocal spectrum condition.

Let us  fix a gauge-invariant, quasi-free state $\omega$ on $\mathcal{A}({\mathcal D}(S)\oplus {\mathcal D}(S), \sigma)$. We denote by $\eta$ its real covariance and $\lambda_{\pm}$ its two complex covariances (see Subsect. \ref{sec0.2}).   We assume that
\[
\eta: {\mathcal D}(S)\oplus {\mathcal D}(S)\to {\mathcal D}'(S)\oplus {\mathcal D}'(S)
\]
(and hence $\lambda_{\pm}$) is continuous. Note that $\lambda_{\pm}$ are ${\mathbb C}-$linear, while $\eta$ is only ${\mathbb R}-$linear.

We associate to $\eta$, $\lambda_{\pm}$ linear operators  $H$, $\Lambda_{\pm}\in L({\mathcal D}(M), {\mathcal D}'(M))$ as follows:
\begin{equation}\label{defino}
\begin{array}{rl}
\langle u_{1}| H u_{2}\rangle& := (\rho\circ E )u_{1}\!\cdot\! \eta  (\rho\circ E) u_{2}, \\[2mm]
 (u_{1}| \Lambda_{\pm} u_{2})&:= \left( (\rho\circ E )u_{1}| \lambda_{\pm}  (\rho\circ E) u_{2}\right).
\end{array}
\end{equation}
Note that $\Lambda_{\pm}$ are ${\mathbb C}-$linear, while $H$ is only ${\mathbb R}-$linear.  From Subsects. \ref{sec0.2}, \ref{sec22.0} we have
$\Lambda_{+}- \Lambda_{-}={\rm i} E$.
\begin{definition}\label{def2}
 Let $\Lambda_{\pm}: {\mathcal D}(M)\to {\mathcal D}'(M)$ be as in (\ref{defino}). Then $\Lambda_{\pm}$ satisfies the {\em microlocal spectrum condition} if
\begin{equation}\label{eq:msc}\tag{$\mu$sc}
{\rm WF}(\Lambda_{\pm})'\subset \{(X_{1}, X_{2})\in {\mathcal N}_{\pm}\times {\mathcal N}_{\pm} : \  \ X_{1}\sim X_{2}\}.
\end{equation}
A gauge invariant quasi-free state is a {\em Hadamard state} if its complex covariances $\Lambda_{\pm}$ satisfy the microlocal spectrum condition.
\end{definition}

\begin{remark}
  The above form of the Hadamard condition for charged fields seems to have been first noticed by Hollands \cite{hollands} for Dirac fields. Let us explain its relationship with the standard form in \cite{radzikowski}: in the standard form, symplectic spaces (real or complex) are considered as real symplectic spaces, and a state is  labelled by its real covariance. The standard Hadamard condition takes then the form
  \begin{equation}\label{radiko}
  {\rm WF}(H_{{\mathbb C}}+ \frac{{\rm i}}{2}({\rm Re} E)_{{\mathbb C}})'\subset  \{(X_{1}, X_{2})\in {\mathcal N}_{+}\times {\mathcal N}_{+} : \  \ X_{1}\sim X_{2}\}.
\end{equation}
Note that if we consider a charged Klein-Gordon equation, $H_{{\mathbb C}}+ \frac{{\rm i}}{2}({\rm Re} E)_{{\mathbb C}}$ is actually a $2\times 2$ matrix of distributions on $M\times M$, its wave front set being then the union of the wave front set of its entries.

We can apply  the discussion of Subsubsect. \ref{chargi} and identify $H_{{\mathbb C}}+ \frac{{\rm i}}{2}({\rm Re} E)_{{\mathbb C}}$ with the matrix $\mat{\Lambda_{+}}{0}{0}{\overline{\Lambda_{-}}}$. Using that ${\rm WF}(\overline{\Lambda_{-}})'= \overline{{\rm WF}(\Lambda_{-})'}$, we see that (\ref{radiko}) is equivalent to (\ref{eq:msc}).

 One can also show that the inclusion in (\ref{eq:msc}) and (\ref{radiko}) can be replaced by an equality.
\end{remark}
\begin{remark}
 Assume now that we consider a neutral Klein-Gordon equation, i.e. that $A_{\mu}=0$. Then the Klein-Gordon operator $P$ is real, and the map $u\mapsto \overline{u}$  is a charge reversal of the symplectic spaces $({\rm Sol}_{\rm sc}(P), \sigma)$, $({\mathcal D}(M)/P{\mathcal D}(M), \varsigma)$.  It follows then from (\ref{reversal}) that $\overline{\Lambda}_{\pm}= \Lambda_{\mp}$, so each of the two conditions in (\ref{eq:msc}) implies the other.
\end{remark}

\section{Background on pseudo-differential  calculus}\label{sec1}
\subsection{Notation}\label{sec1.0}
- If $f:{\mathbb R}_{t}\times {\mathbb R}^{n}_{x}\to {\mathbb C}$ is a function, and $t\in {\mathbb R}$ we denote by $f(t)$ the function:
\[
f(t): \ {\mathbb R}^{n}\ni x\mapsto f(t, x)\in {\mathbb C}.
\]

- We denote by $C^{\infty}_{\rm bd}({\mathbb R}^{n})$ the space of smooth functions on ${\mathbb R}^{n}$ uniformly bounded with all derivatives. We equip $C^{\infty}_{\rm bd}({\mathbb R}^{n})$ with its canonical Fr\'echet space structure. We denote by $H^{m}({\mathbb R}^{d})$  the Sobolev space of order $m\in {\mathbb R}$.

- We denote by ${\mathcal S}({\mathbb R}^{d})$, resp. ${\mathcal S}'({\mathbb R}^{d})$ the space of Schwartz functions, resp. distributions on ${\mathbb R}^{d}$.

- We set
\[
 {\mathcal H}({\mathbb R}^{d}):=\cap_{m\in {\mathbb R}}H^{m}({\mathbb R}^{d}), \ {\mathcal H}'({\mathbb R}^{d}):=\cup_{m\in {\mathbb R}}H^{m}({\mathbb R}^{d}),
 \]
  equipped with their canonical topologies.

- If $E, F$ are two topological vector spaces, we write $A: E\to F$ if $A$ is linear continuous from $E$ to $F$.

- We set $D_{x}= {\rm i}^{-1} \partial_{x}$, $\langle x\rangle= (1+x^{2})^{\frac{1}{2}}$, $x\in {\mathbb R}^{d}$.
\subsection{Symbol classes}\label{sec1.1}
We denote by $S^{m}\!({\mathbb R}^{2d})$, $m\in {\mathbb R}$ the  symbol class
\begin{equation}\label{e1.1}
a\in S^{m}\!({\mathbb R}^{2d}) \ \ \hbox{if} \ \  \partial_{x}^{\alpha} \partial_{k}^{\beta}a(x, k)\in O(\langle k\rangle^{m-|\beta|}), \ \alpha, \beta\in {\mathbb N}^{d}.
\end{equation}
Similarly we will denote by $S^{m}\!({\mathbb R})$ the class
\begin{equation}
\label{e1.1b}
f\in S^{m}\!({\mathbb R}) \hbox{ if } \partial_{\lambda}^{\alpha}f(\lambda)\in O(\langle \lambda\rangle^{m-\alpha}), \ \alpha\in {\mathbb N}.
\end{equation}
We denote by $S^{m}_{\rm h}\!({\mathbb R}^{2d})$ the subspace of $S^{m}\!({\mathbb R}^{2d})$ of symbols homogeneous of degree $m$ in the $k$ variable, (outside a neighborhood of the origin) :
\begin{equation}\label{e1.2}
a\in S^{m}_{\rm h}\!({\mathbb R}^{2d})\hbox{ if }a\in S^{m}\!({\mathbb R}^{2d})\hbox{ and }a(x, \lambda k)= \lambda^{m}a(x, k), \ \lambda\geq 1, \  |k|\geq 1.
\end{equation}
We set
\[
S^{-\infty}({\mathbb R}^{2d}):= \bigcap_{m\in {\mathbb R}}S^{m}({\mathbb R}^{2d}).
\]
If $a_{m-k}\in S^{m-k}\!({\mathbb R}^{2d})$ for $k\in {\mathbb N}$ and $a\in S^{m}\!({\mathbb R}^{2d})$ we write
\[
a\sim \sum_{k=0}^{\infty}a_{m-k}
\]
if
\begin{equation}
\label{e1.3}
a-\sum_{k=0}^{n}a_{m-k}\in S^{m-n-1}\!({\mathbb R}^{2d}), \ \forall n\in {\mathbb N}.
\end{equation}
Note that if $a_{m-k}\in S^{m-k}\!({\mathbb R}^{2d})$ for $k\in {\mathbb N}$, then it is well-known that there exists $a\in S^{m}({\mathbb R}^{2d})$, unique modulo $S^{-\infty}({\mathbb R}^{2d})$, such that $a\sim \sum_{k=0}^{\infty}a_{m-k}$.

We say that a symbol $a\in S^{m}({\mathbb R}^{2d})$ is {\em poly-homogeneous} if $a\sim\sum_{k=0}^{\infty}a_{m-k}$ for $a_{m-k}\in S^{m-k}_{\rm h}\!({\mathbb R}^{2d})$.The symbols  $a_{m-k}$ are then   clearly unique modulo $S^{-\infty}({\mathbb R}^{2d})$.

The   subspace of {\em poly-homogeneous symbols} of degree $m$ will be denoted by $S^{m}_{\rm ph}\!({\mathbb R}^{2d})$.
The space $S_{\rm ph}^{m}({\mathbb R})\subset S^{m}({\mathbb R})$ is defined similarly.

We will often write  $S^{m}_{({\rm ph})}$  for $S^{m}_{({\rm ph})}\!({\mathbb R}^{2d})$. We equip  $S^{m}_{({\rm ph})}\!({\mathbb R}^{2d})$ with the Fr\'echet space topology given by the semi-norms:
\[
\| a\|_{m, N}:= \sup_{|\alpha|+ |\beta|\leq N}| \langle k\rangle^{-m+ |\beta|} \partial_{x}^{\alpha}
 \partial_{k}^{\beta}a|.
\]
We set
\[
S^{\infty}_{({\rm ph})}({\mathbb R}^{2d}):=\bigcup_{m\in {\mathbb R}}S^{m}_{({\rm ph})}({\mathbb R}^{2d}).
\]
\subsection{Principal symbol and characteristic set}\label{sec1.2}

The {\em principal symbol} of $a\in S^{m}$, denoted by $\sigma_{\rm pr}(a)$ is the equivalence class $a+ S^{m-1}$ in $S^{m}/S^{m-1}$. If $a\in S^{m}_{\rm ph}$ then $a+ S^{m-1}$ has a unique representative in $S^{m}_{\rm h}$, namely the function $a_{m}$ in (\ref{e1.3}).
Therefore  in this case the principal symbol is a  function on $ {\mathbb R}^{2d}$, homogeneous of degree $m$ in $k$.

The {\em characteristic set} of $a\in S^{m}_{\rm ph}$ is defined as
\begin{equation}
\label{e1.4}
{\rm Char}(a):= \{(x, k)\in\cooo{{\mathbb R}^{d}} : \  a_{m}(x, k)= 0\},
\end{equation}
it is clearly conic  in the $k$ variable.

A symbol $a\in S^{m}\!({\mathbb R}^{2d})$ is {\em elliptic} if  there exists $C, R>0$ such that
\[
 |a(x, k)|\geq C \langle k\rangle^{m}, \ |k|\geq R.
\]
Clearly $a\in S^{m}_{\rm ph}\!({\mathbb R}^{2d})$ is elliptic iff ${\rm Char}(a)= \emptyset$.

\subsection{Pseudo-differential operators}\label{sec1.3}

In this subsection we collect some well-known results about pseudo-differential  calculus.

For $a\in S^{m}\!({\mathbb R}^{2d})$, we denote by  ${\rm Op}^{\rm w}(a)$ the {\em Weyl quantization} of $a$ defined by:
\[
{\rm Op}^{\rm w}(a)u(x)=  a^{\rm w}(x, D)u(x):=(2\pi)^{-d}\iint{\rm e}^{{\rm i}
(x-y)k}a(\frac{x+ y}{2}, k)u(y){\rm d} y{\rm d} k.
\]
One has ${\rm Op}^{\rm w}(a): {\mathcal H}({\mathbb R}^{d})\to {\mathcal H}({\mathbb R}^{d})$  and
\[
{\rm Op}^{\rm w}(a)^{*}= {\rm Op}^{\rm w}(\overline{a}),
\]
 hence ${\rm Op}^{\rm w} (a):  {\mathcal H}'({\mathbb R}^{d})\to {\mathcal H}'({\mathbb R}^{d})$.

We denote by  $\Psi^{m}_{({\rm ph})}({\mathbb R}^{d})$ the space ${\rm Op}^{\rm w} (S^{m}_{({\rm ph})}\!({\mathbb R}^{2d}))$ and set
\[
\Psi_{({\rm ph})}^{\infty}({\mathbb R}^{d})= \bigcup_{m\in {\mathbb R}}\Psi^{m}_{({\rm ph})}({\mathbb R}^{d}), \ \Psi^{-\infty}({\mathbb R}^{d})= \bigcap_{m\in {\mathbb R}}\Psi^{m}({\mathbb R}^{d}).
\]
We will often write $\Psi^{m}_{({\rm ph})}$ instead of $\Psi^{m}_{({\rm ph})}({\mathbb R}^{d})$. We will equip $\Psi^{m}_{({\rm ph})}({\mathbb R}^{d})$ with the  Fr\'{e}chet space topology obtained from the topology of $S^{m}({\mathbb R}^{2d})$.

If $a= a^{\rm w}\!(x, D_{x})\in \Psi^{m}_{\rm ph}({\mathbb R}^{d})$ the $m-$homogeneous   function
$\sigma_{\rm pr}(a)=:a_{m}(x, k)$ is called the {\em principal symbol} of $a$.

Let $s,m\in {\mathbb R}$. Then   the map
\begin{equation}\label{e0.0}
S^{m}\!({\mathbb R}^{d})\ni a\mapsto {\rm Op}^{\rm w}(a)\in B(H^{s}({\mathbb R}^{d}), H^{s-m}({\mathbb R}^{d}))
\end{equation}
is continuous.

 An operator ${\rm Op}^{\rm w}\!(a)\in \Psi^{m}({\mathbb R}^{2d})$ is {\em elliptic} if its symbol $a(x, k)$ is elliptic in $S^{m}\!({\mathbb R}^{2d})$.  If $a\in \Psi^{m}$ is elliptic
 then there exists $b\in \Psi^{-m}$, unique modulo $\Psi^{-\infty}$ such that $ab=ba=\bbbone$ modulo $\Psi^{-\infty}$. Such an operator $b$ is called a {\em pseudo-inverse} or a {\em parametrix} of $a$. We will denote it by $b=:a^{(-1)}$. As a typical example  $\bbbone + b$ for $b\in \Psi^{-m}$, $m>0$ is elliptic in $\Psi^{0}$.

The following lemma is proven in  \cite[Prop. A.1.1, A.1.2]{S} .

\begin{lemma}\label{1.0}
 Let $u\in {\mathcal D}'({\mathbb R}^{d})$, $(x_{0}, k_{0})\in \cooo{{\mathbb R}^{d}}$. Then $(x_{0}, k_{0})\not \in {\rm WF}(u)$ iff there exists $\chi\in C_0^\infty({\mathbb R}^{d})$  and $a\in S^{0}_{\rm ph}\!({\mathbb R}^{2d})$ with $\chi(x_{0})\neq 0$, $(x_{0}, k_{0})\not\in {\rm Char}(a)$ and ${\rm Op}^{\rm w}\!(a)\chi u\in {\mathcal S}({\mathbb R}^{d})$.
\end{lemma}
\subsection{Functional calculus for pseudo-differential  operators}\label{sec1.4}

We now recall various well-known results about functional calculus for pseudo-differential  operators.
\begin{proposition}\label{1.1}
Let $a\in \Psi^{m}({\mathbb R}^{d})$  for $m\geq 0$ be elliptic in $\Psi^{m}({\mathbb R}^{d})$ and symmetric on ${\mathcal S}({\mathbb R}^{d})$. Then:
\ben \item $a$ is selfadjoint on $H^{m}({\mathbb R}^{d})$,
\item  Denote by $\rho(a)$ the resolvent set of $a$, with domain $H^{m}({\mathbb R}^{d})$. Then for $z\in \rho(a)$ $(z-a)^{-1}\in \Psi^{-m}({\mathbb R}^{d})$,
\item if $f\in S^{p}({\mathbb R})$, $p\in {\mathbb R}$,  then $f(a)$, defined by the functional calculus, belongs to $\Psi^{mp}({\mathbb R}^{d})$.
\item if  $f$ is elliptic in $S^{p}({\mathbb R})$ then $\sigma_{\rm pr}(f(a))= f(\sigma_{\rm pr}(a))$ mod $S^{mp-1}({\mathbb R}^{2d})$.
\item if $a\in \Psi^{m}_{\rm ph}({\mathbb R}^{d})$ and $f\in S^{p}_{\rm ph}({\mathbb R})$, then $f(a)\in \Psi^{mp}_{\rm ph}({\mathbb R}^{d})$, and if $f\in S^{p}_{\rm ph}({\mathbb R})$ is elliptic, then $\sigma_{\rm pr}(f(a))= f_{p}(\sigma_{\rm pr}(a))$, where $f_{p}\in S^{p}_{\rm h}({\mathbb R})$ is the principal symbol of $f$.
\een
\end{proposition}
\proof
We refer the reader for example to \cite[Thm. 5.4]{Ro}, \cite[Corr. 4.5]{Bo} for the proof of similar statements.   Statement (1) follows from the fact that $a+ {\rm i} \lambda \bbbone$ maps $H^{m}({\mathbb R}^{d})$ bijectively onto $L^{2}({\mathbb R}^{d})$ for $|\lambda|$ large enough.

To prove statement (2), the most direct way is to use the Beals criterion (see e.g. \cite{Bo}), which characterizes pseudo-differential  operators by properties of the multi-commutators with the operators $x_{i}$, $D_{j}$:  an operator $a$ belongs to $\Psi^{m}({\mathbb R}^{d})$ iff:
\[
\begin{array}{l}
a: {\mathcal S}({\mathbb R}^{d})\to {\mathcal S}({\mathbb R}^{d}),\\[2mm]
{\rm ad}_{x}^{\alpha}{\rm ad}_{D}^{\beta}a\in B(H^{m- |\alpha|}({\mathbb R}^{d}), L^{2}({\mathbb R}^{d})), \ \forall\ \alpha, \beta\in {\mathbb N}^{d},
\end{array}
\]
where ${\rm ad}_{x_{i}}b= [x_{i}, b]$, ${\rm ad}_{D_{j}}b= [D_{j}, b]$, and ${\rm ad}_{x}^{\alpha}= {\rm ad}_{x_{1}}^{\alpha_{1}}\cdots {\rm ad}_{x_{d}}^{\alpha_{d}}$ and similarly for ${\rm ad}_{D}^{\beta}$.

From (2) one can deduce (3) by expressing $f(a)$ for $f\in S^{p}({\mathbb R})$ using the resolvent $(z- a)^{-1}$ and an almost analytic extension of $f$, see e.g. \cite{HS, D}.  Statement (4) follows from the parametrix construction of $(z-a)^{-1}$, which has $(z- \sigma_{\rm pr}(a))^{-1}$ as principal symbol. Statement (5) can be proved similarly. $\Box$\medskip

We conclude this subsection by stating a useful lemma, which follows from symbolic calculus.

\begin{lemma}\label{idiotic}
 Let $a\in \Psi^{p}({\mathbb R}^{d})$, $p\in {\mathbb R}$,  $f,g\in C^\infty({\mathbb R}^{d})$ with $\nabla f, \nabla g\in C_0^\infty({\mathbb R}^{d})$ and $f\equiv 0$ near ${\rm supp}\,g$. Then  \[
f(x)a g(x)\in \Psi^{-\infty}({\mathbb R}^{d}).
\]
In particular $f(x)a g(x)$ maps ${\mathcal H}'({\mathbb R}^{d})$ into ${\mathcal H}({\mathbb R}^{d})$.
\end{lemma}

\subsection{Propagators}\label{sec1.5}

Let us fix a  map $\epsilon(t)=\epsilon_{1}(t)+ \epsilon_{0}(t)$, where $\epsilon_{i}(t)\in C^{\infty}({\mathbb R}, \Psi^{i}({\mathbb R}^{d}))$ for $i=0,1$. Assume moreover that
$\epsilon_{1}(t)$  is elliptic in $\Psi^{1}({\mathbb R}^{d})$ and symmetric on ${\mathcal S}({\mathbb R}^{d})$. It follows by Prop. \ref{1.1} that  $\epsilon_{1}(t)$
is selfadjoint with domain $H^{1}({\mathbb R}^{d})$, hence $\epsilon(t)$ with domain $H^{1}({\mathbb R}^{d})$ is closed, with non empty resolvent set.

We denote by $\Texp(\int_{s}^{t}{\rm i}\epsilon(\sigma)d\sigma)$ the associated  propagator defined by:

\[
\left\{
\begin{array}{rl}
&\frac{ \partial}{ \partial t}\Texp(\int_{s}^{t}{\rm i}\epsilon(\sigma)d\sigma)= {\rm i}  \epsilon(t)\Texp(\int_{s}^{t}{\rm i}\epsilon(\sigma)d\sigma),\\[2mm]
&\frac{ \partial}{ \partial s}\Texp(\int_{s}^{t}{\rm i}\epsilon(\sigma)d\sigma)= -{\rm i}  \Texp(\int_{s}^{t}{\rm i}\epsilon(\sigma)d\sigma)\epsilon(s),\\[2mm]
&\Texp(\int_{s}^{s}{\rm i}\epsilon(\sigma)d\sigma)=\bbbone.
\end{array}
\right.
\]
Note that   the propagator $\Texp(\int_{s}^{t}{\rm i}\epsilon_{1}(\sigma)d\sigma)$ exists  and is unitary by e.g. \cite[Thm. X.70]{RS2}. Since $\epsilon(t)- \epsilon_{1}(t)$ is locally uniformly bounded in $B(L^{2}({\mathbb R}^{d}))$, one easily deduces  the existence of $\Texp(\int_{s}^{t}{\rm i}\epsilon(\sigma)d\sigma)$, which is strongly continuous in $(t,s)$ with values in $B(L^{2}({\mathbb R}^{d}))$.

\begin{definition}\label{def1}
 Assume in addition that $\epsilon(t)\in \Psi^{1}_{\rm ph}({\mathbb R}^{d})$.  Then we denote by $\Phi_{\epsilon}(t,s): \cooo{{\mathbb R}^{d}}\to \cooo{{\mathbb R}^{d}}$ the symplectic flow associated to  the time-dependent Hamiltonian $ -\sigma_{\rm pr}(\epsilon)(t, x, k)$.

  Equivalently $ \Phi_{\epsilon}(t,s)$  is the restriction to the variables $(x, k)$ of the symplectic flow on $\cooo{{\mathbb R}^{1+d}}$ associated to the Hamiltonian $\tau- \sigma_{\rm pr}(\epsilon)(t, x, k)$.
\end{definition}
Clearly $\Phi_{\epsilon}(t,s)$ is  an homogeneous map of degree $0$.

The following classical result is known as {\em Egorov's theorem}, see for instance \cite[Sec. 0.9]{taylor} for the proof.
\begin{proposition}\label{1.2}
\ben
\item  $\Texp({\textstyle\int_{s}^{t}}{\rm i}\epsilon(\sigma)d\sigma)$ is bounded on ${\mathcal H}({\mathbb R}^{d})$ hence on ${\mathcal H}'({\mathbb R}^{d})$ by duality.

\item  Let $a\in \Psi^{m}({\mathbb R}^{d})$. Then
 \[
a(t,s):=  \Texp({\textstyle\int_{s}^{t}}{\rm i}\epsilon(\sigma)d\sigma)a\Texp({\textstyle\int_{t}^{s}}{\rm i}\epsilon(\sigma)d\sigma)
\]
belongs to $C^{\infty}({\mathbb R}^{2}, \Psi^{m}({\mathbb R}^{d}))$.   Moreover if $\epsilon(t)\in C^{\infty}({\mathbb R}, \Psi^{1}_{\rm ph}({\mathbb R}^{d}))$ and $a\in \Psi^{m}_{\rm ph}({\mathbb R}^{d})$ then $a(t,s)\in C^{\infty}({\mathbb R}^{2}, \Psi_{\rm ph}^{m}({\mathbb R}^{d}))$
 with
 \[
\sigma_{\rm pr}(a)(t,s)= \sigma_{\rm pr}(a)\circ \Phi_{\epsilon}(s,t).
\]
\een
\end{proposition}
From Proposition \ref{1.2} and Lemma \ref{1.0} we obtain the following result (the steps of the proof are explained in \cite[Sec. 0.10]{taylor}).
\begin{proposition}\label{1.2b}
For $u\in {\mathcal H}'({\mathbb R}^{d})$ one has:
\[
{\rm WF}(\Texp({\textstyle\int_{s}^{t}}{\rm i}\epsilon(\sigma)d\sigma)u)= \Phi_{\epsilon}(t,s){\rm WF}(u),
\]
hence
\[
{\rm WF}(\Texp({\textstyle\int_{s}^{t}}{\rm i}\epsilon(\sigma)d\sigma))'= \{(x, k, x', k')  : \ (x, k)= \Phi_{\epsilon}(t,s)(x', k')\}.
\]
\end{proposition}
\begin{lemma}\label{buda}
 Let $\epsilon(t)\in C^\infty({\mathbb R}, \Psi^{1}({\mathbb R}^{d}))$ as above,  $s_{-\infty}(t,s)\in C^{\infty}({\mathbb R}^{2}, \Psi^{-\infty}({\mathbb R}^{d}))$. Then
 \[
\Texp({\textstyle\int_{s}^{t}}{\rm i}\epsilon(\sigma)d\sigma) s_{-\infty}(t,s)\in C^{\infty}({\mathbb R}^{2}, \Psi^{-\infty}({\mathbb R}^{d})).
\]
\end{lemma}
\proof The proof will be given in Subsect. \ref{budaproof}. $\Box$\medskip

\section{Concrete Klein-Gordon equations}\label{sec2}\setcounter{equation}{0}

\subsection{Model Klein-Gordon equation on ${\mathbb R}^{1+d}$}\label{sec2.0b}
In this subsection we describe the model Klein-Gordon equation that will be considered in the sequel.  We take $M= {\mathbb R}^{1+d}$, $x= (t, {\rm x})\in {\mathbb R}^{1+d}$
 and fix a second order differential operator
\begin{equation}\label{e0.3b}
a(t, {\rm x}, D_{{\rm x}})= -\sum_{j,k=1}^{d} \partial_{{\rm x}^{j}}a^{jk}(x) \partial_{{\rm x}^{k}}+ \sum_{j=1}^{d} b^{j}(x) \partial_{{\rm x}^{j}}-  \partial_{{\rm x}^{j}}\overline{b}^{j}(x)+ m(x),
\end{equation}
where
\begin{equation}
\label{e0.3}
\begin{array}{l}
a^{jk}, b^{j}, m \in C^{\infty}({\mathbb R}, C^{\infty}_{\rm bd}({\mathbb R}^{d})), \   m(x)\in {\mathbb R},\\[2mm]
[a^{jk}](x)\geq  c(t)\bbbone\hbox{ uniformly on }{\mathbb R}^{1+d}, \ c(t)>0.
\end{array}
\end{equation}

 We introduce the model Klein-Gordon operator
 \[
P(x, D_{x})=  \partial_{t}^{2}+ a(t, {\rm x}, D_{{\rm x}}),
\]
which is formally selfadjoint for the scalar product $(u_{1}| u_{2})= \int_{{\mathbb R}^{1+d}} \overline{u}_{1}u_{2}dx$.

We will consider the Cauchy problem:
\begin{equation}
\label{e3.1}
\left\{
\begin{array}{rl}
& \partial_{t}^{2}\phi(t)+a(t, {\rm x}, D_{{\rm x}})\phi(t)=0,\\[2mm]
&\phi(s)=f_{0},\\[2mm]
&{\rm i}^{-1} \partial_{t} \phi(s)=f_{1},
\end{array}
\right.
\end{equation}
for $f= (f_{0}, f_{1})\in {\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$. By the well-known method of energy estimates, one obtains the existence and uniqueness of a solution $\phi\in C^\infty({\mathbb R}, {\mathcal H}({\mathbb R}^{d}))$. Similarly if $f\in {\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}$, there exists a unique solution $\phi\in C^\infty({\mathbb R}, {\mathcal H}'({\mathbb R}^{d}))$.

\subsection{Reduction to the model case}\label{ss:modelcase}
Consider a globally hyperbolic space-time $(M, g)$ with a Cauchy hypersurface diffeomorphic to ${\mathbb R}^{d}$.
This implies that we can assume that $M= {\mathbb R}_{t}\times {\mathbb R}^{d}$ and
\begin{equation}\label{metric}
g= -c (x)dt^{2}+ h_{jk}(x)d{\rm x}^{j}d{\rm x}^{k},
\end{equation}
where $x= (t, {\rm x})$, $c(x)>0$ is a smooth function and $h_{jk}(x)d{\rm x}^{j}d{\rm x}^{k}$ is a smooth Riemannian metric on ${\mathbb R}^{d}$.

In this subsection we explain how to reduce the  Klein-Gordon operator (\ref{e0.00}) to the model case considered in Subsect. \ref{sec2.0b}.

Writing $A_{\mu}(x)= (V(x), \rA_{j}(x))$, we have:
\[\begin{array}{rl}
&P(x, D_{x})\\[2mm]
=&\ c^{-\frac{1}{2}}|h|^{-\frac{1}{2}}( \partial_{t}+ {\rm i} V) c^{-\frac{1}{2}}|h|^{\frac{1}{2}}( \partial_{t}+{\rm i} V)\\[2mm]
&- c^{-\frac{1}{2}}|h|^{-\frac{1}{2}}( \partial_{j}+{\rm i} A_{j})c^{\frac{1}{2}}| h|^{\frac{1}{2}}h^{jk}( \partial_{k}+{\rm i} A_{k})+ \rho,
\end{array}
\]
where $|h|= \det[h_{jk}]$, $ [h^{jk}]= [h_{jk}]^{-1}$.

We choose the Cauchy hypersurface $S=\{0\}\times {\mathbb R}^{d}$ so that
\begin{equation}\label{e2.0}
\overline\phi_{1}\sigma \phi_{2}= \int_{S}\left(\overline{( \partial_{t}+{\rm i} V)\phi_{1}}\phi_{2}- \overline{\phi_{1}}( \partial_{t}+ {\rm i} V)\phi_{2}\right)c^{-\frac{1}{2}}| h|^{\frac{1}{2}}d {\rm x}.
\end{equation}
Set:
\[
\begin{array}{rl}
F(t, {\rm x})=&  \int_{0}^{t}V(s, {\rm x})ds, \ \tilde{\rA}=\rA+ \nabla F,\  \tilde{\rho}= c \rho- c^{\frac{1}{4}}|h|^{-\frac{1}{4}} \partial_{t}^{2}(c^{-\frac{1}{4}}|h|^{\frac{1}{4}}),\\[2mm]
a(t, {\rm x}, D_{{\rm x}})= &c^{\frac{1}{4}}|h|^{-\frac{1}{4}}( \partial_{j}+{\rm i} \tilde{\rA}_{j})c^{\frac{1}{2}}|h|^{\frac{1}{2}}h^{jk}( \partial_{k}+{\rm i} \tilde{\rA}_{k})c^{\frac{1}{4}}|h|^{-\frac{1}{4}}+ \tilde{\rho}.
\end{array}
\]
\begin{lemma}\label{2.1}
\ben
\item Let $b= c^{-\frac{1}{4}}|h|^{\frac{1}{4}}$. Then:
\[
P(x, D_{x})= c^{-\frac{1}{2}}|h|^{-\frac{1}{2}} b{\rm e}^{-{\rm i} F}\left( \partial_{t}^{2}+ a(t, {\rm x}, D_{{\rm x}})\right)b{\rm e}^{{\rm i} F},
\]
hence
\[
 E= b^{-1}{\rm e}^{- {\rm i} F}\tilde{E}c^{\frac{1}{2}}|h|^{\frac{1}{2}}b^{-1}{\rm e}^{{\rm i} F},
\]
where $\tilde{E}$ is the Pauli-Jordan function for $\tilde{P}(x, D_{x}):=  \partial_{t}^{2}+ a(t, {\rm x}, D_{{\rm x}})$.
\item The map:
\[
\phi\mapsto \tilde{\phi}:=c^{-\frac{1}{4}}|h|^{\frac{1}{4}}{\rm e}^{-{\rm i} F}\phi,
\]
is symplectic from  $({\rm Sol}_{\rm sc}(P), \sigma)$  to $({\rm Sol}_{\rm sc}(\tilde{P}), \tilde{\sigma})$ for
\[
\overline{\tilde{\phi}_{1}}\tilde{\sigma}\tilde{\phi}_{2}=  \int_{S}\overline{ \partial_{t}\tilde{\phi}_{1}}\tilde{\phi}_{2}-  \overline{\tilde{\phi}_{1}} \partial_{t}\tilde{\phi}_{2}d{\rm x}.
\]
\item Assume the following hypotheses:
\[
(H) \begin{array}{rl}
&\forall \ I\subset {\mathbb R}\hbox{ compact interval }\exists\  C>0\hbox{ such that }\\[2mm]
&C\leq c(x), \ C\bbbone\leq [h_{jk}(x)],\hbox{uniformly for }x\in I\times {\mathbb R}^{d},\\[2mm]
&  h_{jk}(x), \ c(x), \ \rho(x), A_{\mu}(x)\in C^{\infty}({\mathbb R}, C^{\infty}_{\rm bd}({\mathbb R}^{d})).
\end{array}
\]
Then the operator $a(t, {\rm x}, D_{{\rm x}})$ is of the form (\ref{e0.3b}) and the conditions (\ref{e0.3}) are satisfied.
\een
\end{lemma}
\proof
Set
\[
b^{2}= c^{-\frac{1}{2}}|h|^{\frac{1}{2}}, \ a= ( \partial_{j}+{\rm i} A_{j})c^{\frac{1}{2}}| h|^{\frac{1}{2}}h^{jk}( \partial_{k}+{\rm i} A_{k})+ c^{\frac{1}{2}}|h|^{\frac{1}{2}}\rho.
\]
Then
\[
\begin{array}{rl}
P(x, D_{x})= &c^{-\frac{1}{2}}|h|^{-\frac{1}{2}}(( \partial_{t}+ {\rm i} V)b^{2}( \partial_{t}+ {\rm i} V)+a),\\[2mm]
=&c^{-\frac{1}{2}}|h|^{-\frac{1}{2}}b(( \partial_{t}+ {\rm i} V)^{2}+ \tilde{a})b,
\end{array}
\]
for $ \tilde{a}=  b^{-1} ab^{-1}- b^{-1}( \partial_{t}^{2}b)$.    Since $ \partial_{t}+ {\rm i} V= {\rm e}^{- {\rm i} F} \partial_{t}{\rm e}^{{\rm i} F}$ we finally get:
\[
P(x, D_{x})= c^{- \frac{1}{2}}|h|^{-\frac{1}{2}} b{\rm e}^{- {\rm i} F}\left( \partial_{t}^{2}+ a(t, {\rm x}, D_{{\rm x}})\right){\rm e}^{{\rm i} F}b,
\]
  This proves (1). (2) and (3) are  left to the reader. $\Box$\medskip

By Lemma \ref{2.1}, the task of constructing Hadamard states for $P(x,D)$ can be reduced to constructing Hadamard states for the model Klein--Gordon equation $\tilde P(x,D)$. Indeed, suppose we have a Hadamard state with two-point function $\tilde{\Lambda}$ and charge ${\rm i} \tilde E$. Then
\[
\Lambda:= b^{-1}{\rm e}^{- {\rm i} F}\tilde{\Lambda}c^{\frac{1}{2}}|h|^{\frac{1}{2}}b^{-1}{\rm e}^{{\rm i} F}
\]
defines the two-point function of a Hadamard state with charge ${\rm i}  E$, the wave front set being preserved by multiplication by smooth densities.

\section{The model Klein-Gordon equation}\label{sec3}\setcounter{equation}{0}
In this section we consider the model Klein-Gordon operator $P(x, D_{x})=  \partial_{t}^{2}+ a(t, {\rm x}, D_{{\rm x}})$ introduced in Subsect. \ref{sec2.0b}.  The associated symplectic form is:
\[
\overline{\phi}_{1}\sigma\phi_{2}=  \int_{{t}\times {\mathbb R}^{d}}\overline{ \partial_{t}\phi_{1}}\phi_{2}-  \overline{\phi_{1}} \partial_{t}\phi_{2}d{\rm x}.
\]
\subsection{Notation}
If $a\geq 0$ is a selfadjoint operator on a Hilbert space $\mathfrak{h}$ we write $a>0$ if ${\rm Ker}\, a= \{0\}$. Then $a^{-1}$ with domain ${\rm Ran}\,a$ is selfadjoint.

If $a,b\geq 0$ are two selfadjoint operators  on a Hilbert space $\mathfrak{h}$ then we write $a\leq b$ if ${\rm Dom} b^{\frac{1}{2}}\subset {\rm Dom} a^{\frac{1}{2}}$ and $(u| au)\leq (u| b u)$ for $u\in {\rm Dom} b^{\frac{1}{2}}$.

The Kato-Heinz theorem  says that if $0\leq a\leq b$ then $0\leq a^{s}\leq b^{s}$ for all $0\leq s \leq 1$ and if
$0<a\leq b$ then $0\leq b^{-1}\leq a^{-1}$.

We write $a\sim b$ if $c^{-1}a\leq b\leq c a$ for some $c>0$.
\subsection{Bicharacteristic curves}
We denote
$a_{2}(t, {\rm x}, k)= k_{i}a^{ij}(t, {\rm x})k_{j}$ the principal symbol of $a(t, {\rm x}, D_{{\rm x}})$ and by $p(x,\xi)= - \tau^{2}+ a_{2}(t, {\rm x}, k)$ the principal symbol of $P(x, D_{x})=  \partial_{t}^{2}+ a(t,{\rm x}, D_{{\rm x}})$.

We set:
\[
\epsilon_{1}(t, {\rm x}, k)= (k_{i}a^{ij}(t, {\rm x})k_{j})^{\frac{1}{2}}.
\]
As in Def. \ref{def1} we denote by $\Phi_{\pm}(t,s): \cooo{{\mathbb R}^{d}}\to \cooo{{\mathbb R}^{d}}$ the restrictions to  the variables
$({\rm x}, k)$ of the symplectic flows on $\cooo{{\mathbb R}^{1+d}}$ associated to the hamiltonians $\tau\mp \epsilon_{1}(t, {\rm x}, k)$.

The following lemma is immediate:
\begin{lemma}\label{2.3}
Let $X_{i}\in {\mathcal N}_{\pm}$, $i=1,2$ with $X_{1}\sim X_{2}$. Then
 \[
X_{i}= (t_{i}, {\rm x}_{i}, \pm\epsilon_{1}(t_{i}, {\rm x}_{i}, k_{i}), k_{i})\hbox{ with }
({\rm x}_{2}, k_{2})= \Phi_{\pm}(t_{2}, t_{1})({\rm x}_{1}, k_{1}).
\]
\end{lemma}

\subsection{Parametrix for the Cauchy problem}\label{sec3.2}
In this subsection we outline the well-known construction of a parametrix for the Cauchy problem (\ref{e3.1}).
The construction is well-known and belongs to the folklore of microlocal analysis. Usually it is done using  Fourier integral operators.
Our construction relies more on Hilbert space methods.

We start by an auxiliary lemma.
\begin{lemma}\label{2.2}
Assume (\ref{e0.3}) and let $a(t, {\rm x}, D_{{\rm x}})$ be given by (\ref{e0.3b}). Then there exists  smooth maps
\[
\begin{array}{rl}
{\mathbb R}\ni t&\mapsto\epsilon(t)= \epsilon(t, {\rm x}, k)\in S^{1}_{\rm ph}({\mathbb R}^{2d}),\\[2mm]
{\mathbb R}\ni t&\mapsto r_{-\infty}(t)= r_{-\infty}(t, {\rm x}, k)\in S^{-\infty}({\mathbb R}^{2d}),
\end{array}
\]
such that:
\ben
\item $\epsilon(t, {\rm x}, k)$ is real-valued, $\epsilon(t, {\rm x}, k)\in C^{\infty}({\mathbb R}, S^{1}_{\rm ph}({\mathbb R}^{2d}))$ with principal symbol
\[
\epsilon_{1}(t, {\rm x}, k)= (k_{i}a^{ij}(t, {\rm x})k_{j})^{\frac{1}{2}},
\]
\item $\epsilon^{\rm w}(t, {\rm x}, D_{{\rm x}})\geq c(t)(D_{x}^{2}+\bbbone)^{\frac{1}{2}}$ for $c(t)>0$,
\item   $a(t, {\rm x}, D_{{\rm x}})=\epsilon^{\rm w}(t, {\rm x}, D_{{\rm x}})^{2}- r_{-\infty}^{\rm w}(t, {\rm x}, D_{{\rm x}})$.
\een
Moreover $\epsilon(t)$ and $r_{-\infty}(t)$ are unique modulo $C^\infty({\mathbb R}, \Psi^{-\infty}({\mathbb R}^{d}))$.
\end{lemma}
\proof  The proof will be given in Subsect. \ref{sec-app1.1bis}. $\Box$\medskip

The following theorem is the main result of this section. It will be used later on to characterize and construct examples of quasi-free Hadamard states.

\begin{theoreme}\label{3.1}
There exist  operators $b(t)\in C^\infty({\mathbb R}, \Psi^{1}_{\rm ph}({\mathbb R}^{d}))$, $r(t)\in C^\infty({\mathbb R}, \Psi^{-1}_{\rm ph}({\mathbb R}^{d}))$  with
\[
\begin{array}{rl}
(i)& b(t)=\epsilon(t)+ (2\epsilon(t))^{-1} {\rm i}  \partial_{t}\epsilon(t)\hbox{ mod }C^\infty({\mathbb R}, \Psi^{-1}({\mathbb R}^{d})),\\[2mm]
(ii)&r(t)=b^{*}(t)^{(-1)}\hbox{ mod }C^\infty({\mathbb R}, \Psi^{-\infty}({\mathbb R}^{d})),\\[2mm]
(iii)&r(t)= \epsilon(t)^{-1}+ C^\infty({\mathbb R}, \Psi^{-2}({\mathbb R}^{d})),\\[2mm]
(iv)&r(t)+ r(t)^{*}\sim \epsilon(t)^{-1},
\end{array}
\]
such that if
\[
\begin{array}{rl}
u_{+}(t,s):=& \Texp({\rm i}\int_{s}^{t}b(\sigma)d\sigma),\ u_{-}(t,s):= \Texp(-{\rm i} \int_{s}^{t}b^{*}(\sigma)d \sigma)\\[2mm]
d_{+}(t):=&(\bbbone +b^{*}(t)^{(-1)}b(t))^{(-1)},\ d_{-}(t):= d_{+}(t)^{*},\\[2mm]
r_{+}(t):=& r(t), \ r_{-}(t):= r^{*}(t),
\end{array}
\]
the following properties hold:
\ben
\item set for $f\in {\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}$:
\[
\begin{array}{rl}
U(t,s)f= &u_{+}(t,s)d_{+}(s)\left(f_{0}+ r_{+}(s)f_{1}\right)+ u_{-}(t,s)d_{-}(s)\left(f_{0}- r_{-}(s)f_{1}\right),\\[2mm]
=:& U_{+}(t,s)f+ U_{-}(t,s)f.
\end{array}
\]
then
\begin{equation}\label{e.buda5}
\left\{\begin{array}{l}
( \partial_{t}^{2}+ a(t, {\rm x}, D_{{\rm x}}))U(t,s)f= s_{-\infty}(t,s)f,\\[2mm]
U(s,s)f= f_{0}+ r_{-\infty,0}(s)f,\\[2mm]
 {\rm i}^{-1} \partial_{t}U(t,s)f|_{t=s}= f_{1}+ r_{-\infty, 1}(s)f,
\end{array}\right.
\end{equation}
for $s_{-\infty}(t,s)\in C^\infty({\mathbb R}^{2}, \Psi^{-\infty}({\mathbb R}^{d}))\otimes {\mathbb C}^{2}$,
$r_{-\infty, i}(s)\in C^\infty({\mathbb R}, \Psi^{-\infty}({\mathbb R}^{d}))\otimes {\mathbb C}^{2}$.

\item  let $\phi(t)$ be the unique solution of (\ref{e3.1}) for $f\in {\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}$. Then:
\begin{equation}\label{e3.02}
\phi(t)=U(t,s)f { \ \rm mod \ } C^{\infty}({\mathbb R}, {\mathcal H}({\mathbb R}^{d})).
\end{equation}
\een
\end{theoreme}
\proof  the proof will be given in Subsect. \ref{sec-app1.2}. $\Box$\medskip

To simplify notation, in the rest of the paper, we will fix $s=0$, and set:
\[
\begin{array}{rl}
&b:=b(0), \ u_{\pm}(t):=u_{\pm}(t,0), \ U(t):= U(t,0),\\[2mm]
&d:=d(0), \ r:= r(0), \ \epsilon:= \epsilon(0).
\end{array}
\]

The constructions of Hadamard states in Sect. \ref{sec4} will a priori depend on the choice of an operator $r$. To study this dependence it is convenient to introduce the following definition.

\begin{definition}\label{defo}
 We denote by ${\mathcal M}$ the set:
 \[
{\mathcal M}:= \{r\in \Psi^{-1}_{\rm ph}({\mathbb R}^{d})\ : \ r= b^{*(-1)}+ \Psi^{-\infty}({\mathbb R}^{d}), \   r+ r^{*}\sim  \epsilon^{-1}\}.
\]
\end{definition}
\begin{remark}
Note that the operator $b$ in Thm. \ref{3.1} is unique modulo $\Psi^{-\infty}$.   Thm. \ref{3.1} implies  that ${\mathcal M}$ is not empty. Since two elements of ${\mathcal M}$ are equal modulo $\Psi^{-\infty}$, the conclusions of Thm. \ref{3.1} are valid for any $r\in {\mathcal M}$.
\end{remark}
The following corollary is a consequence of (\ref{e3.02}), Prop. \ref{1.2b} and Lemma \ref{2.3}.
\begin{corollary}\label{3.1b}
If $\phi(t)$ is the unique solution of (\ref{e3.1}) for $f\in {\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}$, one has
\[
\phi(t)= U_{+}(t)f+ U_{-}(t)f { \ \rm mod \ } C^{\infty}({\mathbb R}, {\mathcal H}({\mathbb R}^{d})),
\]
and\[\begin{array}{rl}
& {\rm WF}(U_{\pm}(t)f)\\[2mm]
=&\{(x_{2}, \xi_{2})  : \ \exists \ ({\rm x}_{1}, k_{1})\in {\rm WF} (f_{0}\pm r_{\pm}f_{1})\hbox{ with }(s, {\rm x}_{1}, \pm \epsilon({\rm x}_{1}, k_{1}), k_{1})\sim (x_{2}, \xi_{2})\}.
\end{array}
\]
In particular \[
{\rm WF}(U_{\pm}(t)f)\subset {\mathcal N}_{\pm}.
\]
\end{corollary}
\subsection{Symplectic properties of the spaces of positive/negative wavefront set solutions}\label{sec3.3}
We now investigate the properties, with respect to the symplectic form $\sigma$, of the spaces of solutions of the Klein-Gordon equation having wavefront set included in  the positive/negative energy surfaces ${\mathcal N}_{\pm}$.

Of course we cannot work with  solutions in ${\mathcal E}(M)$, since their wavefront set is empty, nor with solutions in ${\mathcal D}'(M)$, since the symplectic form $\sigma$ is not defined on arbitrary distributional solutions. A natural space of solutions is the space of {\em finite energy solutions} defined as follows:
\[
{\rm Sol}_{E}(P):= \{\phi\in C^{0}({\mathbb R}, H^{1}({\mathbb R}^{d}))\cap C^{1}({\mathbb R}, L^{2}({\mathbb R}^{d}))\ :\ P(t, {\rm x}, D_{{\rm x}})\phi=0\}.
\]
It is well-known that  $\phi\in{\rm Sol}_{E}(P)$ iff $f=(\phi(0), {\rm i}^{-1} \partial_{t}\phi(0))\in H^{1}({\mathbb R}^{d})\oplus L^{2}({\mathbb R}^{d})$.
 Moreover the symplectic form $\sigma$ is well defined in ${\rm Sol}_{E}(P)$.

Recall that for  $r\in {\mathcal M}$ one sets $r_{+}= r, \ r_{-}= r^{*}$. We define now
\[
C^{\pm}(r):= \{f\in H^{1}({\mathbb R}^{d})\oplus L^{2}({\mathbb R}^{d})\ : \ f_{0}\mp r_{\mp}f_{1}=0\},
\]
and
\[
{\rm Sol}_{E}^{\pm}(P,r):= \{\phi\in {\rm Sol}_{E}(P)\ : \ (\phi(0); {\rm i}^{-1} \partial_{t}\phi(0))\in C^{\pm}(r)\}.
\]
We call ${\rm Sol}_{E}^{\pm}(P,r)$ the space of {\em  positive/negative wavefront set solutions}.
\begin{theoreme}\label{posit}
 Let $r\in{\mathcal M}$. Then the following properties hold:
 \ben
 \item ${\rm Sol}_{E}(P)= {\rm Sol}^{+}_{E}(P,r)\oplus{\rm Sol}^{-}_{E}(P,r)$,
 \item $\phi\in {\rm Sol}^{\pm}_{E}(P,r)\ \Rightarrow \ {\rm WF}(\phi)\subset {\mathcal N}_{\pm}$,
 \item $\pm{\rm i} \sigma= \pm q>0$ on ${\rm Sol}^{\pm}_{E}(P,r)$,
 \item the spaces $ {\rm Sol}^{\pm}_{E}(P,r)$ are symplectically orthogonal.
   \een
\end{theoreme}
\begin{remark}
 We can interpret Thm. \ref{posit} as follows:  the space of finite energy solutions decomposes as the direct sum of the spaces of positive resp. negative wavefront set solutions. The charge $q$ is positive, resp. negative on these spaces. Moreover these two spaces are symplectically orthogonal.
 This is the exact analogue of the well-known situation in the static case,
where $a(t, {\rm x}, D_{{\rm x}})$ does not depend on $t$ (cf. Subsect. \ref{for-record}).
\end{remark}
For the proof of Thm. \ref{posit}, we will use the following lemma.
\begin{lemma}\label{4.1}
Let $r\in {\mathcal M}$. Then:
 \ben
 \item $r+ r^{*}: L^{2}({\mathbb R}^{d})\to H^{1}({\mathbb R}^{d})$ is invertible  and
 \[
 (r+ r^{*})^{-\frac{1}{2}}= \frac{1}{\sqrt{2}} \epsilon^{\frac{1}{2}}+ \Psi_{\rm ph}^{0}({\mathbb R}^{d}).
 \]
 \item The operator
 \[
T(r) := (r+r^{*})^{-\frac{1}{2}}\mat{\bbbone}{r}{\bbbone}{-r^{*}}: \  {\mathcal H}(\rr^{d})\otimes \cc^{2}\to  {\mathcal H}(\rr^{d})\otimes \cc^{2},
\]
is bounded with bounded inverse:
\[
T(r)^{-1}= \mat { r^{*}}{r}{\bbbone}{-\bbbone}(r+r^{*})^{-\frac{1}{2}}.
\]
\item $T(r): H^{1}({\mathbb R}^{d})\oplus L^{2}({\mathbb R}^{d})\to H^{\frac{1}{2}}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$ is bounded with inverse $T(r)^{-1}$.
\item one has:
\[
\tilde{q}:= (T(r)^{-1})^{*}\circ q\circ T(r)^{-1}= \mat{\bbbone}{0}{0}{-\bbbone}.
\]
 \een
\end{lemma}
\proof
From Thm.  \ref{3.1}  we obtain that $r+ r^{*}\sim  \epsilon^{-1}$, which implies that $r+r^{*}$ is bijective from $L^{2}({\mathbb R}^{d})$ to $H^{1}({\mathbb R}^{d})$. Moreover $(r+r^{*})^{-1}= {\frac{1}{2}} \epsilon+ \Psi^{0}$ and  $(r+r^{*})^{-1}\sim  \epsilon$. By Prop. \ref{1.1} we obtain (1). Statements (2), (3), (4) follow from (1). $\Box$\medskip

\medskip

{\bf Proof of Thm. \ref{posit}.}
For $\tilde{f}\in H^{\frac{1}{2}}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$ we set $\tilde{f}= (\tilde{f}_{+}, \tilde{f}_{-})$.  Since $r_{+}= r$, $r_{-}= r^{*}$ we obtain that by Lemma \ref{4.1} $f\in C^{\pm}$ iff $(Tf)_{\mp}=0$.  The theorem follows then from Lemma \ref{4.1} (4).   $\Box$\medskip
\section{Construction of Hadamard states}\label{sec4}\setcounter{equation}{0}
\subsection{Microlocal spectrum condition}\label{sec4.2}

In this subsection we discuss conditions under which a covariance $c$ on ${\mathcal D}(\rr^{d})\otimes \cc^{2}$ leads by(\ref{defino}) to a covariance $C$ on ${\mathcal D}(M)$ satisfying the microlocal spectrum condition in Def. \ref{def2}.

Recall that  $c:E_{1}\to E_{2}$ means that $c$ is linear continuous between the two topological vector spaces $E_{1}$ and $E_{2}$.

We consider the model Klein-Gordon equation:
\[
 \partial_{t}^{2}\phi + a(t, {\rm x}, D_{{\rm x}})\phi=0,
\]
introduced in Subsect. \ref{sec2.0b}.

Let $c$  be a  bounded hermitian form on  ${\mathcal D}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$. We identify it with the operator:
\begin{equation}\label{e4.001}
c= \mat{c_{00}}{c_{01}}{c_{10}}{c_{11}}: \ {\mathcal D}({\mathbb R}^{d})\otimes {\mathbb C}^{2}\to {\mathcal D}'({\mathbb R}^{d})\otimes {\mathbb C}^{2},
\end{equation}
and associate to it the  bounded hermitian form  $C$ on ${\mathcal D}({\mathbb R}^{1+d})$ given by:
\begin{equation}\label{e4.01}
(u|Cv)= (\rho\circ E u|c \rho\circ Ev), \ C: {\mathcal D}({\mathbb R}^{1+d})\to {\mathcal D}'({\mathbb R}^{1+d}).
\end{equation}
We still denote  by $C\in {\mathcal D}'({\mathbb R}^{1+d}\times {\mathbb R}^{1+d})$ the distribution kernel of $C$ given by:
\[
(u|Cv)= \int C(x, y)\overline{u}(x)v(y)dxdy.
\]

We fix now an operator $r\in {\mathcal M}$ (see Def. \ref{defo}). The map $T(r)$ in Lemma \ref{4.1} will be denoted by $T$ for simplicity.

We would like to define:
\begin{equation}\label{e4.05}
\tilde{c}:= (T^{-1 })^{*}\circ c \circ T^{-1}=: \mat{\tilde{c}_{++}}{\tilde{c}_{+-}}{\tilde{c}_{-+}}{\tilde{c}_{--}}.
\end{equation}
Since $T$, $T^{-1}: {\mathcal H}(\rr^{d})\otimes \cc^{2}\to {\mathcal H}'(\rr^{d})\otimes \cc^{2}$, a natural requirement is that
\[
c: {\mathcal H}(\rr^{d})\otimes \cc^{2}\to {\mathcal H}'(\rr^{d})\otimes \cc^{2},
\]
which implies that $c: {\mathcal D}(\rr^{d})\otimes \cc^{2}\to {\mathcal D}'(\rr^{d})\otimes \cc^{2}$.  In the next theorem we will need to impose stronger conditions on $c$.

We can now state the main result of this subsection.  Recall that the notations $ _{M}\!\Gamma$ and $\Gamma_{M}$  for a conic set $\Gamma$ are defined in Subsect. \ref{sec-1}.
\begin{theoreme}\label{4.2}
Assume that
\[
 \begin{array}{l}
(1a)\ c: {\mathcal H}(\rr^{d})\otimes \cc^{2}\to {\mathcal H}(\rr^{d})\otimes \cc^{2}, \\[2mm]
(1b)\ \ _{{\mathbb R}^{d}}\!{\rm WF}(c)'= {\rm WF}(c)'_{{\mathbb R}^{d}}= \emptyset.
\end{array}
\]
Let $C$ be defined by (\ref{e4.01})  and
\[
 \Lambda_{+}:= C, \Lambda_{-}:= C- {\rm i} E.
\]
Then $\Lambda_{\pm}$ satisfies the microlocal spectrum condition iff:
\[
(2)\ {\rm WF}(\tilde{c}_{--})'= {\rm WF}(\tilde{c}_{+-})'= {\rm WF}(\tilde{c}_{-+})'=  {\rm WF}(\bbbone -\tilde{c}_{++})'= \emptyset.
\]
\end{theoreme}
\begin{remark}\label{remaro}
 Note that condition (2) implies that condition (1b) is satisfied by $\tilde{c}$.
 Using that $T$, $T^{-1}$ are pseudo-differential  operators, it is easy to see that condition (1b) is then also satisfied by $c$. Therefore if conditions (1a), (2) are satisfied, $\Lambda_{\pm}$ satisfies $(\mu{\rm sc})$.
\end{remark}
\begin{remark}
Note  that we strengthen the assumption on the sesquilinear form $c$, since we require in (1a) that $c: {\mathcal H}(\rr^{d})\otimes \cc^{2}\to {\mathcal H}(\rr^{d})\otimes \cc^{2}$ instead of  $c: {\mathcal H}(\rr^{d})\otimes \cc^{2}\to {\mathcal H}'(\rr^{d})\otimes \cc^{2}$ as before.  In fact since the Cauchy surface is not compact, some care is needed with the composition  of operators.

Condition (1b) is satisfied for example if ${\rm WF}(c)'\subset \Gamma$, where $\Gamma$ is the graph of a conic, bijective map on $T^{*}{\mathbb R}^{d}$.
This is the case if the entries of $c$ are pseudo-differential  or even Fourier integral operators.
\end{remark}
\proof

We set  \[
\tilde{\rho}= T\circ \rho=: \tilde\rho_{+}\oplus \tilde\rho_{-},
\]so that:
\begin{equation}\label{e4.04}
C= (\tilde{\rho}\circ E)^{*}\circ \tilde{c}\circ (\tilde{\rho}\circ E)=\sum_{\alpha, \beta\in\{+,-\}}C_{\alpha \beta},
\end{equation}
for
\begin{equation}\label{e4.04b}
 \ C_{\alpha \beta}:= (\tilde\rho_{\alpha}\circ E)^{*}\circ \tilde{c}_{\alpha\beta}\circ (\tilde\rho_{\beta}\circ E).
\end{equation}
Let us first check that we can perform the various compositions in (\ref{e4.04}).

Because of the well-known support properties of $E_{\pm}$  we have:
\begin{equation}\label{shiti}
\begin{array}{rl}
\rho\circ E:&\begin{array}{l}
{\mathcal D}(M)\to {\mathcal D}({\mathbb R}^{d})\otimes {\mathbb C}^{2},\\[2mm]
{\mathcal E}'(M)\to {\mathcal E}'(\rr^{d})\otimes \cc^{2},
\end{array}\\[4mm]
(\rho\circ E)^{*}:&\begin{array}{l}
{\mathcal D}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}\to {\mathcal D}'(M),\\[2mm]
{\mathcal E}(\rr^{d})\otimes \cc^{2}\to {\mathcal E}(M).
\end{array}
\end{array}
\end{equation}
Note also that  $T: {\mathcal D}({\mathbb R}^{d})\otimes {\mathbb C}^{2}\to {\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$ and $T: {\mathcal E}'(\rr^{d})\otimes \cc^{2}\to {\mathcal H}'(\rr^{d})\otimes \cc^{2}$. We obtain that
\begin{equation}\label{shito}
\begin{array}{rl}
\tilde\rho\circ E:&\begin{array}{l}
{\mathcal D}(M)\to {\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2},\\[2mm]
{\mathcal E}'(M)\to {\mathcal H}'(\rr^{d})\otimes \cc^{2},
\end{array}\\[4mm]
(\tilde\rho\circ E)^{*}:&\begin{array}{l}
{\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}\to {\mathcal D}'(M),\\[2mm]
{\mathcal H}(\rr^{d})\otimes \cc^{2}\to {\mathcal E}(M).
\end{array}
\end{array}
\end{equation}

Since we assumed that $c:{\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}\to {\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}$, we have $\tilde{c}: {\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}\to {\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}$, using that $T$, $T^{-1}$ preserve ${\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$ and ${\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}$. Therefore we can perform the compositions in (\ref{e4.04}).

By Lemma \ref{lem:raff}, we have
\[
{\rm WF}(E)'=\{(X_{1}, X_{2}) : \ X_{1}\sim X_{2}, \ X_{1}, X_{2}\in {\mathcal N}\},
\]
and using that the Cauchy surface $\{t=0\}$ is  non-characteristic for the Klein-Gordon equation, we have for $i=0,1$:
\[
\begin{array}{rl}
&{\rm WF}(\rho_{i}\circ E)'\\[2mm]
=& \{\left(({\rm x}_{1}, k_{1}), (x_{2}, \xi_{2})\right)\in \coo{({\mathbb R}^{d}\times M)} : \ (0, {\rm x}_{1}, - \epsilon_{1}(0, {\rm x}_{1}, k_{1}), k_{1})\sim (x_{2}, \xi_{2})\} \\[2mm]
 \cup & \{\left(({\rm x}_{1}, k_{1}), (x_{2}, \xi_{2})\right)\in \coo{({\mathbb R}^{d}\times M)} : \ (0, {\rm x}_{1}, + \epsilon_{1}(0, {\rm x}_{1}, k_{1}), k_{1})\sim (x_{2}, \xi_{2})\}.
\end{array}
\]
Then from Corollary \ref{3.1b}, we obtain that:
\begin{equation}
\label{e4.1}
{\rm WF}(\tilde\rho_{\pm}\circ E)'= \Gamma_{\pm},
\end{equation}
for
\[
\Gamma_{\pm}=\{\left(({\rm x}_{1}, k_{1}), (x_{2}, \xi_{2})\right)\in \coo{({\mathbb R}^{d}\times M)} : \ (0, {\rm x}_{1}, \pm \epsilon(0, {\rm x}_{1}, k_{1}), k_{1})\sim (x_{2}, \xi_{2})\}.
\]
We also have
\begin{equation}\label{e4.1bis}
{\rm WF}((\tilde\rho_{\pm}\circ E)^{*})'= {\rm Exch}(\Gamma_{\pm}).
\end{equation}
We now want to apply the composition rule for the wave front set recalled in Subsect. \ref{sec-1} to the identity (\ref{e4.04b}), in order to bound ${\rm WF}(C_{\alpha\beta})'$.  It clearly suffices to bound ${\rm WF}(\chi C_{\alpha\beta}\chi)'$ for $\chi\in C_0^\infty(M)$.

\medskip

{\em Step 1.}
The first step is to reduce oneself  to the composition of properly supported kernels, modulo smoothing operators.

Because of the support properties of the kernel of $E$, there exists $\psi\in C_0^\infty({\mathbb R}^{d})$ such that (denoting again $\mat{\psi}{0}{0}{\psi}$ by $\psi$):
\[
 \rho\circ E \chi= \psi \rho\circ E \chi.
\]
Let us also fix $\psi_{1}\in C_0^\infty({\mathbb R}^{d})$ with $\psi_{1}\equiv 1$ near ${\rm supp} \psi$.  By Lemma \ref{idiotic} we know that $ (1- \psi_{1})T\psi\in \Psi^{-\infty}({\mathbb R}^{d})$, hence
\[
\tilde\rho\circ E \chi= \psi_{1}\tilde\rho\circ E \chi+ R_{1}\rho\circ E \chi,
\]
for
\begin{equation}\label{shitu}
R_{1}:= (1- \psi_{1})T \psi: {\mathcal D}'(\rr^{d})\otimes \cc^{2}\to {\mathcal H}(\rr^{d})\otimes \cc^{2}.
\end{equation}
Taking adjoints we have:
\[
\begin{array}{l}
\chi(\tilde\rho\circ E)^{*}= (\tilde \rho \circ E \chi)^{*}= \chi(\tilde \rho\circ E)^{*}\psi_{1}+ \chi (\rho\circ E)^{*}R_{1}^{*},\\[2mm]
R_{1}^{*}: {\mathcal H}'(\rr^{d})\otimes \cc^{2} \to {\mathcal D}(\rr^{d})\otimes \cc^{2}.
\end{array}
\]
It follows that
\[
\begin{array}{rl}
& \chi(\tilde \rho\circ E)^{*}\tilde{c}(\tilde \rho\circ E)\chi\\[2mm]
=&  \chi(\tilde \rho\circ E)^{*}\psi_{1}\tilde{c}\psi_{1}(\tilde \rho\circ E)\chi
+  \chi(\tilde \rho\circ E)^{*}\tilde{c}R_{1}(\rho\circ E)\chi+  \chi(\tilde \rho\circ E)^{*}R_{1}^{*}\tilde{c}\psi_{1}(\tilde \rho\circ E)\chi\\[2mm]
=:& \chi(\tilde \rho\circ E)^{*}\psi_{1}\tilde{c}\psi_{1}(\tilde \rho\circ E)\chi+ I_{1}+ I_{2}.
\end{array}
\]
We claim that
\begin{equation}\label{sta}
I_{1}, \ I_{2}: {\mathcal D}'(M)\to {\mathcal D}(M).
\end{equation}
Note that from hypothesis (1a) and the fact that $T$, $T^{-1}$ are (matrices of) pseudo-differential  operators, we know that $\tilde{c}: {\mathcal H}(\rr^{d})\otimes \cc^{2}\to {\mathcal H}(\rr^{d})\otimes \cc^{2}$.

Using then (\ref{shiti}), (\ref{shito}), (\ref{shitu}) we have:
\[
\begin{array}{rl}
I_{1}:&{\mathcal D}'(M):\mathop{\to}\limits^{(\rho\circ E)\chi} {\mathcal E}'(\rr^{d})\otimes \cc^{2}\subset {\mathcal H}'(\rr^{d})\otimes \cc^{2}\mathop{\to}\limits^{R_{1}} {\mathcal H}(\rr^{d})\otimes \cc^{2}\\[2mm]
 &\mathop{\to}\limits^{\tilde{c}}{\mathcal H}(\rr^{d})\otimes \cc^{2}\subset {\mathcal E}(\rr^{d})\otimes \cc^{2}\mathop{\to}\limits^{\chi(\tilde{\rho}\circ E)^{*}} {\mathcal D}(M),
\end{array}
\]
and:
\[
\begin{array}{rl}
I_{2}:&{\mathcal D}'(M):\mathop{\to}\limits^{(\tilde\rho\circ E)\chi} {\mathcal H}'(\rr^{d})\otimes \cc^{2}\mathop{\to}\limits^{\tilde{c}\psi_{1}} {\mathcal H}'(\rr^{d})\otimes \cc^{2}\\[2mm]
 &\mathop{\to}\limits^{R_{1}^{*}}{\mathcal D}(\rr^{d})\otimes \cc^{2}\subset {\mathcal E}(\rr^{d})\otimes \cc^{2}\mathop{\to}\limits^{\chi(\rho\circ E)^{*}} {\mathcal D}(M),
\end{array}
\]
which proves (\ref{sta}).  It follows that if $\psi_{2}\in C_0^\infty({\mathbb R}^{d})$ with $\psi_{2}\equiv 1$ near ${\rm supp} \psi_{1}$ we have:
\[
\chi C_{\alpha\beta} \chi=   \chi (\tilde \rho_{\alpha}\circ E)^{*} \psi_{1}\circ \psi_{2}\tilde{c}\psi_{2}\circ \psi_{1}(\tilde{\rho}_{\beta}\circ E)\chi\hbox{ mod }C^{\infty}(M\times M),
\]
the three operators in the composition above having compactly (hence properly) supported kernels.

{\em Step 2.} We check that we can apply the composition rule for wave front sets.  Note first that since $T$, $T^{-1}$ are (matrices of) pseudo-differential  operators, we obtain using hypothesis (1b) that:
\begin{equation}
\label{wavefr}
_{{\mathbb R}^{d}}\!{\rm WF}(\tilde{c})'= {\rm WF}(\tilde{c})'_{{\mathbb R}^{d}}= \emptyset.
\end{equation}
Let us fix $\alpha, \beta\in \{+,-\}$ and set:
 \[
K_{\alpha}=  \chi (\tilde \rho_{\alpha}\circ E)^{*} \psi_{1},\ K_{\alpha\beta}= \psi_{2}\tilde{c}_{\alpha\beta}\psi_{2}, \ K_{\beta}=  \psi_{1}(\tilde{\rho}_{\beta}\circ E)\chi.
\]  Using (\ref{e4.1}), (\ref{e4.1bis}) we have:
\[
{\rm WF}(K_{\alpha})'\subset {\rm Exch}(\Gamma_{\alpha}), \ {\rm WF}(K_{\alpha\beta})'\subset {\rm WF}(\tilde{c}_{\alpha\beta})', \ {\rm WF}(K_{\beta})'\subset \Gamma_{\beta}.
\]
It follows also from (\ref{wavefr}) that:
\[\begin{array}{rl}
&_{M}\!{\rm WF}(K_{\alpha})'= {\rm WF}(K_{\alpha})'_{{\mathbb R}^{d}}= _{{\mathbb R}^{d}}\!{\rm WF}(K_{\alpha\beta})'\\[2mm]
&= {\rm WF}(K_{\alpha\beta})'_{{\mathbb R}^{d}}= _{{\mathbb R}^{d}}\!{\rm WF}(K_{\alpha\beta})'= {\rm WF}(K_{\beta})'_{M}= \emptyset.
\end{array}
\]
It follows that we can compose $K_{\alpha\beta}$ and $K_{\beta}$ and
\[
{\rm WF}(K_{\alpha\beta}\circ K_{\beta})'\subset {\rm WF}(\tilde{c}_{\alpha\beta})\circ \Gamma_{\beta}.
\]
We can also compose $K_{\alpha}$ and $K_{\alpha\beta}\circ K_{\beta}$ and
\[
{\rm WF}(K_{\alpha}\circ K_{\alpha\beta}\circ K_{\beta})\subset  {\rm Exch}(\Gamma_{\alpha})\circ {\rm WF}(\tilde{c}_{\alpha\beta})\circ \Gamma_{\beta}.
\]
 {\em Step 3.}
 Recalling the definition of $\Gamma_{\alpha}$, $\Gamma_{\beta}$, we obtain from Step 2  that:
\begin{equation}
\label{e4.2}
\begin{array}{rl}
& {\rm WF}(C_{\alpha\beta})'\\[2mm]
\subset&\{(x_{1}, \xi_{1}, x_{2}, \xi_{2})\ : \ (x_{1}, \xi_{1})\in {\mathcal N}_{\alpha}, \ (x_{2}, \xi_{2})\in {\mathcal N}_{\beta}, \ \exists \ ({\rm x}_{1}, k_{1}, {\rm x}_{2}, k_{2})\in {\rm WF}(\tilde{c}_{\alpha\beta})'\\[2mm]
&\hbox{such that }(x_{1}, \xi_{1})\sim (0,{\rm x}_{1}, \alpha\epsilon(0, {\rm x}_{1}, k_{1}),k_{1}), \ (x_{2}, \xi_{2})\sim (0,{\rm x}_{2}, \beta\epsilon(0, {\rm x}_{2}, k_{2}))\}.
\end{array}\end{equation}
Let $S_{\alpha\beta}$ be the set in the r.h.s. of (\ref{e4.2}).  Using (\ref{e4.04}) and the fact that the $S_{\alpha\beta}$ are pairwise disjoint, we obtain that:
\begin{equation}\label{astu}
{\rm WF}(C)'\subset  \bigcup_{\alpha, \beta \in \{+,-\}}S_{\alpha\beta}.
\end{equation}
{\em Step 4.}
Recall that we set $\Lambda_{+}= C$, $\Lambda_{-}= C- {\rm i} E$.  We first consider the condition
\[
(\mu{\rm sc}_{+}):\  {\rm WF}(\Lambda_{+})'\subset \{(X_{1}, X_{2})\ : \ (X_{i})\in {\mathcal N}_{+}, \ X_{1}\sim X_{2}\}.
\]
By (\ref{astu}) $(\mu{\rm sc}_{+})$ is satisfied iff:
\[
\begin{array}{rl}
S_{\alpha\beta}= &\emptyset,\hbox{ for }(\alpha, \beta)\neq (+,+),\\[2mm]
S_{++}\subset &\{(x_{1}, \xi_{1}, x_{2}, \xi_{2})\ : \ (x_{i}, \xi_{i})\in {\mathcal N}_{+}, \ (x_{1}, \xi_{1})\sim (x_{2}, \xi_{2})\}.
\end{array}
\]
This condition is satisfied iff
\begin{equation}\label{plus}
{\rm WF}(\tilde{c}_{\alpha\beta})'=\emptyset \hbox{ for }(\alpha, \beta)\neq (+, +), \ {\rm WF}(\tilde{c}_{++})'\subset \Delta,
\end{equation}
where $\Delta$ is the diagonal in $\cooo{{\mathbb R}^{d}}\times \cooo{{\mathbb R}^{d}}$.
Let us now consider
\[
(\mu{\rm sc}_{-}):\  {\rm WF}(\Lambda_{-})'\subset \{(X_{1}, X_{2})\ : \ (X_{i})\in {\mathcal N}_{-}, \ X_{1}\sim X_{2}\}.
\]
By (\ref{e0.2}), replacing $C$ by $C-{\rm i} E$ amounts to replace  $\tilde{c}$ by $\tilde{c}- \tilde{q}$.  Therefore  $(\mu{\rm sc}_{-})$ is satisfied iff:
\begin{equation}\label{minus}
{\rm WF}(\tilde{c}_{\alpha\beta}- \delta_{\alpha\beta}\bbbone)'=\emptyset \hbox{ for }(\alpha, \beta)\neq (-, -), \ {\rm WF}(\tilde{c}_{--}+\bbbone)'\subset \Delta.
\end{equation}
 Combining (\ref{plus}) and (\ref{minus}) we obtain that $\Lambda_{\pm}$ satisfy $(\mu{\rm sc})$ iff
 \[
{\rm WF}(\tilde{c}_{--})'= {\rm WF}(\tilde{c}_{+-})'= {\rm WF}(\tilde{c}_{-+})'=  {\rm WF}(\bbbone -\tilde{c}_{++})'= \emptyset,
\]
which completes the proof of the theorem.   $\Box$\medskip
\subsection{Construction of Hadamard states}\label{sec4.3}
In this subsection we construct a large class of two-point functions $\lambda$ with pseudo-differential  entries, such that $\Lambda$ is  the two-point function  of a (gauge-invariant) quasi-free Hadamard state.
Beside the microlocal condition in Thm. \ref{4.2}, $\lambda$ should also satisfy the positivity conditions recalled in Subsect. \ref{sec0.1}, i.e. $\lambda\geq 0$, $\lambda\geq q$, where $q={\rm i}\sigma=\left(\begin{smallmatrix}0 & \bbbone \\ \bbbone & 0 \end{smallmatrix}\right)$.

As before we fix $r\in {\mathcal M}$.
\begin{proposition}\label{4.4}
Let $\lambda$ be a two-point function with pseudo-differential  entries.  Let $\tilde{\lambda}_{\alpha \beta}$ for $\alpha, \beta\in \{+,-\}$ be defined as in (\ref{e4.05}).
Then $\lambda$ is a Hadamard charge density iff
\[
\begin{array}{rl}
&(\mu{\rm sc}')\ \tilde \lambda_{--}, \ \tilde \lambda_{+-}, \tilde \lambda_{-+}\in \Psi^{-\infty}({\mathbb R}^{d}),\\[2mm]
&(1)\ \tilde \lambda_{++}\geq \bbbone, \hbox{ on }{\mathcal H}({\mathbb R}^{d}), \ \tilde \lambda_{--}\geq 0\hbox{ on }{\mathcal H}({\mathbb R}^{d}),\\[2mm]
&(2)\ |(u| \tilde \lambda_{+-}v)|\leq (u|\tilde  \lambda_{++}u)^{\frac{1}{2}}(v|\tilde  \lambda_{--}v)^{\frac{1}{2}}, \ u, v\in {\mathcal H}({\mathbb R}^{d}),\\[2mm]
&(3)\ |(u| \tilde\lambda_{+-}v)|\leq (u| (\tilde\lambda_{++} -\bbbone)u)^{\frac{1}{2}}(v|(\tilde\lambda_{--}+\bbbone)v)^{\frac{1}{2}}, \ u, v\in {\mathcal H}({\mathbb R}^{d}).
\end{array}
\]
\end{proposition}
\proof
Since the entries of $\lambda$, $\tilde \lambda$ are pseudo-differential  operators, condition (1a) of Thm. \ref{4.2} is satisfied and condition (1b) as well by Remark \ref{remaro}. Moreover, condition $(\mu{\rm sc}')$ is equivalent to (2) of Thm. \ref{4.2}, hence $(\mu{\rm sc}')$ is equivalent to the microlocal spectrum condition.

From Sect. \ref{sec0} we know that  $\lambda$ is the two-point function of a gauge-invariant quasi-free state iff
\begin{equation}\label{e4.7}
 \lambda\geq 0, \ \lambda\geq q \ \ \hbox{on }{\mathcal D}(\rr^{d})\otimes \cc^{2}
\Leftrightarrow  \lambda\geq 0, \ \lambda\geq q \ \  \hbox{on }{\mathcal H}(\rr^{d})\otimes \cc^{2},
\end{equation}
using that  the entries of $\lambda$ are  pseudo-differential operators.

We recall that  $\tilde{\lambda}= (T^{-1})^{*}\circ \lambda \circ T^{-1}$ and  $\tilde{q}= (T^{-1})^{*}\circ q \circ T^{-1}$.   By Lemma \ref{4.1} we have
\[
\tilde{q}= \mat{\bbbone}{0}{0}{ -\bbbone}: \ {\mathcal H}(\rr^{d})\otimes \cc^{2} \to {\mathcal H}'(\rr^{d})\otimes \cc^{2}.
\]

Since $T$ maps ${\mathcal H}(\rr^{d})\otimes \cc^{2}$ into itself bijectively, (\ref{e4.7}) is equivalent to:
\begin{equation}
\label{e4.3}
\tilde{\lambda}\geq 0, \  \tilde{\lambda}\geq \tilde{q} \ \ \hbox{on }{\mathcal H}(\rr^{d})\otimes \cc^{2}.
\end{equation}
Clearly if $a, b, c$ are linear operators with domain ${\mathcal H}({\mathbb R}^{d})$ one has:
\[
\begin{array}{rl}
&2{\rm Re}(u| bv)+ (u|au)+ (v| cv)\geq 0, \ u, v\in {\mathcal H}({\mathbb R}^{d})\\[2mm]
\Leftrightarrow& |(u| bv)|\leq (u| au)^{\frac{1}{2}}(v|cv)^{\frac{1}{2}}, \ u, v\in {\mathcal H}({\mathbb R}^{d})\hbox{ and }a, c\geq 0\hbox{ on }{\mathcal H}({\mathbb R}^{d}).
\end{array}
\]
Applying this observation and noting that $r+r^{*}\geq 0$, we obtain that condition (\ref{e4.3}) is equivalent to conditions (1), (2), (3). $\Box$\medskip

We now proceed to construct a large class of pseudo-differential  operators $\tilde\lambda_{\alpha\beta}$ satisfying the conditions in Prop. \ref{4.4}.

\begin{theoreme}\label{4.3}
Let us fix pseudo-differential  operators:
\[
a_{-\infty}, b_{-\infty}\in \Psi^{-\infty}({\mathbb R}^{d}), \ a_{0}\in \Psi^{0}({\mathbb R}^{d}) \hbox{ with }\|a_{0}\|\leq 1,
\]
and set:
\[
\begin{array}{rl}
\tilde{\lambda}_{++}=&\bbbone+ b_{-\infty}^{*}b_{-\infty},\\[2mm]
\tilde \lambda_{--}= &a_{-\infty}^{*}a_{-\infty},\\[2mm]
\tilde \lambda_{+-}=&\tilde \lambda_{-+}^{*}= b_{-\infty}^{*}a_{0}a_{-\infty}
.\end{array}
\]
Then the two-point function $\lambda$ given by (\ref{e4.05}) is the two-point function of a Hadamard state.
\end{theoreme}

\proof  We  check the conditions in Prop. \ref{4.4}.
Conditions $(\mu{\rm sc})$ and (1) are clearly satisfied.   From the form of $\tilde{\lambda}_{+-}$ we have
\[
|(u| \tilde{\lambda}_{+-}v)|\leq  (u| b_{-\infty}^{*}b_{-\infty}u)^{\frac{1}{2}}(v| a_{-\infty}^{*}a_{-\infty}v)^{\frac{1}{2}}, \ u,v\in {\mathcal H}({\mathbb R}^{d}),
\]
which implies (2) and (3), using the form of $\tilde{\lambda}_{++}$ and $\tilde{\lambda}_{--}$.  $\Box$\medskip

\subsection{Symplectic transformations}\label{sec4.4}
Recall from Sect. \ref{sec0} that if $({\mathcal Y}, \sigma)$ is a complex symplectic space and $q= {\rm i} \sigma$, then the set of two-point functions of gauge-invariant quasi-free states is invariant under conjugation by elements of $U({\mathcal Y}, q)$. The same is true for the set of two-point functions of pure quasi-free states.

In this subsection we describe a class of symplectic transformations $u\in U({\mathcal D}({\mathbb R}^{d})\otimes {\mathbb C}^{2},q)$ which preserve the  microlocal spectrum condition $(\mu{\rm sc})$. We start with a general result.

\begin{proposition}\label{p4.4}
 Let
 $u$ such that $u, u^{*}: {\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}\to {\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$.
 Set
 \[
\tilde{u}:= TuT^{-1}=\mat{\tilde{u}_{++}}{\tilde{u}_{+-}}{\tilde{u}_{-+}}{\tilde{u}_{--}},
\]
and assume that
\[
{\rm WF}(\tilde{u}_{++}^{*} \tilde{\lambda}_{++}\tilde{u}_{++})'\subset \Delta, \ \tilde{u}_{+-}: {\mathcal H}'({\mathbb R}^{d})\to {\mathcal H}({\mathbb R}^{d}).
\]
Then if $\lambda$ is a two-point function with pseudo-differential  entries satisfying $(\mu{\rm sc})$,  the two-point function
$u^{*}\lambda u$ satisfies also $(\mu{\rm sc})$.
\end{proposition}
\proof
We set  $ c:= u^{*}\lambda u$ and as in Subsect. \ref{sec4.2}:
\[
\tilde{\lambda}:= (T^{-1})^{*} \lambda T^{-1}, \ \tilde{c}:= (T^{-1})^{*} c T^{-1}= \tilde{u}^{*} \tilde{\lambda}\tilde{u}.
\]
Since $\lambda$ has pseudo-differential  entries and satisfies $(\mu{\rm sc})$ we have:
\begin{equation}\label{e.buda7}
\tilde{\lambda}_{\alpha\beta}\in \Psi^{\infty}({\mathbb R}^{d}), \ \tilde{\lambda}_{\alpha\beta}\in \Psi^{-\infty}({\mathbb R}^{d}) \ \hbox{ for }(\alpha, \beta)\neq (+,+).
\end{equation}

We will check that $c$ satisfies the hypotheses of Thm. \ref{4.2}. Since $u, u^{*}, \lambda$ preserve ${\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$ condition (1a) is satisfied. By Remark \ref{remaro}, it remains to check condition (2). We compute $\tilde{c}$ and
obtain using (\ref{e.buda7}) that
\[
\tilde{c}=\mat{\tilde{u}_{++}^{*}\tilde{\lambda}_{++}\tilde{u}_{++}}{\tilde{u}^{*}_{++}\tilde{\lambda}_{++}\tilde{u}_{+-}}{\tilde{u}^{*}_{+-}\tilde{\lambda}_{++}\tilde{u}_{++}}{\tilde{u}^{*}_{+-}\tilde{\lambda}_{++}\tilde{u}_{+-}}+ s,
\]
 where $s: {\mathcal H}'({\mathbb R}^{d})\otimes {\mathbb C}^{2}\to {\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}$ is a smoothing operator.  Since $\tilde{u}_{+-}, \tilde{u}_{+-}:{\mathcal H}'({\mathbb R}^{d})\to {\mathcal H}({\mathbb R}^{d})$ we have
 \[
\tilde{c}=\mat{\tilde{u}_{++}^{*}\tilde{\lambda}_{++}\tilde{u}_{++}}{0}{0}{0}+ s_{1},
\]
for $s_{1}$ as $s$. Therefore condition (2) is satisfied. $\Box$\medskip
\medskip

\begin{definition}
 We denote by $U_{-\infty}({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q)$  the subgroup of $U({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q)$ defined by:
\[
U_{-\infty}({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q):=\{u\in U({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q)\ : \ u-\bbbone\in\Psi^{-\infty}({\mathbb R}^{d})\otimes M_{2}({\mathbb C})\}.
\]
\end{definition}
\begin{corollary}\label{p4.5}
The conjugations by elements of $U_{-\infty}({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q)$ preserve the set of (pure) quasi-free Hadamard states. \end{corollary}

\begin{remark}
It is easy to see that
 if $\mat{a}{b}{c}{d}\in U({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q)$ and  $a$ is invertible, then
\begin{equation}
\mat{a}{b}{c}{d}=\mat{\bbbone}{0}{e}{\bbbone}\mat{g^*}{0}{0}{g^{-1}}\mat{\bbbone}{f}{0}{\bbbone}
\end{equation}
for some $g$ invertible and $e^*=-e$, $f^*=-f$.  Moreover  the matrices
\[
(1) \ \mat{\bbbone}{0}{e}{\bbbone}, \ \ (2) \  \mat{\bbbone}{f}{0}{\bbbone} \mbox{ or } \ (3) \  \mat{g^*}{0}{0}{g^{-1}},
\]
where $e, \ f\in \Psi^{-\infty}({\mathbb R}^{d})$ with $e^*=-e$, $f^*=-f$,  and $g-\bbbone\in \Psi^{-\infty}({\mathbb R}^{d})$ with $g$, $g^{*}$  invertible, belong to $U_{-\infty}({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q)$.
\end{remark}
\

\subsection{Pure Hadamard states}\label{sec4.5}
We now characterize pure Hadamard states with pseudo-differential  entries and discuss some examples.

\begin{theoreme}\label{thm:canonicform}Let $\lambda\in\Psi^{\infty}({\mathbb R}^d)\otimes M_2({\mathbb C})$. Then $\lambda$ is the two-point function of a pure Hadamard state iff
\begin{equation}\label{eq:canonicform}
\begin{array}{rl}
\tilde{\lambda}_{++}=&\bbbone+a_{-\infty}a_{-\infty}^*,\\[2mm]
\tilde \lambda_{--}= &a_{-\infty}^{*}a_{-\infty},\\[2mm]
\tilde \lambda_{+-}=&\tilde \lambda_{-+}^{*}= a_{-\infty}(\bbbone+a_{-\infty}^*a_{-\infty})^{{\frac{1}{2}}}
\end{array}
\end{equation}
for some $a_{-\infty}\in\Psi^{-\infty}({\mathbb R}^d)$.
\end{theoreme}
\proof
Set $\tilde\eta= 2\tilde \lambda-\tilde q$. From Prop. \ref{istu} we see that $\lambda$ is the two-point function of a pure state iff
\begin{equation}
\label{turu}
i) \ \tilde\eta\geq 0, \ ii) \ \tilde \eta \tilde q^{-1} \tilde \eta=\tilde q.
\end{equation}
Writing $\tilde\eta$ as $\mat{a}{b}{b^{*}}{c}$ we obtain that (\ref{turu}) is equivalent to:
\begin{equation}
\label{toro}
\begin{array}{rl}
i')& a\geq0, \ c\geq 0, \ | (u| bv)|\leq  (u| au)^{\frac{1}{2}}(v| cv)^{\frac{1}{2}}, \ u, v\in {\mathcal H}({\mathbb R}^{d}),\\[2mm]
ii')& a^{2}= \bbbone + bb^{*}, \ c^{2}= \bbbone + b^{*}b, \ ab-bc=0.
\end{array}
\end{equation}
Note that if $b$ is a bounded operator on $L^{2}({\mathbb R}^{d})$ then:
\begin{equation}
\label{tara}
bf(b^{*}b)= f(bb^{*})b, \hbox{ for any Borel function }f.
\end{equation}
In fact (\ref{tara}) is immediate for $f(\lambda)= (\lambda-z)^{-1}$, $z\in {\mathbb C}\backslash {\mathbb R}$ and extends to all Borel functions by a standard argument.

Since $a,c\geq 0$ by {\it i')}, the first two equations of {\it ii')} yield
\[
a= (\bbbone + bb^{*})^{\frac{1}{2}}, \ c= (\bbbone + b^{*}b)^{\frac{1}{2}}.
\] The third equation of {\it ii')} then holds using (\ref{tara}). The second condition in {\it i')} is equivalent to $\| (\bbbone + bb^{*})^{\frac{1}{2}}b(\bbbone + b^{*}b)^{\frac{1}{2}}\|\leq 1$, which holds using again (\ref{tara}).

 Going back to $\tilde\lambda$ we obtain
\begin{equation}
\label{tiri}
\tilde\lambda=\frac{1}{2}\mat{(\bbbone + bb^{*})^{\frac{1}{2}}+ \bbbone}{b}{b^{*}}{(\bbbone + b^{*}b)^{\frac{1}{2}}-\bbbone}.
\end{equation}
Let now
\[
a:= \frac{b}{\sqrt{2}}((\bbbone + b^{*}b)^{\frac{1}{2}}+ \bbbone)^{\frac{1}{2}}.
\]
Using (\ref{tara}) we obtain by an easy computation that
\[
\bbbone + a^{*}a= \frac{1}{2}((\bbbone + b^{*}b)^{\frac{1}{2}}+ \bbbone), \ \bbbone + aa^{*}=  \frac{1}{2}((\bbbone + bb^{*})^{\frac{1}{2}}+ \bbbone), \ b= 2a(\bbbone + a^{*}a)^{\frac{1}{2}}.
\]
Hence $\tilde\lambda$ in (\ref{tiri}) can be rewritten as:
\begin{equation}
\label{tere}
\tilde\lambda= \mat{\bbbone + aa^{*}}{a(\bbbone + a^{*}a)^{\frac{1}{2}}}{(\bbbone + a^{*}a)^{\frac{1}{2}}a^{*}}{a^{*}a}.
\end{equation}
By Prop. \ref{4.4} $\lambda$ satisfies $(\mu{\rm sc})$ iff $a^{*}a$, $a(\bbbone + a^{*}a)^{\frac{1}{2}}\in\Psi^{-\infty}$, which is equivalent to $a\in \Psi^{-\infty}$. $\Box$\medskip

\begin{proposition}\label{prop:transinfty}
Let $\lambda_i\in\Psi^{\infty}({\mathbb R}^d)\otimes M_2({\mathbb C})$, $i=1,2$, be two-point functions of pure Hadamard states (for the model Klein-Gordon equation). Then there exists $u\in U_{-\infty}({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q)$ s.t.
\[
\lambda_2=u^* \lambda_1 u.
\]
\end{proposition}
\proof Without loss of generality we can assume that $\tilde\lambda_{1}= \mat{\bbbone}{0}{0}{0}$ and $\tilde{\lambda}_{2}$ is given by (\ref{tere}) for $a\in \Psi^{-\infty}$.  Then
\[
\tilde{u}=\mat{(\bbbone + aa^{*})^{\frac{1}{2}}}{a}{a^{*}}{(\bbbone + a^{*}a)^{\frac{1}{2}}}
\]
belongs to $U_{\infty}({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2},\tilde q)$ and satisfies $\tilde{u}^{*}\tilde{\lambda}_{1}\tilde{u}= \tilde{\lambda}_{1}$. $\Box$\medskip

\subsubsection{Canonical Hadamard state}\label{ssec:canonic}

Once having fixed $r\in{\mathcal M}$, let us consider the two-point function
\[
\lambda(r):=\mat{(r+r^{*})^{-1}}{(r+r^{*})^{-1}r}{r^{*}(r+r^{*})^{-1}}{r^{*}(r+r^{*})^{-1}r}.
\]
An easy computation shows that
\[
\tilde\lambda(r)= (T(r)^{-1})^* \lambda(r) T(r)^{-1} = \mat{\bbbone}{0}{0}{0}.
\]
This is a particular case of Theorem \ref{thm:canonicform} with $a_{-\infty}=0$ and it follows that $\lambda(r)$ is the two-point function of a pure Hadamard state. One can show that it is distinguished among all two-point functions $\lambda:{\mathcal H}(\rr^{d})\otimes \cc^{2}\to{\mathcal H}(\rr^{d})\otimes \cc^{2}$ of pure quasi-free states by the property
\[
{\rm Ran}\, P_{\pm} \subset C_{\pm}(r),
\]
where $P_{\pm}$ is defined on $H^{1}({\mathbb R}^d)\oplus L^{2}({\mathbb R}^d) $ by
\[
P_{\pm}:={\textstyle \frac{1}{2}}\bbbone\pm q\eta, \quad \eta=\lambda-{\frac{1}{2}} q.
\]

We now study the dependence of  $\lambda(r)$ on $r\in{\mathcal M}$.

\begin{proposition}\label{defu}
Let:
 \[
{\mathcal G}:=\{(g, f)\ : \ g-\bbbone, \ f\in \Psi^{-\infty}({\mathbb R}^{d}), \  g, g^{*}:L^{2}({\mathbb R}^{d})\to L^{2}({\mathbb R}^{d})\hbox{ invertible, }f=-f^{*}\}.
\]
We equip ${\mathcal G}$ with the group structure given by:
\[
\begin{array}{rl}
&{\rm Id}:= (\bbbone, 0),\\[2mm]
&G_{2}G_{1}:=(g_{2}g_{1}, (g_{2}^{*})^{-1}f_{1}g_{2}^{-1}+ f_{2}), \hbox{ for } G_{i}= (g_{i}, f_{i}).
\end{array}
\]
Then the following holds:
\ben
\item the map
\[
{\mathcal G}\ni G=(g,f)\mapsto u_{G}:=  \mat{g^*}{0}{0}{g^{-1}}\mat{\bbbone}{f}{0}{\bbbone}\in U_{-\infty}({\mathcal H}({\mathbb R}^{d})\otimes {\mathbb C}^{2}, q)
\]
is a group homomorphism.
\item
\
 ${\mathcal G}$ acts transitively on ${\mathcal M}$ by
 \[
\alpha_{G}(r):= (g^{*})^{-1} r g^{-1}+f, \ r\in {\mathcal M}, \ G= (g,f)\in {\mathcal G}.
\]
\een
\end{proposition}
\proof  Statement (1) is an easy computation. Let us prove (2).

We  first check that $\alpha_{G}$ preserves ${\mathcal M}$. Let $r\in {\mathcal M}$ and $\tilde{r}= \alpha_{G}(r)$ for $G\in {\mathcal G}$. Clearly  $\tilde{r}-r\in \Psi^{-\infty}$ so $\tilde{r}- (b^{*})^{(-1)}= r-(b^{*})^{(-1)}+ \Psi^{-\infty}\in \Psi^{-\infty}$.  It remains to check that
\begin{equation}
\label{etoto.10}
\tilde{r}+ \tilde{r}^{*}\sim  \epsilon^{-1}.
\end{equation}
This is obvious if $G= (\bbbone, f)$, since then $\tilde{r}= r+f$ and $f^{*}=-f$.

Assume now that $G= (g, 0)$, so that $\tilde{r}+ \tilde{r}^{*}= (g^{*})^{-1}(r+ r^{*})g^{-1}$, and $g-\bbbone \in \Psi^{-\infty}$, $g, g^{*}: L^{2}\to L^{2}$ invertible.
It follows that
\begin{equation}\label{etoto.11}
(\tilde{r}+ \tilde{r}^{*})^{-1}= g (r+ r^{*})^{-1}g^{*}= \frac{1}{2} \epsilon + \Psi^{0}.
\end{equation}
Since $r\in {\mathcal M}$ we have $r+ r^{*} \sim  \epsilon^{-1}$, hence $(r+ r^{*})^{-1} \sim  \epsilon$, by the Kato Heinz inequality.
In particular we have $(r+ r^{*})^{-1}\geq c_{0}>0$. Using then (\ref{etoto.11}) we obtain that:
\[
\begin{array}{rl}
(\tilde{r}+ \tilde{r}^{*})^{-1}\geq&  c_{3}>0,\\[2mm]
 (\tilde{r}+ \tilde{r}^{*})^{-1}\geq& \frac{1}{2} \epsilon- c_{4}.
\end{array}
\]
This implies that $(\tilde{r}+ \tilde{r}^{*})^{-1}\geq c \epsilon$ for some $c>0$. On the other hand (\ref{etoto.11}) directly implies that $(\tilde{r}+ \tilde{r}^{*})^{-1}\leq c \epsilon$ for some $c>0$. Therefore we have $(\tilde{r}+ \tilde{r}^{*})^{-1}\sim  \epsilon$, which implies (\ref{etoto.10}) by applying Kato-Heinz theorem once again.  This completes the proof that $\alpha_{G}$ preserves ${\mathcal M}$.

It remains to prove that the action is transitive.

 let $r_{i}\in {\mathcal M}$, $i=1,2$  As we saw above  $(r_{i}+ r_{i}^{*})^{-1}\in \Psi^{1}$  and $(r_{i}+ r_{i}^{*})^{-1}\sim  \epsilon$.  By Prop. \ref{1.1} we obtain that $(r_{i}+r_{i}^{*})^{-\frac{1}{2}}\in\Psi^{\frac{1}{2}}$ and by Kato-Heinz theorem we have $(r_{i}+r_{i}^{*})^{-\frac{1}{2}}\sim  \epsilon^{\frac{1}{2}}$. In particular $(r_{i}+r_{i}^{*})^{-\frac{1}{2}}$ is bijective from $H^{\frac{1}{2}}({\mathbb R}^{d})$ to $L^{2}({\mathbb R}^{d})$. It follows that
\begin{equation}\label{etiti.2}
(r_{2}+r_{2}^{*})^{-\frac{1}{2}}= g(r_{1}+r_{1}^{*})^{-\frac{1}{2}}= (r_{1}+r_{1}^{*})^{-\frac{1}{2}}g^{*},
\end{equation}
where $g$, $g^{*}$ are invertible on $L^{2}({\mathbb R}^{d})$. Using also that $r_{1}-r_{2}\in \Psi^{-\infty}$, we obtain that $g-\bbbone\in \Psi^{-\infty}$.
From (\ref{etiti.2}) we get:
\begin{equation}
\label{etiti.3}
r_{2}+r_{2}^{*}= (g^{*})^{-1}(r_{1}+r_{1}^{*})g^{-1}.
\end{equation}
We set now
\begin{equation}\label{etiti.4}
 r_{2}- r_{2}^{*}=: (g^{*})^{-1}(r_{1}- r_{1}^{*})g^{-1}+2f.
\end{equation}
Clearly $f^{*}= -f$, and since $g-\bbbone$ and $r_{1}-r_{2}$ belong to $\Psi^{-\infty}$, we see that $f\in \Psi^{-\infty}$. From (\ref{etiti.3}), (\ref{etiti.4}) we obtain that $r_{2}= (g^{*})^{-1}r_{1}g^{-1}+f= \alpha_{G}(r_{1})$,  for $G= (g,f)$. This completes the proof of the proposition. $\Box$\medskip

\medskip

The following theorem explains the dependence of the pure quasi-free state with two-point function $\lambda(r)$ on the choice of $r\in {\mathcal M}$.
\begin{theoreme}\label{p4.7}
We have
\[
\lambda(\alpha_{G}(r))= u_{G}^{*}\lambda(r)u_{G}, \ \forall \ r\in {\mathcal M}, \ G\in {\mathcal G}.
\]

\end{theoreme}
\proof  writing $\lambda(r)$ as:
\[
 \lambda(r)= \mat{\bbbone}{0}{r^{*}}{\bbbone}\mat{(r+r^{*})^{-1}}{0}{0}{0}\mat{\bbbone}{r}{0}{\bbbone},
\]
we easily obtain that if $u=\mat{\bbbone}{f}{0}{\bbbone}$, then
\[
u^{*}\lambda(r)u= \lambda(r+f),
\]
and if $u= \mat{g^{*}}{0}{0}{g^{-1}}$, then
\[
u^{*}\lambda(r)u= \lambda((g^{*})^{-1}rg^{-1}).
\]
This completes the proof of the theorem. $\Box$\medskip

\subsection{The static case}\label{for-record}
Let us  illustrate our results in  the {\em static case}, when $a(t, {\rm x}, D_{{\rm x}})$ is independent on $t$. We assume for simplicity that $a({\rm x}, D_{{\rm x}})\geq m^{2}>0$, in order to avoid infrared problems.   We can work in an abstract setting and denote by $a>0$ a selfadjoint operator on a (complex) Hilbert space ${\mathfrak h}$.   We set $\epsilon:= a^{\frac{1}{2}}$.

The solution of the Cauchy problem:
\begin{equation}\label{e5.1}
\left\{
\begin{array}{rl}
& \partial_{t}^{2}\phi(t)+a\phi(t)=0,\\[2mm]
&\phi(0)=f_{0},\\[2mm]
&{\rm i}^{-1} \partial_{t} \phi(0)=f_{1},
\end{array}
\right.
\end{equation}

is:
 \[
\phi(t)= {\textstyle \frac{1}{2}}{\rm e}^{{\rm i} t \epsilon}(f_{0}+ \epsilon^{-1}f_{1})+ {\textstyle \frac{1}{2}}{\rm e}^{-{\rm i} t \epsilon}( f_{0}- \epsilon^{-1}f_{1})=:U(t,0)f.
\]
Therefore, when ${\mathfrak h}=L^2({\mathbb R}^d)$ and $\epsilon\in\Psi^{1}({\mathbb R}^d)$ we can choose
\[
b(t)= \epsilon,  \ u_{\pm}(t,s)= {\rm e}^{\pm {\rm i} (t-s)\epsilon}, \ d_{\pm}(s)= \frac{1}{2}\bbbone, \ r_{\pm}(s)= \epsilon^{-1}.
\]
\begin{remark}Using the reduction to the model case described in Subsect. \ref{ss:modelcase}, one obtains $a(t, {\rm x}, D_{{\rm x}})$ independent on $t$ if the metric is static and the electric field vanishes, i.e. $\partial_i V+\partial_t A_i\equiv 0$, $i=1,\dots, d$.
\end{remark}

For sake of completeness we list below the essential examples of Hadamard states in the static case.
\begin{itemize}
\item The two-point function of the vacuum state is:
\[
(f|\lambda_{\rm vac}f)= \frac{1}{2}(f_{0}+ \epsilon^{-1}f_{1}|\epsilon(f_{0}+\epsilon^{-1}f_{1}))_{{\mathfrak h}}.
\]
The matrix elements of $\tilde{\lambda}_{\rm vac}$ are:
\[
\tilde\lambda_{++}= \bbbone, \ \ \tilde \lambda_{--}= \tilde \lambda_{+-}= \tilde \lambda_{-+}=0.
\]
It follows that $\lambda_{\rm vac}$ equals to $\lambda(\epsilon^{-1})$ with the notation in Subsect. \ref{ssec:canonic}. Setting
\[
\phi_+(t):=U(t,0)P_+ f=U(t,0)(q\lambda ) f
\]
we have
\[
\phi_+(t)={\textstyle \frac{1}{2}}{\rm e}^{{\rm i} t \epsilon}( f_{0}+  \epsilon^{-1}f_{1}).
\]

\item Let us consider a special case of Theorem \ref{4.3}, namely let $\lambda$ be such that the entries of $\tilde\lambda$ are given by
\[
\tilde{\lambda}_{++}=\bbbone+ b, \ \ \tilde \lambda_{+-}=\tilde \lambda_{-+}^{*}= 0, \ \ \tilde \lambda_{--}= a,
\]
where $a, b\in\Psi^{-\infty}({\mathbb R}^d)$ are both assumed to be positive. The corresponding state is not pure unless $a=b=0$. More explicitly, $\lambda$ is given by
\[
\lambda={\frac{1}{2}}\mat{\epsilon+\epsilon^{{\frac{1}{2}}}(a+b)\epsilon^{{\frac{1}{2}}}}{\bbbone+\epsilon^{{\frac{1}{2}}}(b-a)\epsilon^{-{\frac{1}{2}}}}{\bbbone+\epsilon^{-{\frac{1}{2}}}(b-a)\epsilon^{{\frac{1}{2}}}}{\epsilon^{-1}+\epsilon^{-{\frac{1}{2}}}(a+b)\epsilon^{-{\frac{1}{2}}}}
\]
Defining $\phi_+(t)$ as before we get
\begin{eqnarray*}
\phi_+(t)= & {\textstyle \frac{1}{2}}{\rm e}^{{\rm i} t \epsilon}\left((\bbbone+\epsilon^{-{\frac{1}{2}}}b\epsilon^{{\frac{1}{2}}})f_{0}+ \epsilon^{-1}(\bbbone+\epsilon^{-{\frac{1}{2}}}b\epsilon^{{\frac{1}{2}}})f_{1}\right)
 \\ & + {\textstyle \frac{1}{2}}{\rm e}^{-{\rm i} t \epsilon}\left(\epsilon^{-{\frac{1}{2}}}a\epsilon^{{\frac{1}{2}}} f_{0}- \epsilon^{-{\frac{1}{2}}}a\epsilon^{-{\frac{1}{2}}} f_{1}\right).
\end{eqnarray*}
One can show that the thermal state at inverse temperature $\beta$ is obtained by taking
\[
a=b=\frac{{\rm e}^{- \beta \epsilon}}{\bbbone - {\rm e}^{- \beta \epsilon}}.
\]
\end{itemize}

\section{Hadamard states on general space-times}\label{sec7}\setcounter{equation}{0}

\subsection{Space-times with compact Cauchy surfaces}\label{ssec:junker}
The results  in Sects. \ref{sec1}, \ref{sec2}, \ref{sec3}  and \ref{sec4} extend verbatim to the case where ${\mathbb R}^{d}$ is replaced by a compact manifold $S$. It suffices to replace everywhere  ${\mathcal E}({\mathbb R}^{d})$,  ${\mathcal H}({\mathbb R}^{d})$ and ${\mathcal D}({\mathbb R}^{d})$ by ${\mathcal D}(S)$ and similarly for their dual spaces.  The Weyl pseudo-differential  calculus has to be replaced by the  standard calculus on compact manifolds.
This case is related to the results in \cite{junker,JS02}.

\begin{remark}In \cite{junker,JS02}, a different convention is employed for the symplectic form acting on Cauchy data. This amounts to considering $(\phi(s),\partial_t\phi(s))$ as Cauchy data instead of $(\phi(s),{\rm i}^{-1}\partial_t\phi(s))$. A two-point function $\lambda_{\rm Ju}$ in the convention used in \cite{junker,JS02} corresponds in our notation to the two point function $\lambda=v\lambda_{\rm Ju}v^*$, where $v$ is  diagonal  with entries $v_{++}=\bbbone$ and $v_{--}={\rm i} \bbbone$.
\end{remark}

In \cite[Thm. 5.10]{JS02} it is shown how to construct families of operators $J(t)\in C^{\infty}({\mathbb R},\Psi^1(S))$, $R(t)\in C^{\infty}({\mathbb R},\Psi^0(S))$, such that
\[
\lambda={\frac{1}{2}}\mat{RJ^{-1}R+J}{\bbbone-{\rm i} R J^{-1}}{\bbbone+{\rm i} J^{-1}R}{J^{-1}}
\]
is the two-point function of a pure Hadamard state on ${\mathbb R}\times S$. In our approach, this corresponds to setting
\[
R(t)=\frac{{\rm i}}{2} \left(b(t) - b^*(t)\right), \quad J(t)=\frac{1}{2} \left(b(t) + b^*(t)\right),
\]
where $b(t)$ is as in Theorem \ref{3.1}. Using $r(t)=b^*(t)^{(-1)}$ mod $C^{\infty}({\mathbb R},\Psi^{-\infty})$, it is not difficult to check the microlocal spectrum condition by means of Theorem \ref{4.2}. It is worth pointing out that one of the advantages of basing the construction on $r(t)$ (as we do) rather than on $b(t)$ is that the former is more closely related to the operator $\epsilon(t)$, cf. Theorem \ref{3.1}.

\begin{remark}
 Since $S$ is compact we know that ${\rm Op}(a): H^{m}(S)\to H^{p}(S)$ is Hilbert-Schmidt for any $m,p\in {\mathbb R}$  and $a\in \Psi^{-\infty}({\mathbb R}^{d})$.
 By Shale's theorem it follows that the CCR representations obtained from two pure Hadamard states as in Thm. \ref{thm:canonicform} are unitarily equivalent, since two such states are obtained from one another by a symplectic transformation in $U_{\infty}({\mathcal D}(S)\otimes {\mathbb C}^{2},q)$.
\end{remark}
\subsection{General space-times}\label{ssec:general}

In this subsection we give a new construction of quasi-free Hadamard states on an arbitrary globally hyperbolic space-time $M$ and compare it with the classical construction of Fulling, Narcowich and Wald \cite{FNW}. For the reader's convenience let us first outline this construction.
 \subsubsection{The FNW construction}
Let $(M,g)$ a globally hyperbolic space-time, $S$ a Cauchy hypersurface.  We may assume that
\[
M= {\mathbb R}_{t}\times S, \quad g= -c(x)dt^{2}+ h_{jk}(x)d{\rm x}^{j}d{\rm x}^{k}
\]
and we set $S_{t}:= \{t\}\times S$. We fix a real function $r\in C^\infty(M)$ and consider $P= - \nabla^{a}\nabla_{a}+ r(x)$ --- the associated Klein-Gordon operator. The construction of Hadamard states for $P$ is performed using a `deformation argument'. Namely, one chooses an  ultra-static metric
\[
g^{0}= - dt^{2}+ h_{jk}({\rm x})d{\rm x}^{j}d{\rm x}^{k}, \quad r^{0}(x)= m^{2}>0,
\]
and interpolating  metric $g'= - c'(x)dt^{2}+ h_{jk}'(x)d{\rm x}^{j}d{\rm x}^{k}$, and  real function $r'\in C^\infty(M)$ such that  $(g',r')= (g,r)$ near $[-T/2, T/2]\times S$, $(g',r')= (g^{0}, m^{2})$ near ${\mathbb R}\backslash [-T, T]\times S$.
We denote by $P^{0}, P'$ the associated Klein-Gordon operators, and by $E, E^{0}, E'$ the respective Pauli-Jordan commutator functions.

It is easy to show  that $E'(x, y)= E(x,y)$ (resp. $E'(x,y)=E^{0}(x, y)$) for $x, y$ in a neighborhood of $\{0\}\times S$ (resp. $\{-T\}\times S$).

Let now $\omega_{0}$ be a quasi-free Hadamard state for $P^{0}$. Such Hadamard state always exists on an ultra-static space-time, for instance one can take the associated ground state. Parametrizing elements of ${\rm Sol}_{\rm sc}(P^{0})$ by Cauchy data on $S_{-T}$, we can consider its two point function  $\lambda_{0}$, acting on $C_0^\infty(S_{-T})\otimes {\mathbb C}^{2}$.  Then  $\Lambda'_{0}:= (\rho_{-T}\circ E')^{*}\circ \lambda_{0}\circ (\rho_{T}\circ E')$ satisfies the Hadamard condition for $P'$ in a neighborhood of $S_{-T}$, since the kernels of $E'$ and $E^{0}$ coincide on a neighborhood of $S_{-T}$ and $\omega^{0}$ is Hadamard for $P^{0}$.  By a well-known argument $\Lambda^{'}_{0}$ satisfies the Hadamard condition for $P'$ globally, and defines a Hadamard state $\omega_{0}'$ for $P'$.

We now associate to $\omega_{0}'$ its two-point function $\lambda_{0}'$ acting on $C_0^\infty(S_{0})\otimes{\mathbb C}^{2}$, expressed in terms of Cauchy data on $S_{0}$. Then by the same argument $\Lambda:=(\rho_{0}\circ E)^{*}\circ \lambda_{0}'\circ (\rho_{0}\circ E)$  satisfies the Hadamard condition for $P$ locally near  $S_{0}$  hence globally. In this way we obtain a quasi-free Hadamard state $ \omega$ for $P$.

Denoting  by $U'(-T, 0): {\mathcal D}(S_{0})\otimes{\mathbb C}^{2}\to {\mathcal D}(S_{-T})\otimes {\mathbb C}^{2}$ the propagator for $P'$, mapping the Cauchy data on $S_{0}$ to the Cauchy data on $S_{-T}$ we have:
\begin{equation}\label{ascoli}
\lambda_{0}'= U'(-T, 0)^{*}\circ \lambda_{0}\circ U'(-T, 0).
\end{equation}
Since $U'(-T,0)$ is symplectic we see  by Prop. \ref{istu} that $\omega'_{0}$ is pure iff $\omega_{0}$ is pure, and the same argument gives pureness of $\omega$.

\subsubsection{An alternative construction}
In our approach, we reduce the general problem to the special case of space-times considered by us in Sections \ref{sec2}-\ref{sec3} (or simply to the case of a compact Cauchy surface). Namely, using a suitable partition of unity, we glue together Hadamard states on smaller regions of the space-time. \medskip

The steps of the construction are the following:

\medskip

We fix a Cauchy surface $S$, so that we can assume that $M={\mathbb R}\times S$ and the metric $g$ is of the form (\ref{metric}).

We choose  an open set $\Omega$ in $M$ and  for $n\in {\mathbb N}$,  open, pre-compact sets  $U_{n}$, $\tilde{U}_{n}$ in $S$,  constants
 $0<\delta_{n}$    such that:
\begin{equation}
\label{Had.e1}
\begin{array}{rl}
(i)&U_{n}\Subset \tilde{U}_{n}, \ \bigcup_{n}U_{n}=S,\\[2mm]
(ii)&\tilde{U}_{n}\hbox{ are coordinate charts for }S,\\[2mm]
(iii)&y\in \Omega, \ J(y)\cap U_{n}\neq \emptyset \ \Rightarrow \ y\in ]-\delta_{n}, \delta_{n}[\times \tilde{U}_{n}=: \tilde{\Omega}_{n},\\[2mm]
(iv)&\Omega\hbox{ is a neighborhood of }S\hbox{ in }M.
\end{array}
\end{equation}
In {\it (iii)} $J(y)$ denotes the causal shadow of $y\in M$.
This is clearly possible.

We fix a partition of unity $1= \sum_{n}\chi_{n}^{2}$ of $S$, with  $\chi_{n}\in C_0^\infty(U_{n})$ for $n\in {\mathbb N}$. Denoting still by $\chi_{n}$ the map $\chi_{n}\otimes \bbbone$ on $C_0^\infty(S)\otimes {\mathbb C}^{2}$, we note that
\begin{equation}
\label{Had.e2}
q= \sum_{n}\chi_{n}^{*}q\chi_{n}.
\end{equation}

Fix for each $n\in {\mathbb N}$ a coordinate map $\varphi_{n}: \tilde{U}_{n}:\to \tilde{V}_{n}$,
where $V_{n}$ is a neighborhood of $0$ in ${\mathbb R}^{d}$.  The symplectic form $\sigma$
on $C_0^\infty(\tilde{U}_{n})\otimes {\mathbb C}^{2}$ transported to $C_0^\infty(\tilde{V}_{n})\otimes{\mathbb C}^{2}$ will be given by  (\ref{e2.0}).

We also transport  with $\varphi_{n}$ the Klein-Gordon operator $P$ on $\tilde{\Omega}_{n}$ to an operator on $]-\delta_{n}, \delta_{n}[\times V_{n}\subset {\mathbb R}\times {\mathbb R}^{d}$.  We can extend this operator to ${\mathbb R}\times {\mathbb R}^{d}$ so that it satisfies the conditions in Sect. \ref{sec2}. Let us denote by $P_{n}$ the Klein-Gordon operator on ${\mathbb R}\times {\mathbb R}^{d}$ obtained in this way. We choose for each $n\in {\mathbb N}$ a two-point function $c_{n}$ (acting on the  space of Cauchy data) of a quasi-free state, which is Hadamard for $P_{n}$.  We will have in particular
\begin{equation}\label{Had.e3}
c_{n}\geq 0, \ c_{n}\geq q.
\end{equation}
We restrict $c_{n}$ to $C_0^\infty(\tilde{V}_{n})\otimes {\mathbb C}^{2}$ and transport it back to $C_0^\infty(\tilde{U}_{n})\otimes {\mathbb C}^{2}$ by $\varphi_{n}^{-1}$, denoting it by $\lambda_{n}$. Finally we set:
\begin{equation}\label{ascola}
\lambda:= \sum_{n\in {\mathbb N}}\chi_{n}^{*}\lambda_{n}\chi_{n},
\end{equation}
which is well defined as a two-point function on $C_0^\infty(S)\otimes {\mathbb C}^{2}$, since the sum is locally finite.

Since $q= \sum \chi_{n}^{*}q\chi_{n}$, and $c_{n}$ was the two-point function of a gauge-invariant quasi-free state, we deduce from (\ref{Had.e3}) that  $\lambda\geq 0$, $\lambda\geq q$, i.e. $\lambda$ is the two-point function of a gauge-invariant quasi-free state.

It remains to check that $\lambda$ satisfies the Hadamard condition, i.e. that
\[
\Lambda= (\rho\circ E)^{*}\lambda (\rho\circ E)
\]
 satisfies $(\mu{\rm sc})$. By the well-known propagation  property of \cite{FSW}(see also \cite{SV}), it  suffices to check $(\mu{\rm sc})$ in $T^{*}\Omega\times T^{*}\Omega\backslash\{0\}$,  since $\Omega$
is a neighborhood of $S$ in $M$ by (\ref{Had.e1}).  Set:
\[
\Lambda_{n}:= (\rho\circ E)^{*}\chi_{n}^{*}\lambda_{n}\chi_{n} (\rho\circ E),
\]
so that  $\Lambda= \sum_{n}\Lambda_{n}$. It suffices to check that $\Lambda_{n}$
satisfies $(\mu{\rm sc})$ in  $T^{*}\Omega\times T^{*}\Omega\backslash\{0\}$.

Using the support properties of $E$ and condition (\ref{Had.e1}) {\it (iii)}, we obtain that the restriction to $\Omega\times \Omega$ of the distribution kernel of $\Lambda_{n}$ is supported in $\tilde{\Omega}_{n}\times \tilde{\Omega}_{n}$. Therefore, up to diffeomorphisms, $\Lambda_n$ is equal to
\[
C_{n}:= (\rho\circ E_n)^{*}\chi_{n}^{*}\lambda_{n}\chi_{n} (\rho\circ E_n)
\]
on $\tilde{\Omega}_{n}\times \tilde{\Omega}_{n}$, where $E_n$ is the propagator associated to $P_n$. Using the invariance of the wavefront set under diffeomorphisms, it follows that $\Lambda_n$ satisfies the microlocal spectrum condition.  Therefore $\Lambda$ is the two-point function of a gauge-invariant quasi-free Hadamard state.

\subsubsection{Comparison between the two constructions}
 The FNW construction has the important advantage to generate pure Hadamard states for an arbitrary space-time from pure Hadamard states for an ultra-static one.  On the other hand, the relation (\ref{ascoli}) between the two two-point functions expressed in terms of Cauchy data is difficult to control in practice, since it involves the propagator $U'(-T, 0)$ for the intermediate metric $g'$, which is a non-local operator.

In contrast, our construction does  in general not generate a pure state, even if all the local states $\omega_{n}$ are so.  However the two-point function  given by (\ref{ascola}) is easier to control in practice, since it depends only on the two-point functions $\lambda_{n}$ and on the local operators $\chi_{n}$.