\magnification1200
\tolerance=10000
\hsize=17truecm\vsize=23truecm
\multiply\baselineskip by 15\divide \baselineskip by10
\parindent=40pt
\font\eightpoint=cmsl8



%\input mssymb
\def\Bbb{\bf} %if \Bbb not available

% shorthand

\def\<{\langle}
\def\>{\rangle}

\def\header#1#2{\hfil\underbar{#1} \hfil\bigskip
                \headline={{\eightpoint #2\hfill\today}}}

% misc math stuff

\def\slim{\mathop{\hbox{\rm s-lim}}}
\def\Ker{\mathop{\rm Ker}}
\def\Ran{\mathop{\rm Ran}}
\def\Re{\mathop{\rm Re}}
\def\Im{\mathop{\rm Im}}
\def\ess{\hbox{\it ess}}
\def\d{\hbox{\it d}}
\def\c{\hbox{\it c}}
\def\pp{\hbox{\it pp}}

\def\const{{\rm const}}
\def\parderivs#1#2{{\partial#1\over \partial#2}}
\def\parderiv#1{{\partial\over \partial#1}}
\def\secparderivs#1#2{{\partial^2#1\over \partial#2^2}}
\def\secparderiv#1{{\partial^2\over \partial#1^2}}
\def\derivs#1#2{{d #1\over d #2}}
\def\deriv#1{{d\over d #1}}
\def\secderivs#1#2{{d^2#1\over d #2^2}}
\def\secderiv#1{{d^2\over d #1^2}}
\def\nderivs#1#2#3{{d^{#1}#3\over d #2^{#1}}}
\def\nderiv#1#2{{d^{#1}\over d #2^{#1}}}
\def\inprod#1{\langle#1\rangle}
\def\max{{\rm max}}
\def\min{{\rm min}}


%symbols

\def\L{{\rm L}}

\def\R{{\bf R}}
\def\Clow{\lower1pt\hbox{\bf C}}
\def\C{{\bf C}}

\def\O{{\cal O}}
\def\H{{\cal H}}
\def\W{{\cal W}}
\def\J{{\cal J}}
\def\F{{\cal F}}
\def\P{{\cal P}}
\def\V{{\cal V}}
\def\U{{\cal U}}
\def\N{{\cal N}}
\def\Q{{\cal Q}}
\def\B{{\cal B}}
\def\D{{\cal D}}

% misc

%
%       the stuff after \ifdraftmode only appears if \draftmode has previously
%       appeared
%
\def\draftmode{\let\dfmd 1}
\def\ifdraftmode#1{\ifx\dfmd 1{#1}\fi}




%
%\def\today{\ifcase\month\or January\or February \or March \or April\or May\or
%           June\or July\or August\or September\or October \or November
%           \or December\fi
%           \space\number\day, \number\year}
%



\def\today{\number\month\slash \number\day\slash \number\year}


%\def\today{7\slash 9\slash 1994}



\def\quote#1{{\narrower\smallskip\noindent\it#1\smallskip}}

\def\Box{\vbox{\hrule
              \hbox{\vrule height7pt\kern7pt\vrule height7pt}
              \hrule
              }
        }

\def\Picbox{\centerline{
                \vbox{\hrule
                \hbox{\vrule height2.5in\kern5in\vrule height2.5in}
                \hrule
                }
            }
        }



%
%       the following macros read in the files \jobname.eqr
%       and \jobname .thm which
%       contain the names of equations and how they are nested
%       within sections in a format like
%       \\{name1}\bgnscn\\{name2}\bgnscn\\{name3}\ndscn\ndscn
%       In this example \eqlabelname1 will expand to 1
%                       \eqlabelname2 will expand to 1.1
%                       \eqlabelname3 will expand to 1.1.1
%       This is the maximum nesting supported
%       In the theorem list each \\{name} is preceded by \typetheorem,
%       \typelemma or \typeproposition to indicate what it is
%



\def\bgnsctn{
        \advance\secnest by1
        \cntnumber=0
        \ifnum\secnest=1
                \advance\seclevelone by1
                \edef\prenum{\number\seclevelone .}
        \fi
        \ifnum\secnest=2
                \advance\secleveltwo by1
                \edef\prenum{\number\seclevelone .\number\secleveltwo .}
        \fi
        }
\def\ndsctn{\advance\secnest by-1}

\def\thmtype{} %initialize
\def\typetheorem{\def\thmtype{Theorem}}
\def\typelemma{\def\thmtype{Lemma}}
\def\typeproposition{\def\thmtype{Proposition}}
\def\typecorollary{\def\thmtype{Corollary}}
\def\typehypothesis{\def\thmtype{Hypothesis}}




%
%       \beginsection begins a new section
%       writes \bgnsctn to token list in \eqrtoks
%       \endsection ends a section
%       writes \ndsctn to token list in \eqrtoks
%
\def\beginsection{
        \advance\secnumber by1
        \probctr=0
        \toksdef\ta=0
        \ta={\bgnsctn}
        \immediate\write\eqroutfile{\the\ta}
        \immediate\write\thmoutfile{\the\ta}
        }
\def\endsection{
        \toksdef\ta=0
        \ta={\ndsctn}
        \immediate\write\eqroutfile{\the\ta}
        \immediate\write\thmoutfile{\the\ta}
        }

%
%       \num{name} inserts (<number>) where <number> is the equation number
%       referred to by "name"
%
\def\num#1{(\csname eqlabel#1\endcsname )}

%
%       \thmnum is the analogous thing for theorems
%
\def\thmnum#1{\hbox{\csname typethmlabel#1\endcsname}}



%
% first do the references - define \\ to be a macro which defines \eqlabelname
% to be a reference number and then evaluate the input string as in Appendix D
%
\def\\#1{
        \advance\cntnumber by1
        \expandafter\xdef\csname
        eqlabel#1\endcsname{\prenum\number\cntnumber}
        }
\newcount\seclevelone \seclevelone=0
\newcount\secleveltwo \secleveltwo=0
\newcount\cntnumber \cntnumber=0
\newcount\secnest \secnest=0
\def\prenum{}
\newread\eqrinfile
\openin\eqrinfile=\jobname .eqr
\ifeof\eqrinfile
  \def\eqrinput{\relax}
\else\loop
  \read\eqrinfile to\eqrinput
  \eqrinput
  \ifeof\eqrinfile\let\endfile 1\else\let\endfile 0\fi
  \ifx\endfile 0 \repeat
\fi
\closein\eqrinfile

%
% now the same thing for theorems
%
\def\\#1{
        \advance\cntnumber by1
        \expandafter\xdef\csname
        thmlabel#1\endcsname{\prenum\number\cntnumber}
        \expandafter\xdef\csname
        typethmlabel#1\endcsname{\thmtype\ \prenum\number\cntnumber}
        }
\seclevelone=0
\secleveltwo=0
\cntnumber=0
\secnest=0
\def\prenum{}
\newread\thminfile
\openin\thminfile=\jobname .thm
\ifeof\thminfile
  \def\thminput{\relax}
\else\loop
  \read\thminfile to\thminput
  \thminput
  \ifeof\thminfile\let\endfile 1\else\let\endfile 0\fi
  \ifx\endfile 0 \repeat
\fi
\closein\thminfile

%
% theorems etc
% These are numbered by section
%

