First direct metallicity measurement of a lensed star-forming galaxy at z=1.7

Abstract: We present the rest-frame optical spectrum of a strongly lensed galaxy at redshift z =1.7 behind the cluster Abell 1689. We detect the temperature sensitive auroral line [O III] 4363, which allows the first direct metallicity measurement for galaxies at z > 1. Our high signal-to-noise spectrum indicates that the target is an extremely low metallicity star-forming galaxy.We estimate an intrinsic absolute B band magnitude of M_{B}=-18.3 \pm 0.1$, with a stellar mass of 4.4$\pm1.2\times10^{8}$ M$_{\odot}$. This galaxy extends the luminosity-metallicity relation of star-forming galaxies at z > 2 by more than an order of magnitude. Given the double-nuclei like morphology and velocity profile of \ha, we tentatively suggest that it could be a merger or a proto-rotating disk galaxy.

\section{Introduction}
Five years ago,  our knowledge about galaxy evolution still had a glaring gap between redshifts 1.5 and 3, termed the redshift desert. The desert is even more desolate for galaxy chemical evolution studies, where high signal-to-noise (S/N) spectra are needed for reliable measurements.  However, great
progress has been made in the last few years in studying galaxies at z$\sim$ 2
and we are beginning to harvest the ``desert".  To date, metallicities are available for
a few tens of  individual galaxies  at $1.5<z<3$ \citep{Shapley04, Erb03, Erb06, Hayashi09}.
These galaxies are selected using broad-band colors either in the UV  \citep[Lyman Break technique;][]{Steidel96,Steidel03}
 or using $B$, $z$, and $K$-band colors \citep[BzK selection;][]{Daddi04}.
 The Lyman break and $BzK$ selection techniques favor the galaxies that are luminous in the UV or  blue and could therefore be biased against low  luminosity (low-metallicity) galaxies, and dusty (potentially metal-rich) galaxies.

To date, metallicity measurements for high-redshift ($z \geq 1$) galaxies have been made
via indirect methods from the ratios of strong nebular lines.  In local galaxies, the preferred method for determining metallicity
is via the ``direct" $T_{e}$ method which relies on the electron temperature sensitivity of the {[O\,{\sc iii}]}\lam4363 line.
At high redshift,  metallicity measurement via the direct method has not been possible because the
{[O\,{\sc iii}]}\lam4363 line used to derive the electron temperature is very weak.
The highest redshift that the {[O\,{\sc iii}]}\lam4363 line has been reported is at $z \sim 0.8$ by \citet{Hoyos05} and \citet{Kakazu07}.

Strong lensing by galaxy clusters offers an alternative tool to explore the metallicity of galaxies at high-redshift.
One advantage of lensing selected galaxies is that they are immune to color and luminosity biases.
Previously, only five lensed galaxies had been observed with sufficient spectral coverage and S/N to yield metallicities
 \citep{Yee96,Pettini00,Pettini02,LB03,Nesvadba07,Stark08}. Compared with Lyman break selected
galaxies in the same redshift and luminosity range, the five lensed objects
span a substantially broader metallicity range (12\,+\,log(O/H) $\sim$7.1---9.1), indicating that current metallicity history
studies at $z > 1.5$  (12\,+\,log(O/H) $\sim$8.4---8.8) may not be sampling the full metallicity range at this redshift.
It is worth noting that the metallicity range quoted for the lensed galaxies is extremely uncertain, depending
on various assumptions about reddening and metallicity calibrators.
Another great advantage of gravitational lensing is its natural magnification of the source.
The flux of the lensed background galaxy is commonly
 boosted by a factor of 1-3 mag or more around the critical lines of the cluster center.
Optical spectroscopic surveys of high-z galaxies near cluster critical lines have already yielded a large catalog of lensed galaxies beyond z $\sim$1
 \citep[e.g.,][]{Broadhurst05, Richard08, Frye07, Limousin07}. Unfortunately,
spectra from these surveys cannot be directly utilized in metallicity studies because the commonly used metallicity diagnostic emission lines have shifted into
the near infrared (near-IR) at  $z>1$.

With the revolutionary multi-object cryogenic near-IR spectrographs that have been recently installed on 8-10m class telescopes,
high efficiency near-IR spectroscopy  is now available.   We carried out
the first near-IR spectroscopic survey of lensed galaxies behind lensing clusters, aimed at measuring the metallicity of galaxies
at $1.5<z<3$.  This Letter reports the rest-frame optical spectrum
of a strongly lensed star-forming galaxy at $z = 1.7$ behind the Abell cluster 1689 from our ongoing survey.
The galaxy cluster Abell 1689 was chosen as our first target because it has the largest Einstein radius
of $\sim 50\,\!\!^{\prime\prime}$ \citep{Broadhurst05},
and nearly $\sim$ 100 spectroscopically confirmed lensed images \citep{Limousin07}.
The galaxy in question is designated as ID 22.3 in the catalog of \citet{Broadhurst05},
we refer to it as Lens22.3 hereafter.  Throughout we use the standard lambda cold dark matter ($\Lambda$CDM) cosmology with $H_0=70$ km s$^{-1}$ Mpc$^{-1}$,
$\Omega_M=0.3$ and $\Omega_{\Lambda}=0.7$.

\section{Observations and Data Reduction}
\begin{figure}[!ht]
\epsscale{1.1}
\plotone{f1.eps}
\caption{Observed {\it Subaru} MOIRCS spectrum of Lens22.3 at z= 1.7.
The extracted one-dimensinoal spectrum is shown in black, with the identified emission lines labeled.  The sky telluric absorption and OH emission
lines are shown in red and blue, respectively.  The two-dimensional spectrum for each of the detected lines is shown on top.}
\label{fig:spec}
\end{figure}

The data were obtained on 2008 March 26th,  using the
Multi-Object InfraRed Camera and Spectrograph (MOIRCS) at {\it Subaru} \citep{Ichikawa06},
in the spectroscopic mode \citep{Tokoku06}.    We used the  low-resolution grism zJ 500 with a nominal
resolution of 500 and  a wavelength coverage of 0.9---1.78 $\mu$m.
We designed two masks targeting  the spectroscopically confirmed lensed images in the field of Abell 1689.
We  obtained useable data for one mask.
During the observation for this mask,  the seeing varied between 0.$\!\!^{\prime\prime}$3 and 0.$\!\!^{\prime\prime}$6,  and the airmass
 ranged between 1.1 and 1.7. Lens22.3 was placed at one of the center slits of the mask, with a slit width of $\sim 1\, \!\!^{\prime\prime}$ and a slit length of
 $\sim 9\, \!\!^{\prime\prime}$.
The telescope was nodded along the slit in a ABAB mode with a dithering length of 2.$\!\!^{\prime\prime}$5.
The exposure time for each individual frame was 400-600 seconds, resulting in a total integration time of 6800 seconds.

