Topological crystalline insulator states in Pb(1-x)Sn(x)Se

Abstract: Topological insulators are a novel class of quantum materials in which time-reversal symmetry, relativistic (spin-orbit) effects and an inverted band structure result in electronic metallic states on the surfaces of bulk crystals. These helical states exhibit a Dirac-like energy dispersion across the bulk bandgap, and they are topologically protected. Recent theoretical proposals have suggested the existence of topological crystalline insulators, a novel class of topological insulators in which crystalline symmetry replaces the role of time-reversal symmetry in topological protection [1,2]. In this study, we show that the narrow-gap semiconductor Pb(1-x)Sn(x)Se is a topological crystalline insulator for x=0.23. Temperature-dependent magnetotransport measurements and angle-resolved photoelectron spectroscopy demonstrate that the material undergoes a temperature-driven topological phase transition from a trivial insulator to a topological crystalline insulator. These experimental findings add a new class to the family of topological insulators. We expect these results to be the beginning of both a considerable body of additional research on topological crystalline insulators as well as detailed studies of topological phase transitions.

\section*{Methods}
The preparation of the highest quality materials was essential for the studies presented here. Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace single crystals exhibiting a rock-salt structure were grown by the self-selecting vapour growth (SSVG) method under near-equilibrium thermodynamic conditions~\cite{Szczerbakow-ProgCrystalGrowth-2005, Szczerbakow-JCrystalGrowth-1994}, where the specific profile of the temperature field ensures contact between the growing crystal and the crystal's own source material only. This method offers unique advantages in the growth of compositionally uniform solid solutions of high crystal quality and natural (001) facets (Supplementary Information Fig.~\ref{fig:S1}). For our photoemission and magnetotransport studies, parallelepiped samples were cleaved from the bulk crystals along the (001) crystal planes (see the AFM microscopy pictures in Supplementary Information Fig.~\ref{fig:S2}). All samples were prepared such that they possessed $n$-type conductivity, which was performed to achieve both occupied conduction-band and TCI states.

Angle-resolved photoemission spectroscopy studies were carried out on the (001) surface of cleaved single crystals. These studies were performed using linearly polarised light having photon energies between $10.5$~eV and $40$~eV. Temperature-dependent photoemission experiments were performed at the BALTAZAR laser-ARPES facility~\cite{Berntsen-RevSciInstrum-2011} using $10.5$~eV photons. The total energy and crystal-momentum resolution for these measurements were $5$~meV and $0.008$~\AA$^{-1}$, respectively. The single crystal used in the temperature study was cleaved at room temperature under ultra-high-vacuum (UHV) conditions. The photoemission spectra were acquired at selected temperatures in the range of $9$~K to $300$~K. Complementary measurements at higher photon energies and temperatures as low as $110$~K have been conducted using the I4 beam line~\cite{Jensen-NIMA-1997} at the MAX III synchrotron at MAX-lab, Lund University.

\section*{Acknowledgments}
We would like to acknowledge V.~Domukhovski and A.~Reszka for structural and chemical composition analyses of the crystals, J.~Adell for his help during our beamtime at MAX-lab, S.~Safai for her help in numerical calculations, and P.~Kacman for critical reading of the manuscript. In Poland, this work was supported by the European Commission Network SemiSpinNet (PITN-GA-2008-215368) and by the European Regional Development Fund through the Innovative Economy grant (POIG.01.01.02-00-108/09). In Sweden, this work was made possible through support from the Knut and Alice Wallenberg Foundation and the Swedish Research Council.

\input{manuscript-arXiv-20130810.bbl}

\newpage

\renewcommand{\thesection}{S\arabic{section}}

\section{Crystal growth, structural and chemical characterization}
\label{sec:S1}

The method of self-selecting vapour growth (SSVG) was applied to manufacture
Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace monocrystals (Supplementary Fig.~\ref{fig:S1}) of excellent crystal structure as well as
chemical homogeneity. The SSVG method is particularly suitable for optimisation of the
latter because the crystals grow in a near-equilibrium regime. Moreover, the temperature
profile set in the technological furnace ensures growth of a crystal in exclusive contact with
its source material. In the case of the considered alloy, the method allows controlling
deviations from stoichiometry to achieve (without heterodopants) the required $n$-type
conductivity with a moderate electron concentration. The SSVG method was successfully
employed by one of us (A.S.) to grow monocrystals of various II-VI and IV-VI
semiconductor compounds and their solid solutions~\cite{Szczerbakow-ProgCrystalGrowth-2005, Szczerbakow-JCrystalGrowth-1994}.

X-ray diffractometry (XRD) analysis of the crystal structure of Pb\texorpdfstring{$_{1-x}$}{\ifpdfstringunicode{\unichar{"2081}\unichar{"208B}\unichar{"2093}}{1-x}}Sn\texorpdfstring{$_{x}$}{\ifpdfstringunicode{\unichar{"2093}}{x}}Se\xspace was done at room
temperature for both as-grown and freshly cleaved (001) crystal surfaces as well as for
powdered polycrystals. It revealed the rock-salt single-crystal phase with a lattice parameter $a=(6.0979 \pm 0.0006)$~\AA{}
corresponding to a Sn content of $x = 0.228 \pm 0.005$~\cite{Szczerbakow-JCrystalGrowth-1994}. XRD rocking-curve measurements of the (002) diffraction peaks yield a full width at half maximum
parameter $\Delta\omega = (70~\mathrm{to}~150)$~arcsec for various (001) surfaces. The chemical composition of
the crystals was determined by energy-dispersive X-ray spectroscopy (EDX) using a scanning
electron microscope SEM Hitachi SU-70 equipped with a Thermo Scientific NSS7
microanalysis system. Probing different places (over an area of $1~\mathrm{mm}^2$) of a sample freshly
cleaved (in air) we found a Sn content of $x = 0.234 \pm 0.005$.

\begin{figure}[hb]
\centering
\includegraphics[width=.9\columnwidth]{S1.eps}
\caption{{\bf As-grown Pb\texorpdfstring{$_{\boldsymbol{0.77}}$}{0.77}Sn\texorpdfstring{$_{\boldsymbol{0.23}}$}{0.23}Se\xspace bulk monocrystal.} The mirror-like planes at the front part of the crystal are natural (001) facets of a rock-salt-type crystal.}
\label{fig:S1}
\end{figure}

