Physics goals and parameters of photon colliders

Abstract: Linear colliders offer a unique possibility to study gamma gamma and gamma electron interactions at the energies 0.1--2 TeV. This option is now included in design reports of NLC, JLC and TESLA/SBLC. This paper includes: status of photon colliders, new possibilities in study of Higgs boson, ways to achieve high luminosities.

\section{Introduction}

\hspace*{\parindent}\ It is very likely that linear colliders with the c.m.s energies of
0.2--2 TeV will be built sometime, may be in about ten years from
now~\cite{LOEW}. Besides \mbox{e$^+$e$^-$} collisions, linear colliders give an
unique possibility to study \mbox{$\gamma\gamma$}\ and \mbox{$\gamma$e}\ interactions at energies
and luminosities comparable to those in \mbox{e$^+$e$^-$}\
collisions~\cite{GKST81}-\cite{GKST84}.

The basic scheme of a photon collider is shown in
fig.~\ref{ris1}.
\begin{figure}[!hbt]
\hspace*{1cm}\begin{minipage}[b]{0.45\linewidth}
\centering
\vspace*{-0.0cm}
\hspace*{.0cm} \epsfig{file=ris1.eps,height=6.5cm,angle=-90}
\vspace*{0.3cm}
\caption{Scheme of  \mbox{$\gamma\gamma$}; \mbox{$\gamma$e}\ collider.}
\label{ris1}
\end{minipage}%
\hspace*{0.3cm} \begin{minipage}[b]{0.45\linewidth}
\centering
\vspace*{-0.0cm}
\hspace*{0.5cm} \epsfig{file=fig9.eps,width=0.8in,angle=-90}
\vspace*{1.5cm}
\caption{Crab-crossing scheme}
\label{fig9}
\end{minipage}
\end{figure}
Two electron beams
after the final focus system are traveling toward the interaction
point (IP).  At a distance of about 0.1--1 cm upstream from the IP, at
the conversion point (C), the laser beam is focused and Compton
backscattered by electrons, resulting in the high energy beam of
photons.  With reasonable laser parameters one can ``convert'' most of
the electrons into high energy photons. The photon beam follows the
original electron direction of motion with a small angular spread of
the order $1/\gamma $, arriving at the IP  tightly focused, where it
collides with a  similar opposing high energy photon beam or with an
electron beam. The photon spot size at the IP may be almost equal to
that of electrons at IP, and, therefore, the luminosity of \mbox{$\gamma\gamma$}, \mbox{$\gamma$e}\
collisions will be of the same order of magnitude as the ``geometric''
luminosity of the basic $ee$ beams.  The detailed discussion of photon
colliders  can be found in refs~\cite{GKST83}-\cite{TEL95}.
and in the Berkeley Workshop Proceedings~\cite{BERK}.

  As of today, this option is included into the Conceptual Design Reports of
 NLC~\cite{NLC}, TESLA--SBLC~\cite{TESLA}, and JLC~\cite{JLC} linear
colliders. All these  projects foresee the second
interaction regions for \mbox{$\gamma\gamma$}, \mbox{$\gamma$e}\ collisions.

   This paper covers three topics:
\begin{itemize}
\vspace*{-0.3cm}\item  physics at photon
colliders,  paticularly some important remarks on new possibilities
to study the Higgs boson in \mbox{$\gamma\gamma$}\ collisions;
\vspace*{-0.3cm}\item  parameters of photon colliders in `Zero design'
projects of linear colliders;
\vspace{-0.3cm}\item  ways to achieve very high luminosities.
\end{itemize}

\section{Physics}
\subsection{General remarks}
  The physics in high energy \mbox{$\gamma\gamma$}, \mbox{$\gamma$e}\ colliders is very rich. The
total number of papers devoted to this subject exceeds one
thousand. Recent reviews of physics at photon colliders can be found,
for instance, in TESLA/SBLC Conceptual Design Report~\cite{TESLA}
(with many references therein) and in the J.Jikia's talk at
``Photon 97''~\cite{JIKIA97}.  In this paper, besides some general words
on the physics in \mbox{$\gamma\gamma$}, \mbox{$\gamma$e}\ collisions, I would like to discuss some new
important aspects of the Higgs study in \mbox{$\gamma\gamma$}\ collisions.

   For scientific (and political) motivation of building photon colliders the
   following short list of  physics goals can be presented:

\vspace*{-1mm}

\begin{enumerate}

\vspace*{-1mm}
 \item Some phenomena can be better studied at photon colliders
   than with pp or \mbox{e$^+$e$^-$}\ collisions, for example, the
  measurement of the two-photon decay width of the Higgs boson.  Due the
  loop diagrams, all massive  (even ultra-heavy) charged particles
  contribute to this width if their mass is originated by the Higgs
  mechanism.  Some Higgs decay modes and its mass can be
  measured at \mbox{$\gamma\gamma$}\ colliders better than in \mbox{e$^+$e$^-$}\ or pp collisions
  (see the next section).

\vspace*{-1mm}
\item Cross sections for production of charged scalar,
  lepton and top pairs in \mbox{$\gamma\gamma$}\ collisions are larger than those in
  \mbox{e$^+$e$^-$}\ collisions approximately by a factor of 5 (see
  fig.~ \ref{fig16}); for WW production, this factor is even larger,
  about 10--20.

\vspace*{-1mm}
\item In \mbox{$\gamma$e}\ collisions, charged supersymmetric
  particles with masses higher than in \mbox{e$^+$e$^-$}\ collisions can be
  produced (a heavy charged particle plus a light neutral), \mbox{$\gamma\gamma$}\
  collisions also provide higher accessible masses for particles which
  are produced as a single resonance in \mbox{$\gamma\gamma$}\ collisions (such as the Higgs).
\end{enumerate}

These examples together with the fact that the luminosity in \mbox{$\gamma\gamma$}\
 collisions is potentially higher than that in \mbox{e$^+$e$^-$}\ collisions (see
 sect. 4) are very strong arguments in favor of photon colliders.

\begin{figure}[thb]
\centering
\vspace*{-1.4cm}
\epsfig{file=fig16.eps,width=6.5in}
\vspace*{-1.2cm}
\caption{ Comparison of  cross sections for charged pair production
in \mbox{e$^+$e$^-$}\ and \mbox{$\gamma\gamma$}\ collisions. The cross section  $\sigma
=(\pi\alpha^2/M^2)f(x)$, P=S (scalars), F (fermions),
W (W-bosons); M is particle mass, $x=W_{p\bar{p}}^2/4M^2$. The functions
$f(x)$ are shown.}

\vspace{-0.3cm}

\label{fig16}
\end{figure}

  A typical distribution of \mbox{$\gamma\gamma$}\ or \mbox{$\gamma$e}\  luminosity on the  invariant mass
has  high energy peaks at the maximum masses with widths
$\Delta W_{\mbox{$\gamma\gamma$}} / W_{\mbox{$\gamma\gamma$}} \approx 0.15$,  $\Delta
 W_{\mbox{$\gamma$e}}/W_{\mbox{$\gamma$e}} \approx 0.05$\cite{TEL95,TESLA}. Below
the high energy peak, there is usually a flat part of the luminosity
distribution with 1--10 times larger total luminosity (depending
on details of the collision scheme).

Statistics equal to that in \mbox{e$^+$e$^-$}\ collisons, can be achieved in \mbox{$\gamma\gamma$}\
collisions with much smaller (at least by a factor of five) luminosity. Of
course, in any case data in \mbox{$\gamma\gamma$}\ and \mbox{e$^+$e$^-$}\ collisions are
complimentary to each other because  the coresponding diagrams are different).

\subsection{Higgs in {\boldmath$\gamma\gamma$} collisions}
Search and study of the Higgs boson, the key particle of the Standard
model, is one of the primary goals of linear colliders. The LEP-2 will
put the lower limit on its mass ot about 95 GeV.  Indirectly, from
radiative corrections, it follows that the Higgs mass (if the Standard
Model is correct) is $140^{+150}_{-80}$ GeV~\cite{BLONDEL}. The MSSM
predicts a mass of the lightest neutral Higgs less of than 130
GeV. So, the mass of the Higgs most probably lies in the region of
100$<M_H<$300 GeV.  Some parameters of the SM Higgs (total width, \mbox{$\gamma\gamma$}\
width, and main branching ratios) are presented in
fig. \ref{hbranch}$\;$~\cite{BAIL94,DJOUADI}.

  Higgs production in \mbox{$\gamma\gamma$}\ collisions has been considered in many
papers. Unfortunately, different authors have used their own definitions
of the \mbox{$\gamma\gamma$}\ luminosity which often led  to underestimation of the
expected number of events for a given integrated \mbox{$\gamma\gamma$}\ luminosity. In the
present paper, we will directly compare  the Higgs production cross
sections in \mbox{$\gamma\gamma$}\ and \mbox{e$^+$e$^-$}\ collisions. I would like also to pay
attention to some additional possibilities for Higgs studies at photon
colliders connected with a very sharp edge in the \mbox{$\gamma\gamma$}\ luminosity spectra.

\begin{figure}[thb]
\hspace*{1cm}\begin{minipage}[b]{0.45\linewidth}
\centering
\hspace*{-2.0cm}
\scalebox{0.9}[1.]{
\includegraphics[width=3.5in, angle=0, bb=90 140 490 425,clip] {boudjema1.ps}}
\end{minipage}%
\hspace*{0.5cm} \vspace*{0.3cm} \begin{minipage}[b]{0.45\linewidth}
\centering
\scalebox{1.}[1.]{
\includegraphics[width=2.43in, angle=0, bb=308 293 457 441,clip] {jik1.ps}
}
\end{minipage}

\vspace{0.3cm}

\caption{Some of the SM Higgs branching ratios as a function
of its mass~\cite{BAIL94,DJOUADI}. On the left figure~\cite{BAIL94},
the total width and the $\gamma\gamma$ width (thick lines) are also shown.}
\label{hbranch}
\end{figure}

 The total Higgs width at masses below 400 GeV is much smaller than
the characteristic width of \mbox{$\gamma\gamma$}\ luminosity spectra (FWHM $\sim 10-15 \%$
), therefore the production rate is proportional to $d\mbox{$L_{\gamma\gamma}$}/W_{\mbox{$\gamma\gamma$}}$:
\begin{equation}
\dot{N}_{\mbox{$\gamma\gamma$}\ \to H}
=\mbox{$L_{\gamma\gamma}$} \times \frac{d\mbox{$L_{\gamma\gamma}$} M_H}{d\mbox{$W_{\gamma\gamma}$} \mbox{$L_{\gamma\gamma}$}}\frac{4\pi^2\Gamma_{\mbox{$\gamma\gamma$}}
(1+\lambda_1 \lambda_2)}{M^3_H}\equiv \mbox{$L_{\gamma\gamma}$} \times \sigma^{eff}.
\end{equation}
where $\lambda_i$ is the photon helicity.
  Below we  assume that the Higgs search and study is done utilizing the high
energy peak of the \mbox{$\gamma\gamma$}\ luminosity energy spectrum.

The effective cross section for $(d\mbox{$L_{\gamma\gamma}$}/d\mbox{$W_{\gamma\gamma}$}) (M_H/\mbox{$L_{\gamma\gamma}$})=7$ and
$1+\lambda_1 \lambda_2=2$ is presented in fig.~\ref{cross}.
Note that here \mbox{$L_{\gamma\gamma}$}\
is defined as the \mbox{$\gamma\gamma$}\ luminosity at the high energy luminosity peak
($z=\mbox{$W_{\gamma\gamma}$}/2E_e>0.65$ for $x=4.8$), thus ignoring the lower energy part
of the luminosity spectrum,
which is much less valuable for experiment.\footnote{It
is also possible to search for the Higgs at a constant
collider energy utilizing a flat part of the luminosity spectrum
instead of energy scanning. However, in this case $d\mbox{$L_{\gamma\gamma}$}/d\mbox{$W_{\gamma\gamma}$}\ $ is
lower (no peak and a wider energy range) and we have much larger
backgrounds ($\gamma g \to b\bar{b}$, etc), worse polarization degree
(partially due to contribution of multiple Compton scattering and
unpolarized beamstrahlung photons) and, consequently, worse
suppression of backgrounds; therefore, the required total integrated
luminosity will be only larger than in the method of scanning (this
statement should be checked more carefully). If the Higgs mass is
already known, then, obviously, one has to work at $\mbox{$W_{\gamma\gamma}$}_{max}\sim
M_H$.}
 This luminosity is approximately equal to
$0.25k^2\mbox{$L_{ee}$}(geom)$.  For comparison, in the same figure the cross
sections of the Higgs production in \mbox{e$^+$e$^-$}\ collisions is shown. We see
that for $M_H=$ 100--250 GeV the effective cross section in \mbox{$\gamma\gamma$}\
collisions is larger than that in \mbox{e$^+$e$^-$}\ collisions by a factor of
about 5--30. This interesting fact has never been emphasized.

\begin{figure}[!htb]

\vspace{-1cm}

\hspace*{0cm}\begin{minipage}[b]{0.45\linewidth}
\centering
\vspace*{-0.cm}
\hspace*{-1.cm} \epsfig{file=cross.eps,width=8.cm,angle=0}

\vspace{-0.4cm}

\caption{Cross sections for the Standard model Higgs in \mbox{$\gamma\gamma$}\ and
 \mbox{e$^+$e$^-$}\ collisions.}
\label{cross}
\end{minipage}%
\hspace*{1.1cm} \begin{minipage}[b]{0.45\linewidth}
\centering
\vspace*{0.cm}
\hspace*{-1.4cm} \epsfig{file=edge.eps,width=8cm,angle=0}

\vspace{-0.4cm}

\caption{ Shape of the \mbox{$\gamma\gamma$}\ luminosity spectrum near the high energy edge.
Here $z=W_{\mbox{$\gamma\gamma$}}/2E_0.$}
\label{edge}
\end{minipage}
\end{figure}

How do we study the Higgs in \mbox{$\gamma\gamma$}\ collisions? One obvious way is
searching for a peak in the invariant mass distribution  measured by
the detector. This is the only method  utilizing the broad
luminosity spectrum. The mass resolution in this method can be somewhat better
than the width of the high energy luminosity peak.

Another method is the energy scanning where  we can use some important
features of the \mbox{$\gamma\gamma$}\ luminosity distribution: quite narrow width of
high energy peak and a {\it very sharp edge of luminosity distribution}
which is much narrower than the width of the luminosity peak and the detector
resolution, see fig. \ref{edge}. For $x=4.8$ the differential
luminosity ($dL/dz$) reaches  half of its maximum  at $\Delta
z=z_{max}-z=0.8\%$ (at $\rho = b/\gamma a_e=1$).
During the scanning, we will observe a sharp increase in the visible cross
section when the maximum energy of the \mbox{$\gamma\gamma$}\ collider reaches $M_H$ (if
the Higgs is a very narrow resonance). Observation of a step at the same
energy for different decay modes significantly increases confidence of
results, especially for modes with small branchings or with
difficulties in reconstruction. For example, detection of the
$H\to\tau\tau$ decay is a difficult task due to undetected neutrinos. The most
reliable way to see this channel is to observe a step in visible
cross section for events with the following selection criteria: two collinear
low-multiplicity jets with unbalanced transverse momentum and with the energy
of one jet not far from the maximum photon energy.

 The total number of events for 10 fb$^{-1}$ of integrated \mbox{$\gamma\gamma$}\
luminosity (as it was defined above) for the case when the peak of the
luminosity spectrum sits at the Higgs boson mass is 13000 (18000) for
$M_H$ = 130 (150) GeV respectively (follows from cross section given
in fig. \ref{cross}). That is a lot! Due to energy scanning, the
statistics will be somewhat smaller, about $N \approx
0.5N_{peak}$. The number of events for various decay modes (without
any cuts) is given in table \ref{tabl2}.

\begin{table}[htb]
\caption{The number of  Higgs events in various decay modes (no
cuts) for $M_{H}$= 130 and 150 GeV and 10 fb$^{-1}$ of integrated \mbox{$\gamma\gamma$}\
luminosity.}
\begin{center}
\begin{tabular}{l|cccccccc} \hline
       Mode & $b\bar{b}$ &WW$^*$&ZZ$^*$&$\tau\tau$&$c\bar{c}$&gg&
\mbox{$\gamma\gamma$}\ &$\gamma$ Z\\
\hline\hline
Br $M_H(130)$ &0.52&0.29&0.037&0.055&0.025&0.063&0.0022&0.0019\\ \hline

events&3400&1900&240&350&160&410&15&12 \\ \hline\hline

Br $M_H(150)$ &0.18&0.67&0.083&0.019&0.008&0.03&0.0014&0.002\\ \hline
events&1600&6000&750&171&74&260&13&22 \\ \hline\hline
\end{tabular}
\end{center}
\label{tabl2}
\end{table}

  The Higgs production at photon colliders and various backgrounds have been
  studied in many papers for various decay modes:
  $b\bar{b}\;$
\cite{BARKLOW}-\cite{OHGAKI97};
 ZZ~\cite{BBC1,GUNION93,JIKIAZZ};
  WW~\cite{BBC1,GUNION93,MORRIS,GINZ97}.

Below $M_H < 150$ GeV, the SM Higgs decays mainly into $b\bar{b}$
pairs. The main backgrounds here are the QED processes $\mbox{$\gamma\gamma$}\ \to
b\bar{b},\; c\bar{c}$, which can be suppressed using  vertex detectors and
polarized photon beams ($\sigma_{\gamma\gamma \to q\bar{q}} \propto
1-\lambda_1\lambda_2$, while $\sigma_{\gamma\gamma\to H} \propto
1+\lambda_1\lambda_2$). The remaining background is small.

  For $M_H>190$ GeV, the Higgs can be observed in the best way in the ZZ
  decay mode. The main background here is the process $\gamma\gamma \to
  WW$, which can be suppressed by requiring that at least one of the Z's  be
  detected in $l^+l^-$ decay modes. In this channel, the SM Higgs can
  be observed in the range $M_H \sim$ 120 -- 350 GeV. At higher masses,
  the Higgs signal will be much smaller than irreducible background
  $\gamma\gamma \to ZZ$~\cite{JIKIAZZ}.

In the region $M_H > 140$ GeV, the SM Higgs decays mainly
into $WW$ (or $WW^*$) pairs. The main problem here is the background from
$\gamma\gamma \to WW$  with large cross section. Fortunately, this
background is very small for $M_H <2M_W$~\cite{GUNION93,MORRIS,GINZ97},
so that  Higgs with $100 <M_H < 160$ GeV
can be detected in this channel practically without backgrounds.

 In the region $M_H > 160$ GeV, the interference between the Higgs
signal and the background becomes important~\cite{MORRIS,GINZ97}.  The cross
section of $\mbox{$\gamma\gamma$}\ \to $ WW for various Higgs masses is shown in
fig.~\ref{WW} (the left figure is taken from ref.\cite{MORRIS}, the right one
from ref.\cite{GINZ97}) These cross sections correspond to monochromatic photon
beams. In a real situation the photon beam energy spread and the
detector resolution are much larger than the Higgs width, so that the
resonance structures is much broader but nevertheless observable for
not too large $M_H$. For $M_H>2M_W$, there are four additional
constraints in the 4-jet decay mode: two come from the zero transverse
momentum of the system, and two come from the jet-jet mass equal to
$M_W$. Using of these constraints should improve the WW invariant mass
resolution.

Much better sensitivity to narrow structures can be obtained utilizing
the sharp edge of the \mbox{$\gamma\gamma$}\ luminosity spectrum. Using all these
methods, one can see the Higgs in $\mbox{$\gamma\gamma$}\ \to WW$ decay mode from
$M_H\sim 100$ GeV up to $M_H \sim$ 200--250 GeV or even the higher.

\begin{figure}[thb]
\hspace*{1cm}\begin{minipage}[b]{0.45\linewidth}
\centering
\hspace*{-1.6cm}
\scalebox{1.05}[1.05]{
\includegraphics[width=7cm, angle=0, bb=104 451 471 772,clip] {morris.eps}}
\end{minipage}%
\hspace*{-0.cm} \vspace*{-0.cm} \begin{minipage}[b]{0.45\linewidth}
\centering
\vspace*{-0.3cm}
\includegraphics[width=6cm, angle=0, bb=104 229 504 657,clip] {gz1.ps}
\end{minipage}

\vspace{0.2cm}

\caption{The $\gamma\gamma\to WW$ cross section for various Higgs masses,
see comments in the text}
\label{WW}
\end{figure}

Detection of the Higgs in $\tau\tau,cc,gg,\gamma\gamma,\gamma Z$ decays are
also possible at photon colliders. In \mbox{$\gamma\gamma$}\ and $\gamma Z$ modes, the
background is rather small~\cite{JIKIAgg,JIKIAgZ}. Conclusions on
perspectives of Higgs detection  in these modes and required
luminosities can be made only after a detail simulation.

  Photon colliders also allow to measure the Higgs mass with a high
  precision. The very sharp edge of the luminosity spectrum and high
  statistics make it possible to obtain a much better  accuracy  than
  in \mbox{e$^+$e$^-$}\ or pp collisions. In this measurement, it is important to
  remember that nonlinear effects in the conversion region can shift
  the maximum energy of Compton photons and change the shape of the
  luminosity spectrum~\cite{TEL90,TEL95}. To reduce this effect, the
  parameter $\xi^2$ characterizing nonlinearity should be kept small
  enough. The shape of the luminosity spectrum near $z_{max}$ should
  be measured very precisely using the processes $\mbox{$\gamma\gamma$}\ \to \mu\mu\;
  (\mbox{e$^+$e$^-$})$.

As we have seen, the Higgs can be very successfully studied in \mbox{$\gamma\gamma$}\
collisions for $M_H$ = 100--350 GeV,  may be even better than
in \mbox{e$^+$e$^-$}\ collisions.

\section{Luminosity}
\subsection{Current projects}

 Due to the absence of beamstrahlung, beams in \mbox{$\gamma\gamma$}\ collisions can
have much smaller horizontal beam size than  in \mbox{e$^+$e$^-$}\ collisions,
thus the beta functions at the interaction point can be made as
small as possible (some restrictions are imposed by the Oide effect
connected with chromatic aberrations due to synchrotron radiation in
the final quads).  However, even after optimization of the final focusing
system the attainable \mbox{$\gamma\gamma$}\ luminosity in curent LC projects is
determined by the attainable ``geometric'' ee--luminosity.

 The results of simulations for different projects are the
   following. For the ``nominal'' beam parameters (the same as in
   \mbox{e$^+$e$^-$}\ collisions) and the optimum final focus system the
   luminosity $\mbox{$L_{\gamma\gamma}$}(z>0.65) \sim (0.8/1.2/0.7)\times10^{33}$ \mbox{cm$^{-2}$s$^{-1}$}\ for
   NLC/TESLA/SBLC, or by about factor 5 smaller than
   \mbox{e$^+$e$^-$}\ luminosity.

 Obviously, this is not a fundamental limit. Below we will discuss what
 the limit really is and how to approach  it.

\subsection{Ultimate luminosity}

  The only collision effect restricting \mbox{$\gamma\gamma$}\ luminosity at photon
  colliders is the coherent pair creation which leads to the
  conversion of a high energy photon into \mbox{e$^+$e$^-$}\ pair in the field of
  the opposing electron beam~\cite{CHEN,TEL90,TEL95}. There are
  three ways to avoid this effect: a) use flat beams; b) deflect the
  electron beam after conversion at a sufficiently large distance from
  the IP; c) under certain conditions (low beam energy, long bunches)
  the beam field at the IP is below the critical one due to the
  repulsion of electron beams~\cite{TELSH}.  The problem of ultimate
  luminosities for different beam parameters and energies was analyzed
  recently in ref.\cite{TSB2} analytically and by simulation. The  resume
  is  following.

  The maximum luminosity is attained when the conversion point is
  situated as close as possible to the IP, at $b \approx
  3\sigma_z+0.04E[\mbox{TeV}]$ cm (here the second term is equal to the
  minimum length of the conversion region). In this case, the vertical
  radius of the photon beam at the IP is also minimal: $a_{\gamma}\sim
  b/\gamma$ (assuming that the vertical size of the electron beam is
  even smaller). The optimal horizontal beam size ($\sigma_x$) depends
  on the beam energy, the number of particles in a bunch and the bunch
  length.  The dependence of the \mbox{$\gamma\gamma$}\ luminosity on $\sigma_x$ for
  various energies and numbers of particles per bunch is shown in
  fig. \ref{sb3}. The bunch length is fixed at 0.2 mm. The collision
  rate is calculated from the total beam power, which is equal to
  15E[TeV] MW (close to that in current projects). From the
  fig. \ref{sb3} we see that at low energies and small numbers of
  particles, the luminosity curves follow their natural behavior $L
  \propto 1/\sigma_x$, while at high energies and large numbers of
  particles per bunch the curves make a zigzag which is explained by
  $\gamma \to$ \mbox{e$^+$e$^-$}\ conversion in the field of the opposing beam.

 \begin{figure}[!hb]
\centering
\vspace{3mm}
\epsfig{file=sb3.eps,width=10cm,bb= 32 50 550 540 }
\caption{Dependence of the \mbox{$\gamma\gamma$}\ luminosity on the horizontal beam size
for $\sigma_z = 0.2$ \mbox{mm}, see comments in the text.}
\label{sb3}
\end{figure}

   What is so remarkable about these results? First of all, the maximum
   attainable luminosities are huge, above $10^{35}$ \mbox{cm$^{-2}$s$^{-1}$}. At low
   energies, there is no coherent pair creation, even for a  very small
   $\sigma_x$'s, when the field in the beam is much higher than the
   critical field $B_{cr}=\alpha e/\gamma r^2_e$. This is explained by
   the fact that during the collision the beams are repulsing each other
   so that the field on the beam axis (which affects  high energy
   photons) is below the critical field. It means that the
   $\gamma\gamma$ luminosity is simply proportional to the geometric
   electron-electron luminosity (approximately $\mbox{$L_{\gamma\gamma}$}(z>0.65)\sim 0.1
   L_{ee}$) for $\sigma_x,\sigma_y < b/\gamma \sim 3\sigma_z/\gamma +
   0.2$ nm.  For  energies 2E$ <$ 2 TeV, which are in reach of the next
   generation of linear colliders, the luminosity limit is much higher
   than it is usually required  and much higher than in
   \mbox{e$^+$e$^-$}\ collisions (especially at low energies).

  One of the main problems in obtaining such high \mbox{$\gamma\gamma$}\ luminosity is
generation of electron beams with small emittances in both horizontal
and vertical directions. In the following sections we will discuss
various approaches to solving this problems.

\subsection{Ways to  higher luminosities}

There are several possibilities for increasing  luminosity.

1) Reduction of the horizontal emittance by optimizing the damping
rings. For example, at the TESLA, a decrease of \mbox{$\epsilon_{nx}$}\ by a factor of
3.5 leads to an increase in \mbox{$L_{\gamma\gamma}$}\ up to $3\times 10^{33}$
\mbox{cm$^{-2}$s$^{-1}$}. However, it seems that this way is quite difficult.

2) One can use low emittance RF-photoguns instead of damping
rings. Unfortunately, even with the best photoguns the luminosity will
be somewhat lower than that with damping rings. However, there is one
possible solution. The normalized emittance in photoguns is
approximately proportional to the number of particles in the electron
bunch. It seems possible to merge (using some difference in energies)
many ($N_g \sim 5-10$) low current beams with low emittances to one
high current beam with the same transverse emittance. This gives us a
gain in luminosity more than by a factor of $N_g$ in comparison with a
single photogun (``more'' because the lower emittance allows smaller
beta functions due to the Oide effect).  Estimations show that with
ten photoguns, one can achieve \mbox{$L_{\gamma\gamma}$}\ $\sim L_{\mbox{e$^+$e$^-$}}$ in all
considered projects.  There is only one ``small'' problem: such
RF-guns with polarized electrons do not exist yet, though there are no
visible fundamental problems~\cite{CLENDENIN}. Also there is the
``thermal'' limit on emittance in photoguns, which puts some
restrictions on this method. In the AsGa photoguns with polarized
electrons, this limit should be much smaller than for metal
photocathodes~\cite{CLENDENIN}.

3) If one has flat electron beams (\mbox{$\epsilon_{nx}$} $>>$ \mbox{$\epsilon_{ny}$}) from some injector
(a damping ring or, better, an RF-gun), then it is possible, in principle, to
make exchange of the horizonlal and longitudinal emittances using
an off-axis RF cavity (some head-tail kick arizing here can be removed by
additional RF cavities rotated by 90$^0$. This method was never tested
or even discussed. I do not know why. This method is especially
promising for beams from RF-guns, where very small energy spread
allows both horizontal and longitudinal beam compression.
This method as well as the previous one should be  studied more carefully.

4) For a considerable step up in luminosity,
beams with much lower emittances are required. That needs
development of new approaches, such as {\it laser cooling}~\cite{TSB1}
(see next subsection).  Potentially, this method allows to attain a
geometric luminosity by two orders higher than that achievable by
other known methods.

\subsection{Laser cooling}

  Recently~\cite{TSB1}, a new method of beam production with small
emittance was considered --- laser cooling of electron beams --- which
allows, in principle, to reach $\mbox{$L_{\gamma\gamma}$} \geq 10^{35}$ \mbox{cm$^{-2}$s$^{-1}$}.

The idea of laser cooling of electron beams is very simple.  During a
collision with optical laser photons (in the case of strong fields it
is more appropriate to consider the interaction of an electron with an
electromagnetic wave) the transverse distribution of electrons
($\sigma_i$) remains almost the same. Also, the angular spread
($\sigma_i^{\prime}$) is almost constant, because electrons loss  momenta
 almost along their
trajectory (photons follow the initial electron trajectory with a
small additional spread). So, the emittance $\mbox{$\epsilon_i$} = \mbox{$\sigma_i$} \mbox{$\sigma_i^{\prime}$}$ remains
almost unchanged. At the same time, the electron energy decreases from
$E_0$ down to $E$. This means that the transverse normalized
emittances have decreased: $ \mbox{$\epsilon_n$} = \gamma \mbox{$\epsilon$} = \mbox{$\epsilon_n$}_0(E/E_0)$.  One can
reaccelerate the electrons up to the initial energy and repeat the
procedure. Then after N stages of cooling $ \mbox{$\epsilon_n$} /\mbox{$\epsilon_n$} _0 = (E/E_0)^N$
(if \mbox{$\epsilon_n$}\ is far from its limit).

 Some possible set of parameters for  laser cooling is: $E_0 = 4.5$ GeV,
$l_e=0.2 $ mm, $\lambda = 0.5$ \mbox{$\mu$m}, flash energy $A \sim 10 $ J.  The
final electron bunch will have an energy of 0.45 \mbox{GeV}\ with an energy
spread $\sigma_E/E \sim 13 \%$, the normalized emittances \mbox{$\epsilon_{nx}$},\mbox{$\epsilon_{ny}$}\
are reduced by a factor of 10.  A two-stage system with the same parameters
reduces the emittances by a factor of 100. The limit on the final
emittance is $\mbox{$\epsilon_{nx}$}\ \sim\mbox{$\epsilon_{ny}$}\ \sim2\times 10^{-9}\;$ m~rad at
$\beta_i= 1\; \mbox{mm}$.  For comparison, in the TESLA (NLC)
project the damping rings have $\mbox{$\epsilon_{nx}$}\ =14(3)\times 10^{-6}\;$ m~rad,
$\mbox{$\epsilon_{ny}$}\ =25(3)\times 10^{-8}\;$ m~rad.

This method requires a laser system even more powerful than that for e
$\to \gamma$ conversion. However, all the requirements are reasonable,
taking into account the fast progress of laser technique and time
plans of linear colliders. Multiple use of the laser bunch (some
optical resonator) can  considerably reduce the required  average laser power.
This method can be tested already now at a low repeatition rate.

\section{Acknowledgements}

 I would like to thank Clem Heusch for a nice Workshop, which was
 one of the important steps towards \mbox{e$^+$e$^-$}, ee, \mbox{$\gamma$e}, \mbox{$\gamma\gamma$}\ colliders.

High Energy Physics - Experiment (hep-ex); High Energy Physics - Phenomenology (hep-ph); Accelerator Physics (physics.acc-ph); Optics (physics.optics)


The HST Survey of BL Lacertae Objects. II. Host Galaxies

Abstract: We have used the HST WFPC2 camera to survey 132 BL Lac objects comprising seven complete radio-, X-ray-, and optically-selected samples. We obtained useful images for 110 targets spanning the redshift range 0 < z < 1.3. In two thirds of the BL Lac images, host galaxies are detected, including nearly all for z < 0.5 (58 of 63). The highest redshift host galaxy detected is in a BL Lac object at z=0.664. In 58 of the 72 resolved host galaxies, a de Vaucouleurs profile is significantly preferred, at >99% confidence, over a pure exponential disk; the two fits are comparable in the remaining 14 cases. These results limit the number of disk systems to at most 8% of BL Lacs (at 99% confidence), and are consistent with all BL~Lac host galaxies being ellipticals. The detected host galaxies are luminous ellipticals with a median absolute K-corrected magnitude of M_R= -23.7 +- 0.6 mag, at least one magnitude brighter than M* and comparable to brightest cluster galaxies. The galaxy morphologies are generally smooth and undisturbed, with small or negligible ellipticities (<0.2). There is no correlation between host galaxy and observed nuclear magnitude or estimated jet power corrected for beaming. If black hole mass is correlated linearly with bulge mass in general, this implies a large range in Eddington ratio. Present data strongly support the unification picture with FR I galaxies constituting the bulk of the parent population of BL Lac objects.

\section{Introduction}
\label{sec:intro}

The {\it Hubble Space Telescope (HST)} has been used extensively to
study the host galaxies of AGN, primarily quasars and radio galaxies
with relatively high nuclear luminosities
(Disney et al. 1995; McLeod \& Rieke 1995; Bahcall et al. 1997;
Best, Longair \& R\"ottgering 1997; Hooper, Impey \& Foltz 1997;
McCarthy et al. 1997; Ridgway \& Stockton 1997;
Serjeant, Rawlings \& Lacy 1997; Boyce, Disney \& Bleaken 1999;
McLeod, Rieke \& Storrie-Lombardi 1999; McLure et al. 1999).
Its order-of-magnitude better spatial resolution (over
a large field compared to adaptive optics) provides unique and
critical information at sub-arcsecond scales.

Collectively, {\it HST}
and ground-based observations of host galaxies have already led to
interesting results. The idea that radio-loudness is uniquely related
to host galaxy type has been shown to be incorrect:
while radio-loud AGN are almost always found in elliptical
galaxies --- often luminous ellipticals comparable to brightest
cluster galaxies (Smith \& Heckman 1986; Hutchings, Janson \& Neff 1989;
Veron-Cetty \& Woltjer 1990; Taylor et al. 1996; Zirm et al. 1998) ---
radio-quiet AGN are found in both elliptical or disk galaxies (Taylor
et al. 1996, Bahcall et al. 1997, McLure et at. 1999).
Several studies have reported that
the host galaxies of radio-quiet AGN are systematically less luminous
than those of radio-loud AGN (Smith \& Heckman 1986;
Hutchings et al. 1989; Veron-Cetty \& Woltjer 1990;
Lowenthal et al. 1995)
although this effect was not found in other samples matched for
redshift and luminosity (Taylor et al. 1996; Hooper et al. 1997).
It has also been suggested that radio-quiet AGN
have less disturbed morphologies (Hutchings et al. 1989),
and certainly dust, tidal tails, and/or close companion galaxies are
prevalent in radio-loud AGN (Smith \& Heckman 1986; Yee \& Green 1987;
Bahcall et al. 1997; Canalizo \& Stockton 1997; Martel et al. 1997;
Pentericci et al. 1999).

An interesting but controversial issue is the possible relation
between host galaxy magnitude and nuclear brightness.
A trend for the brightest nuclei to lie in the most luminous galaxies
has been found in some
(McLeod \& Rieke 1994, 1995; Hooper et al. 1997; McLeod et al. 1999)
though not all
(Taylor et al. 1996; Wurtz, Stocke \& Yee 1996, hereafter WSY)
host galaxy studies.
Such a correlation would indicate a close connection between
small-scale, black-hole-related phenomena and large-scale galactic
phenomena, possibly related to galaxy formation scenarios (e.g., Small
\& Blandford 1992; Haehnelt \& Rees 1993). Further, where black hole
masses have been reliably estimated in local galaxies they appear to
be proportional to the bulge mass (Kormendy \& Richstone 1995,
Magorrian et al. 1998, van der Marel 1999);
if the efficiency of converting accreting mass into AGN
luminosity does not vary widely, then for AGN this would translate to
AGN luminosity being proportional to host galaxy magnitude (e.g.,
McLure et al. 1999).

A related issue is whether the cosmic evolution of galaxies and AGN
is inextricably linked or essentially independent.
That is, does the central black
hole grow more or less independently of the stellar mass, or is there
significant feedback between small-scale and large-scale systems?
Many have noted the similarity of the cosmic evolution of star-forming
galaxies and of AGN (e.g., Silk \& Rees 1998) --- both peak somewhere
in the redshift range $z=1-3$ --- possibly indicating a close
connection between galaxy and black hole evolution. So far, there is
evidence, at least in radio-loud AGN, that AGN host galaxies contain
old stars, as if the galaxy formed at high redshift (Dunlop et al. 1996;
Ridgway \& Stockton 1997; Best et al. 1998; De Vries et al. 1998),
close to the epoch of quasar dominance (Foltz et al. 1992;
Maloney \& Petrosian 1999). It is not yet clear whether this is the case for all
radio-loud AGN or only for the most luminous (those with the most massive
black holes).

Despite extensive observations, the results to date on host galaxies
are surprisingly mixed, possibly because the samples studied contain
mostly higher luminosity AGN and are often far from complete, in part
because of selection effects. We therefore
undertook an {\it HST} snapshot survey\footnote{Based on observations
made with the NASA/ESA Hubble Space Telescope, obtained at the Space
Telescope Science Institute, which is operated by the Association of
Universities for Research in Astronomy, Inc., under NASA contract
NAS~5-26555.} to investigate the morphology, color, and luminosity of
AGN host galaxies, and the evolution of these properties with cosmic
epoch. Our strategy was to complement existing investigations by
investigating lower luminosity AGN out to moderate redshifts, $z\sim 1$.
We also concentrated on radio-loud AGN, which seem less affected by
dense gaseous environments than radio-quiet AGN and which therefore
may be simpler systems on galactic scales. According to the current
paradigm, radio-loud AGN all have relativistically outflowing jets,
and depending on the orientation of the jet, they present markedly
different appearances to the observer (Urry \& Padovani 1995). The
jet luminosity is a critical parameter, influencing the spectral
energy distributions (Sambruna et al. 1996; Fossati et al. 1997, 1998)
and radio morphological properties (Baum, Zirbel \& O'Dea 1995).
BL~Lac objects have intrinsically lower luminosities than radio-loud quasars
but can be seen to moderately high redshifts thanks to fortuitous
relativistic beaming, hence they form the ideal sample for our study.
BL~Lac nuclei are also relatively weaker than in beamed quasars,
making them even easier targets for host galaxy studies.

Our well-defined survey of more than a hundred BL~Lac objects included the
full range of BL~Lac types (Padovani \& Giommi 1995), at redshifts
out to $z\sim1.3$.
The observations and data analysis are described very briefly in
\S~\ref{sec:obsda}, with details given by Scarpa et al. (2000).
Results are given \S~\ref{ssec:results} and discussed further in
\S~\ref{sec:disc}. Conclusions are given in \S~\ref{sec:conc}. For
ease of comparison to the published literature, we used \cosmo
throughout the paper.

\section{Observations and Data Analysis}
\label{sec:obsda}

\subsection{The BL~Lac Sample}
\label{ssec:sample}

The BL~Lac snapshot survey is based on seven complete samples
selected at radio, optical, and X-ray wavelengths (see Table~1).
We deliberately targeted both radio-selected and X-ray-selected
BL~Lac objects because the sample content varies strongly with
selection wavelength (Ledden \& O'Dell 1985; Stocke et al. 1985).
Specifically, there is a selection effect for BL~Lac ``type''
because the two broad components
in BL~Lac spectral energy distributions (SEDs)
have peak power outputs ($\nu L_\nu$) at wavelengths that
increase systematically with luminosity
(Sambruna, Maraschi \& Urry 1996,;
Fossati et al. 1997,;Ulrich, Maraschi \& Urry 1997).
``Red'' BL~Lacs, also known as low-frequency-peaked BL~Lacs (LBL),
have SEDs peaking at infrared-optical wavelengths and in the
MeV-GeV gamma-ray band, and have luminosities approaching those of quasars.
``Blue'' or high-frequency-peaked BL~Lacs (HBL)
have SEDs peaking at UV-X-ray wavelengths and again at TeV energies,
and are generally less luminous.
Because of these diverse spectral shapes,
``red'' BL~Lacs dominate existing radio-selected samples
and ``blue'' BL~Lacs dominate most X-ray-selected samples
(Padovani \& Giommi 1995),
although there is clearly a continuous distribution of SED shapes
between these extrema (e.g., Laurent-Muehleisen et al. 1998;
Fossati et al. 1997; Perlman et al. 1998).
Since BL~Lac SEDs are dominated by beamed emission from
aligned relativistic jets, from radio through gamma-ray wavelengths
(Urry \& Padovani 1995; Ulrich et al. 1997),
the range of SED shapes sampled in our {\it HST} snapshot survey
corresponds to the full range of jet physics in BL~Lac objects.

The final list of 132 BL~Lac objects (some are in more than one
sample) was approved for snapshot observations in Cycle~6,
and in the end 110 were observed.
These spanned the redshift range $0.027 \leq z \leq 1.34$,
with a median redshift of $\langle z \rangle =0.29$
and 22 having $z>0.5$; the distribution of redshifts is shown
in Figure~1.
Ten of the 132 were also observed with WFPC2 in Cycle~5,
for longer exposures in a different filter
(Falomo et al. 1997; Jannuzi et al. 1997; Yanny et al. 1997;
Urry et al. 1999).

\subsection{HST Observations}
\label{ssec:obs}

The list of observed objects is given in Table~2, along with the
redshift and SED type (HBL or LBL). Scarpa et al. (2000) give a more
detailed journal of the observations, including the BL~Lac position,
date of observation, and exposure information.
Observations were done with the {\it HST} WFPC2 and the F702W filter,
a sensitive, red, broad-band filter that minimizes contamination from
extended emission line gas (which is in any case much less important
in BL~Lacs than in quasars and radio galaxies), dust, and recent star
formation. In the few cases for which WFPC2 F814W images already existed,
we used the F606W or F555W filters instead,
to get a broader baseline for estimating colors.
The scheduling was done in snapshot mode, meaning the observations fit
random holes in the {\it HST} schedule for regular GO observations.
Thus, the final list of 110 observed BL~Lac objects constitutes an
unbiased subset of the original target list.

To obtain for each target a final image well exposed
both in the inner, bright nucleus and in the faintest outer regions
of the host galaxy, we used a series of exposures ranging from a few tens
of seconds to as long as $\sim 1000$ seconds. From 3 to 5
images were obtained for each target, and were
later combined to remove cosmic-ray events and to increase the signal-to-noise
ratio of the final image. The median exposure times are
480~seconds for $z<0.5$ and 840~seconds for $z>0.5$.

\subsection{Data Reduction and Galaxy Surface Brightness Profile Fitting}
\label{ssec:datared}

Data reduction was carried out as summarized by Urry et al. (1999) and
described in more detail by Scarpa et al. (2000), who also show
the final summed images. Magnitudes reported here are in
the Cousins system, transformed from {\it HST} magnitudes assuming colors
appropriate for a redshifted elliptical galaxy (for details, see
Scarpa et al. 2000). For the host galaxies detected here ($z\lesssim 0.7$),
the color corrections are $\lesssim 0.6$~mag because the
Cousin R and WFPC2 F702W filters are similar,
as are the Johnson V and WFPC2 F606W.

We estimated the expected amount of reddening
due to interstellar matter in our Galaxy
from HI column densities, using the
conversion $\log N_H /E(B-V) = 21.83$~cm$^{-2}$ mag$^{-1}$ appropriate
for high latitudes
(Shull \& Van Steenberg 1995), assuming a total-to-selective extinction
$A_R = 2.3 E(B-V)$ (Cardelli, Clayton \& Mathis 1989).
In general, the reddening is
quite small, with median value 0.2~mag at R;
values for each object are in Table~2.
These corrections were applied to the reported absolute magnitudes,
to give our best estimate of the intrinsic physical quantity,
but not to the apparent magnitudes, which reflect directly measured
quantities.
In any case, the reddening corrections are usually comparable to or
smaller than the estimated systematic uncertainties (\S~\ref{ssec:errors}).
In three cases $A_R$ is as high as $\sim 2.5$, but for these BL~Lacs
the redshift is not known, so they do not affect our conclusions about
absolute quantities.
Note that we had no constraints on, and therefore did not correct for,
reddening in the host galaxy or BL Lac nucleus itself.

To evaluate the morphologies and apparent magnitudes of the BL~Lac
host galaxies, most of which are quite smooth and round, we fitted
one-dimensional surface brightness profiles. This is computationally
much simpler than two-dimensional analysis and even for our
well-exposed (2-orbit) Cycle~5 images gave equivalent results (Falomo
et al. 1997; Urry et al. 1999); extensive two-dimensional analysis
for the low-redshift ($z<0.3$) BL~Lacs is described by Falomo et al. (2000).
Azimuthal averaging also improves the
signal-to-noise ratio in the outermost parts of the host galaxy,
allowing us to go to fainter surface brightnesses.

Details of the fitting procedure can be found in Scarpa et
al. (2000), and Urry et al. (1999). Briefly, we fitted the profile
with a galaxy plus point source, convolved with the point spread function (PSF),
adjusting the parameters
simultaneously to determine a best-fit and statistical errors using
the chi-squared statistic. The PSF consisted
of a Tiny Tim model (Krist 1995) in the inner 2~arcsec joined smoothly
to a composite stellar profile for the wings (the pure Tiny Tim PSF
model does not include large angle scattered light, leading
to overestimates of the host galaxy brightness; see Fig.~2 of
Scarpa et al. 2000).

We tested both exponential disk and de~Vaucouleurs $r^{1/4}$ models
for the galaxy. We used an F-test to evaluate which if either of the two galaxy
models
was preferred, at 99\% confidence or better. Our
threshold for formal detection of a host galaxy was that
PSF-plus-galaxy fit be better at 99\% confidence than the PSF-only
fit. For unresolved objects we determined 99\% confidence upper
limits (statistical errors) to the host galaxy magnitudes ($\Delta
\chi^2 = 6.6$ for one parameter of interest, $M_{gal}$), fixing the
half-light radius at $r_e =10$~kpc, slightly larger (to be
conservative) than the median ($<r_e>=8.5$~kpc) for the
72 resolved objects.

We did use two-dimensional analysis to test for decentering of the nucleus
with respect to the host galaxy, for the $\sim30$ well-resolved cases
(Falomo et al. 2000).
With the exception of peculiar cases like double nuclei (Scarpa et
al. 1999), the point sources are well-centered in the host galaxies,
with a tolerance generally better than 0.03~arcsec. Extensive
two-dimensional analysis carried out on the $\sim30$ nearest host
galaxies ($z<0.3$), where spatial resolution and signal-to-noise ratio
are highest, also shows very small ellipticities (generally $\epsilon
< 0.2$) and few cases of isophotal twists or distortions (analysis and
results described fully by Falomo et al. 2000).

\subsection{Systematic Errors}
\label{ssec:errors}

Comparison of our fitted host galaxy magnitudes with those obtained
by other authors, even on the same data (e.g., Jannuzi et al. 1997)
reveals systematic discrepancies of up to several tenths, even for
bright, easily detected elliptical host galaxies. Sources of
systematic error include uncertainty in the PSF shape, variations in
how the PSF is normalized, and uncertainty in the sky background. We
have done extensive simulations to estimate the size of the systematic
uncertainties in our derived magnitudes, as reported by Scarpa et
al. (2000); here we mention the results.

The PSF shape we have adopted and the fitting procedure recover very
accurately the input parameters in simulated data, with uncertainties in
the total galaxy magnitude of less than $0.15$~mag (Gaussian
half-width of the distribution of measurement minus true value; Scarpa
et al. 2000).
The uncertainty in the sky background is more significant, particularly
when the galaxy is only marginally resolved. From
the simulations we estimate that the typical systematic
uncertainty in measured host galaxy magnitude is $\pm0.2$~mag when the
point source brightness is within 2~mag of the host galaxy. This
uncertainty of a few tenths is over and above the statistical
uncertainties quoted in Table~2.

We also note that derived galaxy properties reported in the literature
can differ widely because of different calibrations, different
aperture sizes for photometry, and different fitting assumptions, as
well as the usual differences in cosmology. Conversion to absolute
magnitude can introduce further discrepancies because K~corrections in
the literature vary widely (King \& Ellis 1985; Frei \& Gunn 1994;
Fukugita, Shimasaku \& Ichikawa 1995; Kinney et al. 1996). For elliptical
galaxies, published values differ by 0.1~mag at $z\lesssim0.2$ and as
much as 0.5~mag at $z\sim1$. The range of values for spiral galaxies is
similar or perhaps even larger. The difference between the
K~corrections for E-type and Sb-type spectra is of course much larger
(as much as 2~mag at $z\sim1$). The K~correction values we used are
given in Table~2.

\subsection{Results of Host Galaxy Fits}
\label{ssec:results}

Results of the one-dimensional fitting for de~Vaucouleurs models are
summarized in Table~2, along with the 68\% confidence statistical
uncertainties (in most cases, systematic errors dominate; see
\S~\ref{ssec:errors}). Plots of the radial surface brightness
profiles with best-fit de~Vaucouleurs model and residuals, and of the
chi-squared confidence contours for the two parameters of interest
(host galaxy magnitude and effective radius), are shown by Scarpa et
al. (2000), along with images and discussion of individual sources.

In 72 of the 110 BL~Lac objects, host galaxies are detected.
This is strongly dependent on redshift, since our relatively short
exposures become insensitive to $L^*$ galaxies for $z\gtrsim 0.5$.
Figure \ref{fig:histo} shows detection fraction as a function of
redshift for the observed BL~Lac sample.

Like the host galaxy detection rate, other trends with redshift result
from declining signal-to-noise ratios and/or absolute spatial resolution
with increasing distance of the BL~Lac.
For example, the absolute magnitude of the nucleus increases with
redshift, as expected in a flux-limited sample (Fig.~\ref{fig:Mnuc_z}).
Note that each BL~Lac type spans nearly the whole redshift range,
although there is redshift segregation because the ``red'' BL~Lacs
are systematically more luminous than the ``blue'' BL~Lacs.

The half-light radius of detected host galaxies increases slightly
with redshift (Fig.~\ref{fig:re_z}),
corresponding to the larger sizes of more luminous host
galaxies (less luminous host galaxies being harder to detect at high
redshift). The uncertainties in the fitted values of $r_e$ also increase
with redshift. The median value over all redshifts is
$\langle r_e \rangle = 8.5$~kpc.

Because the flux limits of the input sample were almost entirely unaffected
by optical flux, the host galaxy magnitude is in principle not correlated
with redshift. However, two selection effects affect the magnitude range of
detected host galaxies. First, faint host galaxies will not be detected
around bright nuclei, and second, very bright host galaxies with weak nuclei
will have been classified as galaxies rather than AGN.
Figure~\ref{fig:nuc/host} shows how the ratio of nuclear to host galaxy
luminosity is confined to a relatively narrow range by these two effects.

\section{Discussion}
\label{sec:disc}

\subsection{Luminosities of the BL~Lac Host Galaxies }
\label{ssec:lum}

Table~2 lists the absolute magnitudes of the host galaxy and nucleus
(i.e., point source) for each observed BL~Lac object, calculated from
the best-fit de~Vaucouleurs model parameters using the K~corrections
listed. The 72 detected hosts are very luminous, round galaxies.
Their median absolute magnitude is $\langle M_R \rangle = -23.7$~mag,
with a relatively small dispersion of 0.6~mag. These results are
largely in agreement with previous, mostly smaller, surveys of BL~Lac
objects (Abraham, Crawford \& McHardy 1991;
Stickel, Fried \& K\"{u}hr 1993; Falomo 1996; WSY),
although for individual objects the differences average
$\sim 1$~mag (see Scarpa et al. 2000).
Because we probe
higher redshifts on average than these ground-based surveys, it is not
surprising that our detected host galaxies are also somewhat more
luminous on average.

Figure~\ref{fig:Mr_z} shows the distribution of our absolute host
galaxy magnitudes for the 85 BL~Lacs with known redshifts. The
median value (dashed line) is roughly one magnitude brighter than
$L^*$; $L^*_R =-22.4$~mag at low redshift (Efstathiou, Ellis
\& Peterson 1988, converted from $L^*_B$ assuming $B-R=1.56$). The
BL~Lac host galaxies are similar in luminosity to brightest cluster
galaxies (Taylor et al. 1996; WSY), $M_R = -23.9$~mag (Thuan \&
Puschell 1989, converted from the H band assuming $R-H=2.5$) or to
Fanaroff-Riley type~I radio galaxies (Ledlow \& Owen 1996; cf. WSY),
which are often found in moderate to rich cluster environments.

At high redshifts, we have many upper limits to the host galaxy
magnitudes. Most are uninteresting because the nuclei are quite bright
but a few are faint, $\gtrsim -23$~mag, indicating at least a few high
redshift BL~Lacs ($z\gtrsim 0.3$) have $L^*$-like host galaxies, like
the lower luminosity hosts at $z\lesssim0.1$.

\subsection{Morphologies of the BL~Lac Host Galaxies }
\label{ssec:morph}

In the vast majority of cases, a de~Vaucouleurs $r^{1/4}$ model fit
the data significantly better than an exponential disk model. In only
14 cases, all with relatively low signal-to-noise ratios, were the
two fits even comparable. (In one case, 0446+449, the disk fit was
unequivocally preferred but there was no point source present,
indicating the identification as a BL~Lac is in error;
see Scarpa et al. 2000 for a full
discussion of all dubious classifications.)
We are not biased against finding disks, since their surface brightness
would fall off slower than the $r^{1/4}$ profile.
Given our large sample of resolved host galaxies,
we can say at the 99\% confidence level
that at most $\sim8$\% can be in disk systems, and
our results are consistent with all BL~Lac objects being found
exclusively in elliptical galaxies.

Table~2 lists the half-light radii of the host galaxies, in
kiloparsecs, from the best-fit de~Vaucouleurs model. As well as being
luminous, the host galaxies of BL~Lac objects are large, and the
larger galaxies tend to be more luminous.
Figure~\ref{fig:mue_re} shows the relation between half-light radius, $r_e$,
and surface brightness at that radius, $\mu_e$, for the detected BL~Lac
host galaxies with known redshifts. The data describe a linear
trend such that larger, more luminous galaxies have lower central
surface brightnesses. This has been seen in many samples of
ellipticals, in clusters or out, whether radio-loud or not, and is
basically a projection of the fundamental plane for elliptical
galaxies (Djorgovski \& Davis 1987, Hamabe \& Kormendy 1987). The
best-fit correlation for the BL~Lac host galaxies,
$\mu_e = (3.9\pm 0.9) \log (r_e/{\rm kpc})  + (17.2\pm0.7)$ mag arcsec$^{-2}$,
is consistent with those reported for FR~I
radio galaxies (Govoni et al. 2000), bright cluster
ellipticals (BCE; Ledlow \& Owen 1995), and non-cluster ellipticals
(Kormendy 1977). It also agrees well with the slope determined
for other, more powerful AGN (McLure et al. 1999).

The implication is that BL~Lac host galaxies have absolutely normal
elliptical morphologies and are somewhat more luminous than average.
Because they are very round, there is no obvious
alignment with the more linear radio structures.
The BL~Lac host galaxies also show normal color profiles
(Kotilainen, Falomo \& Scarpa 1998; Urry et al. 1999), and
follow quite well a $r^{1/4}$ law, so the nuclear activity
appears to have markedly little effect on the galaxy properties.
The integrated colors, where available, are
consistent with redshifted emission from a passively
evolving elliptical galaxy with an old stellar population
(rest-frame colors $R-I = 0.70$, $V-I= 1.31$) and
imply an initial star formation epoch for the host galaxies
of at least $\sim6$~billion years ago (Bruzual \& Charlot 1993).
Bluer data are required to constrain new star formation,
to which our WFPC2 F814W images are generally not sensitive
at these low redshifts.

Finally, we note that Figure~\ref{fig:mue_re} includes the 14
morphologically unclassified host galaxies (for which disk and de Vaucouleurs
fits gave similar chi-squared values). That they fit nicely
into the $\mu_e - r_e$ relation for elliptical hosts supports the
idea that these galaxies are indeed ellipticals.

\subsection{Host Galaxies of ``Red'' and ``Blue'' BL~Lacs}
\label{ssec:redblue}

There are systematic differences between ``red'' and ``blue'' BL~Lac
objects --- in luminosity, redshift, and spectral energy distribution.
The ``blue'' objects have less
luminous nuclei and jets with lower kinetic powers (Celotti, Padovani
\& Ghisellini 1997), and dissipate most of their energy in
synchrotron radiation from highly relativistic electrons.
The ``red'' objects, which have systematically higher bolometric
luminosities, are probably redder because the highest energy
electrons cool quickly by Compton scattering ambient UV and X-ray
photons to gamma-ray energies, which can dominate the bolometric
output (Ghisellini et al. 1998).
These two classes of BL~Lac object therefore
reflect two different kinds of jets (probably extrema of a continuous
distribution), which result from different jet formation and/or
evolution.

Despite these strong nuclear trends, we find no differences between
the host galaxies of ``red'' and ``blue'' BL~Lacs,
either in luminosity or size, confirming the earlier result by WSY
for somewhat fewer objects.
This strongly suggests that nuclear
properties, which can strongly influence jet formation and
propagation, do not have a dramatic effect on large-scale host galaxy
properties (or vice-versa).

\subsection{Comparison to Radio Galaxies}
\label{ssec:rg}

According to unified schemes (Barthel 1989), BL~Lac objects are
FR~I radio galaxies whose jets are aligned along the line of sight
(Urry \& Padovani 1995). This implies BL~Lac host galaxies should be
statistically indistinguishable from FR~I host galaxies. It has been
suggested that the parent population of BL~Lacs might instead be
FR~IIs, or a subset thereof (Kollgaard 1987; Urry \& Padovani 1995; WSY;
Laing 1994). Our host galaxy study directly tests this alternative
unification hypothesis.

The original division between FR~I and FR~II galaxies was
morphological --- whether hot spots occurred at the inner or outer
edges of the radio source, respectively --- and the excellent
correlation of morphology with radio luminosity was noted at the same
time (Fanaroff \& Riley 1974). For low-frequency radio luminosities
below (above) $P_{178} = 2 \times 10^{25}$~W~Hz~$^{-1}$~sr$^{-1}$,
almost all radio sources were type I (II). This clean separation in
luminosity disappears at higher radio frequencies, where the overlap
can be several decades in radio power.

Owen and Ledlow (1994) showed that FR~I/II division depends on both
radio power and optical luminosity, with a diagonal line dividing FR~Is
from IIs.
If the observed radio power is a measure of the kinetic power of the
jet, and if optical luminosity correlates with the mass of the host
galaxy, the Fanaroff-Riley division can be explained in (at least) two
ways. In the ``nurture'' scenario, more massive galaxies have a
denser interstellar medium better able to decelerate an outflowing
relativistic jet (Bicknell 1995). For a given jet power, the FR~Is
would be in more luminous galaxies than FR~IIs (i.e., to the right of
the diagonal dividing line); or for a given galaxy mass,
the FR~Is would have less powerful jets than FR~IIs (i.e., below the diagonal
line).
In contrast, in one ``nature'' scenario,
FR~Is and FR~IIs are distinguished at birth because the power
delivered to the jet depends on a magnetic switch that essentially
links higher power jets with more massive, spinning black holes (Meier
1999). A correlation between black hole mass and galaxy mass then
leads to the diagonal FRI/II dividing line.

Figure~\ref{fig:pext_opt} shows a new version of the Owen \& Ledlow
(1994) diagram of radio power versus optical magnitude. Because we
plot extended radio power instead of total radio power, and host
galaxy magnitude rather than total magnitude, relativistic beaming of
BL~Lac nuclei has no effect, and thus a direct comparison between the
host galaxies of BL~Lac objects and radio galaxies is possible.
We took FR~I and II galaxies from the 2~Jy sample (Wall \& Peacock 1985)
because it has similar depth and selection criteria as the 1~Jy BL~Lac
sample (Stickel et al. 1991); morphological classifications are from
Morganti, Killeen \& Tadhunter (1993); values for the extended
radio power are from references listed in Table~2; and we
restricted all samples to $z<0.5$ to avoid luminosity-redshift biases.
Because the sample selection biases still differ, one cannot compare
the distributions in extended radio power and host galaxy magnitude;
rather, unified populations should simply occupy similar regions in
Figure~\ref{fig:pext_opt}.

The BL~Lacs overlap extremely well with the FR~I galaxies, with only a
few in the FR~II region.
Similarly, radio-loud quasars (also restricted to $z<0.5$)
lie in the FR~II region of the diagram.
Thus, the present data strongly support the
unification picture with FR~I galaxies constituting the bulk of the parent
population.

Note that the projection of this plot onto the host galaxy magnitude
axis will give statistically distinguishable distributions for BL~Lacs
and FR~Is; formally, in this one dimension alone, FR~IIs might appear
to be a better match (WSY). This is a misleading approach, however, since it
ignores important information about radio power. As Figure~\ref{fig:pext_opt}
clearly shows, BL~Lac objects are not well matched to FR~II radio galaxies.
Instead, what confuses the one-dimensional approach is that BL~Lacs (so far)
have not been found in host galaxies as luminous as the most luminous
FR~Is, nor have they been found in clusters as rich as those FR~Is.
If the absence of very luminous host galaxies and/or rich cluster
environments is significant (Owen, Ledlow \& Keel 1996), it is possible
that dense intracluster environments or extremely massive host galaxies
completely quench any would-be relativistic jet.

\subsection{Near Environments and Close Companions}
\label{ssec:near}

Further refinement of the unification picture is possible from
consideration of the larger environments of the BL~Lacs (Fried et
al. 1993; Falomo, Pesce, \& Treves 1993, 1995; Pesce, Falomo, \&
Treves 1994, 1995; Smith, O'Dea \& Baum 1995; Wurtz et al. 1997).
Since FR~I radio galaxies commonly occur in clusters, so should BL~Lac
objects. The present data add to previous work by allowing detection
of fainter companions closer to the BL~Lac nucleus. The small field of
view, however, limits the statistics with which the possible excess of
companion galaxies can be assessed.

Preliminary results for the environments of BL~Lacs indicate a large
number of objects with companions, some as close as 5~kpc (projected).
Defining ``companion'' galaxies as
those within 70~kpc of the BL~Lac nucleus and brighter than
1~magnitude below $m^*$ (for those objects without measured redshift,
$z=0.2$ was used), we find companion galaxies in 42\% of the
BL~Lac sample. Without the magnitude limit, so that fainter galaxies
are included, 47\% of the sample has companions.
For comparison, 42\% of a sample of low redshift FR~I
galaxies have companion galaxies within the same radius (Pesce et al.,
in preparation.).

On the larger scale environment, BL~Lacs have been seen to lie in
regions of enhanced galaxy density, on average. Typically,
the clusters around BL~Lacs are poor, of Abell richness class 0-1,
although a few richer clusters are detected. Our preliminary results
for the snapshot survey are similar, with 40\% of the sample showing
regions of enhanced galaxy density. Nonetheless, a
significant number of objects appear to be completely
isolated (i.e., no close companions and no surrounding galaxies
above the average background).
A comprehensive analysis of the environments of BL~Lac objects,
determined from the {\it HST} images, will be given in separate papers
(Pesce et al. and Falomo et al., in preparation).

If mergers are a significant part of the galaxy formation process
(especially ellipticals) then hosts should be more disturbed at high
redshift. The BL~Lac objects in our sample generally appear
undisturbed, with a few exceptions (details of individual sources are
given by Scarpa et al. 1999, 2000).
This is in marked contrast to the case for more powerful radio sources.
Since among BL~Lacs, jet power has no discernible effect on galaxy morphology,
this suggests that instead age may be important.
Radio sources having undergone recent mergers would be more likely to show
dust lanes, tidal tails, and the like, even if most of the stars were formed
at high redshift. In contrast, the relaxed morphologies
of BL~Lac host galaxies suggest they are old, more evolved sources
that have not recently merged.
The observation that radio-quiet AGN have less disturbed morphologies
than radio-loud AGN (Hutchings et al. 1989) could then
be a function more of intrinsic AGN luminosity and/or evolutionary
state than radio-loudness, since our radio-loud sample is markedly undisturbed.
Comparison to observations of BL~Lac objects at higher redshifts ($z>0.5$)
will be very illuminating on this point.

\subsection{Comparison of Nucleus and Host Galaxy }
\label{ssec:nuc_hg}

A correlation or trend between galaxy magnitude and nuclear brightness has been
reported in several host galaxy studies (McLeod \& Rieke 1994, 1995;
Hooper et al. 1997; McLeod et al 1999). This can be interpreted as an
extension of the
correlation between black hole mass and bulge mass in nearby
ellipticals (Kormendy \& Richstone 1995; Magorrian et al. 1998;
van der Marel 1999), providing the
Eddington ratio does not vary widely among the AGN considered
(cf. McLure et al. 1999).

In our sample of BL~Lac objects, there is a slight correlation
between measured nuclear and host galaxy luminosities,
but it becomes insignificant when upper limits are included.
This can be seen in Figure~\ref{fig:hg_nuc},
which shows that the host galaxy magnitudes
cluster near the median value, $\langle M_R \rangle = -23.7$~mag,
independent of the luminosity of the nucleus.
Furthermore, the best-fit slope of a linear relation
is much shallower than implied by the correlation between black hole mass and
bulge
mass for fixed Eddington ratio. If the bulge-black hole correlation
translates to $L_{gal}$-$L_{nuc}$ in our data, then the Eddington ratio
must range over at least two orders of magnitude among otherwise
similar jet sources.

The observed point source magnitude, which is dominated by synchrotron
emission from an unresolved jet, is likely affected by relativistic
beaming. This could cause the points plotted in Figure~\ref{fig:hg_nuc}
to extend across an artificially large range in point source magnitude,
possibly washing out an underlying correlation.
Estimates for the Doppler factor (actually lower limits)
are available for only a fraction of our target sources
(Burbidge \& Hewitt 1987; Xie et al. 1991; Dondi \& Ghisellini 1995),
so wholesale correction of
the observed nuclear magnitudes is not presently feasible.
Instead we considered whether extended radio power --- which
correlates with jet power and is unaffected by beaming ---
was correlated with host galaxy magnitude. Including upper limits,
there is no significant correlation; ignoring upper limits, there is
a marginal correlation, with a slope much shallower than that implied
by the bulge-black hole relation.

It remains to be explained why a nucleus-galaxy correlation is seen in some
other
samples and not in the present sample of BL~Lac objects. One
possibility is that there is a luminosity threshold for the effect and
that it does not appear in low-luminosity AGN, as suggested by
McLeod \& Rieke (1995). Our sample includes
some luminous AGN at $M_R < -25$~mag, the region where
McLeod \& Rieke (1995) found a correlation in their quasar sample,
but few have detected host galaxies and the nuclear magnitudes are
affected by beaming.
Combining our BL~Lac sample with quasar samples matched in redshift,
it should be possible to assess this issue directly.

Two selection effects could in
principle induce a spurious correlation, especially in data with the
low spatial resolution typical of ground-based observations: (1) the
difficulty of finding faint host galaxies around bright nuclei and (2)
the absence of AGN with bright host galaxies and weak nuclei (these
are identified as galaxies rather than AGN). For any given
investigation, simulations can indicate whether these effects are
significant. (We note that because McLeod \& Rieke 1995 detected host
galaxies for 100\% of the AGN in their sample, the correlation they
report should not be influenced by the first effect.)

The correlation of luminosity with redshift in flux-limited samples
could confuse the effects of evolution or steep luminosity functions
with physical effects like true nuclear/host galaxy relations.
To measure the latter effect definitively therefore require spanning
a large range of luminosity at a fixed redshift.
At present, conclusions drawn from flux-limited samples of limited
luminosity range at any one redshift, whether high luminosity (quasars)
or low-luminosity (BL~Lacs), must be considered tentative.

Comparing the positions of nucleus and host galaxy, we are able,
with the high spatial resolution of {\it HST}, to place
tight limits on any de-centering of the BL~Lac nucleus. If any of our
BL~Lac objects were actually background quasars microlensed by stars
in a foreground galaxy (which we are calling the host galaxy), there
could well be an offset between the position of the nucleus (the
amplified background quasar) and the lensing galaxy (Ostriker \&
Vietri 1985). Instead, we find that the nuclei are generally
well-centered in the host galaxy, with deviations typically less than
0.03~arcsec (Falomo et al. 1998, 2000). Thus there is no evidence
for microlensing occurring in a large fraction of our sample.

\section{Conclusions}
\label{sec:conc}

We have shown that with {\it HST} it is easy to detect and characterize
the host galaxies of low-luminosity AGN like BL~Lac objects,
up to moderately high redshifts, $z \sim 0.6$.
We detected host galaxies in almost all cases with $z<0.5$
(58 of 63), and in 6 of 22 with $z>0.5$.
The highest redshift BL~Lac object with a detected host galaxy
is 1823+568 at $z=0.664$ (Falomo et al. 1997).

The detected host galaxies are smooth, round, very luminous ellipticals,
well fitted with de~Vaucouleurs surface brightness profiles.
In most cases --- generally, where the signal-to-noise ratio
is high --- the $r^{1/4}$ law fits significantly better than an
exponential disk; in the remaining cases, neither fit is preferred.
Thus our data are consistent with all BL~Lac host galaxies being
ellipticals.

The median K-corrected absolute magnitude of the detected host
galaxies is $\langle M_R \rangle = -23.7$~mag, with a dispersion of
0.6~mag. This is more than 1~mag brighter than $L^*_R$ galaxy, and
comparable to brightest cluster galaxies, or to
Fanaroff-Riley type~I radio galaxies (Ledlow \& Owen 1996) which are
often found in moderate to rich cluster environments. This strongly
supports the unification of BL~Lac objects with low-luminosity radio
galaxies, and rules out the possibility, at least at these low
redshifts, that the parent population of a substantial fraction of
BL~Lacs is FR~II radio galaxies. Note that there is a decade or so of overlap
between the FR~I and II populations, and between the quasar and BL~Lac
populations, not quite the clean division that was the original paradigm.

The BL~Lac host galaxies follow the same trend in
the $\mu_e$-$r_e$ projection of the fundamental plane as
other luminous elliptical galaxies. By any measure, BL~Lac
host galaxies look like completely normal ellipticals that are
somewhat brighter than average --- as far as the galaxy goes,
there is no evidence of the nuclear activity.

There are no systematic differences in the host galaxies of ``red'' and
``blue'' BL~Lac objects, once the obvious selection effects (on the BL~Lac
nuclei) are taken into account. Thus active nuclei with relativistic jets of
very different kinetic powers can live in very similar galaxies.
Their formation cannot be strongly affected by galaxy mass or morphology,
nor can their effect on the host galaxy be dramatic.

We confirm previous studies that BL~Lacs tend to lie in regions of
enhanced galaxy density, either groups or poor clusters, although the
small WFPC2 field-of-view limits the statistical significance of this
result. In some cases, however, the BL~Lac object appears truly
isolated, with no nearby companions or surrounding cluster galaxies,
to limits several magnitudes below the BL~Lac brightness.

We can rule out that a substantial fraction of BL~Lac objects at (apparently)
low redshift are actually high-redshift quasars microlensed by
intervening galaxies. Were the detected galaxies not hosts but
lensing galaxies,
in at least some cases the nuclei should be displaced from the center of
the galaxy. With our very large sample, we can say with high confidence
that this is not the case.

Contrary to previous studies, we do not find any correlation between nuclear
and host galaxy luminosities, such as might have been expected from the
trend of black hole mass with bulge mass seen in nearby ellipticals.
Although the observed nuclear properties of the BL~Lac objects are clearly
affected by beaming, correction for this effect makes no difference to
the lack of correlation.
Simulations show that selection effects, wherein bright nuclei can obscure
all but the brightest host galaxies, could contribute to spurious correlation.
The lack of an observed correlation for low-luminosity radio-loud AGN
implies a large scatter in Eddington ratio.

The unification of radio-loud AGN is strongly supported by our results.
This means that the properties of BL~Lac host galaxies and near environments
are basically universal to all low-luminosity radio-loud AGN. Just as
FR~I and FR~II radio galaxies span the full range of central engine power,
so do BL~Lacs represent the low-luminosity version of radio-loud quasars.
To understand fully trends in luminosity and/or redshift, samples of
BL~Lac objects and quasars should be combined, as we intend to do in
future work.

\acknowledgements
Support for this work was provided by NASA through grant number GO6363.01-95A
from the Space Telescope Science Institute, which is operated
by AURA, Inc., under NASA contract NAS~5-26555. MG acknowledges support from
the Hubble Fellowship Program through grant number HF 1071.01-94A, awarded by
the Space Telescope Science Institute. We thank Marcella Carollo,
Bob Hanisch, John Hutchings, Marek Kukula, and Susan Ridgway
for useful and stimulating discussions,
and Jochen Heidt, Kim McLeod, Ross McLure, and Chris O'Dea
for helpful comments on the manuscript.
RF thanks the {\it HST} visitor program for hospitality
during several visits to STScI. We thank G. Fasano for the use of the AIAP
package. This research made use of the NASA/IPAC Extragalactic Database (NED),
operated by the Jet Propulsion Laboratory, Caltech, under contract with NASA,
and of NASA's Astrophysics Data System Abstract Service (ADS).

\newpage


Cohesion of BaReH$_9$ and BaMnH$_9$: Density Functional Calculations and Prediction of (MnH$_9)^{2-}$ Salts

Abstract: Density functional calculations are used to calculate the structural and electronic properties of BaReH$_9$ and to analyze the bonding in this compound. The high coordination in BaReH$_9$ is due to bonding between Re 5$d$ states and states of $d$-like symmetry formed from combinations of H $s$ orbitals in the H$_9$ cage. This explains the structure of the material, its short bond lengths and other physical properties, such as the high band gap. We compare with results for hypothetical BaMnH$_9$, which we find to have similar bonding and cohesion to the Re compound. This suggests that it may be possible to synthesize (MnH$_9)^{2-}$ salts. Depending on the particular cation, such salts may have exceptionally high hydrogen contents, in excess of 10 weight %

\section{Introduction}

BaReH$_9$ is the prototypical member of a family of hydrides, based
on (ReH$_9)^{2-}$ and (TcH$_9)^{2-}$ structural units.
\cite{knox,gin1,ginsberg,stetson,stetson2,ryan}
These compounds are of interest from a fundamental point
of view
\cite{pauling}
because of the unusual coordination of the transition metal
atoms, which are surrounded by 9 hydrogen atoms, with relatively
short metal-H bond lengths and are in an high formal
valence state of 7.
These compounds are also of practical interest because
of the high hydrogen to metal ratio of 4.5 in BaReH$_9$.
Although no Mn based examples of these compounds have been
synthesized to date, the hypothetical 3$d$ analogue, MgMnH$_9$,
would have a hydrogen content in excess of 10 weight percent.
The purpose of this paper is to analyze the electronic structure of
BaReH$_9$ in order to understand its bonding and the prospects
for synthesis of Mn based analogues.
This follows a previous electronic structure calculation
by Orgaz and Gupta, \cite{orgaz} done using the X-ray crystal
structure of Ref. \onlinecite{stetson} and
the $X\alpha$ method. The present calculations were
done using a general potential self-consistent method, with
calculated H positions, and allow us to present a more
accurate electronic structure and detailed analysis of the bonding.

The crystal structure (spacegroup $P6_3/mmc$), which was
determined using X-ray diffraction by Stetson and
co-workers, \cite{stetson}
is depicted in Fig. \ref{structure}. It
consists of alternating triangular layers of Ba and Re stacked along the
$c$-axis, such that the Re are at the center of trigonal prisms
formed by the Ba. Each Re is coordinated by three H atoms
in its own plane (the H1) positions, and six other H atoms (the H2)
above and below the plane. These H2 atoms form a trigonal prism. Thus the
H atoms are arranged in tricapped trigonal prisms around the Re atoms.
The unit cell contains two (ReH$_9)^{2-}$ units, stacked so that
the orientation of these units alternates along the hexagonal $c$-axis.

\begin{figure}
\vskip 0.4 cm
\includegraphics[width=3.0in]{struct-small.eps}
\vskip 0.4 cm
\caption{(color online) Structure of BaReH$_9$, showing Ba and ReH$_9$
layers stacked along the $c$-axis (top) and the structure of
the hexagonal planes (bottom). The H1 and H2 sites are shown
by small light and dark red spheres, respectively, Ba by large
light green spheres, and Re by large dark blue spheres. The
atomic coordinates are those obtained by LDA structure
relaxation.}
\label{structure} \end{figure}

\section{Method}

The calculations reported here were done within the local density
approximation (LDA) using the general potential linearized augmented
planewave (LAPW) method. \cite{singh-book}
Local orbitals \cite{singh-lo}
were used for the semicore states (Re 5$s$ and 5$p$,
and Ba 5$s$ and 5$p$), and to relax the linearization of the $d$
bands. LAPW sphere radii of 2.0 $a_0$ and 1.1 $a_0$ were
used for the metal and H atoms, respectively.
We used well converged basis sets consisting of more than 3200
LAPW functions and local orbitals for the two formula unit
primitive cell. We tested larger basis sets, but found
no significant changes in the results.
We also tested various Brillouin zone samplings for the iteration
to self-consistency, but found that the results were already converged
at a sampling of 6 special {\bf k}-points in the irreducible wedge
of the hexagonal zone. No doubt this reflects the large band
gap insulating character of the compound.
The calculations for hypothetical BaMnH$_9$ were done
with similar parameters, except that the Mn and H radii were
1.8 $a_0$ and 1.05 $a_0$, respectively.

\section{Structure}

While the lattice parameters, metal positions,
and basic structure determined by
Stetson and co-workers (Ref. \onlinecite{stetson})
are undoubtedly correct,
it is difficult
to accurately determine H positions in a complex compound containing
heavy atoms without neutron scattering. Accordingly, we used the
experimental lattice parameters ($a$=5.287\AA, $c$=9.323\AA) here,
but determined the internal
coordinates of the H atoms by energy minimization within the LDA.
The calculated structural parameters are given in Table \ref{struct-tab}
along with those of hypothetical BaMnH$_9$, which is discussed
later.
As may be seen, the tricapped trigonal prisms are quite compact and
regular in
that the Re-H1 and Re-H2 bond lengths are short and similar, although
there is a compression along the $c$-axis, as may be seen from the
fact the the H2-H2 neighbor distance on the top of the trigonal
prism is longer than the H1-H2 distance between the prism and
cap hydrogen sites.
It is also notable that the H1-H2 distance
of 1.92\AA$~$ is very short compared with most metal hydrides, which
generally have H - H distances larger than 2.1 \AA. \cite{westlake1,westlake2}
These short H-H distances within the (ReH$_9)^{2-}$ units
suggest that direct H-H interactions may be important in forming
the electronic structure, which in fact is what we find.

\begin{table}[tbp]
\caption{Calculated structural parameters of  $P6_3/mmc$ BaReH$_9$
and hypothetical BaMnH$_9$, obtained
within the LDA. The lattice parameters were fixed at the
values $a$=5.287\AA, and $c$=9.323\AA, as reported
in Ref. \onlinecite{stetson} for BaReH$_9$.
For hypothetical BaMnH$_9$ we used $a$=5.067\AA, and $c$=8.883\AA$~$
(see text).
The coordinates of the
atoms are Ba (2$a$): (0,0,0), Re/Mn (2$c$): (1/3,2/3,1/4),
H1 (6$h$): ($x_1$,2$x_1$,1/4) and H2
(12$k$): ($x_2$,2$x_2$,$z_2$). Distances, $d$, are
the shortest bonds of a given type.}
\begin{tabular}{lcc}
\hline
parameter~~~~~ & ~~BaReH$_9$~~ &~~ BaMnH$_9$~~ \\
\hline
$x_1$ & 0.151 & 0.161 \\
$x_2$ & 0.466 & 0.460 \\
$z_2$ & 0.122 & 0.129 \\
$d$(Re/Mn-H1)   & 1.67\AA & 1.51\AA \\
$d$(Re/Mn-H2)   & 1.71\AA & 1.54\AA \\
$d$(H1-H2)      & 1.92\AA & 1.73\AA \\
$d$(H2-H2)      & 2.10\AA & 1.92\AA \\
$d$(H1-H1)$^*$  & 2.39\AA & 2.44\AA \\
$d$(Ba-H1)      & 2.71\AA & 2.63\AA \\
$d$(Ba-H2)      & 2.89\AA & 2.80\AA \\
\hline
$^*$between different (ReH$_9)^{2-}$ units.
\end{tabular}
\label{struct-tab}
\end{table}

\begin{figure}
\vskip 0.4 cm
\includegraphics[width=2.0in,angle=270]{band.eps}
\vskip 0.4 cm
\caption{Calculated band structure of hexagonal BaReH$_9$,
using the relaxed crystal structure.}
\label{bands} \end{figure}

\begin{figure}
\vskip 0.4 cm
\includegraphics[width=2.4in,angle=270]{dos.eps}
\includegraphics[width=2.4in,angle=270]{dos-h.eps}
\vskip 0.4 cm
\caption{(color online) Calculated electronic density of states
(DOS) of BaReH$_9$ on a per formula unit basis.
The top panel shows the total, Re $d$ and Ba contributions, as defined
by projections onto the LAPW spheres.
The bottom panel shows the projections of H1 $s$ and H2 $s$ character.
Note that for 1.1 $a_0$ H LAPW spheres approximately 2/3 of the
H $s$ charge for a H atom would lie outside the sphere.
}
\label{dos} \end{figure}

\section{Electronic Structure}

Our main electronic structure results are given in
Figs. \ref{bands} - \ref{integrated}.
The calculated band structure and electronic density of states (DOS)
are shown in Figs. \ref{bands} and \ref{dos}, respectively.
Fig. \ref{dos} also shows projections of the DOS onto the
LAPW spheres. Electron counts corresponding to integration
of the DOS and projections are shown in Fig. \ref{integrated}.
As in the calculation of Ref. \onlinecite{orgaz}, we obtain a wide
band gap insulator, but the bands and the structures and positions
of the features in the DOS are significantly different.
We obtain a direct band gap of 3.58 eV at $\Gamma$, which may be
an underestimate, as is common in LDA calculations.
In any case, this large band gap is consistent with
the observed transparent nature of the material, and also
indicates the nature of the band formation. In particular,
it is an extremely large value for a crystal field gap, and also
is not consistent with an ionic gap, since in a scenario
with H$^-$ and Re$^{7+}$ ions, one would expect low
lying Re $d$ levels and therefore a small gap
due to the high Re valence. Such an ionic scenario is also unlikely
considering the short H-H bond lengths in the (ReH$_9)^{2-}$ units.

\begin{figure}
\vskip 0.4 cm
\includegraphics[width=2.4in,angle=270]{int.eps}
\includegraphics[width=2.4in,angle=270]{int-h.eps}
\vskip 0.4 cm
\caption{(color online) Integrated
DOS of BaReH$_9$ on a per formula unit basis.
The top panel shows the total, Re $d$ and Ba contributions, as in Fig.
\ref{dos}.
The bottom panel shows the projections of H1 $s$ and H2 $s$ character.
The integrations are on a per formula unit basis.
}
\label{integrated} \end{figure}

Examining the band structure in greater detail, we note that there
are three manifolds of states. These are a narrow set of two bands
(one band per formula unit) extending from -6.2 eV to -5 eV (relative
to the valence band maximum), a manifold of 16 bands extending
from -3.3 eV to the valence band maximum and a broad set of conduction
bands. While valence bands are often less dispersive than
conduction bands in semiconductors, this electron-hole asymmetry is
particularly noticeable in the BaReH$_9$. As may be seen from
the projections, of the DOS there is no appreciable Ba
contribution to the valence bands, and therefore, as expected, Ba
occurs as Ba$^{2+}$. There is, however, Ba $d$ character in the upper
conduction bands starting at $\sim$ 7 eV, and the Ba $f$ resonance is visible
from $\sim$ 10 - 11 eV.

With the exception of the lowest, split-off peak in the DOS, which
is H $s$ in character, the remainder of the DOS shows a mixture
of Re $d$ and H $s$ character, demonstrating covalency between Re and H.
This is distinct from what is seen in NaAlH$_4$, for example.
\cite{aguayo}

The projected densities of states and electron counts are based
on integration within the LAPW spheres. In the LAPW method these
spheres are constrained to be non-overlapping, and in addition
they must be chosen so that core states are contained within them.
These constraints require the use of rather small spheres at the
H sites. Here we used $r_{\rm H}$ = 1.1 $a_0$. This is significant
for the interpretation of the projected DOS, because most of the
H $s$ charge lies outside this radius. In particular,
for a neutral non-spin-polarized H atom in the LDA only 0.329 electrons
are contained within 1.1 $a_0$,
and 0.305 electrons within 1.05 $a_0$. This factor needs to be kept in
mind when analyzing the DOS. The advantage of the small H LAPW sphere
is that the H wavefunction near the nucleus is much less sensitive to
the environment, which means that charge within this size sphere
is almost entirely due to the H $s$ states, with very little
contamination by charge entering from neighboring atoms. We
verified this by integrating the $p$ projection in the H spheres,
which would come from tails of neighbors,
and found less than 0.02 electrons, supporting this conclusion.
Thus, integration of the H $s$ projection can provide a useful
measure of the H valence, when normalized by the charge of a neutral H.
The charge obtained in this way corresponds to H$^{0.5-}$ for
both the H1 and H2 sites. This analysis cannot be easily applied
to the conduction bands because there is no  plateau at the top of
the Re $d$ derived DOS, presumably because of contributions from the
free electron like states, which would provide H $2s$ and other contributions.
Nonetheless, by this measure, 2 electrons per H would be included
by integrating up to 9 eV, while integrating to 10 eV, where
the Re $d$ character becomes small, yields 2.1 electrons.
In any case, these numbers indicate substantial covalency between
Re and H. This is also consistent with the large bandgap, which
comes from the bonding - antibonding splitting associated with the
covalency.

\section{Bonding}

To understand the bonding, we examined the symmetry of the H $s$ states
when expanded about the Re site, using the projections of the H tails
into the Re sphere. Based on this, the lowest split-off peak
containing 2 electrons comes
from the symmetric combination of the H $s$ states, and is $s$-like
at the Re site. The second peak, which contains 16 electrons per
(ReH$_9)^{2-}$ unit has both $p$-like and $d$-like combinations of
the H $s$ states. Considering the short H-H distance in the H$_9$
cage, and the simple electronic structure of H, the electronic
states of a H shell should be ordered according to the number of nodes,
yielding, in order, a fully symmetric state ($s$-like), three states with
one node ($p$-like) and five states with two nodes ($d$-like).
In fact, calculations for free Be$_9$ and Mg$_9$ clusters,
\cite{kawai,kumar}
which
have 18 electrons and single $s$ orbital ions, show that these
clusters are exceptionally stable and have a tricapped trigonal prismatic
ground state structure, similar to 18 electron (ReH$_9)^{2-}$.
Furthermore, the electronic structure \cite{kumar}
shows clearly separated groups of states, which are, in order,
s-like, p-like and d-like. The energy splittings of the states
within the groups due to the non-spherical, $D_{3h}$, symmetry
are small compared to the splittings between the groups.
This strong shell structure of the orbitals is associated with
the stability of these clusters and their shape, as discussed
by Reimmann and co-workers, \cite{reimann} and may be expected
to be a general tendency in 18 electron, 9 site clusters with only
$s$-orbitals.

This forms the basis for understanding the electronic structure and
bonding of BaReH$_9$,
as schematically illustrated in Fig. \ref{levels}.
Specifically,
the d-like states of the H$_9$
shell have angular character about the Re site that favors bonding
with the Re d-states. This yields five bonding states, which are
$\sigma$ bonding combinations in the radial direction, and five
antibonding states, which have nodes between the Re and the H shell.
This provides an explanation for
the asymmetry of the valence and conduction bands,
since the antibonding combinations will be directed away from the
(ReH$_9)^{2-}$ clusters, leading to greater dispersion than for the
bonding combinations.

The strong bonding of the $d$-like combinations of H $s$ states on the H$_9$
shell and Re $d$ states leads to a large band gap, and shifts the bonding
states so that they are in the same energy region as the $p$-like
combinations of shell states, which are non-bonding with Re $d$ orbitals.
Thus, the 16 electron main peak in the valence DOS comes from
10 electrons in bonding Re $d$ - shell $d$ states, and 6 electrons
in shell $p$-like states. Considering the large bonding - antibonding splitting
of $\sim$ 8 eV as measured by the energy separation of the centers
of the valence and conduction Re $d$ DOS, there should be approximately
equal mixtures of Re $d$ and H $s$ character in both the bonding
and antibonding states. This would yield an electron count of
13 H 1$s$ electrons in the valence bands (including the $s$-like and $p$-like
shell states) and therefore a hydrogen valence of H$^{0.44-}$, in reasonable
accord with the estimate obtained from integrating the projected
DOS.

\begin{figure}
\vskip 0.4 cm
\includegraphics[width=3.5in,angle=0]{levels.eps}
\vskip 0.4 cm
\caption{Schematic depiction of the formation of the BaReH$_9$
electronic structure from the levels of a tricapped trigonal
prismatic H$_9$ cluster (see Ref. \onlinecite{kumar})
and
the Re $d$ orbitals. The levels of the cluster are labeled according
using $D_{3h}$ representations, and also $s$, $p$ and $d$ according to
their approximate angular momentum character when expanded around
the center of the cluster. Note that all the H states shown are derived from
different linear combinations of H 1$s$ orbitals.
}
\label{levels} \end{figure}

\section{Hypothetical Manganese Based Compounds}

The bonding mechanism discussed above is quite general in principle,
depending on a proper electron count
and the formation of strong bonds between $d$ orbitals and
appropriate combinations of H $s$ orbitals.
The $d$ orbitals of 3$d$ elements are much less extended than
those of 4$d$ and 5$d$ elements.
Thus, as is well known, 3$d$ compounds typically have smaller band widths
and reduced covalency with ligands than 4$d$ and 5$d$ analogues.
This might be used to provide
an explanation for why salts containing (MnH$_9)^{2-}$
have not been reported. Nonetheless, considering that the bonding in
BaReH$_9$ appears to be very strong, based on the large bonding - antibonding
splitting and the fact that the phase competes with the very stable hydride
BaH$_2$, it is of interest to study the Mn based analogue. We note
that the binding energy of MgH$_2$ is less than half of that of BaH$_2$,
and that the binding energy of BeH$_2$ is even lower, and so even with
reduced stability of the anionic cluster, salts may be formed with
appropriate cations.

In order to investigate the electronic structure of hypothetical
BaMnH$_9$, we reduced the $a$ and $c$ lattice parameters of
the Re compound based on the difference in the covalent radii
of Re and Mn (0.11\AA), to obtain $a$=5.067\AA, and $c$=8.883\AA.
This corresponds to a volume reduction of 12.5\%
relative to BaReH$_9$. Using these lattice parameters, we
relaxed the H positions in the unit cell. The resulting
H coordinates and bond lengths are give in Table \ref{struct-tab}.
Remarkably, the relaxed structure maintains a cage of nine
H atoms at approximately the same distance from Mn, and with
short H-H distances.

\begin{figure}
\vskip 0.4 cm
\includegraphics[width=2.4in,angle=270]{dos-mn.eps}
\includegraphics[width=2.4in,angle=270]{dos-h-mn.eps}
\vskip 0.4 cm
\caption{(color online) Calculated electronic density of states
(DOS) of hypothetical BaMnH$_9$ on a per formula unit basis,
as in Fig. \ref{dos}. Note that the H LAPW sphere radius
for BaMnH$_9$ was 1.05 $a_0$.
}
\label{dos-mn} \end{figure}

As may be seen from the DOS and projections
(Fig. \ref{dos-mn}), the electronic structure is qualitatively
similar to BaReH$_9$, although the detailed peak structure of the
main valence band DOS differs and the split-off H state moves
to lower energy. This larger splitting between the shell $s$-like
state and the shell $p$-like states is what is expected based on
the smaller H-H distances in the more compact H$_9$ cage of the
Mn compound, since this will yield larger hopping matrix elements
and therefore stronger shell structure.
As in the Re compound, integration of the H $s$ projection of the
DOS, normalized by the amount of charge within a 1.05 $a_0$ sphere
for a neutral non-spin-polarized H atom in the LDA, yields a
hydrogen valence, H$^{0.5-}$. Furthermore, although the band gap
of BaMnH$_9$ is reduced to 3.0 eV, the bonding - antibonding splitting
as measured by the centers of the Mn projection of the DOS is
nearly the same as in the Re compound.

The implication of the above results is that salts containing
(MnH$_9)^{2-}$ may be realizable. Synthesis may, however, be
difficult because of competition from stable hydride phases.
As noted, MgH$_2$ is less stable than BaH$_2$. However, MgH$_2$ still
has a binding energy of 75 KJ/mol H$_2$, and furthermore Mg
forms other stable hydrides with Mn, such as Mg$_3$MnH$_7$. \cite{bortz,gupta}

The LDA often overbinds solids. Nonetheless, it is of interest to
compare the heats of formation of BaReH$_9$ and hypothetical
BaMnH$_9$. We calculated these values by subtracting calculated total
energies of the elemental metals and the H$_2$ molecule from the
total energy of the compounds. For Mn, we used the fcc structure
and made a correction of 5 mRy per atom, which comes from energy
difference between the fcc structure and the ground state non-collinear
antiferromagnetic $\alpha$-Mn, as obtained by Hobbs and Hafner.
\cite{hobbs}
We made no correction for zero point energy.

With the LDA value of the energy of H$_2$, as obtained
in a supercell calculation, $E({\rm H}_2)$ = -2.294 Ry,
the resulting formation energies are -99 KJ/mol H$_2$ for
BaReH$_9$ and -86 KJ/mol H$_2$ for hypothetical BaMnH$_9$.
The main uncertainties are (1) the crystal structure of BaMnH$_9$,
which if different from that assumed would lead to a higher
stability, (2) the neglected zero point motion, which typically
affects the cohesive energies by less than 20 KJ/mol H$_2$ and
here would likely destabilize the compounds as the H-metal and
H-H bonds are short, suggesting stiff phonons, and
(3) the use of the LDA, which based on calculations for
other metal hydrides,
\cite{smithson,halilov,miwa}
overestimates stability typically
by $\sim$ 20 KJ/ mol H$_2$, though this varies from compound to compound.
In any case, based on these results, we conclude that BaReH$_9$ and
hypothetical BaMnH$_9$ have similar stabilities, assuming that the
decomposition is to similar products, {\em e.g.} BaH$_2$ and Re/Mn metal.

\section{summary and conclusions}

Electronic structure calculations show that the high coordination in
BaReH$_9$ is due to bonding between Re 5$d$ states
and $d$-like states formed from combinations
of H $s$ orbitals in the H$_9$ cage. This explains the structure of the
material, its short bond lengths and other physical properties, such
as the high band gap. Similar bonding is found in hypothetical BaMnH$_9$
and in fact both compounds are found to have similar cohesive energies.
While BaMnH$_9$ may not be the most favorable target due to the
stability of BaH$_2$ as a competing phase, our results suggest that
synthesis of (MnH$_9)^{2-}$ salts may be possible and should be
attempted. We note that the chemistry of Mn VII is different
than that of Re VII. For example, KMnO$_4$ is highly oxidizing
in aqueous solution, while
KReO$_4$ is not. Thus, while it may be possible to synthesize (MnH$_9)^{2-}$
salts, the solvents and synthesis path needed may be different from that
used in the (ReH$_9)^{2-}$ salts.
Light element salts of (MnH$_9)^{2-}$ would have very high hydrogen
contents, in excess of 10 weight \% in the cases of Li, Mg or Be.

\acknowledgments

DJS thanks the University of Paris-Sud for their hospitality,
which made this work possible.
Crystal structures were plotted using the XCrySDEN program. \cite{xcrysden}
Work at Oak Ridge National Laboratory is supported by the U.S. Department
of Energy.
We thank the Institut du Developpement et des Ressources en
Informatique Scientifique (IDRIS) for a grant of computer time.