Verifying the Planetary Nebula Luminosity Function Method

Abstract: The planetary nebula luminosity function (PNLF) technique for determining distances to galaxies now has been applied to 34 galaxies, including 6 in the Virgo cluster and 3 in the Fornax cluster. Of these, 16 galaxies are late-type or spirals and presumably contain Cepheid variables useful for verifying the PNLF method. For 7 of these galaxies, Cepheid distances exist; the PNLF distances agree with the Cepheid distances within the dispersion of 8% and within a zero-point offset of 1%. In addition, 3 small groups were studied (NGC 1023, Leo I, and Coma I) where both spiral and elliptical distances were obtained to investigate the magnitude of any systematic dependence on spiral versus elliptical Hubble type. None was found. Since the PNLF method agrees well with the Cepheid system, and there is no measurable dependence on Hubble type, it follows that PNLF distances to the ellipticals in Virgo and Fornax also are on the Cepheid scale. This conclusion is strengthened by the Cepheid distances to several Virgo galaxies and the recent determination of a Cepheid distance to Fornax.

\section{Introduction}

The Planetary Nebula Luminosity Function (PNLF) technique for
determining distances to galaxies is described in \cite{Jacoby92}.
Briefly, the method relies on the constancy of the [OIII] $\lambda5007$
luminosity at which the population of high luminosity planetary nebulae
(PN) declines rapidly.  While the number of PN decreases monotonically
with increasing luminosity, there exists a luminosity beyond which PN
are not observed at all.  This limiting luminosity is a consequence of
several factors (Jacoby (1989), Kaler \& Jacoby (1991), M\'endez \mbox{\it et al.\ }
(1993), Han \mbox{\it et al.\ } (1994), Stanghellini (1995), Richer \mbox{\it et al.\ } (1997)), the
most important being the independence between progenitor
age (or initial mass) and the final white dwarf (or central star) mass
over reasonable age ranges ({e.g.,\ } 3--11 Gyrs).  For the great majority of
progenitors and therefore, for normal populations, the critical PN
luminosity varies by less than 0.1 mag.

The PNLF technique is one of the few that has been demonstrated to
yield consistent results in both elliptical and spiral galaxies at
distances up to 20 Mpc.  Consequently, it offers a path to unify the
population I distances ({e.g.,\ } Cepheids) with population II objects ({e.g.,\ }
elliptical galaxies). The former are needed for calibration while the
latter are needed to define distances to rich clusters.  Other
promising cross-population indicators include SN Ia when corrected for
decline rate variations (\cite{Hamuy96}), surface brightness
fluctuations (SBF), and globular cluster luminosity functions (GCLF).

This paper describes recent results in which the PNLF has been extended
to late-type galaxies in order to improve the zero-point calibration,
and new observations in M87 to assess the legitimacy of the arguments
posed against the use of the PNLF by \cite{Bottinelli91} and
\cite{Tammann92}.

\section{What Makes a Good Distance Indicator?}

Several speakers at this meeting presented their views on the important
qualities of a good distance indicator.  They noted the advantages of
an indicator that is luminous, easily identified, and has easily
measured properties ({e.g.,\ } magnitudes). My list focuses more
on the physical nature of the indicator rather than its technical
properties. Thus, a good distance indicator:

\medskip
\begin{enumerate}
\renewcommand{\theenumi}{\arabic{enumi}}
\item Has a good zero point calibration
\item Has a good prescription for metallicity correction
\item Has a good prescription to correct for effects of stellar ages
\item Can be corrected for effects of foreground and internal extinction
\item Has a good physical rationale
\item Can be tested (a) against other methods and (b) by other investigators
\end{enumerate}
\medskip

\noindent
This paper concentrates on the first of these. It is worth reviewing, though,
how well the PNLF satisfies the remaining 5 criteria.

Metallicity effects were modeled by \cite{Dopita92}. Limited
observational testing was performed by \cite{CJ92}. Most galaxies of
interest have metallicities within a factor of 2 of solar abundances,
and the predicted and observed errors in PNLF distances are smaller
than 5\% over this range.

The effects of population age on distance have been modeled by \cite{Dopita92},
\cite{Mendez93}, and \cite{Stanghellini95} and shown to be
$<$5\% for galaxies having ages in the oberved range between 3 and 11
Gyrs.  Direct tests are complicated by our present inability to measure
population ages accurately, but PNLF distances to young (LMC, SMC,
M101) and old populations (M31's bulge, M81's bulge) where distances
are known from Cepheids, fail to detect any age effect at all.

Internal extinction is less of a problem than intuitively suspected.
In ellipticals, extinction is not an issue since the dust density is
very low. In spirals, significant errors are expected
if internal extinction in the galaxy is ignored. Observations, though,
fail to reveal any measurable distance errors in the 7 late-type
calibrators (Feldmeier, Ciardullo, \& Jacoby 1997).
Feldmeier, Ciardullo, \& Jacoby (1997) modeled the
effects of dust to understand this unexpected situation, and
found that PN, which generally have scale heights well above the
population I disk of a spiral, are either so extincted by heavy dust
that they fall out of the PNLF sample, or they are so little affected
that their magnitudes are not significantly diminished.

The excellent agreement between the PNLF and Cepheids (as well as SBF
and with other methods to a lesser degree) demands that a physical
basis for the PNLF must exist. Before all the comparisons were made,
though, the theory had been described by \cite{Jacoby89}, \cite{Dopita92},
and \cite{Mendez93}; recently, Stanghellini (1995) investigated
the H$\beta$ PNLF. In short, it has proven easy to reproduce the constancy
of the PNLF provided the population age is within the range of 3 to 11 Gyrs.
If the progenitors are as young as 0.5 Gyrs, the PNLF brightens by $\sim$0.3 mag.
And, if all progenitors in a galaxy are much older than 11 Gyrs, they fail to
produce observable PN at all.

The requirement that a distance indicator be testable against another
method is fundamental to the concept of the scientific method.
Because we never know the ``right answer'' in the distance scale business,
we turn to intercomparisons between different methods assuming that
2 independent methods are very unlikely to repeatedly yield the
same wrong answer. If a method is not testable, it relies solely
on the validity of a model and scientists generally agree that models must
be tested. By inference, an untestable indicator is equivalent
to an untestable model. Fortunately, the PNLF can be tested against
numerous methods (see \cite{CJT93}, \cite{Jacoby95}, \cite{Feldmeier97}).

The second component of the last requirement is that multiple
investigators must be able to derive the same answer using the same
technique. This sounds simple enough, and again, is fundamental to the
scientific method.  Results from some methods, however, cannot be
reproduced at a later time should a question arise about their
validity.  The most obvious of these methods is supernovae. A second
observer cannot go back in time to observe a supernova in order to
check on the observational accuracy of a prior observer's measurements,
or to utilize a superior instrument. Thus, SN Ia fail to satisfy this
requirement.

The PNLF technique satisfies the prescription for a good distance indicator
on each count. Cross-testing with other methods (the last and most
important criterion) shows that disagreements between the PNLF and other
reliable methods ({e.g.,\ } Cepheids) are smaller than 8\%. Thus, systematic
errors due to extinction, age, metallicity, or application of the method
are not accumulating beyond the 8\% level. In fact, when consideration
is made for the error contribution from the Cepheid distances, the PNLF
errors must be smaller than $\sim$5\%.

\begin{table}
\begin{center}
\begin{tabular}{llrrl}
{\bf Galaxy} & {\bf Type} & {\bf Nr.~PN} & {\bf (m-M)$_0$} & {\bf Reference} \\
     &    &    &               & \\
\multicolumn{5}{l}{\bf Local Group} \\
LMC &SBm &42 &$18.44 \pm 0.18$ & \cite{JWC} \\
SMC &Im &8 &$19.09 \pm 0.29$  & \cite{JWC} \\
185 &dE3p &4 &... &  \cite{Ciardullo89} \\
205 &S0/E5p &12 &$24.68 \pm 0.35$ & \cite{Ciardullo89} \\
221 &E2 &9 &$24.58 \pm 0.60$ & \cite{Ciardullo89} \\
224 &Sb &104 &$24.26 \pm 0.04$  & \cite{Ciardullo89} \\
     & &    &    &              \\
\multicolumn{5}{l}{\bf NGC 1023 Group} \\
891 &Sb &34 &$29.97 \pm 0.16$ & \cite{CJH91} \\
1023 &SB0 &97 &$29.97 \pm 0.14$ & \cite{CJH91} \\
     & &    &    &               \\
\multicolumn{5}{l}{\bf Fornax Cluster} \\
1316 &S0p &58 &$31.13 \pm 0.07$ & \cite{MCJ93} \\
1399 &E1 &37 &$31.17 \pm 0.08$  & \cite{MCJ93} \\
1404 &E2 &19 &$31.15 \pm 0.10$  & \cite{MCJ93} \\
     &    &    &               & \\
\multicolumn{5}{l}{\bf Leo I Group} \\
3377 &E6 &22 &$30.07 \pm 0.17$  & \cite{CJF89} \\
3379 &E0 &45 &$29.96 \pm 0.16$  & \cite{CJF89} \\
3384 &SB0 &43 &$30.03 \pm 0.16$ & \cite{CJF89} \\
3368 &Sab &25 &$29.91 \pm 0.15$ & \cite{Feldmeier97} \\
     &    &    &               & \\
\multicolumn{5}{l}{\bf Virgo Cluster} \\
4374 &E1 &37 &$30.98 \pm 0.18$  & \cite{JCF90} \\
4382 &S0 &59 &$30.79 \pm 0.17$  & \cite{JCF90} \\
4406 &S0/E3 &59 &$30.98 \pm 0.17$ & \cite{JCF90} \\
4472 &E1/S0 &26 &$30.71 \pm 0.19$ & \cite{JCF90} \\
4486 &E0 &201 &$30.73 \pm 0.19$  & this paper \\
4649 &S0 &16 &$30.76 \pm 0.19$  & \cite{JCF90} \\
     &    &    &               & \\