The surface morphology of Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace crystals was studied using a Dimension Icon
Atomic Force Microscope (AFM). The measurements were carried out for (001) surfaces
prepared by cleaving the crystals in air at room temperature. Supplementary Fig.~\ref{fig:S2} presents an
AFM image for an area of $1\times 1~\mu\mathrm{m}^2$ and the section analysis taken along the indicated line.
The root-mean-square surface roughness parameter for this area equals $6.8$~\AA{} (i.e., about just
one lattice parameter). The freshly cleaved (001) surface consists of large atomically flat
regions with small vertical steps with heights of a few monolayers. Time-dependent AFM
measurements of a cleaved surface exposed to air showed substantial surface contaminations
appearing after $1.5$~hours.

\begin{figure}[hb]
\centering
\includegraphics[width=.9\columnwidth]{S2.eps}
\caption{{\bf Morphology of the (001) surface of a Pb\texorpdfstring{$_{\boldsymbol{0.77}}$}{0.77}Sn\texorpdfstring{$_{\boldsymbol{0.23}}$}{0.23}Se\xspace crystal.}
Atomic force microscopy (AFM) image of a (001)
surface of a bulk crystal freshly cleaved in air at room
temperature. The height differences between the points
marked with green and red arrows equal $6.8$~\AA{} and $11.5$~\AA{},
respectively. The macroscopically large region between
the black arrows is atomically flat.}
\label{fig:S2}
\end{figure}

\section{Angle-resolved photoemission measurements at various photon energies}
\label{sec:S2}

\begin{figure*}[ht]
\centering
\includegraphics[width=.9\textwidth]{S3.eps}
\caption{Selected constant-energy surfaces of the Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace surface state for the photon energies $17.5$~eV, $25$~eV and $40$~eV.}
\label{fig:S3}
\end{figure*}

Supplementary Fig.~\ref{fig:S3} shows three sets of selected constant-energy surfaces for Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace,
in the vicinity of the $\overline{X}$ point, measured at a temperature of $110$~K using photon energies of
$17.5$~eV, $25$~eV, and $40$~eV with corresponding overall energy resolution of $25$~meV to
$40$~meV. The data were acquired at the I4 beam line of the MAX III synchrotron at MAX-lab,
Lund University, Sweden. As seen in the figure, the spectral weight is shifted when changing
the photon energy. For the measurements at $25$~eV and $40$~eV the part of the surface state
located in the second Brillouin zone (BZ) is the most intense. At a photon energy of $17.5$~eV
the most intense part shifts from the first BZ for high binding energies (below the Dirac
points) to the second zone for lower binding energies (above the Dirac points). Apart from
these intensity variations the position of the Dirac-like feature with respect to the surface
Brillouin zone is independent of photon energy (within the limits of experimental accuracy)
suggesting that the observed feature is indeed a surface electronic state.

\begin{figure*}
\centering
\includegraphics[width=.9\textwidth]{S4.eps}
\caption{Energy-momentum slices across the surface Dirac cones in Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace taken at the photon energies $10.5$~eV, $17.5$~eV, $25$~eV and $40$~eV.}
\label{fig:S4}
\end{figure*}

Comparing energy-momentum slices from the data sets taken at different photon energies also
confirms that the feature has virtually no photon-energy dependence. Supplementary Fig.~\ref{fig:S4}
displays such energy-momentum slices through the Dirac point located in the second surface BZ
(analogous to Fig.~\ref{fig:3}a in the main text) for the photon energies $10.5$~eV, $17.5$~eV, $25$~eV and
$40$~eV. The cuts are parallel to the $\overline{M}$--$\overline{X}$--$\overline{M}$ direction. The $10.5$~eV data are taken at the
BALTAZAR laser ARPES facility at a temperature of $100$~K and the $17.5$~eV, $25$~eV and $40$~eV data are acquired at the I4 beam line of MAX III at a temperature of $110$~K. Disregarding
for the moment the intensity differences, a clear X-shaped state is seen for all photon
energies. However, for the $17.5$~eV data the Dirac point is shifted about $55$~meV towards
lower binding energy as compared to the $10.5$~eV measurement. A slight further shift towards
lower binding energy is observed for the other measurements. The $10.5$~eV data were
recorded within $12$~hours after cleaving whereas the $17.5$~eV, $25$~eV and $40$~eV data were
recorded at $80$~hours, $85$~hours and $100$~hours post cleaving, respectively. Thus, the shift in
binding energy is likely related to the aging of the surface due to residual gas adsorption.
Similar observations have previously been made for topological insulators, e.g. Bi$_2$Se$_3$, for
which the Dirac point shifts towards higher binding energy upon gas adsorption and these
shifts can amount to several hundred millielectronvolt~\cite{King-PhysRevLett-2011}.
In the present case, the shift is towards lower binding energy and seems
at first glance to be of smaller magnitude. In the Bi$_2$Se$_3$ case, the adsorbtion of residual gas is
accompanied by the appearance of quantum well states that display a large Rashba split~\cite{King-PhysRevLett-2011}. For Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace, neither the laser-based
data nor the synchrotron data show any clear sign of quantum well states or other
additional surface-related states.

Employing a free-electron model with an inner potential of $15$~eV deduced from normal
emission data, for states close to the Fermi energy at $\overline{X}$ the photon energies $10.5$~eV,
$17.5$~eV, $25$~eV and $40$~eV correspond to $k_z$ values $2.2$, $2.5$, $2.9$ and $3.4$ (in units of $2\pi/a$, where
$a = 6.075$~\AA{}), respectively. However, it should be kept in mind that the $k_z$ resolution in the
studied photon-energy range is very limited as a consequence of the high surface sensitivity.
Nevertheless, different $k_z$ values are probed and bulk-related dispersive features are expected
to differ. The lack of dispersive differences for the Dirac state between the data sets once
more confirms the surface nature of the observed state.

Looking in more detail at the intensity ``inside'' the Dirac state it is clear from Fig.~\ref{fig:3} in the
main text that this intensity in the $10.5$~eV data is partially originating from the ``second
cone'', which has its Dirac point located on the other side of $\overline{X}$. This is not the case for the
higher photon energy measurements since the shift of the chemical potential causes this
contribution to lie above the Fermi level. However, from Supplementary Fig.~\ref{fig:S4} it is apparent
that bulk features contribute to the intensity ``inside'' the cone since it changes as a function of
photon energy. The Fermi level intensity varies with photon energy as does the intensity close
to $k_y = 0$ for high binding energies. The high intensity at the Fermi level in the $17.5$~eV data
and at higher binding energy in the $25$~eV data might indicate that the corresponding points in
reciprocal space are close to a conduction-band minimum and a valence-band maximum,
respectively, although final-state effects and the limited $k_z$ resolution hinder a more accurate
determination.

