Pion and photon production in heavy ion collisions

Abstract: Measurement of neutral pions and direct photons are closely connected experimentally, on the other hand they probe quite different aspects of relativistic heavy ion collisions. In this short review of the $\pi^0$ results from the PHENIX experiment at RHIC our focus is on the $\phi$-integrated nuclear modification factor, its energy and system size dependence, and the impact of these results on parton energy loss models. We also discuss the current status of high $p_T$ and thermal direct photon measurements both in $p$ + $p$ and Au+Au collisions. Recognizing the advantages of measuring not only the "signal", but also all the "references" needed for proper interpretation in the same experiments (with same or similar systematics) we argue that RHIC should regularly include $d$ + A and even $d$ + $d$ collisions into its system size and energy scan.

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

\begin{figure}
\begin{centering}
\includegraphics[width=.6\textwidth]{gpe_RAA.eps}
\caption{``Historic'' nuclear modification factors $R_{AA}$ {\it vs}
  $p_T$ in  200GeV/A Au+Au  collisions.  Squares: direct photons (no
  apparent  suppression~\cite{ppg042}, published 2005).
  Triangles and circles: preliminary results on $\pi^0$ and $\eta$,
  also from 2005. }
\label{fig:gpe_RAA}
\end{centering}
\end{figure}

One of the {\it predicted} signals of the formation of a quark-gluon
plasma (QGP) in relativistic heavy ion collisions was the
in-medium energy loss of hard scattered partons, leading to ``jet
quenching'' and manifesting itself in a drastic reduction of high
$p_T$ hadron production in relativistic heavy ion collisions, as
compared to the production rate in $p$+$p$ scaled by the nuclear thickness
function $T_{AA}$ to account for the increased probability of hard
scattering.  The effect has been characterized by the
{\it nuclear modification  factor} $R_{AA}$ defined as

$$
R_{AA}^{\alpha}(p_T) = \frac{1}{N_{evt}}
\frac{d^2N_{AA}^{\alpha}/dp_Tdy}
     {\left\langle T_{AA} \right\rangle d^2\sigma_{pp}^{\alpha}/dp_Tdy}
$$

for an arbitrary particle type $\alpha$.
In the absence of nuclear effects $R_{AA}=1$, otherwise it is enhanced
($R_{AA}>1$, like in the ``Cronin-effect'') or suppressed ($R_{AA}<1$
like in jet quenching.)  The first major discovery at RHIC was that in
central Au+Au collisions hadrons were strongly
suppressed~\cite{ppg003,ppg014}
(see Fig.~\ref{fig:gpe_RAA}), in accordance with the predicted energy
loss of partons in the hot, dense medium.  Lack of suppression in
d+Au collisions~\cite{ppg028} (where we assume no such medium is formed)
supported the conclusion that the jet suppression was indeed a final
state effect, as opposed to some drastic changes in the initial state
(parton distribution functions).
The first observation that direct
photons were not suppressed, not even in the most central Au+Au
collisions~\cite{ppg042} made this first, qualitative picture
coherent, by showing that the crucial step in forming $R_{AA}$ -
scaling $p$+$p$ cross sections by $\left\langle T_{AA} \right\rangle$ to
obtain the expected Au+Au yields - makes sense: photons will not
interact with the medium ($\alpha_{em} << \alpha_s$), so if the
$\left\langle T_{AA} \right\rangle$ scaling is the right concept, the
direct photon $R_{AA}$ should be unity - and indeed, it is
(see Fig~\ref{fig:gpe_RAA}).