We reduced the data using the standard {\sc iraf} and IDL routines.
Sky OH lines were removed by subtracting the closest dithering pair A and B.
A scaling factor between each pair was determined to optimize the sky subtraction.
The sky subtracted A-B pair was flat-fielded, cleaned for comic rays, and then traced and rectified using a dome flat spectrum.  The final two-dimensional spectrum was created by co-adding the reduced individual pairs.
 Wavelength calibration was carried out using the strong OH lines between 0.9 and 1.8 $\mu$m.
A F5V standard star was observed at approximately the same airmass as the target. The Telluric absorption
was obtained by dividing the star spectrum with a fitted blackbody curve. The resulting stellar spectrum
was then divided into the object spectrum to correct for telluric absorption.

\section{Analysis and Results}
\subsection{Metallicity}

\begin{deluxetable}{lccccc}
\tabletypesize{\footnotesize}
\tablewidth{0pt}
\tablenum{1}
\tablecolumns{6}
\tablecaption{Observed Line Intensities and Ratios\label{tab:obs}}
\tablehead{
\colhead{Ion Line} &
\colhead{$z_{obs}$}&
\colhead{Flux} &
\colhead{Flux (corrected)} &
\colhead{1$\sigma$} &
\colhead{S/N} \\
\cline{1-6}\\
\colhead{(1)} &
\colhead{(2)} &
\colhead{(3)} &
\colhead{(4)} &
\colhead{(5)} &
\colhead{(6)} \\
}
\startdata
{[N\,{\sc ii}]}\lam6584& \nodata&   $<$0.009 &$<$0.025 &\nodata & \nodata \\
{H$\alpha$}\lam6563 & 1.705  &3.12 &5.03&$\pm$0.4 & $13$  \\
{[O\,{\sc iii}]}\lam5007& 1.706 & 6.45 &6.45& $\pm$ 0.3 & $15$  \\
{[O\,{\sc iii}]}\lam4959 &1.705 &  1.98 &1.98& $\pm$0.3 & $10$  \\
{H$\beta$}\lam4861 &1.705   &1.0 &1.0& $\pm$0.1 & $5$  \\
{[O\,{\sc iii}]}\lam4363 & 1.696&0.27 &0.27&  $\pm$0.1 & $3$  \\
{[O\,{\sc ii}]}\lam3727 &1.708   &0.95 &1.11& $\pm$0.3 & $8$  \\
\enddata
\tablecomments{
(1) Detected emission lines and their rest frame wavelength.
(2) Redshift calculated from the observed wavelength.
(3) Measured flux, relative to {H$\beta$}, detector response uncorrected.
(4) Measured flux, relative to {H$\beta$},  detector response corrected.
(5) 1 $\sigma$ error of the line flux from Gaussian fitting.
(6) Signal-to-noise of the line flux.}
\end{deluxetable}

At redshift 1.7, all the rest-frame optical metallicity diagnostic lines fall into the spectral coverage
 of MOIRCS zJ500 grism. We have detected all the rest frame nebular emission lines from {[O\,{\sc ii}]}\ to {H$\alpha$},
 including the usually very weak {[O\,{\sc iii}]}\lam4363 line. Figure~\ref{fig:spec} shows our final two-dimensional and one-dimensional rest-frame optical spectrum of Lens22.3.   Table~\ref{tab:obs} lists the measured line fluxes relative to {H$\beta$},  the 1 $\sigma$ error,
and S/N.  Since the metallicity studies use only line ratios, we did not  attempt to flux calibrate the spectrum,
which is nontrivial in near-IR observations.
 The instrument response is essentially constant for the wavelength range from
{[O\,{\sc iii}]}\lam4363  to  {[O\,{\sc iii}]}\lam5007.  The {H$\alpha$}\  line lies in the region where the
response has dropped by $\sim 35\%$. The response on the expected {[N\,{\sc ii}]}\ line  position has dropped by
$\sim 65\%$, which may be the reason that there is no detection of {[N\,{\sc ii}]}\  within 3$\sigma$.
After correcting for instrumental response,  we obtain an upper limit for {[N\,{\sc ii}]}/{H$\alpha$}.
Metallicities were calculated using the following two methods.

{\it 1. Strong line diagnostics.}  In the absence of direct measurements from electron temperature,
 the ratio of strong nebular lines is used to determine the oxygen abundance. The
well-known discrepancies between the different line ratio indicators can be removed by converting
 to the same base calibration following the methods of \citet{Kewley08}.
We use the indicator $R_{23}$\,=\,({[O\,{\sc ii}]}\lam3727 + {[O\,{\sc iii}]}\llam4959, 5007)/{H$\beta$}\/
 to calculate metallicity.  In our case, the upper and lower branch degeneracy of $R_{23}$
can be broken by the upper limit of {[N\,{\sc ii}]}/{H$\alpha$}\,$\,<\,$0.005.
In addition, the detection of the weak {[O\,{\sc iii}]}\lam4363 line almost certainly means the adoption of the lower branch.  For the actual calculation,  we use
 the lower branch calibration of \citet{McGaugh91} (M91 method) in the analytical form given by \citet{Kobulnicky99}.
 The resulting metallicity is  \mbox{12\,+\,log(O/H)\,=\,}\ 8.0 in the M91 method.   Converted to the \citet{KD02} method (KD02),
 we have \mbox{12\,+\,log(O/H)\,=\,}\ 8.1 assuming zero extinction, or 8.3 assuming a high extinction value of E(B-V)\,=\,0.4 \citep{Cardelli89,Osterbrock89}.   The high extinction value agrees with the
 direct estimate from the {H$\alpha$}/{H$\beta$}\ in Table~\ref{tab:obs}. However,  because  the {H$\alpha$}\  line lies near the edge where the response has dropped by $\sim 35\%$,  the ratio of {H$\alpha$}/{H$\beta$}\ can be largely uncertain because it depends sensitively on the precision of correction for the detector response.  Therefore instead of using only the high extinction from {H$\alpha$}/{H$\beta$}, we use two extreme values of extinction to quantify the range of uncertainty caused by extinction.
 Note that the logarithm in the definition of 12\,+\,log(O/H)\ means that 10\% of change in the line ratio  causes only  $\sim$ 0.04 dex change in metallicity.  The uncertainty for our derived metallicity is therefore $<$ 0.1 dex prior to
 extinction correction.