Returning to the Dirac state, looking at the intensity variations between the different photon
energies one observes that the expected increase in intensity at the Dirac point (c.f. Fig.~\ref{fig:3}a in
the main text), due to the band crossing, is not present. This indicates that the gap is not yet
fully closed and that the transition temperature is slightly below $100$~K. Some of the data seem to
indicate a closed gap but caution is called for since intensity variations as a function of photon
energy might mask a small gap. Such variations are clearly present since at $10.5$~eV the main
intensity is in the part above the Dirac point whereas the opposite is true for $40$~eV. Also, the
intensity at the Dirac point is varying across the photon energies used. The difference in
energy resolution between the laser ARPES and synchrotron measurements provides another
potential explanation for the apparently closed gap in parts of the data.

\section{Electrical properties}
\label{sec:S3}

\begin{figure*}
\centering
\includegraphics[width=.3\textwidth]{S5a.eps}\hfill
\includegraphics[width=.3\textwidth]{S5b.eps}\hfill
\includegraphics[width=.3\textwidth]{S5c.eps}
\caption{{\bf Electrical characterisation of Pb\texorpdfstring{$_{\boldsymbol{0.77}}$}{0.77}Sn\texorpdfstring{$_{\boldsymbol{0.23}}$}{0.23}Se\xspace crystals.}
Temperature dependence of the resistivity (a), the Hall electron concentration (b),
and the Hall electron mobility (c) determined from low-field classical transport measurements.}
\label{fig:S5}
\end{figure*}

The samples for electron transport measurements were prepared from as-grown bulk crystals
by cleaving along three mutually orthogonal $(001)$ planes to get standard parallelepiped Hall-bars
of typical dimensions $2\times 2\times 7~\mathrm{mm}^3$. Directly after cleaving, six contacts were made by
soldering with indium. The low-field ($B = 0.5$~T) Hall effect and resistivity measurements
were performed in a continuous-flow helium cryostat with a DC current of $50$~mA applied
along the $[100]$ crystallographic direction. The electron transport parameters (resistivity, Hall
electron concentration and Hall mobility) presented in Supplementary Fig.~\ref{fig:S5} were obtained
considering one transport channel only in the analysis, as usually done for characterisation
purposes.

\section{High-field magnetotransport measurements}
\label{sec:S4}

\begin{figure*}
\centering
\includegraphics[width=.3\textwidth]{S6a.eps}\hfill
\includegraphics[width=.3\textwidth]{S6b.eps}\hfill
\includegraphics[width=.3\textwidth]{S6c.eps}
\caption{{\bf High-field magnetotransport measurements}
The high-field longitudinal resistivity $\rho_{xx}$ (a) and
transverse resistivity $\rho_{xy}$ (b) of the Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace sample at several temperatures. (c) presents
the $\rho_{xy}$ data after subtraction of the linear term.}
\label{fig:S6}
\end{figure*}

The high-field measurements of the resistivity tensor components $\rho_{xx}$ and $\rho_{xy}$ were carried out
in a liquid-helium cryostat equipped with a $B = 13$~T superconducting solenoid. The
magnetoresistance presented in Supplementary Fig.~\ref{fig:S6}a is of usual orbital origin with
characteristic two-carrier conductivity contributions. Supplementary Fig.~\ref{fig:S6}b displays the as-measured
experimental data for the transverse $\rho_{xy}$ tensor component, which after the
subtraction of linear terms were transformed to the $\Delta\rho_{xy}$ curves shown in Supplementary Fig.~\ref{fig:S6}c.

The experimentally observed magnetic field nonlinearities in the $\rho_{xy}$ tensor component and
the non-quadratic magnetoresistance indicate the presence of at least two parallel conduction
channels. The dominant conduction channel can be identified as a bulk conduction band
contribution. The identification of the other electronic channel needs further magnetotransport
studies, in particular for samples with lower bulk electron concentration.

\begin{figure}
\centering
\includegraphics[width=.9\columnwidth]{S7.eps}
\caption{{\bf Shubnikov-de Haas oscillations}
The high-field magnetoresistance data recorded for a Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace sample at the temperature
$T = 4.2$~K in tilted magnetic fields. The tilting angle $\beta$ is the angle between the
external magnetic field and the normal of the $(001)$ sample plane.}
\label{fig:S7}
\end{figure}

The good electrical quality of the Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace crystals reflected in the high electron mobility
permits the experimental observation at very low temperatures of quantum oscillations of the
magnetoresistance (the Shubnikov-de Haas effect, SdH). The SdH oscillations in $n$-Pb\texorpdfstring{$_{0.77}$}{0.77}Sn\texorpdfstring{$_{0.23}$}{0.23}Se\xspace displayed in
Supplementary Fig.~\ref{fig:S7} are of bulk (3D) origin. This was confirmed
in a tilted-field experiment, where the magnetic-field positions of the oscillations did not
change with the tilting angle. The electron concentration (in a single energy valley) obtained
from the analysis of the period of the SdH oscillations ($n_{\mathrm{SdH}} = 8 \times 10^{17}~\mathrm{cm}^{-3}$)
perfectly corresponds to the total concentration obtained from the low-field Hall effect
($n_{\mathrm{H}} = 3.2 \times 10^{18}~\mathrm{cm}^{-3}$). The factor of $4$ difference is due to the $L$ valley
4-fold degeneracy in Pb\texorpdfstring{$_{1-x}$}{\ifpdfstringunicode{\unichar{"2081}\unichar{"208B}\unichar{"2093}}{1-x}}Sn\texorpdfstring{$_{x}$}{\ifpdfstringunicode{\unichar{"2093}}{x}}Se\xspace. Since the bulk contribution dominates the total conductivity at low temperatures,
the SdH oscillations due to topological surface states are not visible in our crystals. The very
small contribution of surface states to the conductivity at low temperatures is typically
observed in other well established topological insulators~\cite{Qu-Science-2010, Taskin-PhysRevLett-2011, Analytis-NatPhys-2010, Ren-PhysRevB-2010}.

\section{Tight-binding (TBA) calculations: the crystal-slab-thickness dependence}
\label{sec:S5}