\multicolumn{5}{l}{\bf Coma I Group} \\
4278 &E1 &23 &$30.04 \pm 0.18$  & \cite{JCH96} \\
4494 &E1 &101 &$30.54 \pm 0.14$ & \cite{JCH96} \\
4565 &Sb &17 &$30.12 \pm 0.17$  & \cite{JCH96} \\
     &    &    &               & \\
\multicolumn{5}{l}{\bf NGC 5128 Group} \\
5102 &S0p &19 &$27.47 \pm 0.22$  & \cite{MCJ94} \\
5128 &S0p &224 &$27.73 \pm 0.04$ & \cite{Hui93} \\
5253 &Amorph &16 &$27.80 \pm 0.29$ & \cite{Phillips92} \\
     &    &    &               & \\
\multicolumn{5}{l}{\bf Other} \\
Bulge &Sbc &22	&$14.54 \pm 0.20$  & \cite{Pottasch90} \\
300  &Sc &10 &$26.78 \pm 0.40$  & \cite{Soffner96} \\
3031 &Sb &88 &$27.72 \pm 0.25$  & \cite{JCFB} \\
3109 &Sm &7 &$26.03 \pm 0.30$  & \cite{RM92} \\
3115 &S0 &52 &$30.11 \pm 0.20$  & \cite{CJT93} \\
4594 &Sa &204 &$29.76 \pm 0.13$ & \cite{Ford96} \\
5194/5 &Sbc/SB0 &38 &$29.56 \pm 0.15$  & \cite{Feldmeier97} \\
5457 &Sc &27 &$29.36 \pm 0.15$  & \cite{Feldmeier97} \\
\end{tabular}
\end{center}
\medskip
\caption{All distances must be increased by 0.06 mag to be placed on
the recent M31 Cepheid scale of Freedman \& Madore (1990).
A metallicity correction (Ciardullo \& Jacoby 1992) has been applied
to NGC 5253 only, because it has SMC-like abundances.}
\end{table}

\section{Some New Insights}

If we accept that the PNLF method yields accurate distances, it seems
odd, at first, that the errors from population differences and extinction
aren't larger. A potentially dominant population effect is age, as
discussed already. The key point is that intermediate age populations
all produce nearly identical central star masses. This follows from
the initial-to-final mass relation (\cite{Weidemann87}). That is,
for progenitor masses between 1 and 2 M$_\odot$ corresponding to ages
of about 1 to 10 Gyr, the central star mass will be in the narrow
range of $\sim 0.58\pm0.02$ M$_\odot$.  This narrow range is close
to that observed for white dwarfs (\cite{McMahon89}).

Another important effect arises in young ($<$0.5 Gyr) populations to
inhibit [OIII] luminous PN from forming. \cite{KJ91} showed, and
\cite{Dopita96} confirmed, that for young progenitors producing central
stars more massive than 0.65 M$_\odot$, the surface abundances are
strongly altered such that nitrogen is greatly enhanced.  The added
nitrogen competes with oxygen in cooling the nebula, to the detriment
of the [OIII] luminosity. Thus, PN deriving from young, massive
progenitors fail to populate the high luminosity end of the PNLF and
the effect of a young population on the PNLF is lost.

Similarly, metallicity seems like it ought to play a large role. A
competition exists, though, between the efficiency of the nebula to
radiate in [OIII] and the luminosity input from the central star.
Higher metallicity values enhance the nebula's ability to
radiate at $\lambda5007$.  The central star, however, is predicted to
have a lower mass and luminosity as a consequence of experiencing
higher mass loss prior to leaving the AGB. The reduced heating
compensates to first order for the enhanced radiative efficiency
as metallicity increases (\cite{Dopita92}).

\section{Summary of PNLF Distances}

Table 1 summarizes the available PNLF distances.
Distances are referenced to a zero-point based on M31 having
a distance of 710 kpc and a reddening of E(B--V)$ = 0.11$. This is
the baseline zero-point used in all our papers. To place
these distances on the recent scale where the distance to M31
is 770 kpc (\cite{Freedman90}) and E(B--V)$ = 0.08$, all distances
should be increased by 0.06 mag (3\%).

Ten galaxies of the 34 in this list are obvious spirals.
Six more have a significant late-type component for which Cepheid
distances either exist already or could be determined.

\section{Spirals}

Extending the PNLF technique to spiral galaxies requires extra care due
to 3 factors. First, potential confusion exists between PN and HII
regions. Second, spiral arms contain obvious dust lanes that could
reduce observed luminosities.  Third, a young population of stars
must exist in spiral arms.  Nevertheless, the advantage of working in
spirals is tremendous because their Cepheids provide the most accepted
reference distances to test the PNLF method for systematic errors.
Also, spirals are where zero-point calibrations are most believed.

To identify PN in spiral galaxies, we complement the $\lambda5007$
on-band/off-band imaging technique with on-band/R-band images at
H$\alpha$.  We define the following criteria for PN candidates in spirals.
A PN candidate must:

\medskip
\begin{enumerate}
\renewcommand{\theenumi}{\arabic{enumi}}
\item have a stellar PSF
\item be on the [OIII] image and absent on a continuum image
\item be absent in an R-band image
\item be absent or extremely weak in an H$\alpha$ image
\item not be in a spiral arm
\end{enumerate}
\medskip

\noindent
These criteria have been applied recently to observations of M51 (NGC 5194),
M96 (NGC 3368), and M101 (NGC 5457) by \cite{Feldmeier97}.  M96 and M101
distances exist already from HST Cepheid surveys.

\begin{figure}
\epsfysize 432pt \epsfbox{fig1.eps}
\caption{A comparison between PNLF and Cepheid distances. Solid circles
represent direct galaxy comparisons; the solid triangle is the calibrator
galaxy, M31, which, by definition, falls exactly on the dashed 1:1 line.
The points labeled N1023, Virgo, and Fornax represent comparisons
between elliptical (PNLF) and spiral (Cepheid) galaxies within the
same cluster. The lower panel  clarifies the level of disagreement
by plotting the relative differences in distances. Only M96 deviates
by more than 1$\sigma$; at this time, it is unclear if the disagreement
is significant, and if so, whether the PNLF or Cepheid distance is
discrepant.}
\end{figure}

Figure 1 illustrates the excellent agreement between the PNLF (adjusted to
the modern M31 distance and extinction) and Cepheid distance scales.
The PNLF distances are based solely on M31 as the zero-point calibrator,
yet the mean offset for the additional 6 Cepheid galaxies differs by
only 1\%. The random scatter of 8\% about the ridge line is consistent
with, or smaller than, the combined errors of the Cepheid and PNLF distances.
Note that some uncertainties are common to both methods (e.g., foreground
extinction) which slightly reduce the apparent uncertainty in the
combined distance error.

Figure 1 also includes comparisons between PNLF distances to elliptical
galaxies and Cepheid distances for different galaxies in the same
group.  These indirect comparison points provide further support, but
do not confirm the PNLF distances to the Virgo, Fornax, Leo, and NGC
1023 groups.

The lack of any evidence for systematic errors, either among the direct
or indirect comparisons, strongly implies that any residual systematic
errors for zero-points, population effects (age and metallicity),
internal extinction, and methodology (adopted PNLF shape, sizes of PNLF
samples) must be smaller than $\sim$8\%. Alternatively, one can insist
that a conspiracy exists among these parameters such that both Cepheid
and PNLF distances have errors that correlate.

\section{The Spiral--Elliptical Connection}

To a limited degree, the 3 indirect distance comparisons shown in
Figure 1 illustrate that the PNLF distances to ellipticals are on
the same system as the spiral distances. In addition, tests relying
purely on PNLF distances have been described for the NGC 1023 group
(\cite{CJH91}), the Coma I group (Jacoby, Ciardullo, and Harris 1996),
and the Leo group (\cite{Feldmeier97}). Figure 2 summarizes the 4 PNLF
and 2 Cepheid distances to the Leo group.

\begin{figure}
\epsfysize 216pt \epsfbox{fig2.eps}
\caption{A comparison between PNLF and Cepheid distances in the Leo Group.
NGC 3368 (M96) has been measured using both the PNLF (Feldmeier \mbox{\it et al.\ } (1997))
and with Cepheids (Tanvir \mbox{\it et al.\ } (1995)). In addition, group member NGC 3351
has been measured with Cepheids (Graham \mbox{\it et al.\ } (1996)).}
\end{figure}

Tests to different galaxies cannot be conclusive, since the two
galaxies are not guaranteed to be at exactly the same distance.
Nevertheless, the deviations in distance (spiral minus elliptical) for
the 3 purely PNLF tests are as follows. NGC 1023 group: $0.00\pm0.21$
mag (NGC 1023); Coma I group: $+0.08\pm0.15$ mag (NGC 4278) and
$-0.42\pm0.12$ mag (NGC 4494); Leo group: $-0.12\pm0.24$ mag (NGC
3377), $-0.11\pm0.24$ mag (NGC 3379), and $-0.18\pm0.24$ (NGC 3384).
At face value, assuming that in each of these 3 cases the spiral and
associated ellipticals are at exactly the same distance, the average
offset is $-0.16$ mag for the 3 cases (0.00 for NGC 1023, $-0.34$ for Coma
I, and $-0.13$ for Leo) in the sense that our elliptical distances are
too large.

This conclusion is too simplistic because we have not addressed
the assumption that the spiral galaxy is always at the same distance
as the comparison galaxy. For the Coma I group, for example, the GCLF
(\cite{Fleming95}) and the SBF (\cite{SP94}) methods concur that NGC 4494
is beyond NGC 4565, and curiously, the deviance for this galaxy is the
largest we see. If removed from the sample, the spiral--elliptical offset
is reduced to $-0.02$ mag, a level that is too small to consider reliable
within the assumptions. Until this problem can be addressed more
thoroughly, we apply no correction to our elliptical galaxy
distances.

