\section{introduction}
There is a growing interest\cite{list} on the subject of non-extensive
thermodynamics introduced ten years ago by
 Tsallis \cite{CONSTANTINO88}, by means of the entropy
 \begin{equation}
 S_{q}= \frac{1 - \sum_{i = 1}^{W}p_{i}^{q}}{q-1}.
 \label{entropy}
 \end{equation}
 Note that this entropy is characterized by the index $q$ whose
 departure from the
 conventional value $q = 1$ signals the thermodynamic effects of
 either long-range correlations in fractal dynamics
 \cite{heuristic1,heuristic2,heuristic3} or the non-local
 character of quantum mechanics\cite{luigi,rajagopal}.The growing interest for
 Tsallis' non-extensive entropy
 is testified by the exponentially
 growing list of publications on this hot issue\cite{list}.
 Of remarkable interest for the subject of fractal dynamics is the
 discovery recently made
 in Refs.\cite{heuristic1,heuristic2,heuristic3}
 that
 the entropic index $q$ also determines the specific analytical
 form, more general than the exponential form, adopted by two
 trajectories moving from infinitelly close initial conditions, to
 depart from one another.

  This important result is based
 on the generalization of the Kolmogorov-Sinai (KS) entropy \cite{K,S}
 and, consequently, of the important theorem of Pesin \cite{pesin}.
In Refs.~\cite{heuristic1,heuristic2,heuristic3} it was shown that the
adoption of
the Tsallis entropy yields a natural generalization of the
Kolmogorov-Sinai (KS)
entropy \cite{ruelle}. For the sake
of simplicity, henceforth
we shall be referring ourselves
to this form
of entropy
as Kolmogorov-Sinai-Tsallis (KST) entropy.
The adoption of the KST entropy yields a generalization of the
Pesin theorem \cite{heuristic1,heuristic2,heuristic3} and, with it,  the
important result:
\begin{equation}
\delta(t) \equiv \lim_{\Delta y(0) \rightarrow 0} \frac{\Delta
y(t)}{\Delta y(0)}
= [1 + \lambda_{q}t(1-q)]^{1/(1- q)}.
\label{delta(t)}
\end{equation}
The symbol $\Delta y(t) $ denotes the distance at time $t$ between
two trajectories departing from two close initial conditions.
Note that this
prescription for $q>1$ implies a divergence at finite time, and,
consequently, a
faster than exponential departure of two
initially very close trajectories from one to the other.
The case $q<1$ yields a power law
dependence on time, which can be interpreted as being slower than the
exponential departure. It is well known that the case $q=1$ results
in ordinary Brownian diffusion. In Ref.~\cite{condmatt} it was shown
that the case $q>1$ leads to superdiffusion. Consequently, the case
$q<1$, with a departure of the trajectories slower than exponential,
ought to result,
in general, in subdiffusion. It has to be pointed out that
this dynamical
argument fits the theoretical prediction of Ref.~\cite{tsallisbukman96}. The
authors of this interesting paper assign an entropic index $q$
to a non-linear Fokker-Planck equation, and prove
that the condition $q<1$ generates a subdiffusion process.
We are convinced that there is an intimate relation between
dynamic and probabilistic approach\cite{polish,bianucci}; as a
consequence of this conviction we are led to believe that also a dynamic
process characterized by $q<1$ must produce a form of
subdiffusion.

 As well known\cite{ruelle,pincus}, the concept of KS entropy
 implies the assumption of an invariant distribution. This condition, which
  must
 be extended to the KST entropy
 is ,  in turn, closely related to the stationary assumption recently
 adopted\cite{mannella} to establish a dynamic derivation of
 the diffusion processes. We want to prove that this assumption is
 the key physical ingredient necessary to make the dynamic treatment
 of
 anomalous diffusion compatible with a treatment based on
 non-extensive thermodynamics. This has also the effect of making
 subdiffusion as well as superdiffusion\cite{condmatt} compatible with
 the extended thermodynamics of Tsallis.

 On the other hand, in Section II we will show that
 the only form of subdiffusion compatible with the
 stationary condition is a localization process taking
 place with a finite time scale. For all these reasons
 we are led to make the conjecture that a generator of diffusion
 characterized by the entropic index $q < 1$ must result
 in a localization process, whose time scale is finite,
 and it is expected to tend to infinity with $q \rightarrow 1$.

 We shall not produce a general proof of this conjecture.
 We shall limit ourselves to
 a numerical investigation
 of a dynamical generator
 with $q<1$ and we shall show, using in fact a numerical approach,
 that it results in a localization process.
  The dynamic generator of diffusion studied in this paper
 is the logistic map.
 The rationale for this choice is that the logistic map is one of the
 most used prototypes to illustrate the transition to chaos
 \cite{chinesebook}. On top of that, there are already
 some papers\cite{heuristic1,heuristic2,heuristic3}
 devoted to studying the sensitivity to initial condition in the case
 of the logistic map, and claiming that the proper entropic
 index in this case is $q <1$. Thus, to prove our conjecture
 that the stationary dynamics with $q<1$ yield localization,
 we are only left with the problem of proving, with numerical
 calculations, that the
 diffusion process generated by the logistic map
 is in fact characterized by localization.

 The paper is organized as follows. In Section II we review the
 dynamic approach to diffusion, a picture based on very simple
 mathematical arguments. We show that the condition that anomalous
 diffusion is  slower than normal diffusion, and at the same time
 compatible with the stationary assumption, yields localization
 occurring within a finite time scale. Section III is devoted to
 discussing the time evolution of the KST entropy of the logistic map
 at the onset of chaos, and to the detection of the entropic index
 $q<1$.
 In Section IV we check the theoretical conjecture that
 $q <1$ yields localization with a finite time scale.
 We study also the localization breakdown caused
 by the control parameter approaching the condition of full chaos.
 Finally, Section V is devoted to concluding remarks.

\section{a dynamic approach to diffusion}
The simplest dynamical approach to diffusion is established by the
equation of motion\begin{equation}
\frac{dx}{dt} = \xi(t).
\label{mainequation}
\end{equation}
The variable $x$ is the diffusing variable and at a given time $t>0$
we have\begin{equation}
x(t) = \int_{0}^{t}\rm d$t$^{\prime}\xi(t^{\prime}) + x(0).
\label{x(t)}
\end{equation}
Following \cite{mannella} we make first of all the assumption that the
time evolution of $\xi(t)$ corresponds to a stationary statistical
condition, which, in the specific case of the logistic map studied in
this letter, implies the existence of an invariant
distribution, whose structure will be studied numerically in Section
IV. In this specific condition it is straightforward to
prove\cite{mannella} that\begin{equation}
\frac{d\langle x^{2}(t)\rangle}{dt}
 = 2 D(t),
\label{2D(t)}
\end{equation}
where $D(t)$ can be interpreted as being a time dependent diffusion
coefficient defined by\begin{equation}
D(t) \equiv
\langle \xi(t)^{2}\rangle_{eq}\int_{0}^{t} \Phi_{\xi}(t^{ \prime})
dt^{ \prime}
\label{defineD(t)}
\end{equation}
and the symbol $\Phi_{\xi}(t)$ denotes the autocorrelation
function of the fluctuation $\xi(t)$ defined by\begin{equation}
\Phi_{\xi}(t) \equiv
\frac{\langle \xi(0)\xi(t)\rangle_{eq}}
{\langle \xi^{2} \rangle_{eq}}.
\label{correlationfunction}
\end{equation}
Note that the stationary property implied by the one-time
correlation function of
Eq.(\ref{correlationfunction}) is generated by averaging over an equilibrium
distribution, or invariant measure, of initial condition.
Ordinary statistical mechanics is ensured by the condition that the
correlation time $\tau$ defined by\begin{equation}
\tau \equiv
\int_{0}^{t} \Phi_{\xi}(t^{ \prime})
dt^{ \prime}
\label{tau}
\end{equation}
exits and is finite, namely, that the important
condition\begin{equation}
0 < \tau < \infty
\label{importantcondition}
\end{equation}
is fulfilled.

The important condition of
Eq.(\ref{importantcondition}) can be violated in different ways.
A first example is given by a physical condition described by
the autocorrelation function $\Phi_{\xi}(t)$
 with the time
asymptotic property
\begin{equation}
\lim_{t\rightarrow\infty}\Phi_{\xi}(t )
= \frac{const}{t^{\beta }}
\label{asymptotics}
\end{equation}
and