In the case of the critical Sn composition $x_{\mathrm{c}} = 0.381$ the calculated fundamental bulk bandgap
at the $L$ point of the Brillouin zone is closed. However, in our calculations performed for
crystal slabs of finite thickness we observe a small gap at the $\overline{X}$ point. This effect is due to
the electron confinement as verified by the calculation of the bandgap dependence on the
crystal-slab thickness. We have calculated the bandgaps at various points of the surface
Brillouin zone for Pb\texorpdfstring{$_{1-x}$}{\ifpdfstringunicode{\unichar{"2081}\unichar{"208B}\unichar{"2093}}{1-x}}Sn\texorpdfstring{$_{x}$}{\ifpdfstringunicode{\unichar{"2093}}{x}}Te\xspace crystal slabs of varying thickness up to about $160$~nm.
Supplementary Fig.~\ref{fig:S8} presents the results in two cases. For the fundamental gap at the $\overline{X}$
point in the crystal with critical Sn content ($x = x_{\mathrm{c}}$) the finite-size effect is clearly observed.
The bandgap is inversely proportional to the slab thickness (blue line), as expected, and
decreases asymptotically to zero with increasing slab thickness. In contrast, the surface states
bandgap in the $\overline{X}$--$\overline{M}$ direction in the crystal with a Sn content $x = 0.6$ (red line) remains
constant for thick slabs thus proving its intrinsic origin. In this case the gap is not determined
by the electron confinement. Its apparent dependence on the thickness for very thin slabs
stems from the surface wave functions overlap.

\begin{figure}[hb]
\centering
\includegraphics[width=.9\columnwidth]{S8.eps}
\caption{{\bf TBA calculations of the crystal-slab-thickness dependence of energy gaps in Pb\texorpdfstring{$_{\boldsymbol{1-x}}$}{\ifpdfstringunicode{\unichar{"2081}\unichar{"208B}\unichar{"2093}}{1-x}}Sn\texorpdfstring{$_{\boldsymbol{x}}$}{\ifpdfstringunicode{\unichar{"2093}}{x}}Te\xspace.}
The blue line shows the bandgap at the $\overline{X}$ point in Pb\texorpdfstring{$_{1-x}$}{\ifpdfstringunicode{\unichar{"2081}\unichar{"208B}\unichar{"2093}}{1-x}}Sn\texorpdfstring{$_{x}$}{\ifpdfstringunicode{\unichar{"2093}}{x}}Te\xspace with a Sn content $x = x_{\mathrm{c}} = 0.381$, whereas the red line
displays the bandgap on the $\overline{X}$--$\overline{M}$ line for a Sn content $x = 0.6$. For the Sn content range investigated,
the thickness of one monolayer is about $3.2$~\AA{}.}
\label{fig:S8}
\end{figure}

\end{document}


On the origin of very high energy gamma-rays from radio galaxies

Abstract: Radio galaxies have emerged as a new gamma-ray emitting source class on the extragalactic sky. With their jets misaligned, i.e. not directly pointing towards us, they offer a unique tool to probe some of the fundamental (and otherwise hidden) non-thermal processes in AGN. This contribution briefly summarizes the observed characteristics of the four radio galaxies detected so far at very high energies (VHE). Given its prominence, particular attention is given to the origin of the variable VHE emission in M87. We discuss some of the theoretical progress achieved for this source within recent years highlighting, amongst others, the relevance of magnetospheric particle acceleration and emission models.

\section{Introduction}
The extragalactic sky as observed at very high energies (VHE; $>100$ GeV) is
highly dominated by Active Galactic Nuclei (AGNs) of the blazar-type, i.e., AGNs
with relativistic jets pointing almost directly towards us. Flux enhancement
of the intrinsic non-thermal jet emission by relativistic Doppler boosting effects
then naturally favours their detection at VHE. On the other hand, misaligned (i.e.,
non-blazar) AGNs, characterized by jets substantially inclined with respect to
the observer, could allow us to probe some of the fundamental non-thermal
VHE processes in AGN that are otherwise swamped by boosted jet emission.
Nearby radio galaxies (RGs) are especially attractive as their proximity often
makes it possible to resolve their radio jets down to sub-parsec scales, and to
investigate the behavior of the emission at different wavelengths.

\section{Radio Galaxies at gamma-ray energies}
At high gamma-ray energies (HE; $>100$ MeV) Fermi-LAT has reported the
detection of about ten misaligned RGs, mostly of the Fanaroff-Riley-type~I
(FR~I) \cite{Abdo10a,Ackermann11,Nolan12}, see also Fig. 1.
At very high energies (VHE), only four RGs have been clearly identified so far:
Cen~A ($d\simeq 3.8$ Mpc), M87 ($\simeq 16.7$ Mpc) and the Perseus Cluster
($d\sim77$ Mpc, $z\sim 0.018$) sources NGC~1275 and IC~310. A detection
of the RG 3C66B at VHE was reported by MAGIC (2007 observations \cite{Aliu09}),
but the VHE emission seems not sufficiently disentangled from the nearby
(separation $\theta \sim 0.12^{\circ}$) IBL blazar 3C66A to consider it here as a
confirmed VHE emitter.\\
\begin{figure}
 \includegraphics[height=0.5\textheight]{fig1.pdf}
  \caption{ Radio Galaxies (RGs) detected at gamma-ray energies.
  Out of $\geq$ 886 AGN ("clean" sample) detected by Fermi at high
  energies (HE) less than a few percent are RGs. Only four of these
  RGs are currently confirmed to be VHE emitters.}
\end{figure}