\section{Details of New M87 PNLF Studies}

Bottinelli \mbox{\it et al.\ } (1991) and Tammann (1992) argued that PNLF distances to
Virgo ellipticals were underestimated because the luminosity depth of the PN
surveys was not adequate to sample beyond the brightest (0.5 mag) edge
of the PNLF.  Since that edge is nearly linear in the logarithmic PNLF,
the method becomes insensitive to distance modulus. In addition, those
authors challenged the PNLF distances on the basis that a shallow
survey of a large galaxy will suffer from a sample size bias. The sense
of this argument is that $N$ objects are more likely to be drawn from
the low probability bright tail of the large elliptical galaxy PN
sample than are $N$ objects from the smaller sample in M31's bulge.

While it is true that PN surveys must extend deep enough to sense the
curvature of the PNLF reliably with statistical methods, the required
depth is only 0.8 mag.  With the exception of NGC 4649 which
was observed under poor conditions, Jacoby, Ciardullo, \& Ford (1990)
estimated the depth of their surveys for 6 Virgo galaxies to be $\sim
1.0$ mag. Thus, it seemed unlikely that a serious systematic error was
contaminating those distances.

An independent assessment of the likelihood of a serious systematic
error is provided by recent Cepheid distances to Virgo galaxies.
The average PNLF distance to Virgo, based on 6
galaxies, is 15.3 Mpc (using the modern M31 distance and extinction for
the zero-point). This result agrees very well with the
\cite{Ferrarese96} distance of 15.8 Mpc to M100 based on HST Cepheids,
and the \cite{Pierce94} distance of 14.7 Mpc to NGC 4571 based on CFHT
Cepheids. In addition, \cite{Sandage96} reports HST Cepheid distances
to three near-Virgo galaxies: NGC 4496 at 16.6 Mpc, NGC 4536 at 16.6
Mpc, and NGC 4639 at 25.1 Mpc. Thus, four galaxies are reported in the
range 14.7 to 16.6 Mpc, and these are very comparable to the PNLF range
of distances (14.3 to 16.2 Mpc). One galaxy, though, is behind
all of these. It is unclear which, if any, of these spirals
represents the distances to the ellipticals,
but it is evident that most of the spirals (four out of five) have
distances that support the PNLF distances.

A direct resolution of the challenges to the PNLF distances lies in a
short observing project. Deep PN observations in M87 can push well into
the plateau region of the PNLF. Data were obtained in April 1995 with
the KPNO 4-m telescope to examine the claim that the earlier Virgo
data were not deep enough.  A total of 7 hours of on-line integration
were devoted to detecting fainter PN. This survey also extends to large
radial (10 arcmin) distances from M87's nucleus. The survey results are
described below.

\subsection[]{A Deep PNLF Distance Contradicts the Challenges}

Figure 3 shows the new PNLF for M87. A total of 320 PN were identified,
but many are fainter than the completeness limit. A total of 201 PN
are in the complete sample which extends $\sim 1.2$ mag down the PNLF.

\begin{figure}
\epsfysize 216pt \epsfbox{fig3.eps}
\caption{The 1995 PNLF for M87 based on 7 hours at the KPNO 4-m.
This PNLF extends 1.2 mag down the PNLF and clearly reaches beyond
the linear portion of the bright edge ($25.9 < m_{5007} < 26.3$)
of the PNLF. The scaled and shifted
PNLF from M31 is superposed to illustrate the agreement with the reference
galaxy's PNLF and demonstrates that the M87 PNLF is not a power law
over this regime. Note the single very luminous
object at $m_{5007}=25.6$ which was found first in the 1990 survey.
See Jacoby \mbox{\it et al.\ } (1996) for a discussion of what this luminous object may be.}
\end{figure}

Figure 4 shows a curious effect, though, which has not been fully
evaluated at the time of this conference. That is, when the sample of
PN is divided in half such that the PN drawn from M87's halo are
separated, the PNLF for the inner half (out to 4 arcmin in radius) are
systematically fainter by 0.3 mag than the halo sample. I return to
this point below, but note here that the more reliable distance is
derived from the inner sample because it is less likely to be contaminated
by intracluster PN. That distance, $14.4\pm1.3$ Mpc (on the modern
M31 distance scale), is nearly identical to the \cite{JCF90} value of
$14.9\pm1.2$ Mpc.

\begin{figure}
\epsfysize 216pt \epsfbox{fig4.eps}
\caption{The 1995 PNLF for M87 where open circles represent those PN
found in the inner 4 arcmin and triangles show those PN found beyond 4 arcmin
and out to 10 arcmin.}
\end{figure}

Thus, a deeper PNLF argues {\it against} the contentions of \cite{Bottinelli91}
and \cite{Tammann92} that the PNLF distance to Virgo is underestimated
as a consequence of inadequate survey depth.

\subsection[]{Sample Size Effects are Small}

With the new larger sample of PN, it is possible to investigate the
effects that different sample sizes have on the final PNLF distance.
Subsets of PN were drawn randomly from the sample of 201 PN.  Distances
for these subsamples of 20, 30, 40, 50, and 100 PN were derived
following our standard procedures and compared to the distance based on
201 PN. Figure 5 shows the magnitude of the effect of sample size
differences.

\begin{figure}
\epsfysize 216pt \epsfbox{fig5.eps}
\caption{Results of a Monte Carlo experiment to estimate the effect
of deriving distances with different sample sizes. The error bars
represent the scatter in the multiple attempts to derive distances with
a given number of PN in the sample.}
\end{figure}

In the worst case, for a sample of 20 PN, there is a slight tendency to
overestimate the distance to a galaxy by up to 3\%. Again, this
contradicts the challenges of \cite{Bottinelli91} and \cite{Tammann92}
who claim that our distances would be {\it underestimated}. The
reason that the effect is small is that the statistical process described
by \cite{Ciardullo89} is cognizant of the sample size and adjusts the
derived distances to the most likely one for a given sample. That is,
a statistical correction for sample size has always been applied to the
PNLF results.

\subsection[]{Bright PN in M87's Halo}

As noted above, M87's halo PNLF is $\sim 0.3$ mag brighter than
the central PNLF. We can all agree that the halo of M87
is {\it not} 15\% closer than its core! Thus, something must be artificially
enhancing the brightness of the PN in the outer regions. Since we have
not seen this effect in the two other samples that permit a similar
radial test (Cen A, \cite{Hui93}; NGC 4494, \cite{JCH96}), we
consider what could cause such an effect here. Five possibilities
come to mind:

\medskip
\begin{enumerate}
\renewcommand{\theenumi}{\arabic{enumi}}
\item Metallicity decrease in halo
\item Age decrease in halo
\item Dust near center
\item Instrumental effect
\item Intracluster PN contamination
\end{enumerate}
\medskip

The color gradient in ellipticals is such that the outer halos are
bluer than the inner regions. The first 3 possibilities above
have been suggested as possible causes of gradients. The first,
a metallicity decrease, reduces the luminosity of PN, in contradiction
to the observed effect. The second, the presence of a young population
can enhance the PNLF luminosity if the ages of the halo stars are $<$0.5 Gyrs,
provided the central population is $>$3 Gyrs. We can
neither dismiss nor confirm this possibility.

The third option, dust in the central regions, has been suggested by
Wise \& Silva (1996), \cite{GD95}, and \cite{WTC92} to explain the
color gradients.  To explain the PNLF enhancement, central reddening
of E(B--V)$\sim$0.06 is needed. Again, this is plausible.

Although the effect has not been seen before, a brightness enhancement
in the halo could be caused by instrumental effects. Thus far,
flat-fielding errors and filter transmission variations have been dismissed.

The fifth option initially seems highly speculative. The key point is
that a large population of intracluster stars exists and can produce PN
having a range of distances representing the full depth extent of the
Virgo cluster. Thus, some PN may be foreground to M87
and appear brighter than M87's own PN, while other PN
will be in the background and be lost in the faint end
of the PNLF.  Since the number of intracluster PN found is proportional
to the area surveyed, there is a survey bias against finding foreground
PN in the smaller central region, while simultaneously, there is a bias
in favor of finding true M87 PN at the center where the stellar density
is high.

What is the likelihood of finding intracluster PN? \cite{Arnabaldi}
found that 3 of 19 PN in the \cite{JCF90} NGC 4406 sample are intracluster.
Since that sample is biased against intracluster PN due to velocity
rejection in the survey filter (NGC 4406 was sampled at its systemic
velocity of $-220$ km/s, or about 1500 km/s from the Virgo systemic velocity),
there are likely to be many intracluster PN. \cite{Arnabaldi} discuss
the intracluster population in detail. To zeroth order, the number density
of intracluster PN is not a problem.

A definitive source of the enhanced halo PNLF in M87 is not possible
at this time. Kinematics can be used to investigate the likelihood
that intracluster PN are contaminating the halo sample.
The other likely causes, extinction and very young populations, seem
less secure at this time because we don't know that they exist while
we do know that intracluster PN do exist.

\section[]{Can M87 Really Be At 22 Mpc?}

Sandage and Tammann (e.g., \cite{Sandage90}) have long argued that for
values of H$_0 < $55, the Virgo Cluster must be at a distance of about
22 Mpc since the cosmic velocity of Virgo is $\sim 1170$ km/s.
Additionally, M87 is usually claimed to be representative of the Virgo
distance since its recessional velocity is close to the Virgo average,
it is projected near the center of the cluster, and it appears to be
at the bottom of the potential well as evidenced by
an x-ray cloud typical of x-ray clusters.

A reasonable question to ask is ``what must one give up to
derive a PNLF distance of 22 Mpc for M87?'' The simple answer is
that the constancy of the maximum PN luminosity must be abandoned. And,
it must be abandoned in a big way -- by $\sim$0.9 mag.

The astrophysical question becomes ``can the PNLF be enhanced by 0.9 mag,
and if so, what are the observable manifestations?'' The latter turn out to
be severe and obvious. Models of the PNLF by \cite{Jacoby89}, \cite{Mendez93},
and \cite{Stanghellini95} all show that a shift in the PNLF as large
as 0.9 mag is almost impossible because there is inadequate UV flux
being radiated by any reasonable collection of central stars.  A 0.9
mag enhancement demands that the central stars be very massive,
originating from 4--5 M$_\odot$ progenitors. A population this young
($<100$ Myrs) should be evident.

The presence of a very young population in ellipticals is not
unreasonable; many ellipticals exhibit evidence for gas and dust
(\cite{Goudfrooij}). Usually, though, the mass fraction of the
elliptical involved in the young population is tiny. To enhance the
PNLF by 0.9 mag, though, requires that $\sim$40\% of the total
luminosity of the galaxy comes from the young population.  (The 40\%
value derives from the PN production rate for bright PN in M87.  This
rate is 40\% of the maximum value that is attainable under the
assumption that all stars produce PN. Since some stars may not produce
PN, the actual fraction of luminosity from young stars could be higher.)

A trivial observational test for a young population is galaxy color.
A young population ($<$100 Myr) has a color (B--V) $<$ 0.0 while an old
population ($>$11 Gyr) has a color (B--V) $\sim$ 1.0. Thus a mix where
40\% of the luminosity derives from a young population has a color
(B--V) $\sim$ 0.6. This color is the direct consequence of insisting that
M87 be at 22 Mpc, combined with stellar evolution predictions
(e.g., \cite{Ciotti}); \cite{Vass})
for the initial-to-final mass relation that was devised to explain
the observed relation (\cite{WK83}).