While very useful and providing many early insights, causing a flurry
of theoretical activities rarely seen before, the above picture
turned out to be only rudimentary.  As RHIC experiments moved
``from the discovery phase to the exploration phase'' more precise data
became available, with higher statistics, extended $p_T$ ranges
and more complex physics quantities (like two- and three-particle
correlations, azimuthal anisotropies, {\it etc}.
Equally important, RHIC launched a systematic energy and colliding
species scan.  This made measuring excitation functions in the same
collider and within the same experiments possible.  The collider
environment ensures that multiplicities at midrapidity rise only
slowly with $\sqrt{s_{NN}}$, and measuring in the same experiment
means that most systematic errors will be similar - even if individual
measurements are not very precise, the evolution with energy and/or
system size can be established quite accurately.  Increased accuracy
of the data made comparisons with theories more
meaningful, put some constraints on free parameters~\cite{ppg079} and
gave the first experimental indications at what region of the phase
diagram the transition to the strongly coupled quark gluon plasma
might occur.  Almost needless to say that in the process our initial,
relatively simple picture became richer (much more complicated).

\section{Pion suppression systematics}
\label{sec:pionsuppression}

So far four aspects of the $\phi$-integrated  $R_{AA}$ systematics
have been studied at RHIC.  First, the magnitude and shape of
$R_{AA}(p_T)$ at the highest
available transverse momenta - namely, whether it remains constant,
consistent with a constant fractional energy loss picture, or rises
with increasing $p_T$, as predicted by most parton energy loss
models~\cite{ppg079}.  These studies gave the first quantitative
constraints on free parameters of models, like initial gluon density
or the transport coefficient of the PQM~\cite{PQM2007} model
$\left\langle \hat{q}=\mu^2/\lambda \right\rangle =
13.2^{+2.1}_{-3.2}$GeV$^2$/fm, where $\mu$ is the average transverse
momentum transfer to the medium per mean free path $\lambda$.
PQM - which includes radiative loss, a static medium, no initial state
multiple scattering and unmodified PDFs - predicts a slow rise
of $R_{AA}$ with $p_T$, as do other models with very different
assumptions~\cite{ppg079}.  In light of this it is interesting, that
the data are best described by a constant fit which in turn for the
power law spectra of high $p_T$ pions would imply a constant fractional
energy loss (constant $S_{loss}=\Delta p_T / p_T$ fractional shift
of the spectrum).  These statements are based on data taken in
RHIC Run-4 (2004), which still have large errors on the high $p_T$
points, results from the much higher statistics 2007 data should be
available soon.

Second, the evolution of $R_{AA}$ with collision centrality has been
studied.  One very interesting and non-trivial observation was that the
shape of the high $p_T$ $\pi^0$ spectra is virtually unchanged from
$p$+$p$ to the most central Au+Au collisions.  PHENIX found~\cite{ppg080}
that the exponent of the power law function fitting the high $p_T$
part of the spectra ($\propto p_T^n$) varies from
$n=-8.22\pm0.09$ in $p$+$p$ to $n=-8.00\pm0.012$ in the most central
(0-5\%) Au+Au collisions.  This combined with the fact that $R_{AA}$
in the most central collisions is essentially flat (constant) means
that $R_{AA}$ is constant at {\it all} centralities to the best of our
current knowledge which in turn means that $R_{AA}$ integrated above a
certain $p_T$ is a useful quantity to characterize the onset and
evolution of jet quenching.  This is shown on Fig.~\ref{fig:int_raa}.
Fitting the points with a function motivated by the constant
fractional energy loss ({\it i.e} the energy loss increases with $p_T$)
$R_{AA} = (1 - S_0 N^a_{part})^{n-2}$ the exponent $a$ is found to be
$a=0.58\pm0.07$ for $p_T>5$GeV/c and $a=0.56\pm0.10$ for $p_T>10$GeV/c
which is in reasonable (albeit not perfect) agreement with the
$a \approx 2/3$ predictions of the PQM~\cite{PQM2007} and GLV~\cite{GLV2000}
 models.  Note that in the ``constant fractional
energy loss'' picture $R_{AA}$ depends on the slope (exponent $n$) of
the power-law $p_T$ spectrum which in turn depends on the per-nucleon
collision energy $\sqrt{s_{NN}}$.

\begin{figure}
\begin{centering}
\includegraphics[width=.6\textwidth]{int_raa.eps}
\caption{Integrated $R_{AA}$ for $\pi^0$ as a function of collision
  centrality expressed in terms of $N_{part}$.  The last two points
  correspond to overlapping centrality bins 0-10\% and 0-5\%.  See
  text for an explanation of the fits.}
\label{fig:int_raa}
\end{centering}
\end{figure}

Therefore, the third issue is what can be learned from the energy
dependence of $R_{AA}$, and the first answer is given on
Fig.~\ref{fig:raa_CuCu} which shows the nuclear modification factors
measured at RHIC in the most central Cu+Cu collisions at
$\sqrt{s_{NN}}=$22.4,62.4 and 200GeV.  Instead of suppression the
lowest energy data show a substantial {\it enhancement} of $R_{AA}$ in
line with the well-known ``Cronin-effect'' (and within errors this
enhancement is {\it independent} of centrality~\cite{ppg084}) whereas at
$\sqrt{s_{NN}}=$62.4GeV the pions are already significantly
suppressed.  Of course it does {\it not} follow that the onset of jet
quenching - presumably due to QGP formation - happens necessarily
between $\sqrt{s_{NN}}=$22.4 and 62.4GeV; all we can say is that the
quenching mechanism overtakes the Cronin-effect in this region.
Being a consequence of the initial $k_T$ smearing of partons the
Cronin-effect is present already in cold nuclei and could be unfolded
from the $\sqrt{s_{NN}}$-dependence of $R_{AA}$ if $p$+A ($d$+A) data were
available from RHIC in the entire available collision energy range.
In the author's opinion taking $d$+A data (along with A+A and $p$+$p$ for
reference) should be part of the planned ``energy scan'' at RHIC.
- Since pions are strongly suppressed in Cu+Cu already at
$\sqrt{s_{NN}}=$62.4GeV, it is not surprising that they are also
suppressed in Au+Au, in fact, almost as strongly as at 130 or 200GeV.