\noindent The nearest AGN {\it Cen A}, hosting a black hole of mass
$M_{\rm BH} \simeq (0.5-1) \times 10^8\,M_{\odot}$, was the second RG
to be detected at VHE energies in a deep ($>120$h) exposure by H.E.S.S.
\cite{Aharonian09}. Cen~A is a weak VHE source with an integral flux
$I(>250$ GeV) above 250 GeV of about $0.8\%$ of the Crab Nebula
(equivalent to an apparent isotropic VHE luminosity of $L(>250$ GeV)
$\simeq 2 \times 10^{39}$ erg/s). The VHE spectrum extends up to
$\sim5$ TeV and is compatible with a single power law of photon index
$2.7\pm0.5$. No significant VHE variability has been found so far, so
that an extended origin of the VHE emission \cite{Hardcastle} cannot
be simply excluded yet.  At high energies, Fermi-LAT has also detected
HE $\gamma$-rays up to 10 GeV from the core of Cen~A \cite{Abdo10b}.
The HE light curve (15d bins) of the core is consistent with no variability
and the HE $\gamma$-ray spectrum can be described by a photon index
comparable to the one in the VHE. A simple power law extrapolation of
the HE spectrum, however, cannot account for the observed TeV flux. This
could be indicative of an additional emission contribution to the VHE domain
beyond the common synchrotron-Compton emission \cite{Chiaberge01},
emerging at the highest energies \cite{Rieger12a}. At GeV energies, Fermi-LAT
has also detected $\gamma$-ray emission from the giant radio lobes of Cen A.
A recent analysis of a larger HE data set revealed significant evidence for
a spatial extension of the $\gamma$-ray lobes beyond the radio (WMAP)
image \cite{Zhang12}, suggesting the need for a more complex ($\geq$
two-zone) lobe emission modeling than done so far. Note that the angular
resolution of H.E.S.S. ($\sim0.1^{\circ}$) clearly allows to exclude the lobe
regions as source of the detected TeV emission.\\
The Virgo Cluster galaxy {\it M87} was the first radio galaxy detected at
TeV energies \cite{Aharonian03}. Classified as an FR I, M87 is known to
host one of most massive black holes with $M_{\rm BH} \simeq(2-6)\times
10^9\,M_{\odot}$ and to exhibit a relativistic jet misaligned by $\theta
\simeq(15-25)^{\circ}$. These latter values are consistent with modest
Doppler boosting $D=1/[\Gamma_j (1-\beta\cos\theta)] \simlt 3$ (for review,
see e.g. \cite{Rieger12b}). M87 has shown particularly interesting VHE
features within recent years, including rapid day-scale variability (observed
flux doubling time scales $\Delta t_{\rm obs} \sim 1$) above 350 GeV
during active source states (with flux levels sometimes exceeding $10\%$
of the Crab Nebula), and a hard spectrum (compatible with a power law of
photons index $\Gamma \simeq 2.2$) extending beyond $10$ TeV
\cite{Aharonian06,Albert10,Aliu12,Abramowski12,Beilicke12}.
Both the hard VHE spectrum and the observed rapid VHE variability are
remarkable features for a misaligned AGN. M87 is the only RG where
clear evidence for day-scale TeV variability (see Fig. 2) has been found,
although there appear to be some hints for day-scale variability in IC~310
(and for month-type in NGC~1275) as well (see below). The observed
variations in M87 seen at VHE are the fastest compared to any other
waveband so far, and imply a compact $\gamma$-ray emitting region
($R < c\Delta t_{\rm obs} D$), comparable to the Schwarzschild radius
$r_s=(0.6-1.8) \times10^{15}$ cm of its black hole. This could point to the
black hole vicinity as the most likely origin of the observed, variable VHE
radiation, a consideration supported by high-resolution radio (VLBA)
observations during the 2008 (February) VHE high state of M87, showing
that the VHE outburst was followed by an ejection of a new radio component
from very close to the black hole \cite{Acciari09}.
\begin{figure}
  \includegraphics[height=0.4\textheight]{fig2.pdf}
  \caption{The 2010 April VHE flare of the radio galaxy M87 as observed
  by different Cherenkov Telescopes. Significant day-scale activity is
  evident. The curve shows a fit of an exponential function to the data.
  From Ref.~\cite{Abramowski12}.}
\end{figure}

\noindent The Perseus Cluster galaxy {\it IC~310} has been detected by
MAGIC above 300 GeV in about 21h of data (taken in 2009/2010) at an
average level of $\sim 3\%$ of the Crab Nebula \cite{Aleksic10}. The
source has been originally classified in the literature as a head-tail RG,
but a recent radio analysis finds little evidence for jet bending \cite{Kadler12}
(see also Ref. \cite{Rector99} for the proposal that IC~310 may instead be a
weakly beamed blazar). The VHE spectrum measured between 150 GeV
and 7 TeV is very hard (even harder than in M87) and compatible with a
single power law of photon index $\Gamma \simeq 2.0$. There is clear
evidence for VHE variability on yearly and monthly time scales, with some
indications for day-scale activity found in a new analysis \cite{Eisenacher},
features that are reminiscent of the VHE activity seen in M87. Note however,
that due to its larger distance, IC~310 needs to be intrinsically much more
luminous.\\
In addition to IC~310, MAGIC has also recently reported the detection of
the central dominant (FR I) cluster galaxy {\it NGC~1275} above $\sim
100$ GeV during enhanced high energy (Fermi-LAT) activity in 46h of
data (taken between 08/2010-02/2011) \cite{Aleksic12}. This source is
known to have radio jets misaligned by $\simgt30^{\circ}$. Fermi-LAT
data indicate significant flaring activity above 0.8 GeV down to time
scales of days \cite{Brown11}. At VHE energies the situation is less
evident, with no evidence for VHE variability found in the noted data set,
while a recent, improved analysis of an earlier (10/2009-02/2010) VHE
data set provides hints for a possible month-type variability \cite{Colin}.
NGC1275 shows a very steep VHE spectrum (power law photon index
$\Gamma \simeq 4.1$) extending up to $\sim500$ GeV. When compared
with the hard HE (Fermi-LAT) spectrum (photon index $\Gamma\simeq
2.1$), this implies the existence of a break or cut-off in the SED around
some tens of GeV.

\begin{figure}[ht]
\includegraphics[height=.5\textheight]{fig3.pdf}
\caption{Possible theoretical scenarios for the origin of the variable
VHE emission as proposed for M87. For details, including some pros
and cons of the models, see Ref. \cite{Rieger12b}.}
\end{figure}

\section{On the origin of variable TeV gamma-rays from M87}
Being one of the best studied extragalactic object, M87 has become
a prominent benchmark for theoretical research. A multitude of VHE
emission scenarios, operating on different scales and relying on
either hadronic and/or leptonic radiation processes, have been
developed and applied within recent years (cf. Fig. 3). The observed
day-scale variability favors an origin of the variable TeV emission on
scales of the inner jet and below, although outer-jet scenario may
possibly contribute to the overall quiescent VHE source state (e.g.,
\cite{Rieger12b} for discussion).\\
(1) One recent theoretical development in the context of leptonic models
considers relativistic (Petschek-type) reconnection as acceleration
mechanism of electrons to high energy \cite{Giannios10}, see Fig. 4
for illustration. Compton up-scattering of external photons could then
possibly account for the production of the observed, variable VHE
gamma-rays \cite{Cui12}. Efficient reconnection in a highly magnetized
($\sigma\sim 100$) electron-proton jet could lead to an additional
relativistic velocity component of the ejected plasma ($\Gamma_r
\simeq \sqrt{\sigma}$, at some angle to the main outflow direction)
with respect to the mean bulk flow of the jet. This makes strong,
differential Doppler boosting effects for the emitting zone possible,
circumventing the problem of a modest Doppler factor for the general
bulk flow. At present, there seem to be a number of open issues that
need to be clarified within such an approach in order to assess its
explanatory potential. This includes the fact that (i) typically much
lower magnetizations ($\sigma\sim 10$) are expected for electron-proton
(disk-driven) jets in AGN (in contrast to BZ-driven electron-positron
jets), that (ii) at the anticipated scale of VHE production a non-negligible
magnetic guide field may persist allowing for weak dissipation only,
and that (iii) the achievability of power-law electron acceleration much
beyond the thermal Lorentz factor $\sqrt{\sigma} m_p/2m_e$ is not
clear.\\