\def\begintheoremlabel#1{
        \toksdef\ta=0
        \ta={\typetheorem\\{#1}}
        \immediate\write\thmoutfile{\the\ta}
        \medbreak\par\noindent{\ifdraftmode{\mcomment{#1}}\bf Theorem
\csname thmlabel#1\endcsname :}
        \begingroup \it
}
\def\endtheoremlabel{\par\medbreak\endgroup}

\def\begintheorem{
        \advance\nonamectr by1
        \edef\temp{{noname\the\nonamectr}}
        \expandafter\begintheoremlabel\temp
}
\def\endtheorem{\endtheoremlabel}

\def\beginlemmalabel#1{
        \toksdef\ta=0
        \ta={\typelemma\\{#1}}
        \immediate\write\thmoutfile{\the\ta}
        \medbreak\par\noindent{\ifdraftmode{\mcomment{#1}}\bf Lemma \csname
thmlabel#1\endcsname :}
        \begingroup \it
}
\def\endlemmalabel{\par\bigbreak\endgroup}

\def\beginlemma{
        \advance\nonamectr by1
        \edef\temp{{noname\the\nonamectr}}
        \expandafter\beginlemmalabel\temp
}
\def\endlemma{\endlemmalabel}



\def\beginpropositionlabel#1{
        \toksdef\ta=0
        \ta={\typeproposition\\{#1}}
        \immediate\write\thmoutfile{\the\ta}
        \medbreak\par\noindent{\ifdraftmode{\mcomment{#1}}\bf Proposition
\csname thmlabel#1\endcsname :}
        \begingroup \it
}
\def\endpropositionlabel{\par\medbreak\endgroup}

\def\beginproposition{
        \advance\nonamectr by1
        \edef\temp{{noname\the\nonamectr}}
        \expandafter\beginpropositionlabel\temp
}
\def\endproposition{\endpropositionlabel}

\def\begincorollarylabel#1{
        \toksdef\ta=0
        \ta={\typecorollary\\{#1}}
        \immediate\write\thmoutfile{\the\ta}
        \medbreak\par\noindent{\ifdraftmode{\mcomment{#1}}\bf Corollary
\csname thmlabel#1\endcsname :}
        \begingroup \it
}
\def\endcorollarylabel{\par\medbreak\endgroup}

\def\begincorollary{
        \advance\nonamectr by1
        \edef\temp{{noname\the\nonamectr}}
        \expandafter\begincorollarylabel\temp
}
\def\endcorollary{\endcorollarylabel}

\def\beginhypothesislabel#1{
        \toksdef\ta=0
        \ta={\typehypothesis\\{#1}}
        \immediate\write\thmoutfile{\the\ta}
        \medbreak\par\noindent{\ifdraftmode{\mcomment{#1}}\bf Hypothesis
\csname thmlabel#1\endcsname :}
        \begingroup \it
}
\def\endhypothesislabel{\par\medbreak\endgroup}

\def\beginhypothesis{
        \advance\nonamectr by1
        \edef\temp{{noname\the\nonamectr}}
        \expandafter\beginhypothesislabel\temp
}
\def\endhypothesis{\endhypothesislabel}

%
% The following allows one to give a theorem or lemma name
% from a different part of the paper
% to a result quoted here.
%

\def\begintheoremref#1{
        \bigbreak\par\noindent{\bf \hbox{\csname typethmlabel#1\endcsname}:}
        \begingroup \it
                       }
\def\endtheoremref{\par\medbreak\endgroup}



%
%
%       open files for writing
%
%
\newwrite\eqroutfile
\immediate\openout\eqroutfile=\jobname .eqr
\newwrite\thmoutfile
\immediate\openout\thmoutfile=\jobname .thm




%
%       \be{name} and \ee begin and end a numbered displayed equation
%       the tokens in the name are added to the token list in \eqrtoks
%       if \eqlabelname has been defined uses this as a number
%       \lastdisplayline{name}{...} gives a referenced number  to the last
%       line in a displaylines series
%
\def\be#1{
        \def\temp{#1}
        \toksdef\ta=0
        \ta={\\{#1}}
        \immediate\write\eqroutfile{\the\ta}$$
        }
\def\ee{
        \eqno{(\csname eqlabel\temp\endcsname)
        \ifdraftmode{\rlap{\ {\smallfont\temp}}}}$$
        }
\def\eel#1{
        \eqno{(\csname eqlabel\temp\endcsname #1)
        \ifdraftmode{\rlap{\ {\smallfont\temp}}}}$$
        }
\def\lastdisplayline#1#2{
        \hfill#2\hfill\llap{(\csname eqlabel#1\endcsname)
        \ifdraftmode{\rlap{\ {\smallfont#1}}}}\cr
        \toksdef\ta=0
        \ta={\\{#1}}
        \immediate\write\eqroutfile{\the\ta}
        }



\def\beginproof{\bigbreak\par\noindent{\bf proof:\ }}
%\def\endproof{\thinspace\Box\par\medbreak}
\def\endproof{\hfill\vrule height .75em width .75em depth 0pt\bigbreak}

\def\beginproofof#1{\bigbreak\par\noindent{\bf proof of #1:\ }}
%\def\endproofof{\thinspace\Box\par\medbreak}



\def\begindefn{\medbreak\par\noindent{\it Definition:\/\ }}
\def\enddefn{\par\medbreak}

\def\beginexample{\medbreak\par\noindent{\it Example:\/\ }}
\def\endexample{\par\medbreak}




\def\chead#1{\medbreak\centerline{#1}\medskip}

\def\bea#1{$$\eqalign\begingroup}
\def\eea{\endgroup $$}
\def\beann{$$\eqalign\begingroup}
\def\eeann{\endgroup $$}


%
%       These go at the beginning  of input
%
\global\newtoks\eqrtoks
\global\newtoks\thmtoks
\global\newcount\secnumber
\global\newcount\nonamectr
\global\newcount\probctr

\message{reference macros need to run twice to get references right}

\def\la{\langle}
\def\ra{\rangle}

\hfill\today

\ \vskip1cm

\centerline{\bf Stimulated particle emission -- a simple example.}

\vskip1cm
\centerline{Christopher King$^{(1)}$ and Roger Waxler$^{(2)}$}
\bigskip
\centerline{\it $^{(1)}$ Department of Mathematics,}\par
\centerline{\it \hfill Northeastern University, Boston, Mass. 02115 \hfill}\par
\centerline{\it $^{(2)}$ Department of Mathematics, S.U.N.Y. at Buffalo,
Buffalo, N.Y., U.S.A. 14214}


\bigskip\bigskip
{\narrower\bigskip\noindent{\bf Abstract:}\/
We find evidence of stimulated photon emission
in a simple model of the atom plus radiation
system, namely a two state atom interacting with a quantised
scalar field in the rotating wave approximation.
We use an initial state where the atom is excited and
there is one photon present, and the photon energy is
close to the energy difference between the two atomic states.
By taking the long-time, weak-coupling limit, we compute in closed
form the probability that the atom remains in its excited state
as a function of the rescaled time. Our results indicate that
stimulated emission does occur for certain initial photon states,
but only at quadratic order in the rescaled time.
We also find partial recurrence of the excited state.
\bigskip}


\vfill\eject


\noindent{\bf 1. \underbar{Introduction.}}
\medskip

The broad subject of this paper is the interaction of matter
and radiation, in particular the absorption and emission of
photons from a one-electron atom. The literature on this subject
is vast -- we refer the reader to [CDG], [D], [St], [WW] for the basic
physics.
Recently there has been considerable progress in the mathematical
analysis of the Hamiltonian which describes this interaction
([BFS1], [BFS2]). In particular it has been shown how
an excited state of the atom turns into a resonance for the
atom-radiation system. The evolution of initial states in a
quantum mechanical model is
strongly influenced by the existence of such resonances, as they
generally determine the long-time asymptotics of solutions of
Schr\"odinger's equation. In the atom-radiation model
they lead to decay of excited states, along with
emission of photons. For additional background on the mathematics
of quantum radiation, resonances and
decay of states, we refer the reader to [A1], [A2], [Da], [DE],
[Fro], [HSp], [Ho], [Hu], [JP], [K1], [OY],
[Sp2].

The influence of
such a resonance on the time evolution operator can be studied
in a simplified model where the atom is replaced by a two-state system,
the radiation field is replaced by a massless scalar field
and the rotating wave approximation is used in the interaction.
This model and its generalizations have been widely studied --
some examples are to be found in [Da], [EY], [Fri], [GT], [JC], [K1], [K2],
[Pf], [Sp1], [W].
In particular the issue of spontaneous emission of a photon has been
examined; here one follows the time evolution of the initial state
${\Psi}_{0}$ where the atom is excited
and there are no photons present.
In this paper we will consider initial states
${\Psi}_{1}$ for this same model where the atom is excited
and there is {\it one} photon present. Our goal is to find
evidence for ``stimulated emission", namely an increase in the
rate at which the atom relaxes to its ground state,
due to the presence of the photon.

The precise formulation of such a result requires some choices. We follow
the lead of Davies [Da] and explore a long-time, weak-coupling
limit which can be used to define the rate of spontaneous emission. Namely we
introduce a parameter $\beta$ which governs the size of the coupling
between atom and field, and we rescale the time coordinate $t$
by defining $\tau = {\beta}^{2}t$. We then
take the limit $\beta \rightarrow 0$, $t \rightarrow \infty$ with
$\tau$ fixed. As Davies showed ([Da]), when the initial state is ${\Psi}_{0}$
the probability that the atom remains in its excited state can be
computed exactly in this limit. The result is
$\exp [- 2 \Gamma \tau]$ where $\Gamma$ is the usual decay rate.
As mentioned above, in this paper we choose an initial state where the
atom is excited, and there is one photon present. In order to obtain
an interesting result (that is, something other than $\exp [- 2 \Gamma
\tau]$),
we must choose the photon state to depend
on the parameter $\beta$, in such a way that its energy is localised to
within $O({\rm {\beta}^{2}})$ of the atom's energy gap.
Let $p(\beta,t)$ be the resulting probability that the atom remains excited at
time $t$, with this initial state. It turns out that we can explicitly
compute $\alpha(\tau)=\lim_{\beta \rightarrow 0} p(\beta, {\beta}^{-2}\tau)$
as a function of $\tau$ (see equations $\num{result2}$ --
$\num{result4}$ in section 2.5). It shows interesting
behavior, depending on the precise form of the initial one-photon
state. Initially $\alpha(\tau)$ decreases more rapidly than the corresponding
probability with no initial photon, and goes to zero
exponentially as $\tau \rightarrow \infty$. However it behaves as
$1 - 2 \Gamma \tau$ for small $\tau$, which means that the {\it linear}
decay rate is unaffected by the initial one-photon state. In addition,
$\alpha(\tau)$ may actually increase at some later times. Presumably this
records some re-excitation of the atom by the initial one-photon state.
Such a ``rebound" or recurrence effect is well-known in models
with discrete spectrum, for example the Jaynes-Cummings model ([JC]),
but occurs less often in models with only continuous spectrum,
as is the case here.


We prove the result by using a resummed perturbation theory,
based on a resolvent expansion. Specifically we separate the interaction
into its direct and exchange parts, and we resum the direct
part. The technical part of our work
is concerned with obtaining bounds
for the remainder terms in this expansion, and one of the chief
tools here is dilation analyticity (for background on this topic,
see [AC], [BC], [Hu], [Si]).
It is worth noting that our results
cannot be obtained using the standard formulation of time-dependent
perturbation theory for this model. Indeed repeating the calculation
which leads to Fermi's Golden Rule for this model shows that the
rate of decay is independent of the initial one-photon state.
The reason is that time-dependent perturbation theory
yields a useful expansion for times which are less than
$O({\beta}^{-1})$, while our results show that
the effects of stimulated emission
occur at times which are $O({\beta}^{-2})$.


The paper is organised as follows. In section 2 we define our model and
state the main results in $\thmnum{first theorem}$ and $\thmnum{second
theorem}$.
In section 3 we use spectral theory to
derive our resummed perturbation theory, and prove $\thmnum{first theorem}$.
Section 4 contains the
results about stimulated emission, and the proof of $\thmnum{second theorem}$.
In section 5 we prove convergence
of our resummed perturbation theory, and the technical estimates
needed in sections 3 and 4.

Concerning notation, we are cavalier about numerical factors in the
estimates throughout the paper. In particular we use $``C"$ to denote a
positive number which is independent of $t$, $\beta$ and any variables.
Its value may change from line to line.

\vfill\eject

\noindent{\bf 2. \underbar{Definitions and Results.}}


\medskip
\noindent{\it 2.1 The model.}
\smallskip

First we define the model. The state space is ${\bf C}^{2} \otimes {\cal F}$,
where $\cal F$ is the symmetric Fock space over $L^{2}({\bf R})$
(see [RS], sections II.4 and VIII.10 for relevant definitions).
For any $g \in L^{2}({\bf R})$, let $a(g)$ and $a^{\dagger}(g)$ denote
respectively the annihilation and creation operators acting on $\cal F$.
They are defined by their action on ${\cal F}_{n}$, the $n$-particle
subspace of $\cal F$. Let ${\phi}_{n}(k_{1},\dots,k_{n})$ be a state
in ${\cal F}_{n}$; then the state $a^{\dagger}(g) {\phi}_{n}$ lies in
${\cal F}_{n+1}$, and is given by
\be{cr.op}
\big(a^{\dagger}(g) \phi_n\big)(k_1, \dots, k_{n+1})=
{1 \over \sqrt{n+1}} \sum_{j=1}^{n+1}\phi_n(k_1, \dots, \hat k_j,
\dots , k_{n+1})g(k_j)
\ee
where the argument $\hat k_j$ on the right side is omitted.
Similarly, for $n \geq 1$, $a(g) \phi_n$ lies in $\F_{n-1}$, and is given by
\be{ann.op}
\big(a(g) \phi_n\big)(k_1, \dots, k_{n-1})=
\sqrt{n} \int \phi_n(k_1, \dots, k_n) g(k_n) dk_n
\ee
The subspace $\F_0 \simeq {\bf C}$ is spanned by the Fock vacuum $\Omega$,
and $a(g)\Omega=0$.

Let $h$ denote the second quantised Hamiltonian for the free scalar field
on $\cal F$ defined by the single particle dispersion relation
$\omega(k) = |k|$. Since $h$ preserves each subspace ${\cal F}_{n}$,
it is defined by its restriction to ${\cal F}_{n}$, which we call
$h_{n}$; for ${\phi}_{n}(k_{1},\dots,k_{n})$ in ${\cal F}_{n}$,
\be{en.op}
h_{n} {\phi}_{n}(k_{1},\dots,k_{n}) = \sum_{j=1}^{n}
|k_{j}| {\phi}_{n}(k_{1},\dots,k_{n})
\ee

The full Hamiltonian $H(\beta)$ for our model
is most conveniently written as a $2 \times 2$ matrix,
with respect to the decomposition ${\bf C}^{2} \otimes {\cal F} =
{\cal F} \oplus {\cal F}$:
\be{full.Ham}
H(\beta) = \left(\matrix{
h+E & \beta a({\bar f}) \cr & \cr \beta a^{\dagger}(f) & h \cr}\right)
\ee
Here $E > 0$ is the energy gap between the ground state and excited
state of the atom, $\beta$ is the (real) coupling parameter,
and $f(k)$ is the atomic form factor. We shall make strong and explicit
regularity assumptions for $f$ below.
Note that for all $n \geq 0$, the subspace
\be{inv.sub}
{\cal H}_{n} = {\cal F}_{n} \oplus {\cal F}_{n+1}
\ee
is invariant under $H(\beta)$. Hence under the decomposition
${\bf C}^{2} \otimes {\cal F} = \oplus_{n=0}^{\infty} {\cal H}_{n}$,
$H(\beta)$ has block diagonal form. We let $H_{n}(\beta)$ denote its
restriction to ${\cal H}_{n}$.

We shall be concerned with
the dynamics generated by $H_{1}(\beta)$ on ${\cal H}_{1}$.
Explicitly,
\be{H1.def}
H_{1}(\beta) = \left(\matrix{
h_{1}+E & \beta V \cr & \cr \beta V^{\dagger} & h_{2}\cr}\right)
\ee
where the operators $h_{1},h_{2},V,V^{\dagger}$ act as follows:
let ${\phi}_{1} \in {\cal F}_{1} = L^{2}(R)$ and
${\phi}_{2} \in {\cal F}_{2} =
\big(L^{2}(R) \otimes L^{2}(R)\big)_{\rm sym}$;
\be{H1.detail}\eqalign{
(h_{1}{\phi}_{1})(k) &= |k| {\phi}_{1}(k)\cr
(h_{2}{\phi}_{2})(k_{1},k_{2}) &= (|k_{1}| + |k_{2}|)
{\phi}_{2})(k_{1},k_{2})\cr
(V^{\dagger}{\phi}_{1})(k_{1},k_{2}) &= {1 \over \sqrt{2}}
[{\phi}_{1}(k_{1})f(k_{2}) + {\phi}_{1}(k_{2})f(k_{1})]\cr
(V{\phi}_{2})(k_{1}) &= \sqrt{2} \int {\bar f}(k_{2})
{\phi}_{2}(k_{1},k_{2}) d k_{2}\cr}
\ee

\medskip
\noindent{\it 2.2 Dilation analyticity.}
\smallskip

We say that a function $g(k)$ is ${\cal F}_{1}$-dilation analytic
in the set ${\cal S} \subset  {\bf C}$, if $g_{\theta}$ is analytic
as an ${\cal F}_{1}$-valued function of $\theta$ for all $\theta \in {\cal S}$,
where
\be{dil.an}
g_{\theta}(k) = e^{-i \theta/2} g(e^{-i\theta}k)
\ee
Henceforth the statement ``$g(k)$ is dilation analytic" will mean
that $g$ is ${\cal F}_{1}$-dilation analytic in a set
${\cal S} \subset  {\bf C}$ which contains some strip
$\{ \theta \,|\, -a \leq
{\rm Re}\theta \leq b \}$, for some $a,b > 0$.
Note in this case it is also true that the function ${\bar g}(k)$
is dilation analytic, and $({\bar g})_{\theta}(k)
= \overline{g_{- \bar \theta}(k)}$.
Furthermore the product of dilation analytic functions is
dilation analytic, as is the inverse of a dilation analytic
function which is bounded away from zero.

Suppose $f(k)$ and $g(k)$ are both
dilation analytic in some strip containing
${\rm Re}\theta = 0$. Then the inner product
$\la f_{- {\bar \theta}}, g_{\theta} \ra$ is independent of
$\theta$. This follows because it is analytic in
$\theta$, and when $\theta$ is pure imaginary the
$\theta$-dependence can be removed by a
change of variables. Next, recall the multiplication operator
$h_{1}$ defined in $\num{H1.detail}$.
Let $h_{1,\theta}$ be the following multiplication operator:
\be{dil.h1}
(h_{1,\theta}{\phi}_{1})(k) = e^{-i\theta} |k| {\phi}_{1}(k).
\ee
Suppose that ${\rm Im}z>0$. Then the inner product
$\la f_{- {\bar \theta}}, (h_{1,\theta}-z)^{-1} g_{\theta} \ra$
is independent of $\theta$ in the region
$0 \leq {\rm Re}\theta \leq a$ for some $a > 0$.
Again this follows by noting
that the inner product is analytic as a function of $\theta$
in some open set containing the strip
$0 \leq {\rm Re}\theta \leq a$,
and that the $\theta$-dependence can be removed by a change of
variables when ${\rm Re}\theta=0$. Similarly if ${\rm Im}z<0$,
$\la f_{- {\bar \theta}}, (h_{1,\theta}-z)^{-1} g_{\theta} \ra$
is independent of $\theta$ in the region
$-a < {\rm Re}\theta \leq 0$.

Finally we say that a bounded operator $A$ on ${\cal F}_{1}$ is
dilation analytic if $A\psi$ is dilation analytic
for every $\psi$ which is dilation analytic, and if there
exists $A_{\theta}$ bounded on ${\cal F}_{1}$ such that
\be{dil.op}
(A \psi)_{\theta} = A_{\theta} {\psi}_{\theta}.
\ee


\medskip
\noindent{\it 2.3 Assumptions on $f(k)$.}
\smallskip

For comparison we recall some results for the spin-boson Hamiltonian.
This Hamiltonian has the same form as $\num{full.Ham}$, with
the off-diagonal terms both replaced by
$\beta \big(a({\bar f})+a^{\dagger}(f)\big)$. For this model,
it is known that a sufficient condition for essential self-adjointness
and existence of a unique ground state is
$\int |f(k)|^{2} (1 + |k|^{-2})d k < \infty$ ([Sp1], [HSp]).
We now state our assumptions on the form factor.
\medskip
\par\noindent i) $f(k)$ is real, and $f(k) = f(-k)$ for all $k$.
\par\noindent ii) $f(E) \neq 0$.
\par\noindent iii) $f(k)$ is ${\cal F}_{1}$-dilation analytic
in some strip $-{\theta}_{0} \leq {\rm Re} \theta \leq {\theta}_{0}$,
with ${\theta}_{0} > 0$.
\par\noindent iv) For some $\delta >0$, $C>0$,
$|f_{\theta}(k)|^{2} \leq C |k|^{1+\delta}$
for all $|k| \leq 1$, and for all
$-{\theta}_{0} \leq {\rm Re} \theta \leq {\theta}_{0}$.
\par\noindent v) $|f(k)|^{2} \leq C \,(1 + k^{2})^{-1}$
for all $k$.

\medskip
\noindent Some comments on these assumptions are in order:
(i) is unessential, and is chosen for convenience;
(ii) is essential -- otherwise the transition is forbidden
in leading order;
(iii) is essential for our method to work;
(iv) is necessary for our estimates --
it is slightly stronger than
the conditions stated above for the spin-boson model;
(v) gives a useful shortcut to our result, and is
physically reasonable.


\medskip
\noindent{\it 2.4 Special initial states.}
\smallskip

In order to find evidence of stimulated emission, we consider
a family of initial states which depend on $\beta$. Specifically, we choose
states which are localised at energy $E$ to within
O(${\beta}^{2}$). For ${\rm Im}w \neq 0$
we define
\be{special}
{\psi}^{(w)}(k) = {\kappa(w) \over |k| -E + w}, \quad\quad\quad
{\Psi}^{(w)} = \left(\matrix{{\psi}^{(w)} \cr 0\cr}\right)
\ee
where $\kappa(w)$ is chosen so that $||{\psi}^{(w)}||_{2} = 1$.
Then our special initial states will be ${\Psi}^{({\beta}^{2}w)}$,
and we will take ${\rm Im}w > 0$ for definiteness.
Note that
\be{asymp.norm}
\lim_{\beta \rightarrow 0} {\beta}^{-2} |\kappa({\beta}^{2}w)|^{2}
= {|{\rm Im}w| \over 2 \pi}
\ee


\medskip
\noindent{\it 2.5 Results.}
\smallskip

We define
\be{def.u}
u = \lim_{\eta \rightarrow 0}\
\la f, (h_{1} - E - i \eta)^{-1} f \ra \,\,
{\buildrel \rm def \over =}\,\, E_{1} + i \Gamma
\ee
It follows from section 2.3 that
\be{}
\Gamma = 2 \pi |f(E)|^{2} >0.
\ee
Note that $2 {\beta}^{2} \Gamma$ is the rate of decay of the excited state
which follows by applying the Fermi Golden Rule to this model.

We write $P_{1}$ and $P_{2}$ for the orthogonal projections onto
the $1$- and $2$-particle subspaces of ${\cal H}_{1}$ respectively:
\be{def.P}
P_{1} = \left(\matrix{1 & 0\cr 0 & 0\cr}\right),\quad\quad
P_{2} = \left(\matrix{0 & 0\cr 0 & 1\cr}\right).
\ee
Let $\psi \in {\cal F}_{1}$, and
$\Psi = \left(\matrix{{\psi} \cr 0\cr}\right)$. We are interested
in the time evolution of the initial state $\Psi$, projected onto
the $1$-particle subspace. Accordingly let
\be{def.p1}
p_{1}(\beta, t; \psi) =
||P_{1} \big(e^{-iH_{1}(\beta)t} \Psi \big)||^{2}
\ee
We examine the bahavior of this probability as a function of the
rescaled time $\tau = {\beta}^{2}t$, in the limit
$\beta \rightarrow 0$.

\medskip
\begintheoremlabel{first theorem} Let $\psi(k)$ be dilation analytic.
Then for $\tau \geq 0$,
\be{result1}
\lim_{\beta \rightarrow 0} \,p_{1}(\beta, {\beta}^{- 2}\tau; \psi)
= e^{-2 \Gamma \tau}
\ee

\endtheoremlabel

\medskip
\noindent In order to state the next result we need some notation. Define
\be{def.K}
K(u,w;\tau,\alpha;v) =
{e^{i \tau u} \over v+w}\bigg((\alpha + 1) e^{-i \tau v}
- e^{i \tau w}\bigg) -
{e^{i \tau w} \over v+u}\bigg(e^{-i \tau v}
- e^{i \tau u}\bigg)
\ee
Recall the definition of ${\psi}^{(w)}$ in $\num{special}$, and of
$u$ in $\num{def.u}$.


\medskip
\begintheoremlabel{second theorem}
For ${\rm Im}w > 0$,
\be{result2}
\lim_{\beta \rightarrow 0} \,p_{1}(\beta,
{\beta}^{-2} \tau; {\psi}^{({\beta}^{2}w)}) =
{4 {\Gamma}^{2} {\rm Im}w \over \pi |u-w|^{2}}
\int_{-\infty}^{\infty}\bigg|
K(u,w;\tau,{u-w \over 2i\Gamma};v)\bigg|^{2} d v.
\ee

\endtheoremlabel


\medskip

The integral in $\num{result2}$ can be computed using Cauchy's theorem; we
record
the result in section 4.6. One case of interest is when
${\rm Re}(w) = E_{1} = {\rm Re}(u)$, so we state this as a separate
result below. We rescale time by setting $s = \Gamma \tau$, and we write
$w=E_{1}+ic\Gamma$.

\medskip
\begincorollarylabel{first cor}
For any $c > 0$,
\be{result3}\eqalign{
\lim_{\beta \rightarrow 0} \,p_{1}(\beta,
{\beta}^{-2} {\Gamma}^{-1}s;
{\psi}^{\big({\beta}^{2}(E_{1}+ic\Gamma)\big)}) = &
\big({3-c \over 1-c}\big)^{2} e^{-2s}
+ {4c \over (1-c)^{2}} e^{-2sc} \cr
& - {8c(3-c) \over (1-c)^{2}(1+c)}
e^{-s-sc} - {8 \over 1+c} e^{-2s-2sc}\cr}
\ee
\endcorollarylabel

\medskip

As one would expect, the right side
of $\num{result3}$ approaches the right side of $\num{result1}$
as $c \rightarrow \infty$.
In fact this also happens in the limit $c \rightarrow 0$.
So the interesting behavior occurs for $0 < c < \infty$.
Perhaps the
most interesting case is when $w=u$. We present this as another special case
below.

\medskip
\begincorollarylabel{second cor}
\be{result4}
\lim_{\beta \rightarrow 0} \,p_{1}(\beta,
{\beta}^{-2}{\Gamma}^{-1}s; {\psi}^{({\beta}^{2}u)}) =
[5 - 8 s + 4 s^{2}] e^{-2s} - 4 e^{-4s}
\ee

\endcorollarylabel

\medskip
\noindent{\it 2.6 Discussion.}
\smallskip

We present $\thmnum{second theorem}$ and its corollaries
as evidence for stimulated emission in this model.
The initial one photon state ${\psi}^{({\beta}^{2}w)}$
is designed to mimic the spontaneously
emitted photon state, so we expect that it should produce an enhanced decay
rate for the atom. In position space it is centered at the atom at
time zero. As the results show, enhancement does occur
at some values of the rescaled time ${\beta}^{-2}t$.
Notice that for $\tau$ small, the right sides of $\num{result1}$ and
$\num{result2}$ are both
approximately $1-2\Gamma \tau$, so they have the same linear decay rate.
Therefore the effects of stimulated emission show up at quadratic order
in $\tau = {\beta}^{2}t$, which gives another way of seeing that
application of standard
time dependent perturbation theory does not yield an enhanced decay
rate.

As noted in the introduction, the probability that the atom remains excited may
increase at later times. This is a little unexpected in our
model, and we take this interesting feature as motivation
for the continued study of this most
`elementary' quantum field theory model.


\vfill\eject


\noindent{\bf 3. \underbar{Stone's formula and the resolvent expansion.}}

\medskip
\noindent{\it 3.1 The setup.}
\smallskip

To prove our results we use a resolvent expansion
which is essentially a resummation of the usual perturbative
expansion. We separate the interaction between atom and field into
a direct part and an exchange part, and we resum the direct part.

For any $\Psi \in {\cal H}_{1}$, Stone's formula [RS] gives
\be{stone}
e^{-iH_{1}(\beta)t} \Psi = \lim_{\eta \downarrow 0} {1 \over 2 \pi i}
\int_{-\infty}^{\infty} e^{-itx} \bigg[\big(H_{1}(\beta) - x - i \eta\big)^{-1}
\Psi - \big(H_{1}(\beta) - x + i \eta\big)^{-1} \Psi \bigg] d x
\ee
The resolvent $\big(H_{1}(\beta) - z\big)^{-1}$ can be expressed as
a $2 \times 2$ matrix with respect to the direct sum decomposition
${\cal H}_{1} = {\cal F}_{1} \oplus {\cal F}_{2}$
(see $\num{inv.sub}$); for ${\rm Im}z \neq 0$,
\be{eval.resol}
\big(H_{1}(\beta) - z\big)^{-1} =
\left(\matrix{G(\beta,z) & -\beta G(\beta,z)V(h_{2}-z)^{-1}\cr
& \cr
-\beta (h_{1}-z)^{-1}V^{\dag}G(\beta,z) & M(\beta,z)\cr}\right)
\ee
where
\be{def.G}
G(\beta,z) = \big(h_{1}+E-z-{\beta}^{2}V(h_{2}-z)^{-1}V^{\dag}\big)^{-1}
\ee
\be{def.M}
M(\beta,z) = \big(h_{2}-z-{\beta}^{2}V^{\dag}(h_{1}+E-z)^{-1}V\big)^{-1}
\ee
For our purposes it will be enough to apply Stone's formula to  a state
$\Psi = \left(\matrix{\psi \cr 0\cr}\right)$, and to project the
result onto the first component. Therefore we will study
\be{proj.stone}
P_{1}(e^{-iH_{1}(\beta)t}\Psi) = \lim_{\eta \downarrow 0} {1 \over 2 \pi i}
\int_{-\infty}^{\infty} e^{-itx} \bigg[G(\beta,x + i \eta)
\psi - G(\beta,x - i \eta) \psi \bigg] d x
\ee

Our resummed expansion is obtained by splitting the term
$V(h_{2}-z)^{-1}V^{\dag}$ in $\num{def.G}$ into two factors, corresponding to
the direct and exchange parts of the interaction. For ${\rm Im}z \neq0$,
we define the following operators on ${\cal F}_{1}$:
\be{def.S}
(S(z) \psi)(k) = \bigg(\int_{-\infty}^{\infty}
{|f(p)|^{2} \over |p| + |k| - z} d p \bigg) \psi(k)
\ee
\be{def.T}
(T(z) \psi)(k) = \bigg(\int_{-\infty}^{\infty}
{{\bar f(p)} \psi(p) \over |p| + |k| - z} d p \bigg) f(k)
\ee
Then for all $\psi \in {\cal F}_{1}$,
\be{S+T}
V(h_{2}-z)^{-1}V^{\dag}\psi = S(z)\psi + T(z)\psi.
\ee

We also define
\be{def.G0}
G_{0}(\beta,z) = \big(h_{1}+E-z-{\beta}^{2}S(z)\big)^{-1}
\ee
\be{def.G1}
G_{1}(\beta,z) = {\beta}^{2}G_{0}(\beta,z)T(z)G_{0}(\beta,z)
\ee
\be{def.GR}
G_{R}(\beta,z) = {\beta}^{4}G_{0}(\beta,z)T(z)G(\beta,z)T(z)G_{0}(\beta,z)
\ee
Using $\num{S+T}$, and applying the second resolvent formula to
$\num{def.G}$, we obtain
\be{expand.G}
G(\beta,z) = G_{0}(\beta,z) + G_{1}(\beta,z) + G_{R}(\beta,z)
\ee

\medskip
\noindent{\it 3.2 Leading order parts and remainder.}
\smallskip

On substituting the right side of $\num{expand.G}$ into Stone's formula,
we obtain an expansion for the time evolved state. Define
\be{def.T*}
T_{\sigma}(\psi,t) = \lim_{\eta \rightarrow 0} {1 \over 2 \pi i}
\int_{-\infty}^{\infty} e^{-itx} \big[
G_{\sigma}(\beta,x+i\eta) \psi - G_{\sigma}(\beta,x-i\eta) \psi \big] d x
\ee
where $\sigma \in \{0,1,R\}$. Then we obtain
\be{expand.in.T}
P_{1}\bigg(e^{-iH_{1}(\beta)t} \left(\matrix{{\psi} \cr 0\cr}\right)
\bigg) = T_{0}(\psi,t) + T_{1}(\psi,t) + T_{R}(\psi,t)
\ee
The following lemmas give bounds for the three terms on the right
side of $\num{expand.in.T}$. The first result, $\thmnum{bound T0,1}$,
will be proved in section 4.5. Recall the definition of $u$ in
$\num{def.u}$.

\medskip
\beginlemmalabel{bound T0,1}
Let $|\beta|$ be sufficiently small and
let $\psi$ be dilation analytic. There is $z_{0} \in {\bf C}$,
with $|z_{0}-E+{\beta}^{2}u|={\rm O}({\beta}^{4})$,
such that for all $t \geq 0$,
\be{boundforT0,1}
||T_{1}(\psi,t)|| \leq C\, {\beta}^{2}, \quad\quad
||T_{0}(\psi,t) - e^{-it(h_{1} + z_{0})} \psi||
\leq C\, {\beta}^{2}
\ee
\endlemmalabel

The number $z_{0}$ appearing in $\num{boundforT0,1}$ will
be defined below in $\num{props.g}$.
The next lemma gives a bound for the remaining term on the right side
of $\num{expand.in.T}$; the proof is in section 5. Recall that
$\delta$ was defined in section 2.3. It will be convenient to
use the following additional notation:
\be{def.cases}
{\beta}^{\{2,\delta\}} = \cases{
{\beta}^{2 \delta} & if $0 < \delta < 1$\cr
\ln ({\beta}^{-2}) {\beta}^{2}& if $\delta=1$\cr
{\beta}^{2}& if $\delta>1$\cr}
\ee

\medskip
\beginlemmalabel{bound TR}
For $|{\beta}|$ sufficiently small, for all $t \geq 0$ and all
$\psi \in {\cal F}_{1}$,
\be{boundforTR}
||T_{R}(\psi, t)|| \leq C\, {\beta}^{\{2,\delta\}} \, ||\psi||
\ee

\endlemmalabel


\medskip
\noindent{\it 3.3 Proof of $\thmnum{first theorem}$.}
\smallskip

Combining the results of $\thmnum{bound T0,1}$ and $\thmnum{bound TR}$
we see that for
$|\beta|$ small enough, for any dilation
analytic state $\psi$ with $||\psi||=1$ there is $C$ such that
\be{est1}
||P_{1}\bigg(e^{-iH_{1}(\beta)t} \left(\matrix{{\psi} \cr 0\cr}\right)
\bigg) - e^{-it(h_{1} + z_{0})} \psi||
\leq C\,{\beta}^{\{2,\delta\}}
\ee
for all $t \geq 0$. Since
\be{est2}
||e^{-it(h_{1} + z_{0})} \psi|| = e^{-t ({\beta}^{2}\Gamma +
{\rm O}({\beta}^{4}))},
\ee
the result follows immediately.
\endproof

\medskip
\noindent{\it 3.4 Propagators.}
\smallskip

We require some detailed properties of the propagators
$G_{0}(\beta,z)$ and $G_{1}(\beta,z)$. To this end, we define
the following functions for ${\rm Im}z \neq 0$:
\be{def.Pi}
\Pi(z) = \la f, (h_{1} - z)^{-1} f \ra = \int_{-\infty}^{\infty}
{|f(p)|^{2} \over |p| - z} d p
\ee
\be{def.g}
g(z) = E - z - {\beta}^{2} \Pi(z).
\ee
It follows that $S(z)\psi(k) = \Pi(z-|k|)\psi(k)$, so we have
\be{G.relate.g}
\big(G_{0}(\beta,z)\psi\big)(k) = g(z - |k|)^{-1} \psi(k).
\ee
Following the comments in Section 2.2, for ${\rm Im}z > 0$
and $0 \leq {\rm Re}\theta \leq {\theta}_{0}$, and for
${\rm Im}z < 0$ and $0 \geq {\rm Re}\theta \geq -{\theta}_{0}$, we have
\be{dilate.Pi}
\Pi(z) = \la f_{- {\bar \theta}}, (h_{1,\theta}-z)^{-1} f_{\theta}\ra
= \int_{-\infty}^{\infty} \overline{f_{- {\bar \theta}}(p)}
f_{\theta}(p) (e^{-i\theta}|p| - z)^{-1} d p
\ee
Define the two subsets of the complex plane:
\be{def.calD}\eqalign{
{\cal D}_{+} &= \{z \,:\,  |z| > 0, -{\theta}_{0} < {\rm Arg}z < \pi \} \cr
{\cal D}_{-} &= \{z \,:\,  |z| > 0, -\pi < {\rm Arg}z < {\theta}_{0} \} \cr}
\ee
Then letting $\theta = \pm {\theta}_{0}$
(see section 2.3) $\num{dilate.Pi}$
provides the analytic continuations of $\Pi(z)$ from
the region ${\rm Im}z > 0$
to ${\cal D}_{+}$, and from the region
${\rm Im}z < 0$ to
${\cal D}_{-}$. We denote the former by ${\Pi}_{+}(z)$ and
the latter ${\Pi}_{-}(z)$. These in turn define the
corresponding analytic continuations of $g(z)$ to
${\cal D}_{+},{\cal D}_{-}$ which we denote by
$g_{+}(z)$ and $g_{-}(z)$.

We list some properties of ${\Pi}_{+}(z)$:
\par\noindent (i) $
{\rm Im}(z) \geq 0 \quad \Rightarrow \quad
{\rm Im}\big({\Pi}_{+}(z)\big) \geq 0
$
\par\noindent (ii) $\lim_{\eta \downarrow 0}{\Pi}_{+}(i \eta)
= \int |p|^{-1} |f(p)|^{2} d p < \infty$.
\par\noindent (iii) for all $z \in {\cal D}_{\pm}$,
$$
|{\Pi}_{\pm}(z)| \leq \int_{-\infty}^{\infty}
|f_{\mp  {\theta}_{0}}(p)|\,|f_{\pm  {\theta}_{0}}(p)|
\, |p|^{-1}\,d p \,< \, \infty
$$
\par\noindent (iv) ${\Pi}_{+}(E) = u = E_{1} + i \Gamma$.

Since by assumption ${\rm Im}({\Pi}_{+}(E)) = \Gamma > 0$,
it follows that there is $\epsilon > 0$ such that
${\rm Im}\big({\Pi}_{+}(z)\big) > 0$ for all $z$ satisfying
$|z - E| \leq \epsilon$.

For $|\beta|$ sufficiently small, $g_{+}(z)$ has a unique zero in
${\cal D}_{+}$, which we denote $z_{0}$. It satisfies the following:
\be{props.g}
g_{+}(z_{0})=0, \quad
g_{+}'(z_{0}) = -1 + O({\beta}^{2}), \quad
z_{0} = E - {\beta}^{2}u + O({\beta}^{4})
\ee
The following bound is easily verified; for $|\beta|$ sufficiently
small and for all $z \in {\cal D}_{+}$,
\be{bound1.g}
|g_{+}(z)| \geq {|z - z_{0}| \over 2}
\ee
Together with $f(E) \neq 0$ this in turn implies that there
is $C > 0$ such that for all $z$ with ${\rm Im}z \geq 0$,
\be{bound2.g}
|g_{+}(z)| \geq C(|z - E|^{2} + {\beta}^{4})^{1/2}
\ee
Similar bounds hold for $g_{-}(z)$:
for all $z \in {\cal D}_{-}$,
\be{bound3.g}
|g_{-}(z)| \geq {|z - {\bar z}_{0}| \over 2};
\ee
for all $z$ with ${\rm Im}z \leq 0$,
\be{bound4.g}
|g_{-}(z)| \geq C(|z - E|^{2} + {\beta}^{4})^{1/2}.
\ee

Returning to $\num{G.relate.g}$ we deduce that for any function
$\psi(k) \in {\cal F}_{1}$,
and any $z$ with ${\rm Im}z \neq 0$,
\be{bound1.G0}
|(G_{0}(\beta,z)\psi)(k)| \leq C
\big(|z - |k| - E|^{2} + {\beta}^{4}\big)^{-1/2} |\psi(k)|.
\ee

We shall also need some estimates for the dilated version of the operator
$G_{0}$ (see section 2.2). Since the operator $G_{0}$
is multiplication by a function of $|k|$,
its dilation is obtained by replacing $|k|$ everywhere by
$e^{-i\theta}|k|$. So we define
\be{def.G.theta}
(G_{0,\theta}(\beta,z)\psi)(k) =
g(z - e^{-i\theta}|k|)^{-1} \psi(k).
\ee
For ${\rm Im}z > 0$, $g(z - e^{-i\theta}|k|)$ is analytic
in $\theta$ for $\theta$ in some open set containing the
strip $0 \leq {\rm Re}\theta < {\theta}_{0}$.
By considering separately the regions
$|e^{-i\theta}|k| + E -z| \leq \epsilon$
and $|e^{-i\theta}|k| + E -z| > \epsilon$
for $\epsilon$ sufficently small, we get the
bounds
\be{bound5.g}\eqalign{
|g(z - e^{-i\theta}|k|)| = &
\big|e^{-i\theta}|k| + E -z - {\beta}^{2}
\Pi(z-e^{-i\theta}|k|)\big| \cr
\geq & C \, \max
({\beta}^{2}\Gamma, {\rm Re}(\theta)|k| + {\rm Im}z)}
\ee
This shows that
the operator $G_{0,\theta}(\beta,z)$ is analytic in $\theta$
and $z$ in the region $0 \leq {\rm Re}\theta \leq {\theta}_{0}$
and ${\rm Im}z >0$, and satisfies the bound
\be{bound.G.theta}
|(G_{0,\theta}(\beta,z)\psi)(k)| \leq C \,
\min \{{1 \over {\beta}^{2}}, {1 \over |k||{\rm Re}\theta|},
{1 \over |{\rm Im}z|} \}
\ee
By similar reasoning the same bound holds for ${\rm Im}z<0$
and $-{\theta}_{0} \leq {\rm Re}\theta \leq 0$.

\vfill\eject


\noindent{\bf 4. \underbar{Computation of the leading order terms.}}

\medskip
\noindent{\it 4.1 The setup.}
\smallskip


As shown by $\thmnum{bound TR}$ in section 3.2, the leading order contribution
to the time evolved state is provided by the two terms $T_{0}(\psi,t)$
and $T_{1}(\psi,t)$. In this section we will extract and compute the
most important parts of these terms when the state $\psi$ is a special
initial state as defined in $\num{special}$, and estimate the error terms.


\medskip
\noindent{\it 4.2 Leading order part of $T_{0}$.}
\smallskip

The leading order part of $T_{0}(\psi,t)$ can be extracted for
any state $\psi$.
For any $\psi \in {\cal F}_{1}$, we have
\be{eval.T0}
T_{0}(\psi,t)(k) = {1 \over 2\pi i} \lim_{\eta \rightarrow 0}
\int_{-\infty}^{\infty} e^{-itx} \big[g(x+i \eta -|k|)^{-1} \psi(k)
- g(x-i \eta -|k|)^{-1} \psi(k) \big] d x
\ee
Recall the definition of $z_{0}$ in $\num{props.g}$. Define:
\be{def.L0}
L_{0}(\psi, \beta, t;k) = e^{-it(|k|+z_{0})} {\psi}(k).
\ee

\medskip
\beginlemmalabel{T0-L0}
For all $t \geq 0$, and $\psi$ with $||\psi|| = 1$,
\be{bound.T0-L0}
\bigg(\int_{-\infty}^{\infty} \Big|T_{0}({\psi},t)(k)
 - L_{0}(\psi, \beta, t;k)\Big|^{2} \bigg)^{1/2}
\leq C {\beta}^{2}
\ee
\endlemmalabel

\beginproof
First, it is sufficient to prove that for all $k$,
\be{est3}
\Big|{1 \over 2\pi i}
\int_{-\infty}^{\infty} e^{-itx} \big[g_{+}(x -|k|)^{-1}
- g_{-}(x -|k|)^{-1} \big] d x -
e^{-it(|k|+z_{0})}\Big| \leq C {\beta}^{2}
\ee
Second, the bounds $\num{bound2.g}$, $\num{bound4.g}$ imply that for
$|x-|k|-E| \geq E/2$, we have
\be{est4}
|g_{\pm}(x-|k|)|^{-1} \leq C |x - |k| - E|^{-1}.
\ee
Combining this with the formula
\be{resol.expand}
g_{+}(z)^{-1} - g_{-}({z})^{-1} =
g_{+}(z)^{-1} g_{-}({z})^{-1}[-
{\beta}^{2} {\Pi}_{-}({z})
+ {\beta}^{2} {\Pi}_{+}(z)],
\ee
and recalling that ${\Pi}_{\pm}(z)$ is uniformly bounded,
it is sufficient to show that
\be{est5}
\Big|{1 \over 2\pi i}
\int_{|k|+E/2}^{|k|+3E/2} e^{-itx} \big[g_{+}(x -|k|)^{-1}
- g_{-}(x -|k|)^{-1} \big] d x -
e^{-it(|k|+z_{0})}\Big| \leq C {\beta}^{2}
\ee

Let $S$ be the semicircular arc in the lower half plane
which begins at $|k|+E/2$ and ends at $|k|+3E/2$.
For $|\beta|$ sufficiently small the pole of
$g_{+}(z -|k|)^{-1}$ lies above $S$. Therefore by
deforming the contours of integration we get
\be{int.S}\eqalign{
{1 \over 2\pi i}
\int_{|k|+E/2}^{|k|+3E/2}  e^{-itx} & \big[g_{+}(x -|k|)^{-1}
- g_{-}(x -|k|)^{-1} \big] d x \cr
 =
{1 \over 2\pi i}
\int_{S} e^{-itz} & \big[g_{+}(z -|k|)^{-1}
- g_{-}(z -|k|)^{-1} \big] d z \cr
& - {\rm Res} \big[e^{-itz} g_{+}(z -|k|)^{-1}; z=z_{0}+|k|\big]\cr}
\ee
where Res($F(z);z=w$) denotes the residue of the meromorphic
function $F(z)$ at $z=w$. Using
$\num{resol.expand}$ again, and the bounds $\num{bound1.g}$,
$\num{bound3.g}$,
we deduce that the integral along $S$ is ${\rm O}({\beta}^{2})$.
Finally the evaluation of the residue and the estimate
$g_{+}'(z_{0})=-1+ {\rm O}({\beta}^{2})$
yield the stated result.

\endproof

\medskip
\noindent{\it 4.3 Leading order part of $T_{1}$.}
\smallskip

Unlike the estimate in section 4.2 for $T_{0}$, we must now use the
particular form of the initial state.
Recall that the special initial state is
$$
\psi(k) = {\psi}^{({\beta}^{2}w)}(k) =
{\kappa({\beta}^{2}w) \over |k|-E+{\beta}^{2}w}
$$
where Im$(w)>0$, and $\kappa$ is chosen so that
$||{\psi}^{({\beta}^{2}w)}|| = 1$.
Define for ${\rm Im}z \neq 0$,
\be{def.B}
B(z,w,|k|) = \int_{-\infty}^{\infty}{{\bar f}(p) \over
|p| + |k| - z} g(z-|p|)^{-1} {\psi}^{({\beta}^{2}w)}(p)
d p
\ee
Then we can write
\be{eval.T1.B}\eqalign{
T_{1}({\psi}^{({\beta}^{2}w)},t)(k)=
{{\beta}^{2}f(k) \over 2 \pi i} \lim_{\eta \rightarrow 0}
\int_{-\infty}^{\infty} &e^{-itx} \big[g(x+i \eta -|k|)^{-1}
B(x+i \eta,w,|k|)\cr
& - g(x-i \eta -|k|)^{-1} B(x-i \eta,w,|k|)
\big] d x\cr}
\ee
Using $f(p)=f(-p)$ real we can rewrite $\num{def.B}$ as
\be{def2.B}
B(z,w,|k|) = 2 \int_{0}^{\infty}{{f}(p) \over
p + |k| - z} g(z-p)^{-1}
{\kappa({\beta}^{2}w) \over p-E+{\beta}^{2}w}
d p
\ee
For ${\rm Im}(z)>0$, the integrand in
$\num{def2.B}$ has only one pole below the real axis,
namely $p=E-{\beta}^{2}w$. Its residue is
\be{def.b}
b(z,w,|k|) = - 2 \pi i \kappa({\beta}^{2}w)
{\rm Res}\, \bigg({2 {f}(p) \over p+|k|-z} g(z-p)^{-1};
p = E - {\beta}^{2} w \bigg)
\ee
where ${f}$ now denotes the analytic continuation
off the real axis. To estimate the contribution of $\num{def.b}$ to the
integral $\num{def2.B}$, we will rotate the contour of integration.
Accordingly we define
\be{def.B.theta}
B_{\theta}(z,w,|k|)=
2 \int_{0}^{\infty}{{f}(pe^{-i\theta}) \over
pe^{-i\theta} + |k| - z} g(z-pe^{-i\theta})^{-1}
{\kappa({\beta}^{2}w) \over pe^{-i\theta}-E+{\beta}^{2}w}
e^{-i\theta} d p
\ee

\medskip
\beginlemmalabel{bound B}
There exists $C$ such that for all
$\eta>0$, all $\theta$ with ${\theta}_{0}/2 \leq
{\rm Re}\theta \leq {\theta}_{0}$,
\be{bound.B.theta}
\int_{-\infty}^{\infty}\big|
B_{\pm\theta}(x\pm i \eta,w,|k|) \big|^{2} d x
\leq C {\beta}^{2}
\ee
\endlemmalabel

\beginproof
Let $\theta \geq {{\theta}_{0}/2}$ and $\eta>0$.
We use the bounds (recall $\num{bound1.g}$):
\be{bounds}\eqalign{
|g(x+i \eta-pe^{-i\theta})| & \geq
{1 \over 2}|x+i \eta - p e^{-i\theta} -z_{0}|
\geq C [p^{2} {\sin }^{2}\theta
+ (x- p\cos \theta -{\rm Re}z_{0})^{2}]^{1/2} \cr
|pe^{-i\theta}+|k|-x-i \eta| & \geq C [p^{2} {\sin }^{2}\theta
+ (x- p\cos \theta -|k|)^{2}]^{1/2} \cr
|pe^{-i\theta}-E+{\beta}^{2}w| & \geq
[p^{2} {\sin }^{2}\theta + (E- p\cos \theta)^{2}]^{1/2}
\geq C \cr}
\ee
Hence an application of Cauchy-Schwarz shows that
\be{est7}\eqalign{
\big| B_{\theta}(x+i \eta,w,|k|) \big|^{2} \leq & C
|\kappa({\beta}^{2}w)|^{2} \bigg(\int_{0}^{\infty}
|{f}(pe^{-i\theta})|[p^{2} {\sin }^{2}\theta
+ (x-p\cos \theta-{\rm Re}z_{0})^{2}]^{-1} d p \bigg)\cr
& \bigg(\int_{0}^{\infty}
|{f}(qe^{-i\theta})|[q^{2} {\sin }^{2}\theta
+ (x-q\cos \theta-|k|)^{2}]^{-1} d q \bigg)\cr}
\ee
We now integrate over $x$. We interchange the $x$ integral
with the integrals over $p$ and $q$, and use Cauchy-Schwarz
to estimate the integral over $x$. The result is
\be{est8}
\int_{-\infty}^{\infty}
\big| B_{\theta}(x+i \eta,w,|k|) \big|^{2} d x
\leq C |\kappa({\beta}^{2}w)|^{2} \bigg(\int_{0}^{\infty}
|{f}(pe^{-i\theta})| p^{-3/2} d p\bigg)^{2}
\ee
By assumption the integral is finite (see section 2.3, (iv)), and
$|\kappa({\beta}^{2}w)|^{2} = {\rm O}({\beta}^{2})$,
and hence the result follows.
The analysis for $B_{-\theta}(x-i \eta,w,|k|)$ is
identical.
\endproof

\medskip

We now observe that for ${\rm Im}z<0$,
$B_{\theta}(z,w,|k|)$ is independent of
$\theta \leq 0$, since in this case the integrand on the right side of
$\num{def2.B}$ has no poles
above the real axis. For ${\rm Im}z>0$, and
$\theta \geq {\theta}_{0}/2$, an application of
Cauchy's theorem gives
\be{Cauchy}
B(z,w,|k|) = b(z,w,|k|) + B_{\theta}(z,w,|k|)
\ee
Hence for ${\theta}_{0}/2 \leq \theta \leq {\theta}_{0}$,
\be{rotate.T1}\eqalign{
T_{1}({\psi}^{({\beta}^{2}w)},t)(k) =
{{\beta}^{2}f(k) \over 2 \pi i} \lim_{\eta \downarrow 0}
\int_{-\infty}^{\infty} & e^{-itx} g(x+i \eta -|k|)^{-1}
b(x+i \eta,w,|k|) d x \cr
+ {{\beta}^{2}f(k) \over 2 \pi i} \lim_{\eta \downarrow 0}
\int_{-\infty}^{\infty} & e^{-itx} \big[g(x+i \eta -|k|)^{-1}
B_{\theta}(x+i \eta,w,|k|) \cr
& - g(x-i \eta -|k|)^{-1}
B_{-\theta}(x-i \eta,w,|k|) \big] d x \cr}
\ee
We now use $\num{bound2.g}$, $\num{bound4.g}$ to deduce
\be{est9}
\int_{-\infty}^{\infty} \big|
g(x \pm i \eta - |k|) \big|^{-2} d x\leq C {\beta}^{-2}.
\ee
It follows by Cauchy-Schwarz and $\num{bound.B.theta}$ that
\be{est10}
\bigg| \int_{-\infty}^{\infty} e^{-itx}
g(x \pm i \eta - |k|)^{-1}
B_{\pm \theta}(x \pm i \eta,w,|k|) d x \bigg| \leq C,
\ee
and hence the second term on the right side of $\num{rotate.T1}$ is bounded by
$C {\beta}^{2} |f(k)|$. Therefore the dominant contribution
to $T_{1}({\psi}^{({\beta}^{2}w)},t)$ comes from the first term
on the right side of $\num{rotate.T1}$. This integral has three poles,
all located at distance
${\rm O}({\beta}^{2})$ below the real axis (see $\num{eval.b}$ below
for the evaluation of this residue): these are
\be{poles}\eqalign{
x_{1} = & |k| + z_{0}\cr
x_{2} = & E + |k| - {\beta}^{2}w\cr
x_{3} = & E + z_{0} - {\beta}^{2}w\cr}
\ee
Given $a,b,c$ with nonzero imaginary parts, we define
\be{def.J}
J(a,b,c;t) = \int_{-\infty}^{\infty} e^{-ity}
(y-a)^{-1} (y-b)^{-1} (y-c)^{-1} d y
\ee
The next lemma shows that the leading order
part of $\num{rotate.T1}$ is the following expression
(recall $\num{def.u}$ $u = {\Pi}_{+}(E)$):
\be{def.L1}
L_{1}(w,\beta,t;k) = 2 {\beta}^{2} \kappa({\beta}^{2}w)
{f}(E) f(k)
 J(|k|+E-{\beta}^{2}u,
|k|+E-{\beta}^{2}w, 2E-{\beta}^{2}u
-{\beta}^{2}w;t)
\ee

\medskip
\beginlemmalabel{T1-L1}
\be{bound.T1-L1}
\bigg(\int_{-\infty}^{\infty} \Big|
T_{1}({\psi}^{({\beta}^{2}w)},t)(k) -
L_{1}(w,\beta,t;k) \Big|^{2} d k \bigg)^{1/2}
\leq C |\beta|.
\ee

\endlemmalabel

\beginproof
It is enough to establish the bound
$\Big|T_{1}({\psi}^{({\beta}^{2}w)},t)(k) -
L_{1}(w,\beta,t;k) \Big| \leq C |\beta| |f(k)|$ for each $k$.
Also using the estimate $\num{est10}$,
it is sufficient to establish such a
bound for the difference between the first term on the right side
of $\num{rotate.T1}$ and the function $L_{1}(w,\beta,t;k)$. The evaluation
of the residue $b$ gives
\be{eval.b}
b(x+i \eta,w,|k|) = - 4 \pi i \kappa ({\beta}^{2}w)
{f}(E-{\beta}^{2}w) (E+|k|-x-i \eta -{\beta}^{2}w)^{-1}
g(x+i \eta -E+{\beta}^{2}w)^{-1}
\ee
where $f$ again denotes the analytic continuation.
Inserting this in $\num{rotate.T1}$ we see that the integrand has three
poles at the points $x_{1},x_{2},x_{3}$ defined in $\num{poles}$.
In the resulting expression
we want to replace $g_{+}(x-|k|)^{-1}$ by
$(x-x_{1})^{-1}$ and $g_{+}(x-E+{\beta}^{2}w)^{-1}$ by
$(x-x_{3})^{-1}$, and replace ${f}(E-{\beta}^{2}w)$
by ${f}(E)$. This can be accomplished by noting that
(i) only the neighborhoods of the poles at
$x_{1},x_{2},x_{3}$ contribute to the integral -- the integral over the
remaining intervals is bounded by $C {\beta}^{2} |f(k)|$;
(ii) in a neighborhood of ${\rm Re}x_{1}$, we have $g_{+}(x-|k|)^{-1}
= (x-x_{1})^{-1} (-1 + {\rm O}({\beta}^{2}))$ for real $x$
(see $\num{props.g}$), and in a
neighborhood of ${\rm Re}x_{3}$ we have
$g_{+}(x-E+{\beta}^{2}w)^{-1}=
(x-x_{3})^{-1}(-1 + {\rm O}({\beta}^{2}))$ for real $x$;
(iii) $f(E-{\beta}^{2}w)=f(E)+{\rm O}({\beta}^{2})$.
Now the following crude estimate
\be{crude.est}\eqalign{
\int_{-\infty}^{\infty}& |x-x_{1}|^{-1}
|x-x_{2}|^{-1} |x-x_{3}|^{-1} d x \cr
& \leq
\bigg(\int |x-x_{1}|^{-2}  d x \bigg)^{1/2}
\bigg(\int |x-x_{2}|^{-2}  |x-x_{3}|^{-2} d x \bigg)^{1/2}\cr
& \leq \bigg(\int |x-x_{1}|^{-2}  d x \bigg)^{1/2}
\bigg(\int |x-x_{2}|^{-4}  d x \bigg)^{1/4}
\bigg(\int |x-x_{3}|^{-4}  d x \bigg)^{1/4}\cr
& \leq C {\beta}^{-4} \cr}
\ee
shows that
\be{est11}
\bigg|{{\beta}^{2}f(k) \over 2 \pi i} \int_{-\infty}^{\infty}
e^{-itx} g_{+}(x-|k|)^{-1} b(x,w,|k|) d x \bigg| \leq
C {\beta}^{-1} |f(k)|.
\ee
Therefore if we replace $g_{+}(x-|k|)^{-1}$
by $(x-x_{1})^{-1}$ as discussed above, the error is bounded
by $C|\beta||f(k)|$. For the same reason we can
replace $g_{+}(x-E+{\beta}^{2}w)^{-1}$ by
$(x-x_{3})^{-1}$, and
$f(E-{\beta}^{2}w)$ by $f(E)$, obtaining the same error.
Putting it together, we have shown that
\be{est12}
\Big|T_{1}({\psi}^{({\beta}^{2}w)},t)(k)
- 2 {\beta}^{2} \kappa({\beta}^{2}w)
{f}(E) f(k) J(x_{1},x_{2},x_{3};t)\Big|
\leq C\,|\beta|\,|f(k)|
\ee
Now we use the estimate $\num{crude.est}$ again to deduce that
for any ${\epsilon}_{1},{\epsilon}_{2},{\epsilon}_{3}$,
\be{est13}
\Big|J(x_{1},x_{2},x_{3};t) -
J(x_{1}+{\beta}^{4}{\epsilon}_{1},
x_{2}+{\beta}^{4}{\epsilon}_{2},
x_{3}+{\beta}^{4}{\epsilon}_{3};t)\Big| =
{\rm O}({\beta}^{-2})
\ee
Since $z_{0}=E-{\beta}^{2}u+{\rm O}({\beta}^{4})$, $\num{est13}$
implies that we can replace $J(x_{1},x_{2},x_{3};t)$ in
$\num{est12}$ with $J(|k|+E-{\beta}^{2}u,|k|+E-{\beta}^{2}w,
2E-{\beta}^{2}u-{\beta}^{2}w;t)$, with error
$C|\beta||f(k)|$, thus proving the result.

\endproof



\medskip
\noindent{\it 4.4 Proof of $\thmnum{second theorem}$.}
\smallskip

Recall $\num{def.L0}$. With some abuse of notation,
we write
\be{newdef.L0}
L_{0}(w,\beta,t;k) = L_{0}({\psi}^{({\beta}^{2}w)},\beta,t;k)
\ee
Then by the results of the previous sections, we have
\be{est14}
||P_{1}(e^{-iH_{1}(\beta)t} {\Psi}^{({\beta}^{2}w)})||^{2}
= \int_{-\infty}^{\infty}|L_{0}(w,\beta,t;k) +
L_{1}(w,\beta,t;k)|^{2} d k  + {\rm O}(|\beta|)
\ee
We make the change of variables
$|k|=E+{\beta}^{2}v$, and recall the definition
$u={\Pi}_{+}(E)$. Then the integral can be written
\be{int.L0+L1}
\int_{-\infty}^{\infty}|L_{0}(w,\beta,t;k) +
L_{1}(w,\beta,t;k)|^{2} d k =
2 {\beta}^{2} \int_{-E/{\beta}^{2}}^{\infty}
|F(w,\beta,t;v)|^{2}  d v
\ee
where
\be{def.F}\eqalign{
F(w,\beta,t;v) = & e^{-it(z_{0}-E+{\beta}^{2}v)}
{\kappa({\beta}^{2}w) \over {\beta}^{2}v+{\beta}^{2}w} \cr
 + & 2 {\beta}^{2} \kappa({\beta}^{2}w){f}(E)
f(E+{\beta}^{2}v)
J(-{\beta}^{2}u+{\beta}^{2}v, {\beta}^{2}v-{\beta}^{2}w,
 -{\beta}^{2}u -{\beta}^{2}w;t),\cr}
\ee
and we have removed an overall phase factor
$e^{-2iEt}$. By rescaling the integration variable
in the definition of $J$, we can rewrite this as
\be{rewrite.F}
F(w,\beta,t;v) = {\kappa({\beta}^{2}w) \over {\beta}^{2}}
\bigg[{e^{-it({\beta}^{2}v+z_{0}-E)} \over v+w} +
2 {f}(E) f(E+{\beta}^{2}v) J(-u+v,v-w,-u-w;
{\beta}^{2}t) \bigg]
\ee
Let $\tau = {\beta}^{2}t$; evaluating $J$ gives
\be{eval.J}
{-1 \over 2 \pi i}J(-u+v,v-w,-u-w;\tau) =
{e^{-i \tau (v-u)} \over (-u+w)(v+w)} +
{e^{-i \tau (v-w)} \over (u-w)(v+u)} +
{e^{-i \tau (-u-w)} \over (-w-v)(-u-v)}
\ee
Substituting this in $\num{int.L0+L1}$ we see that the integral
over $v$ is dominated by the contributions from the
neighborhoods of the poles. Therefore we can
replace $f(E+{\beta}^{2}v)$ by $f(E)$ with an error
O(${\beta}^{2}$), and we can replace the lower limit
of integration $-E/{\beta}^{2}$ by $-\infty$ with the same
error.
Finally, replacing $e^{-it({\beta}^{2}v+z_{0}-E)}$ by
$e^{-i{\beta}^{2}t(v-u)}$ produces an error which is
${\rm O}({\beta}^{4}t) = {\rm O}({\beta}^{2}\tau)$.

Therefore since $\Gamma = 2 \pi |f(E)|^{2}$ and using
$\lim_{\beta \rightarrow 0} {\beta}^{-2}
|\kappa({\beta}^{2}w)|^{2} = {\rm Im}(w)/2 \pi$, we get
\be{est15}\eqalign{
\lim_{\beta \rightarrow 0}
||P_{1}\big(e^{-iH_{1}(\beta)({\beta}^{-2}\tau)}
& {\Psi}^{({\beta}^{2}w)}\big)||^{2} \cr
& =
\lim_{\beta \rightarrow 0} 2 {\beta}^{2} \int_{-E/{\beta}^{2}}
^{\infty} |F(w,\beta,{\beta}^{-2}\tau;v)|^{2} d v \cr
& =
2 \Big({{\rm Im}w \over 2 \pi}\Big) \int_{-\infty}^{\infty}
 \bigg| {e^{-i\tau (v-u)} \over v+w} + {\Gamma \over \pi}
J(-u+v,v-w,-u-w;\tau) \bigg|^{2} d v \cr}
\ee
After some rearrangement this produces $\num{result2}$.
\endproof

\medskip
\noindent{\it 4.5 Proof of $\thmnum{bound T0,1}$.}
\smallskip

Suppose that $\psi$ in $\num{def.B}$, $\num{eval.T1.B}$ is dilation analytic.
Then the first term on the right side of $\num{rotate.T1}$ is absent, and hence
$\num{est10}$ immediately implies the bound
\be{est25}
\Big|T_{1}({\psi}^{({\beta}^{2}w)},t)(k)\Big|
\leq C {\beta}^{2} |f(k)|,
\ee
which together with $\num{bound.T0-L0}$ gives $\num{boundforT0,1}$.
\endproof

\medskip
\noindent{\it 4.6 Evaluation of limit.}
\smallskip

Here we write down the explicit expression for the integral
on the right side of $\num{result2}$. For convenience, define
\be{def.alpha}
\alpha = {u-w \over 2i\Gamma}
\ee
Then Cauchy's theorem gives the following:
\be{eval.int.K}\eqalign{
{1 \over 2 \pi i} \int_{-\infty}^{\infty}
\Big| K(u,w;\tau,\alpha;v) \Big|^{2} d v = &
{e^{i \tau (u -{\bar u})} \over w -{\bar w}}\big[|1+\alpha|^{2}-
e^{i \tau (w -{\bar w})}(1+\alpha+{\bar \alpha})\big] \cr
& + {e^{i \tau (w -{\bar w})} \over u -{\bar u}}\big[1 -
e^{i \tau (u -{\bar u})}\big] \cr
& - {(1+{\bar \alpha})e^{i \tau (w-{\bar u})}
\over u -{\bar w}}\big[1 -
e^{i \tau (u - {\bar w})}\big] \cr
& - {(1+\alpha)e^{i \tau (u-{\bar w})}
\over w -{\bar u}}\big[1 -
e^{i \tau (w - {\bar u})}\big] \cr}
\ee





\vfill\eject


\noindent{\bf 5. \underbar{Convergence of the expansion.}}

\medskip
\noindent{\it 5.1 Proof of $\thmnum{bound TR}$.}
\smallskip

Our first goal is to prove $\thmnum{bound TR}$ which shows that
$T_{R}(\psi,t)$ in $\num{expand.in.T}$ is uniformly bounded in $t$. Using
Stone's
formula $\num{def.T*}$ this follows immediately from the following
estimate (recall the definition of $G_{R}$ in $\num{def.GR}$,
and the definition of $\delta$ in section 2.3 (iv)).

\medskip
\beginlemmalabel{bound GR}
For $|\beta|$ sufficiently small, and for all $z$ with
${\rm Im}z \neq 0$,
\be{bound.GR}
||G_{R}(\beta,z)|| \leq {C {\beta}^{\{2,\delta\}}
\over ({\rm Re}z-E)^{2} + 1}
\ee
\endlemmalabel
\medskip
We will deduce this result from $\thmnum{bound A}$ below.
For ${\rm Im}z \neq 0$,
and $p,q$ real, define
\be{def.A}
{\cal A}(z;p,q) = \la f,(h_{1}+|p|-z)^{-1}G(\beta,z)
(h_{1}+|q|-z)^{-1}f \ra
\ee
Then the integral kernel of the operator
$T(z)G(\beta,z)T(z)$ is given as follows:
\be{kernel1}
\big(T(z)G(\beta,z)T(z)\big)(p,q) = f(p) {\cal A}(z;p,q) {f}(q)
\ee
Therefore the integral kernel of the operator $G_{R}(\beta,z)$
in $\num{def.GR}$ is the following:
\be{kernel2}
(G_{R}(\beta,z))(p,q) = {\beta}^{4}
g(z-|p|)^{-1} f(p) {\cal A}(z;p,q) g(z-|q|)^{-1} {f}(q)
\ee
Our main technical result is the following.

\medskip
\beginlemmalabel{bound A}
For $|\beta|$ sufficiently small, for all $z$ with
${\rm Im}z \neq 0$, and for all $p,q$,
\be{bound.A}
|{\cal A}(z;p,q)| \leq C {\beta}^{-2} {\beta}^{\{2,\delta\}}
\ee
\endlemmalabel

\medskip

\noindent{\it Proof of $\thmnum{bound GR}$.}\par
\noindent We use the Hibert-Schmidt norm to bound the operator norm:
\be{est16}
||G_{R}(\beta,z)||^{2} \leq \int \int
|(G_{R}(\beta,z))(p,q)|^{2}  d p d q
\ee
Using $\thmnum{bound A}$ and also the bounds $\num{bound2.g}$ we deduce
that for ${\rm Im}z \neq 0$
\be{est17}\eqalign{
& ||G_{R}(\beta,z)||^{2} \cr
& \leq C {\beta}^{4} ({\beta}^{\{2,\delta\}})^{2}\int \int
|f(p)|^{2} |f(q)|^{2}
\big(|z-|p|-E|^{2} + {\beta}^{4} \big)^{-1}
\big(|z-|q|-E|^{2} + {\beta}^{4} \big)^{-1}  d p d q \cr}
\ee
Recall from section 2.3(v) that $|f(k)|^{2} \leq C (k^{2}+1)^{-1}$;
combined with the estimate
\be{est18}
\int_{-\infty}^{\infty} \big((|z-|p|-E)^{2} + {\beta}^{4}\big)^{-1}
(p^{2}+1)^{-1} d p \leq {C {\beta}^{-2} \over
({\rm Re}z-E)^{2}+1}
\ee
\be{est19}
|f(k)|^{2} \leq C (k^{2}+1)^{-1},
\ee
we deduce the result.
\endproof

So the proof reduces to $\thmnum{bound A}$ for the kernel ${\cal A}$.
Our proof of this uses dilation analyticity in an essential way.
First we dilate the operator $T(z)$ defined in $\num{def.T}$,
using the prescription $\num{dil.op}$.
We assume below that ${\rm Im}z>0$; the proof for ${\rm Im}z<0$
is the same.
If ${\rm Im}z >0$, then
the operator $T(z)$ is dilation analytic for $\theta$ satisfying
$0 \leq {\rm Re}\theta \leq {\theta}_{0}$, and for these values
\be{def.T.theta}
(T_{\theta}(z)\psi)(k) =
f_{\theta}(k) \la f_{-{\bar \theta}},
(h_{1,\theta} + e^{-i\theta}|k| - z)^{-1} \psi \ra.
\ee
Furthermore since the integral kernel of this operator is
\be{kernel3}
(T_{\theta}(z))(p,q) =
{f_{\theta}(p) \overline{f_{-{\bar \theta}}(q)} \over
e^{-i\theta}|p| + e^{-i\theta}|q| - z}
\ee
we deduce that for ${\rm Im}z>0$ and
${\theta}_{0} \geq {\rm Re}\theta \geq 0$,
\be{est19}
|(T_{\theta}(z))(p,q)| \leq C
|f_{\theta}(p)| |f_{-{\bar \theta}}(q)|
\min \Big({1 \over {\rm Im}z},
{1 \over {\rm Re}(\theta)\sqrt{|p|\,|q|}}\Big)
\ee
Combining this with $\num{bound.G.theta}$ produces the bound:
\be{est20}
\bigg|\bigg( G_{0,\theta}^{1/2} T_{\theta}(z)
G_{0,\theta}^{1/2} \bigg) (p,q) \bigg| \leq
C \,|f_{\theta}(p)| |f_{-{\bar \theta}}(q)|
\min \Big({1 \over ({\rm Im}z)^{2}},
{1 \over ({\rm Re}(\theta))^{2}|p|\,|q|}\Big)
\ee
Using the Hilbert-Schmidt norm to bound the operator norm,
and using the fact that $|k|^{-r} |f_{\theta}(k)|^{2}$
is integrable for $0 \leq r \leq 2$, we get
\be{est21}
\bigg|\bigg|\bigg( G_{0,\theta}^{1/2} T_{\theta}(z)
G_{0,\theta}^{1/2} \bigg)\bigg|\bigg| \leq
C \min \Big({1 \over ({\rm Im}z)^{2}},
{1 \over ({\rm Re}(\theta))^{2}}\Big)
\ee


We now define the dilated version of the operator
$G(\beta,z)$, by the following series:
\be{def.series}
G_{\theta}(\beta,z) = G_{0,\theta}(\beta,z)
\sum_{N=0}^{\infty} \bigg( {\beta}^{2} T_{\theta}(z)
G_{0,\theta}(\beta,z) \bigg)^{N}
\ee
For $|\beta|$ sufficiently small, the bound $\num{est21}$
shows that (i) this series converges uniformly
in $z$ and $\theta$ for ${\rm Im}z>1$ and
$0 \leq {\rm Re}\theta < {\theta}_{0}$, and
(ii) if ${\rm Re}\theta$ is bounded away from zero, then
this series converges uniformly in $z$ for ${\rm Im}z>0$.
Hence $G_{\theta}(\beta,z)$ defined by $\num{def.series}$
is analytic in $\theta$, $z$ in both of these regions.

Returning to ${\cal A}(z)$ in $\num{def.A}$,
it follows that for ${\rm Im}z>1$
and for $0 \leq {\rm Re}\theta \leq {\theta}_{0}$,
\be{rotate.A}
{\cal A}(z;p,q) = \la f_{-{\bar \theta}}, (h_{1,\theta}+|p|-z)^{-1}
G_{\theta}(\beta,z) (h_{1,\theta}+|q|-z)^{-1} f_{\theta} \ra
\ee
The independence of $\theta$ follows as usual by noting that
the right side is analytic in $\theta$, and is independent of
$\theta$ when ${\rm Re}\theta = 0$.

Now fix $\theta={\theta}_{0}$. Since $G(\beta,z)$ is a projection of the
resolvent of a self-adjoint operator, ${\cal A}(z;p,q)$ is an
analytic function of $z$ for all ${\rm Im}z>0$. As noted above,
the right side of $\num{rotate.A}$ is also analytic for ${\rm Im}z>0$.
Thus by the uniqueness of analytic continuation,
\be{continue.A}
{\cal A}(z;p,q) = \la f_{-{\theta}_{0}}, (h_{1,{\theta}_{0}}+|p|-z)^{-1}
G_{{\theta}_{0}}(\beta,z)
(h_{1,{\theta}_{0}}+|q|-z)^{-1} f_{{\theta}_{0}} \ra
\ee
for all $z$ with ${\rm Im}z>0$.

We now bound ${\cal A}(z)$.
>From $\num{est21}$ and $\num{def.series}$ it follows that
\be{est22}\eqalign{
|{\cal A}&(z;p,q)|\cr
&\leq C  \|\big(G_{0,-{\theta}_{0}}(\beta,\bar z)\big)^{1\over2}
(h_{1,-{\theta}_{0}}+|p|-\bar z)^{-1}f_{-{\theta}_{0}}\|
\|\big(G_{0,{\theta}_{0}}(\beta,z)\big)^{1\over2}
(h_{1,{\theta}_{0}}+|q|-z)^{-1}f_{{\theta}_{0}}\|.
}\ee
Further, using $\num{bound5.g}$ and section 2.3(iv),
\be{est23}\eqalign{
\|\big(G_{0,{{\theta}_{0}}}(\beta,z)\big)^{1\over2}
(h_{1,{{\theta}_{0}}}+&|q|-z)^{-1}
f_{{\theta}_{0}}\|^2\cr
&=\int {|f_{{\theta}_{0}}(p)|^2\over \big|e^{-i{{\theta}_{0}}}|p|+|q|-z\big|^2
\big|e^{-i{{\theta}_{0}}}|p|+E-z-\beta^2\Pi(z-e^{-i{{\theta}_{0}}}|p|)\big|}\, d
p\cr
&\le C \int_{0}^{1} {|p|^{1+\delta} \over |p|^{2} (p + {\beta}^{2})}\, dp +
C \cr
& \leq C {\beta}^{-2} {\beta}^{\{2,\delta\}} \cr
}\ee
The same estimate holds for
$\|\big(G^{(1)}_{0,-{{\theta}_{0}}}
(\beta,\bar z)\big)^{1\over2}(h_{1,-{{\theta}_{0}}}+|p|-\bar z)^{-1}
f_{-{{\theta}_{0}}}\|^2$.
The bound claimed in $\num{bound.A}$ now follows.\endproof\medskip


\vfill\eject


\noindent{\bf \underbar{Bibliography.}}
\medskip



{\noindent\frenchspacing

\item{[AC]} J. Aguilar, J. M. Combes, {\it A class of analytic perturbations
for one-body Schr\"odinger Hamiltonians},
Commun. Math. Phys. {\bf 22}, 269-279 (1971).\bigskip

\item{[A1]} A. Arai, {\it Selfadjointness and spectrum of Hamiltonians
in nonrelativistic quantum electrodynamics},
Journal of Mathematical Physics {\bf 22}, 534-537 (1981).\bigskip

\item{[A2]} A. Arai, {\it On a model of a harmonic oscillator coupled
to a quantized, massless scalar field},
Journal of Mathematical Physics {\bf 22}, 2539-2548, 2549-2552 (1981).\bigskip

\item{[BFS1]} V. Bach, J. Fr\"ohlich, and I. M. Sigal,
{\it Mathematical theory of nonrelativistic matter and
radiation}, Lett. Math. Phys. {\bf 34},  183-201 (1995).\bigskip

\item{[BFS2]} V. Bach, J. Fr\"ohlich, and I. M. Sigal,
{\it Quantum Electrodynamics of Confined Non-Relativistic
Particles}, preprint mp-arc 96-346 (1996).\bigskip

\item{[BC]} E. Balslev, J. M. Combes, {\it Spectral properties of
many-body Schr\"odinger operators with dilatation analytic interactions},
Commun. Math. Phys. {\bf 22}, 280-294 (1971).\bigskip

\item{[CDG]} C. Cohen-Tannoudji, J. Dupunt-Roc, G. Grynberg,
{\it Photons and Atoms -- Introduction to Quantum Electrodynamics},
John Wiley, New  York (1991).\bigskip

\item{[Da]} E.B. Davies, {\it Dynamics of a Multilevel Wigner-Weisskopf Atom},
Journal of Mathematical Physics {\bf 15}, 2036-2041 (1974).\bigskip

\item{[D]} P.A.M. Dirac, {\it Quantum Mechanics},
Oxford University Press (1978).\bigskip

\item{[DE]} J. Dittrich, P. Exner, {\it A non-relativistic
model of two-particle decay I. Galilean invariance}, Czech.
Jour. Phys. {\bf B37}, 503-515 (1987); {\it II. Reduced
resolvent}, Czech. Jour. Phys. {\bf B37}, 1028-1034 (1987);
{\it III. The pole approximation},  Czech. Jour. Phys. {\bf B38},
591-610 (1988); {\it IV. Relation to the scattering theory,
spectral concentration, and bound states},  Czech. Jour. Phys.
{\bf B39}, 121-138 (1989).\bigskip

\item{[EY]} J. H. Eberly, H. Yoo, {\it Dynamical
theory of an atom with two
or three levels interacting with quantized cavity fields},
Physics Reports {\bf 118}, 239 (1985).\bigskip

\item{[Fri]} K.O. Friedrichs, {\it On the perturbation
of continuous spectrum}, Comm. Pure Appl. Math {\bf 1}, 361-406
(1948).\bigskip

\item{[Fro]} J Fr\"ohlich, {\it On the infrared problem in a model
of scalar electrons and massless, scalar bosons},
Ann. Inst. Henri Poincare {\bf 19}, 1-103 (1973).\bigskip

\item{[GT]} A. Grossmann, A. Tip, {\it Hydrogen atoms interacting with a
quantised radiation mode},
Jour. Phys. A {\bf 13}, 3381-3397 (1980).\bigskip

\item{[Ho]} J. Howland, {\it Spectral Concentration and Virtual Poles, II},
Transactions of the A.M.S. {\bf 162}, 141-156 (1971);\hfill\break
{\it The Livsic Matrix in Perturbation Theory}, Journal of Mathematical
Analysis and Applications {\bf 50}, 415-437  (1975).\bigskip

\item{[HSp]} M. H\"ubner, H. Spohn, {\it Radiative decay:
nonperturbative approaches},
Rev. Math. Phys. {\bf 7}, 363-387 (1995).\bigskip

\item{[Hu]} W. Hunziker, {\it Resonances,
Metastable States and Exponential Decay Laws in Perturbation Theory},
Commun. Math. Phys. {\bf 132}, 177-188 (1990).\bigskip

\item{[JP]} V. Jaksic, C. A. Pillet, {\it On a model for quantum friction
I: Fermi's golden rule and dynamics at zero temperature},
Ann. Inst. Henri Poincare {\bf 62}, 47 (1995).\bigskip

\item{[JC]} E. T. Jaynes, F. W. Cummings, {\it },
Proc. IEEE {\bf 51}, 98 (1963).\bigskip

\item{[K1]} C. King, {\it Exponential Decay Near Resonance,
without Analyticity}, Letters in Math. Phys. {\bf 23},
215-222 (1991).\bigskip

\item{[K2]} C. King, {\it Scattering Theory for a Model of an Atom in a
Quantized Field},
Letters in Math. Phys. {\bf 25}, 17-28 (1992).\bigskip

\item{[OY]} T. Okamoto, K. Yajima, {\it Complex scaling technique in
non-relativistic massive QED},
Ann. Inst. Henri Poincare {\bf 42}, 311-327 (1985).\bigskip

\item{[Pf]} P. Pfeifer, {\it Radiative decay versus photonic Jahn-Teller
distortion of molecular states},
Phys. Rev. A {\bf 26}, 701-704 (1982).\bigskip

\item{[RS]} M. Reed and B. Simon, {\it Methods of Mathematical
Physics vol. 1: Functional Analysis},
Academic Press (1980).\bigskip

\item{[Si]} B. Simon, {\it Resonances in N-Body Quantum
Systems with Dilation Analytic Potentials and the Foundations of
Time-Dependent Perturbation Theory}, Annals of Math. {\bf 97},
247-274 (1973).\bigskip

\item{[Sp1]} H. Spohn, {\it Ground state(s) of the spin-boson Hamiltonian},
Communications in Mathematical Physics {\bf 123}, 277-304 (1989).\bigskip

\item{[Sp2]} H. Spohn, {\it Ground state of a quantum particle coupled to
a scalar Bose field},
preprint mp-arc 97-174 (1997).\bigskip

\item{[St]} S. Stenholm, {\it Quantum theory of electromagnetic fields
interacting with atoms and molecules},
Physics Reports {\bf 6}, 1-122 (1973).\bigskip

\item{[W]} D. F. Walls, {\it Higher order effects in the single atom field
mode approximation},
Phys. Lett. A {\bf 42}, 217-218 (1972).\bigskip

\item{[WW]} V. Weisskopf and E. Wigner,
{\it Berechnung der Nat\"urlichen Linienbreite auf
Grund der Diracshen Lichttheorie}, Zeitschrift f\"ur Physik
{\bf 63}, 54-73 (1930).\bigskip}




\end
