%Paper: 
%From: bf08@ibm3090.rz.uni-karlsruhe.de
%Date: Thu, 20 Oct 94 15:13

\documentstyle{ichep}
\hoffset 0mm
\voffset 3mm
% Use suitable offset values to centre the output on your own system.
%
\begin{document}

\title{
\vskip -2.8cm
\hfill {\normalsize TTP94--19}\\
\vskip -0.2cm
\hfill {\normalsize September 1994}\\
\vskip 1.0cm
QCD Corrections to Inclusive Distributions of Leptons\\
in Decays of Polarised Heavy Quarks*}
%
\author{
\vskip-0.5cm
Marek Je\.zabek}

\affil{Institute of Nuclear Physics,
        ul.Kawiory 26a, PL-30055 Cracow, Poland\\
and\\
Institut f. Theoretische Teilchenphysik, Univ. Karlsruhe
D-76128 Karlsruhe, Germany}
%
\abstract{
Compact analytic expressions have been obtained for the first order
perturbative QCD corrections to the inclusive spectra of the leptons
in the semileptonic decays of polarised heavy quarks.
\cabs
Charmed and beautiful $\Lambda$ baryons from $Z^0$ decays can be viewed
as sources of highly polarised charm and bottom quarks. Charged leptons
and neutrinos from $\Lambda_b$ and $\Lambda_c$ decays can be used in
the polarisation studies for the corresponding heavy quarks. Thus our
results are applicable for the b quark polarisation measurements at LEP.
\cabs
Short lifetime  enables polarisation studies for the top quark.
The angular-energy spectra
of the charged leptons are particularly useful in this respect whereas
the distributions of the neutrinos are sensitive to deviations from
the V-A structure of the charged weak current in the decay.
}

\twocolumn[\maketitle]

%\fnm{7}{Work partly supported by KBN under contract
%            26 and by DFG under contract 436POL173193S.}
%\fnm{2}{E-mail: jezabek@chopin.ifj.edu.pl}
\fnm{7}{Presented at XXVII International Conference on High Energy
Physics, 20-27 July 1994, Glasgow, Scotland; to appear in the
proceedings.}

\section{Introduction}
Inclusive semileptonic decays
of polarised charm and bottom quarks play important
role in present day particle physics. With
increasing statistics at LEP and good prospects for $B$-factories
quantitative description of
these processes may offer the most
interesting tests of the standard quantum theory
of particles. In fact the first measurement of $b$ quark
polarisation at LEP has been presented by the ALEPH collaboration
at this conference \cite{Roudeau}.
At the high energy frontier semileptonic decays
of the top quark will be instrumental in establishing its
properties \cite{Kuehn1,Kuehn2,Kuehn3,teupitz}.

In this article I present the results of calculations
of the first order perturbative
QCD corrections to semileptonic
decays of polarised heavy quarks. Some of these results
have been published in \cite{CJK,CJKK}.
In \cite{CJ} compact analytic formulae have been obtained
for the distributions of the charged lepton and the neutrino.
These formulae  agree with those given in \cite{CJK}
for the joint angular and energy distribution
of the charged lepton in top quark decays and are much simpler.

\section{The formula and cross checks}
The QCD corrected triple differential distribution
of the charged lepton
for the semileptonic decay  of the polarised quark
with the weak isospin $I_3=\pm 1/2$
can be written in the following way \cite{CJ}:
\begin{eqnarray}
{{\rm d}\Gamma^{\pm} \over {\rm d}x\,{\rm d}y\,{\rm d}\cos\theta  }
&& \sim\qquad \left[\,
{\rm F}^\pm_0(x,y) + S\cos\theta\,{\rm J}^\pm_0(x,y)\,
\right]
\nonumber\\
&& -\; {2\alpha_s\over3\pi}\;
  \left[\,
  {\rm F}^\pm_1(x,y) + S\cos\theta\,{\rm J}^\pm_1(x,y)\,
  \right]
\nonumber\\
\label{eq:1}
\end{eqnarray}
In the rest frame of the decaying heavy quark
$\theta$ denotes the angle between the
polarisation vector $\vec s$ of the heavy quark and the
direction of the charged lepton,
$S=|\,\vec s\,|$,
$x= 2Q\ell/Q^2$ and $y= 2\ell\nu/Q^2$ where
$Q$, $\ell$ and $\nu$ denote the four-momenta of the decaying
quark, charged lepton and neutrino. Eq.(\ref{eq:1}) describes also
the triple differential
distribution of the neutrino for  $I_3=\mp 1/2$. In this case,
however, $x= 2Q\nu/Q^2$  and
$\theta$ denotes the angle between
$\vec s$ and the three-momentum of the neutrino.
The functions ${\rm F}^\pm_0(x,y)$ and ${\rm J}^\pm_0(x,y)$
corresponding to Born approximation read:
\begin{eqnarray}
{\rm F}^+_0(x,y) &=& x (x_m-x)
\label{eq:5}\\
{\rm J}^+_0(x,y) &=& {\rm F}^+_0(x,y)
\label{eq:6}\\
{\rm F}^-_0(x,y) &=& (x-y) (x_m-x+y)
\label{eq:7}\\
{\rm J}^-_0(x,y) &=& (x-y) (x_m-x+y-2y/x)
\label{eq:8}
\end{eqnarray}
where $x_m=1-\epsilon^2$, $\epsilon^2= q^2/Q^2$ and $q$
denotes the four-momentum of the quark originating from
the decay.
The functions ${\rm F}^\pm_1(x,y)$ and
${\rm J}^+_1(x,y)$ correspond to
the first order QCD corrections and
are given in \cite{CJ}.

Non-trivial cross checks are fulfilled by the polarisation
independent parts of the distributions (\ref{eq:1}):
\begin{itemize}
\item
the distributions
${\rm d}\Gamma^{\pm} / {\rm d}x\,{\rm d}y$
agree with the results
for unpolarised decays which were
obtained in \cite{JK}. The present formulae are simpler.
\item
in the four-fermion (Fermi) limit integration over $y$
can be performed numerically. The resulting distributions
${\rm d}\Gamma^{\pm} / {\rm d}x$
also agree with those of \cite{JK}.
Recently the results of \cite{JK}
have been confirmed \cite{CCM1}. Thus an old conflict
with other calculations \cite{CCM0} is solved and
the agreement with \cite{JK} can be considered as a
non-trivial cross check. Moreover, the
analytic result of \cite{JK} for
${\rm d}\Gamma^{+} / {\rm d}x$
and $\epsilon=0$
has been also confirmed \cite{GWF}.
\item
$${{\rm d}\Gamma^{+} / {\rm d}y} =
{{\rm d}\Gamma^{-} / {\rm d}y}$$
and the analytic formula for this distribution
exists \cite{JK0} which at the same time describes the lifetime
of the top quark as a function of its mass.
This formula has been confirmed by a few groups, c.f. \cite{JK93}
and references therein.
\item
in the four-fermion limit the result for the total rate $\Gamma$
derived from
eq.(\ref{eq:1})
agrees with the results of \cite{CM} and
the analytical formula of \cite{nir}.
\end{itemize}

\section{Applications}
\subsection{Polarised bottom and charm quarks}