\begin{equation}
0 < \beta < 1.
\label{anotherimportantcondition}
\end{equation}
This corresponds to the time dependent diffusion coefficient of
Eq.(\ref{defineD(t)}) becoming infinite in the asymptotic time limit.
Another way of violating this condition, on wich we focus our attention in
this paper, is when the time $\tau$ of
Eq. (\ref{tau}) vanishes. In this case, in the asymptotic time
limit  the time dependent diffusion
coefficient of Eq.({\ref{defineD(t)}) vanishes,
thereby implying that also the localization
process might take place asymptotically
in time.

The methods of nonequilibrium statistical mechanics
can be extended to deal with
processes violating the important condition of
Eq.(\ref{importantcondition}).
Examples can be found in the recent
literature, see for instance \cite{giorgio}, and in interesting review
papers\cite{geisel,klafter}.  Furthermore, as an effect of
the increasing popularity of the non-extensive entropy proposed
ten year ago by Tsallis\cite{CONSTANTINO88} it is becoming clear that
also the regime of anomalous diffusion can be given a thermodynamical
as well a dynamical significance\cite{condmatt}.

We aim at establishing a unified mechanical and
 thermodynamical picture ranging from normal to anomalous diffusion,
 namely,
 including both super and subdiffusion. For this
 ambitious task to be succesfully completed
 it is necessary to turn our attention to the
 stationary property. First of all, we have to point out that the
 analysis made in Ref.~\cite{condmatt} was limited to the case where
 the fluctuating variable $\xi(t)$ is dichotomous, and the map of
 Geisel and Thomae \cite{geiselthomae84} is used to establish the
 distribution of sojourn times in one or the other of the velocity
 states \cite{allegrini}. This waiting time distribution, called
 $\psi(t)$, has
 the inverse power law structure:
 \begin{equation}
 \lim_{t \rightarrow \infty}
 \psi(t)
= \frac{const}{t^{\mu}}
\label{psi(t)}
\end{equation}
with
 \begin{equation}
 2 < \mu <3 .
\label{restriction}
\end{equation}
The restriction to values of $\mu > 2$ is enforced by the need of
establishing a contact with the stationary theory earlier developed.
In fact, as recently pointed by Zaslavsky\cite{zaslavsky}, in the
specific case when the
waiting time distribution has an inverse power behavior, the asymptotic
time regime of the distribution of Poincar\'{e} recurence time is
dominated by the distribution of times of sojourn at the border between
the chaotic sea and the stability islands. On the other hand the Kac
theorem\cite{KAC} establishes that ergodicity implies that the
distribution of Poincar\'{e} recurrence times is characterized by a
finite first moment, thereby yielding $\mu >2$. According to
Ref.\cite{allegrini}, $\psi(t)$
is proportional to the second time
derivative of the correlation function $\Phi_{\xi}(t)$ of
Eq.( \ref{asymptotics}). Thus, $\mu = \beta +2$ and the condition $\mu >3$
implies $\beta>1$. The values of $\mu>3$, corresponding to a condition where
also the second moment of $\psi(t)$ is finite, means that the ensuing
diffusion process collapses into the attraction basin of the central
limit theorem, and consequently into the dominion of ordinary
statistical mechanics. In conclusion, the stationary condition and our
wish to study anomalous condition force us to consider the interval
defined by Eq.(\ref{anotherimportantcondition}).

The stationary condition involved in the dynamical
theory of diffusion behind Eqs.(\ref{2D(t)})
and (\ref{defineD(t)}) is shared by the mathematical foundation
of the KS entropy, resting on the assumption
that the dynamical system under study
has an invariant measure \cite{pincus}. This
makes it possible to settle the apparent paradox
associated with the use of the map of
Geisel and Thomae \cite{geiselthomae84}
for subdiffusion. It is known \cite{bettin} that
when the map of Geisel and Thomae is adopted
to establish the distribution of sojourn times
on space states rather than on velocity states,
it results in subdiffusion, with the following
time asymptotic expression of the space second moment:
\begin{equation}
 \lim_{t \rightarrow \infty} \langle x^{2}(t)\rangle
 = const \cdot t^{\mu -1},
 \label{bettinsubdiffusion}
 \end{equation}
where $\mu$ is the inverse power index of the function
$\psi(t)$ with
\begin{equation}
1 < \mu < 2.
\label{bettinrange}
\end{equation}
On the other hand, it has been recently
established \cite{condmatt} that according
to the prescriptions of Eq.(\ref{delta(t)})
the map of Geisel and Thomae must be associated
with the entropic index $q$ given by
\begin{equation}
q = 1 + \frac{1}{\mu}.
\label{annamarco}
\end{equation}
This prescription and Eq.(\ref{bettinsubdiffusion})
seem to contradict the simple rule that $q>1$ and
$q<1$ are associated to super and subdiffusion, respectively.
This apparent paradox is settled by noticing that Bettin
et al. \cite{bettin} have demonstrated that the
subdiffusion property of Eq.(\ref{bettinsubdiffusion})
is incompatible with the stationary approach
yielding Eqs.(\ref{2D(t)}) and (\ref{defineD(t)}).

After proving that the stationary assumption is an ingredient
necessary to
satisfactorily relate mechanics to thermodynamics along the lines
established by
Tsallis\cite{CONSTANTINO88,heuristic1,heuristic2,heuristic3},
 we now establish an
important property of subdiffusion: A dynamic
process fulfilling the stationary condition yields
localization within a finite time scale.  We have
seen that, in addition to superdiffusion, another way of breaking the
condition for normal diffusion is given by the property
\begin{equation}
\lim_{t \rightarrow \infty} D(t) = 0.
\label{vanishingdiffusion}
\end {equation}
Note that at $t=0$ the diffusion coefficient $D(t)$  vanishes, due to its
own definition. Thus we expect $D(t)$ to reach the maximum value $D_{max}$
at a finite time
$\tau$, corresponding to the duration of the transition from
 microscopic dynamics to the ``macroscopic'' diffusion regime.
Let us define ($ t > \tau$)
\begin{equation}
R(t) \equiv  \frac{D(t)}{D_{max}}
\label{R(t)}
\end{equation}
and the corresponding localization time:
\begin{equation}
T\equiv \int_{\tau}^{\infty} dt R(t).
\label{localizationtime}
\end{equation}
The case $T =\infty$ corresponds to the asymptotic property of
Eq.(\ref{bettinsubdiffusion}). This is a case of no interest for the
discussion of this paper, focusing on the stationary condition.
Note, however, that the stationary condition implies that use is
made of the correlation
function $\Phi_{\xi}(t)$. The asymptotic behavior of
this correlation function can be derived from Eq.(\ref{defineD(t)}) by
time differentiation. Let us also adopt the following
subdiffusion condition
\begin{equation}
 \lim_{t \rightarrow \infty} \langle x^{2}(t)\rangle
 = const \cdot t^{2 - \beta},
 \label{stationarysubdiffusion}
 \end{equation}
 with
\begin{equation}
1 < \beta < 2.
\label{integrable}
\end{equation}
This is the same subdiffusion regime as that of
Eq.(\ref{bettinsubdiffusion}). The change of notation with the
adoption of the new symbol $\beta$ is due to the fact that
we want now to establish a connection with
a stationary condition and with the correlation
function
$\Phi_{\xi}(t)$. From Eqs. (\ref{2D(t)}) and (\ref{defineD(t)}) we see that
the asymptotic behavior of the correlation function
$\Phi_{\xi}(t)$ is proportional to the second time derivative
of $x^{2}(t)$ of Eq. (\ref{stationarysubdiffusion}). Thus we see that
asymptotic limit of this correlation function
is characterized by negative
tail:
\begin{equation}
\lim_{t\rightarrow\infty}\Phi_{\xi}(t )
= - \frac{const}{t^{ \beta }}.
\label{negativetail}
\end{equation}
This
negative tail would yield a negative contribution to the diffusion
coefficient balancing, in part or totally, the positive contribution
stemming from the short times, thus leading to either normal
diffusion or localization at finite times.
If we want to set the condition of anomalous diffusion
with a vanishing diffusion coefficient,
we are left with a localization process
taking place with the
finite time scale $T$.

In conclusion, the main result of this Section is
that \emph{a dynamic process
with $q<1$ must involve localization with finite time scale}.
The demonstration proceeds as follows. First we assume
that the property found by the authors of Ref.
Ref.\cite{tsallisbukman96} that $q<1$ implies subdiffusion is of
general validity. The authors base their arguments on
Fokker-Planck-like
equations, and consequently on a probabilistic
view, which is historically dictated by
the need itself of understanding thermodynamics.
We are convinced that a unifying perspective
directly relating thermodynamics to dynamics is possible. On the
basis of this conviction we make the conjecture
that $q<1$ implies subdiffusion even when a totally dynamics
perspective is adopted. We note that the extension of KS entropy
rests on the existence of an invariant distribution and consequently
is compatible with a dynamic approach to diffusion
resting on the stationary assumption. Finally, the adoption of the
dynamic approach to diffusion\cite{mannella} shows that stationary
subdiffusion means localization occurring within a finite time scale.

\section{thermodynamics of fractal dynamics}
The authors of Ref.\cite{jin} expressed the Kolmogorov-Sinai-Tsallis
(KST) entropy in terms of an average over the invariant distribution
$p(x)$ by means of the following expression:
 \begin{equation}
  H_{q}(t) = [1 - k(q) \int dx p(x)^{q} \delta(t,x)^{1-q}]/(q-1).
  \label{fundamentalresult}
  \end{equation}
  Note that the departure point of the theory leading to this
  interesting result is given by a repartition of the phase space of
  the system under study into cells of small size. This is
  mirrored by the fact that $k(q) \equiv  (2/l)^{q-1}$, where
  $2/l$ denotes the size of these cells. The function $\delta(t,x)$
  is defined by Eq.(\ref{delta(t)}) with the dependence on the
  initial condition $x$ made explicit.

 It is worth noticing that the interesting result of
 ref.\cite{jin}, as expressed by
 Eq.(\ref{fundamentalresult}),
 fits the predictions of the more heuristic treatment of
 Refs.\cite{heuristic1,heuristic2,heuristic3}. Let us illustrate that
 this aspect in detail. Let us assume that the function
 $\delta(t)$ of Eq.(\ref{delta(t)}) is expressed by
 \begin{equation}
 \delta(t) = [1 + (1-q) \lambda_{Q} t]^{1/(1-Q)}.
 \label{typicalsensitivity}
 \end{equation}
Note that we are using the capital letter $Q$ rather than
$q$ on purpose. This is because we are assuming that
that only one power law exists, that this is characterized
by the power index $\beta= \frac{1}{(1-Q)}$ and that the resulting $Q$ is
the magic value that the analysis made in terms of the KST ought to
discover. Plugging Eq.(\ref{typicalsensitivity}) into the KST entropy
of Eq.(\ref{fundamentalresult}) we immediately assess that the
KST entropy become a linear function of time only when we assign to the mobile
entropic index $q$ the magic value $q = Q$.

 Note that
 the threshold of chaos is characterized by fractal dynamics. This
 means that system's dynamics significantly depart from the condition
 of full chaos, where\cite{zaslavsky} a nice connection between
 thermodynamics and mechanics can be established using the
 KS entropy. The non-extensive thermodynamics of Tsallis makes it
 possible to extend a thermodynamic
 perspective to the case of fractal dynamics .
We want to point out again that the onset of thermodynamics is
associated, in a full accordance with the traditional
 wisdom\cite{zaslavsky}, to the presence of a time region where
 the generalized form of Kolmogorov entropy grows linearely in time.

 Actually, as earlier pointed out by the authors of
 Refs.\cite{politi,mori}, in the real dynamic cases, the sensivity to
initial condition
 does not result in a
 simple expression as that of Eq. (\ref{typicalsensitivity}). Rather,
  a distribution of power indexes is found. Work in progress of our
  group aims at assessing which value of $Q$ will emerge from the
  search for
  linear increase of $H_{q}(t)$ of Eq.(\ref{fundamentalresult}), if a
  proper average over this distribution of power indices is made.
 Although a final conclusion
 has not yet been reached, we
 are convinced that the prediction  $Q <1$ of
 Refs.\cite{heuristic1,heuristic2,heuristic3} is correct.
  Thus we can restate our convition that the family of logisticlike maps
 studied in Section IV has to produce subdiffusion under the specific
 form of localization with a finite localization time.

\section{The logistic map and the process of dynamical localization
with a finite time scale}
We devote this Section to checking the conjecture that
$Q < 1$ yields a localization process with a finite localization time.
For the reasons explained in Section I, we use
the logistic map  as a diffusion generator. More precisely,
we follow \cite{heuristic3} and we
study
the logisticlike family of maps:
\begin{equation}
y_{n+1}= 1 - \mu | y_{n}|^{z},
\label{logisticmap}
\end{equation}
with
\begin{equation}
z> 1; 0 < \mu \le 2; n = 0,1,2,\ldots; y \in [-1,1].
\label{range}
\end{equation}

First of all, we evaluate numerically the invariant distribution,
a sample of which is illustrated in Fig. 1 for the case $z = 2, \mu =
\mu_{c}(2) = 1.4011551\ldots$. Then we define the fluctuation
\begin{equation}
\xi_{n} \equiv y_{n} - <y_{n}>_{eq},
\label{fluctuation}
\end{equation}
where $<y_{n}>_{eq}$ denotes an average over the invariant distribution.
This makes it possible for us to identify the fluctuation of
Eq.(\ref{fluctuation}) with the diffusion generator of
Eq.(\ref{mainequation}). In fact, it is evident that with
$n\rightarrow \infty$ the discrete variable of Eq.(\ref{fluctuation})
becomes virtually continuous as the variable $\xi$ of
Eq.(\ref{mainequation}).
Then we evaluate $x(t)$ using the discrete counterpart of
Eq.(\ref{x(t)}) with different values of z,
we determine $x(t)^{2}$ and we average it over the set of
initial conditions that in the case $z = 2$ corresponds to the invariant
distribution of
Fig.1.

Actually, we do our calculations for different values of $z$, and for
each $z$ we assign to $\mu$ the corresponding critical value
$\mu_{c}(z)$.
 According to the authors of
 Ref\cite{heuristic3} the best fitting to their numerical result is
 given by:
 \begin{equation}
 Q(z) = 1 - a_{0}/(z-1)^{a_{1}},
 \label{bestfitting}
 \end{equation}
 with $a_{0}= 0.75$ and $a_{1}= 0.60$. This means that with
 increasing $z$ the value of $Q$ tends to the value $Q = 1$,
 signalling the
 condition of standard statistical
 mechanics\cite{list,CONSTANTINO88,heuristic1,heuristic2,heuristic3}.
 This implies that for $z \rightarrow \infty$, the diffusion process
 must tend to become ordinary Brownian diffusion again. Let us see
 how the time evolution of the second moment $M_{2}(t) \equiv
 <x^{2}(t)>$ mirrors this earlier entropic analysis\cite{heuristic3}.
The time evolution of $M_{2}(t)$ is illustrated by Fig.2. We see that
the time evolution of $M_{2}(t)$
shows a scenery distinctly different from that of
ordinary statistical mechanics, corresponding
to a function $M_{2}(t)$  increasing linearly in time. At $z= 2$ the
function $M_{2}(t)$
immediately
settles in a state characterized by wild fluctuations but
 not exceeding a finite upper value. We also see that
 increasing $z$ from  $z =2$ to
 $z = 100$ has the effect of resulting
 in a visible localization onset time. In fact in the case $z=100$
 it takes the system about
 $25$ time steps to settle in another state characterized by wild
 fluctuations around a time indipendent mean value.
 Although the intensity of the fluctuations is larger than in the case
 $z = 2$, the transient nature of this transition to the localized
 state is clear.
 Exploring the time evolution of the logisticlike map of
 Eq.(\ref{logisticmap}) for larger values of $z$ so as to make
 more extended the time necessary to produce localization
 involves some technical problems\cite{capel}. However, on the basis
 of the results illustrated in Fig.2 we make the conjecture that the
 effect of increasing $z$ must be that of making increasingly larger
 the localization time so that the condition of ordinary statistical
 mechanics $Q(\infty) = 1$ is recovered when the localization time
 becomes infinite.

 We explored another path to full chaos, namely to the entropic
 condition characterized by $Q = 1$. This
 is obtained by keeping $z=2$ and moving $\mu(2)$ from
 the critical value $\mu = \mu_{c}(2) = 1.4011551\ldots$ to the value
 $\mu = 2$, in correspondence of which the logistic map is
 characterized by full chaos\cite{chinesebook}.
 First of all we adopt the same numerical method as that used in
 Refs.\cite{heuristic1,heuristic2,heuristic3} to establish the
 entropic index $Q$ corresponding to a given dynamical condition.
 We use this same method with a changing  value of the control parameter
 $\mu$: we select some
 values of $\mu$ in the interval $[1.4011551\ldots, 2]$ that do not
 correspond to the periodic windows.
 The behavior of $\delta(t)$ in these cases is illustrated by Fig.3 and the
 second moment of the diffusing distribution is illustrated by Fig.4.

 These numerical results are somewhat unexpected. In fact, we would
 have
 expected the localization time to become increasingly larger till to
 diverge at the condition of full chaos. On the contrary, we see that there
 exists a localization state regime. This is of infinite time
 duration at the critical threshold $\mu = \mu_{c}$. As the control
parameter $\mu$
 exceeds this critical value the localization state become unstable.
The closer and the closer the control parameter $\mu$ to the full chaos
condition, the shorter and the shorter the lifetime of this localized state.
 In the literature of localization processes it is possible to find a
 discussion of the effect produced by the environmental fluctuations
 on the phenomenon of localization exhibited by the kicked rotator
 (see, for instance, Ref.\cite{graham}). This is a kind of
 localization breakdown taking place as soon as the interaction with
 environment is switched on. The diffusion coefficient of the
 ensuing process of ordinary diffusion is proportional to
 the square of the intensity of the coupling between system and environment.
 Here, on the contrary, we see that the localization breakdown occurs
 in time, at times of increasingly larger duration as the control parameter
 approaches the threshold value $\mu_{c}$.

In conclusion, it seems that the effect revealed by the numerical
calculations of this paper is related to the conjecture of Ref.\cite{jin}
 that the entropic index Q is time dependent. As a result of a proper average,
made on the
the crowding index distribution of Ref.\cite{beck} rather than on the
coordinate $x$, we expect that Eq.(\ref{fundamentalresult}) results in a
linear in time entropy increase at
a "magic" entropic value Q. This value, in turn, is not the same at all
times. This is equivalent to saying that, in a sort of effective sense, the
form of Eq.(\ref{typicalsensitivity}) is replaced by

 \begin{equation}
 \delta(t) = [1 + (1-Q(t)) \lambda_{Q(t)} t]^{1/(1-Q(t))}.
 \label{agingeffect}
 \end{equation}

In accordance with the main result of this paper, let us assume that $ Q(t) <
1$ means localization occurring within a finite time scale, and, as shown by
Fig.2, almost instantaneously at $z=2$.  Then the results of Fig.4 can be
immediately accounted for if we assign to $Q(t)$ the following time behavior:
$Q(t)$ keeps a value smaller than the ordinary prescription $Q = 1$ for an
extended time interval corresponding to the time duration of the localized
states of Fig.4, then $Q(t)$ makes an abrupt jump to the traditional value $Q
= 1$.
Note that Fig.3, and other figures as well, concerning different
values of system's parameters and not shown here, exhibit an interesting
property worth of a comment. The overall time behavior of $\delta(t)$ looks
like exponential in accordance with the fact that the localization
breakdown implies that the standard form of sensitivity to the initial
conditions is recovered. However, a more careful analysis of the numerical
results shows that the exponential increase of $\delta(t)$ and the
multifractal oscillations of Refs. \cite{heuristic1,heuristic2,heuristic3}
 coexist in the same curve. The
multifractal oscillations last for a finite amount of time at a given
level, with a behavior very similar to that illustrated in Refs.
\cite{heuristic1,heuristic2,heuristic3}. At
a given time a sort of abrupt transition to another level, with the same
kind of multifractal oscillations, take place. The overall behavior looks
like that exponentiallike behavior of Fig.3. We think that the abrupt
transition from a localized state condition to the regime of Brownian
diffusion, shown in Fig.4, is  a consequence of the coexistence of these
two properties. Thus we make the following attractive conjecture: The
occurrence of a localization process with a finite time scale, as a result
of the stationary condition and of an entropic index $Q <1$, is closely
connected to the multifractal nature of the oscillations exhibited by
$\delta(t)$, as revealed by the analysis of Ref.
\cite{heuristic1,heuristic2,heuristic3}.
In these papers no significant attention was devoted to the
physical meaning of these
oscillations. We hope that the numerical results of this Section might
trigger further research work aiming at establishing the nature of this
attractive connection.

\section{Concluding remarks}
The so called KST entropy shares with the conventional KS entropy the
stationary condition, which is made evident by
Eq.(\ref{fundamentalresult}), resting in fact on an average over the
invariant distribution. The findings of Ref.\cite{tsallisbukman96}
refer to a probabilistic level of description, preceding the first
attempts at establishing a dynamic foundation of non-extensive
thermodynamics\cite{heuristic1,heuristic2,heuristic3,condmatt}.
We make the plausible assumption that the unification of dynamics and
thermodynamics\cite{polish,bianucci} can be extended. As ambitious as
this task is, this assumption is plausible, and the first results
obtained\cite{condmatt} are very encouraging. Consequently, it is
reasonable to make the conjecture that the results of
Ref.\cite{tsallisbukman96} have a general validity and apply also to
merely dynamic processes, described with no use of probabilistic
ingredients. This means that a dynamic process characterized by $Q <1$
is expected to yield subdiffusion.
On the other hand, according to the compelling results of Section III a
process of subdiffusion generated by stationary fluctuations must
lead to localization within a finite time scale. All these arguments
lead us to make the plausible conjecture that dynamics with $Q<1$ generate
localization, and that the localization time is finite.

This is not a rigorous demonstration. However, the numerical results
illustrated in Section IV fully support this conjecture and, in our
opinion, should serve the useful purpose of stimulating a search for a
more rigorous demonstration.


\section{Introduction}

The conventional picture of luminous early type cluster galaxies (E and
S0) is that their stars formed in a short period of time nearly a
Hubble time ago (Tinsley \& Gunn 1976).  Since the suggestion of
Sandage \& Visvinathan (1978) that the UV-optical color-magnitude
relation for cluster early type galaxies could be used to constrain
their star formation histories, studies of the cluster color-magnitude
relations have almost universally confirmed the conventional picture
(Bower, Lucey, \& Ellis 1992, BLE; Ellis et al.\ 1997); the luminous,
early type populations of rich clusters are remarkably homogeneous and
the bulk of their stellar populations appear to be old.

Dressler et al.\ (1997, D97), on the other hand, suggest that the
luminous S0 galaxies are less abundant
in distant clusters relative to the fractions of these types seen in present
day clusters.  They suggest that there must have been strong
evolution in the S0 population across the redshift range $z=0$--0.5 to
account for the large S0 populations seen in local clusters.  If S0
galaxies are being produced in clusters in recent times then they
should appear characteristically younger than the elliptical
population.  Indeed, Kuntschner \& Davies (1998), find such an effect
with the majority of the S0 galaxies in a sample drawn from the Fornax
cluster having luminosity-weighted ages several Gyrs younger than the
ellipticals.   Moreover, the post-starburst galaxies (PSG) found in
large numbers in the distant clusters (e.g.\ Couch \& Sharples 1987,
CS87; Dressler et al.\ 1999, D99) would provide an obvious intermediate
phase in the transformation of star forming disk galaxies into passive
S0's (Poggianti et al.\ 1999, P99).  However, this model then leads to
a potential conflict with those results confirming the conventional
model of old, coeval E and S0 galaxies.

In this Letter, we apply the spectroscopic technique of Kuntschner \&
Davies (1998) in clusters out to $z\sim 0.5$ to place new constraints
on the evolution of the S0 population in rich clusters by a direct
comparison of the relative ages of the cluster S0 and E populations at
three epochs spanning the range where strong evolution is claimed in
the proportion of S0 galaxies within rich clusters (D97).  In \S2, we
discuss the morphological and spectroscopic data as well as describing
the method we used to separate age and metallicity effects in the
spectra, then in \S3, we present the derived constraints on S0 evolution,
and finally in \S4, we suggest a resolution to the apparent conflict.

\section{Observations and Analysis}

Our analysis makes use of existing moderate resolution spectroscopy and
{\it HST} morphologies for galaxies in three clusters at $z=0.31$ and a
further 9 clusters between $z=0.37$--0.55. Existing spectroscopic and
morphological data for Coma are included as a local ($z\sim$0.0), rich
cluster comparison.  In particular we highlight
the usefulness of three clusters at $z=0.31$ from our sample in this
analysis -- for the evolution claimed by D97, almost half the S0
population of these clusters should have been produced in the 1.5\,Gyrs
since $z=0.5$.  This relatively high fraction coupled with the
reasonable numbers of S0 galaxies in these systems provides a strong
test of the D97 scenario.

\subsection{Morphological Samples}

For Coma, the morphologies are taken from Godwin, Metcalfe, \& Peach
(1983).  There are 63 elliptical and 121 S0 cluster members for which
we also have high quality spectra.

For the three clusters at $z=0.31$ (AC\,103, AC\,114, and AC\,118), we
have used the morphological classifications from Couch et al.\ (1998, C98)
and the details of the observations are contained therein.  The F702W
imaging consists of 3-orbit exposures of the cores of AC\,118 and
AC\,103, and a mosaic of four, 6 orbit exposures on AC\,114.  The
morphologies from Couch et al.\ (1998) are robust down to
$R_{702}\sim22.5$ yielding classifications for 25 spectroscopically-confirmed cluster ellipticals and 29 S0 cluster members.

Our more distant cluster sample comes from observations by the MORPHS
collaboration, consisting of deep {\it HST} imaging in the F702W, F814W, or
F555W
filters of 11 fields in 10 clusters in the redshift range $0.37 \leq z
\leq 0.56$.  Details of these observations can be found in Smail et
al.\ (1997).  The spectroscopic follow-up of this survey provides
membership for 15 elliptical and 10 S0 galaxies brighter than $R$=23.5
or $I$=23.0 (D99).  Smail et al.\ (1997) and D97 have
shown that the incidence of round and/or face-on S0 galaxies likely to be
misclassified as ellipticals in the distant clusters
is no worse than that seen in local samples (15\%).

\subsection{Spectroscopic Observations}

The Coma cluster data have been published by Caldwell et al.\ (1993)
and Caldwell \& Rose (1997, CR97), and the reduced galaxy spectra were kindly
provided by the authors.  The spectral resolution of their spectra is
$\sim4$\,\AA\ FWHM with S/N ranging from 30 to 100 -- substantially
above that achieved for individual galaxies in the distant samples.

Spectroscopic observations for the three clusters at $z=0.31$ are taken
from two sources, CS87 and C98.  The CS87 spectra were
collected at the 3.9-m Anglo-Australian Telescope (AAT) with the FOCAP
multi-fiber system and the IPCS detector, and have a spectral
resolution of $\sim 4$\,\AA\ FWHM.  The C98 spectra
were also collected at the AAT with LDSS-1, a multi-slit spectrograph, and
have a spectral resolution of $\sim 9$\,\AA\  FWHM.  The spectra in
both samples range in signal-to-noise ratio (S/N) from 3 to 10, with
the median being $\sim$8.

The spectroscopic observations of galaxies in the MORPHS sample is
presented in D99 and these spectra were kindly
provided by the authors.  The majority of these spectra come from the
COSMIC multi-object spectrograph on the 5.1-m Hale Telescope.  The
spectral resolution of these observations is 8\,\AA\ FWHM and the
median S/N is 25 for the spectra analysed here (the lowest having
S/N=5).

The majority of the spectra provide rest-frame wavelength coverage which includes
3900\,$\leq\lambda\leq$\,5300\,\AA, i.e.\ from the Ca{\sc ii} H \& K
lines through to the Mg triplet.  All of the galaxy spectra have been
smoothed to the 9 \AA\ resolution of the LDSS-1 data.

\subsection{Spectral Indices}

The main obstacle to an analysis of ages and metallicities of galaxies
is the degeneracy of those two effects on the galaxy spectral energy
distribution.  Focusing on individual spectral features has allowed
marked progress in breaking that degeneracy (c.f.\ Worthey 1994, W94;
Worthey \& Ottaviani 1997; Jones \& Worthey 1995), and for this paper
we have chosen seven of the commonly used Lick spectral indices.  The
H$\delta_{A}$ index acts as an age indicator (Worthey \& Ottaviani 1997),
while the Ca4227, Fe4383, Ca4455, Fe4531, C$_{2}$4668, and Fe5015 indices
trace metallicity effects (W94). This analysis relies upon high
quality spectra to robustly separate the effects of age and metallicity,
and given that we are looking for evidence of evolution in the characteristics
of the overall population of S0 galaxies {\it relative} to ellipticals,
we have chosen to co-add the S0 and elliptical spectra at each redshift
and perform a {\it differential} (E versus S0) analysis in each case
(z=0.0, 0.31, $\sim$0.50). Co-addition of the spectra was achieved using
the SCOMBINE routine in IRAF after having first transformed them onto
the rest wavelength scale. This provided final spectra with S/N$=40$ and
43 for the E and S0's in the $z=0.31$ clusters, and S/N$=97$ and 79 for
the same respective types in the $z\sim 0.5$ clusters. Seven indices were
measured in the co-added spectra, with those pertaining to the same
element being combined to form average values.

Figure 1 shall be used to present our results. The grid lines
represent lines of constant age (dotted lines) and constant metallicity
(solid lines) as taken from the single stellar population models of
W94; they are labelled in Panel c to indicate age differences in Gyr
and metallicity differences in dex units of [Fe/H]. Each redshift is
represented by a different shaped symbol with open and filled versions
being used to differentiate between the E and S0 measurements, respectively.
Regarding the separation of age and metallicity, we see that while the
age and metallicity lines are not orthogonal, they do indeed provide
discrimination between these two effects.  Note that any galaxy `ages'
discussed in this Letter will be luminosity-weighted mean ages, and are
not meant to indicate the ages of the oldest stars in the galaxies.
Also, because the three different samples of spectra have been
collected with different telescope--instrument--detector combinations,
and because we have not attempted to reconcile any offsets between
spectral index measurements of the different systems and the models, we
stress that there are zero point offsets relative to the age--metallicity
model grid associated with each redshift.  Accordingly, we cannot simply
compare `ages' between redshifts; only the relative ages are meaningful at
a given redshift.  This limitation does not bear on our analysis since we are
concerned only with a differential comparison between the E and S0 populations
at a single redshift. The following sections address the results derived
from Fig.~1, and their implications.

\section{Results and Discussion}

\subsection{Main Result}

Panels a--c of Fig.~1 show the H$\delta_{A}$ index plotted against the
mean of the three Fe indices, the mean of the two Ca indices, and the
C$_{2}4668$ index, respectively. The C$_{2}4668$ has been proposed as a
possible metallicity indicator and below we discuss its usefulness in the
context of our results. The error bars for H$\delta_{A}$ are approximately
the size of the points for the two higher redshift samples and one fifth
the size of the points for Coma, so the separation between points in
H$\delta_{A}$ is less than 1$\sigma$ for the E and S0 spectra at a given
redshift. Furthermore, in Panels~a \& b, no pair of points is separated by
more than 1$\sigma$ in the metallicity index. In contrast, we see for
the C$_{2}4668$ index good agreement between the E and S0 populations in
Coma, some consistency between the populations in the MORPHS sample (less
than 2$\sigma$ difference), but a 7$\sigma$ difference between the
apparent metallicities of the E and S0 population in the $z=0.31$
clusters from CS87/C98. We note that the CS87/C98 elliptical and S0
spectra were obtained in the same observing runs and reduced without
differentiating between the morphological types, so it is difficult to
imagine that a systematic error in the reduction could be responsible for
this disagreement. The near orthogonality of the age and metallicity
sequences with the index axes means that the  difference in the C$_{2}4668$
index for elliptical and S0 galaxies in the CS87/C98 sample translates
into an age difference of around 1--2 Gyrs.

The discrepent results from the C$_{2}4668$ index compared to the broad
agreement between our other two age/metallicity indicators requires an
explanation.  The C$_{2}4668$ index is extremely sensitive to changes
in the photospheric carbon abundance in stars (Tripicco \& Bell 1995)
and the most striking change in the carbon abundance of a star is in
the asymptotic giant branch (AGB) phase of evolution where stars reach
the thermally pulsing AGB and become Carbon stars (Iben 1993).  The
brightest {\em and} coolest\footnote{They must be cool
(K and M) stars because they would otherwise not have formed molecular
carbon.} of these stars have high initial masses and will be extinct by
about 1 Gyr after a burst of star formation suggesting that a recent
burst is indeed  the cause of the difference between the E and S0
C$_{2}4668$ indices. However, while the integrated number of Carbon
stars produced over the first 1 Gyr after a burst is enough to produce
the observed change in the C$_{2}4668$ index, the number of such stars
present {\em at any one time} is too small (Buzzoni 1989).  Moreover,
the similarity between the E and S0 galaxies at other redshifts (unless
the three CS87/C98 clusters are conspiratorily pathological), suggests that
we cannot appeal to global differences in carbon abundances due to
environment.  The lack of a simple explanation of the discrepent
C$_{2}4668$ indices leads us to concentrate on our other
age/metallicity indicators and we suggest that the C$_{2}4668$ index,
despite its empirical connection with metallicity demonstrated by
W94, requires more study before it can be reliably used to
investigate the age/metallicity of integrated stellar populations.

In summary, with the exception of the CS87/C98 C$_{2}4668$ index in the $z=0.31$ clusters, we find that
the luminosity weighted average age of elliptical and S0 galaxies in
rich clusters derived from this data are indistinguishable out to
$z\sim0.5$.

\subsection{Star Formation Histories}

The distribution of H$\delta_{A}$ values for E and S0 galaxies in the
different clusters can also be used to place limits on recent star
formation in the galaxies, and hence we have measured this spectral
index in the individual galaxy spectra as well.  Comparing
the E and S0 H$\delta_{A}$ distributions for the three different cluster
samples, we find using the KS test that they are indistinguishable
with the formal probabilities that they are drawn from different
populations being: P(KS)$_{\rm Coma}$=53\%, P(KS)$_{\rm CS87/C98}$=34\%,
and P(KS)$_{\rm MORPHS}$=0\%.

The D97 results suggest that the number of S0 galaxies
in the more distant MORPHS clusters would have to increase four-fold to
reproduce the E:S0 ratio in local clusters, while 50\% of the S0's in
the $z=0.31$ samples would need to have formed since $z\sim0.5$.  We
can directly evaluate the sensitivity of our analysis to this form of
evolution by selecting a random sample of 50\% of the $z=0.31$ S0's and
adding a contribution from a model spectrum of a young stellar
population from the library of Bruzual \& Charlot (1993) to each S0
spectrum and then recomputing the spectral indices.  We have done this
with two solar metallicity models, taken 0.5 and 1.0\,Gyrs after the end
of a  burst of star formation.  We then use the new indices to
recalculate the H$\delta_{A}$ distribution of the entire S0 sample, and compare
the new S0+model H$\delta_{A}$ distribution to the original E
H$\delta_{A}$ distribution.  This allows us to determine the minimum
contribution from a young population in the S0 sample which our
analysis of the $z=0.31$ clusters would detect (assuming the broad
framework of the evolution claimed by D97). The resulting constraints,
for the two different models, are that the H$\delta_{A}$ distribution
for the S0+model galaxies is statistically distinguishable from the
ellipticals at $\ge 96$\% (92\%) significance level when 10\%
(10\%) of the $V$-band light in half of the S0 sample is contributed by a
stars from a 0.5 (1.0)\,Gyr old post-starburst stellar
population.\footnote{This corresponds to a change in the average S0
H$\delta_{A}$ index of $\sim$0.5\AA, which, given that only 50\% of the
S0 population has the  young star light added, means the average index
change we are sensitive to in the post-starburst S0 galaxies is
$\sim$1.0\AA.}

Converting this constraint to the mass fraction of a galaxy converted to stars
in the burst, we conclude that the S0 galaxies at $z=0.31$ have not
experienced bursts of more than 6\% (11\%) by mass in the 0.5\,Gyr
(1\,Gyr) prior to the time of observation (Bruzual \& Charlot 1993).
These modest burst strengths are consistent with the
strength of the star bursts found by Barger et al.\ (1996) in these
same clusters.  The equivalent star formation rates are typical of
those seen in normal spiral galaxies.  We conclude that for our analysis
to detect a change in the characteristic age of S0 galaxies relative to
ellipticals at the same redshift, at least half of the S0's in these
clusters would have to have undergone particularly intense bursts of star
formation ($>11$\% of their current stellar mass) in the 1\,Gyr prior to
observation, i.e.\ between $z=0.3$ and $z\sim 0.5$.  {\it This limit
suggests that the bulk of the cluster S0 population is not being built up
through a process of intense starbursts where spirals are converting
their entire gas reservoir into stars in a single event.}

\section{Conclusions}

By combining high resolution {\it HST} imagery and existing
spectroscopy, we have been able to investigate the evolution of early
type galaxies in moderate redshift clusters.  We have used the stellar
population models of Worthey (1994) to disentangle the effects
of age and metallicity in our spectra. We have two main conclusions: Firstly,
we find no evidence for an age difference between E and S0 galaxies in rich
clusters out to $z\sim$0.5. Secondly, the precision with which our
spectral index method rules out an age difference would imply that any
recent star formation in the cluster S0's at $z=0.31$ could have involved
no more than 50\% of these galaxies undergoing bursts in which no more
than $\sim 11$\% of their mass was converted into stars.

Poggianti et al.\ (1999) have analysed the spectroscopic and
morphological catalogs gathered for the $z=0.37$--0.57 clusters by
D99.  They find evidence for a wide-spread
suppression of star formation in cluster members, leading to large
populations of PSG and passive spirals which are much more prevelant in
the cluster environment. The progenitors of these galaxies are
identified with dusty, star
burst galaxies infalling from the surrounding field.  They suggest
that the resulting passive, disk population are the progenitors of the S0
galaxies.  At lower redshifts, Caldwell et al.\ (1996, C96) have shown
that in the Coma cluster, Balmer line strengths of PSG early type
galaxies cannot be explained by a previously constant star formation
rate, typical of spiral galaxies, which has been recently truncated,
but that they require truncation to occur after a brief burst of
$\sim$10\% by mass, and well within the above limits placed on
starbursts at higher redshifts.  In CR97 they conclude that these
final, modest bursts of star formation in Coma, as well as other nearby
rich clusters, are indeed associated with the cluster environment.
Also, in both C96 and CR97, they show that most of the starburst and
post-starburst galaxies in their clusters are, morphologically, early
type disk galaxies, either early type spirals or S0's.

In light of the evidence from D97 that the S0 population in rich
clusters needs to strongly increase between $z\sim 0.5$ and today, the
results of P99, C96 and CR97, together with the results of this Letter,
suggest a mostly unspectacular origin for the missing S0 galaxies.  We
conclude from the assembled evidence that the missing S0 galaxies are
the remnants of old, early-
type spiral galaxies which, as a result of
their interaction with the rich cluster environment, have had their
star formation truncated after possibly undergoing a final, modest
burst of star formation.  This conclusion suggests that with more
sensitive age indicators and spectra of higher signal-to-noise and
higher spectral resolution, we would be able to detect the effects of
those bursts and test this hypothesis.  Such observations are possible
from current 4- and 8-meter class telescopes.

\acknowledgements

The authors would like to thank the MORPHS collaboration for the use of their
spectra in advance of publication, and Jim Rose and Nelson Caldwell for
provding the Coma spectra.  This work was supported by the Australian Research
Council.  IRS acknowledges additional support from PPARC and a travel grant
from the Royal Society.

\clearpage


\section{Introduction}

Applications of dimensional analysis to  mechanical problems are
well known in physics. \par\indent
From a historical point of view, the first seeds of that method can be
found in Newton's ``\emph{Principia}'' ${}^{1}$ and in Fourier's
``\emph{Th\'eorie Analytique de la Chaleur}'' ${}^{2}$.
Successively, important contributions to the modern treatment of
dimensional theory have been made by R.C.Tolman who defined the
principle of dimensional homogeneity and stated the $\Pi$
theorem ${}^{3}$. \par\indent
It is worth noticing that dimensional analysis was employed by
Albert Einstein too in the calculus of the eigenfrequencies in
solid bodies  ${}^{4}$. He set forth the pendulum example and, on
obtaining $T= C \sqrt{\dfrac {l}{g}} $ as usual, said: ``\emph{One can,
as is known, get a little more out from dimensional considerations,
but not with complete rigor. The dimensionless numerical factor
(like $C$ here), whose magnitude is only given by a more or less
detailed mathematical theory, are usually of the order of unity. We
cannot require this rigorously, for why shouldn't a numerical factor
like $(12\pi )^{3}$ appear in a mathemathical-physical deduction? But
without doubt such cases are rarities \ldots }''  ${}^{5}$.\par\indent
For a more complete view of the subject, we refer the reader to the
book of Bridgman~${}^{6}$.\par\indent
Our aim here is deducing the form of the Stephan-Boltzman law and of
the Planck's formula in a  space-time with one time and $n$ space
coordinates, by employing the methods of dimensional analysis only.
The minimum knowledge of physics we need getting started is given by
the definitions of length $L$, mass $M$, time $T$ and absolute
temperature $\theta$; we also accept the validity of Gauss's law  in
all dimensions and know the definitions of force and energy. Finally,
in the particular case $n=3$ we shall recover the known results
on the basis of little physical information.\par\indent

\section{The universal constants and the $\Pi$ theorem}

Let us recall here some elements of physical dimensions theory.
For instance, if we aim at a physical theory of gravity, the three
fundamental units $L,M$ and $T$ are sufficient, but they must be
associated as usual to some universal constants. In the gravitational
case two universal constants are needed: an invariant velocity $c$ with
dimensions $[c] = [LT^{-\, 1}]$ and a gravitation constant $G_{n}$,
where the subscript $n$ stands for the number of space coordinates,
with physical dimensions $[G_{n}] = [L^{n}M^{-\, 1}T^{-\, 2}]$. The
latter dimensional equality follows immediately from Gauss's law for
the gravitational field $\boldsymbol{g}_{n}$ produced by a
$n$-dimensional distribution of mass $\rho_{n}$:
\begin{equation}
\boldsymbol{\nabla \cdot g}_{n} = -\, \dfrac{2\pi^{n/2}}
{\Gamma (n/2)}\, G_{n}\rho_{n}
\end{equation}
and from the relation
\begin{equation}
\boldsymbol{F}_{n} =  m \boldsymbol{g}_{n}
\end{equation}
which gives the force $\boldsymbol{F}_{n}$ (with dimensions
$[LMT^{-\, 2}]$) acting on a test mass $m$ located in the gravitational
field $\boldsymbol{g}_{n}$.\par\indent
If the source is a point particle of mass $M$ we have, in a
$n$-dimensional space:
\begin{equation} |\boldsymbol{g}_{n}| =
\begin{cases}
\dfrac{G_{n}M}{r^{n-1}} \hspace{2in} &\text{for} \quad n>1 \\ {}\\
\dfrac{G_{1}M}{2} \hspace{2in} &\text{for} \quad n=1
\end{cases}\end{equation}
In dealing with electromagnetism, we need again $c$ and another
fundamental constant.We do not introduce any additional physical
quantities besides $L,M$ and $T$; so,
Gauss's law for the electrostatic field $\boldsymbol{E}_{n}$
produced by a $n$-dimensional distribution of charge $\rho_{n}$ will
be written, in Gaussian units, as:
\begin{equation}
\boldsymbol{\nabla \cdot E}_{n} = \dfrac{2\pi^{n/2}}{\Gamma (n/2)}\,
\rho_{n}
\end{equation}
Acting as in the gravitational case, the electrostatic field
$\boldsymbol{E}_{n}$ of a point charge $Q_{n}$ (whose dimensions
turn out to be $[L^{n/2}M^{1/2}T^{-\, 1}]$) is :
\begin{equation} |\boldsymbol{E}_{n}| =
\begin{cases}
\dfrac{Q_{n}}{r^{n-1}} \hspace{2in} &\text{for} \quad n>1 \\ {}\\
\dfrac{Q_{1}}{2} \hspace{2in} &\text{for} \quad n=1
\end{cases}\end{equation}
Before going on, we notice that in the three dimensional space
Coulomb's law for two point charges $Z_{1}e$ and $Z_{2}e$
\begin{equation}
|\boldsymbol{F}| = \dfrac {Z_{1}Z_{2}e^{2}}{r^{2}}
\end{equation}
can be written, by introducing into it the fine structure
constant\linebreak $\alpha = e^{2}/(\hbar c)$, as
\begin{equation}
|\boldsymbol{F}| = \dfrac {Z_{1}Z_{2}\alpha\hbar c}{r^{2}}
\end{equation}
So the other fundamental constant we need can be either the electron
charge $e$ or the Planck's constant $\hbar$. Therefore when dealing
within both gravity and electromagnetism, we can associate $(L,M,T)$
with $(c,G,e)$ or with $(c,G,\hbar )$: this agrees with the choices
made, respectively by Stoney in $1874$ ${}^{7}$ and by Planck in
$1899$ ${}^{8}$. In our approach the appearance of the Planck's
constant $\hbar$ is clearly related to the quantization of the electric
charge, however.\par\indent
Let us recall that, according to Planck, ${}^{8}$ the system of
``natural units'' must include $\sqrt{(\hbar c)/G}$ as unit of mass,
$\sqrt{(\hbar G)/c^{3}}$ as unit of length and $\sqrt{(\hbar G)/c^{5}}$
as unit of time; while according to Stoney ${}^{7}$ the ``natural
units'' for mass, length and time are $e/\sqrt{G}$,
$(\sqrt{G} e)/c^{2}$ and $(\sqrt{G}e)/c^{3}$ respectively.
One sees immediately that Stoney's units are
given by Planck's  times the square root of the fine structure
constant. By the way, one may notice that Stoney's conception of the
electron is somewhat different from that resulting from the subsequent
discovery by Thomson.\par\indent
Going back to our task, it is suitable, in $n$ spatial dimensions,
to pass from $G$ to $G_{n}$ and from $\alpha$ to
$\alpha_{n} = e_{n}^{2}/(\hbar_{n} c)$; so, being $\alpha_{n}$
dimensionless by definition, the replacement of $e$ with $e_{n}$ and
of $\hbar$ with $\hbar_{n}$ makes their respective dimensions to be
$[e_{n}] = [L^{n/2}M^{1/2}T^{-\, 1}]$ and
$[\hbar_{n}] = [L^{n-\, 1}MT^{-\, 1}]$. Of course,
the quantities labelled with $n$ are defined only formally and not
operationally, but this procedure will provide interesting results.\par\indent
Let us fix our attention at this point on ordinary thermodynamics.
In this case we need a fourth  quantity, namely, the absolute
temperature $\theta$ which then must be associated with another
universal constant. In a quite natural way, we choose the Boltzmann
constant $k$ whose dimensions are $[L^{2}MT^{-\, 2}\theta^{-\, 1}]$.
\par\indent Consequently in the dimensional equations which refer to the
thermodynamics of the electromagnetic radiation  factors of the type
$c^{\alpha}\hbar_n^{\beta}k^{\gamma}$ must appear.\par\indent
Let us finally recall the so called $\Pi$ theorem ${}^{6}$ which we
state as follows: ``Denote by $P_{1},P_{2},\ldots P_{r}$ the
magnitudes of quantities which may be physical magnitudes or
experimental constants. Suppose that a functional relation
$f(P_{1},P_{2},\ldots P_{r}) = 0$ holds,
and furthermore that this is a complete equation
and is the only relation between the magnitudes. Then, if there are
$s$ fundamental magnitudes the relation can be expressed in the form
$F(\Pi_{1},\Pi_{2},\ldots \Pi_{r-\, s}) = 0$ where the $\Pi$'s are
the $r-\, s$ independent products of the arguments $P_{1},P_{2},
\ldots P_{r}$ which are dimensionless in the fundamental magnitudes''.

\section{The Stefan-Boltzmann law}

Now we are ready to establish the Stefan-Boltzmann law in a
$n$-dimensional space. Denoting by $u_{n}(\theta )$ the energy
density of the electromagnetic radiation enclosed in a cavity and
in equilibrium at the temperature $\theta$, the
functional relation quoted in the $\Pi$ theorem is $f(u_{n},\theta ,c,
\hbar_{n},k) = 0$. In this case $r-\, s=5-\, 4=1$, so there is only
one product $\Pi$ which can be written in the form:
\begin{equation}
\Pi_{1} = u_{n}^{-1} c^{\alpha_1}
\hbar_{n}^{\beta_1}k^{\gamma_1} \theta^{\delta_1}
\end{equation}
Remembering that $[u_{n}] = [L^{2-n}MT^{-2}]$ and that consequently
we shall have to discard the values $n=1$ and $n=2$ by physical
reasons, we obtain the algebraic system:
\begin{equation}
\left\{ \begin{split} &\beta_1 + \gamma_1 - 1 = 0  \\
&\alpha_1 + (n-1)\beta_1 + 2\gamma_1 + n -2 = 0 \\
&\alpha_1 + \beta_1 + 2\gamma_1 -2 = 0  \\ &\gamma_1 -
\delta_1 =0 \end{split} \right.
\end{equation}
For $n \geq 3$ the solutions are:
\begin{equation}
\left\{ \begin{split} &\alpha_1 = -\, \dfrac {n}{n-2}
\\ &\beta_1 = -\, \dfrac {n}{n-2}  \\ &\gamma_1 =\dfrac
{2(n-1)}{n-2}  \\ &\delta_1 = \dfrac {2(n-1)}{n-2}
\end{split} \right.
\end{equation}
and therefore
\begin{equation}
u_n \propto \, \dfrac {(k\theta)^{2(n-1)/(n-2)}}
{(\hbar_n c)^{n/(n-2)}}
\end{equation}
When $n=3$ we recover the Stephan-Boltzmann law:
\begin{equation}
u \propto \, \dfrac{(k\theta)^4}{(\hbar c)^3}
\end{equation}
Obviously the numerical factor $\pi^{2}/15$ which appears in that
law cannot be determined by dimensional analysis.
\section{The Planck's formula}

Let us now obtain the Planck's formula in a $n$-dimensional space.\par\indent
Calling $u_{n,\omega}(\omega ,\theta )$ the spectral energy density
of the black-body radiation in the frequency interval
between $\omega$ and $\omega + d\omega$, the functional relation is
now $f(u_{n,\omega },\omega ,\theta ,c,\hbar,k)=0$. In this case
$r-s = 6-4 = 2$, so there are two products $\Pi$'s which will be
written as
\begin{align}
\Pi_2 &= u_{n,\omega }^{-1}c^{\alpha_2}\hbar^{\beta_2}k^{\gamma_2}
\omega^{\delta_2} \\
\Pi_3 &=\theta^{-1}c^{\alpha_3}\hbar^{\beta_3}k^{\gamma_3}
\omega^{\delta_3}
\end{align}
and satisfy a relation of the form $F(\Pi_2,\Pi_3) = 0$. \par\indent
The two algebraic systems one obtains are:
\begin{equation} \hspace{-0.5in}
\left\{ \begin{split} &\beta_2 + \gamma_2 - 1 =0  \\
&\alpha_2 + (n-1)\beta_2 + 2\gamma_2 + n - 2 = 0  \\
&\alpha_2 + \beta_2 + 2\gamma_2 + \delta_2 - 1 = 0  \\
&\gamma_2 = 0  \end{split} \right.
\hspace{0.3in}
\left\{ \begin{split} &\beta_3 + \gamma_3 = 0  \\
&\alpha_3 + (n-1)\beta_3 + 2\gamma_3 = 0 \\ &\alpha_3 +
\beta_3 + 2\gamma_3 + \delta_3 = 0  \\ &\gamma_3 + 1 = 0
\end{split} \right.
\end{equation}
with the respective solutions:
\begin{equation} \hspace{-0.5in}
\left\{ \begin{split} &\alpha_2 = -\, (2n - 3)  \\
&\beta_2 = 1  \\ &\gamma_2 = 0  \\ &\delta_2 = 2n - 3
\end{split} \right.
\hspace{1.5in}
\left\{ \begin{split} &\alpha_3 = -\, (n - 3)  \\
&\beta_3 = 1  \\ &\gamma_3 = -\, 1  \\ &\delta_3 = n - 2
\end{split} \right.
\end{equation}
Therefore from the relation $F(\Pi_2,\Pi_3) = 0$ it follows that
\begin{equation}
u_{n,\omega} = \dfrac {\hbar_n \omega^{2n-3}}{c^{2n-3}}\, \varphi
\left( \dfrac {\hbar_n \omega^{n-2}}{c^{n-3}k\theta}\right)
\end{equation}
On the other hand it results, by definition, that $\int_{0}^{\infty}
u_{n,\omega}(\omega ,\theta )\, d\omega = u_n(\theta )$; so, by using
the expressions already found for $u_{n,\omega}(\omega ,\theta )$ and
for $u_n(\theta )$, we obtain the following equation to be
satisfied by the unknown function $\varphi$:
\begin{equation}
\int_0^\infty \dfrac {\hbar_n \omega^{2n-3}}{c^{2n-3}}\, \varphi
\left( \dfrac {\hbar_n \omega^{n-2}}{c^{n-3}k\theta}\right)
d\omega \propto \dfrac {(k\theta)^{2(n-1)/(n-2)}}
{(\hbar_n c)^{n/(n-2)}}
\end{equation}
Passing to the new variable  $\varepsilon_n = (\hbar_n
\omega^{n-2})/c^{n-3}$, whose meaning as a photon energy is apparent,
the last equation becomes:
\begin{equation}
\dfrac {1}{(n-\, 2)}\,\int_0^\infty
\varepsilon_n^{2(n-1)/(n-2)\,-\, 1} \varphi \left(\dfrac
{\varepsilon_n}{k\theta}\right)\, d\varepsilon_n  \propto
(k\theta)^{2(n-1)/(n-2)}
\end{equation}
As to the function $\varphi$, some considerations are in order.\par\indent
First of all $\varphi$ cannot contain monomial terms in
$(\hbar_n \omega^{n-2})/(c^{n-3}k\theta)$, otherwise it would
modify the correct depedendence in Eq.(17) of $u_{n,\omega}$ on the
term $\omega^{2n-\, 3}$ which comes by multiplying the
factor $\omega^{n-\, 2}$ contained in the photon energy
$\varepsilon_n$ by the factor $\omega^{n-\, 1}$ due to the evaluation
of the number of photon states in the interval $d\omega$.\par\indent
Second, when $n=3$ we know the behavior of
$u_{n,\omega}$ in the limits $\omega \to 0$ (Rayleigh-Jeans's law)
and $\omega \to \infty$ (Wien's law), so, when $n=3$, one must have
\begin{equation}
\varphi \left( \dfrac{\hbar \omega }{k\theta } \right) \sim
\begin{cases}
\left( \dfrac {\hbar\omega }{k\theta }\right)^{-\, 1} \hspace{1.25in}
{\text for} \quad \omega \to 0  \\ {} \\
\exp\left(-\,\dfrac {\hbar \omega }{k\theta }\right) \hspace{1in}
{\text for} \quad \omega \to \infty
\end{cases}\end{equation}
From a strictly mathematical point of view, the function $\varphi$
which satisfies Eq.(19) is not unique, but, taking into account our
requirements, the choice is very restricted and, discarding too
involved expressions, one is left only with the function
\begin{equation}
\varphi \left(\dfrac {\varepsilon_n}{k\theta }\right) \,\propto
\,\dfrac {1}{\exp{\left(\dfrac {\varepsilon_n}{k\theta }\right)}- 1}
\end{equation}
Once again, we notice that for $n=3$ Eq.(17) becomes
\begin{equation}
u_{n,\omega} \propto \dfrac {\hbar_n \omega^{2n-3}}{c^{2n-3}}\,
\dfrac {1}{\exp{\left(\dfrac {\varepsilon_n}{k\theta }\right)}- 1}
\end{equation}
so that one recovers, apart from a numerical factor $1/\pi^{2}$, the
Planck's formula and hence the Bose-Einstein statistics.\par\indent
Had we ignored the behavior of the function $\varphi$ as $\omega
\to 0$, the choice made in Eq.(21) might have been changed to $\varphi
(\varepsilon_{n}/(k\theta )) \propto
\exp (-\,\varepsilon_{n}/(k\theta ))$, which is proper for the
Boltzmann statistics, or even to $\varphi (\varepsilon_{n}/(k\theta ))
\propto (\exp (\varepsilon_{n}/(k\theta )) + 1)^{-\, 1}$, which is
proper for the Fermi-Dirac statistics.

\section{Conclusions}

We would like to focus attention on some interesting properties of
dimensional analysis we made use of in this
paper.\par\indent
First of all, in treating by the dimensional method a quite general
problem of physics, it appeared  of capital importance to associate
the four fundamental quantities ($L,M,T,\theta $) with as
many universal constants, namely ($c,G_{n},\hbar_{n},k$). In less
general situations, one needs only some of these costants; for example,
when restricting to the three-dimensional space, one can select $c$
and $G$ for treating gravity, $c$ and $\hbar$ (or $e$) for treating
electromagnetism or, as in our case, $c,\hbar$ and $k$ for treating
the electromagnetic radiation thermodynamics.\par\indent
On the other hand, such a prescription was alredy contained in the
works of Stoney ${}^{7}$ and of Planck ${}^{8}$. In Planck's words
natural units are such because\linebreak ``\emph{\ldots they are
independent of particular bodies or substances, necessarily retain
their meanings in all times and in every culture even if
extra-terrestrial or extra-human.}'' ${}^{9}$ \par\indent
Secondly, let us comment on the power of the physical dimensions
theory to test the soundness of a physical law.
To make an example let us consider the Rayleigh-Jeans's formula for
black-body radiation. We can obtain it if, contrary to what we just
stated, we  accept the  widespread belief that electromagnetism
can be treated by employing as universal constant only
$c$, disregarding $\hbar$ (or $e$). Acting as before, we would write,
in three space dimensions, the new functional relation among the
quantities of interest as
$f(u_{,\omega},\omega ,\theta ,c,k) = 0$. Now we have $r-s = 5 - 4 =1$
so there is only the following product
\begin{equation}
\Pi = u_{,\omega}^{-\, 1} \omega^{\alpha} \theta^{\beta} c^{\gamma}
k^{\delta}
\end{equation}
The algebraic system is
\begin{equation}
\left\{ \begin{split} &\alpha +\gamma + 2\delta - 1 = 0  \\
&\beta - \delta = 0 \\
&\gamma + 2\delta + 1 = 0  \\ &\delta - 1 = 0
\end{split} \right.
\end{equation}
and its solutions are:
\begin{equation}
\left\{ \begin{split} &\alpha = 2\\ &\beta = 1 \\ &\gamma = -\, 3\\
&\delta = 1
\end{split} \right.
\end{equation}
The final result
\begin{equation}
u_{,\omega}(\theta ) \propto (k\theta \omega^{2})/c^{3}
\end{equation}
reproduces indeed, apart from numerical factors, the Rayleigh-Jeans's
formula. As to the Stephan-Boltzmann law one easily realizes that,
when removing $\hbar$, Eq.(8) leads to $r-s = 4-4 = 0$ and hence no
physical law is obtainable, in accordance with the fact that now the
energy density $u(\theta )$ blows up to infinity.\par\indent
Lastly, we wish to remind the reader that a gravitational constant
$G_{n}$ with dimensions $L^{n}M^{-\, 1}T^{-\, 2}$ is currently used in
higher dimensional theories of gravity; so, it seems quite consistent
that other universal constants such as $e$ or $\hbar$ be
extended to $n$ dimensions! One must take however into account the
following. We have seen that the photon energy is
$\varepsilon_{n} = \hbar_{n}\omega^{n-2}/c^{n-\, 3}$ and hence its
momentum is $|\boldsymbol{p}_{n}| =\varepsilon_{n}/c = \hbar_{n}
(\omega /c)^{n-2} = \hbar_{n} |\boldsymbol{k}|^{n-2}$; if now we pass
to  matter waves and make use of the de Broglie's relations, it
happens that only for $n = 3$ the connection is relativistically
invariant.\par\indent
Finally, we would like to add a brief comment about the
so-called ``black-hole thermodynamics'', invented by authors so as
Bekenstein ${}^{10}$ and Hawking~${}^{11}$ by generalizing the
ordinary thermodynamics.\par\indent
For instance, a black-body absorbs by definition every
electromagnetic radiation impinging on it, and was shown to emit a
characteristic (Planck) electromagnetic spectrum. By contrast, a
black-hole is defined as a body absorbing every kind of particle and
radiation impinging on it, and is predicted to emit a characteristic
spectrum for each kind of particle or radiation. Such spectra are
expected to be Planck-like, except that in their denominator the
minus sign (valid in general for boson emissions) has to be replaced
by a plus sign in the case of fermion emissions.\par\indent
These beautiful results have been reached by a clever, simultaneous
recourse to general relativity, quantum mechanics and thermodynamics.
They still lack, therefore, a definite theoretical foundation. From
this point of view, it may be interesting that they get some support
by simple dimensional considerations.

\begin{flushleft}
{\large{\textbf{Acknowledgments}}}\\ \vspace{0.2in}
We are grateful to Massimiliano Badino, Decio Pescetti and
Nadia Robotti for their comments and suggestions.
\end{flushleft}

\newpage


\section{INTRODUCTION}

A systematic examination of the manner in which data is described by
Monte Carlo models can give valuable insight into the underlying
physics involved.  However, many such comparisons are to a single
result, limiting the extent to which physics conclusions can be drawn
and raising the possibility that agreement is achieved in one cross
section at the expense of another. The aim of this study is to tune,
for the first time, general purpose Monte Carlo models to the existing
and expanding set of jet photoproduction data. To that end, we
have compared model predictions to published inclusive jet, di-jet and
3-jet data from the ZEUS
\cite{zeus1,zeus2,zeus3,zeus4,zeus5,zeus6,zeus7} and H1 \cite{h11,h12}
experiments at HERA, the OPAL \cite{opal1,opal2} experiment at LEP,
and the TOPAZ \cite{topaz} and AMY \cite{amy} experiments at TRISTAN.
The inclusion of both $\gamma p$ and $\gamma\gamma$ data (see figure
\ref{fig:lephera}) is particularly significant for constraining the
tuning.

Our tuning currently focuses on the role of the so-called
``underlying event'' in hadronic jet production and the extent to
which this can be described by perturbative QCD inspired models. We
constrain the multiparton interaction (minijet) models contained
within the Monte Carlos described below, and use these tunings to
estimate backgrounds at future colliders.

Such a tuning can also lead to a number of practical benefits.  The
method employed allows us to check the consistency of the data
contained in the various publications, which ranges across different
colliders, experiments, years, energies, kinematic regions etc. A
better description of the data by the models can be expected to lead
to a reduction in the systematic errors due to detector corrections in
future measurements.

\begin{figure}[htb]
\vspace{9pt}
{\includegraphics[width=7.5cm]{lephera.eps}}
\caption{Illustration of the correspondence between the different $\gamma\gamma$ and $\gamma p$ event classes.}
\label{fig:lephera}
\end{figure}

\section{THE MODELS}

The models tuned were HERWIG 5.9 \cite{hwg}, interfaced to the JIMMY
\cite{mi} eikonal model for multiparton interactions, and PYTHIA 6.125
\cite{pyt}, which contains a new (since version 5.7) mode aiming to
simulate the virtuality of the bremsstrahlung
photons \cite{friberg}. For HERWIG, a previous tuning to DIS data
(CLMAX=5.5, PSPLT=0.65) \cite{distune} was used.  In addition to the
default mode, preliminary investigations were carried out on a
modification to HERWIG, whereby the intrinsic transverse momentum
distribution of the partons in the photon takes the form of a power
law rather than a gaussian. This modified version, referred to as
HERWIG+$k_t$, was motivated by studies at HERA \cite{kt1} and at LEP
\cite{kt2} which indicate that it can lead to a better description of
the data.

The main free parameters investigated in the tuning are the
description of the ``underlying event'' and the the choice of photon
structure. Multiple hard scatters, illustrated in figure \ref{fig:mi},
can be enabled or not, and it is expected that some of the
distributions studied will exhibit great sensitivity to this
phenomenon. The other facet of the underlying event is the minimum
transverse momentum of the hard scatter(s), hereafter referred to as
$\hat{p}_T^{\rm min}$.  Sensitivity to this parameter is greatly increased when
multiple interactions are turned on. As one goes to lower $\hat{p}_T$, higher
parton densities are being probed, leading to an enhanced probability
that more than one hard scatter will occur. It should also be noted at
this point that, whilst HERWIG
employs a cluster fragmentation scheme, PYTHIA uses string hadronization.

A number of different parameterisations of the photon structure
function were investigated. Only leading order sets were used since
the matrix elements in the Monte Carlo models are leading order.  The
sets used were GRV 92 \cite{grv}, WHIT2 \cite{whit}, SaS1d and SaS2d
\cite{sas}, and LAC1 \cite{lac}. Throughout this work the GRV 94 LO
\cite{proton} proton pdf was used in the simulation of HERA data.

Finally, the overall normalisation of the Monte Carlo was treated as a
tunable parameter. This is justified within a range of around a factor
of two or less because of the uncertainty in the scale of $\alpha_s$.

\begin{figure}[htb]
\vspace{9pt}
\begin{center}
{\includegraphics[width=4.5cm]{mi.eps}}
\end{center}
\caption{An example of a multiple scattering.}
\label{fig:mi}
\end{figure}

\section{FITTING METHOD}

This work was carried out within the framework of a package developed
by the HERA community to permit easy comparison of data to Monte Carlo
\cite{hztool}. Extensions were made to facilitate the inclusion of the
$\gamma\gamma$ data.

The procedure for finding the best fit, for a given set of parameters,
of the Monte Carlo to the data was as follows. An overall $\chi^2$ per
degree of freedom across all the distributions (some 50 in all) was
defined as:
\begin{equation}
  \chi^2 = \frac{1}{n-1}\sum_{i=1}^{n}\frac{(\textrm{MC}(i) - \textrm{Data}(i))^2}{\sigma_{\textrm{\scriptsize{MC}}(i)}^2 + \sigma_{\textrm{\scriptsize{Data}}(i)}^2}
\end{equation}
where Data($i$) and MC($i$) are the values of the distribution in a given bin
$i$ for the data and Monte Carlo respectively. The sum runs over the total
number of bins $n$ in all the distributions.

The aim was to minimise this $\chi^2$, and a good fit is when it is
approximately one or less. The normalisation of the Monte Carlo was
varied to find the best fit (The sum in equation (1) is divided by
$n-1$ rather than $n$ to take into account the resulting loss of a
degree of freedom). For the HERA distributions, each data plot was
allowed to float within the quoted correlated error (typically of
15-20\%) and the number of plots was then subtracted from $n$ as
well. Otherwise, all systematic and statistical errors were treated as
uncorrelated and added in quadrature.

The distributions from TRISTAN were not included in the fitting
procedure. There appears to be a discrepancy between the TOPAZ and AMY
results. Indeed, this has been observed in a previous study
\cite{drees}.

\section{RESULTS AND DISCUSSION}

\begin{table}[hbt]
\newlength{\digitwidth} \settowidth{\digitwidth}{\rm 0}
\catcode`?=\active \def?{\kern\digitwidth}
\caption{HERWIG results (multiparton interactions on, {\bf not} $k_t$ version). The $\chi^2$ per degree of freedom for a particular parameter set is given first for the combined fits (see text), and then for the separate HERA and LEP fits.}
\label{tab:hrw}
\begin{tabular*}{7.5cm}{ccccc}
\hline
$p.d.f.$ & $\hat{p}_T^{\rm min}$ & $\chi^2/D.o.F.$ & HERA & LEP \\
\hline

WHIT 2 & 1.6 & 3.29 & 4.24 & 1.06 \\
       & 1.8 & 1.89 & 2.27 & 0.99 \\
       & 2.0 & 1.08 & 1.11 & 0.99 \\
       & 2.2 & 1.08 & 1.07 & 1.08 \\
       & 2.4 & 1.47 & 1.60 & 1.14 \\
GRV 92 & 1.6 &      &      & 1.03 \\
       & 1.8 & 1.47 & 1.71 & 0.89 \\
       & 2.0 & 1.35 & 1.46 & 1.09 \\
       & 2.2 &      &      & 1.17 \\
SaS 1d & 1.4 & 2.01 & 2.30 & 1.31 \\
       & 1.6 &      &      & 1.65 \\
       & 1.8 & 1.37 & 1.16 & 1.81 \\
SaS 2d & 1.6 & 1.57 & 1.81 & 0.98 \\
       & 1.8 & 1.28 & 1.40 & 0.97 \\
       & 2.0 &      &      & 1.31 \\
       & 2.2 & 1.27 & 1.25 & 1.31 \\
LAC 1  & 2.2 & 8.70 & 10.54& 4.35 \\
       & 2.6 &      & 2.14 &      \\

\hline
\end{tabular*}
\end{table}

\begin{table}[hbt]
\caption{PYTHIA results (multiparton interactions on). The $\chi^2$ per degree of freedom for a particular parameter set is given first for the combined fits (see text), and then for the separate HERA and LEP fits.}
\label{tab:pyt}
\begin{tabular*}{7.5cm}{ccccc}
\hline
$p.d.f.$ & $\hat{p}_T^{\rm min}$ & $\chi^2/D.o.F.$ & HERA & LEP \\
\hline
WHIT 2 & 2.4 & 1.17 & 0.91 & 1.86 \\
       & 2.6 & 1.13 & 0.99 & 1.50 \\
       & 2.8 & 1.18 & 1.07 & 1.44 \\
       & 3.0 & 1.24 & 1.12 & 1.52 \\
GRV 92 & 2.4 & 1.21 & 0.96 & 1.89 \\
       & 2.6 & 1.24 & 1.08 & 1.66 \\
       & 2.8 & 1.24 & 1.02 & 1.72 \\
SaS 1d & 1.4 & 1.15 & 1.15 & 1.15 \\
       & 1.6 & 1.15 & 1.19 & 1.03 \\
       & 1.8 & 1.35 & 1.38 & 1.23 \\
SaS 2d & 1.4 & 1.48 & 1.32 & 1.90 \\
       & 1.6 & 1.23 & 1.21 & 1.27 \\
       & 1.8 & 1.09 & 1.08 & 1.09 \\
       & 2.0 & 1.08 & 1.07 & 1.08 \\
       & 2.2 & 1.14 & 1.11 & 1.18 \\
LAC 1  & 2.8 & 6.18 & 5.08 & 10.0 \\

\hline
\end{tabular*}
\end{table}

The results of the fits using HERWIG ({\bf not} the $k_t$ version) are
shown in table \ref{tab:hrw}, and those for PYTHIA in table
\ref{tab:pyt}. Full details of all the fits, including the individual
plots, can be found at
http://www.hep.ucl.ac.uk/\~ ~jmb/HZTOOL/. Figure \ref{fig:results}
shows a summary of the results in graphical form. In most cases there
is a clear, favoured value of $\hat{p}_T^{\rm min}$, and overall - at least where
reasonable fits are found - there appears to be a favoured range of
$\sim1.6-2.2$ GeV. Note that although we do not constrain the $\hat{p}_T^{\rm min}$
to have the same value for both LEP and HERA fits, when we combine to
give an overall $\chi^2$ we have taken those runs with the same
$\hat{p}_T^{\rm min}$. This is motivated by the idea that, as one moves toward low
transverse momentum, two effects can render perturbative QCD
inapplicable.  Not only does $\alpha_s$ become large, but $x$ becomes
small, potentially leading to large $\ln(x)$ corrections. Our
conjecture is that if $\alpha_s$ effects are independent of the centre
of mass energy, then $\hat{p}_T^{\rm min}$ can be considered to be a universal
parameter. The low $x$ effects are modelled by multiple interactions,
and so the minijet model is ``universal'', but different
effects are seen
depending on the cms energy and beam particle type. All the results shown
have multiple
interactions enabled.  The agreement without multiple interactions is
very poor. A particularly striking example of this is shown in figure
\ref{fig:mi2}, and it is low $E_T^{Jet}$ measurements such as this which
are especially valuable for constraining multiple interaction models.

\begin{figure}[htb]
\vspace{9pt}
{\includegraphics[width=7.5cm]{fits2.eps}}
\caption{$\chi^2$/DoF for combined LEP and HERA data samples (see text).}
\label{fig:results}
\end{figure}

\begin{figure}[htb]
\vspace{9pt}
{\includegraphics[width=7.5cm]{mi2.eps}}
\caption{An example of the effect of multiple interactions. The data shown is ZEUS 1993 dijets data [2], along with the HERWIG predictions, using the GRV 92 photon and $\hat{p}_T^{\rm min} = 2.0$, with the JIMMY multiple interactions model either enabled or not.}
\label{fig:mi2}
\end{figure}

The scaling of the overall normalisation generally agrees between LEP and
HERA for the good
fits, and is in the range $1.2-1.3$. In addition, there is general
consistency between different datasets for both LEP and HERA. With
regard to pdf sets, WHIT and GRV lead to better fits for HERWIG,
whereas for PYTHIA there are no good fits to LEP data for sets other
than the SaS sets. In no circumstances does it appear possible to
obtain anything close to a good description of the data with LAC1.

The DIS tuned parameters for HERWIG do help, if only marginally, but
there is no firm evidence yet that the HERWIG+$k_t$ modification
improves matters.

\section{LINEAR COLLIDER}

The results of the tuning can be used to estimate photon-induced
minijet backgrounds at future colliders. The best fits found using
HERWIG have been extrapolated to linear collider energies ($\sqrt s =
500$ GeV), with the TESLA beamstrahlung and bremsstrahlung spectra
included \cite{circe}. The parameters thus used were the WHIT 2 photon
pdf and a $\hat{p}_T^{\rm min}$ of 2.0 GeV with multiple interactions enabled. The
overall normalisation was scaled by a factor of 1.4. The estimate for
the minijet transverse momentum cross-section is shown in figure
\ref{fig:lc1}.

The results are consistent at the 50\% level with independent studies
using PYTHIA and the DELPHI Monte Carlo \cite{sitges}. They indicate
that the backgrounds will not be a concern for detector occupancy and
dosage, but that minijets potentially present a very significant
background to physics owing to the high $\hat{p}_T$ tail.

\begin{figure}[htb]
\vspace{9pt}
{\includegraphics[width=7.5cm]{hrwlc_whit_2.0b.eps}}
\caption{Jet transverse energy differential cross section in the
barrel (top) and endcap (lower) calorimeters of the TESLA conceptual
design report detector \cite{CDR}, for the high luminosity TESLA option.}
\label{fig:lc1}
\end{figure}

\section{CONCLUSIONS}

The general purpose Monte Carlo models HERWIG and PYTHIA have been
tuned to jet photoproduction data from LEP and HERA. An optimal range
of hard scales appears to emerge, and some level of multiparton
interactions is found to be necessary in these models in order that an
adequate description of the data be achieved. Favoured sets of parton
distribution functions are established for each generator but
currently no generator independent conclusion on photon structure is
drawn. The HERWIG tuned parameters are used to estimate photon-induced
backgrounds at the linear collider. This work is ongoing and
up-to-date results can be found on our web page
http://www.hep.ucl.ac.uk/\~ ~jmb/HZTOOL/.