Comparing the predicted color of 0.6 to the observed (B--V) $\sim$ 1.0
(which depends slightly on radial position), it is evident that a sizable
young population is untenable. Thus, we are forced to accept that either
M87 is much closer than 22 Mpc in order to alleviate the pressure to
enhance the PNLF luminosity, or that the initial-to-final mass relation
predicted by stellar evolution models and observed in young Galactic clusters
is seriously flawed in such a way that old stars can produce massive
central stars. Since there is no evidence for the latter, the more
likely conclusion is that the existence of PN in M87 at the observed
apparent magnitudes demands a distance much smaller than 22 Mpc.
(It is possible to push the distance as large as 17 Mpc
before the implied population color becomes inconsistently blue.)

\section{Some Important Distances}

Based on the comparisons between PNLF and Cepheid distances, between
PNLF and SBF (\cite{CJT93}) distances, and between PNLF spiral and PNLF
elliptical distances within a cluster, there is repeated evidence that
PNLF distances are accurate to $\sim$8\%.  Similar claims are being
made for SN Ia.  Yet, PNLF distances suggest values for H$_0$ near 75
km/s/Mpc while SN Ia distances are between 55 and 67 km/s/Mpc
(\cite{Sandage96}; Riess \mbox{\it et al.\ } (1996); Hamuy \mbox{\it et al.\ } (1996)) and both
methods are based on Cepheids for calibration. While an H$_0$ of 67 is
within the combined error budget of the methods, an H$_0$ of 55
stresses the PNLF method (and others) severely.

The PNLF and SN Ia distance comparison shown in Table 2 for galaxies in
common in the Virgo and Fornax clusters summarizes the issue.  For
Virgo, these include a single SN Ia in each of NGC 4374, NGC 4382, and
NGC 4486, although none was observed very well.  Although SN1991bg in NGC 4374
was observed very well, we reject it for being anomalously
underluminous (\cite{Phillips93}; \cite{Leib93}).  For Fornax, 2 SN Ia
occurred in NGC 1316, and both were reasonably well observed
(\cite{Hamuy91}).  Distances are given as the average of the 3 galaxies
in Virgo.

The SN Ia zero-point is set using the well observed SN Ia SN1981B and
SN1990N. No others (SN1895B, SN1937C, SN1960F, SN1972E) were adequately
observed for their peak magnitudes and decline rates to be determined
at the same level of reliability as these 2 primary calibrators. From
\cite{Sandage96}, we have M$_{\rm B}^{max} = -19.3\pm0.08$ for SN1981B
and SN1990N. The apparent peak magnitudes for the Virgo and Fornax SN
Ia are taken from \cite{Leib91}.

\begin{table}
\begin{center}
\begin{tabular}{lccc}
Cluster & PNLF Distance & SN Ia Distance & Discrepancy     \\
        &   Mpc         &     Mpc        &   in $\sigma$'s \\
        &               &                &                 \\
Virgo   & $15.1\pm0.9$  &  $17.4\pm2.6$  &    0.8$\sigma$  \\
Fornax  & $17.7\pm0.5$  &  $22.9\pm1.5$  &    3.3$\sigma$  \\
\end{tabular}
\end{center}
\caption{Comparison Between PNLF and SN Ia Distances to Virgo and Fornax}
\end{table}

Evidently, the discrepancy in Virgo is not significant while the
discrepancy in Fornax is significant. Thus, either the PNLF distance to
Fornax is incorrect (although 2 other galaxies yield the same distance
and the Cepheid distance is nearly identical), or the 2 well observed SN
Ia in Fornax are underluminous by 0.5 mag.  There is evidence that the
Fornax SN Ia are, in fact, somewhat fast declining ones, but this can
only explain about 0.15 mag of the discrepancy. If the decline rate is
considered, using the \cite{Hamuy96} slope of 0.78, then the
discrepancy is reduced to a 2.2$\sigma$ event, which still is
marginally significant.  Because the discrepancies between the Fornax
SN Ia distance and Fornax distances from the HST Cepheids, SBF, and
PNLF are all very similar, the SN Ia Fornax inconsistency
cannot be solved by appealing solely to errors in the PNLF method.
For now, this issue remains open.

\section{Conclusions}

No distance indicator is perfect. Confidence in the results, though,
is built by extensive testing against other distance indicators,
and through the development of the astrophysical theory supporting
the indicator's use.

The PNLF technique has been well tested against Cepheids, SBF, and against
itself in clusters. Constancy of the bright end of the PNLF has been
easy to reproduce theoretically if the ages of the stars producing PN are
within a plausibly wide range. Thus, PNLF distances
appear reliable at this time. The main points of this paper are:

\medskip
\begin{enumerate}
\renewcommand{\theenumi}{\arabic{enumi}}
\item PNLF distances to spirals are consistent with Cepheid distances.
\item There are now 7 PNLF calibrators; their dispersion is 8\%.
\item Good agreement with Cepheid and SBF distances must mean that
systematic errors (due to age, metallicity, extinction, PN sample sizes,
PNLF methodology) are under control to the limits of the deviations
of these methods (typically, 8\%).
\item The M87 deep survey demonstrates that the bright end of the PNLF
is not a power law as some have maintained and that PNLF distances are
not sensitive to the number of PN in the sample.
\item Both the Virgo and Fornax distances agree with the Cepheid-based
distances, but the Fornax SN Ia distance is $\sim$25\% larger. The Virgo
SN Ia distance is not discrepant, within the errors.
\end{enumerate}
\medskip

In terms of the future application of the PNLF, it is somewhat costly
in telescope time (exposure time is proportional to distance$^4$).
So, while the method yields results as reliable as the Cepheids or SBF,
its application may be reserved for galaxies that present problems
for other methods. For example, it can be used to measure distances
to SN Ia hosts such as Sa/S0 galaxies where Cepheids won't be found easily
and where galaxy structure compromises the SBF.

\medskip
\medskip

\noindent
{\bf Acknowledgements}
\begin{acknowledgments}
This paper derives from other papers in press. The discussions
on spirals was based on information provided by John Feldmeier; preliminary
results on the deep M87 survey were provided by Robin Ciardullo.
I am grateful to John Graham who gave me his HST Cepheid distance to
NGC 3351 prior to publication.
\end{acknowledgments}

\clearpage


Where to look for radiatively inefficient accretion flows (and find them)

Abstract: We have studied the nuclear emission detected in high-resolution radio and optical data of carefully selected samples of low luminosity AGN (LLAGN) in the local universe. When the Eddington ratio is plotted against the nuclear ``radio-loudness'' parameter, sources divide according to their physical properties. It is thus possible to disentangle between nuclear jets and accretion disks of different radiative efficiencies. If this simple interpretation is correct, we now have a powerful tool to investigate the nature of the nuclear radiation, and identify radiatively inefficient accretion flows (RIAF) candidates. Our results show that the best chance of investigating RIAF processes in the IR-to-UV spectral region is to observe (at the resolution provided by HST) the nuclei of unobscured Seyferts of the lowest luminosity, as well as a sub-class of LINERs. In all other objects other radiation processes dominate. In a sample of 132 LLAGN we identify 8 objects in which we predict the radiation from a RIAF can be directly detected.

\section{Introduction}