Finally, the fourth issue is how does $R_{AA}$ depend on the
(size of the) colliding systems at any given energy.  A comparison of
$\sqrt{s_{NN}}=$200GeV central Cu+Cu collisions to mid-central Au+Au
(such that the number of participating nucleons $N_{part}$ is the same
in both cases albeit the geometry is not) shows similar values
of $R_{AA}$.  It would be interesting to know which is the lightest,
symmetric heavy ion system where jet quenching at
$\sqrt{s_{NN}}=$200GeV would disappear.

So far we discussed only the $\phi$-integrated $R_{AA}$.  For any
given dataset neglecting azimuthal dependence results in higher
precision of the measurement itself (better statistics and no
systematic errors related to the determination of the reaction plane)
but it also hides any information encoded in azimuthal anisotropies
of $R_{AA}$ due to the overlap geometry of the colliding nuclei.
Such anisotropies clearly exist: after all, some of the most
spectacular observations at RHIC are the large ``elliptic flow'' of
hadrons and the quark number scaling thereof.  Indeed, when studying
$R_{AA}(\Delta\Phi)$ where $\Delta\Phi$ is the azimuthal angle
with respect to the reaction plane, determined event-by-event,
we observe a rich structure~\cite{ppg054}: jet quenching varies
substantially with the average pathlength of the parton in the medium.
Back-to-back dihadron (or more ambitiously, dijet) correlations, and
particularly the ``golden channel'', photon-hadron (photon-jet)
correlations have also been advocated - and now measured -
as tools of ``jet tomography''.  It has been argued~\cite{wang2007}
that single hadron spectra at high $p_T$ are dominated by surface
emission, therefore, $R_{AA}$ is barely sensitive to the medium deep
inside the collision volume.  On the other hand back-to-back
correlations force one of the partons to go through the medium (with
the exception of ``tangential'' emissions) thus the dihadron
suppression $I_{AA}$ is more sensitive to the input parameters of
theories like initial energy or gluon density.  This is certainly true
as long as the uncertainties of both theory and experimental data are
negligible: a multidifferential quantity always provides more
constraints.  However, theories usually have uncertainties (even if
often unstated, regrettably), and experimental data always have them
and it is hard to imagine a situation in which
{\it using a given dataset} a multidifferential quantity could be
measured more precisely than a single differential.  Just the
opposite.  The {\it quantitative} question then is: does the increased
sensitivity on the theory side compensate for the inevitable loss of
precision in the data?  The answer varies from case to case.

\begin{figure}
\begin{centering}
\includegraphics[width=.6\textwidth]{raa_CuCu_energy_Cc08.eps}
\caption{Nuclear modification factors for the most central Cu+Cu
  collisions~\cite{ppg084} for $\sqrt{s_{NN}}=$22.4,62.4 and 200GeV
  compared to theoretical calculations~\cite{vitev2006}.}
\label{fig:raa_CuCu}
\end{centering}
\end{figure}

\section{Direct photons}
\label{sec:photons}

Direct photons - defined as those not originating from final state hadron
decays - offer unique possibilities to study heavy ion collisions.
First, there are physics mechanisms that produce direct photons at
each stage of the collision, thus they have the potential to provide
information on the entire evolution, ``history'' of the collision.
Equally important, once produced, they leave the interaction region
mostly unaltered ($\alpha_{em}<<\alpha_s$), so the information they
carry is unbiased.  Unfortunately, the very same properties make them
challenging to measure (photons from hadron decays provide a large
background, particularly at low $p_T$) and their ``message'' is hard
to decypher because it is difficult (if not impossible) to deconvolute
the contributions of individual mechanisms to the observed direct
photon spectrum.

On Fig.~\ref{fig:gpe_RAA} we already have seen the first measurement
of direct photon $R_{AA}$ in central Au+Au collisions at
$\sqrt{s_{NN}}=$200GeV, showing that in the region where photons from
primordial hard scattering were expected to dominate
($p_T>$5-6GeV/$c$) $R_{AA}$ is indeed consistent with unity.  This was
consistent with the picture that both the PDFs and the production
cross section are unchanged, the only difference between $p$+$p$ and
Au+Au as far as hard scattering is concerned is the increased
probability of an otherwise rare process; this increased
probability is properly described by a straightforward geometric
overlap integral (the nuclear thickness function $T_{AA}$), and the
photons produced escape the collision region unaltered even if a
dense, strongly interacting medium is formed which quenches parton
jets.  Although we have higher quality data now (in an extended range)
and our understanding of the processes generating photons is much more
nuanced, we should emphasize that at a very basic level, in first
approximation, qualitatively the above picture is still correct.