{\it 2. Direct metallicity from electron temperature $T_{e}$.}
The auroral line {[O\,{\sc iii}]}\lam4363 is often very weak, even in  low metallicity environments.  It is usually not observed in
high metallicity galaxies without very sensitive, high S/N spectra \citep[e.g.,][]{Garnett04}.
The 3 $\sigma$ detection of  {[O\,{\sc iii}]}\lam4363 in Lens22.3  strongly suggests a low metallicity environment.
The ratio of {[O\,{\sc iii}]}\lam5007, 4959 and  {[O\,{\sc iii}]}\lam4363
allows a direct measurement of the oxygen abundance via electron temperature $T_{e}$.
We use the {\tt nebular} package in {\sc iraf} to derive the electron temperature
$T_e$\,=\,2.3557 $\times$ 10$^{4}$ K, assuming an electron density n$_{e}$ of 100 cm$^{-3}$ ($T_e$ is
relatively insensitive to $n_{e}$, e.g.,  for $n_{e}$\,=\,1000 cm$^{-3}$, $T_e$\,=\,2.3508 $\times$ 10$^{4}$ K).  Following the procedure of \citet{Izotov06}, we obtain \mbox{12\,+\,log(O/H)\,=\,} 7.5 $\pm$ 0.1 assuming zero extinction; or \mbox{12\,+\,log(O/H)\,=\,} 7.3 $\pm$ 0.1 assuming
a high extinction value of E(B-V)\,=\,0.4.

The 0.6 dex discrepancy between the KD02 and  $T_{e}$ method is not surprising since the metallicity measured from the direct $T_e$ method is systematically lower than other methods;  such offsets can be as large as 0.7 dex \citep{Liang06, Yin07b,Nagao06}.  Regardless of the offset between the two methods, both metallicity estimates indicate that
 Lens22.3 is  low in oxygen abundance.  The evidence that Lens22.3 is a star-forming galaxy instead of a narrow line AGN comes from its position on the traditional
 Baldwin---Phillips---Terlevich (BPT) diagram \citep{Baldwin81}.  Using the upper limit of  {[N\,{\sc ii}]}/{H$\alpha$}\  and the ratio of {[O\,{\sc iii}]}/{H$\beta$},  Lens22.3 is  located on the leftmost region of the star-forming branch of the BPT  {[N\,{\sc ii}]}\ diagram, excluding its possibility of being an AGN \citep{Kewley06}.

\subsection{Photometry}
We searched the {\em Hubble Space Telescope} ({\it{HST}}) archive and obtained the photometry for Lens22.3.
The observed AB magnitudes are 23.2 $\pm$ 0.1, 23.6 $\pm$ 0.2, 23.3 $\pm$ 0.1, 23.4 $\pm$ 0.1 and 23.1 $\pm$ 0.1 in the {\textit{ACS}} {\ensuremath{g_{F475W}}}, {\ensuremath{r_{F625W}}}, {\ensuremath{i_{F775W}}}, {\ensuremath{z_{F850LP}}}, and {\textit{NICMOS}}  {\ensuremath{J_{F110W}}} bands, respectively.
After correcting for the lensing magnification of 15.5 $\pm$ 0.3 (provided by J. Richard 2009, private communication),
the absolute  $B$ band magnitude is $M_{B}$\,=\,-18.3 $\pm$ 0.1 (Vega System).
A simple spectral energy distribution (SED) fitting using the \verb+kcorrect v4_1_4+ code of \citet{Blanton07} yields a stellar mass of
1.6 $\times$ 10$^{9}$ M$_{\odot}$.  Since the \verb+kcorrect+ may overestimate the stellar masses of high-redshift galaxies \citep{Blanton07},  we recalculated the stellar mass using the
advanced version of the code
\verb+LE PHARE+ \citep[photometric redshift and simulation package;][]{Ilbert09} based on
population synthesis models of \citet{BC03}.
The resulting stellar mass is 4.4$\pm$1.2$\times$10$^{8}$ M$_{\odot}$.

\subsection{Morphology and Kinematics}
\begin{figure}[!ht]
\epsscale{1}
\plotone{f2.eps}
\caption{{\it HST} rescaled F775W image for the field of Lens22.3  The image has been rescaled and registered onto the MOIRCS chip. The slit on Lens22.3 has a width of 1\,$\!\!^{\prime\prime}$ and a length of 9\,$\!\!^{\prime\prime}$.  Slit $X$ (dispersion) and $Y$ (spatial) direction are marked.  A zoomed view
of Lens22.3 with intensity contours (0.3, 0.2, 0.1, 0.04 of the peak intensity) is shown on the lower left corner.}
\label{fig:lens}
\end{figure}
\begin{figure}[!ht]
\epsscale{1}
\plotone{f3.eps}
\caption{Top two panels: velocity profile of {H$\alpha$}\ and {[O\,{\sc iii}]}\ from spectra extracted within the spatial apertures
 A, B, C, D (bottom right panel).  Bottom left panel:  velocity profile for the sky line nearest to {H$\alpha$}\ using the
 same apertures. Bottom right panel:   two-dimensional spectrum of {H$\alpha$}, spatially cut into sub-apertures A, B, C, D.
The velocity zero point is chosen to be the averaged center of the cuts.}
\label{fig:lineprof}
\end{figure}
Thanks to the excellent seeing ($\sim$  0.$\!\!^{\prime\prime}$4)
 and the powerful magnification of gravitational lensing,  the target is spatially resolved.
 In Figure~\ref{fig:lens}, we show the {\textit{HST}} F775W image in the vicinity of Lens22.3.
The image has been registered onto the MOIRCS chip, which has a pixel scale of 0.$\!\!^{\prime\prime}$117 pixel$^{-1}$. The MOIRCS slit layout is shown, where
$X$ is the dispersion direction and $Y$ is the spatial direction which is also along the elongated side of the galaxy.
On the lower right corner is a zoomed image of Lens22.3, with various intensity contours overlaid.
It has a clear double-nuclei-like morphology.  The  nuclei are separated by $\sim$ 0.$\!\!^{\prime\prime}$4,
corresponding to a physical scale of $\sim$3 kpc  at the redshift of $z = 1.7$.
The intrinsic separation will be smaller than this value because the image has been elongated by lensing effects.
 We divide the two-dimensional spectrum into subsections  A, B, C, D along the spatial direction $Y$,  using a 3 pixel bin that roughly matches the average seeing of the observation.
We extract the spectrum for apertures A, B, C, D using the strongest line {[O\,{\sc iii}]}\lam5007 and {H$\alpha$}.  The velocity profile relative to the center of the subsections are given in Figure~\ref{fig:lineprof}.
For comparison,  the velocity profile for the sky OH line nearest to {H$\alpha$}\  extracted from the same apertures is also shown. There seems to be a systematic rise in the  velocity from A to D in {H$\alpha$},  but not in {[O\,{\sc iii}]}\lam5007.
The misalignment of the line profile in the spatial direction and the double-nuclei like morphology may indicate that
the system is rotating or merging.   Future follow up with integral field spectroscopy with adaptive optics is highly desirable to reveal the kinematics
of the system.

