%PARAM.1
\magnification=\magstep1
\pageno=0
%FONTS.2
\font\huge=cmbx10 scaled\magstep2
\font\xbold=cmbx10 scaled\magstep2
\font\cbold=cmsy10 scaled\magstep2
\font\mbold=cmmi10 scaled\magstep2
\font\sub=cmbx10 scaled\magstep1
\font\small=cmr7
%TITLES.2
\count5=0
\count6=1
\count7=1
\count8=1
\def\proof{\medskip\noindent{\bf Proof.\ }}
\def\qed{\hfill\smallskip
	 \line{\hfill\vrule height 1.8ex width 2ex depth +.2ex
	       \ \ \ \ \ \ }\bigskip}
\def\remark{\medskip\noindent{\bf Remark.\ }}
\def\abstract#1\par{\noindent{\sub Abstract.} #1 \par}
\def\equ(#1){\csname e#1\endcsname}
\def\clm(#1){\csname c#1\endcsname}
\def\equation(#1){\eqno\tag(#1)}
\def\tag(#1){(\number\count5.
	      \number\count6)
    \expandafter\xdef\csname e#1\endcsname{
    (\number\count5.\number\count6)}
    \global\advance\count6 by 1}
\def\claim #1(#2) #3\par{
    \vskip.1in\medbreak\noindent
   {\bf #1\ \number\count5.\number\count7.}
    \expandafter\xdef\csname c#2\endcsname{
    #1~\number\count5.\number\count7}
    \global\advance\count7 by 1
    {\sl #3}\par
    \ifdim\lastskip<\medskipamount
    \removelastskip\penalty55\medskip\fi}
\def\section#1\par{\vskip0pt plus.3\vsize\penalty-75
    \vskip0pt plus -.3\vsize\bigskip\bigskip
    \global\advance\count5 by 1
    \message{#1}\leftline
     {\xbold \number\count5.\ #1}
    \count6=1
    \count7=1
    \count8=1
    \nobreak\smallskip\noindent}
\def\subsection#1\par{\vskip0pt plus.2\vsize\penalty-75
    \vskip0pt plus -.2\vsize\bigskip\bigskip
    \message{#1}\leftline{\sub
    \number\count5.\number\count8.\ #1}
    \global\advance\count8 by 1
    \nobreak\smallskip\noindent}
%MACROS.11
\footline{\ifnum\pageno=0\hss\else\hss\tenrm\folio\hss\fi}
\let\cl=\centerline
\def\mean{{\rm I\kern-.18em E}}
\def\natural{{\rm I\kern-.18em N}}
\def\integer{{\rm Z\kern-.32em Z}}
\def\real{{\rm I\kern-.2em R}}
\def\complex{\kern.1em{\raise.47ex\hbox{
	    $\scriptscriptstyle |$}}\kern-.40em{\rm C}}
\def\I{|\hskip-0.08em |\hskip-0.08em |}
\def\Id{{\rm I}}
\def\Int{\int\limits_{-\infty}^\infty\!}
\def\HT{_{\scriptscriptstyle HT}}
\def\UV{_{\scriptscriptstyle UV}}
\def\IR{_{\scriptscriptstyle IR}}
\def\AA{{\cal A}}
\def\BB{{\cal B}}
\def\CC{{\cal C}}
\def\DD{{\cal D}}
\def\EE{{\cal E}}
\def\FF{{\cal F}}
\def\GG{{\cal G}}
\def\HH{{\cal H}}
\def\II{{\cal I}}
\def\JJ{{\cal J}}
\def\KK{{\cal K}}
\def\LL{{\cal L}}
\def\MM{{\cal M}}
\def\NN{{\cal N}}
\def\OO{{\cal O}}
\def\PP{{\cal P}}
\def\QQ{{\cal Q}}
\def\RR{{\cal R}}
\def\SS{{\cal S}}
\def\TT{{\cal T}}
\def\UU{{\cal U}}
\def\VV{{\cal V}}
\def\WW{{\cal W}}
\def\XX{{\cal X}}
\def\YY{{\cal Y}}
\def\ZZ{{\cal Z}}
%ABSTRACT.2
\vskip.4in
\cl{{\huge A Nontrivial Renormalization Group Fixed Point}}
\cl{{\huge for the Dyson--Baker Hierarchical Model}}
\vskip.6in
\cl{Hans Koch
\footnote{$^1$}
{{\small Supported in Part by the
National Science Foundation under Grant No. DMS--9103590.}}}
\cl{Department of Mathematics, University of Texas at Austin}
\cl{Austin, TX 78712}
\vskip.25in
\cl{Peter Wittwer
\footnote{$^2$}
{{\small Supported in Part by the
Swiss National Science Foundation.}}}
\cl{D\'epartement de Physique Th\'eorique, Universit\'e de Gen\`eve}
\cl{Gen\`eve, CH 1211}
\vskip.8in
\abstract
We prove the existence of a nontrivial Renormalization Group (RG) fixed point
for the Dyson--Baker hierarchical model in $d=3$ dimensions.
The single spin distribution of the fixed point is shown to be
entire analytic, and bounded by $\exp(-{\rm const}\times t^6)$
for large real values of the spin $t$.
Our proof is based on estimates for the zeros of a 
RG fixed point for Gallavotti's hierarchical model.
We also present some general results for the heat flow
on a space of entire functions, including an order preserving
property for zeros, which is used in the RG analysis.

\par\vfill\eject

\def\rDi{1}
\def\rW{2}
\def\rBak{3}
\def\rMore{4}
\def\rG{5}
\def\rKWi{6}
\def\rKWii{7}
\def\rKWiii{8}
\def\rBo{9}
%INTRO.4
\section Introduction and Main Results

One of the basic assumptions in the modern theory of critical phenomena
is the existence of nontrivial renormalization group (RG) fixed points,
associated with certain universality classes of interactions.
Within the framework of statistical mechanics,
this assumption is actually a conjecture,
and it should be possible to either prove it or disprove it.
However, even for the simplest classes of ``realistic'' interactions,
such as the one represented by the three--dimensional Ising model,
the rigorous construction of a nontrivial RG fixed point
seems beyond the reach of presently known methods.
In addition, there is a lack of good numerical results in this area,
which indicates that even at a quite fundamental level
there are still gaps in our understanding of RG transformations.

The traditional approach in such situations is to try
to first understand some simpler,
and thus necessarily less realistic, class of interactions.
In this case, the model with the longest history
is Dyson's hierarchical model [\rDi--\rMore].
For the Dyson--type hierarchical analogue
of the nearest--neighbor (continuous spin) Ising model in $d=3$ dimensions,
which was first considered by Baker [\rBak],
the full RG transformation $\TT$ reduces to the third power of the following
nonlinear operator $\RR$.
$$
\bigl(\RR(h)\bigr)(t)=K\Int dx\
e^{-2\sigma x^2}h(\alpha t+x)h(\alpha t-x)\,,\qquad t\in\real\,.
\equation(1)
$$
Here, $h$ is the density of the single spin distribution,
$K$ and $\sigma$ are arbitrary but fixed positive real numbers, and 
$$
\alpha=2^{-{d-2+\eta\over 2d}}=2^{-{1\over 6}},\qquad (d=3\,,\eta=0)\,.
\equation(2)
$$

\noindent
More details about the connection between the transformations $\TT$
and $\RR$ will be given in Section 5.
Our main result is the following.

\claim Theorem(1) $\RR$ has a fixed point which is
the restriction to $\real$ of an even entire analytic function $h\IR\,$.
Given any $\gamma>6$, the function $h\IR$ satisfies the bounds
$$
0 < h\IR(t) < c_1 e^{-c_2 t^6}\,,\qquad\forall t\in\real\,,
\equation(3)
$$
$$
|h\IR(t)|< c_3 e^{c_4 |t|^\gamma}\,,\qquad\forall t\in\complex\,.
\equation(3b)
$$
for some positive constants $c_1,\ldots,c_4\,$.

In what follows, we choose the values $K=2\alpha$ and
$\sigma=(2\alpha^2-1)/(4\alpha^2-1)$ for the normalization
constants that appear in the definition of $\RR$.
This can be done without loss of generality,
since any two transformations of the type\equ(1)
are conjugate via a scaling of the form $h\mapsto ah(b.)$.

Let now $h\IR$ be the fixed point of $\RR$ described in\clm(1).
Then we can define an entire analytic function $f\IR$ by the equation
$$
h\IR(t)=e^{\sigma t^2}\,{\sigma\over\pi}\Int ds\ 
e^{-2i\sigma st}f\IR(is)\,,\qquad t\in\real .
\equation(4)
$$
A straightforward calculation shows that $f\IR$ satisfies
the fixed point equation for the following transformation $\NN$.
$$
\bigl(\NN(f)\bigr)(t)
={1\over\sqrt{(1-\beta^2)\pi}}
\Int ds\;e^{-{1\over 1-\beta^2}s^2}f(\beta t+s)^2\,,\qquad t\in\complex\,,
\equation(5)
$$
with $\beta={1\over 2\alpha}=2^{-5/6}$.
Conversely, if $f\IR$ is some given fixed point for $\NN$, one can try
to use equation\equ(4) in order to obtain a fixed point $h\IR$ for $\RR$.
This is exactly the strategy which we have adopted here.
The main problem is to estimate the function $f\IR$ along the imaginary axis.
The following lemma, together with equation\equ(4), implies\clm(1).

\claim Lemma(2) $\NN$ has a fixed point $f\IR$ with the following properties.
$f\IR$ is entire analytic, $f\IR(i.)$ is a function of positive type
(i.e., the Fourier transform of a finite positive measure),
and for every positive $q<6/5$ there are constants $b_1,\ldots,b_4>0$
such that
$$
\;|f\IR(t)| < b_1 e^{b_2 |t|^{6/5}}\,,\qquad\forall t\in\complex\,,
\equation(7)
$$
$$
|f\IR(it)| < b_3 e^{-b_4 |t|^q}\,,\qquad\forall t\in\real\,.
\equation(6)
$$

We note that the same result,
but without the bound\equ(6), was already obtained in [\rKWii].
To prove only that part of\clm(2),
the analysis can be restricted, modulo a finite dimensional problem,
to a neighborhood of the high temperature fixed point $f\HT\equiv 1$.
The result is sufficient to show that there exists a nontrivial
``weak'' solution of the fixed point equation for $\RR$.
But in order to establish the regularity of this solution,
as implied by the bound\equ(6), it becomes necessary to analyze
the transformation $\NN$ on a more {\it global} scale.
Roughly speaking,
the infrared fixed point $h\IR$ of $\RR$ inherits its asymptotic behavior
(fast decay) from the high temperature fixed point $h\HT(t)=\delta(t)$,
and its local regularity
from the ultraviolet fixed point $h\UV={\rm const}$.
This is indicated e.g. by the fact that the fixed point $f\IR$
is ``close'' to $f\HT$ only with respect to its relatively
slow growth; its behavior\equ(6) along the imaginary axis
is qualitatively much closer to that of the ultraviolet fixed point
$f\UV(t)=\exp(\sigma t^2+{\rm const})$.
This may be regarded as a consequence of the generally believed fact
(which we think could be proved with the methods developed here)
that the ultraviolet fixed point lies on the boundary
of the stable manifold of the infrared fixed point.

As was already shown in [\rKWii], the transformation $\NN$ has
the following interesting property $P$:
If all zeros of a polynomial $f$ lie on the imaginary axis,
then the same is true for $\NN(f)$.
In addition, we have the bound\equ(7) on the growth of $f\IR$.
Thus, since it seems likely that there exists at least one polynomial
(e.g. close to the ultraviolet fixed point) which has only imaginary zeros,
and which lies on the stable manifold of the fixed point $f\IR$,
it is natural to conjecture that the infrared fixed point of $\NN$
can be written as a canonical product
$$
f\IR(t)=f\IR(0)\prod_{k=0}^\infty
\Bigl[1+\Bigl({t\over \nu_k}\Bigr)^{\!2}\,\Bigr]
\,,\qquad t\in\complex\,,
\equation(8)
$$
with purely imaginary zeros $\pm i\nu_0,\pm i\nu_1,\ldots$.
This conjecture is indeed correct, but it appears very hard to prove.
In fact, our first attempt contained an error [\rKWii, Lemma 4.2].
A new proof, which follows the argument given above
(not for $\NN$, but for a contraction with the same property $P$
and the same fixed point $f\IR$)
can be found in [\rKWiii].

\claim Theorem(3) $\NN$ has a fixed point $f\IR$ of the form\equ(8),
with $0<\nu_0<\nu_1<\ldots$ and $f\IR(0)>0$.
Furthermore, given any $p<{1\over 6}$,
there are positive constants $a_1$ and $a_2$ such that for $k=1,2,\ldots$
$$
\nu_k > a_1 k^{5/6}\,,\qquad
\nu_k-\nu_{k-1} < a_2 k^{-p}\,.
\equation(9)
$$

This theorem implies\clm(2),
even though for a general product of the form\equ(8), the bounds\equ(9)
are not sufficient to guarantee decay along the imaginary axis.
In order to get\equ(6),
we also need to use that $f\IR$ is a fixed point of $\NN$.

Our proof of\clm(3) uses input (see Theorem 3.2)
from a computer--assisted analysis, which is described in detail in [\rKWiii].
In particular, it requires upper bounds on the first $80$ gaps
$\gamma_k=\nu_k-\nu_{k-1}$ for the fixed point $f\IR$ of $\NN$.
Given these bounds,
we use an order preserving property (for zeros) of the heat flow,
and estimates on the evolution of polynomial approximations
to the hyperbolic functions, in order to inductively estimate
the gaps $\gamma_{2k}$ and $\gamma_{2k+1}$
in terms of the gaps near $\nu_k$.
This method is suggested by a connection between the heat flow
and the evolution of a system of interacting charges on a line.

The results stated so far will be proved in Section 3.
We shall now describe some general results on the one--dimensional
heat flow, which we use in our RG analysis,
but which may be of independent interest.
The proofs and further details will be given in Section 4.

Denote by $\EE$ the set of all entire analytic functions $f$ that satisfy
$$
\mathop{\overline{\rm lim}}_{|z|\to\infty}|z|^{-2}\ln |f(z)| =0\,,
\equation(10)
$$
and whose zeros all lie on a straight line $\mu(f)+i\real$
parallel to the imaginary axis.
For every $f\in\EE$ and for every $\lambda>0$,
define a function $H_\lambda f$ by the equation
$$
\bigl(H_\lambda f\bigr)(z)={1\over\sqrt{4\pi\lambda}}\Int ds\,
e^{-{1\over 4\lambda}s^2}f(z-s)\,,\qquad z\in\complex\,,
\equation(11)
$$
and let $H_0f=f$.

\claim Lemma(4) Let $f\in\EE$. Then for all $\lambda>0\,$,
the function $H_\lambda f$ lies in $\EE$ and has no multiple zeros.
The set of solutions $(\lambda,z)$ of the equation $(H_\lambda f)(z)=0$
is a union of curves $\lambda\mapsto z(\lambda)$
which are continuous on $[0,\infty)$ and real analytic on $(0,\infty)$.

In what follows, we will assume that the zeros of a function $f\in\EE$
have been arranged into an indexed set $\{\mu(f)+i\nu_k(f)\}_{k\in I}\,$,
where multiple zeros are repeated according to their multiplicities.
We also assume that $I$ is a set of consecutive integers,
and that the zeros are labeled in such a way that the gaps
$$
\gamma_k(f)=\cases{
\nu_k(f)-\nu_{k-1}(f)\,,\quad &if $\{k-1,k\}\subset I$ ;\cr
+\infty\,, &otherwise ;\cr}
\equation(12)
$$
are nonnegative, for all integers $k$.
When indexed this way, the zeros of $f$ determine
a double--sided sequence $\gamma(f)$ of numbers in $[0,\infty]$.
On the set of such sequences,
one has the following canonical partial order relation.

\claim Definition(5)
$\gamma\prec\gamma^\prime\quad\Longleftrightarrow\quad
\gamma^{\phantom\prime}_k\le\gamma^\prime_k\,,\ \forall k\in\integer.$

\noindent
We note that by\clm(4),
if $f$ is a function in $\EE$ with zeros indexed by $I$,
then there exists a canonical (unique) way of indexing the zeros
of the entire family $\{H_\lambda f\}_{\lambda\ge 0}$ by the same set $I$,
in such a way that the curves $\lambda\mapsto\nu_k(H_\lambda f)$
are continuous for all $k\in I$.

\claim Theorem(6) Let $f$ and $g$ be functions in $\EE$ with indexed zeros.
If the functions $H_\lambda f$ and $H_\lambda g$ are indexed canonically,
then $\gamma(f)\prec\gamma(g)$ implies
$\gamma(H_\lambda f)\prec\gamma(H_\lambda g)\,$, for all $\lambda>0$.

\medskip\noindent
{\bf Remarks.}

\smallskip\noindent$\bullet$\hskip1.5em
In [\rKWiii] it is shown that the fixed point $f\IR$ described here
is unique in a small open subset $\OO$ of some Banach space.
Previously, in [\rKWii, Section 3], we proved the existence
of a nontrivial fixed point for $\NN$ in some other set $G\oplus H$.
We have no doubt that $\OO$ is contained in $G\oplus H$,
but we did not try to prove this,
since we have no uniqueness result on $G\oplus H$.
On the other hand, the fixed point properties derived in [\rKWii, Section 4]
hold for any non--constant fixed point of $\NN$ in a function space
that is much larger than (and contains) the spaces considered here.
Thus, all these results apply to $f\IR$.
This includes the following stronger version of the estimate\equ(7).
If we define $\ell_n(t^2)=(\beta^{-n}t)^{-6/5}\ln f\IR(\beta^{-n}t)$,
then the limit $\ell(x)=\lim_{n\to\infty}\ell_n(x)$ exists and is nonzero
for all $x\in\real_+$, and the function $s\mapsto\ell(e^s)$ is periodic
with period $\ln(\beta^{-2})$. In fact, it can be proved that
the same holds if $\real_+$ is replaced by $\complex\setminus (-\infty,0]$.

\smallskip\noindent$\bullet$\hskip1.5em
We would like to stress that the hierarchical model discussed here
is {\it not} a mean--field model,
and that the infinite--volume measure associated with the function $h\IR$
is an {\it exact} fixed point for the Wilson--Kadanoff transformation,
without any approximations.
The ``only'' flaw of this model is the lack of translation invariance
and its consequences (no anomalous dimension).
Apart from that, the model fully supports the general RG picture,
down to the numerical value of the critical index of the free energy
($\nu=0.6495704\ldots$ from non--rigorous numerical computation),
which differs only by $2\%$ from what is believed to be the correct
value for translation--invariant models such as the Ising model
in $d=3$ dimensions ($\nu\approx 0.638$).

%PRODUCTS.2
\section Canonical Products and {\cbold N}

In this section, we will use\clm(6) in order to show that
the gaps for an $\NN$--invariant canonical product of the form\equ(8)
decrease like $\nu_k-\nu_{k-1}={\scriptstyle\OO}(k^{-p})$,
if a finite number of these gaps satisfy a certain bound.

In addition to the convolution operators $H_\lambda$ defined in
equation\equ(11), consider now also the map
$Q_\beta: f\mapsto f^2(\beta .)$.
Let $f\in\EE$, and assume that the zeros $\omega_k(f)=\mu(f)+i\nu_k(f)$
of $f$ have been indexed by a set $I$ of consecutive integers
(the intersection of $\integer$ with a connected subset of $\real$),
such that the gaps defined by equation\equ(12) are non--negative.
Then we will index the zeros of $Q_\beta(f)$ by setting
$$
\omega_{2k}\bigl(Q_\beta(f)\bigr)=\omega_{2k+1}\bigl(Q_\beta(f)\bigr)
=\beta^{-1}\omega_k(f)\,,\qquad\forall k\in I.
\equation(201)
$$
With this convention, it is clear that $Q_\beta$ is order--preserving,
in the same sense as the operators $H_\lambda$.
Thus, since $\NN=H_\lambda\circ Q_\beta\,$ for $\lambda=(\beta^{-2}-1)/4$,
we obtain the following as a consequence of\clm(4) and\clm(6).

\claim Corollary(201) Let $f$ and $g$ be functions in $\EE$ with indexed
zeros. Then $\NN(f)$ and $\NN(g)$ are also functions in $\EE$, and if their
zeros (which are all simple) are indexed as described above, then
$\gamma(f)\prec\gamma(g)$ implies
$\gamma\bigl(\NN(f)\bigr)\prec\gamma\bigl(\NN(g)\bigr)$.

In order to take advantage of this property of $\NN$, we need
a class of functions $g$ such that the gaps for $\NN(g)$ can be estimated.
The following functions $h^\pm_\ell\,$, and their translates,
turn out to be ideal for this purpose.

\claim Definition(202) For any fixed $\kappa>0$,
define $h^+=\cosh(\kappa\pi .)$, $h^-=\sinh(\kappa\pi .)$,
and define $h_0^\pm, h_1^\pm,\ldots$
by the equation 
$$
\eqalign{
h^+_\ell(z)&=\prod_{k=0}^\ell\left[1+{\kappa^2z^2\over(k+1/2)^2}\right]\,,
\qquad z\in\complex,\cr
h^-_\ell(z)&=\kappa\pi z\prod_{k=1}^{\ell+1}
\left[1+{\kappa^2z^2\over k^2}\right]\,,\ \ 
\qquad z\in\complex.\cr}
\equation(202)
$$

Notice that $h^\pm_\ell\to h^\pm$ pointwise, as $\ell$ tends to infinity.
The image under $\NN$ of the limit can be computed explicitly:
$$
\bigl(\NN(h^\pm)\bigr)(z)
={\textstyle{1\over 2}}e^{(1-\beta^2)\kappa^2\pi^2}
\left[\cosh(2\kappa\pi\beta z)\pm e^{-(1-\beta^2)\kappa^2\pi^2}\right]\,.
\equation(203)
$$
The important observation here is that
for large values of $\kappa$, the gaps between the zeros of
$\NN(h^\pm)$ are approximately a factor of $2\beta=2^{1/6}$
smaller than the gaps between the zeros of $h^\pm$.
The following proposition will be used to show
that the same is true for the ``middle'' gaps for $\NN(h^\pm_\ell)$,
if $\ell$ is sufficiently large. Define
$$
a(\kappa)=(1-\beta^2)
\exp\!\left({\textstyle{\beta^2\over 4\kappa^2(1-\beta^2)}}\right)
\left[\kappa^{-2}+2(1-\beta^2)\pi^2+\kappa^{-2}
e^{-(1-\beta^2)\kappa^2\pi^2}\right]\,.
\equation(204)
$$

\claim Proposition(203) Let $\kappa>0$ and $\ell\in\natural$.
If $t$ is a real number that satisfies $|t|\le{1\over 2\kappa}$, then
$$
\left\vert\bigl(\NN(h^\pm)\bigr)(it)-\bigl(\NN(h^\pm_\ell)\bigr)(it)\right\vert
\le{a(\kappa)\kappa^4\over\ell+1}\cdot{\textstyle{1\over 2}}
e^{(1-\beta^2)\kappa^2\pi^2}.
\equation(205)
$$

\proof Let $\sigma=``+$''.
We estimate first the difference between the function
$\phi_\ell=[h^\sigma_\ell/h^\sigma]^2$ and the constant function $1$.
For real values of $s$, we obtain
$$
\eqalign{
\left\vert\phi_\ell^\prime(s)\right\vert
&=\left\vert {d\over ds}\prod_{k=\ell+1}^\infty
\left[1+{\kappa^2s^2\over(k+1/2)^2}\right]^{-2} \right\vert\cr
&=\left\vert\sum_{k=\ell+1}^\infty{-4\kappa^2s\over(k+1/2)^2}
\left[1+{\kappa^2s^2\over(k+1/2)^2}\right]^{-1}\phi_\ell(s)\right\vert\cr
&\le 4\kappa^2|s|\sum_{k=\ell+1}^\infty{1\over(k+1/2)^2}
\le{4\kappa^2|s|\over\ell+1}\,,\cr}
\equation(206)
$$
and thus
$$
\left\vert\phi_\ell(s)-1\right\vert
=\left\vert\int\limits_0^sdt\,\phi_\ell^\prime(t)\right\vert
\le{4\kappa^2\over\ell+1}\int\limits_0^{|s|}dt\,t
={2\kappa^2s^2\over\ell+1}\,.
\equation(207)
$$
By using this last bound, the left hand side of\equ(205)
can be estimated as follows.
$$
\eqalign{
{\rm l.h.s.}\,
&=\left\vert{1\over\sqrt{(1-\beta^2)\pi}}\Int ds\,
e^{-{1\over 1-\beta^2}(s-i\beta t)^2}
\bigl(h^\sigma(s)\bigr)^2\bigl[1-\phi_\ell(s)\bigr]\right\vert\cr
&\le{1\over\sqrt{(1-\beta^2)\pi}}\Int ds\,
\left\vert e^{-{1\over 1-\beta^2}(s-i\beta t)^2}\right\vert
\bigl(\cosh(\kappa\pi s)\bigr)^2\,{2\kappa^2s^2\over\ell+1}\cr
&={1\over 2}e^{(1-\beta^2)\kappa^2\pi^2}
(1-\beta^2){\kappa^2\over\ell+1}e^{\beta^2t^2\over 1-\beta^2}
\left[1+2(1-\beta^2)\kappa^2\pi^2+e^{-(1-\beta^2)\kappa^2\pi^2}\right].\cr}
\equation(208)
$$
After substituting the maximum value ${1\over 2\kappa}$
for $|t|$, we get precisely the bound given in\equ(205).
The proof for the case $\sigma=``-$'' is similar.
\qed

The resulting bound on the middle gaps for $\NN(h^\pm_\ell)$
will now be formulated more generally, for functions that have
a sufficiently long sequence of gaps of size $\le\kappa^{-1}$.

\claim Definition(204) If $g$ is a function in $\EE$
with indexed zeros and gaps $\gamma_j(g)$,
we define for every $k\in\integer$ and $\ell\in\natural$,
$$
\hat\gamma_k(g,\ell)=\max\{\gamma_j(g): k-\ell\le j\le k+\ell+1\}\,.
\equation(209)
$$

\claim Definition(205) Define $\UU$ to be the set of all triples
$(\kappa,\theta,\ell)$ in
$\real_+\!\times\!\real_+\!\times\!\natural$ such that
$$
{1\over 2\beta}+\theta<1\,,\qquad
e^{-(1-\beta^2)\kappa^2\pi^2}+{a(\kappa)\kappa^4\over\ell+1}
< \sin(\pi\beta\theta)\,.
\equation(210)
$$

\claim Proposition(206) Let $(\kappa,\theta,\ell)\in\UU$.
Then for every $f\in\EE$ with indexed zeros, and for every $k\in\integer$,
the following holds.
$$
\hat\gamma_k(f,\ell)\le\kappa^{-1}\quad\Longrightarrow\quad
\hat\gamma_{2k}(\NN(f),0)
\le\bigl({\textstyle{1\over 2\beta}}+\theta\bigr)\kappa^{-1}\,.
\equation(211)
$$

\proof Consider a fixed $k\in\integer$. Then, for any given $\ell\in\natural$,
we index the zeros of $h^+_\ell$ and $h^-_\ell$
in increasing order of their imaginary parts,
by successive integers
ranging from $k-\ell-1$ up to $k+\ell$ and $k+\ell+1$, respectively.
Define $t={1\over 2\kappa}\bigl({1\over 2\beta}+\theta\bigr)$.

We claim that
under the given assumptions,
$$
\gamma_{2k}\bigl(\NN(h^+_\ell)\bigr)\le 2t\,,\qquad
\gamma_{2k+1}\bigl(\NN(h^-_\ell)\bigr)\le 2t\,.
\equation(212)
$$
Notice that the gaps considered in\equ(212)
are the ones between the two zeros of $h^\pm_\ell$ that lie
closest to the origin.
By assumption, we have 
$$
\cos(2\kappa\pi\beta t)+e^{-(1-\beta^2)\kappa^2\pi^2}
\le -{a(\kappa)\kappa^4\over\ell+1}\,,
\equation(213)
$$
and $|t|\le{1\over 2\kappa}$.
Thus, it follows from equation\equ(203) and\clm(203) that
$$
\bigl(\NN(h^\pm)\bigr)(it)
\le -{a(\kappa)\kappa^4\over\ell+1}\cdot{\textstyle{1\over 2}}
e^{(1-\beta^2)\kappa^2\pi^2}
\le\bigl(\NN(h^\pm)\bigr)(it)-\bigl(\NN(h^\pm_\ell)\bigr)(it)\,,
\equation(214)
$$
which implies that $\bigl(\NN(h^\pm_\ell)\bigr)(it)\le 0$.
Since the function $\bigl(\NN(h^\pm_\ell)\bigr)(i.)$ is even
and takes a positive value at the origin,
we conclude that it has at least two zeros in the interval $[-t,t]$.
This proves the bounds\equ(212).

Assume now that $\hat\gamma_k(f,\ell)\le\kappa^{-1}$.
Then, since the (finite) gaps for the functions $h^\pm_\ell$
are all of size $\kappa^{-1}$, we have $\gamma(f)\prec\gamma(h^\pm_\ell)$.
But by\clm(201), this implies that
$\gamma_j\bigl(\NN(f)\bigr)\le\gamma_j\bigl(\NN(h^\pm_\ell)\bigr)$ for all $j$.
Thus, the assertion follows from the bounds\equ(212).
\qed

The statement\equ(211) will now be used to inductively estimate
the gaps for a fixed point $f$ of $\NN$. Given the hypotheses of
the lemma below, we can assume that the zeros $i\nu_k(f)$
of $f$ are indexed by $\integer$,
in such a way that
$$
\nu_{k-1}(f)\le\nu_k(f)=-\nu_{-k-1}(f)\,,\qquad \forall k\in\integer.
\equation(215)
$$

\claim Lemma(207) Let $f$ be an a non--constant even function
in $\EE$ that satisfies the fixed point equation for $\NN$,
and let $(\kappa_0,\theta,\ell_0)$ be an element of $\,\UU$. If
$$
\gamma_k(f)<\kappa_0^{-1}\,,\qquad n_0-\ell_0\le k<2n_0\,,
\equation(216)
$$
for some $n_0>\ell_0+1$, then for every positive $p<1/6$
there exists a constant $b$ such that
$$
\gamma_k(f)<b k^{-p}\,,\qquad k=1,2\ldots\ .
\equation(217)
$$

\proof Assuming that the hypotheses of this lemma are satisfied,
we shall first show that the conclusion holds
for some (as opposed to all) positive $p<1/6$.

Depending on a choice of $p$, define two sequences
$n\mapsto\kappa(n)$ and $n\mapsto\ell(n)$ as follows.
$$
\kappa(n)=\kappa_0\left({n\over n_0-\ell_0}\right)^p\,,\qquad
\ell(n)=\left\lfloor\ell_0\left({n\over n_0}\right)^{5/6}\right\rfloor\,,
\equation(218)
$$
where $\lfloor r\rfloor$ denotes the integer part of a real number $r$.
Below we will prove that
$$
\bigl(\kappa\bigl(n-\ell(n)\bigr),\theta,\ell(n)\bigr)\in\UU
\equation(219)
$$
for every $n\ge n_0\,$, provided that $p>0$ is sufficiently small.
Our goal now is to show that for arbitrary $n\ge n_0\,$, if
$$
\gamma_j(f)<\kappa(j)^{-1}\,,\qquad n-\ell(n)\le j<2n\,,
\equation(220)
$$
then
$$
\gamma_j(f)<\kappa(j)^{-1}\,,\qquad 2n\le j\le 2n+1\,.
\equation(221)
$$
This is sufficient to prove the assertion,
since for $n=n_0\,$, the bound\equ(220) follows from the assumption\equ(216),
if $p>0$ is sufficiently small.

Let now $n\ge n_0$ be fixed, and assume that\equ(220) holds.
Then $\hat\gamma_n\bigl(f,\ell(n)\bigr)$ is less than or equal to
$\kappa\bigl(n-\ell(n)\bigr)^{-1}$,
and we obtain the bound
$$
\hat\gamma_{2n}(f,0)\le\bigl({\textstyle{1\over 2\beta}}+\theta\bigr)
\kappa\bigl(n-\ell(n)\bigr)^{-1}
\equation(222)
$$
from\clm(206). Thus, in order to prove\equ(221), it suffices to show that
$$
{1\over 2\beta}+\theta<{\kappa(n-\ell(n))\over\kappa(2n+1)}\,,
\qquad\forall n\ge n_0\,.
\equation(223)
$$
To this end, choose an arbitrary positive real number $\varepsilon$.
Then for every $n\ge n_0\,$, we have
$$
\eqalign{
{\kappa(n-\ell(n))\over\kappa(2n+1)}
&=\left({n-\ell(n)\over 2n+1}\right)^p
=2^{-p}\left(1-{1\over 2n+1}\right)^p \left(1-{\ell(n)\over n}\right)^p\cr
&\ge 2^{-p}\left(1-{1\over 2n_0+1}\right)^p
\left(1-{\ell_0\over n_0}\right)^p > 2^{-p}(1-\varepsilon)\,,\cr}
\equation(224)
$$
if $p>0$ is sufficiently small.
Thus, since the left hand side of\equ(223) is less than $1$ by assumption,
and since $\varepsilon$ was arbitrary, it follows that\equ(223) holds
for sufficiently small $p>0$. This shows that\equ(220) implies\equ(221).

In order to prove\equ(219), choose again $\varepsilon>0$,
and denote by $n_1$ the largest value of $n$ for which
$\ell_0(n/n_0)^{1/2}<\ell_0+1$.
Since the inequalities\equ(210) in the definition of $\UU$ are strict,
the condition\equ(219) holds for $n_0\le n<n_1\,$,
if $p$ is sufficiently small. And if $n\ge n_1$, we have the bound
$$
\eqalign{
{\kappa(n-\ell(n))^4\over\ell(n)+1}
&\le{\kappa_0^4\over\ell_0+1}\cdot
{\ell_0+1\over(\ell(n)+1)(n/n_0)^{-4p}}
\left({1\over 1-\ell_0/n_0}\right)^{4p}\cr
&\le{\kappa_0^4\over\ell_0+1}\cdot
{\ell_0+1\over\ell_0(n/n_0)^{5/6}(n/n_0)^{-4p}}
\left({1\over 1-\ell_0/n_0}\right)^{4p}\cr
&\le{\kappa_0^4\over\ell_0+1}
\left({1\over 1-\ell_0/n_0}\right)^{4p}
\le{\kappa_0^4\over\ell_0+1}
(1+\varepsilon)\,,\cr}
\equation(225)
$$
for sufficiently small $p>0$.
Thus, since the factor $a(\kappa)$ in\equ(210)
is a decreasing function of $\kappa$, we find that the condition
\equ(219) also holds for $n\ge n_1\,$, if $p$ is sufficiently small.

After having proved\equ(217) for some $p>0$,
we will now choose new values for the parameters $(\kappa_0,\theta,\ell_0)$,
and then proceed as above,
but with ``fixed $p$ and sufficiently large $n_0$''
instead of ``fixed $n_0$ and sufficiently small $p$''.

Given some positive $p<1/6$, let $\theta$ be a positive real number satisfying
$$
{\textstyle{1\over 2\beta}}+\theta<2^{-p}\,.
\equation(226)
$$
Then choose $\kappa_0$ and $\ell_0$ such that
$(\kappa_0,\theta,\ell_0-1)\in\UU$, and define two sequences
$n\mapsto\kappa(n)$ and $n\mapsto\ell(n)$ as in\equ(218).
Given any $\varepsilon>0$, a bound similar to\equ(225)
can be obtained for all $n\ge n_0\,$, if $n_0$ is chosen sufficiently large:
$$
\eqalign{
{\kappa(n-\ell(n))^4\over\ell(n)+1}
&\le{\kappa_0^4\over\ell_0}\cdot
{\ell_0\over\ell_0(n/n_0)^{5/6}(n/n_0)^{-4p}}
\left({1\over 1-\ell_0/n_0}\right)^{4p}\cr
&\le{\kappa_0^4\over\ell_0}
\left({1\over 1-\ell_0/n_0}\right)^{4p}
\le{\kappa_0^4\over\ell_0}\,(1+\varepsilon)\,.\cr}
\equation(227)
$$
Thus, if $n_0$ is sufficiently large,
we have again\equ(219) for all $n\ge n_0\,$.
Under the same condition, the bound\equ(220) holds in the case $n=n_0\,$,
since $\gamma_j(f)\to 0$ as $j\to\infty$.
And as shown before,\equ(220) implies\equ(221),
provided e.g. that\equ(223) is satisfied.
But for sufficiently small $\varepsilon>0$,
the inequality\equ(223) follows from\equ(226) and\equ(224);
and the latter holds for all $n\ge n_0\,$,
if $n_0$ is chosen sufficiently large.
\qed

%PROOF.2
\section Proof of Theorems 1.1 and 1.3

The proofs given in this section use \clm(207),
together with the following two theorems.
The first theorem is a reformulation of two results from [\rKWii, Section 4].

\claim Theorem(501)
Let $f\not\equiv 0$ be an even entire analytic
fixed point for $\NN$, whose restriction to $\real$ is real--valued,
and whose Taylor coefficients (at the origin) are bounded
in absolute value by those of the function $z\mapsto K\exp({1\over r}z^2)$,
for some constants $K>0$ and $r>(4\alpha^2-1)/(2\alpha^2-1)$.
Then $f$ is the Fourier transform of a positive measure on $\real$
whose moments are all finite,
and there are constants $b_1$ and $b_2$ such that
$$
|f(z)|\le b_1e^{b_2|z|^{6/5}}\,,\qquad\forall z\in\complex.
\equation(501)
$$

\claim Theorem(502) The transformation $\NN$ has a fixed point $f^*$
with the following properties.
\item{$(a)$} $f^*$ is an even entire analytic function
which takes real values when restricted to $\real$.
\item{$(b)$} The Taylor coefficients at zero of $f^*$
are bounded by the corresponding coefficients
of the function $z\mapsto K\exp\bigl({1\over 10}z^2\bigr)$,
for some positive constant $K$.
\item{$(c)$} All zeros of $f^*$ of lie on the imaginary axis.
\item{ } \hskip-\parindent In addition, there are positive real numbers
$y_0<y_1<\ldots<y_{79}\,$, such that
\item{$(d)$} $f^*(iy_k)=0$ for $k=0,1,\ldots,79\,$.
\item{$(e)$} the function $f^*(i.)$ has exactly $5$ zeros
in the interval $[0,y_4]$.
\item{$(f)$} $y_k-y_{k-1}<{4\over 3}\,$, for $k=5,6,\ldots,79$.

\noindent A proof of\clm(502) is given in [\rKWiii].
In addition, we verify in [\rKWiii] that
$$
\bigl({\textstyle{3\over 4}},{\textstyle{1\over 10}},35\bigr)\in\UU\,,
\equation(502)
$$
where $\UU$ is the set specified in\clm(205).

\bigskip\noindent
{\bf Proof of Theorem 1.3.\ }
First, we note that $(4\alpha^2-1)/(2\alpha^2-1)<4$.

Let $f\IR=f^*$,
where $f^*$ is the fixed point of $\NN$ described in\clm(502).
Since the properties $(a)$ and $(b)$ of this fixed point
(and $(e)$, which implies that $f\IR\not\equiv 0$) verify the assumptions
of\clm(501), we obtain for $f=f\IR$ the bound\equ(501)
and the inequality $f(0)>0$.
Thus, given that $f\IR$ is even and satisfies $(c)$ as well,
it follows from Hadamard's factorization theorem that $f\IR$ can be
represented as a convergent canonical product of the form\equ(8).
The lower bound on $\nu_k$ in\equ(9) is obtained from\equ(501),
by using the standard inequality which bounds the number of zeros of
an entire function $f$ in the disk $|z|\le r$
by the logarithm of the maximum modulus of $f$ on the disk $|z|\le er$,
if $f(0)=1$.

Finally, if we set $(\kappa_0,\theta,\ell_0)=
\bigl({\textstyle{3\over 4}},{\textstyle{1\over 10}},35\bigr)$
and $n_0=40$, then the last three statements in\clm(502),
together with\equ(502), imply that the hypotheses of\clm(207)
are satisfied. But the conclusion of this lemma is precisely
the upper bound on the gaps $\gamma_k=\nu_k-\nu_{k-1}$ in\equ(9).
Thus,\clm(3) is proved.
\qed

In order to show that the function $f\IR$ decreases along the imaginary axis
as claimed in\equ(6), we need the following simple fact.

\claim Proposition(503) If $h$ is a ${\rm C}^n$ function
that has $n>0$ zeros in an interval $[0,x]$, then
$$
|h(x)|\le{x^n\over n!}\sup_{t\in[0,x]}\left\vert h^{(n)}(t)\right\vert\,.
\equation(503)
$$

\proof Let $t_0=x_{-1}=x$. For $k=0,1,\ldots,n-1$,
if we define $x_k$ to be the largest zero of $h^{(k)}$ in $[0,x_{k-1}]$,
then $h^{(k+1)}$ has at least $n-k-1$ zeros in $[0,x_k]$. Thus,
$$
|h(t_0)|=\left\vert\,
\int\limits_{x_0}^{t_0}\!dt_1\!
\int\limits_{x_1}^{t_1}\!dt_2\,\cdots\!
\int\limits_{x_{n-1}}^{t_{n-1}}\!dt_n\,h^{(n)}(t_n)\right\vert
\le\int\limits_0^{t_0}\!dt_1\!
\int\limits_0^{t_1}\!dt_2\,\cdots\!
\int\limits_0^{t_{n-1}}\!dt_n\,
\sup_{t\in[0,x]}\left\vert h^{(n)}(t)\right\vert\,,
$$
and the assertion follows.
\qed

\bigskip\noindent
{\bf Proof of Lemma 1.2.\ }
Let $f\IR$ be the fixed point of $\NN$ described in\clm(3).
Then the hypothesis (and hence the conclusion) of\clm(501)
is clearly satisfied for the function $f=f\IR\,$.
Thus, the only thing that remains to be proved is the bound\equ(6).

Given that $f\IR(i.)$ is the Fourier transform of some finite
measure $\mu$, we have $|f\IR(it)|\le|f\IR(0)|$ for every $t\in\real$.
The same inequality holds for all even derivatives of $f\IR\,$,
since the moments of $\mu$ are finite. Thus, by Cauchy's formula
and\equ(501), there are constants $d_1,d_2>0$ such that
$$
\left\vert{1\over n!}f\IR^{(n)}(it)\right\vert
\le\left\vert{1\over n!}f\IR^{(n)}(0)\right\vert
\le d_1(d_2n)^{-5n/6}\,,
\qquad t\in\real,\ n\in 2\natural.
\equation(504)
$$

Let us now choose an arbitrary $q\in(1,{6\over 5})$.
Let $q-{1\over 5}<r<1$, and define $p=(q-r)/(1+q-r)$.
Since $p<1/6$, it follows from\clm(3)
that $\nu_k\le d_3 k^{1-p}+d_4$ for some positive constants $d_3,d_4\,$.
By combining this inequality with the second bound in\equ(9),
we find $d_5>0$ such that
if $t$ is sufficiently large and $\nu_k>t-t^r$, then
$$
\nu_k-\nu_{k-1}
\le a_2\left({\nu_k-d_4\over d_3}\right)^{-p/(1-p)}
\le{t^{-q+r}\over 2d_5}\,.
\equation(505)
$$
Thus, the function $f\IR(i.)$
has $n\ge d_5\,t^q$ zeros in the interval $[t-t^r,t]$ for large $t$.
Now we can apply\clm(503) and the bound\equ(504) in order
to estimate the value of the function $h=f\IR(i(t-t^r+.))$ at $x=t^r$.
The result is that
$$
|f\IR(it)|\le d_6\exp\bigl(-d_7\,t^q\ln t\bigr)\,,\qquad t>0,
\equation(506)
$$
for some positive constants $d_6$ and $d_7\,$.
\qed

\bigskip\noindent
{\bf Proof of Theorem 1.1.\ }
Let $f\IR$ be the function described in\clm(2). By analytic continuation,
we can write the fixed point property of $f\IR$ in the form
$$
{1\over\sqrt{(1-\beta^2)\pi}}
\Int ds\;e^{-{1\over 1-\beta^2}s^2}f\IR\bigl(i\beta^{-1}(t+s)\bigr)
=f\IR(it)^2\,,\qquad t\in\complex\,.
\equation(507)
$$
To be more precise, consider the operator $S_i$ which maps an entire
function $f$ to the function $f(i.)$. A group property of the heat
flow $\lambda\mapsto H_\lambda\,$, which will be proved in the next section,
implies that $S_i^{-1}H_\lambda S_i H_\lambda f=f$,
for every entire function $f$ that satisfies\equ(10).
Thus, given the bound\equ(7), it follows that\equ(507)
is equivalent to the equation $\NN(f\IR)=f\IR\,$.

Let now $h\IR$ be the function defined by equation\equ(4).
From\equ(6) it is clear that this function is entire analytic,
and that it satisfies the bound\equ(3b).
In order to prove\equ(3),
consider the Fourier transform $\psi$ of the function $t\mapsto f\IR(it)^2$.
Below we will show that $\psi$ satisfies the bound
$$
|\psi(t)|\le d_1 e^{-d_2\,t^6}\,,\qquad t\in\real,
\equation(508)
$$
for some constants $d_1,d_2>0$. By equation\equ(507),
this leads to an analogous estimate for the Fourier transform of $f\IR\,$,
from which the upper bound in\equ(3) follows.
Given now that $\RR(h\IR)$ is well defined, it is easy to check
that the function $h\IR$ satisfies the fixed point equation for $\RR$.
In addition, for every $t\in\real$
we have $h\IR(t)\ge 0$, since $f\IR$ is a function of positive type;
and $h\IR(t)=0$ is excluded by the fact that $h\IR$ is a fixed point for $\RR$.

We will now turn to the proof of the bound\equ(508).
From the fact that $f\IR$ is the Fourier transform of a positive measure,
it follows that the function $\psi$
takes only non--negative values on $\real$. In addition,
since $\psi$ and $-\psi^{\prime\prime}$ are themselves
functions of positive type, we have the bounds
$$
|\psi(t)|\le|\psi(0)|=b^2\,,\qquad
|\psi^{\prime\prime}(t)|\le|\psi^{\prime\prime}(0)|=2c^2\,,
\equation(509)
$$
for all $t\in\real$. The constants $b,c>0$ are defined by\equ(509).
Let now $t$ be an arbitrary real number larger than $2b/c$,
and define $a=\psi(t)^{1/2}$. Then the even moments of $\psi$
can be bounded from below as follows.
$$
\eqalign{
\Int dx\,x^{2n}\psi(x)
&\ge\Bigl(t-{a\over c}\Bigr)^{2n}
\!\int\limits_0^{a/c}dx\,\bigl[\psi(t+x)+\psi(t-x)]\cr
&\ge\Bigl({t\over 2}\Bigr)^{2n}
\int\limits_0^{a/c}dx\,2\bigl[a^2-c^2x^2]
=\Bigl({t\over 2}\Bigr)^{2n}\,{4a^3\over 3c}\,.\cr}
\equation(510)
$$
On the other hand, we have a bound on the Taylor coefficients of
the function $f\IR^2$, which follows from\equ(7), and which implies that
$$
\Int dx\,x^{2n}\psi(x)
=\left\vert\bigl(f\IR^2\bigr)^{(2n)}(0)\right\vert
\le\left(d_3 n^{1/6}\right)^{2n}\,,\qquad n=1,2,\ldots\;,
\equation(511)
$$
for some constant $d_3>0$.
The two inequalities\equ(510) and\equ(511)
can now be combined to yield an upper bound on $\psi(t)$:
For any $\delta>0$ with the property that $n=\delta t^6$
is a positive integer, we get
$$
{4\over 3c}\psi(t)^{3/2}
={4a^3\over 3c}\le\left(2d_3\delta^{1/6}\right)^{2\delta t^6}.
\equation(512)
$$
Thus, if $\delta$ is chosen appropriately, the bound\equ(508) follows.
This completes the proof of\clm(1).
\qed

%HEAT.2
\section The Heat Flow on {\cbold E}

Denote by $\FF$ the vector space of all entire analytic functions $f$
that satisfy the bound\equ(10). On $\FF$ we define the following
directed family of seminorms.
$$
\| f\|_\rho=\sup_{z\in\complex}e^{-|z|^2/\rho}|f(z)|\,,\qquad\rho>0\,.
\equation(801)
$$
Equipped with the natural topology defined by these seminorms,
$\FF$ is a Fr\'echet space. We note that a linear operator $L$ on $\FF$
is continuous (bounded) if and only if for every $\rho>0$
there exist $c,r>0$ such that $\|Lf\|_\rho\le c\|f\|_r$
for every function $f\in\FF$.

Examples of continuous linear operators on $\FF$ are
differentiation $D: f\mapsto f^\prime$, translations
$T_\lambda: f\mapsto f(.-\lambda)\,$, and, as shown below,
the operators $K_{n,m}\,$, defined by the equation
$$
(K_{n,m}f)(z)={1\over 2\pi i}\int_{\Gamma_z}
{z^n f(\zeta) d\zeta\over \zeta^m(\zeta-z)}\,,
\qquad n\le m\,,\ z\in\complex\,.
\equation(802)
$$
Here, $\Gamma_z$ is a closed curve in $\complex$ with winding number $1$
with respect to the points $0$ and $z$. The proof for the continuity
of $D$ and $T_\lambda$ is similar to the proof of the following proposition.

\claim Proposition(801) Let $n\le m$ be nonnegative integers.
Then the operator $K_{n,m}$
is continuous on $\FF$, and for every $\rho>0$, it satisfies the bound
$$
\|K_{n,m}f\|_\rho\le 2^{-m+1}e^{1/\rho}\|f\|_{4\rho}\,,\qquad f\in\FF.
\equation(803)
$$

\proof Let $f\in\FF$ and $\rho>0$. Given any $z\in\complex$,
let $\Gamma_z$ be the positively oriented circle
of radius $2(|z|^2+1)^{1/2}$ around the origin.
Then for every $\zeta\in\Gamma_z$ we have $|z^n/\zeta^m|\le 2^{-m}$,
$|\zeta-z|\ge{1\over 2}|\zeta|\,$, and
$$
|f(\zeta)|\le \|f\|_{4\rho}e^{{1\over 4\rho}|\zeta|^2}
=\|f\|_{4\rho}e^{1/\rho+|z|^2/\rho}.
\equation(804)
$$
Thus, we can bound the right hand side of\equ(802)
by $\exp(|z|^2/\rho)$ times the right hand side of\equ(803), as claimed.
Since $f$ and $\rho$ were arbitrary, it follows that $K_{n,m}$ is continuous.
\qed

We shall now consider the heat flow $\lambda\mapsto H_\lambda$
defined in\equ(11), but extended to complex values of $\lambda$.

\claim Definition(802) For every $f\in\FF$ and for every $\lambda\in\complex$,
define a function $H_\lambda f$ by the equation
$$
\bigl(H_\lambda f\bigr)(z)=\hat f(\lambda,z)\equiv
{1\over\sqrt{4\pi}}\Int ds
e^{-s^2/4}f(z-\sqrt{\lambda}\,s)\,,\qquad z\in\complex .
\equation(805)
$$

\noindent
Notice that this definition is independent of the choice
of the square root function, and that $\hat f$ is analytic
in both of its arguments.

\claim Proposition(803) Let $\lambda\in\complex$, $\rho>0$, $f\in\FF$,
and define $r=\rho+4|\lambda|$.
Then $H_\lambda$ is a continuous linear operator on $\FF$,
and it satisfies the following two bounds.
$$\ \ \ \ 
\|H_\lambda f\|_\rho
\le(r/\rho)^{1/2}\|f\|_r\,,
\equation(806)
$$
$$
\|H_\lambda f-f-\lambda f^{\prime\prime}\|_\rho
\le{\textstyle{1\over 2}}|\lambda|^2(r/\rho)^{1/2}
\|f^{\prime\prime\prime\prime}\|_r\,.
\equation(807)
$$

\proof The inequality\equ(806) is obtained by replacing
$f(z-\sqrt{\lambda}\,s)$ in equation\equ(805)
by the bound $\|f\|_r\exp(|z-\sqrt{\lambda}\,s|^2/r)$,
and then computing the resulting Gaussian integral.
The second inequality can be obtained from the first one by using the identity
$$
\bigl(H_\lambda f\bigr)(z)
=f(z)+\lambda f^{\prime\prime}(z)+\lambda^2\!\int\limits_0^1\!dt\,
(1-t)\bigl(H_{t\lambda}f^{\prime\prime\prime\prime}\bigr)(z)\,,
\qquad z\in\complex .
\equation(808)
$$
\qed

The following proposition shows, among other things,
that the heat flow on the space $\FF$ is invertible.

\claim Proposition(804)
$H_{\lambda_1}H_{\lambda_2}f=H_{\lambda_1+\lambda_2}f\,$,
for every $f\in\FF$ and $\lambda_1\,,\lambda_2\in\complex$.

\proof Given any positive integer $n$, consider the vector space
$\PP_n$ of all polynomials of degree $\le n$, equipped with some norm.
On this space, differentiation $D$ is a bounded linear operator.
An explicit calculation shows that $\lambda\mapsto H_\lambda\,$,
$\lambda\ge 0$, is a semigroup of linear operators on $\FF$ which
leave $\PP_n$ invariant.
Thus, it follows from\equ(807) that
$H_\lambda f=\exp(\lambda D^2)f$, for every $f\in\PP_n$ and $\lambda\ge 0$.
By analyticity, the same holds for every $\lambda\in\complex$.
This proves the assertion in the case where $f$ is a polynomial.
But polynomials are dense in $\FF$. This follows e.g. from the fact
that $\Id-K_{n,n}$ projects $\FF$ onto $\PP_{n-1}$, and that
$K_{n,n}f$ converges to zero as $n\to\infty$, for every $f\in\FF$.
The latter is a consequence of\clm(801). Thus, the assertion is proved.
\qed

Let us now consider the $\lambda$--dependence of the zeros for $H_\lambda f$.
We start with properties that hold for every function $f\in\FF$.
The next proposition is an immediate consequence
of the implicit function theorem, given the fact the the function $\hat f$,
defined by equation\equ(805), is entire analytic.

\claim Proposition(805) Let $f\in\FF$
and $(\lambda_0,z_0)\in\complex\!\times\!\!\complex$.
If $H_{\lambda_0}f$ has a simple zero at $z_0\,$,
then in some open polydisk $U\times V$ containing $(\lambda_0,z_0)$,
the set of solutions $(\lambda,z)$ of the equation $(H_\lambda f)(z)=0$
is the graph of an analytic function from $U$ to $V$.

\claim Proposition(806) Let $m\ge 1$,
and denote by $h_m$ the $m^{\rm th}$ Hermite polynomial.
Then for every $\rho>0$ there exists a constant $K>0$
such that the following holds.
If $f$ is a function in $\FF$ that has a zero of order $m$ at the origin,
and if $\varepsilon$ is a complex number of modulus $\le 1$, then
$$
\left\Vert\hat f({\textstyle{\varepsilon^2\over 2}},\varepsilon .)
-{\textstyle{\varepsilon^m\over m!}}f^{(m)}(0)h_m\right\Vert_\rho
\le K|\varepsilon|^{m+1}\|f\|_{8\rho+16}\,.
\equation(809)
$$

\proof Let $|\varepsilon|\le 1$, and
for every $f\in\FF$ define $S_\varepsilon f=f(\varepsilon .)$.
In order to prove the bound\equ(809), we may
assume that $f$ satisfies $z^{-m}f(z)\to 1$ as $z\to 0$.
Then we can write
$$
\bigl(S_\varepsilon f\bigr)(z)=\varepsilon^mz^m+\varepsilon^{m+1}z^{m+1}
\bigl(K_{0,m+1}f\bigr)(\varepsilon z)\,,
\equation(810)
$$
where $K_{0,m+1}$ is the operator defined by the equation\equ(802).
For every $n\in\natural$ and $z\in\complex$, define $p_n(z)=z^n$.
An explicit calculation shows that $H_{1/2}p_m=h_m\,$.
Thus, by using\equ(810) and\clm(801), we obtain the bound
$$
\eqalign{
\Vert S_\varepsilon H_{\varepsilon^2/2}f-\varepsilon^mh_m\Vert_\rho
&=\Vert H_{1/2}S_\varepsilon f-\varepsilon^mh_m\Vert_\rho
=\left\Vert\varepsilon^{m+1}H_{1/2}
\bigl[p_{m+1}\cdot\bigl(S_\varepsilon K_{0,m+1}f\bigr)\bigr]\right\Vert_\rho\cr
&\le(1+2/\rho)^{1/2}|\varepsilon|^{m+1}\left\Vert p_{m+1}
\cdot\bigl(S_\varepsilon K_{0,m+1}f\bigr)\right\Vert_{\rho+2}\cr
&\le(1+2/\rho)^{1/2}|\varepsilon|^{m+1}
\|p_{m+1}\|_{2\rho+4}\Vert S_\varepsilon K_{0,m+1}f\Vert_{2\rho+4}\cr
&\le K|\varepsilon|^{m+1}\|f\|_{8\rho+16}\,,\cr}
\equation(811)
$$
for some positive constant $K$.
\qed

\claim Proposition(810) Let $\omega$ and $r$ be continuous functions
on some compact set $\Lambda\subset\complex$,
with values in $\complex$ and $[0,\infty)$, respectively.
Furthermore, let $f\in\FF$, and denote by $m(\lambda)$
the number of zeros of $H_\lambda f$
in the disk $|z-\omega(\lambda)|\le r(\lambda)$.
Here, and in what follows, zeros of order $m$ are counted as $m$ zeros.
Then, given any sufficiently small $\varepsilon>0$,
there exists an open neighborhood $U$ of $f$ in $\FF$,
such that for every $g\in U$ and for every $\lambda\in\Lambda$,
the function $H_\lambda g$ has exactly $m(\lambda)$ zeros
in the disk $|z-\omega(\lambda)|<r(\lambda)+\varepsilon$.

\proof For every $\delta>0$
define $U(\delta)=\{g\in\FF: \|f-g\|_\rho<\delta\}$,
where $\rho$ is the maximum value of $1+4|\lambda|$ on $\Lambda$.
The assertion follows by a standard argument, if we can show that
for every sufficiently small $\varepsilon>0$
there exists a $\delta>0$, such that a bound of the form
$$
\left\vert\bigl(H_\lambda f\bigr)(z)-\bigl(H_\lambda g\bigr)(z)\right\vert
< {\textstyle{1\over 2}}c_\varepsilon
< \left\vert\bigl(H_\lambda f\bigr)(z)\right\vert
\equation(826)
$$
holds on the circle $C_\varepsilon(\lambda)$
of radius $r(\lambda)+\varepsilon$ centered at $\omega(\lambda)$,
for every $g\in U(\delta)$ and for every $\lambda\in\Lambda$.
Consider first the second inequality in\equ(826).
Since the functions $\omega$ and $r$ are continuous, it follows that the set
$C_\varepsilon=\{(\lambda,z):
\lambda\in\Lambda,\, z\in C_\varepsilon(\lambda)\}$
is compact for every $\varepsilon\ge 0$.
Thus, the continuous function $|\hat f|$ has a minimum value $c_\varepsilon$
on $C_\varepsilon$, and $c_\varepsilon$ is positive for sufficiently small
values of $\varepsilon>0$. In order to get the first inequality in\equ(826),
we can use the bound
$$
\left\vert\bigl(H_\lambda f\bigr)(z)-\bigl(H_\lambda g\bigr)(z)\right\vert
\le\delta\sqrt{\rho}\,e^{|z|^2}\,,\qquad g\in U(\delta)\,,
\equation(827)
$$
which follows from\clm(803).
Since $|z|$ is uniformly bounded on the compact set $C_\varepsilon\,$,
it suffices to take $\delta$ small enough.
\qed

We shall now work towards a proof of\clm(6), in the case where
the functions involved have only a finite number of zeros.
Here, only real values of the heat flow parameter $\lambda$
will be considered, unless specified otherwise.

\claim Definition(807) Denote by $\EE_f$ the set of all functions in $\EE$
that have only a finite number of zeros, and denote by $\EE_p$ the
set of all polynomials whose zeros all lie on the imaginary axis.
In order to simplify notation, define $f_\lambda=H_\lambda f$.

\noindent
Let $f$ be a function in $\EE_f\,$. Then by a standard argument,
there exists a polynomial $g\in\EE_p$ such that
$$
f(z)=e^{bz}g(z-\mu)\,,\qquad z\in\complex ,
\equation(812)
$$
for some constants $b\in\complex$ and $\mu\in\real$.
A short calculation shows that
$$
f_\lambda(z)=e^{b^2\lambda}e^{bz}g_\lambda(z-\mu+2b\lambda)\,,
\qquad z\in\complex ,
\equation(813)
$$
which implies that the zeros sets for $f_\lambda$ and $g_\lambda$ are 
related to each other by a translation, for every $\lambda$.
Thus, if we prove the assertion of\clm(6) for polynomials in $\EE_p\,$,
then the same holds for every function in $\EE_f\,$.

Let now $g$ be a nonzero polynomial of degree $d\ge 2$.
Then the functions $g_\lambda$ are polynomials
of the same degree $d$, and we may write them in the canonical form
$$
g_\lambda(z)=a(\lambda)\prod_{k\in I}
\bigl[z-i\nu_k(g_\lambda)\bigr]\,,\qquad z\in\complex\,.
\equation(814)
$$
Here, the zeros $i\nu_k(g_\lambda)$ of $g_\lambda$
have been indexed by some set $I$ of $d$ consecutive integers,
in such a way that the curves $\lambda\mapsto\nu_k(g_\lambda)$ are continuous
(not necessarily real--valued, at this point).
That this is possible follows e.g.
from the implicit function theorem and\clm(806).
Assume now that the $k^{\rm th}$ zero of $g_\lambda$ is simple,
for some given value $\lambda_0$ of the parameter $\lambda$.
Then, by differentiating the equation
$0=g_\lambda(i\nu_k(g_\lambda))$ with respect to $\lambda$,
and using that ${d\over d\lambda}g_\lambda=g_\lambda^{\prime\prime}\,$,
we obtain the identity
$$
{d\over d\lambda}\,\nu_k(g_\lambda)
=i\,{g_\lambda^{\prime\prime}\bigl(i\nu_k(g_\lambda)\bigr)\over
g_\lambda^\prime\bigl(i\nu_k(g_\lambda)\bigr)}
=\sum_{j\in I\setminus\{k\}}{1\over \nu_k(g_\lambda)-\nu_j(g_\lambda)}\,,
\equation(815)
$$
for every $\lambda$ near $\lambda_0$.
Notice that if all numbers $\nu_k(g_\lambda)$
are real and distinct from each other at $\lambda=0$,
then the same remains true for all $\lambda>0$.

The following definition makes precise
what we mean (e.g. in \clm(6)) by a function with indexed zeros.

\claim Definition(808) A function in $\EE$ with indexed zeros
is a nonzero function $f\in\EE$, together with a pair $(I,\varphi)$,
consisting of a set $I$ of integers
and an onto map $\varphi$ from $I$ to $f^{-1}(0)$,
with the following properties:
For all $j$ and $k$ in $I$, if $j<k$ then $k-1\in I$
and ${\rm Im}\,\varphi(k-1)\le{\rm Im}\,\varphi(k)$;
and for every zero $z$ of $f$, the cardinality of $\varphi^{-1}(z)$
is equal to the order of $z$.

\noindent
{\bf Remark.\ }To simplify notation,
if a quantity $q$ depends on a function $f$ with indexed zeros $(I,\varphi)$,
we will write $q(f)$ instead of $q((f,(I,\varphi)))$.
In particular, the real and imaginary parts of $\varphi(k)$
will be denoted by $\mu(f)$ and $\nu_k(f)$, respectively.
We note that if $f$ has only a finite number of zeros, then
the function $\varphi$ is uniquely determined by the choice of $I$.
In addition, $I$ can be reconstructed from the corresponding sequence
of gaps, defined in equation\equ(12).

\claim Proposition(809) For every $\lambda\ge 0$, 
$H_\lambda$ maps $\EE_f$ to $\EE_f\,$, and the following holds.
Let $g$ and $h$ be functions in $\EE_f$ with indexed zeros,
and assume that the zeros of $H_\lambda g$ and $H_\lambda h$ have been
indexed by the same set of integers as the zeros of $g$ and $h$, respectively.
Then $\gamma(g)\prec\gamma(h)$
implies $\gamma\bigl(H_\lambda g\bigr)\prec\gamma\bigl(H_\lambda h\bigr)$.

\proof As was mentioned earlier, it suffices to prove the assertion
for polynomials $g,h\in\EE_p$.
We may also assume that $g$ is of degree $\ge2$. 
Consider first the case where all zeros of $g$ and $h$ are simple.
Then by equation\equ(815), $g_\lambda$ and $h_\lambda$
lie in $\EE_p\,$ for all $\lambda>0$.
From the same equation it follows that
$$
{d\over d\lambda}\,\gamma_k(g_\lambda)={2\over\gamma_k(g_\lambda)}
+\gamma_k(g_\lambda)S_k(g_\lambda)\,,
\equation(816)
$$
where 
$$
S_k(g_\lambda)=-\!\!\sum_{j\in I\setminus\{k-1,k\}}
\Bigl\vert\bigl[\nu_j(g_\lambda)-\nu_k(g_\lambda)\bigr]
\bigl[\nu_j(g_\lambda)-\nu_{k-1}(g_\lambda)\bigr]\Bigr\vert^{-1}.
\equation(817)
$$

Assume now that $\gamma(g)\prec\gamma(h)$.
Consider the set $\Lambda$ of all $x\ge 0$
such that $\gamma(g_\lambda)\prec\gamma(h_\lambda)$ for all $\lambda\le x$,
and assume that this set contains a maximum value $\lambda_0$.
Then there exists a non--empty set $K\subset\integer$ such that
$\gamma_k(g_{\lambda_0})=\gamma_k(h_{\lambda_0})<\infty$ for all $k\in K$.
At $\lambda=\lambda_0$ we have for all $k\in K$
$$
{d\over d\lambda} \bigl[\gamma_k(g_\lambda)-\gamma_k(h_\lambda)\bigr]
=\gamma_k(g_\lambda)\bigl[S_k(g_\lambda)-S_k(h_\lambda)\bigr].
\equation(818)
$$
But $S_k(u)\le S_k(v)$ whenever $\gamma(u)\prec\gamma(v)$,
and the inequality is strict whenever
$\gamma_k(u)=\gamma_k(v)$ and $\gamma(u)\ne\gamma(v)$.
Thus, it follows that $\gamma(g_\lambda)\prec\gamma(h_\lambda)$
for all $\lambda$ in some interval $(\lambda_0,\lambda_0+\varepsilon)$.
This contradicts the assumption that the set $\Lambda$ has a maximum.

In the case where $g$ has multiple zeros, the same argument can
be used first to prove that
$\gamma\bigl(g^{(n)}_\lambda\bigr)\prec\gamma\bigl(h^{(n)}_\lambda\bigr)\,$,
where $g^{(n)}$ and $h^{(n)}$ are defined such that
$$
\nu_i\bigl(g^{(n)}\bigr)=\nu_i(g)+i/n\,,\qquad
\nu_j\bigl(h^{(n)}\bigr)=\nu_j(h)+j/n\,,
\equation(819)
$$
for all zeros of $g$ and $h$, respectively.
Taking $n\to\infty$, the assertion now follows from\clm(810).
\qed

Consider now an arbitrary function $f$ in $\EE$
with indexed zeros $(I,\varphi)$. By Hadamard's factorization theorem
and a theorem of Lindel\"of [\rBo],
$f$ can be represented as a convergent product
$$
f(z)=ae^{bz}\prod_{k\in I}G\bigl(\mu(f)+i\nu_k(f),z\bigr)\,,
\qquad z\in\complex\,,
\equation(820)
$$
where
$$
G(\omega,z)=\cases{
e^{z/\omega}[1-z/\omega]\,,\quad &if $\omega\not=0$ ;\cr
z\,, &if $\omega=0$ .\cr}
\equation(820a)
$$
For every $n\in\natural$, denote by $I_n$
the intersection of $I$ with the set $\{-n,-n+1,\ldots ,n\}$,
and define the $n^{\rm th}$ partial product $f_n$ for $f$ by the equation
$$
f_n(z)=ae^{bz}\prod_{k\in I_n}G\bigl(\mu(f)+i\nu_k(f),z\bigr)\,,
\qquad z\in\complex\,.
\equation(821)
$$
In what follows,
the zeros of $f_n$ are always assumed to be indexed by the set $I_n\,$.

\claim Proposition(8100) Let $f_1, f_2, \ldots$
be the partial products for a function $f\in\EE$.
Then $f_n\to f$ in the topology of $\FF$.

\proof The task is to show that $\|f-f_n\|_\rho\to 0$ for every $\rho>0$.
We shall only consider the case $I=\{0,1,2,\ldots\}$ here;
the other cases are either similar or trivial (if $I$ is finite).
Let now $\rho>0$ be fixed,
and let $N$ be some positive integer such that
$\omega_k=\mu(f)+i\nu_k(f)$ is nonzero, for all $k\ge N$.
By using that $|e^z(1-z)|\le\exp(2|z|^2)$ for all $z\in\complex$,
we get for all $n\ge N$ the bound
$$
|f_n(z)/f_N(z)|=\prod_{k=N+1}^n|G(\omega_k,z)|
\le\exp\Bigl(2|z|^2\sum_{k=N+1}^n|\omega_k|^{-2}\Bigr)\,.
\equation(822)
$$
From the convergence of the product in\equ(820),
it follows that the sum over all $k$ of $\omega_k^{-2}$ converges.
But since all $\omega_k$'s have the same real part,
this implies that the sum $\sum_{k=0}^\infty|\omega_k|^{-2}$ converges as well.
Consequently, if $N$ is chosen sufficiently large,
we have $\|f_n/f_N\|_{4\rho}=1$ whenever $n\ge N$, and thus
$$
\|f_n\|_{2\rho}=\|f_N\cdot f_n/f_N\|_{2\rho}
\le\|f_N\|_{4\rho}\|f_n/f_N\|_{4\rho}=\|f_N\|_{4\rho}\,.
\equation(823)
$$
Under the same condition on $N$ and $n$, we get now the following bound.
$$
\eqalign{ \|f-f_n\|_\rho
&\le\sum_{k=n+1}^\infty\|f_k-f_{k-1}\|_\rho
=\sum_{k=n+1}^\infty
\left\Vert\bigl(G(\omega_k,.)-1\bigr)f_{k-1}\right\Vert_\rho\cr
&\le\sum_{k=n+1}^\infty\left\Vert G(\omega_k,.)-1\right\Vert_{2\rho}
\|f_{k-1}\|_{2\rho} \le\|f_N\|_{4\rho}
\sum_{k=n+1}^\infty\left\Vert G(\omega_k,.)-1\right\Vert_{2\rho}\,.\cr}
\equation(824)
$$
The last sum in\equ(824) converges and tends to zero as $n\to\infty$,
since we can extract a factor $\omega_k^{-2}$ from the term
$$
G(\omega_k,z)-1
=\omega_k^{-2}\bigl[z^2\bigl(K_{0,2}G(1,.)\bigr)(z/\omega_k)\bigr]\,,
\equation(825)
$$
and bound the $2\rho$--norm of the remaining factor $[\cdots]$
uniformly in $k$.
The latter follows from\clm(801), by using the fact that 
$|z/\omega_k|<|z|$ for sufficiently large $k$.
\qed

We are now ready to prove the first part of\clm(4).

\claim Proposition(811) Let $f$ be a function in $\EE$, and let $\lambda>0$.
Then the function $H_\lambda f$ lies in $\EE$,
has the same number (cardinality) of zeros as $f$,
and all its zeros are simple.

\proof Let $f$ be a function in $\EE$ with indexed zeros,
and let $f_1, f_2, \ldots$ be the partial products for $f$, defined
in equation\equ(821). By\clm(809), the functions $H_\lambda f_n$
all lie in $\EE_f\,$. And from\clm(803) and\clm(8100) it follows
that $H_\lambda f_n\to f$ in $\FF$, and that the zeros of $H_\lambda f$
are limits of zeros of the functions $H_\lambda f_n\,$.
Thus, the difference between any two zeros of $H_\lambda f$ lies in $i\real$,
which proves that $H_\lambda f\in\EE$.

The discussion after\clm(807) shows that the following holds
for every nonnegative integer $d$, and for every $\lambda\in\real\,$: 
If $f$ has $d$ zeros then $H_\lambda f$ has $d$ zeros.
The converse is also true since $H_\lambda^{-1}=H_{-\lambda}\,$,
by\clm(804). Thus, the same holds for $d=\infty$ as well.

Assume now for contradiction that $H_\lambda f$ has a zero $z_0$
of multiplicity $m>1$, for $\lambda=\lambda_0>0$.
Then by\clm(806), the function $H_\lambda f$ has $m$ zeros close to $z_0$
for $\lambda$ near $\lambda_0\,$,
and $2\lfloor m/2\rfloor$ of these zeros approach $z_0$ on curves
that are tangent to the line $z_0+\real$, as $\lambda\uparrow\lambda_0\,$.
This contradicts the fact that $H_\lambda\in\EE$ for all $\lambda>0$.
Thus, the zeros of $H_\lambda f$ are simple for all $\lambda>0$.
\qed

Let $f$ be a function in $\EE$ with indexed zeros $(I,\varphi)$.
By using\clm(806), together with the implicit function theorem,
we can find, for every $k\in I$, a positive real number $\lambda_k$ and
a continuous function $\lambda\mapsto\omega_k(f,\lambda)$ from
$[0,\lambda_k)$ to $\complex$, such that the following holds.

\item{$(a)$}
$\omega_k(f,0)=\varphi(k)\,$, for every $k\in I$.
\item{$(b)$}
$\bigl(H_\lambda f\bigr)\bigl(\omega_k(f,\lambda)\bigr)=0\,$,
for every $k\in I$ and for every $\lambda<\lambda_k\,$.
\item{$(c)$}
${\rm Im}\,\omega_{k-1}(f,\lambda)\le{\rm Im}\,\omega_k(f,\lambda)\,$,
whenever $\{k-1,k\}\subset I$ and $\lambda<\min(\lambda_{k-1},\lambda_k)$.

\noindent
The functions $\omega_k(f,.)$ will be referred to
as the local zero curves for $f$.

\claim Proposition(812) Let $f$ be a function in $\EE$ with indexed zeros,
let $\omega_k(f,.)$ be one of the local zero curves for $f$,
and let $f_1, f_2, \ldots$ be the partial products for $f$.
Then there exists $\varepsilon_k>0$
such that the zero curves $\omega_k(f_n,.)$ for $f_n$
converge to $\omega_k(f,.)$ uniformly on $[0,\varepsilon_k]$.

\proof 
By\clm(806) there exists $\varepsilon_k>0$ and $n_k\in\natural$, such that
for every positive $\varepsilon\le\varepsilon_k$ and for every $n\ge n_k$,
the zero $\omega_k(f,\varepsilon^2/2)$ lies within a distance
$\varepsilon^{3/2}$ of $\omega_k(f_n,\varepsilon^2/2)$,
and all other zeros of $H_{\varepsilon^2/2}f$ are at a distance
$3\varepsilon^{3/2}$ or more from $\omega_k(f,\varepsilon^2/2)$.
Given any positive $\varepsilon\le\varepsilon_k\,$, we can now use\clm(810),
with $\Lambda=[\varepsilon^2/2,\varepsilon_k]$ and $r\equiv 0$,
in order to find $N\ge n_k$
such that $|\omega_k(f,\lambda)-\omega_k(f_n,\lambda)|\le\varepsilon^{3/2}$
for every $\lambda\le\varepsilon_k$ and for every $n\ge N$.
\qed

The proof of\clm(4) is completed with the following proposition.

\claim Proposition(813) Let $f$ be a nonzero function in $\EE$,
and assume that $\bigl(H_{\lambda_0}f\bigr)(z_0)=0\,$,
for some $z_0\in\complex$ and some $\lambda_0>0$.
Then there exists a unique continuous function $\omega: [0,\infty)\to\complex$
such that $\omega(\lambda_0)=z_0\,$, and such that
$\bigl(H_\lambda f\bigr)\bigl(\omega(\lambda)\bigr)=0$ for all $\lambda\ge 0$.
 
\proof First, we note that the uniqueness follows from the fact that
all zeros of $H_\lambda f$ are simple, if $\lambda$ is positive.
Let now $V_0$ be the largest subinterval of $[0,\infty)$
containing $\lambda_0\,$,
on which there exists a function $\omega$ with the desired properties.
By\clm(805), the intersection of $V_0$ with $(0,\infty)$ is open.

Assume for contradiction that $0\not\in V_0\,$,
and denote by $\lambda_1$ the largest real number that is smaller
than any element of $V_0\,$. Since $\hat f$ is analytic,
we must have $|\omega(\lambda)|\to\infty$ as $\lambda\downarrow\lambda_1\,$.
Consider the function $g=H_{\lambda_1}f$.
Given that $g$ has zeros (by\clm(811)), we can index them and choose
a local zero curve $\omega_k(g,.)$ for $g$.
According to\clm(812), there exists a positive $\delta<\lambda_0-\lambda_1$
such that for any given $\varepsilon>0$,
the $k^{\rm th}$ zero curve for the partial product $g_n$ satisfies
the bound $|\omega_k(g,\lambda)-\omega_k(g_n,\lambda)|<\varepsilon$,
whenever $n$ is larger than some number $n_0(\varepsilon)$.
Define $\lambda_2=\lambda_1+\delta$.
Then, by\clm(810), there exists an $\varepsilon_1>0$ such that
$H_{\lambda_2}f$ has no zero other than $\omega(\lambda_2)$
within a distance $\varepsilon_1$ of $\omega(\lambda_2)$,
and such that for every positive $\varepsilon\le\varepsilon_1$,
the function $H_\delta g_n$ has a unique zero
$\zeta_n$ in the disk $|z-\omega(\lambda_2)|<\varepsilon$,
if $n$ is larger than some number $n_1(\varepsilon)\ge n_0(\varepsilon)$.
In addition, we have $\zeta_n=\omega_{j_n}(g_n,\delta)$ for some index $j_n$,
since all zero curves of $g_n$ are defined on $[0,\infty)$.
In order to show that $j_{n+1}=j_n$ for large $n$,
consider an open disk $D\in\complex$,
such that both $D$ and its closure contain the straight line segment $L$
connecting the two points $\omega(\lambda_2)$ and $\omega_k(g,\delta)$,
but none of the zeros of $H_\delta g$ that lie outside $L$.
By\clm(810), there exists $n_2\in\natural$ such that the function
$H_\delta g_n$ has the same number of zeros in $D$
as the function $H_\delta g$, for all $n\ge n_2\,$.
But since the (imaginary parts of the) zero curves of $g_n$ cannot cross,
it follows that $j_n=j$ for some fixed index $j$, whenever
$n$ is larger than $n_1(\varepsilon_1)$ and $n_2\,$.
Now, we can use that by\clm(809),
$\gamma(H_\lambda g_{n+1})\prec\gamma(H_\lambda g_n)$ for every $n$
and for every $\lambda\ge 0$. This shows that the sequence
$n\mapsto|\omega_j(g_n,\lambda)-\omega_k(g_n,\lambda)|$
is non--increasing as $n\to\infty$,
for every $\lambda$ in the interval $[0,\delta]$, i.e.,
that the functions $\omega_j(g_n,.)$ remain bounded on this interval.
On the other hand, if $\lambda=\delta$ then
$\omega_j(g_n,\lambda)$ converges to
$\omega(g,\lambda)=\omega(f,\lambda_1+\lambda)$, as $n\to\infty$;
and by\clm(810), the same holds for every positive $\lambda\le\delta$.
But this contradicts the assumption that $0\not\in V_0\,$,
which implied that $|\omega(\lambda)|\to\infty$
as $\lambda\downarrow\lambda_1$.

% prevent a page break right before \qed
\vskip0pt plus.1\vsize\penalty-75\vskip0pt plus -.1\vsize
A similar argument shows that the set $V_0$ cannot have a finite upper bound.
\qed

\smallskip\noindent
{\bf Proof of Theorem 1.6.\ }
Let $f$ and $g$ be functions in $\EE$ with indexed zeros
$(I,\varphi)$ and $(J,\psi)$, respectively.
As we just proved, there exists a family (indexed by $I$)
of continuous function $\lambda\mapsto\omega_k(f,\lambda)$
from the interval $[0,\infty)$ to $\complex$,
such that $\bigl(H_\lambda f\bigr)(z)=0$ if and only if
$z=\omega_k(H_\lambda f)$ for some $k\in I$.
An analogous statement holds for $g$,
and for the partial products $f_n$ and $g_n$ for $f$ and $g$, respectively.
In addition, we have $H_\lambda f_n\to H_\lambda f$
and $H_\lambda g_n\to H_\lambda g$ for every $\lambda\ge 0$,
as a consequence of\clm(8100) and\clm(803).

Let now $\lambda$ be a fixed positive real number,
and assume that $\gamma(f)\prec\gamma(g)$.
Then for every $n$ we have $\gamma(f_n)\prec\gamma(g_n)$,
and $\gamma(H_\lambda f_n)\prec\gamma(H_\lambda g_n)$
follows from\clm(809).
But by\clm(812) and\clm(810), the zeros $\omega_k(H_\lambda f_n)$
and $\omega_k(H_\lambda g_n)$ converge to $\omega_k(H_\lambda f)$
and $\omega_k(H_\lambda g)$, respectively, as $n$ tends to infinity.
Thus, we conclude that $\gamma(H_\lambda f)\prec\gamma(H_\lambda g)$.
\qed

%MODEL.2
\section The Dyson--Baker Hierarchical Model

Following parts of [\rBak], we show in this section
how the transformation $\RR$ derives from the full RG transformation $\TT$
for the Dyson--Baker hierarchical model in 3 dimensions,
and how the latter relates to a general RG transformation
for continuous--spin lattice models.
We do not consider any specific properties of the IR fixed point for $\TT$,
since the results of [\rKWi, Section 4] are easy to adapt
to the present situation, as the following discussion will show.

A statistical mechanics model of real--valued spins
on a set $\Lambda\subset\integer^d$ may be represented
by a parameterized family of measures $\mu$ on the space
of all possible spin configurations $\phi\in\real^\Lambda$,
such that the quantities of interest can be obtained
from the partition function $\int d\mu$
by differentiating it with respect to the parameters of the family.
The RG transformations considered below are maps $\mu\mapsto\tilde\mu$
with the property that $\int d\tilde\mu\equiv\int d\mu\,$,
which can be used for an iterative computation of partition functions.
Here, $\tilde\mu$ represents a model on a smaller set $\tilde\Lambda$,
whose spins are averages of the original spins.

One of the standard RG transformations is obtained by choosing
the following averaging operator $A$ on $\ell^2(\integer^d)$.
$$
\bigl(A\phi\bigr)(y)=N^{-d/2}\!\!\!
\sum_{x:\lfloor x/N\rfloor =y}\!\!\!\phi(x)\,,
\equation(301)
$$
where $N$ is some integer larger than $1$, and where $\lfloor z\rfloor$
is the point in $\integer^d$ obtained from $z\in\real^d$
by taking the integer part of its coordinates.
In order to simplify notation, assume now that $\Lambda=\integer^d$.
This will lead to integrals that are not well defined,
but the cure to this is well known; see e.g. [\rKWi].
In addition to $A$, we can also choose a Gaussian measure $\mu_C\,$,
with mean zero and (some given) covariance $C$,
which will become a trivial fixed point of the RG transformation.
Whenever possible, a given spin model will now be represented by
the function $F$, defined by the equation $d\mu(\phi)=d\mu_C(\phi)F(\phi)$.
Let $\alpha>0$ be such that the operator $\Gamma$, determined by the equation
$$
C=\Gamma+(\alpha^2N)^{-d}\bigl(C^{-1}AC\bigr)^*C\bigl(C^{-1}AC\bigr)\,,
\equation(302)
$$
has no negative eigenvalues.
Then the following defines a Wilson--Kadanoff type RG transformation $\TT$.
$$
\bigl(\TT(F)\bigr)(\phi)=\tilde F(\phi)
\equiv\int\! d\mu_\Gamma(\psi)
F\bigl((\alpha^2N)^{-d/2}(C^{-1}AC)^*\phi+\psi\bigr)\,,
\equation(303)
$$
where $\mu_\Gamma$ is the Gaussian measure on $\real^\Lambda$
with mean zero and covariance $\Gamma$.
It is easy to check that $\int d\mu_C\tilde F$ is equal to
$\int d\mu_C F$, as required.
Thus, the partition function for $F$ may be computed by iterating the
transformation $\RR$. This is particularly useful 
in cases where $C$ is unbounded but $\Gamma$ is bounded.

In hierarchical models, the covariance $C$ is chosen in such a way
that the integral $\int d\mu_\Gamma$ factorizes over blocks
$B(y)=\{x\in\Lambda : \lfloor x/N\rfloor=y\}$.
And if one wants to mimic the short--range Ising model,
the matrix elements $C(x,y)$ should decay roughly like the
two--point function of the Ising model; or like the kernel of the
inverse Laplacean, if $\eta=0$.
In order to define such a choice, consider the operators
$U$ and $A_0$ on $\ell^2(\Lambda)$, given by
$$
\eqalign{
\bigl(U\phi\bigr)(x_1,x_2,\ldots,x_d)&=\phi(x_2,x_3,\ldots,x_d,x_1)\,,\cr
\bigl(A_0\phi\bigr)(y_1,y_2,\ldots,y_d)
&=N^{-1/2}\!\!\!\!\!\sum_{x_d : \lfloor x_d/N\rfloor=y_d}\!\!\!\!
\phi(y_1,y_2,\ldots,y_{d-1},x_d)\,,\cr}
\equation(302a)
$$
and for $n=1,2,3\ldots$ let
$$
A_n=U^{-n+1}A_0U^{n-1},\qquad P_n=\bigl(A_1^*U^{-1}\bigr)^{n-1}
\bigl(1-A_1^*A_1\bigr)\bigl(UA_1\bigr)^{n-1}.
\equation(302b)
$$
Notice that $A_nA_n^*=1$ for all $n$, which implies that
the operators $A_n^*A_n$ are orthogonal projections.
By using this fact, it is easy to check that the $P_n$'s
are also a family of orthogonal projections,
and that $P_mP_n=\delta_{mn}P_n$ for all $m$ and $n$.
Let now $C=C_\infty\,$, where
$$
C_k=(4\sigma)^{-1}\sum_{n=1}^k(\alpha^2N)^{n-1}P_n\,,
\equation(304)
$$
for some fixed $\sigma>0$.
Since $A_dA_{d-1}\cdots A_1=A$ and $A^*P_nA=P_{n+d}\,$,
the covariance $C$ satisfies the equation
$$
C=\Gamma+(\alpha^2N)^dA^*CA\,,
\equation(304a)
$$
with $\Gamma=C_d\,$.
A straightforward calculation shows that $C^{-1}AC=(\alpha^2N)^dA$,
i.e., that $C=C_\infty$ and $\Gamma=C_d$ also satisfy\equ(302).
Using the decomposition of $A$ and $\Gamma$ given above,
we can now write $\TT=\TT_d\circ\TT_{d-1}\circ\ldots\circ\TT_1$, where
$$
\bigl(\TT_n(F)\bigr)(\phi)
=\int d\mu_{\Gamma_n}(\psi)F(\alpha N^{1/2}A_n^*\phi+\psi)\,,
\equation(305)
$$
with $\Gamma_n={1\over 4\sigma}(1-A_n^*A_n)$.
We note that if $\alpha=N^{1/d-1/2}$, then the spectrum of $C$,
as given by\equ(304), has a scaling property analogous to that
of the inverse Laplacean.

The property that distinguishes hierarchical models
from more realistic models is the fact that Gibbs factors of the form
$$
F(\phi)=\prod_{x\in\Lambda}\alpha
\sqrt{2\pi/\sigma}\,h\bigl(\phi(x)\bigr)
\equation(306)
$$
are mapped to Gibbs factors of a similar form by $\TT$.
In the case considered here, the integral in\equ(305)
factorizes into a product of $(N-1)$--dimensional integrals
associated with blocks of size $N$. In particular,
it is easy to check that for $N=2$
$$
\bigl(\TT_n F\bigr)(\phi)=\prod_{x\in\Lambda}\alpha
\sqrt{2\pi/\sigma}\,\bigl(\RR h\bigr)\bigl(\phi(x)\bigr)\,,
\equation(307)
$$
where $\RR$ is the transformation given in\equ(1), with $K=2\alpha$.

%REF.1
\bigskip\noindent
{\bf References}\hfill\break

\item{[\rDi]} F.J.~Dyson, {\it
Existence of a Phase Transition in a
One-Dimensional Ising Ferromagnet.}
Commun. Math. Phys. {\bf 12}, 91 (1969).

\item{[\rW]} K.~Wilson, {\it
Renormalization Group and Critical Phenomena. II.
Phase Space Cell Analysis and Critical Behavior.}
Phys. Rev. {\bf B4}, 3184 (1971).

\item{[\rBak]} G.A.~Baker, {\it
Ising Model with a Scaling Interaction.}
Phys. Rev. {\bf B5}, 2622 (1972).

\item{[\rMore]} For a more extensive list of references see [\rKWi].

\item{[\rG]} G.~Gallavotti, {\it
Some Aspects of the Renormalization Problems in Statistical Mechanics.}
Memorie dell' Accademia dei Lincei {\bf 15}, 23 (1978).

\item{[\rKWi]} H.~Koch, P.~Wittwer, {\it
A Non--Gaussian Renormalization Group Fixed
Point for Hierarchical Scalar Lattice Field Theories.}
Commun. Math. Phys. {\bf 106}, 495 (1986).

\item{[\rKWii]} H.~Koch, P.~Wittwer, {\it
On the Renormalization Group Transformation for Scalar Hierarchical Models.}
Commun. Math. Phys. {\bf 138}, 537 (1991).

\item{[\rKWiii]} H.~Koch, P.~Wittwer, {\it
Bounds on the Zeros of a Renormalization Group Fixed Point}.
In Preparation.

\item{[\rBo]} See e.g. Theorem (2.10.3) in:\hfill\break
R.P.~Boas, {\it
Entire Functions.}
Academic Press (1954).

\bye