\begin{figure}
\includegraphics[angle=270, width=7.5cm]{fig4a}
\includegraphics[angle=270, width=7.5cm]{fig4b}
\caption{LEFT: Illustration of a jets-in-jet scenario \cite{Giannios10},
where a variety of 'mini-jet' features are thought to be induced by
relativistic magnetic reconnection events within the general jet flow.
This could lead to an additional velocity component ($\Gamma_r$)
relative to the main flow and allow a favorable orientation with
respect to the observer. RIGHT: Illustration of a hadronic model,
where interactions of the jet with a massive gas cloud leads to
Fermi (shock)-type proton acceleration and introduces a sufficient
target density to allow for efficient $pp$-collisions \cite{Barkov12}.}
\end{figure}
\noindent(2) In the context of hadronic models, on the other hand, one new
development considers inelastic proton-proton (pp) collisions (and $\pi^0$-decay)
in the jet for the origin of the observed VHE emission. While AGN jets have
been often considered not to carry (on average) a sufficient amount ($n_p$)
of target matter \cite{Celotti98,Wardle98} to ensure the efficiency of this (pp)
process (given its characteristic timescale $t_{pp}' \simeq 10^{15}/n_p'$ sec),
it seems possible that interactions of a red giant star or a massive, dense gas
cloud (size $r_c$) with the base of the jet (radius $r_j$) could occasionally
introduce the required amount of matter and thereby potentially drive fast VHE
activity \cite{Barkov10,Barkov12}. Recent calculations suggest that such a
scenario could account for the observed temporal and spectral properties if
the jet is powerful enough and a sufficiently large fraction ($\propto r_c^2/
r_j^2$) of it can be channeled into VHE $\gamma$-ray production. Given the
observed large opening angle and transversal dimension of the radio jet in
M87, models of this type require a spine-shear-type jet configuration with
most of the energy flux concentrated into a narrow core.\\
(3) An alternative explanation has been pursued within magnetospheric models
(see e.g. \cite{Rieger11} for review), motivated by the observed day-scale VHE
variability and the observed radio-VHE connection in 2008. For a black hole
mass in M87 expected to be in the range $M_{\rm BH} \simeq (2-6) \times 10^9\,
M_{\odot}$ \cite{Rieger12b}, anticipated variability  timescales could possibly be
as small as $\sim r_g/c =0.1$ d, allowing to accommodate the observed rapid
VHE variability.
One of the latest developments in this context analyses the possible link between
VHE $\gamma$-ray activity and jet formation \cite{Levinson11}. A schematic
representation is shown in Fig. 5:\\
\begin{figure}
\includegraphics[angle=270, width=.95\textwidth]{fig5}
\caption{VHE $\gamma$-ray production within a black hole-magnetospheric
scenario \cite{Levinson11}: Electrons and positrons are accelerated in a
vacuum gap to high Lorentz factors. The gap is exposed to soft radiation field
of the disk. Curvature emission and inverse Compton up-scattering of ambient
disk photons yields VHE $\gamma$-rays and triggers $e^+e^-$-cascades.
The resultant large pair multiplicity establishes the force-free outflow seen
in VLBA images.}
\end{figure}
In this scenario, a maximally rotating black hole, which is embedded in an external
magnetic field $B$ generates an electric potential difference across a (vacuum) gap
of height $h$, given by $\Delta V_e=1.7\times10^{20} B_3 M_9(h/r_s)^2$ Volts,
where $B_3=B/10^3$ G and $M_9=M_{\rm BH}/10^9~M_{\odot}$. Pair creation in
a hot accretion flow (i.e. by $\gamma\gamma$-annihilation of MeV photons in an
ADAF or RIAF) then leads to the injection of primary electrons and positrons into the
gap, with the charged particles being quickly accelerated to high Lorentz factors.
The number density of ADAF bremsstrahlung MeV photons is $n_{\gamma} \propto
\dot{m}^2$, so that the $e^+e^-$-injection rate inside the magnetosphere roughly
becomes
$n_{\pm} \sim \sigma_{\gamma\gamma} n_{\gamma}^2 r_s \propto \dot{m}^4$.
When expressed in terms of the Goldreich-Julian (GJ) density $n_{\rm GJ} =
\Omega B/(2\pi e c) \propto \dot{m}^{1/2}$, the ratio becomes \cite{Levinson11}
\[
n_\pm/n_{\rm GJ}\simeq 0.06~\dot{m}_{-4}^{7/2} M_9^{1/2}.
\]  This numerical estimate has some uncertainty as the bremsstrahlung spectrum
is dependent on $T_e$, which itself depends on the accretion rate $\dot{m} \sim
10^{-4}$ (expressed in terms of the Eddington rate). The sensitive dependence on
$\dot{m}$ suggests, however, that a gap (with $n_\pm/n_{\rm GJ} < 1$) can be
formed during periods of low accretion. Gap-type particle acceleration and
subsequent curvature emission and, in particular, inverse Compton (IC) up-scattering
of disk photons then produces $\gamma$-rays with a spectrum extending up to $10^4$
TeV. Different contributions dominating at different energy bands could then lead
to a different variability behavior as e.g. seen by MAGIC during the 02/2008-VHE
high state in M87 \cite{Albert10}. For low accretion rate, VHE photons with
energies below several TeV can escape $\gamma\gamma$-absorption (see e.g.,
\cite{Levinson11,Rieger12b} for details). On the other hand, IC photons with energies
above $\sim10$ TeV interact with the ambient radiation field initiating pair cascades
just above the gap. This leads to a large pair multiplicity ($M=n_{2\pm}/n_{\rm GJ}$
up to $10^3$). A force-free outflow becomes established above the pair formation
front and accounts for the jet feature seen in radio VLBA. Given the sensitive
dependence of gap formation (primary pair injection) and photon escape (optical
depth for $\gamma\gamma$-absorption) on accretion rate, any emission from the
gap is likely to be intermittent. Modest changes in accretion rate, for example, may
naturally give rise to the variability of the TeV emission and to fluctuations of the
resultant force-free outflow.