\begin{figure*}[!ht]
\epsscale{0.9}
\plotone{f4.eps}
\caption{ Left: mass---metallicity relation.
Dash line:  local relation from the SDSS sample \citep{Tremonti04}.
Filled black circle:  relation for galaxies z $\geq$ 2 \citep{Erb06}.
Filled red square: Lens22.3, without extinction correction.
Empty red square: Lens22.3, after extinction correction.
Right: luminosity---metallicity relation.
Dash line:  local relation from SDSS sample \citep{Tremonti04}.
Filled black circle:  relation for galaxies z $\geq$ 2 \citep{Erb06}.
Filled red square: Lens22.3, without extinction correction.
Empty red square: Lens22.3, after extinction correction.
All metallicities have been recalculated to the same base of KD02 \citep{Kewley08}.
}
\label{fig:ml}
\end{figure*}

\subsection{Mass-Metallicity Relation}
In Figure~\ref{fig:ml},  we show the mass/luminosity---metallicity relation for the local \citep{Tremonti04} and high-redshift
\citep{Erb06}  star-forming galaxies.  All metallicities in Figure~\ref{fig:ml} have been rescaled to the same KD02 method
using the calibration of \citet{Kewley08}.
The local sample of \citet{Tremonti04} is based on $\sim$ 53,000 Sloan Digital Sky Survey (SDSS) star-forming galaxies and
the \citet{Erb06} sample is based on composite spectra of 87 rest-frame UV-selected star-forming galaxies at z $\geq$ 2.   Lens22.3 is consistent with the mass/luminosity---metallicity relation of the high z sample.  We evaluate the effect of extinction
 assuming a high extinction value of E(B-V)\,=\,0.4.
   Extinction causes  metallicity to rise by $\sim$ 0.1 dex,  but
 it does not significantly change the position of Lens22.3 on the mass---metallicity diagram.
Lens22.3 is currently the lowest mass system with direct measurement on the mass---metallicity relation diagram;  Lens22.3 extends the luminosity---metallicity relation to the fainter end by more than an order of magnitude, demonstrating the capability of using gravitationally lensed samples
to probe intrinsically fainter and less massive systems than are found in existing samples.

\section{Summary}
We reported the discovery of a low metallicity star-forming galaxy at $z = 1.7$.
We detect the temperature sensitive auroral line {[O\,{\sc iii}]}\lam4363
which allows us to derive  metallicity using the direct $T_{e}$ method for the first time at high redshift.
The low metallicity and low luminosity/mass of the $z = 1.7$ star-forming
galaxy place it at the lowest end of the mass/luminosity---metallicity relation, extending the
$z>2$ luminosity-metallicity relation by more than an order of magnitude.
Our on-going near-IR spectroscopic survey of lensing clusters aims to obtain
metallicities for more objects like Lens22.3
 in the redshift range 1.5 $<$ z $<$ 3.

\acknowledgments
We want to thank the referee for her/his valuable comments and suggestions that significantly
improved this Letter.  We thank Johan Richard for providing us with the lensing magnification of Lens22.3, and Youichi Ohyama for his help with MOIRCS data reduction.  We also thank Katelyn Allers for useful discussion on Telluric correction methods and C.-J. Ma for running the \verb+LE PHARE+  code for us.   We recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community.  We appreciate the opportunity to conduct observations from this sacred mountain.

{\it Facilities:} {\it Subaru} (MOIRCS)

Cosmology and Nongalactic Astrophysics (astro-ph.CO)


Charm spectroscopy and exotics at ZEUS

Abstract: Light and charmed hadrons are produced copiously in ep collisions with a centre-of-mass energy of 318 GeV at HERA. Results of the ZEUS Collaboration on pentaquark searches, deuteron and antideuteron production and charmed-meson spectroscopy, obtained using the HERA I data, are summarised.

\section{Introduction}
Light and charmed mesons and baryons are produced copiously in $ep$
collisions with a centre-of-mass energy of $318\,$GeV at HERA.
During the first phase of the HERA operation (1992-2000),
the ZEUS Collaboration accumulated a data sample corresponding
to an integrated luminosity of $\sim120\,$pb$^{-1}$.
Results
on pentaquark searches,
deuteron and antideuteron production
and charmed-meson spectroscopy,
obtained by the ZEUS Collaboration
using the HERA I data, are summarised
in this note.

\section{Strange pentaquarks}
A peak in the $K^0_s p ({\bar{p}})$ invariant mass spectrum around $1520\,$MeV
was observed in deep inelastic scattering
(DIS) by the ZEUS Collaboration~\cite{zeus_theta}.
In Fig.~1, the spectrum is shown for
$Q^2>20\,$GeV$^2$,
where $Q^2$ is
the exchanged-photon virtuality
The statistical significance of the signal varies between 3.9$\,\sigma$
and 4.6$\,\sigma$ depending upon the treatment of the background.
The candidate $\Theta^+$ signal was found to be produced predominantly
in the forward (proton) direction in the laboratory
frame~\cite{zeus_theta_study}.
This is unlike the case for the $\Lambda(1520)$ and can indicate that
$\Theta^+$ may have an unusual production
mechanism related to the proton-remnant fragmentation.
Using a smaller integrated luminosity,
the H1 Collaboration did not observe
the candidate $\Theta^+$ signal
and set an upper limit
on its production~\cite{h1_theta}.
The production cross section for the $\Theta^+$ candidate, measured by the
ZEUS Collaboration~\cite{zeus_thexs}, is
$125\pm27({\rm stat.})^{+37}_{-28}({\rm syst.})\,{\rm pb}$
for $Q^2>20\,$GeV$^2$. The value is larger than, but consistent with, the upper
limit obtained by the H1 Collaboration.

\begin{figure}[b]
\centerline{\psfig{file=gladilin_fig1.eps,width=2.6in}}
\caption{
Invariant mass spectrum for the $K^0_s p ({\bar{p}})$ combinations
obtained by the ZEUS Collaboration.}
\label{fig1}
\end{figure}

The ZEUS Collaboration performed also a search
for two strange pentaquarks, reported by
the NA49 Collaboration~\cite{na49_pq},
and observed no signal in the $\Xi \pi$ invariant mass spectrum~\cite{zeus_xi}.
Upper limits on the ratio of a possible $\Xi_{3/2}^{--}(\Xi_{3/2}^0)$
signal to the $\Xi^0(1530)$ signal were set in the mass range
$1.65-2.35\,$GeV.

\section{Deuterons and antideuterons}

A first observation of deuteron and antideuteron production
in $ep$ collisions in the DIS regime has been recently reported by the ZEUS
collaboration~\cite{zeus_deuteron}.
To identify the particles,
the measurement of
the energy losses, $dE/dx$, in the central tracking detector
has been used.
\begin{figure}[b]
\centerline{\psfig{file=gladilin_fig2.eps,width=2.7in}}
\caption{
The mass spectra for positive and negative tracks
with $dE/dx>2.5\,$mips.
The arrows show cuts applied for selection of (anti)proton and
(anti)deuteron candidates.}
\label{fig2}
\end{figure}