Models  of   radiatively  inefficient  accretion   flows  (RIAF)  were
originally developed  to reproduce the low state  of activity observed
in Cyg  X-1 \citep{ichimaru77}. They were later  discussed in relation
to AGN \citep{rees82,fabianrees95} in order to account for the lack of
activity  not  even  at   the  levels  predicted  by  Bondi  accretion
\citep{bondi52}.  The models are based  on the possibility that at low
accretion rates  a physical solution  for the accreting  flow involves
low  density,  hot, optically  thin  gas,  which  is also  radiatively
inefficient.   This regime would  occur instead  of that  of optically
thick, geometrically  thin disk accretion  (`Shakura--Sunyaev disks').
While  the theoretical effort  to characterize  accretion at  very low
rates   has    been   impressive   \citep[and    many   papers   since
then]{narayanyi},   direct observations  of RIAFs
still lack behind.

Low luminosity AGN  (LLAGN) are the class of  extragalactic objects in
which RIAFs  are most likely  to be at  work.  They show some  level of
activity  (compact radio emission  with high  brightness temperatures,
line emission,  nuclear X-ray emission)  which in the majority  of the
objects cannot be  explained in terms of processes  other than a faint
active  nucleus  powered  by  the  central  supermassive  black  hole.
However, because of their  intrinsic low luminosities, it is extremely
difficult to  study the emission from  RIAFs in LLAGN.   This holds in
particular  for  the wavelengths  at  which  the spectral  differences
between RIAF's and ``standard''  accretion disk models are expected to
differ the most, i.e. the IR-to-UV spectral region.  RIAFs should lack
both the "big blue-bump" and  the IR (reprocessed) bump, which instead
characterize  optically  thick,   geometrically  thin  accretion  disk
emission and the surrounding heated dust \citep[e.g][]{elvis94}.

While it is now possible to routinely disentangle the nuclear emission
from      that      of      the      host      thanks      to      HST
\citep[e.g.][]{hopeng,pap1,papllagn,sandrobarbara05},     the     RIAF
emission cannot be seen when  it is swamped by other nuclear radiation
processes, or  obscured by dust.  M~87  \citep{dimatteo03} is a
clear example of  that, where the nuclear emission  at all wavelengths
is  dominated by  non-thermal  synchrotron radiation.   \citet{Maoz05}
have  shown  that  among  a  sample  of 17  LLAGN,  15  of  them  show
variability over a timescale of  a few months, which demonstrate their
non--stellar origin.   However, it was  still unclear from  their work
whether  the nuclear radiation  is from  a jet  or from  the accretion
flow.

Recently,   among  LLAGN   which  show   very  low   Eddington  ratios
$L_{bol}/L_{Edd} << 10^{-3}$, \citet{papllagn} have found that a class
of  LLAGN,  mainly composed  by  LINERs  and low-luminosity  Seyfert~1
galaxies, show faint unresolved optical  nuclei in HST images that may
be  interpreted  as  direct  radiation  from  a  very  low  efficiency
accretion flow.  In fact, that  work shows that when the radio-optical
properties of LLAGN are considered, Seyfert, LINERs and low luminosity
(FR~I) radio  galaxies separate  into different regions  of diagnostic
planes,  according  to  the  properties  of  their  nuclei.   If  this
interpretation is correct, we now have a powerful tool to identify the
nature   of   the    nuclear   radiation   (i.e.    jet-dominated   or
accretion-dominated).   This picture  appears to  be confirmed  by the
nuclear SED of the low luminosity Seyfert galaxy NGC~4565, which shows
very  low obscuration  and  lacks the  typical  signature of  standard
accretion \citep{pap4565}.

In this Letter, by expanding the  sample of LLAGN, we discuss what is,
in our view, the most efficient way of identifying a sample of objects
that are  most likely to host  a RIAF and in  which IR-to-UV radiation
from the accretion process can be directly studied.

\section{Is the ``fundamental plane of black hole activity?''
a diagnostic plane?}

LLAGN have been  found to lie on the  so-called ``fundamental plane of
black hole activity''  \citep{merloni03,falcke04}, which has the great
merit of attempting  a unification of all sources  associated to black
holes, over  a large range  of masses and luminosities,  from Galactic
sources to powerful  quasars.  But the origin of  such a ``fundamental
plane''  and its  relationship  with the  physical  properties of  the
source is still a matter of debate \cite[e.g][]{bregman05,koerding06}.
Here  we want to  test whether  such a  plane can  be easily  used for
discriminating  between different  radiation  processes associated  to
inflows and/or  outflows around black  holes, with the aim  of finding
RIAF candidates.

We  perform a  simple test,  to check  whether objects  with different
emission  processes and of  different physical  origin are  located in
distinct regions.  The ``fundamental plane'' is based on the radio and
X-ray  luminosity, combined  with the  black hole  mass, which  is the
ultimate ``energy source'' for the radiation we observe in the form of
inflows (accretion)  or outflows (jets).   However, if we  check where
other sources  of radio  and X-rays that  are not associated  to black
holes are  located in  this plane, we  find something  intriguing.  We
take the radio and X-ray fluxes of solar system objects, together with
their masses,  and we over plot  them onto the  ``fundamental plane of
black hole  activity''.  Unexpectedly, the Sun, the  Moon, Jupiter and
Saturn do fall on the same  plane, although they are tens of orders of
magnitudes fainter  and less massive  than all ``black  hole'' sources
(Fig.~\ref{funda}).  In  fact, the radiation  processes that originate
the emission from  the solar system objects are  variegated.  While it
is well  known that we observe  thermal radio emission  from the quiet
Sun,  Saturn,  the  Moon  and  Jupiter (below  3cm),  scattered  solar
radiation  is  the origin  for  the X-ray  emission  of  the Moon  and
probably Saturn, high ionization  emission lines dominate in the X-ray
Sun                and                Jupiter               \citep[see
e.g.][]{manson77,kraus86,schmitt91,ness04,branduardi06}.

This implies  that such  a ``fundamental plane''  might be  telling us
something profound  associated with the emitting power  of all sources
of  radio  and  X-ray  radiation,  but  cannot be  easily  used  as  a
diagnostic to  discriminate between different  radiation processes. On
the other  hand, the  apparent correlation might  only be a  result of
mere plotting luminosity  vs. luminosity, and ``artificially'' rescale
the quantities  with the mass  \citep{bregman05}. We speculate  that a
closer look at  the scatter of the ``correlation''  might lead to more
physical  information than just  trying to  reproduce it  assuming ``a
priori''  different models.   However,  as we  show  in the  following
section, in  order to discriminate different  radiation processes (and
eventually  find  RIAF  candidates)  we should  use  diagnostics  that
enhance the differences.

\begin{figure}
\epsscale{1.2} \plotone{f1.ps}
\caption{The so-called  fundamental plane  of black hole  activity, to
which   we   over-plotted    solar   system   sources   \citep[adapted
from][]{merloni03}.  The location of  NGC~4565 (which  we refer  to in
Sect.~\ref{fund}) along the correlation is also shown.}
\label{funda}
\end{figure}

\section{How to find RIAFs: nuclear radio loudness and Eddington ratio}
\label{fund}

In order  to explore the nuclear  properties of all kinds  of LLAGN in
the local universe  with the aim of searching  for RIAF candidates, we
consider the following samples:

{\bf 1)} 33  FR~I Radio Galaxies (i.e. low  luminosity radio galaxies)
from the 3CR catalog (radio selected) \citep{pap1};

{\bf 2)} 25   Seyferts from the Palomar Survey  of nearby galaxies
and from the CfA sample, clearly limiting ourselves to those belonging
to Type~1, i.e.  unobscured, \citep[optically selected][]{hopeng};

{\bf 3)}  a complete sample  of 21 LINERs  from the Palomar  Survey of
nearby galaxies \citep[optically selected][]{ho97};

{\bf 4)} 51  nearby early-type galaxies (E+S0) with  radio emission $>
1$ mJy  at 5 GHz  (optical $+$ radio selection)  \citep[and references
therein]{sandrobarbara05}. The  large majority of the  galaxies in the
sample are spectrally classified as either LINER or Seyfert.

{\bf 5)} the 12 broad-line radio galaxies with  $z<0.3$ included the 3CR
catalog \citep[and references therein]{pap4}.

Samples 1, 2  and 3 have been examined  in detail in \citep{papllagn}.
Note  that  being  selected  according to  different  criteria,  these
objects do not constitute  a complete sample.  However, they perfectly
represent the  overall properties  of low power  active nuclei  in the
local universe, where RIAFs are most likely to be found. The sample of
nearby  ellipticals   partially  overlaps  with  samples   1)  2)  and
3). However, there are only 10  objects in common, so the total number
of objects considered is 132.

In \citet{papllagn} we showed that the nuclear properties of LLAGN are
best  understood  when  the   ratio  between  optical  luminosity  and
Eddington  luminosity $L_o/L_{Edd}$ is  plotted against  the ``nuclear
radio-loudness parameter'' $R$, defined  as the ratio between the {\it
nuclear} radio (core) luminosity at 5GHz $(F_{5GHz})$ as measured from
high-resolution (VLA or higher  resolution) data and the {\it nuclear}
optical luminosity  $L_o$ as measured  from HST images.  The  power of
this  diagnostic plane  resides  in the  fact  that we  only need  two
measured  quantities, plus  an estimate  of  the black  hole mass,  to
clearly separate sources of  different physical origin. This allows us
to  discriminate  between  jets   and  accretion  disks  of  different
radiative efficiency (and/or different accretion rates).

\begin{figure}
\epsscale{1.2} \plotone{f2.ps}
\caption{The diagnostic plane for  the nuclear emission of LLAGN. The
ratio  between  the  optical  nuclear  luminosity  and  the  Eddington
luminosity is plotted against  the radio loudness parameters.  Symbols
are as follows: Seyfert~1 are  marked as triangles, pentagons are BLRG
and QSO  from the  3CR catalog, are  FR~I radio galaxies,  squares are
LINERs,  hexagons  are  early-type  ``core-galaxies'' and  stars  are
early-type ``power-law''  galaxies. The dashed lines  are only plotted
to guide  the eye in  discriminating between different  ``species'' of
nuclei.}
\label{ledd}
\end{figure}

Since here  we are  looking for object  harboring RIAF  candidates, in
Fig~\ref{ledd} we  plot galaxies for which the  ``active nucleus'' has
been detected  both in the radio  and the optical (or  near IR) bands.
The lack of detection of  compact radio emission is considered as lack
of evidence for AGN activity. On the other hand, the lack of detection
of the optical nucleus would set a double upper limit in the plot and,
most  importantly,   the  location  of   the  data  point   along  the
radio-loudness axis would  remain undetermined.  Therefore, we discard
such objects.   We also only plot  the objects for which  a black hole
mass estimate  is available, either  through gas dynamics or  by using
the relation with the stellar velocity dispersion.

Sources  separate into  4  quadrants: Seyfert~1s  occupy the  top-left
quadrant, FR~I radio galaxies are in the bottom-right quadrant.  These
two classes  define the regions characterized  by ``radio-quiet'' high
radiative   efficiency  accretion   and   ``radio-loud''jet  emission,
associated to low radiative efficient accretion, respectively.  LINERs
split  into   two  sub-samples:  both   are  in  the  region   of  low
$L_o/L_{Edd}$, some  of them  being radio quiet  and some  radio loud.
Nuclei  of ellipticals  separate according  to the  properties  of the
radial  brightness profile  of the  host  \citep{sandrobarbara06}, and
behave as the LINERs: core-galaxies are radio-loud, power-law galaxies
are  radio-quiet  \citep[for  a   definition  of  core  and  power-law
galaxies, see  e.g.][]{faber97}.  Broad line  radio galaxies, instead,
are found  at the top  of the plot,  at high values  of $L_o/L_{Edd}$,
while their location on the radio loudness axis most likely depends on
the relative  importance of the jet  and disk emission  in the optical
band  (the radio  core being  dominated  by radiation  from the  jet).
Although it is not easy  to identify any ``bi-modality'' when plotting
all the  samples together as  done in Fig.~\ref{ledd}, we  draw dashed
lines  just  to  guide  the  eye in  discriminating  between  physical
processes  that are  clearly different.   Clearly, there  is a  bit of
ambiguity for the nuclei that fall close to those ``dividing lines''.

Let us  now focus on the  right hand side  of the panel. Here  we find
jet-dominated  nuclei  (LLRG and  a  subsample  of  LINERs). In  these
objects,  both the  radio and  the  optical radiation  is most  likely
synchrotron  emission  from  the base  of  the  jet,  as it  has  been
established      for      low      luminosity      radio      galaxies
\citep[e.g.][]{pap1,gijs}.  The  radiation from the  accretion flow is
swamped by the jet emission,  and cannot be studied directly (at least
in the observing bands considered  here).  We should point out that in
this   case   we   have   an   ``upper  limit''   to   the   Eddington
ratio. Nevertheless, for  several objects, this upper limit  is as low
as $L_o/L_{Edd} \sim 10^{-8}$. On the other hand, none of the detected
objects on the left side of  the plot reaches such low values.  On the
other hand,  these extremely low  efficient accreting black  holes are
still capable of producing a  ``jet''.  This means that although RIAFs
(or whatever  these objects are!)   are inefficient from the  point of
view of producing radiation, under certain circumstances (i.e. for the
``radio-loud''  nuclei)  they  can  be  very  efficient  in  producing
outflows  or  even  (relativistic)  jets,  as  in  the  case  of  LLRG
\citep[see e.g.][]{rees82}.   In \citet{papllagn} we  pointed out that
``radio-loud'' (i.e.  jet dominated) nuclei tend to be associated with
the most massive  black holes.  Thus it is  tempting to speculate that
the enhanced efficiency in  producing powerful jets is somehow related
to the  higher black hole mass.  An alternative picture is  that it is
the   black    hole   spin   that   efficiently    powers   the   jets
\citep{blandfordsaasfee,sikora06},  or possibly  a combination  of the
two proposed scenarios.

Let us now  focus on the left  side of the plot, which  is central for
the   purposes  of   this   Letter.   As   discussed   in  detail   in
\citet{papllagn}, in that region we  find nuclei that show a low value
of  $L_r/L_o$, i.e.   an optical  excess with  respect to  the optical
counterpart  of the  radio (synchrotron)  radiation.   For Seyfert~1s,
such an  optical excess  is readily interpreted  as emission  from the
accretion disk.  A  few bright Seyferts for which  the nuclear SED has
been     derived     support     this    interpretation     \citep[see
e.g.][]{alonso03,pap4565}.   For the brightest  nuclei, the  limits on
the Eddington  ratio are  still compatible with  radiatively efficient
accretion.   Furthermore, we should  stress that  for those  objects a
bolometric correction of  a factor $\sim 10$ \citep{elvis94,marconi04}
should be  performed, because we know  their SED and we  know that the
optical $R$ or  $V$ band do not correspond to  the peak of radiatively
efficient accretion energy output.  On  the other hand, for objects of
the  lowest  Eddington  ratios,  we  don't know  what  the  bolometric
correction should  be, but  a factor of  10 (or even  slightly higher)
would not  change our conclusions.   Objects with Eddington  ratios as
low  as  $10^{-6}$  cannot  be  anyhow  reconciled  with  ``standard''
accretion  disk models.  Therefore,  they are  the best  candidates to
host radiatively inefficient accretion disks.

In  fact,  we have  shown  in  \citet{pap4565}  that NGC~4565,  a  low
luminosity  Seyfert (not  included in  the sample  considered  in this
Letter) that resides in the lower-left region of the diagnostic plane,
shows an unusual SED which  does not display the typical signatures of
radiatively efficient accretion disks (i.e.  lacks both the UV and the
IR bumps).  This strongly supports our prediction that objects in that
part of the plot are not ``standard'' disks.

Since we want  to identify the objects in which  radiation from a RIAF
can be detected and studied  in the IR-to-UV spectral region, in order
to derive  the SED and  set constraints to  the models, our  best RIAF
candidates are  those in which a  nuclear source has  been detected in
archival  HST images.   This reduces  the number  of objects  to eight
candidates, plus  NGC~4565 which  has been already  studied elsewhere.
Summarizing, we predict that RIAFs  can be detected and studied in the
following  objects:  M~81,  NGC~3245,  NGC~3414,  NGC~3718,  NGC~3998,
NGC~4143, NGC~4203, NGC~4565 and NGC~4736.

\section{Summary and conclusions}

By combining the Eddington  ratio as measured from detected unresolved
nuclei in HST images  and the nuclear ``radio-loudness'' parameter, we
have found  a straightforward way of  distinguishing radiation emitted
by  the   accretion  process  from   jets  in  LLAGN  and   find  RIAF
candidates. At the end of  the selection process, which started from a
large sample of 132 AGN in  the local universe, we are left with eight
RIAF candidates for  which the nuclear source is  detected in archival
HST images. These are the best targets for follow-up studies, not only
in  the IR-to-UV  spectral region,  but also  in other  bands  such as
far-IR, mm, and  X-rays. Studying those objects in  more detail is the
best  chance to provide  constraints for  RIAF models  associated with
supermassive black holes.

The method to find RIAF candidates outlined in this Letter is based on
the radio and optical emission.  However, it is possible to extend our
ability of discovering  new candidates using other bands  in which the
nuclear luminosity can be estimated. For example, nuclear IR radiation
between $\sim 10$ and $\sim 20 \mu$m can be isolated in nearby objects
using  8m  class  telescopes,  as done  by  \citet{whysong}.   Another
possible way  of estimating  the nuclear flux  is to use  the emission
line strength, assuming a correlation between these two quantities, as
\citet[e.g.][]{sandrobarbara06}  pointed  out  that  when  the  [OIII]
luminosity is used instead of  the optical nuclear luminosity from HST
images,  similar   results  are  achieved.   However,   in  all  these
alternative methods  we are  not sure that  the optical source  can be
detected in order to derive the SED and constrain the models.

We should also  stress that the diagnostic plane  discussed here holds
for LLAGN  in the local universe.   When high-z more  powerful AGN are
over-plotted,  cosmic evolution  may play  a role  in  determining the
position of  the nuclei  in the  plane. For example,  since it  is not
clear  whether the  local  relations between  properties  of the  host
galaxy and the supermassive black  hole mass can be easily extended to
high-redshift galaxies, the estimate of the black hole mass may change
substantially  with cosmic  time.  Furthermore,  the  overall spectral
properties of jets and disks, and their relative importance in the SED
may  change as  well, and  all this  should be  carefully investigated
before including ``all known AGNs'' in the plane and draw conclusions.
However, RIAFs are  so faint and so difficult to  detect, that at this
stage it is best to look for them in local objects.

\acknowledgments M.C.   wishes to thank  A.  Capetti, A.   Celotti, D.
Macchetto,  R.   Gilli and  W.B.  Sparks  for  useful discussions  and
insightful comments.


Opposite magnetic polarity of two photospheric lines in single spectrum of the quiet Sun

Abstract: We study the structure of the photospheric magnetic field of the quiet Sun by investigating weak spectro-polarimetric signals. We took a sequence of Stokes spectra of the Fe I 630.15 nm and 630.25 nm lines in a region of quiet Sun near the disk center, using the POLIS spectro-polarimeter at the German VTT on Tenerife. The line cores of these two lines form at different heights in the atmosphere. The 3$\sigma$ noise level of the data is about 1.8 $\times 10^{-3} I_{c}$. We present co-temporal and co-spatial Stokes-$V$ profiles of the Fe I 630 nm line pair, where the two lines show opposite polarities in a single spectrum. We compute synthetic line profiles and reproduce these spectra with a two-component model atmosphere: a non-magnetic component and a magnetic component. The magnetic component consists of two magnetic layers with opposite polarity: the upper one moves upwards while the lower one moves downward. In-between, there is a region of enhanced temperature. The Stokes-$V$ line pair of opposite polarity in a single spectrum can be understood as a magnetic reconnection event in the solar photosphere. We demonstrate that such a scenario is realistic, but the solution may not be unique.

\section{Introduction}

Stokes polarimetry provides detailed information about the magnetic field  and its interaction
with the plasma in the solar photosphere and chromosphere~\citep{stenflo_94,landi_landolfi}.
While most of the observed Stokes-$V$ profiles in active regions and the network
are close to antisymmetric with a low degree of asymmetry, abnormal and strongly asymmetric $V$ profiles are
common in the inter-network~\citep{sigwarth01,Lites_02}.
The classification of \cite{sigwarth_etal_99} presents abnormal $V$ profiles
with a single lobe, two lobes with identical polarities, and
$Q$-like, or pathological profiles with four or more lobes.
There are indications from observations and 3-D simulations that the
degree of asymmetry and the fraction of abnormal $V$ profiles
increase with decreasing magnetic flux~\citep{sigwarth_etal_99, khomenko_etal_05}.
While most of the pathological profiles can be reconstructed with models consisting of two or
more magnetic components, \cite{grossman_etal_00} and \cite{steiner_2000} have shown that
one-component models can also account for a large variety of $V$ profiles provided that magnetic,
velocity, and temperature gradients are large enough.

An opposite polarity (OP) Stokes-$V$ profile is a set of two $V$ profiles of two
different spectral lines, recorded in a strictly co-temporal and co-spatial observation,
that shows different polarities.
\cite{almeida_cerdena_kneer03} report opposite polarity $V$ profiles
between the visible~(630\,nm) and the infrared~(1.56\,$\mu$m) neutral iron lines observed at
 the French-Italian solar telescope TH\'EMIS and the
German Vacuum Tower Telescope (VTT), respectively.
This work was questioned later by \cite{khomenko_etal_05}, who showed that different seeing
conditions for the two data sets can spuriously produce OP profiles.

In this paper, we present observations of a quiet Sun region close to the disk center with
the POlarimetric Littrow Spectrograph~\citep[POLIS,][]{schmidt03,beck05b} that show
a few OP Stokes-$V$ profiles.
Each set of OP profiles consists of the two iron lines at $\lambda\lambda$ 630.15 and 630.25\,nm
that are part of  a {\it single spectrum} recorded strictly simultaneously.
The formation heights of these two lines span different layers in
the photosphere~\citep[][and references therein]{cabrera_etal_05,khomenko_collados}.
We study one set of these profiles in detail and argue that it hints at
a magnetic reconnection event in the solar photosphere.

\begin{figure*}
\sidecaption
\resizebox{10.8cm}{!}{
\includegraphics*{7441fig1.eps}}
\caption[]{Co-spatial and co-temporal data of the speckle
reconstructed continuum image at 500\,nm (left), $V_{\mathrm{tot}}$ map (middle), and
Stokes-$I$ continuum map (right). Contours of the $V_{\mathrm{tot}}$ map are overplotted on all panels.
The position of the opposite polarity profile is marked with a rectangle, which also shows
the POLIS pixel size. The horizontal is the scan direction and the vertical is the slit direction.
In the $V_{\mathrm{tot}}$ map, white and black indicate negative and positive
polarities, respectively (for definition of positive/negative polarity, see Fig.~\ref{fig:op}).
Tickmarks indicate arcsec.}
\label{fig:maps}
\end{figure*}

\section{Observations and data reduction}
A sequence of spectra taken in a quiet Sun region close to the disk center
($\cos\theta$\,=\,0.99), was observed with the VTT in Tenerife, July 07, 2006.
The seeing was good and stable during the observation.
The Kiepenheuer Adaptive Optics System  was used for maximum spatial
resolution and image stability~\citep{luhe_etal_03}.
Each of the 38 scans consists of 16 slit positions.
The scanning step size and spatial sampling along the slit were  0.48 and 0.29 arcsec, respectively.
The scanning cadence was about 97~s.

Full Stokes profiles of the neutral iron lines at 630.15 and 630.25\,nm
and the Stokes-$I$ profile of the \ion{Ca}{ii}\,H line
were observed strictly simultaneously with the red (630\,nm) and blue (396.8\,nm)  channels of POLIS.
The average continuum contrast (rms) of the POLIS Stokes-$I$ maps
is $\approx\,3.4\,\%$ of $I_{\mathrm{c}}$.
An absolute velocity calibration was performed using the telluric O$_2$ line
at 630.20\,nm \citep{reza_etal_1}. The spectro-polarimetric data of the red channel
were corrected for instrumental effects and telescope
polarization with the procedures described by \cite{beck05b, beck05a}.
The rms noise level of the Stokes parameters in the continuum
is  $\sigma$\,=\,6.0\,$\times 10^{-4}$\,$I_{\mathrm{c}}$.

Simultaneously, a continuum speckle channel in POLIS recorded a larger field of view at 500\,nm.
The speckle reconstruction was performed using the Kiepenheuer-Institut Speckle
Imaging Package~\citep{mikurda_luehe,fwoeger_phd}.
The spatial resolution of the reconstructed image  is about 0.3~arcsec (cf. Fig.~\ref{fig:maps}, left panel).
We used the POLIS intensity map and the reconstructed image to align the data.

We define the signed $V_{\mathrm{tot}}$ as follows:
\begin{equation}
V_{\mathrm{tot}} = \frac{\int^{\lambda_0}_{\lambda_{\mathrm b}} V(\lambda)\, \mathrm{d}\,\lambda
- \int^{\lambda_{\mathrm r}}_{\lambda_0} V(\lambda)\,
\mathrm{d}\,\lambda}{I_{\mathrm{c}}\,(\lambda_{\mathrm r} - \lambda_{\mathrm b})}{\rm ,}
\end{equation}
where $\lambda_0$ is the zero-crossing wavelength of the Stokes-$V$ profile
and $\lambda_{\mathrm r}$ and $\lambda_{\mathrm b}$ denote  fixed wavelengths in the red and blue
continuum of the lines~\citep{lites_etal_99}.
The continuum  speckle image and the POLIS Stokes-$I$ map are shown in the left
and right panels of Fig.~\ref{fig:maps}, respectively.
The middle panel of Fig.~\ref{fig:maps} shows the $V_{\mathrm{tot}}$ map in which
the position of the OP profile is marked with a rectangle that also shows the POLIS pixel size.
It is located between two patches of opposite polarities.
The OP profile (exposure time 5\,s) was recorded within the time window of 15~s used for the speckle
burst~(left panel, Fig.~\ref{fig:maps}).
This image
and the continuum map (right panel, Fig.~\ref{fig:maps}) indicate that the lower patch
was co-spatial with an intergranular vertex (white patch in the $V_{\mathrm{tot}}$ map).
Therefore, the $V$ profiles in the lower patch show a redshift, which is also
the case for the 630.25\,nm OP profile.

Figure 2 shows the OP profile. The positive (negative) magnetic polarity
in this figure corresponds to black (white) in the $V_{\mathrm{tot}}$ map.
Both lines have positive amplitude and area asymmetries.
The 630.15\,nm line of the OP profile has
a blueshifted zero-crossing of $v_{\mathrm zc}\!\!\approx\!\!-1\,$km\,s$^{-1}$, while
the 630.25\,nm line has a redshifted zero-crossing of $v_{\mathrm zc}\!\approx\!+2\,$km\,s$^{-1}$.

The OP profile along with its neighboring profiles are shown in Fig.~\ref{fig:array}.
In this figure, the slit direction is vertical and the OP profile is at the center.
The Stokes-$V$ profiles above and below the OP profile
show strong $V$ signals of normal shape
which shows that the OP profile was located at the center of a polarity reversal.
The polarity of the 630.25\,nm line of the OP profile is the same as that of the lower pixels
and the polarity of the 630.15\,nm line is identical to that of the upper ones.
The $V$ profile left to the OP profile (Fig.~\ref{fig:array}, second row, left) is also strange:
a pathological profile for the 630.15\,nm line and a regular profile of the 630.25\,nm line
with a polarity similar to the 630.25\,nm line in the OP profile.
The same is true for one scan-step after the OP profile (Fig.~\ref{fig:array}, second row, right column):
a normal $V$ profile at 630.15\,nm  with a polarity identical
to the 630.15\,nm line of the OP profile and no signal for the 630.25\,nm line.

Note that the OP profile can not be reproduced by seeing effects as
the adjacent profiles above and below the central OP profile in
Fig.~\ref{fig:array} have incompatible shifts, line widths, and
amplitudes: A least-squares fit to a superposition of the adjacent
upper and lower profiles yields residuals
exceeding 2$\sigma$. Independently, we
measure from the speckle burst that the standard deviation of the
image motion during the
exposure was only about 0.1\,arcsec.
Therefore the OP profile is not due to seeing effects.

\begin{figure}
\begin{center}
\resizebox{8.9cm}{!}{\includegraphics*{7441fig2.eps}}
\caption[]{Stokes-$V$ line pair that shows opposite polarity in the 630.15 and 630.25\,nm lines.
The position of the Stokes-$I$ line-core and the Stokes-$V$ zero-crossing are indicated by
vertical solid and dashed lines, respectively.
The line-core (v$_\mathrm{\begin{scriptsize}I\end{scriptsize}}$)
and zero-crossing (v$_{\mathrm{\begin{scriptsize}zc\end{scriptsize}}}$) velocities are in m\,s$^{-1}$.
}
\label{fig:op}
\end{center}
\end{figure}
\section{The model profile}
The basic properties of the observed OP profile are: (i) the opposite
polarities of the $V$ profiles at 630.15\,nm and 630.25\,nm, (ii) the
velocity of the deeper forming line (630.25\,nm) is positive
(downward) and the velocity of the higher forming line
(630.15\,nm) is negative (upward), and (iii) both
lines show positive area and amplitude asymmetries. The vicinity of
the OP profile is characterized by two patches of opposite polarity
that are separated by only one resolution element as indicated by the
rectangle in the middle panel of Fig.~\ref{fig:maps} and in the
profile array of Fig.~\ref{fig:array}. This configuration suggests
that we witness an electric current sheet in the intermediate
atmospheric layers with magnetic field of one polarity in deeper
layers and of the opposite polarity in higher layers.
The field lines should be inclined with respect to the line-of-sight
since there is a $U$ signal in the
observed Stokes parameters~(Fig.~\ref{fig:fit}). The finding
of upward velocity in the higher layer and downward velocity in the
deeper layer presents evidence for magnetic reconnection, since
such an event produces a bipolar jet.
In the central region in-between the two magnetic
polarities one expects an enhanced temperature due to Joules
heating~(gray scale, Fig~\ref{fig:model}).
This basic picture is sketched in the left panel of Fig.~\ref{fig:model}.
Based on these ideas, we construct in the following a quantitative model
which allows us to compute a set of synthetic Stokes profiles for a
direct comparison with the observed OP profile.

We use the SIR code~\citep{sir92} to synthesize the Stokes profiles of
the \ion{Fe}{i}\,630\,nm pair.  In order to fit the degree of
polarization and to account for straylight, the model atmosphere must
contain a field free (non-magnetic) component beside the magnetic
component. The non-magnetic component is taken to be the HSRA model
atmosphere~\citep{hsra}.  This thermal stratification is also used for
the magnetic component with a slight modification: in an intermediate
layer around $\log\tau_{500\,\mathrm{nm}}\!\approx\!-1.2$, we
introduce a temperature bulge.  At the peak of this hot bulge, it is
270\,K hotter than in the HSRA atmosphere.  The hot bulge is partly
field free and partly overlaps with the magnetic field of the deeper
layer.  Above and below this hot bulge, the velocity and magnetic
field have opposite orientations.  The ranges of non-vanishing
magnetic field strength and velocity overlap, but do not have their
peak values at the same optical depth.  Constant values for the
inclination of the magnetic field with respect to the vertical
(35\degr) and the azimuth (135\degr) are used to fit the ratios
between $Q$, $U$, and $V$.  The model atmosphere that was finally used
to reproduce the OP profile is shown in Fig.~\ref{fig:model}~(right).
All structure elements of this model atmosphere are necessary
ingredients for a successful reproduction of the OP profile.
In particular, the temperature bump is essential to reproduce the different V-profile
polarities when using only one line-of-sight along which two magnetic
polarities are present.
This demonstrates that OP profiles can be synthesized with a
realistic model, but note that this solution may not be unique.

To obtain the filling factor of the magnetic atmosphere,
we perform an inversion of the regular Stokes profiles  in the surrounding pixels of the OP profile.
The inversion set up is similar to that of \cite{luis_beck05}, except that
we assume a constant value for the stray light and allow for a linear gradient of the line-of-sight
velocity in the magnetic component.
From these inversions we obtain magnetic filling factors of 10\,--\,30\,\%, so
we assume a filling factor of 20\,\% for our model atmosphere.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics*{7441fig3.ps}}
\caption[]{The OP Stokes-$V$ profile (center) and its surrounding profiles.
Abscissae show the wavelength ($\lambda - 630$\,nm). The slit direction is vertical.
Other parameters are like in Fig.~\ref{fig:op}.}
\label{fig:array}
\end{figure}

The upper layer of magnetic field along with the bump of negative velocity leads
to a very asymmetric $V$ profile for both lines. This profile has the polarity of the
630.15\,nm line of the OP profile and shows a higher amplitude in 630.15\,nm than in 630.25\,nm.
The lower-layer magnetic field of opposite polarity produces an almost
antisymmetric $V$ profile at  630.25\,nm, but the profile at 630.15\,nm is asymmetric.
The combination of these two components gives the final fit to the data
that is shown in Fig.~\ref{fig:fit}.
The quality of the fit is satisfactory considering that we only used one single line-of-sight
for this complex topology. With a slight shift of the velocity peaks we are also able to reproduce
the pathological profile to the left to the OP profile (Fig.~\ref{fig:array}).
\begin{figure}
\centerline{\resizebox{8cm}{!}{\includegraphics*{7441fig4.ps}}}
\caption[]{Left: schematic configuration of the magnetic field, the line-of-sight velocity,
and the temperature enhancement~(gray scale).
Right: stratification of the atmosphere as a function of optical depth for the
magnetic\,(solid) and non-magnetic\,(dashed) components.
From top to bottom: temperature, magnetic field strength
(with the the sign denoting the polarity), and
line-of-sight velocity. The non-magnetic atmosphere has a convective blueshift of
300\,m\,s$^{-1}$.}
\label{fig:model}
\end{figure}

\section{Discussion}

The magnetic area with polarity rendered in dark in the $V_{\mathrm{tot}}$ map
of Fig.~\ref{fig:maps} was present throughout the observing time
of 64 minutes and was persistently
bright in \ion{Ca}{ii}\,H, which suggests that it belonged to network
magnetic fields.
On the contrary, the opposite polarity, white patch, visible in the
center of the $V_{\mathrm{tot}}$ map, was a transient feature -- it assembled from
diffuse flux of white polarity, concentrated and intensified to the white patch seen in
Fig.~\ref{fig:maps}, before it rapidly weakened. Within the scanning cadence of
97~s it virtually disappeared, not without influencing the opposite polarity (black)
neighboring patches to the upper and lower sides, which weakened during this time.
Such events must frequently take place considering that the magnetic flux in the quiet
network is permanently replaced by and in interaction with mixed polarity flux
\citep{schrij_title_etal_97}. When flux concentrations
of opposite polarity collide, flux cancellation and magnetic field reconnection are likely
to be involved \citep{zwaan87}. At the reconnection site we expect conversion from
magnetic to thermal energy and the formation of a bipolar jet: processes that are
compatible with the temperature bulge in between the two layers of oppositely directed
velocities and magnetic fields of the atmosphere of Fig.~\ref{fig:model}  and with the observed flux
cancellation and weakening.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics*{7441fig5.ps}}
\caption[]{Comparison between the observed Stokes profiles (pluses) and the
synthesized ones (solid lines). The horizontal dotted lines show the 3$\sigma$ noise level.}
\label{fig:fit}
\end{figure}

For a check of compatibility with the thermal, kinetic, and magnetic
energy fluxes required in a reconnection event, we first consider a
control volume given by the size of one POLIS pixel of
$A_p = \Delta_x\times\Delta_y = 350\,\mbox{km}\times 210\,\mbox{km}$
and the height range spanned by the FWHM of the temperature bulge
of $h = 135$\,km, which encompasses a mass of $m = 1.3\times 10^{12}$\,kg.
With a bulge peak of $\delta T = 270$\,K we obtain from Newton's
law of cooling a radiative heat loss from this volume of
\begin{equation}
\dot{Q} = 16\, \sigma\, \kappa\, T^3\, \delta T\, m\, f = 9\times 10^{16}\,{\mathrm W,}
\end{equation}
where $\sigma$ is the Stefan-Boltzmann constant,
$\kappa\!=\!1.1\times 10^{-2}\,\mbox{m}^2\,\mbox{kg}^{-1}$ the
opacity corresponding to the mean density of
$1.3\!\times\!10^{-4}\,\mbox{kg}\,\mbox{m}^{-3}$,
$T\!=\!5100$\,K the corresponding temperature of the unperturbed
background atmosphere, and $f=0.2$ the filling factor.
This radiative energy flux should be of the same order of magnitude
as the kinetic energy flux,
\begin{equation}
F_\mathrm{kin} = (1/2)\,\rho\, v^3\, A_p\, f,
\end{equation}
carried by the bipolar reconnection jets. Using the values of the
lower, heavier (downward moving) layer with
$v_\mathrm{LOS} = 3\,\mbox{km\,s}^{-1}$ and
$\rho = 3.0\times 10^{-4}\,\mbox{kg\,m}^{-3}$, and assuming
the velocity to be directed parallel to the magnetic field we obtain
$F_\mathrm{kin} = 1\times 10^{17}$\,W. The kinetic energy of the upwards moving
jet is negligible, because of the lower density.
The kinetic and thermal energy fluxes need to be sustained by
a corresponding influx of magnetic energy of
\begin{equation}
F_\mathrm{mag} = v\, (B^2/2\mu_0)\, \Delta_x\, h\, \sqrt{f}.
\end{equation}
With $B = 0.1\,\mbox{T}$ we obtain
$F_\mathrm{mag} = 1\times 10^{17}\,\mbox{W}$ if the
influx velocity is $1.2\,\mbox{km\,s}^{-1}$, which is a reasonable
value for the horizontal velocity of photospheric magnetic flux concentrations.

Photospheric magnetic reconnection has previously been studied by
\citet{litvinenko99} and \cite{takeuchi_shibata_01}. Since we observe a
flux intensification prior to flux canceling for the white polarity
patch of Fig.~\ref{fig:maps}, the scenario of \citet{takeuchi_shibata_01} consisting of
reconnection induced by convective intensification may apply to this event.

We would like to stress that the model presented here is not necessarily unique and that
other model atmospheres may also reproduce the data.
However, we expect any alternative model to show strong gradients in velocity and
magnetic field strength.

\section{Conclusion}
We have discovered several sets of opposite polarity Stokes-$V$ profiles.
Each set is a single spectrum of the
630.15 and 630.25\,nm neutral iron lines pertaining to a single resolution element.
We use a model atmosphere consisting of a non-magnetic and one magnetic component
for synthesizing the observed profiles.
The magnetic component contains strong gradients along the line-of-sight in the magnetic field
strength and the velocity. In particular, the higher and deeper layers have opposite
magnetic polarity.
In accordance with the observed red- and blue-shifts of the spectral lines,
the proposed model atmosphere contains a bipolar flow of material
along the line-of-sight and it contains a temperature enhancement in-between.
These atmospheric elements provide evidence of magnetic reconnection in the solar photosphere.
Our estimated values for the thermal, kinetic, and magnetic energy fluxes are also
compatible with a reconnection event.\vspace{-.1cm}

\begin{acknowledgements}
We wish to thank C. Beck and F. W\"oger for their
assistance during observation and data reduction.
We are grateful to J. Bruls, O.~von~der~L\"uhe, and E.~Khomenko for their useful comments.
The POLIS instrument has been a joint development of the High
Altitude Observatory (Boulder, USA) and the Kiepenheuer-Institut.
Part of this work was supported by the Deutsche Forschungsgemeinschaft (SCHM 1168/8-1).\vspace{-.2cm}
\end{acknowledgements}