\section{Conclusions}
Radio galaxies have emerged as a new and particularly interesting VHE source
class on the extragalactic sky. Evidence for $\gamma$-ray variability, as seen in
most of the sources (apart from Cen~A) suggests that the observed $\gamma-$ray
emission originates in a compact region. In reality, the situation may be somewhat
more complex with extended regions contributing to the average/low states of the
source, and different (additional) zones dominating during high VHE source states.
In the case of M87, the observed VHE flaring characteristics (such as, e.g., the hard
$\gamma$-ray spectrum extending beyond 10 TeV and the day-scale VHE activity)
are of fundamental interest as they allow us to probe the near black hole environment
and to study the possible link between VHE activity and jet formation in AGN. Radio
galaxies like M87 thus appear to be particularly promising targets for future \& more
sensitive instruments like CTA.

 \begin{theacknowledgments}
{Discussions with Felix Aharonian, Amir Levinson, Maxim Barkov, Matthias Beilicke
and Martin Raue, and financial support by a LEA Fellowship are gratefully
acknowledged.}
\end{theacknowledgments}

\bibliographystyle{aipproc}

\input{references.bbl}

\end{document}

High Energy Astrophysical Phenomena (astro-ph.HE)


The Starburst Intensity Limit And Its Ultraviolet Implications

Abstract: Our recent work on starbursts, particularly in the ultraviolet (UV), is summarized. The intrinsic UV fluxes of UV selected starbursts can be derived from UV data alone because, to first order, their dust behaves like a foreground screen. This allows a comparison of the bolometric effective surface brightness S_e of UV selected starbursts to other starburst samples. Starbursts have a robust (90th percentile) upper limit S_e <~ 2.0e{11} L_{sun} kpc^{-2}, which strongly suggest that their global star formation intensities are regulated. The mechanism(s) involved in the regulation are not yet clear. The dust attenuation corrections for high-z starbursts are significant. Calculations of the rate of evolution in the early universe based on zero dust interpretations are probably underestimated by about an order of magnitude. Hence the early universe was not quiescent, but obscured.

\section*{Introduction}

The pace of cosmic evolution is marked out by the condensation of
stars from the ISM.  Starbursts, regions of intense massive star
formation which can dominate the bolometric output of galaxies, are
crucial to our understanding of this evolution at all redshifts.  In
the local universe about 25\%\ of high mass star formation is
contained within starburst galaxies (Heckman, 1997; Gallego {\it et~al.\/}\
1995).  At high redshift ($z > 2$), starbursts are the the easiest
galaxies to detect, and can be used to directly trace the rate of
evolution (Madau {\it et~al.\/}\ 1996). So, to understand the cosmological pace
of evolution, we must understand starbursts.  We would like to know
how to derive intrinsic star formation rates and the physics that
determines these rates.  Can the starburst ``knob'' be cranked up to
arbitrary levels or is its amplitude limited?

Recently we imaged a small sample of starbursts in the
ultraviolet (UV) using the HST (Meurer {\it et~al.\/}\ 1995; hereafter M95).
The basic UV structure of a starburst consists of prominent compact
star clusters embedded in a more diffuse cloud.  In the nearest
starbursts the diffuse light starts to resolve into stars, hence star
formation within starbursts occurs in both clustered and diffuse
modes. In M95 we compared the effective surface brightness, $S_e$ (the
surface brightness within the radius containing half the total
emission), of starbursts and concluded that there is an upper limit to
$S_e$.  This suggests that the intensity of starbursts is limited.
Lehnert and Heckman (1996) found a similar surface brightness limit in
their combined FIR/\Hline{\mbox{$\alpha$}}\ study.  This prompted our more
comprehensive investigation of starburst surface brightnesses (Meurer
{\it et~al.\/}\ 1997; hereafter M97) which we report on here.

\section*{Attenuation and Reddening Due to Dust}

Vacuum-UV imaging is indispensable for understanding starbursts since
it allows the direct detection of the hot-high mass stars that power
starbursts.  In contrast other tracers of high-mass stars,
such as \Hline{\mbox{$\alpha$}}\ (and other emission lines), radio continuum, and
far-infrared (FIR) emission originate in the ISM and
only indirectly trace star formation.  Dust presents the biggest
obstacle to understanding UV observations, since dust extinction and
scattering are strongest in the UV resulting in a net attenuation of
the UV flux.  Fortunately, our work has proven that UV observations
can reveal essential properties of high mass star formation, even in
the presence of modest amounts of dust.

\begin{figure}[t]
\centerline{\hbox{\epsfig{figure=meurer_uvuni_fig1.ps,height=7cm}}
\parbox[b]{7cm}{\footnotesize {\bf Figure 1.} The ratio of far
infrared to UV fluxes compared to UV spectral slope $\beta$ (flux
density $f_\lambda \propto \lambda^\beta$) which is equivalent to a UV
color.  The data correspond to UV-selected starbursts.  The fluxes are
measured in in the FIR by IRAS and in the UV ($\lambda \approx
2300$\AA) by HST or IUE.  The hatched region shows the expected
$\beta$ for naked ionizing populations.  The dashed line shows the
expected relationship for a starburst having an intrinsic $\beta_0 =
-2.5$ behind a foreground screen of dust obeying the Calzetti {\it et~al.\/}\
(1994) attenuation law. \vspace*{9mm}}}
\end{figure}

We are saved because dust is not a sink for photons; the absorbed
radiation is reemitted thermally in the FIR.  Hence the ratio of FIR
to UV fluxes is an indication of net attenuation.  Figure~1 shows a
positive correlation between UV spectral slope $\beta$ and
$\log(F_{\rm FIR}/F_{\rm UV})$ for a sample of strongly star forming
galaxies as indicated by their strong recombination emission line
spectra.  {\em A priori\/} we expect their UV spectra to be that of an
ionizing population, and hence to have an intrinsic color in the
hatched region of the plot.  This correlation is readily understood in
terms of a simple foreground screen dust geometry, as shown by the
model line.  Although the dust distribution, in exact detail, may not
be a {\em homogeneous\/} foreground screen, and other quantities may
be involved in the correlation, this is largely unimportant for
determining the attenuation correction.  So long as the FIR emission
results from dust reradiating the light absorbed in the UV and
optical, any model that reproduces the correlation in Fig.~1 will
recover the UV flux of ionizing populations to within a factor of
$\sim$3 (corresponding to the spread in $\log(F_{\rm FIR}/F_{\rm
UV})$).

\section*{The Starburst Intensity Limit}