What are the currently known/assumed sources of direct photons in
$p$+$p$ and Au+Au collisions?
At leading order two processes generate photons, the
quark-quark annihilation ($q\bar{q} \rightarrow g\gamma$), which is
suppressed at high $p_T$ due to the lack of valence antiquarks, and
the dominant quark-gluon Compton-scattering ($qg \rightarrow q\gamma$).
Note that in the latter the photon balances the momentum of the
original quark (modulo the initial $k_T$ smearing of the partons)
and it comes out freely from the medium, whereas the quark will lose a
significant fraction of its energy. Therefore, in back-to-back
photon-jet correlations the photon  ``calibrates'' the energy scale of
the jet, which is important in ``jet tomography'' as well as in
studying possible modifications of the fragmentation function in the
presence of a medium.  The probe itself is very clean and well
understood - unfortunately the rates are quite low, due to the
$\alpha_{em}$ coupling.

\begin{figure}
\begin{centering}
\includegraphics[width=.6\textwidth]{raa_qm06_cent0010_compArleo.eps}
\caption{Direct photon $R_{AA}$ in central Au+Au collisions.  The
  denominator is a fit to the PHENIX preliminary Run-5 $p$+$p$ data.
  Theoretical curves are LO calculations from~\cite{arleo2006}.  The
  solid line shows the pure isospin effect, the dash-dotted line is
  the combined effect of isospin and (anti)shadowing (EKS), finally
  the band at the botoom combines isospin and antishadowing with
  photon quenching.  }
\label{fig:raa_qm06}
\end{centering}
\end{figure}

At higher order fragmentation (or Bremsstrahlung) photons contribute
both in $p$+$p$ and heavy ion collisions.  In principle Compton and
fragmentation photons are distinguishable (at least in $p$+$p$) because
fragmentation photons should have enhanced, jet-like activity in a
small cone around them
(tracks or energy deposit that exceeds the level
justified by the underlying event),
whereas LO Compton photons should have none
in their neighborhood: they are
``isolated''.  There are predictions that at high $p_T$ fragmentation
photons will be a larger fraction of all direct photons in A+A
collisions than in $p$+$p$.  This is another reason why is it so
important to understand $p$+$p$ processes as well as possible: they
will serve as reference when we try to find new phenomena (new sources
of direct photons) in A+A.  First direct photon cross-sections in
$\sqrt{s}=$200GeV $p$+$p$ collisions at RHIC were published
in~\cite{ppg049} followed by~\cite{ppg060} with extended $p_T$ range,
small systematic errors and detailed comparison to NLO pQCD
calculations.  Above $p_T>5$GeV/$c$ where the systematic errors of
theory and experiment are comparable, the data are well described by
the theory, including the fraction of isolated photons

Is $p$+$p$ the right reference at all?  For pions it certainly is due to
the isospin symmetry of the protons and neutrons of the nucleus.  The
same is not true for direct photons: due to the $\alpha_{em}$ coupling
the photon cross section is proportional to the sum of the squared
quark charges $\Sigma e^2_q$ which is different for $p$ and $n$.  This
``isospin effect'' causes a trivial shift of the direct photon
$R_{AA}$~\cite{arleo2006},
depending on the nucleus and to some small extent on
centrality; in the absence of any other medium or nuclear effect
$R_{AA}$ should go about 10-15\% below unity in Au+Au collisions
(see the solid curve on Fig.~\ref{fig:raa_qm06}, which shows the
latest, still preliminary PHENIX results for photon $R_{AA}$ in the
most central Au+Au collisions).  Note that this
trivial effect may either mask or enhance other mechanisms affecting
$R_{AA}$.  The clean way to generate direct photon $R_{AA}$ would be
to compare to a proper mixture of photon cross sections in
$p$+$p$, $p$+$n$ and $n$+$n$ collisions,
which could be readily available at RHIC from
$d+d$ collisions, since these three types of collisions could easily
be tagged.

\begin{figure}
\begin{centering}
\includegraphics[width=.6\textwidth]{thermal.eps}
\caption{The fraction of the direct photon component as
  a function of $p_T$ in (a) $p$+$p$ and (b) Au+Au (min. bias).
  The error bars and the error band represent the statistical and
  systematic uncertainties, respectively. The curves are
  from a NLO pQCD calculation.}
\label{fig:thermal}
\end{centering}
\end{figure}