Figure~2 shows the mass spectra for positive and negative tracks
with $dE/dx>2.5\,$mips (minimum ionising particles)
in DIS events with $Q^2>1\,$GeV$^2$.
For each track, the mass has been calculated from the track momentum
and $dE/dx$ value
using a parametrisation based on the Bethe-Bloch equation.
The number of deuteron and antideuteron candidates in the mass window
$1.5<M<2.5\,$GeV is 309 and 62, respectively. No antitritons
have been observed, Given the small number of the triton candidates,
no statement on the triton observation in DIS has been made.

Small contributions to the observed signals from the secondary interactions
in the beam-pipe and detector material have been subtracted using side bands of
the distribution on the distance of the closest approach of a track to
the beam spot in the transverse plane.
Possible contamination of the observed deuteron signal
due to interactions of the proton or electron beam with the residual
gas in the beam-pipe has been found to be below $20\%$.
\begin{figure}[b]
\centerline{\psfig{file=gladilin_fig3.eps,width=2.7in}}
\caption{
The corrected (top) ${\bar d}/{\bar p}$ and
(bottom ) ${\bar p}/p$ production ratios as a function
of $p_T/M$.
Open circles show the ${\bar d}/{\bar p}$ ratios obtained
by the H1 Collaboration in the photoproduction regime.}
\label{fig3}
\end{figure}

Figure~3 shows the ${\bar d}/{\bar p}$ and
${\bar p}/p$ production ratios, corrected for the detector and trigger effects,
as a function of $p_T/M$.
The measured ${\bar p}/p$ ratio is consistent with unity.
The ${\bar d}/{\bar p}$ ratios are in fair agreement with
those obtained in $pp$ interactions~\cite{pp_deuteron},
in hadronic $\Upsilon$ decays~\cite{argus_deuteron},
and by the H1 Collaboration at HERA in
the photoproduction regime ($Q^2<1\,$GeV$^2$)~\cite{h1_deuteron}.

No antideuteron candidates have been found in the current region of
the Breit frame~\cite{breit}. Since the region is similar to a single
hemisphere of the $e^+e^-$ annihilation process, the antideuteron observation
in DIS does not contradict to the small antideuteron rates observed
at LEP~\cite{lep_deuteron}.

\section{Charmed pentaquark}

An observation of a candidate
for the charmed pentaquark state, $\Theta^0_c = uudd{\bar c}$,
decaying to $D^{*\pm}p^\mp$
was reported by the H1 Collaboration~\cite{h1_ch5q}.
A fit of the signal in DIS yielded $50.6\pm11.2$ signal events and
the mass of $3099\pm3({\rm stat.})\pm5({\rm syst.})\,$MeV.
The observed resonance was reported to contribute
$(1.46\pm0.32)\%$ of the $D^{*\pm}$ production rate
in the kinematic range studied in DIS~\cite{h1_thetac_study}.
The measured differential distributions were found to be
consistent with the fragmentation modelling of the candidate $\Theta^0_c$
production.

\begin{figure}[b]
\centerline{\psfig{file=gladilin_fig4.eps,width=2.7in}}
\caption{
The distributions of $M(D^{*\pm}p^\mp)$
obtained by the ZEUS Collaborations.}
\label{fig4}
\end{figure}

The observation of the H1 Collaboration was challenged
by the ZEUS Collaboration~\cite{zeus_ch5q}.
Using a larger sample of $D^{*\pm}$ mesons,
ZEUS observed no signature of the narrow resonance in
the $M(D^{*\pm}p^\mp)$ spectra shown in Fig.~4.
The Monte Carlo $\Theta^0_c$ signals normalised to $1\%$
of the number of reconstructed $D^{*\pm}$ mesons are shown
on top of the fitted backgrounds.
The upper limit on the fraction of $D^{*\pm}$ mesons
originating from $\Theta^0_c$ decays was evaluated to be
$0.23\%$ ($95\%$ C.L.). The upper limit for DIS with $Q^2>1\,$GeV$^2$
is $0.35\%$ ($95\%$ C.L.).

\section{Studies of excited $D$ mesons}

Sizeable
production of the excited charmed, $D^0_1$ and $D_2^{*0}$,
and charmed-strange, $D_{s1}^+$,
mesons has
been observed in $ep$ interactions
by the ZEUS collaboration~\cite{zeus_d1d2,zeus_d1s}.
Figure~5 shows the distribution
of the $M(\dspm K^0_s)$
for the $\dspm$ candidates.
A clear signal is seen
at the nominal value of $M(\dsone )$.
The measured $\done$, $\dtwo$ and $\dsone$ rates were converted to
the fractions of $c$ quarks hadronising as a particular excited charmed meson.
These fragmentation fractions agree with previous measurements
in $e^+e^-$ annihilations.

\begin{figure}[b]
\centerline{\psfig{file=gladilin_fig5.eps,width=2.8in}}
\caption{
The distributions of the
$M(\dspm K^0_s)$
for the $\dspm$ candidates.}
\label{fig5}
\end{figure}

The helicity angular distribution of the \dsonepm\ signal
was fitted
to a form $1+R\cos^2\alpha$, where $R$ is a helicity
parameter and
the helicity angle ($\alpha$) is defined as the angle between
the $K_S^0$ and $\pi_S$ momenta in the \dspm\ rest frame.
The unbinned likelihood fit yielded the value
$R=-0.53\pm0.32({\rm stat.})^{+0.05}_{-0.14}({\rm syst.})$.
This value is consistent with the helicity parameter value
$-0.23^{+0.40}_{-0.32}$ obtained
by the CLEO Collaboration for $\dsonepm$ meson
in the $D^{*0}K^{\pm}$ final state~\cite{cleo_d1s}.
The present measurement does not contradict the conclusion
of the CLEO Collaboration that $R=0$ and, thus, the spin-parity of
the \dsone\ is $1^+$. However, this measurement is also consistent
with $R=-1$, expected for spin-parities
$1^-$ and $2^+$~\cite{godfrey_kokoski}.
Recently, a precise measurements of the $\dsone$ helicity parameter,
$R=-0.70\pm0.03$, has been reported by the Belle
Collaboration~\cite{belle_d1s}.
The measured value has been interpreted as a result of the $D$ and $S$
waves mixture due to an interference of the $\dsone$ meson
with the recently discovered $D_{sJ}(2460)^+$ meson.

The ZEUS search for the radially excited $D^{*\prime\pm}$ meson,
reported by the DELPHI Collaboration~\cite{delphi_drad},
revealed
no signal~\cite{zeus_d1d2}.
The upper limit on
the product of the fraction of $c$ quarks hadronising as
a $D^{*\prime +}$ meson
and the branching ratio of the $D^{*\prime +}$ decay to
$D^{*+}\pi^+\pi^-$
was estimated to be $0.7\%$ ($95\%$ C.L.).
This limit is somewhat stronger than
the $0.9\%$ limit obtained by the OPAL Collaboration~\cite{opal_drad}.

\balance

High Energy Physics - Experiment (hep-ex)