Polarisation studies for heavy flavors at LEP are a new
interesting field of potentially fundamental significance,
see \cite{Mele,Roudeau} for recent reviews.
According to the Standard Model
$Z^0\to b\bar b$ and
$Z^0\to c\bar c$ decays
can be viewed as sources of highly polarised heavy quarks.
The degree of longitudinal polarisation is fairly large,
amounting to $\langle P_b\rangle = -0.94$ for $b$
and $\langle P_c\rangle = -0.68$ for $c$ quarks~\cite{Kuehn1}.
The polarisations depend weakly on the production angle.
QCD corrections to Born result are about 3\% \cite{KPT}.
The real drawback is that due to hadronisation the net longitudinal
polarisation of the decaying  $b$ and $c$ quarks is drastically
decreased. In particular these $b$ quarks become depolarised
which are bound in $B$ mesons
both produced directly and from $B^*\to B\gamma$ transitions.
The signal is therefore significantly reduced.
Only those $b$'s (a few percent) which fragment directly
into $\Lambda_b$  baryons retain information on the original
polarisation \cite{Bjo}.
Polarisation transfer from a heavy
quark $Q$ to the corresponding $\Lambda_Q$ baryon is 100\%
\cite{CKPS} at least in the limit $m_Q\rightarrow\infty$.
Thus, a large net polarisation is expected for heavy quarks
in samples enriched with these heavy baryons.

It has been proposed long ago \cite{zerwas}
that distributions of charged leptons from semileptonic
decays of beautiful hadrons can be used in polarisation
studies for $b$ quarks.
Some advantages of neutrino distributions have been
also pointed out \cite{CJKK,BR,dittmar}.
Recently there has been considerable progress in the theory of
the inclusive semileptonic decays of heavy flavor hadrons.
It has been shown that in the leading order of an expansion
in inverse powers of heavy quark mass $1/m_Q$ the
spectra for hadrons coincide with those for the decays of
free heavy quarks \cite{CGG} and there are no $\Lambda_{QCD}/m_Q$
corrections to this result away from the energy
endpoint.
$\Lambda^2_{QCD}/m_Q^2$ corrections have been
calculated in \cite{Bigi,wise} for $B$ mesons and
in \cite{wise} for polarised $\Lambda_b$ baryons.
For some decays the results are similar to those of the well-known
$ACCMM$ model \cite{altar}.
The corrections to charm decays are larger than for bottom
and convergence of
$1/m_Q$  expansion is poorer \cite{Shifman}.
Perturbative first order QCD corrections contribute 10-20\%
to the semileptonic decays and for bottom are much larger
than the nonperturbative ones.

\subsection{Polarised top quarks}
The analysis of polarised top quarks and their decays
has recently attracted considerable attention,
see \cite{Kuehn3,teupitz} and references
cited therein.
The reason is that this analysis
will result in determination of the top
quark coupling to the $W$ and $Z$ bosons either
confirming the predictions
of the Standard Model
or providing clues for physics beyond.
The latter possibility is particularly intriguing for the top quark
because $m_t$
plays an exceptional role in the fermion mass spectrum.

A number of mechanisms have been suggested that will
lead to polarised top quarks.
Studies at a linear electron-positron collider are particularly
clean for precision tests. However, also
$\gamma\gamma$ collisions with circular polarised photons
and subsequent spin analysis of top quarks might reveal new
information.
Related studies may be performed in hadronic collisions
which in this case are mainly based on the correlation
between $t$ and $\bar t$ decay products.
However, single top production through $Wb$ fusion at LHC
may also be a useful source of polarised top quarks.
Electron-positron collisions are
the most efficient and flexible
reactions producing
polarised top quarks. A small
component of polarisation transverse to the production plane
is induced by final state interactions.
The longitudinal polarisation $P_L$
is large.  $P_L$ varies strongly with the production angle.
Averaging over the production angle
leads therefore to a significant reduction of $P_L$ with typical
values of $\langle P_L\rangle$ around -0.2 \cite{KPT}.\\
All these reactions lead to sizable polarisation and can be used
to obtain information on the production mechanism.
However, two drawbacks are evident:
production and decay are mixed in an intricate manner, and
furthermore the degree of polarisation is relatively small
and depends on the production angle. Top quark production
with longitudinally polarised electron beams and close to
threshold provides one important exception: the restricted
phase space leads to an amplitude which is dominantly S-wave
such that the electron (and positron) spin is directly
transferred to the top quark. Close to threshold and with
longitudinally polarised electrons
one can study decays of polarised top quarks
under particularly convenient conditions:
large event rates, well identified rest frame of the top quark,
and large degree of polarisation.  Moreover, short lifetime
of top quark practically eliminates nonperturbative corrections
due to hadronisation.

In the rest frame of the decaying $t$ quark distributions
of the decay products are sensitive to its polarisation.
Eq.(\ref{eq:6}) implies that in Born approximation
the double differential angular-energy distribution
of the charged lepton is the product of the energy
distribution and the angular distribution. The latter distribution
is of the following form
\begin{equation}
{ {\rm d}N\over{\rm d}\cos\theta} =
{1\over 2}\, \left[\, 1\,+
\,S\cos\theta \right]
\label{eq:elec1}
\end{equation}
QCD corrections essentially do not spoil factorisation
of the charge lepton distribution \cite{CJK}.
It is noteworthy that for $S$=1 the angular
dependence in (\ref{eq:elec1}) is maximal because
any larger coeffecient multiplying $\cos\theta$ would be
in conflict with positivity of the decay rate.
Thus the polarisation analysing power of
the charged lepton energy-angular distribution
is maximal \footnote{This is reversed
for $b$ decays and the polarisation analysing power is maximal
for the neutrino distributions
because the formulae for the neutrino distributions
in down-type quark decay describe the charged lepton distributions
for an up-type quark.}.

It follows from eqs.
(\ref{eq:7}) and (\ref{eq:8}) that already in Born approximation
there is no factorisation
for the neutrino energy-angular distribution.
Neutrino distributions are therefore less sensitive to the
polarisation of the decaying top quark than charge lepton
distributions
On the other hand it has been shown \cite{JK94} that
the angular-energy distribution of neutrinos
from the polarised top quark decay will allow for a particularly
sensitive test of the V-A structure of the weak charged current.
%\footnote{A test of V-A charge current in $b$ decays based
%on simultaneous analysis of charged leptons and neutrinos
%has been discussed in \cite{dittmar}.}.
The effect of QCD correction can mimic a small admixture
of V+A interaction. Therefore, inclusion of the radiative
QCD correction to the decay distributions  is necessary for
a quantitative study.

\vskip0.5cm\noindent
{\bf Acknowledgements}\\
I thank  Andrzej Czarnecki, Hans K\"uhn
and J\"urgen K\"orner for collaborations
on research reported in this article.
I would like to gratefully acknowledge helpful
correspondence with Professors N. Cabibbo,
G. Corbo and L. Maiani.
\par\noindent
This
work is partly supported by KBN under contract
26 and by DFG under contract 436POL173193S.