A recent overview of the possible (additional) photon production
mechanisms in heavy ion collisions is given in~\cite{turbide2008}.
At medium $p_T$ (corresponding to the $0.1<x<0.2$ region)
anti-shadowing should enhance the photon yield.  The large parton
energy loss (deduced from $\pi^0$ suppression) means that
fragmentation photons could be suppressed in the entire known $p_T$
range, {\it unless} the fragmentation functions themselves change in
the presence of a medium (which is hard, but not impossible to measure
in Au+Au).  On the other hand the jet-photon conversion mechanism may
increase high $p_T$ photon production.  In this process a hard
scattered quark interacts with a gluon or antiquark of the medium and
the outcoming photon carries essentially the full momentum of the
original parton (hence the name).  Disentangling the individual
contributions of these processes will not be easy, although we should
point out that some mechanisms produce isolated, others give
non-isolated photons, and their azimuthal
asymmetries are also different (for instance jet-photon conversion
should produce isolated photons with a negative ``elliptic flow''
$v_2<0$ which is quite unique since all other particles tested so far
at RHIC exhibited $v_2>0$).

Low $p_T$ real photons (thermal region) are very hard to measure due
to the overwhelming hadron decay background, however, ``quasi-real''
photons are accessible as low mass $e^+e^-$ pairs.  Using this
technique PHENIX measured the fraction of direct photons in the total
inclusive photon spectrum~\cite{thermal_photons}
(see Fig.~\ref{fig:thermal}) and found that while in $p$+$p$ the
results are consistent with NLO pQCD calculations (there are no
unexpected sources), in Au+Au there is a statistically significant
10\% excess at low $p_T$.

In summary, after the early, exciting discoveries RHIC experiments
entered the era of exploration.  Data became more precise enabling
detailed comparisons to theories and setting meaningful constraints on
their free parameters.  Taking advantage of the flexibility of the
accelerator a deliberate study of energy- and system-size dependence
of the new phenomena is underway.  The lively cooperation of theorists
and experimentalists (strengthened by some new initiatives) will
certainly lead to our much deeper understanding of the strongly
coupled quark-gluon plasma.


New Near-Infrared Spectroscopy of the High Redshift Quasar B 1422+231 at z=3.62

Abstract: We present new near-infrared (rest-frame UV-to-optical) spectra of the high redshift, gravitationally lensed quasar B 1422+231 (z=3.62). Diagnostic emission lines of FeII, [OIII]5007, and Hb, commonly used to determine the excitation, ionization, and chemical abundances of radio-quiet and radio-loud quasars, were detected. Our new data show that the ratio FeII(UV)/Hb=18.1+-4.6 and FeII(optical)/Hb=2.3+-0.6 are higher than those reported by Kawara et al. (1996) by factors of 1.6 and 3.3, respectively, although the ration [OIII]5007/Hb=0.19+-0.02 is nearly the same between the two measurements. The discrepancy of the line flux ratios between the measurements is likely due to improved data and fitting procedures rather that to intrinsic variability. While approximately half of the high-z quasars observed to date have much more extreme FeII(optical)/Hb ratios, the line ratio measured for B 422+231 are consistent with the observed range of FeII(optical) ratios of low-z quasars.

\section{INTRODUCTION}