Charm spectroscopy and exotics at ZEUS


The Detection of Silicate Emission from Quasars at 10 and 18 Microns

Abstract: We report the spectroscopic detection of silicate emission at 10 and 18 microns in five PG quasars, the first detection of these two features in galaxies outside the Local Group. This finding is consistent with the unification model for Active Galactic Nuclei (AGNs), which predicts that an AGN torus seen pole-on should show a silicate emission feature in the mid-infrared. The strengths of the detected silicate emission features range from 0.12 to 1.25 times the continuum at 10 microns and from 0.20 to 0.79 times the continuum at 18 microns. The silicate grain temperatures inferred from the ratio of 18-to-10 micron silicate features under the assumption of optically thin emission range from 140 to 220K.

\section{Introduction}

Active galactic nuclei (AGNs) are broadly classified in two
types.  Type 1 AGNs display broad hydrogen emission lines in
the optical, while Type 2 AGNs do not.  The AGN unification
model \citep[e.g.][]{ant93} ties these two types together.  An
optically and geometrically thick dusty torus surrounds a
central black hole, accretion disk and broad-line emission
region.  Sources viewed face-on are recognized as Type 1,
while edge-on sources are Type 2.  Observations at many
wavelengths support this scenario \citep[e.g.][]{am85}, but
until now there was a lack of constraining data in the
mid-infrared.

Silicate dust is a major component of the interstellar medium
in the Milky Way and other galaxies.  It produces two spectral
features in the infrared, a ``10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ feature'' which arises
from a Si$-$O stretching mode, and an ``18\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ feature'' from
an O$-$Si$-$O bending mode \citep[e.g.][]{kt73}.  Emission or
absorption from silicate dust is a dominant feature of the
infrared spectrum produced by most mass-losing evolved stars,
interstellar clouds, H{\sc ii} regions, and circumstellar disks
around young stellar objects (YSOs).  Silicate absorption is
also detected in dusty starburst galaxies \citep{rig99},
ultraluminous infrared galaxies \citep[ULIRGs;][]{gen98,
tra01, spo04, arm04}, and a sample of mid-infrared detected,
optically invisible, high-luminosity galaxies with redshifts
of 1.7 $<$ z $<$ 2.8 \citep{hou05}.

Type 2 AGNs show silicate dust in absorption as expected, but the
mid-infrared spectroscopic data for Type 1 AGNs have been scarce.
Recently, \cite{sie05} have presented spectra of two
quasars with just the 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ feature in emission.  While radiative
transfer models clearly predict silicates to be in emission for Type 1
sources \citep{ld93, pk92, pk93, gd94, err95, nie00, nie01, nie02,
vbd03, dvb05}, the absence of the silicate emission features at 10 and
18\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ has led to models without strong silicate emission.
\cite{ld93} and \cite{vbd03} have suggested that larger grains
dominate the grain-size distribution.  \cite{nie02}, instead,
proposed that if the torus is clumpy, the 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ silicate emission
can be sufficiently suppressed, although \cite{dvb05} have challenged
this suggestion.

Here, we report the first spectroscopic detection of both
silicate features at high signal to noise in five PG quasars.
Quasars are the high luminosity counterparts of Seyfert 1 galaxies,
showing the same broad emission lines that define Type 1 AGNs.

\section{Observations and data reduction}

The five sources presented in this letter (PG 0804+761, PG 1211+143,
PG 1351+640, I Zw 1 [=PG 0050+124] and 3C~273 [=PG 1226+023]) were
selected from a sample of 12 AGNs based
upon their prominant silicate emission features at both 10 and
18\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi. The observations were part of the Guaranteed Time
Observation (GTO) program of the Infrared Spectrograph
(IRS)\footnote{The IRS was a collaborative venture between Cornell
University and Ball Aerospace Corporation funded by NASA through the
Jet Propulsion Laboratory and the Ames Research Center}\citep{hou04}
on the Spitzer Space Telescope \citep{wer04}.  Table \ref{tab1} lists
the basic properties of these targets, along with their observation
dates and on-source integration times.  All of them have broad
H$\beta$ lines \citep{bg92} and therefore are classified as Type 1
AGNs. Two of them, 3C~273 and PG 1211 are radio-loud quasars and the
rest are radio-quiet \citep{san89}.  All five quasars are very
luminous, as can be seen from the bolometric luminosities in
Table\,\ref{tab1} \citep[from][]{san89}.

\begin{table}
\caption{Properties of Sources\label{tab1}}
\begin{tabular}{lccccc}
\colrule
\colrule
 & PG 1351 & I Zw 1 & 3C 273 & PG 1211 & PG 0804 \\
\colrule
AOR key                              & 3760640   & 3761920  & 4978176  & 3760896  & 9074944 \\
Date observed                        & 2004/4/17 & 2004/1/7 & 2004/1/6 & 2004/1/7 & 2004/3/1\\
Integration time (s)                 & 408       & 224      & 336      & 408      & 408\\
Redshift                             & 0.0881    & 0.0611   & 0.158    & 0.0808   & 0.0999 \\
Lum. distance (Mpc)\tablenotemark{a} & 398       & 271      & 749      & 363      & 455 \\
log($L_{\rm bol}/L_{\odot}$)         & 12.14     & 12.20    & 13.44    & 12.44    & 12.56\\
F$_{15\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi}$ (Jy)\tablenotemark{b} & 0.21      & 0.62     & 0.36     & 0.21     & 0.13\\
\colrule
\end{tabular}
\tablenotetext{a}{assuming H$_0$=71\,km\,s$^{-1}$\,Mpc$^{-1}$,
$\Omega_M$=0.27, $\Omega_{\Lambda}$=0.73, $\Omega_K$=0}
\tablenotetext{b}{rest frame}
\end{table}

The observations were made with the Short-Low (SL) and Long-Low (LL)
modules of the {\em IRS}.  The spectra were extracted from the
flatfielded images
provided by the Spitzer Science Center (pipeline version S11.0.2). The
images were background-subtracted by differencing the two SL apertures
and for LL, by differencing the two nod positions.  Spectra were then
extracted and calibrated using the {\em IRS} standard star HR 6348 for
SL and the stars HR 6348, HD 166780, and HD 173511 for LL
\citep{slo05}.

After extraction the orders were stitched, requiring order-to-order
scaling adjustments of less than 5\%.  Finally, the stitched spectra
were scaled to match the observed {\em IRAS} fluxes at 12 and
25\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ or {\em IRS} blue peak-up flux at 16~\ifmmode{\rm \mu m}\else{$\mu$m}\fi.  For 3C 273, we
did not attempt any scaling, since the source varies by up to 0.4
magnitudes at 10.6\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ on time scales of up to several years
\citep{neu99}.  The scaling factors we applied average 6\%, with
the largest scaling factor of 16\%. The observed sources are point
sources for both IRAS and IRS, therefore scaling to IRAS flux is
justified.  Figure \ref{fig1} presents the calibrated spectra.  Table
\ref{tab1} lists the 15\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ rest frame fluxes.