\begin{figure}[t]
\centerline{\hbox{\epsfig{figure=meurer_uvuni_fig2.ps,width=7cm}
\epsfig{figure=meurer_uvuni_fig3.ps,width=7cm}}}
\parbox{140mm}{\footnotesize {\bf Figure 2} (left). Bolometric
luminosity $L_{\rm bol}$ and effective surface brightness $S_e$
plotted against effective radius $R_e$ for UV-selected starbursts and
star clusters.  The dotted and dashed lines are the 90th and 50th
$S_e$ percentiles for the starbursts.  The circles corresponds to
local starbursts ($D < 75$ Mpc), triangles to starbursts at $z = 0.4$,
stars to high redshift ($2.2 < z < 3.5$) starbursts, and crosses
to clusters within starbursts having $D < 10$ Mpc.  \\
{\bf Figure 3} (right).  $S_e$ distributions for star clusters (from
Fig.~2), UV selected starbursts (from Fig.~2), FIR selected starbursts
measured in the FIR and \Hline{\mbox{$\alpha$}}, FIR selected starbursts observed at
21cm, and normal disk galaxies observed at \Hline{\mbox{$\alpha$}}\ (see M97 for
sample details). The dotted lines show the 90th percentiles of the
distributions. In the bottom panel we show preliminary $S_e$
measurements of an UV ($\lambda \approx 1500$\AA) image of M81 at
obtained with the UIT (Hill {\it et~al.\/}\ 1992). }
\end{figure}

The evidence for a starburst intensity limit is summarized in Figs 2
and 3.  In Fig.~2 we show measurements of our rest frame UV sample
after dust attenuation correction as outlined above.  For comparison
with the other samples, shown in Fig.~3, luminosities and surface
brightnesses are converted to bolometric units following the models of
Leitherer \&\ Heckman (1995).  Note that the shape of the low $S_e$
end of the distributions in Fig.~3 are determined by (largely
arbitrary) selection effects.  However there is no selection against
high surface brightness.  The three starburst $S_e$ distributions have
90th percentile upper limits within a factor of three of $S_{\rm lim}
\approx 2 \times 10^{11}\, L_\odot\, {\rm kpc^{-2}}$, which we call
the starburst intensity limit.  This limit is very robust, applying to
starbursts with $R_e \sim 0.1$ to 10 kpc, to both local starbursts,
and systems out to redshift $z \sim 3$, to both UV/optically selected
starbursts and FIR selected starbursts, and to observations obtained
in the UV, \Hline{\mbox{$\alpha$}}, FIR, and radio continuum.  The lack of a strong
$\lambda$ effect, indicates that it corresponds to a limit on star
formation intensity rather than being an opacity effect, and is
further vindication of the UV attenuation estimates. This limit
strongly indicates that the global intensity of star formation is
regulated within starbursts.

In M97 we explored two physical mechanisms that may be responsible for
the intensity limit: galactic winds (Heckman {\it et~al.\/}\ 1990), and
gravitational stability of the inner disk (following Kennicutt, 1989).
Although we find evidence that the intensity limit is related to both
mechanisms, neither predicts the value of the intensity limit, nor its
robustness. Hence {\em determining the physics of the starburst
intensity remains a major theoretical challenge\/}.

Although robust, this limit does not hold for all
star forming scales.  In particular it does not hold
for the star clusters embedded within starbursts.  These are small
$R_e \lesssim 10$ pc (M95) and can have $S_e$ two orders of magnitude
or more intense than $S_{\rm lim}$.  This is not a contradiction
because $S_{\rm lim}$ refers to a global quantity dominated by diffuse
light, whereas individual clusters usually do not dominate the total
luminosity of starbursts.

Star formation in normal disk galaxies is typically three orders of
magnitude less intense than in starbursts, as seen in Fig.~3.
However, this may largely reflect the different star formation
patterns.  Typically star formation fills the central regions of
starburst galaxies, whereas in normal galaxies it usually consists of
\mbox {H\thinspace{\footnotesize II}}\ regions, often arranged in spiral arms or rings, but otherwise
fairly isolated.  We illustrate how large a difference covering factor
can make in the bottom panel of Fig.~3 using a UIT UV image of M81
(Hill {\it et~al.\/}\ 1992).  Its disk averaged $S_e$, marked ``D'', is the
surface brightness within the elliptical aperture containing half the
UV light (the same technique used to measure the $S_e({\rm H}\alpha)$
values shown in the bottom panel).  The {\em isophotal\/} $S_e$,
marked ``I'', is the surface brightness within the isophote containing
half the total flux (in other words, pixels having this or higher
surface brightness comprise half the total UV flux). In effect, this
value does not include any of the ``empty space'' between \mbox {H\thinspace{\footnotesize II}}\
regions. The $S_e$ of one of M81's brightest \mbox {H\thinspace{\footnotesize II}}\ regions is marked
``H''.  This comparison suggests that star formation in normal
galaxies, averaged over only the actively star forming regions, may
occur with the same intensity distribution as starbursts.  More
measurements of normal galaxies are needed to test this hypothesis.

\section*{Implications}

Our results have far reaching implications. We demonstrated that
intrinsic UV fluxes of starbursts can be recovered from UV data alone
(Fig.~1).  This confirms the importance of UV observations to the
study of extragalactic high mass star formation.  The dust attenuation
corrections are significant in the UV.  The median $\beta = -1.1$ of
our UV sample corresponds to an attenuation factor of $\sim$ 8 at
$\lambda \approx 2300$\AA\ and $\sim 15$ at $\lambda = 1600$ \AA.  The
latter value corresponds to the rest $\lambda$ of the Lyman drop out
galaxies observed at high redshift (Madau {\it et~al.\/}\ 1996).  These have
properties (including $\beta$) consistent with those of local UV
selected starbursts.  Madau {\it et~al.\/}\ (1996) use the UV luminosity
density of these systems to infer that the universe was fairly
quiescent (low star formation rate density, or equivalently metal
production rate) at redshift $z \approx 3$.  However they applied no
dust extinction corrections.  After applying our UV attenuation
corrections, and adopting the same cosmology as Madau {\it et~al.\/}\ ($H_0 =
50\, {\rm km\, s^{-1}\, Mpc^{-1}}$, $q_0 = 0.5$) we find a lower limit
of $\sim 3 \times 10^{-3} \, {\rm M_\odot\, yr^{-1}\, Mpc^{-1}}$ for
the metal production rate at $z \approx 3$.  This is about six
times higher than their estimate of the Hubble time averaged metal
production rate.  We find that rather than being quiescent we are
observing an active but moderately obscured early universe.

{\footnotesize {\bf Acknowledgements:} I thank my collaborators for
their work in the M95 and M97 papers: Don Garnett, Anne Kinney, Matt
Lehnert, Claus Leitherer, James Lowenthal, Carmelle Robert, and
especially Tim Heckman.  Those papers and this benefited from many
stimulating discussions with Daniela Calzetti.}