Since near-infrared (NIR) spectroscopy of high-redshift ($z > 2$) quasars
provides information about their rest-frame optical spectroscopic properties,
it is now possible to systematically compare the spectroscopic  properties
(e.g., excitation, ionization, and chemical abundances) of
high-$z$ and low-$z$ quasars.
Although  the first NIR spectroscopic  observations of high-$z$ quasars
were made nearly two decades ago (Hyland, Becklin, \& Neugebauer 1978;
Puetter et al.\ 1981; Soifer et al.\ 1981), high-quality
NIR spectra have been obtained for only $\sim${\thinspace}10 high-$z$ quasars to date
(e.g., Carswell et al.\ 1991;
\markcite{Hill93}Hill et al.\ 1993; \markcite{Elston94}Elston et al.\  1994;
Kawara et al.\ 1996 [hereafter K96]; Taniguchi et al.\ 1997a;
Murayama et al.\ 1998; see also Taniguchi et al.\ 1997b).
However, despite the relatively small number of objects observed,
several very interesting properties of high-$z$ quasars have emerged.
In particular, Hill et al.\ (1993) and Elston et al.\ (1994; hereafter ETH)
suggested that unusually strong optical \ion{Fe}{2} emitter may be
common in the high-$z$ universe ($2<z<3.4$); to date,
five out of eight high-$z$ quasars
observed have very strong optical \ion{Fe}{2} emission with
EW(\ion{Fe}{2})/EW(H$\beta$)$\gtrsim 1$ (see Murayama et al.\ 1998).
By comparison, a much
smaller fraction of far-infrared (FIR) selected AGN ($L_{\rm FIR} \gtrsim
10^{11} L_{\sun}$) have strong optical \ion{Fe}{2} emission
(cf.\ L\'{\i}pari et al.\ 1993).

Recently, we obtained NIR spectra of the  radio-loud, flat-spectrum, high-$z$ quasar
B{\thinspace}1422+231 (Patnaik et al.\ 1992; Lawrence et al.\ 1992)
using the Mayall 4{\thinspace}m telescope at Kitt Peak National Observatory (KPNO) (K96).
Although this quasar is at $z=3.62$, its optical magnitude is
sufficiently bright ($m_{\rm r}$ =15.5; Yee \& Bechtold 1996), due to
gravitational lensing (Patnaik et al.\ 1992; Lawrence et al.\ 1992),
to allow it to be studied by NIR spectroscopy.
The NIR spectra show that the flux ratio
$F(\mbox{\ion{Fe}{2} }{\thinspace}\lambda\lambda\mbox{4434--4684})/F({\rm H}\beta)$
is much less than that of the other high-$z$ quasars (K96),
and in fact is  similar to those of radio-loud, flat-spectrum, low-$z$ quasars
with ``normal'' optical \ion{Fe}{2} emission.
This suggested that high-$z$ quasars may exhibit a range of values of
$F(\mbox{\ion{Fe}{2} }{\thinspace}\lambda\lambda\mbox{4434--4684})/F({\rm H}\beta)$
similar to what has been observed for low-$z$ quasars.
In order to confirm this result, we have obtained new NIR spectra using
the University of Hawaii (UH) 2.2{\thinspace}m telescope. In this paper,
we present our new NIR spectroscopy of B{\thinspace}1422+231 and compare it with our previous
measurements.

\section{OBSERVATIONS AND DATA REDUCTION}

We observed B{\thinspace}1422+231 on 1996, April 7 (UT) using the K-band
spectrograph (KSPEC; \markcite{Hodapp94}Hodapp et al.\ 1996)
at the Cassegrain focus (f/31) of the UH 2.2{\thinspace}m
telescope in combination with the UH tip-tilt system
(\markcite{Jim98}Jim et al.\ 1998).
The cross-dispersed echelle design of KSPEC
provided simultaneous coverage of the entire
1--2.5 \micron{} wavelength region.
The projected pixel size of the HAWAII 1024 $\times$ 1024 array
was 0\farcs{}167 along the slit and $\simeq$ 5.6 \AA{} at
2 \micron{} along the dispersion direction.  We used a 0\farcs{}96
wide slit oriented East-West and centered on the intensity peak
of component B (Patnaik et al.\ 1992) of B{\thinspace}1422+231 (see Figure 1).
However, comparing the two images given in Figure 1, we find that
actual light from component B is 59{\thinspace}\% of the total measured flux
and the remaining 41{\thinspace}\% is associated with light contamination
from components A and C.
Therefore, the flux of our spectra is 1.7 times as bright as real component B
flux.

Thirty exposures, each of 180 sec
integration, taken under photometric conditions, were obtained by
shifting the position of the object along the slit
at intervals of 5\arcsec{} between each integration.
The total integration time was 5400 sec.
An A-type standard star, HD{\thinspace}106965
(\markcite{Elias82}Elias et al.\ 1982),
was observed for flux calibration.
Another A-type star, HD{\thinspace}136754
(\markcite{Elias82}Elias et al.\ 1982), was also observed
before and after observing
B{\thinspace}1422+231  to correct for atmospheric absorption.
Spectra of an incandescent lamp and an argon lamp were taken for
flat-fielding and wavelength calibration, respectively.
Typical widths of the spatial profiles of the
standard star spectra were $\sim${\thinspace}0\farcs{}5 (FWHM) throughout the night.
Note that our previous NIR spectroscopy of this quasar at KPNO
was obtained using a 1.44 arcsec slit
under 1.4 -- 2.1 arcsec seeing conditions (K96).

\placefigure{fig-1}

Data reduction was performed with IRAF\footnote{%
Image Reduction and Analysis Facility (IRAF) is
distributed by the National Optical Astronomy Observatories,
which are operated by the Association of Universities for Research
in Astronomy, Inc., under cooperative agreement with the National
Science Foundation.} using standard procedures as outlined in
\markcite{Hora96}Hora \& Hodapp (1996).
Sky and dark counts were removed by subtracting the average of
the preceding and following exposures, and then the resulting
frame was divided by a normalized dome flat.
The target quasar was not bright enough to trace
its position in each frame with sufficient accuracy.
Therefore, we first fit the spectral positions of the standard star
spectra with a third-order polynomial function
which properly traces the echelle spectrum.
These fitting results were then applied to the quasar spectra.
Using this procedure, we extracted the quasar spectra with an aperture
of 3\arcsec{} which was determined to be the typical width
where the flux was $\sim${\thinspace}10{\thinspace}\% of
the peak flux along the spatial profile of the standard star.
In order to subtract the sky emission, we used the data just adjacent
to the 3\arcsec{} aperture.
The wavelength scale of each extracted spectrum was calibrated to an
accuracy of 18 km s${}^{-1}$ at 1.1 \micron, 20 km s${}^{-1}$ at
1.6 \micron, and 19 km s${}^{-1}$ at 2.0 \micron,
based on both the argon emission lines of the calibration lamp and
on the telluric OH emission lines.
The spectral resolutions (FWHM) measured from the argon
lamp spectra were $\simeq$ 500 km s${}^{-1}$ at 1.1 \micron, $\simeq$
450 km s${}^{-1}$ at 1.6 \micron, and $\simeq$ 500 km s${}^{-1}$ at 2.0 \micron.
Finally, the spectra were median combined in each band.
Atmospheric absorption features were removed
using the spectra of the A-type star HD{\thinspace}136754
because A-type stars are best suited for correcting for atmospheric absorption
features.
However, since A-type stars inherently have hydrogen recombination
absorption lines (e.g., the Brackett series in $H$ and $K$ bands
and the Paschen series in $I$ and $J$ bands),
we removed these features using Voigt profile fitting before the correction.
In order to check whether this procedure worked
we also applied the same atmospheric correction for the spectra of M-type stars
whose data were obtained on the same night.
Comparing our corrected spectra of the M-type stars with their
published spectra (Lan\c{c}on \& Rocca-Volmerange 1992),
we found that our correction procedure works appropriately.
Finally, in order to calibrate the flux scale, we used
the spectrum of the standard star HD{\thinspace}106965 (A2, $K$=7.315)
divided by a 9000{\thinspace}K blackbody spectrum, which fits
the $J\!H\!K\!L$ magnitude of the standard
\markcite{Elias82} (Elias et al.\ 1982)
with only 1.2{\thinspace}\% deviation.
Photometric errors were determined to be $<${\thinspace}10{\thinspace}\%
over all observed wavelengths.

\section{RESULTS AND DISCUSSION}

Figure \ref{fig-2} shows the spectra of B{\thinspace}1422+231 (red line)
in the $I\!H\!K$ bands
together with both our previous measurement at KPNO
(blue line; K96) and
the Large Bright Quasar Survey (LBQS) composite
spectrum shifted to $z=3.62$ (black dashed line; \markcite{Francis91}
Francis et al.\ 1991).
Since the efficiency of KSPEC in the $J$-band is not high,
we have not used the $J$-band data in this paper.

\placefigure{fig-2}

Although our new measurement was made with a narrower slit (0\farcs{}96 wide)
than that used by K96 (1\farcs{}44 wide),
the $H$- and  $K$-band fluxes are slightly higher
than those of the previous measurement.
K96 tried to carefully correct for the effect
of seeing on the relative flux calibration among the spectral bands
because they could take only one band spectrum at a time.
However,
the time variation of the seeing in their observations made it difficult
to perform the correction very accurately.
Since our new  $I$ to $K$-band NIR spectra were obtained simultaneously,
our new measurement should be more reliable than
that of K96.
Our new $I$-band spectrum detects \ion{C}{3}]$\lambda$1909 emission clearly.
Further, some unidentified emission features at 2000{\thinspace}\AA~
and 2080{\thinspace}\AA ~
as well as the dip at 2200{\thinspace}\AA~ in the rest frame, which are
shown in the average spectrum of LBQS quasars (Francis et al.\ 1991),
are also seen.
The presence of [\ion{O}{3}]$\lambda$5007 emission has also been confirmed.