\Bibliography{9}
\bibitem{Roudeau} P.\ Roudeaud, {\it Heavy Quark Physics}, in
these proceedings.
\bibitem{Kuehn1}
J.H. K\"uhn and P.M. Zerwas, in {\it Heavy Flavours}, eds. A.J. Buras
and M. Lindner, (World Scientific, Singapore, 1992), p.434.
\bibitem{Kuehn2}
J.H. K\"uhn et al., DESY Orange Report 92-123A (1992), vol.I,p.255.
\bibitem{Kuehn3}
J.H. K\"uhn, "Top Quark at a Linear Collider", in {\it Physics
and Experiments with Linear $e^+e^-$ Colliders}, eds. F.A. Harris
et al., (World Scientific, Singapore, 1993), p.72.
\bibitem{teupitz}
M. Je\.zabek, {\it Top Quark Physics}, in proceedings
of Zeuthen workshop {\it Physics at LEP 200 and Beyond}, to appear
in \np{B}{94}{Suppl.};\ Karlsruhe\ preprint TTP94-09.
\bibitem{CJK}
A.~Czarnecki, M.~Je{\.z}abek and J.H. K{\"u}hn,
\np{B351}{91}{70}.
\bibitem{CJKK}
A. Czarnecki, M. Je\.zabek, J.G. K\"orner and J.H. K{\"u}hn,
\prl{73}{94}{384}.
\bibitem{CJ}
A. Czarnecki and M. Je\.zabek,
preprint TTP 93-40, Karlsruhe, 1994, 
\np{B}{94}{in print}.
\bibitem{JK}
M.~Je{\.z}abek and J.H. K{\"u}hn,
\np{B320}{89}{20}.
\bibitem{CCM1}
N. Cabibbo, G. Corbo and L. Maiani, private communication.
\bibitem{CCM0}
N. Cabibbo, G. Corbo and L. Maiani,
\np{B155}{79}{93};\
G. Corbo, \np{B212}{83}{99}.
\bibitem{GWF}
C. Greub, D. Wyler and W. Fetscher,
\pl{B324}{94}{109}.
\bibitem{JK0}
M.~Je{\.z}abek and J.H. K{\"u}hn,
\np{B314}{89}{1}.
\bibitem{JK93}
M. Je\.zabek and J.H. K\"uhn, \prev{D48}{93}{R1910}.
\bibitem{CM}
N. Cabibbo and L. Maiani,
\pl{B79}{78}{109}.
\bibitem{nir} Y. Nir,
\pl{B221}{89}{184}.
\bibitem{Mele}
B. Mele, preprint n.1009, Rome, 1994.
\bibitem{KPT}
J.G. K\"orner, A. Pilaftsis and M.M. Tung, preprint MZ-TH/93-3,
\zp{C}{94}{in print}.
\bibitem{Bjo}
J.D. Bjorken, \prev{D40}{89}{1513}.
\bibitem{CKPS}
F.E. Close, J.G. K\"orner, R.J.N. Phillips and
D.J. Summers, \jpg{18}{92}{1716}.
\bibitem{zerwas}
G. K\"opp, L.M. Sehgal and P.M. Zerwas,
\np{B123}{1977}{77};\\
B. Mele and G. Altarelli, \pl{B299}{93}{345}.
\bibitem{BR}
G. Bonvicini and L. Randall,
\prl{73}{94}{392}.
\bibitem{dittmar}
M. Dittmar and Z. W{\c a}s,
\pl{B332}{94}{168};\\
M. Dittmar, in these proceedings.
\bibitem{CGG}
J. Chay, H. Georgi and B. Grinstein,
\pl{B247}{90}{399}.
\bibitem{Bigi}
I.~Bigi, M.~Shifman, N.~Uraltsev and A.~Vainshtein,
\prl{71}{93}{496}.
\bibitem{wise}
A.V. Manohar and M.B. Wise, \prev{D49}{94}{1310}.
\bibitem{altar}
G. Altarelli et al, \np{B208}{82}{365}.
\bibitem{Shifman}
M. Shifman, in these proceedings.
\bibitem{JK94}
M. Je{\.z}abek and J.H. K{\"u}hn,
\pl{B329}{94}{317}.
%
\end{thebibliography}
%
\end{document}
%
%  ichep.sty
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%                                                                      %%%
%%%    INSTITUTE OF PHYSICS PUBLISHING                                   %%%
%%%                                                                      %%%
%%%    ICHEP.STY  LaTeX CRC style file for ICHEP94 conference            %%%
%%%                                                                      %%%
%%%    Copyright 1994 IOP Publishing Ltd                                 %%%
%%%                                                                      %%%
%%%    Permission is hereby given to use this file for                   %%%
%%%    material to be submitted to or published by                       %%%
%%%    Institute of Physics Publishing                                   %%%
%%%                                                                      %%%
%%%    Version 1.2  1 July 1994                                          %%%
%%%                                                                      %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% First we have a character check
%
% ! exclamation mark    " double quote
% # hash                ` opening quote (grave)
% & ampersand           ' closing quote (acute)
% $ dollar              % percent
% ( open parenthesis    ) close paren.
% - hyphen              = equals sign
% | vertical bar        ~ tilde
% @ at sign             _ underscore
% { open curly brace    } close curly
% [ open square         ] close square bracket
% + plus sign           ; semi-colon
% * asterisk            : colon
% < open angle bracket  > close angle
% , comma               . full stop
% ? question mark       / forward slash
% \ backslash           ^ circumflex
%
% ABCDEFGHIJKLMNOPQRSTUVWXYZ
% abcdefghijklmnopqrstuvwxyz
% 1234567890
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
\typeout{Document Style `ichep.sty'. IOP camera-ready copy
style file for ICHEP Conference Proceedings}

\let\reset@font\empty

\def\refname{References}
\def\figurename{Figure}
\def\tablename{Table}
\def\abstractname{Abstract}

\def\@ptsize{0}
\@namedef{ds@11pt}{\def\@ptsize{0}}
\@namedef{ds@12pt}{\def\@ptsize{0}}
\def\ds@twoside{\@twosidetrue
           \@mparswitchtrue}

\def\ds@draft{\overfullrule 5\p@}

\newif\if@titlepage \@titlepagefalse
\def\ds@titlepage{\@titlepagetrue}

\def\ds@twocolumn{\@twocolumntrue}

%
\newdimen\mathindent
\newlength{\digitwidth}
\newlength{\indentedwidth}
\newcounter{firstpage}
\newbox{\captionbox}
\newcounter{eqnval}
%
\@twosidetrue
%\@mparswitchtrue
\def\ds@draft{\overfullrule 5\p@}
\@options
%
\def\hexnumber@#1{\ifcase#1 0\or 1\or 2\or 3\or 4\or 5\or 6\or 7\or 8\or
 9\or A\or B\or C\or D\or E\or F\fi}
%
\lineskip 1pt \normallineskip 1pt
\def\baselinestretch{1}
\def\@normalsize{\@setsize\normalsize{12pt}\xpt\@xpt
\abovedisplayskip 10pt plus2pt minus5pt
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus3pt
\belowdisplayshortskip 6pt plus3pt minus3pt}
\def\small{\@setsize\small{11pt}\ixpt\@ixpt
\abovedisplayskip 8.5pt plus 3pt minus 4pt
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus2pt
\belowdisplayshortskip 4pt plus2pt minus 2pt
\def\@listi{\topsep 4pt plus 2pt minus 2pt\parsep 0pt plus 1pt
\itemsep \parsep}}
\def\footnotesize{\@setsize\footnotesize{9.5pt}\viiipt\@viiipt
\abovedisplayskip 6pt plus 2pt minus 4pt
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus 1pt
\belowdisplayshortskip 3pt plus 1pt minus 2pt
\def\@listi{\topsep 3pt plus 1pt minus 1pt\parsep 0pt plus 1pt
\itemsep \parsep}}
\def\scriptsize{\@setsize\scriptsize{8pt}\viipt\@viipt}
\def\tiny{\@setsize\tiny{6pt}\vpt\@vpt}
\def\large{\@setsize\large{14pt}\xiipt\@xiipt}
\def\Large{\@setsize\Large{18pt}\xivpt\@xivpt}
\def\LARGE{\@setsize\LARGE{22pt}\xviipt\@xviipt}
\def\huge{\@setsize\huge{25pt}\xxpt\@xxpt}
\def\Huge{\@setsize\Huge{30pt}\xxvpt\@xxvpt}
\normalsize
\oddsidemargin -3pc
\evensidemargin -3pc
\marginparwidth .75in
\marginparsep 7\p@
\topmargin=-72\p@
\headheight=12\p@
\headsep=12\p@
\footheight=12\p@
\footskip=25\p@

\textheight 55pc
\textwidth 18cm
\indentedwidth 15.9cm
\columnsep 1cm
\columnseprule 0\p@
\mathindent = 2pc

\newcommand{\onecol}{\parfillskip=0pt\par\eject
   \onecolumn\parfillskip=0pt plus1fil\noindent}
\newcommand{\twocol}{\parfillskip=0pt\par\eject
   \twocolumn\parfillskip=0pt plus1fil\noindent}
\footnotesep 6.65\p@
\skip\footins 9\p@ plus 4\p@ minus 2\p@
\floatsep 12\p@ plus 2\p@ minus 2\p@
\textfloatsep 18\p@ plus 2\p@ minus 4\p@
\intextsep 12\p@ plus 2\p@ minus 2\p@
\@maxsep 20\p@
\dblfloatsep 12\p@ plus 2\p@ minus 2\p@
\dbltextfloatsep 18\p@ plus 2\p@ minus 4\p@
\@dblmaxsep 20\p@
\@fptop 0\p@
\@fpsep 8\p@ plus 1fil
\@fpbot 0\p@ plus 1fil
\@dblfptop 0\p@
\@dblfpsep 8\p@ plus 1fil
\@dblfpbot 0\p@
\marginparpush 5\p@

\parskip 0\p@
\parindent 16\p@
\topsep 4\p@ plus 2\p@ minus 2\p@
\partopsep 2\p@ plus 1\p@ minus 1\p@
\itemsep 0\p@ plus 2\p@
\@lowpenalty 51
\@medpenalty 151
\@highpenalty 301
\@beginparpenalty -\@lowpenalty
\@endparpenalty -\@lowpenalty
\@itempenalty -\@lowpenalty

\@noskipsecfalse   % new version

%
\def\section{\@startsection{section}{1}{\z@}{-3.5ex plus -1ex minus
 -.2ex}{2.3ex plus .2ex}{\noindent\reset@font\normalsize\bf\raggedright}}
\def\subsection{\@startsection{subsection}{2}{\z@}{-3.25ex plus -1ex minus
 -.2ex}{1.5ex plus .2ex}{\noindent\reset@font
  \normalsize\it\raggedright\nohyphens}}
\def\subsubsection{\@startsection{subsubsection}{3}{\z@}{-3.25ex plus
-1ex minus -.2ex}{-1em}{\reset@font\normalsize\it\nohyphens}}
\def\paragraph{\@startsection
 {paragraph}{4}{\z@}{3.25ex plus 1ex minus .2ex}{-1em}
                                                {\reset@font\normalsize\it}}
\def\subparagraph{\@startsection
 {subparagraph}{4}{\parindent}{3.25ex plus 1ex minus
 .2ex}{-1em}{\reset@font\normalsize\it}}

\def\@sect#1#2#3#4#5#6[#7]#8{\ifnum #2>\c@secnumdepth
     \let\@svsec\@empty\else
     \refstepcounter{#1}\edef\@svsec{\csname the#1\endcsname.\hskip 1em}\fi
     \@tempskipa #5\relax
      \ifdim \@tempskipa>\z@
        \begingroup #6\relax
          \noindent{\hskip #3\relax\@svsec}{\interlinepenalty \@M #8\par}%
        \endgroup
       \csname #1mark\endcsname{#7}\addcontentsline
         {toc}{#1}{\ifnum #2>\c@secnumdepth \else
                      \protect\numberline{\csname the#1\endcsname}\fi
                    #7}\else
        \def\@svsechd{#6\hskip #3\relax  %% \relax added 2 May 90
                   \@svsec #8\csname #1mark\endcsname
                      {#7}\addcontentsline
                           {toc}{#1}{\ifnum #2>\c@secnumdepth \else
                             \protect\numberline{\csname the#1\endcsname}\fi
                       #7}}\fi
     \@xsect{#5}}
%
\def\@ssect#1#2#3#4#5{\@tempskipa #3\relax
   \ifdim \@tempskipa>\z@
     \begingroup #4\noindent{\hskip #1}{\interlinepenalty \@M #5\par}\endgroup
   \else \def\@svsechd{#4\hskip #1\relax #5}\fi
    \@xsect{#3}}


\setcounter{secnumdepth}{3}

\def\appendix{\@@par
 \setcounter{section}{0}
 \setcounter{subsection}{0}
 \setcounter{subsubsection}{0}
 \setcounter{equation}{0}
 \setcounter{figure}{0}
 \setcounter{table}{0}
 \def\thesection{Appendix \Alph{section}}
 \def\theequation{\ifnumbysec
      \Alph{section}.\arabic{equation}\else
      \Alph{section}\arabic{equation}\fi}
 \def\thetable{\ifnumbysec
      \Alph{section}\arabic{table}\else
      A\arabic{table}\fi}
 \def\thefigure{\ifnumbysec
      \Alph{section}\arabic{figure}\else
      A\arabic{figure}\fi}}



\labelsep 4\p@

\leftmargini 16\p@
\leftmarginii 18\p@
\leftmarginiii 16\p@
\leftmarginiv 14\p@
\leftmarginv 10\p@
\leftmarginvi 10\p@
\leftmargin\leftmargini
\labelwidth\leftmargini\advance\labelwidth-\labelsep
\parsep 0\p@ plus 1\p@
\def\@listI{\leftmargin\leftmargini \parsep 4\p@ plus2\p@ minus\p@
\topsep 8\p@ plus2\p@ minus4\p@
\itemsep 4\p@ plus2\p@ minus\p@}

\let\@listi\@listI
\@listi

\def\@listii{\leftmargin\leftmarginii
 \labelwidth\leftmarginii\advance\labelwidth-\labelsep
 \topsep 3\p@ plus 1\p@ minus 1\p@
 \parsep 0\p@ plus 1\p@
 \itemsep \parsep}
\def\@listiii{\leftmargin\leftmarginiii
 \labelwidth\leftmarginiii\advance\labelwidth-\labelsep
 \topsep 2\p@ plus 1\p@ minus 1\p@
 \parsep \z@ \partopsep 1\p@ plus 0\p@ minus 1\p@
 \itemsep \topsep}
\def\@listiv{\leftmargin\leftmarginiv
 \labelwidth\leftmarginiv\advance\labelwidth-\labelsep}
\def\@listv{\leftmargin\leftmarginv
 \labelwidth\leftmarginv\advance\labelwidth-\labelsep}
\def\@listvi{\leftmargin\leftmarginvi
 \labelwidth\leftmarginvi\advance\labelwidth-\labelsep}

\pretolerance=5000
\tolerance=8000
\hbadness=5000
\vbadness=5000
%
\def\labelenumi{\theenumi}
\def\theenumi{\arabic{enumi}}
\def\labelenumii{\theenumii}
\def\theenumii{\alpha{enumii}}
\def\p@enumii{\theenumi.}
\def\labelenumiii{\theenumiii.}
\def\theenumiii{\arabic{enumiii}}
\def\p@enumiii{\p@enumii.\theenumii}
\def\labelenumiv{\theenumiv.}
\def\theenumiv{\arabic{enumiv}}
\def\p@enumiv{\p@enumiii.\theenumiii}

\def\labelitemi{$\m@th\bullet$}
\def\labelitemii{\bf --}
\def\labelitemiii{$\m@th\ast$}
\def\labelitemiv{$\m@th\cdot$}

\def\verse{\let\\=\@centercr
 \list{}{\itemsep\z@ \itemindent -1.5em\listparindent \itemindent
 \rightmargin\leftmargin\advance\leftmargin 1.5em}\item[]}
\let\endverse\endlist
\def\quotation{\list{}{\listparindent 1.5em
 \itemindent\listparindent
 \rightmargin\leftmargin\parsep 0\p@ plus 1\p@}\item[]}
\let\endquotation=\endlist
\def\quote{\list{}{\rightmargin\leftmargin}\item[]}
\let\endquote=\endlist

\def\descriptionlabel#1{\hspace\labelsep \bf #1}
\def\description{\list{}{\labelwidth\z@ \itemindent-\leftmargin
 \let\makelabel\descriptionlabel}}
\let\enddescription\endlist
%
\def\enumerate{\ifnum \@enumdepth >3 \@toodeep\else
      \advance\@enumdepth \@ne
      \edef\@enumctr{enum\romannumeral\the\@enumdepth}\list
      {\csname label\@enumctr\endcsname}{\usecounter
        {\@enumctr}\def\makelabel##1{##1\hss}}\fi}
%
\def\itemize{\ifnum \@itemdepth >3 \@toodeep\else \advance\@itemdepth \@ne
\edef\@itemitem{labelitem\romannumeral\the\@itemdepth}%
\list{\csname\@itemitem\endcsname}{\def\makelabel##1{##1\hss}\topsep=3pt
  \parsep=0pt\listparindent=0pt\itemsep=0pt\partopsep=0pt\rightmargin=0pt
  }\fi}
%
\newenvironment{leqnarray}{\begin{leqnarray}}{\end{leqnarray}}
\def\leqnarray{\stepcounter{equation}\let\@currentlabel=\theequation
\global\@eqnswtrue
\global\@eqcnt\z@\tabskip\mathindent\let\\=\@eqncr
\abovedisplayskip\topsep\ifvmode\advance\abovedisplayskip\partopsep\fi
\belowdisplayskip\abovedisplayskip
\belowdisplayshortskip\abovedisplayskip
\abovedisplayshortskip\abovedisplayskip
$$\halign to
\columnwidth\bgroup\@eqnsel$\displaystyle\tabskip\z@
 {##{}}$&\global\@eqcnt\@ne
                    $\displaystyle{{}##{}}$\hfil    %\hfil delete before 2nd $
 &\global\@eqcnt\tw@ $\displaystyle{{}##}$\hfil
 \tabskip\@centering&\llap{##}\tabskip\z@\cr}
%
\def\endleqnarray{\@@eqncr\egroup
 \global\advance\c@equation\m@ne$$\global\@ignoretrue }
%
\arraycolsep 5\p@
\tabcolsep=6\p@
\arrayrulewidth .4\p@
\doublerulesep 2\p@
\tabbingsep \labelsep
\skip\@mpfootins = \skip\footins
\fboxsep = 3\p@
\fboxrule = .4\p@
\def\titlepage{\@restonecolfalse\if@twocolumn\@restonecoltrue\onecolumn
     \else \newpage \fi \thispagestyle{myheadings}\c@page\z@}

\def\endtitlepage{\if@restonecol\twocolumn \else \newpage \fi}

\newcounter {section}
\newcounter {subsection}[section]
\newcounter {subsubsection}[subsection]
\newcounter {paragraph}[subsubsection]
\newcounter {subparagraph}[paragraph]



\def\thesection {\arabic{section}}
\def\thesubsection {\thesection.\arabic{subsection}}
\def\thesubsubsection {\thesubsection .\arabic{subsubsection}}
\def\theparagraph {\thesubsubsection.\arabic{paragraph}}
\def\thesubparagraph {\theparagraph.\arabic{subparagraph}}
\def\@chapapp{Section}


\def\@pnumwidth{1.55em}
\def\@tocrmarg {2.55em}
\def\@dotsep{4.5}
\setcounter{tocdepth}{2}


\def\tableofcontents{\@restonecolfalse\if@twocolumn\@restonecoltrue
 \onecolumn\fi\section*{Contents}{}\thispagestyle{empty}
 \@starttoc{toc}\if@restonecol\twocolumn\fi}
%
\def\l@section{\@dottedtocline{1}{1.5em}{2.3em}}
\def\l@subsection{\@dottedtocline{2}{3.8em}{3.2em}}
\def\l@subsubsection{\@dottedtocline{3}{7.0em}{4.1em}}
\def\l@paragraph{\@dottedtocline{4}{10em}{5em}}
\def\l@subparagraph{\@dottedtocline{5}{12em}{6em}}
\def\listoffigures{\@restonecolfalse\if@twocolumn\@restonecoltrue\onecolumn
 \fi\section*{List of Figures\@mkboth
 {LIST OF FIGURES}{LIST OF FIGURES}}\@starttoc{lof}\if@restonecol\twocolumn
 \fi}
\def\l@figure{\@dottedtocline{1}{1.5em}{2.3em}}
\def\listoftables{\@restonecolfalse\if@twocolumn\@restonecoltrue\onecolumn
 \fi\section*{List of Tables\@mkboth
 {LIST OF TABLES}{LIST OF TABLES}}\@starttoc{lot}\if@restonecol\twocolumn
 \fi}
\let\l@table\l@figure
%
% Redefinition to remove dotted lines from \@dottedtocline
%
\def\@dottedtocline#1#2#3#4#5{\ifnum #1>\c@tocdepth \else
  \vskip \z@ plus .2\p@
  {\leftskip #2\relax \rightskip \@tocrmarg \parfillskip -\rightskip
    \parindent #2\relax\@afterindenttrue
   \interlinepenalty\@M
   \leavevmode
   \@tempdima #3\relax \advance\leftskip \@tempdima
   \hbox{}\hskip -\leftskip
    #4\nobreak\hfill \nobreak \hbox to\@pnumwidth{\hfil
   \rm #5}\@@par}\fi}

\def\footnoterule{}%
\setcounter{footnote}{0}
\@addtoreset{footnote}{page}
\long\def\@makefntext#1{\parindent 1em\noindent
 \makebox[1em][l]{\footnotesize\rm$\m@th{\fnsymbol{footnote}}$}%
 \footnotesize\rm #1}
\def\@makefnmark{\hbox{${\fnsymbol{footnote}}\m@th$}}
\def\@thefnmark{\fnsymbol{footnote}}
\def\footnote{\@ifnextchar[{\@xfootnote}{\stepcounter{\@mpfn}%
     \begingroup\let\protect\noexpand
       \xdef\@thefnmark{\thempfn}\endgroup
     \@footnotemark\@footnotetext}}
\def\@fnsymbol#1{\ifcase#1\or \dagger\or \ddagger\or \S\or
   \|\or \P\or ^{+}\or ^{\tsty *}\or \sharp
   \or \dagger\dagger \else\@ctrerr\fi\relax}
\newcommand\ftnote[1]{\setcounter{footnote}{#1}%
   \addtocounter{footnote}{-1}\footnote}
\newcommand{\fnm}[1]{\setcounter{footnote}{#1}\footnotetext}
%
\def\center{\trivlist\topsep=0\p@\partopsep=0\p@
   \parsep=0\p@\itemsep=0\p@\centering\item[]}
%
\newenvironment{indented}{\begin{indented}}{\end{indented}}
\def\indented{\list{}{\itemsep=0\p@\labelsep=0\p@\itemindent=0\p@
   \labelwidth=0\p@\leftmargin=1.5cm\rightmargin=1.5cm
   \topsep=0\p@\partopsep=0\p@
   \parsep=0\p@\listparindent=0\p@}\rm}

\let\endindented=\endlist
%
\def\catchline{\hfill}

\def\cpyrtline{\hfill}
%
\def\maketitle{\thispagestyle{myheadings}%
   \vspace*{1.8cm}
   \begin{center}\@title\end{center}
   \vspace*{1.1cm}
   \normalsize\rm
   \begin{center}\@author\end{center}
   \begin{center}\@address\end{center}
   \@collab
   \@abstract}
%
%  Title
%
\def\title#1{\def\@title{\exhyphenpenalty=10000\hyphenpenalty=10000
    \Large\bf#1\par}}
\def\shortitle#1{\def\@shorttitle{#1}}
\let\paper=\title
%
% Authors
%
\renewcommand{\author}[1]{\def\@author{{\large #1\par}}}
%
% Affiliation
%
\newcommand{\address}[1]{\def\@address{\rm #1\par}}
\let\affil=\address
%
\newcommand{\collab}[1]{\def\@collab{\begin{center}%
   {\large\rm #1}\par
   \end{center}}}
%
% Default
%
\def\@collab{}
%
% Abstract
%
\def\abstract#1{\def\@abstract{\begin{center}
{\bf\abstractname}\end{center}%
\begin{indented}
\item[]#1\par
\end{indented}
\vspace{2cm minus1cm}}}%
%
\def\endabstract{}
%
% Command for second and subsequent paragraphs in abstract
%
\def\cabs{\\\hspace*{16\p@}}
%
\def\nosections{\vspace{30\p@ plus12\p@ minus12\p@}
    \noindent\ignorespaces}
%
\def\ack{\ifletter\bigskip\noindent\ignorespaces\else
    \section*{Acknowledgments}\fi}
\def\ackn{\ifletter\bigskip\noindent\ignorespaces\else
    \section*{Acknowledgment}\fi}
%
\newif\ifnumbysec
\def\theequation{\ifnumbysec
      \arabic{section}.\arabic{equation}\else
      \arabic{equation}\fi}
\def\eqnobysec{\numbysectrue\@addtoreset{equation}{section}}
%
\def\eqalign#1{\null\vcenter{\def\\{\cr}\openup\jot\m@th
  \ialign{\strut\hfil$\displaystyle{##{}}$&$\displaystyle{{}##}$\hfil
      \crcr#1\crcr}}\,}
%
\def\eqalignno#1{\displ@y \tabskip\z@skip
  \halign to\if@twocolumn\columnwidth\else\displaywidth\fi
   {\hfil$\@lign\displaystyle{##}$%
    \tabskip\z@skip
    &$\@lign\displaystyle{{}##}$\hfill\tabskip\@centering
    &\llap{$\@lign\hbox{\rm##}$}\tabskip\z@skip\crcr
    #1\crcr}}
%
\def\numparts{\addtocounter{equation}{1}%
     \setcounter{eqnval}{\value{equation}}%
     \setcounter{equation}{0}%
     \def\theequation{\ifnumbysec
     \arabic{section}.\arabic{eqnval}{\it\alph{equation}}%
     \else\arabic{eqnval}{\it\alph{equation}}\fi}}

\def\endnumparts{\def\theequation{\ifnumbysec
     \arabic{section}.\arabic{equation}\else
     \arabic{equation}\fi}%
     \setcounter{equation}{\value{eqnval}}}
%
\def\cases#1{%
     \left\{\,\vcenter{\def\\{\cr}\normalbaselines\openup1\jot\m@th%
     \ialign{\strut$\displaystyle{##}\hfil$&\tqs
     \rm##\hfil\crcr#1\crcr}}\right.}%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  Floats
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  \c@topnumber            : Number of floats allowed at the top of a column.
%
\setcounter{topnumber}{4}
%
%  \topfraction            : Fraction of column that can be devoted to floats.
%
\def\topfraction{1}
%
%  \c@dbltopnumber, \dbltopfraction : Same as above, but for double-column
%                          floats.
%
\setcounter{dbltopnumber}{4}
\def\dbltopfraction{1}
%
%  \c@bottomnumber, \bottomfraction : Same as above for bottom of page.
%
\setcounter{bottomnumber}{2}
\def\bottomfraction{.8}
%
%  \c@totalnumber          : Number of floats allowed in a single column,
%                          including in-text floats.
%
\setcounter{totalnumber}{5}
%
%  \textfraction         : Minimum fraction of column that must contain text.
%
\def\textfraction{0}
%
%  \floatpagefraction    : Minimum fraction of page that must be taken
%                          up by float page.
%
\def\floatpagefraction{.8}
%
%  \dblfloatpagefraction : Same as above, for double-column floats.
%
\def\dblfloatpagefraction{.8}
%
\newcounter{figure}
\def\thefigure{\@arabic\c@figure}
\def\figure{\let\@makecaption\@makeonecolcaption\@float{figure}}
\let\endfigure\end@float
%
\@namedef{figure*}{\let\@makecaption\@makewidecaption
      \@dblfloat{figure}}
\@namedef{endfigure*}{\end@dblfloat}
%
\def\@makewidecaption#1#2{\vspace{10\p@}%
     \sbox{\captionbox}{\noindent\footnotesize\rm\raggedright{\bf #1.} #2}%
     \ifdim\wd\captionbox > \indentedwidth
     \begin{indented}
     \item[]\footnotesize\rm\raggedright{\bf #1.} #2\par
     \end{indented}%
     \else
     \hbox to \hsize{\hfil\box\captionbox\hfil}\fi}
%
\def\@makeonecolcaption#1#2{\vspace{10pt}%
     \parbox{\columnwidth}{\noindent
     \footnotesize\rm\raggedright{\bf #1.} #2}\par}
%
%
% The document style must define the following.
%
%    \fps@TYPE   : The default placement specifier for floats of type TYPE.
%
\def\fps@figure{tb}
\def\fps@table{tb}
%
%    \ftype@TYPE : The type number for floats of type TYPE.
%
\def\ftype@figure{1}
\def\ftype@table{2}
%
%    \ext@TYPE   : The file extension indicating the file on which the
%                  contents list for float type TYPE is stored.  For example,
%                  \ext@figure = 'lof'.
%
\def\ext@table{aux}
\def\ext@figure{aux}
%
%    \fnum@TYPE  : A macro to generate the figure number for a caption.
%                  For example, \fnum@TYPE == Figure \thefigure.
%
\def\fnum@table{\tablename~\thetable}
\def\fnum@figure{\figurename~\thefigure}
%
%    \@makecaption{NUM}{TEXT} : A macro to make a caption, with NUM the value
%                  produced by \fnum@... and TEXT the text of the caption.
%                  It can assume it's in a \parbox of the appropriate width.
%
\newcommand{\Figure}[2]{\def\figspace{\vspace*{#1}}%
    \def\figcap{\caption{#2}}%
    \futurelet\next\@figplace}
\def\@figplace{\ifx\next[\let\next=\@figpl
                 \else\let\next=\@fignopl\fi\next}
\def\@figpl[#1]{\begin{figure}[#1]
   \figspace
   \figcap
   \end{figure}}
\def\@fignopl{\begin{figure}
   \figspace
   \figcap
   \end{figure}}
%
\newcommand{\widefigure}[2]{\def\figspace{\vspace*{#1}}%
    \def\figcap{\caption{#2}}%
    \futurelet\next\@wfigplace}
\def\@wfigplace{\ifx\next[\let\next=\@wfigpl
                 \else\let\next=\@wfignopl\fi\next}
\def\@wfigpl[#1]{\begin{figure*}[#1]
   \figspace
   \figcap
   \end{figure*}}
\def\@wfignopl{\begin{figure*}
   \figspace
   \figcap
   \end{figure*}}
%
% \@float{TYPE}[PLACEMENT] : This macro begins a float environment for a
%     single-column float of type TYPE with PLACEMENT as the placement
%     specifier.  The default value of PLACEMENT is defined by \fps@TYPE.
%     The environment is ended by \end@float.
%     E.g., \figure == \@float{figure}, \endfigure == \end@float.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  Tables
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
\newcounter{table}
%
\def\thetable{\@arabic\c@table}
\def\table{\let\@makecaption\@makeonecolcaption
    \footnotesize\rm\@float{table}}
\let\endtable\end@float
%
\@namedef{table*}{\let\@makecaption\@makewidecaption
   \footnotesize\rm
   \@dblfloat{table}}
\@namedef{endtable*}{\end@dblfloat}
%
\def\tabular{\def\@halignto{}\@tabular}
\def\endtabular{\crcr\egroup\egroup $\egroup}
\expandafter \let \csname endtabular*\endcsname = \endtabular
%
\newsavebox{\tablebox}
%
\newcommand{\Table}[2]{\begin{center}
    \lineup
    \begin{tabular}{#1}%
    \hline
    #2
    \hline
    \end{tabular}
    \end{center}}
%
\newcommand{\tabnote}[1]{\begin{indented}
     \item[]\footnotesize\rm\raggedright #1\par
     \end{indented}}
%
% Definitions for centring headings over several columns
% \centre{4}{Results for helium} will centre
% Results for helium over four columns
% \crule{4} will produce a rule centred over four columns
% to go below a centred heading
%
\newcommand{\centre}[2]{\multicolumn{#1}{c}{#2}}
\newcommand{\crule}[1]{\multispan{#1}{\hrulefill}}
%
\def\lineup{\def\0{\hbox{\phantom{\footnotesize\rm 0}}}%
    \def\m{\hbox{$\phantom{-}$}}%
    \def\-{\llap{$-$}}}
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% References
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
\newcommand{\Bibliography}[1]{\section*{References}\par\numrefs{#1}}
\newcommand{\References}[1]{\section*{References}\footnotesize\rm}
%
\def\thebibliography#1{\list
 {\hfil[\arabic{enumi}]}{\topsep=0\p@\parsep=0\p@
 \partopsep=0\p@\itemsep=0\p@
 \labelsep=5\p@\itemindent=0\p@                %-10
 \settowidth\labelwidth{\footnotesize[#1]}%
 \leftmargin\labelwidth
 \advance\leftmargin\labelsep
% \advance\leftmargin -\itemindent
 \usecounter{enumi}}%
 \def\newblock{\ }
 \sloppy\clubpenalty4000\widowpenalty4000
 \sfcode`\.=1000\footnotesize\rm\relax}
\let\endthebibliography=\endlist
%
\def\numrefs#1{\begin{thebibliography}{#1}}
\def\endnumrefs{\end{thebibliography}}
\let\endbib=\endnumrefs

\mark{{}{}}

\def\ps@headings{\let\@mkboth\markboth
 \def\@oddfoot{}%
 \def\@evenfoot{}%
 \def\@evenhead{\makebox[\mathindent][l]{\normalsize\rm \thepage}%
  \normalsize\it\rightmark\hfill}%
 \def\@oddhead{\makebox[\mathindent][r]{\hfill}{\normalsize\it\leftmark}\hfill
  \normalsize\rm\thepage}%
}%

\def\ps@myheadings{\let\@mkboth\markboth
 \def\@oddhead{\catchline}%
 \def\@oddfoot{\cpyrtline}%
 \def\@evenhead{}%
 \def\@evenfoot{}%
}


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

\def\@begintheorem#1#2{\rm \trivlist \item[\hskip \labelsep{\it #1\ #2.}]}
\def\@opargbegintheorem#1#2#3{\rm \trivlist
      \item[\hskip \labelsep{\it #1\ #2\ (#3).}]}

\let\scap=\sc
\renewcommand{\sc}{\protect\scriptsize}
\newcommand{\itsc}{\protect\scriptsize\it}
\newcommand{\bfsc}{\protect\scriptsize\bf}
\def\p@LaTeX{{L\kern-.3em\lower.1em\hbox{$^{\rm A}$}\kern-.15em%
    T\kern-.1667em\lower.7ex\hbox{E}\kern-.125emX}}
%
\newcommand{\nohyphens}{\hyphenpenalty=10000\exhyphenpenalty=10000}
\newcommand{\fl}{\hspace*{-\mathindent}}
\newcommand{\Tr}{\mathop{\rm Tr}\nolimits}
\newcommand{\tr}{\mathop{\rm tr}\nolimits}
\newcommand{\Or}{\mathop{\rm O}\nolimits}
\newcommand{\lshad}{[\![}
\newcommand{\rshad}{]\!]}
\newcommand{\case}[2]{{\textstyle\frac{#1}{#2}}}
\def\pt(#1){({\it #1\/})}
\newcommand{\dsty}{\displaystyle}
\newcommand{\tsty}{\textstyle}
\newcommand{\ssty}{\scriptstyle}
\newcommand{\sssty}{\scriptscriptstyle}
\def\lo#1{\llap{${}#1{}$}}
\def\eql{\llap{${}={}$}}
\def\lsim{\llap{${}\sim{}$}}
\def\lsimeq{\llap{${}\simeq{}$}}
\def\lequiv{\llap{${}\equiv{}$}}
%
\def\;{\protect\psemicolon}
\def\psemicolon{\relax\ifmmode\mskip\thickmuskip\else\kern .3333em\fi}
%
\newcommand{\eref}[1]{(\ref{#1})}
\newcommand{\sref}[1]{section~\ref{#1}}
\newcommand{\fref}[1]{figure~\ref{#1}}
\newcommand{\tref}[1]{table~\ref{#1}}
\newcommand{\Eref}[1]{Equation~(\ref{#1})}
\newcommand{\Sref}[1]{Section~\ref{#1}}
\newcommand{\Fref}[1]{Figure~\ref{#1}}
\newcommand{\Tref}[1]{Table~\ref{#1}}

\newcommand{\opencirc}{\raisebox{2\p@}{\;\circle{5}}}
\newcommand{\opensqr}{\mbox{$\Box$}}
\newcommand{\fullcirc}{\raisebox{-2\p@}{\Large$\bullet$}}
\newcommand{\fullsqr}{\mbox{\vrule height6pt width6pt}}
\newcommand{\dotted}
                 {\mbox{${\mathinner{\cdotp\cdotp\cdotp\cdotp\cdotp\cdotp}}$}}
\newcommand{\dashed}{\mbox{-\; -\; -\; -}}
\newcommand{\broken}{\mbox{-- -- --}}
\newcommand{\longbroken}{\mbox{--- --- ---}}
\newcommand{\chain}{\mbox{--- $\cdot$ ---}}
\newcommand{\dashddot}{\mbox{--- $\cdot$ $\cdot$ ---}}
\newcommand{\full}{\mbox{------}}
%
\newcommand{\etal}{{\it et al\/}\ }
\newcommand{\nonum}{\item[]}
%
% abbreviations for IOPP journals
%
\newcommand{\CQG}{{\em Class. Quantum Grav.} }
\newcommand{\HPP}{{\em High Perform. Polym.} }              % added 4/5/93
\newcommand{\IP}{{\em Inverse Problems\/} }
\newcommand{\JHM}{{\em J. Hard Mater.} }                    % added 4/5/93
\newcommand{\JPA}{{\em J. Phys. A: Math. Gen.} }
\newcommand{\JPB}{{\em J. Phys. B: At. Mol. Phys.} }      %1968-87
\newcommand{\jpb}{{\em J. Phys. B: At. Mol. Opt. Phys.} } %1988 and onwards
\newcommand{\JPC}{{\em J. Phys. C: Solid State Phys.} }   %1968--1988
\newcommand{\JPCM}{{\em J. Phys.: Condens. Matter\/} }    %1989 and onwards
\newcommand{\JPD}{{\em J. Phys. D: Appl. Phys.} }
\newcommand{\JPE}{{\em J. Phys. E: Sci. Instrum.} }
\newcommand{\JPF}{{\em J. Phys. F: Met. Phys.} }
\newcommand{\JPG}{{\em J. Phys. G: Nucl. Phys.} }         %1975--1988
\newcommand{\jpg}{{\em J. Phys. G: Nucl. Part. Phys.} }   %1989 and onwards
\newcommand{\MSMSE}{{\em Modelling Simulation Mater. Sci. Eng.} }
\newcommand{\MST}{{\em Meas. Sci. Technol.} }              %1990 and onwards
\newcommand{\NET}{{\em Network\/} }
\newcommand{\NL}{{\em Nonlinearity\/} }
\newcommand{\NT}{{\em Nanotechnology} }
\newcommand{\PAO}{{\em Pure Appl. Optics\/} }
\newcommand{\PM}{{\em Physiol. Meas.} }                        % added 4/5/93
\newcommand{\PMB}{{\em Phys. Med. Biol.} }
\newcommand{\PPCF}{{\em Plasma Phys. Control. Fusion\/} }      % added 4/5/93
\newcommand{\PSST}{{\em Plasma Sources Sci. Technol.} }
\newcommand{\QO}{{\em Quantum Opt.} }
\newcommand{\RPP}{{\em Rep. Prog. Phys.} }
\newcommand{\SLC}{{\em Sov. Lightwave Commun.} }               % added 4/5/93
\newcommand{\SST}{{\em Semicond. Sci. Technol.} }
\newcommand{\SUST}{{\em Supercond. Sci. Technol.} }
\newcommand{\WRM}{{\em Waves Random Media\/} }
%
% Other commonly quoted journals
%
\newcommand{\AC}{{\em Acta Crystallogr.} }
\newcommand{\AM}{{\em Acta Metall.} }
\newcommand{\AP}{{\em Ann. Phys., Lpz.} }
\newcommand{\APNY}{{\em Ann. Phys., NY\/} }
\newcommand{\APP}{{\em Ann. Phys., Paris\/} }
\newcommand{\CJP}{{\em Can. J. Phys.} }
\newcommand{\JAP}{{\em J. Appl. Phys.} }
\newcommand{\JCP}{{\em J. Chem. Phys.} }
\newcommand{\JJAP}{{\em Japan. J. Appl. Phys.} }
\newcommand{\JP}{{\em J. Physique\/} }
\newcommand{\JPhCh}{{\em J. Phys. Chem.} }
\newcommand{\JMMM}{{\em J. Magn. Magn. Mater.} }
\newcommand{\JMP}{{\em J. Math. Phys.} }
\newcommand{\JOSA}{{\em J. Opt. Soc. Am.} }
\newcommand{\JPSJ}{{\em J. Phys. Soc. Japan\/} }
\newcommand{\JQSRT}{{\em J. Quant. Spectrosc. Radiat. Transfer\/} }
\newcommand{\NC}{{\em Nuovo Cimento\/} }
\newcommand{\NIM}{{\em Nucl. Instrum. Methods\/} }
\newcommand{\NP}{{\em Nucl. Phys.} }
\newcommand{\PL}{{\em Phys. Lett.} }
\newcommand{\PR}{{\em Phys. Rev.} }
\newcommand{\PRL}{{\em Phys. Rev. Lett.} }
\newcommand{\PRS}{{\em Proc. R. Soc.} }
\newcommand{\PS}{{\em Phys. Scr.} }
\newcommand{\PSS}{{\em Phys. Status Solidi\/} }
\newcommand{\PTRS}{{\em Phil. Trans. R. Soc.} }
\newcommand{\RMP}{{\em Rev. Mod. Phys.} }
\newcommand{\RSI}{{\em Rev. Sci. Instrum.} }
\newcommand{\SSC}{{\em Solid State Commun.} }
\newcommand{\ZP}{{\em Z. Phys.} }
%

\def\ap#1#2#3 {Ann. Phys. (NY) {\bf#1} (19#2) #3}
\def\apj#1#2#3 {Astrophys. J. {\bf#1} (19#2) #3}
\def\apjl#1#2#3 {Astrophys. J. Lett. {\bf#1} (19#2) #3}
\def\app#1#2#3 {Acta. Phys. Pol. {\bf#1} (19#2) #3}
\def\ar#1#2#3 {Ann. Rev. Nucl. Part. Sci. {\bf#1} (19#2) #3}
\def\cpc#1#2#3 {Computer Phys. Comm. {\bf#1} (19#2) #3}
\def\err#1#2#3 {{\it Erratum} {\bf#1} (19#2) #3}
\def\ib#1#2#3 {{\it ibid.} {\bf#1} (19#2) #3}
\def\jmp#1#2#3 {J. Math. Phys. {\bf#1} (19#2) #3}
\def\ijmp#1#2#3 {Int. J. Mod. Phys. {\bf#1} (19#2) #3}
\def\jetp#1#2#3 {JETP Lett. {\bf#1} (19#2) #3}
\def\jpg#1#2#3 {J. Phys. G. {\bf#1} (19#2) #3}
\def\mpl#1#2#3 {Mod. Phys. Lett. {\bf#1} (19#2) #3}
\def\nat#1#2#3 {Nature (London) {\bf#1} (19#2) #3}
\def\nc#1#2#3 {Nuovo Cim. {\bf#1} (19#2) #3}
\def\nim#1#2#3 {Nucl. Instr. Meth. {\bf#1} (19#2) #3}
\def\np#1#2#3 {Nucl. Phys. {\bf#1} (19#2) #3}
\def\pcps#1#2#3 {Proc. Cam. Phil. Soc. {\bf#1} (#2) #3}
\def\pl#1#2#3 {Phys. Lett. {\bf#1} (19#2) #3}
\def\prep#1#2#3 {Phys. Rep. {\bf#1} (19#2) #3}
\def\prev#1#2#3 {Phys. Rev. {\bf#1} (19#2) #3}
\def\prl#1#2#3 {Phys. Rev. Lett. {\bf#1} (19#2) #3}
\def\prs#1#2#3 {Proc. Roy. Soc. {\bf#1} (19#2) #3}
\def\ptp#1#2#3 {Prog. Th. Phys. {\bf#1} (19#2) #3}
\def\ps#1#2#3 {Physica Scripta {\bf#1} (19#2) #3}
\def\rmp#1#2#3 {Rev. Mod. Phys. {\bf#1} (19#2) #3}
\def\rpp#1#2#3 {Rep. Prog. Phys. {\bf#1} (19#2) #3}
\def\sjnp#1#2#3 {Sov. J. Nucl. Phys. {\bf#1} (19#2) #3}
\def\spj#1#2#3 {Sov. Phys. JEPT {\bf#1} (19#2) #3}
\def\spu#1#2#3 {Sov. Phys.-Usp. {\bf#1} (19#2) #3}
\def\zp#1#2#3 {Zeit. Phys. {\bf#1} (19#2) #3}
%
\ps@headings \pagenumbering{arabic} \onecolumn