\begin{figure}
\begin{center}
\resizebox{\hsize}{!}{\includegraphics{f1.ps}}
\caption{The 5--35\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ IRS low resolution spectra of five PG
quasars showing silicates in emission.  The pixel-to-pixel
uncertainty in the data is less than 3\%.  Fringing effects with
amplitudes of up to 5\% are visible in LL order 1 (see text). For
each spectrum, three choices for the silicate-graphite continuum
are indicated by the dashed line and grey shaded area.
The spectra have been scaled for plotting purposes, with scaling
factors of:  8.0 for PG\,1351, 1.0 for I\,Zw1, 0.77 for 3C\,273
and PG\,1211 and 0.50 for PG\,0804.
\label{fig1}}
\end{center}
\end{figure}

In some spectra, most notably in IZw1, residual fringing, which
are due to minor pointing errors, appears between 21 and 30\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ (the
1st order of the LL module).  We have not defringed the spectra.
Residual artifacts also appear in some spectra in SL order 1
(7.5--14\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi).  These take the form of a roughly sinusoidal
deviation from the actual spectrum and are most obvious in the
spectrum of 3C 273 at 8--10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi, where they have a maximum strength
of $\sim$3\% of the total flux.  Because 3C 273 has the weakest
silicate emission feature in our sample, this artifact is even more
pronounced in figures showing the continuum-subtracted spectra.  We
estimate uncertainties by comparing the spectra in the two nod
positions, and these artifacts generally manifest themselves as a
disagreement between the nod positions.  The uncertainties only become
apparant in Figures \ref{fig2} and \ref{fig3} shortward of 11\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\
(dark shaded area). None of these artifacts change the conclusions
reported below.

\begin{figure}
\begin{center}
\resizebox{\hsize}{!}{\includegraphics{f2.ps}}
\caption{Continuum-subtracted silicate emission spectra of the
five PG quasars. The spectra have been normalized and offset
from each other for plotting purposes.  Fringing effects are
visible in LL order 1 (see text).  The shaded areas indicate the
uncertainty in the silicate emission spectra as introduced by the
uncertainty in the observed spectrum and choice of continuum.
Vertical dotted lines indicate the rest wavelengths of emission
features discussed in the text.
\label{fig2}}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\resizebox{\hsize}{!}{\includegraphics{f3.ps}}
\caption{
The opacity spectra for the five PG quasars compared to the
scaled synthetic silicate opacity curve of \cite{wein01} and
\cite{li01}.  The opacity spectra from the quasars are derived
using the adopted silicate grain temperature and scaled and offset
for plotting purposes.  The shaded areas indicate the uncertainty
in the silicate profile as introduced by the uncertainty in the
observed spectrum and choice of continuum.
\label{fig3}}
\end{center}
\end{figure}

\section{Analysis}

The spectra in Figure \ref{fig1} are dominated by two broad emission
features, which we identify as the Si--O stretching mode and the
O--Si--O bending mode
in silicate dust grains.  The strength of the features
is best appreciated by comparing the peak fluxes to the underlying
silicate-graphite ``continuum''. For PG 1351, I Zw 1 and PG 1211, we
define the continuum using a spline interpolation between the
5--8\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ spectrum and a point at the extreme red end of the IRS
spectral range, where a change in slope indicates the end of the
18\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ silicate feature.  For 3C 273 and PG 0804, the spectral
structure at the red end of the IRS spectral range does not show a
similar change in slope.  This may indicate that dust emission is
still contributing at these wavelengths.  Therefore, we have assumed
that the continuum at the red end runs somewhat below the observed
spectrum.  Figure \ref{fig1} plots the resulting spline-interpolated
continua (dashed lines), while Figure \ref{fig2} shows the resulting
continuum-subtracted spectra.

Comparison of the silicate emission spectra in Figure \ref{fig2}
reveals that the ratio of the 18 and 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ features varies
significantly within the sample.  It is highest for 3C 273 and lowest
for PG 0804 (see Table \ref{tab2}).  Also the onset and center of the
10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ features shift significantly from one source to the next, 3C
273 and PG 1211 seem to appear at longer wavelengths. Interestingly,
these two sources happen to be the only radio loud quasars in our
sample. In order to test whether these differences can be attributed
to our choice of continuum, for each source we compared the silicate
emission spectra for two other choices of underlying silicate-graphite
continuum: one above and one below the adopted continuum, enclosing
the shaded areas in Figure \ref{fig1}.  As illustrated by the shaded
areas in Figure \ref{fig2}, this results in an uncertainty of up to
30\% in the peak flux of the 18\,$\mu$m feature and in the ratio of
the 18 to 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ features. In contrast, the onset and center of the
10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ features are not affected by different choices of continuum.

\begin{table}
\caption{Measured quantities\label{tab2}}
\begin{tabular}{lccccc}
\colrule
\colrule
 & PG 1351 & I Zw 1 & 3C 273 & PG 1211 & PG 0804 \\
\colrule
10$\mu$m-feature-to-cont. ratio\tablenotemark{a}       & 1.25  & 0.38  & 0.12  & 0.55  & 0.60 \\
18$\mu$m-feature-to-cont. ratio\tablenotemark{a}       & 0.79  & 0.20  & 0.36  & 0.56  & 0.56 \\
18$\mu$m\,/\,10$\mu$m feature ratio\tablenotemark{a}   & 1.75  & 1.31  & 3.34  & 1.89  & 1.08 \\
F(30$\mu$m)/F(5.5$\mu$m)              & 11.0  & 6.86  & 2.53  & 2.99  & 1.95 \\
L$_{\rm sil}$/L$_{\rm bol}$           & 0.087 & 0.040 & 0.009 & 0.023 & 0.021\\
10/18$\mu$m cont. color temp (K)\tablenotemark{b}      & 250   & 270   & 320   & 310   & 340  \\
Silicate temperature (K)              & 200   & 220   & 140   & 180   & 210  \\
\colrule
\end{tabular}
\tablenotetext{a}{evaluated at the peak of the feature}
\tablenotetext{b}{the color temperature is computed from the adopted continuum at 10 and 18\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ }
\end{table}

The presence of the silicate features at 10 and 18\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\
allows us to estimate the mean tempeature of the emitting grains by
performing an unweighted $\chi^2$ fit to the continuum-subtracted
spectrum. This temperature is only physically meaningful under the
assumption that the feature emitting region is optically thin. Using
the continuum-subtracted silicate opacity (see top panel in Figure
\ref{fig3}) derived from \cite{wein01} and \cite{li01}, we force the
product of a blackbody and the silicate feature opacities to fit the
observed continuum-subtracted silicate features. This reduces the
problem to a single free parameter: the grain temperature, and avoids
concerns related to a possible nonthermal contribution to the continuum.
The resulting silicate grain temperatures (T$_{\rm sil}$) range between
140\,K for 3C 273 and 220\,K for I Zw 1, see Table \ref{tab2}.  These
temperatures will change by as much as 15\,K depending on the choice
for the underlying silicate-graphite continuum.  Uncertainties in the
fitting process may amount to as much as another 10\,K.