The continuum flux of component B between 1330{\thinspace}\AA~ and
1380{\thinspace}\AA, which was obtained by Impey et al.\ (1996)
13 months before our observations,
was assumed for the continuum level on the short
wavelength side of our spectra.
We applied the factor 1.7 to this continuum level
in order to correct the contaminated light from component A and C for
our spectra (see previous section).
The continuum level on the long wavelength side of our spectra
was chosen by fitting a power-law
continuum plus a \ion{Fe}{2} template simultaneously to minimize the residual
between 4450{\thinspace}\AA~ and 4750{\thinspace}\AA~ and beyond 5100{\thinspace}\AA.  We used the Balmer continuum
template that was generated to approximate the emission-line strengths
seen in 0742+318 (see Figure 3d of Wills et al.\ 1985).  The \ion{Fe}{2}
emission profile of the very strong \ion{Fe}{2}-emitting low-ionization
BAL quasar PG{\thinspace}0043+039 (Turnshek et al.\ 1994)
was used as the \ion{Fe}{2} template.
Due to a relatively poor fit over the entire wavelength range
when using a single
template, two \ion{Fe}{2} templates were used: one for shortward of
3000{\thinspace}\AA \ and the other for longward of 3000{\thinspace}\AA.
The power-law continuum that we derive is given by $f_\nu \propto \nu^{-0.88}$.
This spectral index is steeper than the
$f_\nu \propto \nu^{-0.54}$ power-law derived in K96.
Yee \& Bechtold (1996) report that B{\thinspace}1422+231 had become brighter by
0.12 mag during 13 months (see also Kundi\'c et al.\ 1997),
but variability of 0.12 mag would change
the spectral index to $-0.88 \pm 0.08$ at most.
Therefore, it is unlikely that
the difference of the spectral index between K96 and our current work
is due to the variability inherent in this quasar.
Thus, we conclude that the
change in the calculated spectral index is due to
our improved absolute NIR spectrophotometry
by simultaneous observations of $I\!H\!K$ bands using KSPEC.

Figure 3 shows the rest frame spectrum of B{\thinspace}1422+231 with the power-law
subtracted and the best fit synthetic spectrum comprised of \ion{Fe}{2} and Balmer
continuum emission as well as other broad-lines.
The fluxes, equivalent widths, and line widths
of the detected emission lines are summarized in Table 1
together with the previous measurements of K96.
Our new data show that the ratio \ion{Fe}{2}(UV)/H$\beta$ and
\ion{Fe}{2}(optical: $\lambda\lambda$3500--6000)/H$\beta$ are
higher than those of K96
by  factors of 1.6 and 3.3, respectively.
These differences are mainly due to adopting the different power-law continuum
described above.
However,  our new \ion{Fe}{2}(optical)/H$\beta$ ratio
for B{\thinspace}1422+231 is still in the
observed range for low-$z$, radio-loud quasars
(Boroson \& Green 1992; see Taniguchi et al.\ 1997a; Murayama et al.\ 1998).
The [\ion{O}{3}]$\lambda$5007/H$\beta$ ratio is nearly the same
between the two measurements.
We also note that there appears to be evidence for excess emission
at $\lambda <$ 2100{\thinspace}\AA ~ in the rest frame
(see the second panel of Figure 3).

Our new measurements yield a flux ratio of
\ion{Fe}{2}(UV)/\ion{Fe}{2}(optical) $\simeq$ 8.0.
Since Wills et al.\ (1985) give a range of
4 $<$  \ion{Fe}{2}(UV)/\ion{Fe}{2}(optical) $<$ 12
for low-$z$ quasars, B{\thinspace}1422+231 is typical of low-$z$ quasars
in this respect.
We also obtain a flux ratio of
\ion{Fe}{2}(optical: $\lambda\lambda$4484--4684)/H$\beta
\simeq 0.53$ which is significantly smaller than the range of values
($\sim${\thinspace}1.5 -- 2) found by Hill et al.\ (1993) for quasars
with $z \simeq$ 2 -- 2.5{\thinspace}.
This demonstrates that it is perhaps dangerous to attempt a measurement of
the iron abundance using solely
optical \ion{Fe}{2} emission features as suggested by Wills et al.\ (1985).

Finally we speculate about the formation epoch of the host galaxy of
B{\thinspace}1422+231. Our new measurement has confirmed that
the \ion{Fe}{2}(total)/H$\beta$
ratio for B{\thinspace}1422+231 is higher-than-normal with respect to
what is found for low-$z$ quasars.
Since the low-$z$ quasars are believed to be associated with nuclei
of massive galaxies, their chemical abundances are expected to be higher than
or roughly equal to solar.
Therefore, the observed higher \ion{Fe}{2}(total)/H$\beta$ ratio suggests that
the iron abundance of B{\thinspace}1422+231 is at least comparable to the solar value.
If this is the case, it is expected that the majority of iron would come from
Type Ia supernovae.
Yoshii, Tsujimoto, \& Nomoto (1996) derived $\sim$\ 1.5 Gyr
for the lifetime of SN Ia progenitors from an analysis of
the O/Fe and Fe/H abundances in solar neighborhood stars.
If the Fe enrichment started at 1.5 Gyr after
the onset of the first epoch of star formation,
the host galaxy of B{\thinspace}1422+231 would
have formed at $z \sim 9$ or earlier for $q_{0}$ = 0.0 and
$H_{0}$ = 75 km s$^{-1}$ Mpc$^{-1}$ (K96).
\placefigure{fig-3}

\begin{table}
\dummytable\label{tbl-1}
\end{table}

\placetable{tbl-1}

\acknowledgments

We are very grateful to the staff of the UH 2.2 m telescope.
In particular, we would like to thank
Andrew Pickles for his technical support and assistance with the observations.
This work was financially supported in part by Grants-in-Aid for
Scientific Research (Nos.\ 07044054 and 09640311)
from the Japanese Ministry of Education, Science, Sports, and Culture
and by the Foundation for Promotion of Astronomy, Japan.
TM is thankful for support from a Research Fellowship from the Japan
Society for the Promotion of Science for Young Scientists.
This research has made use of the NASA/IPAC Extragalactic Database
(NED) and the NASA Astrophysics Data System Abstract Service.