Using the derived silicate grain temperature (T$_{\rm sil}$),
it is possible to compute the silicate opacity spectrum by
dividing the observed silicate emission spectrum by the
blackbody spectrum with temperature T$_{\rm sil}$.
Figure \ref{fig3} presents the resulting opacity spectra
for our sample.

The top panel in Figure \ref{fig3} shows the synthetic silicate
opacity spectrum after subtraction of the silicate-graphite continuum
component.  The lower panels show the opacity spectra for the five
quasars overlaid on the synthetic profile.  The quasar profiles
generally show good agreement at wavelengths beyond 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi.  Below
10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi, only PG 1351 and I Zw 1 match the synthetic profile well.
In the other three spectra, the blue wing of the 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ appears
weak, especially in 3C 273 and PG 1211.  Apart from the artifacts in
the 8--10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ region of the spectrum of 3C 273 (indicated by the
large error bars in Figure \ref{fig2}; see \S 2) and the presence of
atomic lines and PAH emission in some spectra, the 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ band has
a smooth appearance, showing no sign of a departure from amorphous
grain structure.  Evidence for some crystalline grains do appear at
longer wavelengths, most notably in the spectrum of PG 1351 at
23\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi.

The weakness of the blue wing of the 10\,$\mu$m silicate opacity
profile in three of our sources appears to be significant and cannot
be explained by the uncertainties in the observed spectrum (see the
error bars in Figure \ref{fig3}). Grain size and composition,
geometry, optically thick radiative transfer or a combination of these
factors may explain the observed weakness. If our optically thin
assumption is invalid, then the temperatures we have inferred cannot
be interpreted as real dust temperatures. \cite{jaf04} reported a
similar deficiency in the blue wing of the 10um absorption feature in
NGC1068, which they attribute to a different silicate grain
composition.

A close inspection of the spectra in Figures \ref{fig1} and \ref{fig2}
reveals numerous other features.  Those visible in at least one
spectrum in the sample include the PAH emission features at 6.2, 7.7
and 11.2\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi, and emission lines at 10.5\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ ([S{\sc iv}]),
12.8\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ ([Ne{\sc ii}]), 14.3\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ ([Ne{\sc v}]), 15.6\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\
([Ne{\sc iii}]), 18.7\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ ([S{\sc iii}]), and 25.9\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ ([O{\sc
iv}]).

\section{Comparison with Models}

Direct comparisons of our observations to the predictions of published
AGN radiative transfer models are difficult given the limited resolution
of model spectra in the 5--40\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ range of the IRS.  Additionally,
the model output is not tuned to be compared to the observed
properties presented in Table\,\ref{tab2}. However, we can measure the
feature-to-continuum ratio at 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ from the figures presented by
\cite{ld93, gd94, nie02, dvb05}.  For the optically thin dust models
of \cite{ld93} the values range from 0.1 for models where the
grain-size distribution is dominated by larger grains to 1.5 for
models with a standard grain-size distribution.  In the face-on cases
modeled by \cite{gd94}, the values range from near zero for very thick and
compact configurations to 1.0 for more extended torus models.  These
results are similar to those by \cite{nie02} and \cite{dvb05}, but
their highest 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ feature-to-continuum ratios are around 0.7 and
0.4, respectively.  Given the observed range of 0.12 to 1.25 within
our limited sample (Table \ref{tab2}), the model predictions for the
10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ feature-to-continuum ratio are in good agreement. A similar
comparison for the 18\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ feature-to-continuum ratio could not be
undertaken, since the features are not as apparent as the 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\
features in the AGN models.

\section{Discussion and Conclusion}

We have detected the 10 and 18\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ silicate emission
features in five PG quasars --- the first detection of both
silicate features in emission in galaxies outside the Local
Group and in AGNs in particular.  Given the existence of
published mid-infrared photometric and spectroscopic data for
these quasars, our finding may appear surprising at first.
\cite{roc91} obtained 8--13\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ spectra of 3C 273 and
I Zw 1, but the low contrast of the silicate feature to the
strong continnum and the lack of coverage outside the
N-band window prevented them from identifying silicate emission.
The {\it PHT-S} spectrometer on {\em ISO}\ \citep{lem96} obtained
spectra of 3C 273, I Zw 1 and PG 0804 \citep{rig99},
but the limited wavelength coverage of {\it PHT-S} (only to
11.6\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi) and the redshifts of the quasars allowed only the
blue side of the 10~\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ feature to be visible.  Because of
these difficulties, silicate emission had remained unidentified
in these five sources until now.

Our detection of silicate emission in Type 1 AGNs provides support for
the unified AGN model \citep{ant93}.  An AGN torus viewed pole-on
should offer an unobstructed line of sight to the torus dust, the
surface of which is heated to temperatures of several hundred to a
thousand Kelvin by the radiation from the central engine.  Emission
from this hot surface should produce a silicate feature in emission
under most circumstances.

Assuming a simple optically thin geometry, we find that the silicate
grain temperature in our sources ranges from 140 to 220\,K, which is
well below the dust sublimation temperature. Detailed comparison of
the derived opacity profiles for our quasars indicates that our simple
model cannot satisfactorally explain the absence of a blue wing in the
10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ feature for three of our quasars. This likely points to a
more complex geometry, involving radiative transfer within a
(partially) optically thick environment -- probably an AGN torus. It
may also point to differences in the grain size or composition.  Our
derived silicate grain temperatures should, therefore, only be used as
a characterization of the observed spectra, not as a means to infer
the actual temperature and geometry of the AGN.

From our limited sample of PG quasars (three in addition to the sample
discussed in this Letter) and a few Seyfert 1 galaxies, it is clear
that 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ silicate emission is not a ubiquitous feature in
mid-infared spectra of Type 1 AGNs.  Only 5 out of 12
continuum-dominated AGNs show this feature prominently. Interestingly,
the quasars which show a 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ silicate emission feature all have
bolometric luminosities higher than 10$^{12}$ L$_{\odot}$, while the
other Type 1 AGNs are less luminous. Models which can
explain the presence of the 10\,\ifmmode{\rm \mu m}\else{$\mu$m}\fi\ silicate emission feature in
some type 1 AGNs will also need to explain their absence in others.

\acknowledgments

The authors wish to thank Bill Forrest, Elise Furlan and Terry Herter
for their useful discussions.  Support for this work was provided by
NASA through Contract Number 1257184 issued by the Jet Propulsion
Laboratory, California Institute of Technology under NASA contract
1407.  HWWS was supported under this contract through the
Spitzer Space Telescope Fellowship Program.