Super-Solar Metallicity in Weak Mg II Absorption Systems at z ~ 1.7

Abstract: Through photoionization modeling, constraints on the physical conditions of three z ~ 1.7 single-cloud weak Mg II systems (W_r(2796) < 0.3A) are derived. Constraints are provided by high resolution R = 45,000, high signal-to-noise spectra of the three quasars HE0141-3932, HE0429-4091, and HE2243-6031 which we have obtained from the ESO archive of VLT/UVES. Results are as follows:(1) The single-cloud weak Mg II absorption in the three z ~ 1.7 systems is produced by clouds with ionization parameters of -3.8 < logU < -2.0 and sizes of 1-100 pc.(2) In addition to the low-ionization phase Mg II clouds, all systems need an additional 1-3 high-ionization phase C IV clouds within 100 km/s of the Mg II component. The ionization parameters of the C IV phases range from -1.9 < logU < -1.0, with sizes of tens of parsecs to kiloparsecs.(3) Two of the three single-cloud weak Mg II absorbers have near-solar or super-solar metallicities, if we assume a solar abundance pattern. Although such large metallicities have been found for z < 1 weak Mg II absorbers, these are the first high metallicities derived for such systems at higher redshifts.(4) Two of the three weak Mg II systems also need additional low-metallicity, broad Lya absorption lines, offset in velocity from the metal-line absorption, in order to reproduce the full Lya profile.(5) Metallicity in single-cloud weak Mg II systems are more than an order of magnitude larger than those in Damped Lya systems at z ~ 1.7. In fact, there appears to be a gradual decrease in metallicity with increasing N(HI), from these, the most metal-rich Lya forest clouds, to Lyman limit systems, to sub-DLAs, and finally to the DLAs.

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

Many single-cloud weak \ion{Mg}{2} absorbers ($W_r$(2796) $\leq$
0.3\AA) at 0.4 $<$ $z$ $<$ 1.0 have been found to have solar or even
super-solar metallicities \citep{rig02,cha03}. This is quite puzzling
in view of the fact that they tend to be found at distances of 40 --
100~kpc from star-forming galaxies \citep{chu05,mil06}, and not near
detected sites of star formation.  These objects have a cross-section
similar to the absorption cross section of luminous galaxies, thus
they occupy a significant volume of the universe \citep{rig02}.
Because weak \ion{Mg}{2} systems only contain a small mass, they do
not trace the cosmic metal abundance, however, it is still important
to understand local star-formation activities in the full range of
environments over cosmic time.

\citet{nar07} presented the results of a search for weak \ion{Mg}{2}
absorption in VLT/UVES quasar spectra over a cumulative redshift path
of 77.3 at 0.4 $<$ $z$ $<$ 2.4.  A total of 116 \ion{Mg}{2} absorbers
with 0.02 $<$ $W_r(2796)$ $<$ 0.3~\AA\ were detected, with $\sim$ 60\%
of these having just a single-cloud, narrow component in \ion{Mg}{2}.
It has been argued that a large fraction of the population of
multiple-cloud weak \ion{Mg}{2} absorbers (i.e., absorbers with
multiple-clouds, which are resolved by high-resolutional ($R$ $\sim$
40,000) spectroscopies) have a different physical origin from the
single-cloud weak \ion{Mg}{2} absorbers \citep{rig02,mas05,din05}.
\citet{nar07} found that the redshift path density, $dN/dz$, of weak
\ion{Mg}{2} absorbers increases from $z$ $=$ 2.4 to $z$ $=$ 1.2, where
it peaks, and subsequently declines until the present.

The paper of \citet{nar07} followed a weak \ion{Mg}{2} survey of
\citet{lyn06}, which used a small subset of high quality VLT/UVES
quasar spectra.  From that survey, \citet{lyn07} selected the nine $z$
$>$ 1.4 weak \ion{Mg}{2} absorbers and applied photo-ionization
models.  Of these, four were single-cloud weak \ion{Mg}{2} absorbers
with coverage of the corresponding Ly$\alpha$\ transitions so that
metallicities could be constrained.

In general, \citet{lyn07} found that the physical conditions of the
$z$ $>$ 1.4 weak \ion{Mg}{2} absorbers were similar to those of the
same class of absorbers at lower redshift, suggesting that the same
production mechanisms are responsible over this time period.  However,
in contrast to the situation at low redshift \citep{rig02,cha03},
\citet{lyn07} were only able to derive lower limits to the
metallicity, ranging from 1/100$^{th}$ -- 1/10$^{th}$ of the solar
value.  Single-cloud weak \ion{Mg}{2} absorbers, in addition to the
cloud that produced the \ion{Mg}{2} absorption, have separate, lower
density \ion{C}{4} clouds that are spread over tens to hundreds of
{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ \citep{cha03,lyn07}.  The reason that only lower limits could be
derived for the metallicity in these cases was because of the
contributions to the Ly$\alpha$\ profile from the kinematically spread
\ion{C}{4} clouds.  On the other hand, \citet{lyn07} also suggested
that there could be a build-up of metals in the population of
single-cloud weak \ion{Mg}{2} absorbers from $z$ $\sim$ 2 to $z$
$\sim$ 0.4.

In view of the large spread of metallicity constraints at any
redshift, a sample of only four $z$ $>$ 1.4 single-cloud weak
\ion{Mg}{2} absorbers was not sufficient to resolve the issue of
whether metallicity evolution occurs in the population.  We therefore
identified three additional $z$ $>$ 1.4 candidates for
photo-ionization modeling from the larger survey of \citet{nar07}.
These systems, the $z$ $=$ 1.78169 absorber toward HE~0141$-$3932, the
$z$ $=$ 1.68079 absorber toward HE~0429$-$4091, and the $z$ $=$
1.75570 absorber toward HE~2243$-$6031, have high $S/N$-ratio VLT/UVES
coverage of key constraining metal-line transitions, as well as Ly$\alpha$.
We apply Cloudy \citep{fer98} photoionization models to constrain the
ionization parameters and metallicities of the phases\footnote{Phases
are defined as regions of gas with similar metallicity, temperature,
and volume density.} of gas that are required to fit the data.  Our
focus is on the question of the metallicity evolution of the
single-cloud weak \ion{Mg}{2} absorbers, and on comparing the result
to the metallicity evolution of other classes of quasar absorption
line systems.

We will begin in \S~\ref{sec:data} with a summary of the data that we
have used to constrain the properties of the three $z$ $\sim$ 1.7
single-cloud weak \ion{Mg}{2} absorbers.  Next, in \S~\ref{sec:photo},
we describe our procedure for photo-ionization modeling with Cloudy,
which is based on comparing synthetic model profiles to the data.
\S~\ref{sec:systems} presents the detailed results for the three
absorbers.  Finally, in \S~\ref{sec:discussion}, we summarize those
results, compare the metallicities of single-cloud weak \ion{Mg}{2}
absorbers to DLAs, to sub-DLAs, and to Lyman limit absorbers, and
discuss the implications of this comparison.  Throughout this paper,
we use a cosmology with $H_{0}$=72 {\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi} Mpc$^{-1}$, $\Omega_{m}$=0.3,
and $\Omega_{\Lambda}$=0.7.

\section{Data}
\label{sec:data}

All three weak \ion{Mg}{2} systems in our sample were found in the
survey of \citet{nar07} who studied the statistical properties of weak
\ion{Mg}{2} systems, using 81 VLT/UVES spectra.  Data were retrieved
from the ESO archive, and reduced using the ESO MIDAS pipeline. In the
case of multiple epoch observations, the individual exposures were
combined to enhance the $S/N$ ratio.  The detailed description of the
original data analysis is presented in \citet{nar07}.

The main purpose of this paper is to place strict constraints on
physical properties (such as metallicity, ionization
parameter\footnote{The ionization parameter is defined the ratio of
ionizing photons ($n_{\gamma}$) to the number density of hydrogen in
the absorbing gas ($n_{H}$).}, gas temperature, and size) of single
cloud, weak \ion{Mg}{2} absorbers at high redshift. Among the 116 weak
\ion{Mg}{2} systems in \citet{nar07}, 16 systems at relatively high
redshift ($z$ $>$ 1.5) have Ly$\alpha$\ absorption lines covered with the
observed spectra. Ly$\alpha$\ is quite important for determining
metallicities. Of the 16 systems, 7 are classified as single-cloud
systems. Four of these seven, the $z$ = 1.65146 system toward
HE0001$-$2340, the $z$ = 1.70849 system toward HE0151$-$4326, the $z$
= 1.79624 system toward HE2347$-$4342, and the $z$ = 2.17455 system
toward HE0940$-$1050 had already been modeled as described in
\citet{lyn07}, since they were first identified in the earlier survey
of \citet{lyn06}.  The other three systems, the $z$ = 1.78169 system
toward HE~0141$-$3932, the $z$ = 1.68079 system toward HE~0429$-$4091,
and the $z$ = 1.75570 system toward HE~2243$-$6031, are modeled in the
present paper.

In Table~\ref{tab1}, we give an observation log for the three quasars
toward which the single-cloud weak \ion{Mg}{2} absorbers were
detected.  Column (1) is the quasar name, columns (2) and (3) are the
coordinates of the quasars, columns (4) and (5) are the emission
redshift and optical magnitude of the quasars, column (6) is observed
wavelength range, column (7) is the central wavelength for the blue
and red CCDs of VLT/UVES, columns (8) and (9) are total exposure time
and observing date, and column (10) is the proposal ID.  By chance,
all three quasars were observed as part of programs with the same
P.I., Sebastian Lopez.

The absorption profiles for key transitions used to constrain the
three single-cloud weak \ion{Mg}{2} absorbers are displayed in Figures
\ref{fig1}, \ref{fig2}, and \ref{fig3}.  Rest-frame equivalent widths
and $5\sigma$ rest-frame equivalent width limits for these same
transitions are listed in Table~\ref{tab2}.

\section{Photoionization Models}
\label{sec:photo}

In this section we briefly summarize the strategy for photoionization
modeling, which is similar to the procedure practiced in previous
studies \citep[e.g.,][]{chu99,cha03,lyn07}.  At first, we determine
the line parameters (column density, Doppler parameter, and redshift)
for the \ion{Mg}{2} doublets, using a Voigt profile fitting code ({\sc
MINFIT}; \citealt{chu03}).  For each of the three absorbers modeled
here, a single absorption component provided an adequate fit to the
\ion{Mg}{2} doublet.  The best fit parameters are listed in
Table~\ref{tab3}: column (1) is the name of the transition, columns
(2), (3), and (4) are the fit parameters of velocity shift from the
system central redshift\footnote{The system central redshift is
defined as the flux-weighted center of the \ion{Mg}{2}~$\lambda$2796
absorption line.}, column density, and Doppler parameter with their
1$\sigma$ errors.  For completeness, we also applied formal Voigt
profile fits to other detected transitions, and listed the results in
Table~\ref{tab3}.  We emphasize, however, that we did not use these
values (except for the optimized transitions) to determine model
constraints.  Instead we compared directly to the shapes of the
observed absorption profiles as described below.

Assuming that the absorbing clouds are in photoionization equilibrium,
their ionization conditions are derived using the photoionization code
Cloudy, version 07.02.00 \citep{fer98}, optimizing on the observed
column density of \ion{Mg}{2} and comparing model predictions to other
observed low ionization transitions. The clouds in each phase are
modeled as plane-parallel structures of constant density.  The
ionization parameter ($\log U$ $=$ $\log[n_{\gamma}/n_H]$) and
metallicity ($Z/Z_{\odot}$) of the clouds serve as free parameters in
the modeling. The elemental abundance pattern is initially assumed to
be solar, and variations are explored when they are suggested by the
data.

Following \citet{haa96,haa01}, we take the combined flux from quasars
and star forming galaxies at $z \sim 1.7$ (with a photon escape
fraction of 0.1) as an extragalactic background radiation (i.e.,
incident flux on the absorbers). We will consider the possible effects
of a nearby stellar contribution to the ionizing radiation field in
\S~\ref{sec:assumptions}, although it should be safe to ignore them,
since weak \ion{Mg}{2} absorbers are not often within 40~kpc of bright
galaxies \citep{rig02,chu05}.

The metallicity and ionization parameter are initially varied in steps
of 0.5 dex, and then fine-tuned to select the model that corresponds
best with the observed data in steps of 0.1 dex. The column densities
for various ionization stages of each element and the equilibrium gas
temperature are provided by the Cloudy model.  The Doppler parameter
for each element can then be calculated from the expression for total
line width, $b_{tot}^{2}$ $=$ $b_{ther}^{2}$ $+$ $b_{turb}^{2}$, where
$b_{ther}$ ($=$ $\sqrt{2kT/m}$) is the thermal contribution to the
line width corresponding to a gas temperature $T$ and $b_{turb}$, the
contribution from internal gas turbulence. This latter quantity, which
is uniform across elements, is estimated using the gas temperature
from the Cloudy model and the observed Doppler parameter of the
element on which the model is optimized. From the derived column
densities and Doppler parameters, a synthetic spectrum is generated
and compared to the observed spectrum after convolving with a Gaussian
instrumental spread function ($R$ $=$ 45,000 in the case of our
VLT/UVES spectra).

Although the observed and the synthesized spectra are compared ``by
eye'', the method generally gives good results in the sense that other
models whose metallicity and ionization parameter are (sometimes only
slightly) different than the best values would make the model spectrum
deviate greatly from the observed spectrum.  It is not practical to
apply formal procedures to assess the quality of fit because of the
effect of blends and data defects on a $\chi^2$ statistic.  If all
transitions of a certain element are over/under-produced compared to
the observed spectrum in the best model, we consider a specific
abundance pattern for the system.  Our method is illustrated for a
specific example, the $v=1$~{{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}} cloud for the $z=1.78169$ system
toward HE0141-3932 (System 1) in Figure~\ref{fig4}.  We show, for this
component, that only values within the range $-0.7 < \log Z < 0.1$ and
within 0.1~dex of $\log U = -2.3$ produce acceptable fits to Ly$\alpha$,
\ion{Si}{2}, and \ion{Si}{4}, which are the main constraints for this
case.

Usually, low-ionization phase clouds that produce the \ion{Mg}{2}
absorption lines also produce other low-ionization transitions such as
\ion{Fe}{2}, \ion{Si}{2}, \ion{C}{2}, and \ion{Al}{2}. However, the
higher-ionization transition, \ion{C}{4}, and sometimes
intermediate-ionization transitions such as \ion{Si}{4} are not fully
produced by the \ion{Mg}{2} clouds. Therefore, once we find the best
model for the \ion{Mg}{2} clouds, we repeat a similar Cloudy analysis
for the \ion{C}{4} clouds. For the high-ionization phase, multiple
components are sometimes required to fit \ion{C}{4} absorption lines
detected at different velocities.  The results of Voigt profile fits
to these \ion{C}{4} components are also listed in Table~\ref{tab3}.
Finally, the combined models, with contributions from both low- and
high-ionization phases, are compared to the observed spectrum. Since
Ly$\alpha$\ absorption lines have contributions from both the low- and
high-ionization phases, determination of metallicities for the low-
and high-ionization phases can be degenerate.  However, in the three
single-cloud weak \ion{Mg}{2} absorbers that we are considering, it is
possible to place limits on metallicity of both phases, so as not to
over-produce the Ly$\alpha$\ absorption.  Also, in some cases, models with
high metallicity (leading to low temperatures) can be favored when
comparing model profiles to some observed low ionization transitions.

\section{Individual Systems}
\label{sec:systems}

For photoionization models, our procedure is adjusted for each system
based on the available constraints in various transitions.  Some
transitions are blended with other lines, and others may be
saturated. These issues are considered carefully during
photoionization modeling. In this section, we present the spectra of
our three weak \ion{Mg}{2} systems and then describe the results of
photoionization models. We can place strong constraints on metallicity
and ionization parameter because multiple transitions (including Ly$\alpha$)
have been detected in high S/N-ratio regions of the spectra (see
Figures~\ref{fig1} -- \ref{fig3}). We will, hereafter, refer to the
three weak \ion{Mg}{2} systems as Systems~1, 2, and 3.  After
presenting the constraints, we discuss in \S~\ref{sec:assumptions} the
effect of varying model parameters and assumptions such as
photoionization equilibrium, the shape of the incident radiation
field, and the abundance of dust.

\subsection{System~1 (HE~0141$-$3932; $z$=1.78169)}
\label{sec:sys1}

This system is detected in Ly$\alpha$, \ion{Mg}{2}, \ion{Si}{2},
\ion{Al}{2}, \ion{C}{2}, \ion{Al}{3}, \ion{Si}{3}, \ion{Si}{4}, and
\ion{C}{4}, as shown in Figure~\ref{fig1}.  The observed spectrum
covers \ion{Mg}{1}, \ion{O}{1}, and \ion{N}{5}, but these transitions
are not detected to a 5$\sigma$ detection limit.  Although
\ion{Fe}{2}~$\lambda$2383 is formally detected at only 3$\sigma$ and
may be contaminated by spurious features, \ion{Fe}{2}~$\lambda$2600 is
also detected at 2.5$\sigma$ and both are precisely aligned with the
\ion{Mg}{2}.  We therefore are confident in applying the \ion{Fe}{2}
as a constraint.  Rest-frame equivalent widths and $5\sigma$
equivalent width limits are listed in Table~\ref{tab2}. Although the
Ly$\alpha$\ profile is complex, it provides a very strong constraint on
metallicity because the flux at the position of the component centered
on \ion{Mg}{2} begins to recover on its blue side, which enables us to
fit this component of Ly$\alpha$\ effectively. All transitions except for
high-ionization lines such as \ion{C}{4} and \ion{Si}{4} have only
single components detected.

We begin by fitting a Voigt profile to the \ion{Mg}{2} doublet.  The
best fit parameters are listed in Table~\ref{tab3}.  A lower limit to
the ionization parameter is provided by \ion{Fe}{2}, which is
over-produced at $\log U$ $<$ $-$3.8.  \ion{Mg}{1} is also
over-produced at $\log U$ $<$ $-$6.0.  (\ion{Mg}{1} was not detected
in the observed spectrum, and \ion{Fe}{2} is only detected in the
$\lambda$2383 transition at a level of $\sim$ 3 $\sigma$.)  There is
also a strict upper limit on the ionization parameter of $\log U$ $<$
$-$2.8 above which the high-ionization transitions \ion{C}{4} and
\ion{Si}{4} are over-produced.  Based on these considerations, models
with $\log U$ $=$ $-$3.8 -- $-$2.8 seem acceptable, however
\ion{Si}{2} and \ion{Al}{2} provide additional constraints.  An
ionization parameter, $-$3.8 $<$ $\log U$ $<$ $-$3.7, is consistent
with all constraints.  At lower values of $\log U$ \ion{Fe}{2} is
over-produced, and at higher values \ion{Al}{2} and \ion{Si}{2} are
over-produced.  Independent of ionization parameter, the metallicity
must be $\log Z$ $>$ $-$0.7 in order that Ly$\alpha$\ is not over-produced.
An upper limit of $\log Z$ $<$ $-$0.5 would also apply, for solar
abundance pattern, if $\log U$ $=$ $-$3.7, in order that \ion{Fe}{2}
absorption is not over-produced.  The lower limit on metallicity of
$\log Z$ $>$ $-$0.7 applies for the solar abundance pattern, and would
be adjusted downwards for an $\alpha$-enhanced pattern.  Similarly,
$\log U$ could be lower in the $\alpha$-enhanced case.  If dust
depletion was important, the lower limit on metallicity would be
increased.

In addition to the \ion{Mg}{2} cloud, three more higher-ionization
clouds are needed to reproduce the three detected \ion{C}{4}
components, the two \ion{Si}{4} components, and the right wing of the
Ly$\alpha$\ profile. There are no \ion{N}{5} lines detected, which is useful
to place an upper limit on the ionization parameter.  The results of
Voigt profile fits to the \ion{C}{4} doublet are given in
Table~\ref{tab3}.  Beginning with these \ion{C}{4} fits, we adjust the
ionization parameter and metallicity to satisfy constraints from other
transitions.  In order to avoid over-production of \ion{Si}{4} and low
ionization transitions (at low values) or \ion{N}{5} (at high values),
the ionization parameters for the additional 3 clouds at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ $=$
$-$77, 1, and 31~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ from the system center should be $-$1.5 $<$
$\log U$ $<$ $-$1.0, $-$2.35 $<$ $\log U$ $<$ $-$2.25, and $-$1.9 $<$
$\log U$ $<$ $-$1.8, respectively.  For these ionization parameters,
lower limits on metallicity of $\log Z$ $>$ $-1.8$, $\log Z$ $>$
$-0.7$, and $\log Z$ $>$ $-0.5$ will apply for the three \ion{C}{4}
components in order that Ly$\alpha$\ is not over-produced.  For the first
two components, at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ $=$ $-$77 and 1~{{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}} upper limits on $\log
Z$ of 0.5 and 0.1 in order that low ionization transitions are not
over-produced.  The third \ion{C}{4} component, at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ $=$ 31~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}
, has a strong constraint on its metallicity, assuming that it is
responsible for the Ly$\alpha$\ absorption in the red wing of the profile
(since the \ion{Mg}{2} cloud cannot account for that absorption
without over-producing Ly$\alpha$\ in the blue wing).  For that \ion{C}{4}
we find that $-$0.5 \lower0.3em\hbox{$\,\buildrel <\over\sim\,$}\ $\log Z$ \lower0.3em\hbox{$\,\buildrel <\over\sim\,$}\ $-$0.4 in order to match
the Ly$\alpha$\ profile.  We emphasize that the three components have
different ionization parameters from each other, although these values
are all higher than that of the \ion{Mg}{2} cloud.

The one \ion{Mg}{2} and three \ion{C}{4} clouds reproduced all
transitions except for Ly$\alpha$. There are absorption features at both
sides of the Ly$\alpha$\ profile, at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ $<$ $-$50 and at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ $>$
70~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ (see Figure~\ref{fig1}). Neither of them are Ly$\beta$\ profiles
from higher redshift systems, because the corresponding Ly$\alpha$\
absorption lines would be located redward of the quasar Ly$\alpha$\ emission
line.  The \ion{C}{4} component at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ = $-$77~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ cannot itself
give rise to all of the observed Ly$\alpha$\ around that velocity because it
is too narrow.  Thus, we fit the broad Ly$\alpha$\ profile ($b(H)$ $=$
42~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ as listed in Table~\ref{tab3}), and placed upper limits on
its metallicity: ($\log Z$ $<$ 0.0 for $\log U$ $\sim$ $-$3 and $\log
Z$ $<$ $-$2.0 for $\log U$ $\sim$ $-$1.5) so as not to produce other
transitions at this velocity.  Ionization parameter is not constrained
for this broad Ly$\alpha$\ cloud, although high values of $\log U$ are ruled
out unless the metallicity is extremely small.

We conclude that, for a solar abundance pattern, the metallicity of
the $z$ $=$ 1.78169 weak \ion{Mg}{2} system in the HE~0141-3932
spectrum is constrained to be $-$0.7 \lower0.3em\hbox{$\,\buildrel <\over\sim\,$}\ $\log Z$ \lower0.3em\hbox{$\,\buildrel <\over\sim\,$}\ $-$0.5.  A
strict lower limit on metallicity of $\log Z$ $>$ $-$0.7 applies in
order that Ly$\alpha$\ absorption is not over-produced.  Two of the three
\ion{C}{4} clouds related to this system are constrained to have
similar or higher metallicities, and for the other no significant
metallicity constraint is available.  The \ion{Mg}{2} cloud has a
sub-parsec size, much smaller than the three \ion{C}{4} phase
absorbers, with sizes of $>$10~pc. All components are optically thin,
the \ion{Mg}{2} cloud giving rise to $\log N_{\rm HI}$\ $\sim$ 15.9.

A summary of how specific transitions were used to constrain $\log U$
and $\log Z$ for this system is given in Table~\ref{tab4}.  Ranges of
acceptable model parameters for each model cloud are listed in
Table~\ref{tab5}, in which column (1) is the absorption redshift,
column (2) is the optimizing transition used in Cloudy, column (3) is
the velocity shift from the system center, and column (4) is the
acceptable range of ionization parameter.  Column (5) of the table
presents the best constraint on metallicity, considering both
comparison to the observed Ly$\alpha$\ and the fit to metal-line
transitions, assuming a solar abundance pattern.  Column (6) is a more
conservative constraint on the metallicity of the cloud, using only
the requirement that Ly$\alpha$\ absorption is not over-produced.  Column
(7) is the depth of absorber along the line of sight, column (8) is
the Doppler parameter of the optimizing transition (listed in column
(2)), and column (9) is the column density of neutral hydrogen gas in
the absorbers.  An example of a best fit model and its specific
parameters is presented in Figure~\ref{fig1} and Table~\ref{tab6}.

\subsection{System~2 (HE~0429$-$4091; $z$=1.68079)}
\label{sec:sys2}

This is a very weak \ion{Mg}{2} system whose rest-frame equivalent
width ($W_{r}$(2796) = 0.015\AA) is just above the lower limit of the
survey presented in \citet{nar07}. In this system, Ly$\alpha$, \ion{Si}{2},
\ion{Al}{2}, \ion{C}{2}, \ion{Al}{3}, \ion{Si}{3}, \ion{Si}{4}, and
\ion{C}{4} are clearly detected, but \ion{O}{1}, \ion{Fe}{2}, and
\ion{N}{5} are not detected at 5$\sigma$, as shown in
Figure~\ref{fig2}.  The rest-frame equivalent widths and $5\sigma$
limits are tabulated in Table~\ref{tab2}.

Again, we begin by fitting Voigt profiles to the \ion{Mg}{2} doublet
(see Table~\ref{tab3}).  Low ionization parameters, smaller than $\log
U$ $=$ $-$3.5, are not acceptable because of over-production of
\ion{Fe}{2}. At $\log U$ $>$ $-$2.0, \ion{C}{4} is over-produced.  For
this system, the Ly$\alpha$\ absorption is strong, and it probably arises
from offset \ion{C}{4} components rather than from the \ion{Mg}{2}
cloud.  Because of the strong Ly$\alpha$\ absorption we cannot place a
meaningful constraint on the metallicity of the system; the blue wing
of Ly$\alpha$\ only requires $\log Z$ $>$ $-$3.0.  However, at $-$3.5 $<$
$\log U$ $<$ $-$2.0, \ion{Al}{2} is over-produced at low metallicity
(such that $\log Z$ $>$ $-$1.0 for $\log U \sim$ $-$3.0 and $\log Z$
$>$ $-$0.3 for $\log U \sim$ $-$2.0).  Because the acceptable range of
ionization parameter is large for this system, we consider two extreme
cases: (i) the \ion{Mg}{2} cloud has the highest acceptable ionization
parameter ($\log U$ $=$ $-$2.0) that produces the maximum possible
\ion{C}{4} and \ion{Si}{4} absorption, and (ii) the \ion{Mg}{2} cloud
has the lowest acceptable ionization parameter ($\log U$ $=$ $-$3.5)
so that it does not significantly contribute to the \ion{C}{4} and
\ion{Si}{4} profiles. We will refer to these below as the maximal and
minimal \ion{C}{4} models.

In the maximal \ion{C}{4} model, we need an additional three
high-ionization components to fit \ion{C}{4} profile at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ $=$
$-$47, $-$28, and 5~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ from the system center. The fitting
parameters, after removing the contributions to \ion{C}{4} from the
low-ionization clouds, are listed in Table~\ref{tab3}.
Figure~\ref{fig2} also shows that the \ion{Mg}{2} cloud at 0~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ is
offset and is not broad enough to explain all of the \ion{C}{4}
absorption at a similar velocity.  Only the \ion{C}{4} cloud at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\
$=$ $-$47~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ has its metallicity constrained to be $\log Z$ $=$
$-$1.1 -- $-$1.0, assuming that it accounts for the blueward edge of
the saturated Ly$\alpha$\ profile.  The metallicity cannot be lower than
this value or Ly$\alpha$\ would be overproduced.  However, a constraint on
the ionization parameters of the three \ion{C}{4} clouds can be
obtained by the requirement that they also fit the \ion{Si}{3} and
\ion{Si}{4} doublet profile.  For the three \ion{C}{4} clouds, we find
$\log U$ in the range between $-$1.8 and $-$1.6.  The only
disagreement between this maximal \ion{C}{4} model and the observed
spectrum is an over-production of \ion{N}{5}. Adjusting the abundance
pattern so that nitrogen is reduced by at least 0.5 dex compared to
the other elements gives adequate results (see Figure~\ref{fig2}).  In
fact, the production mechanisms for nitrogen are poorly understood
\citep[e.g.,][and references therein]{rus92}. Namely, it is not well
understood yet whether massive stars or intermediate-mass stars (or
both) contribute to the {\it primary} nitrogen production before the
{\it secondary} production starts through the CNO cycle
\citep{spi05}. This nitrogen ``problem'' has also been reported
frequently in DLA systems \citep[e.g.,][]{pet02} and for the
Magellanic Bridge \citep[e.g.,][]{leh01}, as well as for
multiple-cloud weak \ion{Mg}{2} systems \citep{zon04}. Once the best
model parameters for the three \ion{C}{4} components are found, we
adjust the parameters for \ion{Mg}{2} again to compensate, and get
$\log U$ $=$ $-$2.0 and $\log Z$ $=$ $-$0.2 -- 0.1. The lower limit on
metallicity was placed in order not to over-produce \ion{Al}{2}, and
the upper limit in order not to under-produce \ion{Si}{4} and
\ion{C}{4}.  Acceptable ranges of model parameters for one \ion{Mg}{2}
and three \ion{C}{4} components (the maximal \ion{C}{4} model) are
summarized in Table~\ref{tab5}.  An example of the best fit model is
overplotted on the observed spectrum in Figure~\ref{fig2}. Here, we
emphasize that the \ion{Mg}{2} component has a relatively high
ionization parameter similar to those of the \ion{C}{4} components.

We also consider a minimal \ion{C}{4} model, for which we need four
more components to reproduce \ion{C}{4} and \ion{Si}{4} because there
is no contribution to the \ion{C}{4} absorption from the \ion{Mg}{2}
phase. We fit the \ion{C}{4} profile with four high-ionization phase
clouds at {\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ of $-$47, $-$31, $-$12, and 3~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ from the system
center. As with the maximal \ion{C}{4} model, the metallicity of the
{\ifmmode{\Delta v}\else{$\Delta v$}\fi}\ = $-$47~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}\ component is constrained to be $\log Z$ $=$
$-$1.0 -- $-$0.9.  We found a range of acceptable ionization
parameters for the first three \ion{C}{4} components, with $\log U$
$=$ $-$1.9 -- $-$1.7.  However, there are no acceptable solutions for
the fourth component (that overlaps with the \ion{Mg}{2}
component). The best model, with $\log Z$ $=$ $-$0.5 and $\log U$ $=$
$-$1.9, reproduces observations for most all transitions, but it gives
rise to additional \ion{Mg}{2} absorption.  These double contributions
from low- and high-ionization phase clouds overproduce the observed
\ion{Mg}{2} profile. This result suggests that the \ion{Mg}{2}
absorption is partially produced by the same cloud as the \ion{C}{4},
and supports the maximal \ion{C}{4} model above. Thus, we conclude
that the maximal \ion{C}{4} model with three \ion{C}{4} components is
more appropriate to describe this system.

As with System~1, the right wing of the Ly$\alpha$\ profile is not produced
either by the \ion{Mg}{2} cloud or by the three \ion{C}{4} clouds.
Therefore, we placed an additional Ly$\alpha$\ line in this region of the
profile.  The unexplained absorption cannot be Ly$\beta$\ because this
region is not in the Ly$\beta$\ forest, as was the case with System~1.  We
fit the region using one additional Ly$\alpha$\ component with $\log N_{\rm HI}$\
$\sim$ 14.6 and $b$ $\sim$ 34~{\ifmmode{{\rm km~s}^{-1}}\else{km~s$^{-1}$}\fi}. This component is not
well-constrained, but for $\log U$ $>$ $-$2.0 the \ion{C}{4} and
\ion{N}{5} are over-produced unless $\log Z$ $<$ $-$2.0.  For $\log U$
$<$ $-$2.5, the metallicity can be as large as $\log Z$ $=$ 0.0 before
low ionization transitions are over-produced.

We conclude that the $z$ $=$ 1.68079 system in the HE~0429-4091
spectrum has three \ion{C}{4} components as well as one \ion{Mg}{2}
component.  The transitions that provided constraints on $\log U$ and
$\log Z$ based on our Cloudy models are listed in Table~\ref{tab4}.
All components have similar ionization parameters, $\log U$ $\sim$
$-$1.9.  It is difficult to constrain the metallicity of the low
ionization phase for this system since there are kinematically
distributed \ion{C}{4} components that contribute to the Ly$\alpha$\
absorption, however metal line absorption is best reproduced with
about solar metallicity for the \ion{Mg}{2} component.  The blueward,
offset \ion{C}{4} cloud is constrained by the Ly$\alpha$\ profile to have
$\log Z$ $>$ $-$1.1.  Because of its higher ionization state, the size
of the \ion{Mg}{2} cloud of this system ($\sim$ 100~pc) is two orders
of magnitude larger than that of the \ion{Mg}{2} cloud in System~1.
All components are optically thin. Ranges of acceptable model
parameters are listed in Table~\ref{tab5}. An example of the best
models and its model parameters are presented in Figure~\ref{fig2} and
Table~\ref{tab6}. For the model curves shown in Figure~\ref{fig2}, the
nitrogen abundance is decreased by 0.5 dex from the solar abundance.

\subsection{System~3 (HE~2243$-$6031; $z$=1.75570)}
\label{sec:sys3}

This weak \ion{Mg}{2} system has simple, single component profiles
detected in Ly$\alpha$, the low-ionization lines, \ion{Mg}{2} and
\ion{Si}{2}, and in the higher-ionization lines, \ion{Al}{3},
\ion{Si}{3}, \ion{Si}{4} and \ion{C}{4}. \ion{Fe}{2} and \ion{Al}{2}
lines are not detected to a 5$\sigma$ equivalent width limit.  See
Table~\ref{tab2} for the rest-frame equivalent widths and $5\sigma$
limit values.  The regions of spectrum where \ion{O}{1}, \ion{C}{2},
and \ion{N}{5} would appear are contaminated by blends, so that they
only serve to provide weak upper limits.  It is immediately apparent
that the equivalent width of \ion{Mg}{2} is very large for this system
compared to its relatively weak Ly$\alpha$\ profile. This implies that the
system has a high metallicity.  Velocity plots of various transitions
are displayed in Figure~\ref{fig3}.

As usual, we first optimize on the column density and Doppler
parameter of \ion{Mg}{2} obtained from Voigt profile fitting (see
Table~\ref{tab3}). A lower limit on the ionization parameter of $\log
U$ $>$ $-$3.5 can be determined, since \ion{Fe}{2} is over-produced at
smaller values.  On the other hand, higher ionization conditions with
$\log U$ $>$ $-$2.0 are ruled out because of overproduction of
\ion{Si}{4}.  The observed \ion{C}{4} absorption is also over-produced
if $\log U$ $>$ $-$1.5.

More stringent constraints on the ionization parameter can be found if
we consider the observed absorption in \ion{Al}{2} and \ion{Si}{2}.
For metallicities $\log Z$ $\lower0.3em\hbox{$\,\buildrel <\over\sim\,$}$ 0.0, the constraints on ionization
parameter are independent of metallicity.  For larger metallicities,
cooling leads to different constraints.  We will first consider the
constraints for $\log Z$ $<$ 0.0.  Considering \ion{Al}{2} and
\ion{Si}{2}, we found that an ionization parameter at the higher end
of the range $-$3.5 $<$ $\log U$ $<$ $-$2.0 produces the best
agreement, but these ions are still overproduced by the model.  To
reconcile a $\log U$ $=$ $-$2.0 model with the data for $\log Z$ $<$
0.0, we can reduce the abundance of aluminum by 0.5 dex, and that of
silicon by 0.2 dex.  Further variations in abundance pattern could
lead to different preferred values of $\log U$.

Independent of the constraint on the ionization parameter, the Ly$\alpha$\
profile can be used to place a strong lower limit on the metallicity
of the \ion{Mg}{2} cloud.  The \ion{Mg}{2} component is centered
exactly on the Ly$\alpha$\ profile, so that both the blue and the red side
of the Ly$\alpha$\ profile provide the same lower limit, $\log Z$ $>$ 0.0,
in order that Ly$\alpha$\ is not overproduced.
However, there are still some uncertainties on this metallicity
constraint. On one hand, lower values would be permitted if the system
is in $\alpha$-enhanced condition because this constraint does assume
a solar abundance pattern. On the other hand, the required metallicity
would be higher if magnesium is depleted onto dust grains, although we
do not expect a large amount of dust in weak \ion{Mg}{2} systems as
discussed below.
A cloud with $\log U$ $=$ $-$2.0 and $\log Z$ $=$ 0.0 has a
temperature of $T$ $=$ 9870~K and a size of $1.5$~kpc.  A smaller
$\log U$, which would be possible with further reductions in the
abundances of aluminum and silicon, would lead to a smaller cloud
size.

Metallicities greater than the solar value are also possible for the
\ion{Mg}{2} cloud.  In fact, a solar abundance pattern can be
consistent with the data because the cooling at high metallicity leads
to less \ion{Al}{3} and \ion{Si}{4} production.  If we accept the
constraint on ionization parameter of $\log U$ $\sim$ $-$2.0, as
discussed above, and assume a solar abundance pattern, we find that
\ion{Si}{4} will be slightly over-produced unless $\log Z$ $>$ 0.4.
It is over-produced on the blue side of the \ion{Si}{4}~$\lambda$1403
profile.  For $\log Z$ $>$ 0.4, \ion{Si}{4} is under-produced, but in
this case \ion{Si}{4} absorption can arise from a second, higher
ionization phase.  Similarly, \ion{Al}{3} is over-produced unless
$\log Z$ $>$ 0.9.  (The $\lambda$1863 transition is overproduced; the
$\lambda$1855 transition is contaminated by a blend.)  Due to cooling
at extreme metallicities, a $\log Z$ $=$ 1.0 cloud would have a
temperature of $\sim$ 500~K and a size of $\sim$25~pc.  These
constraints are quite dependent on abundance pattern.  As we noted
above, if the abundance of aluminum was decreased by 0.5 dex relative
to the solar pattern and the silicon abundance was decreased by 0.2
dex, $\log Z$ $=$ 0.0 would be consistent with the data for $\log U$
$=$ $-$2.0.  Both super-solar metallicity models and small changes in
abundance patterns are reasonable ways to reconcile models with the
data.  However, we emphasize that the metallicity constraint of $\log
Z$ $>$ 0.0 is directly determined based on comparison with the Ly$\alpha$\
profile.

Next, we added a high-ionization phase because the broad line profile
of \ion{C}{4} cannot be reproduced by the low-ionization \ion{Mg}{2}
component (see the dotted line in Figure~\ref{fig4} on the \ion{C}{4}
profile).  \ion{Si}{4} is also under-produced by an \ion{Mg}{2}
component with $\log U$ $=$ $-$2.0 and $\log Z$ $=$ 1.0.  In order to
find constraints on the \ion{C}{4} phase, we optimized on the
\ion{C}{4} column density from our Voigt profile fit. In this case,
the contribution to the \ion{C}{4} absorption from the low-ionization
phase was minimal, so that removing it made no significant difference
to the fit, listed in Table~\ref{tab3}.  From Cloudy photoionization
models, we derive a constraint of $\log U$ $=$ $-$1.5 -- $-$1.4 on the
ionization parameter, in order to simultaneously fit the \ion{C}{4}
and \ion{Si}{4} profiles.  If the metallicity of the \ion{Mg}{2} cloud
is lower so that it gives rise to more \ion{Si}{4} absorption, then
the ionization parameter of the \ion{C}{4} cloud would need to be
higher to compensate.  Similarly, if the ionization parameter of the
\ion{Mg}{2} cloud is lower, that of the \ion{C}{4} cloud would also
need to be lower in order to produce more \ion{Si}{4} absorption.

Since the Ly$\alpha$\ profile does not have a contribution from
kinematically separate \ion{C}{4} components, it can be used to place
a lower limit on the metallicity in the high-ionization phase.  In
order that Ly$\alpha$\ is not over-produced by the \ion{C}{4} cloud, a
metallicity of $\log Z$ $>$ $-$0.2 is required.  No firm upper limit
can be placed since the Ly$\alpha$\ absorption would be fully accounted for
by the \ion{Mg}{2} cloud if $\log Z$ $=$ 0.0.  However, if the
\ion{Mg}{2} cloud has $\log Z$ $=$ 1.0 in order to produce the
observed \ion{Al}{3} absorption, it does not account for all of the
observed Ly$\alpha$\ absorption.  In this case, $\log Z$ $=$ $-$0.1 -- 0.0
for the \ion{C}{4} cloud would account for the majority of the Ly$\alpha$\
absorption.  A cloud with $\log U$ $=$ $-$ 1.4 and $\log Z$ $=$ 0.0
has a size of $\sim$ 5~kpc.

We conclude that the system can be reproduced with only two
components: one low-ionization and one high-ionization phase.  The
constraints on $\log U$ and $\log Z$ and the transitions they were
derived from are presented in Table~\ref{tab4} for two models.  Model
1 assumes a solar abundance pattern and relies on a super-solar
metallicity for the low-ionization phase to fit the data.  Model 2
allows an adjustment of the abundances of aluminum and silicon
relative to the solar pattern, but still requires at least a solar
metallicity.  Our favored two-phase model has a relatively high
ionization parameter of $\log U$ $=$ $-$2.0 for the \ion{Mg}{2} cloud,
and $\log U$ $=$ $-$1.4 for the \ion{C}{4} cloud.  However, with a
different abundance pattern, the ionization parameter for the
\ion{Mg}{2} cloud could be lower.  All components of acceptable models
are optically thin in {\rm H}~{\sc i}.  Ranges of the acceptable model parameters
are listed in Table~\ref{tab5}, including both models with $\log Z$
$=$ 0.0 and with $\log Z$ $>$ 0.9 for the \ion{Mg}{2} cloud. An
example of best fit models and their model parameters are presented in
Figure~\ref{fig3} (only for the model with $\log Z$ $>$ 0.9) and
Table~\ref{tab6} (for both).

\subsection{Effect of Photoionization Modeling Assumptions}
\label{sec:assumptions}

In using the Cloudy code to derive constraints for the three systems,
we have assumed that photoionization equilibrium applies.  This is
almost certain to be the case for hydrogen since the photoionization
timescale is short, on the order of $10^4$~years.  For metals,
photoionization equilibrium might not apply if the gas has cooled from
a higher temperature, in which case the gas would still be more
ionized than our calculations would suggest.  This would clearly
affect our model constraints.  However, our most important conclusion
of high metallicities for weak \ion{Mg}{2} absorber would only be
stronger if a larger fraction of the magnesium was in higher
ionization states.

Another important assumption we applied in our modeling was that of a
solar abundance pattern, unless a deviation was required by the data.
We note that deviations from the solar pattern are to be expected, and
may apply whether the data require them or not.  In other words, there
is a degeneracy between the parameters $\log U$ and $\log Z$ and the
abundance pattern.
One common type of abundance pattern deviation would be an
$\alpha$-enhanced condition. In this case, we tend to over-estimate
metallicity because our metallicity estimations are based on
magnesium, an $\alpha$-element. Another source of the deviation is
depletion onto dust grains.  We do not really expect a large amount of
dust in weak \ion{Mg}{2} absorbers because they are not typically
close to sites with large {\ifmmode{N_{\rm H\;I}}\else{$N_{\rm H\;I}$}\fi}\ or with current star formation.
However, it is important to note that if magnesium is depleted by some
factor, our inferred metallicity would {\it increase} by roughly that
same factor.

Finally, we consider the effect of changing the shape of the ionizing
radiation field.  For the previous calculations we had assumed a
Haardt and Madau model for the ionizing radiation from quasars and
star forming galaxies, with an escape fraction of 0.1 of ionizing
photons from galaxies \citep{haa96,haa01}.  We have also considered
the opposite extreme, in which only quasars contribute to the ionizing
radiation field.  For a given ionization parameter, the absence of a
stellar contribution leads to an increase in the relative number of
high energy photons.  With that change of spectral shape, we found
negligible difference in the properties derived for the low ionization
phase.  For the high ionization phase, we found that a lower
ionization parameter (e.g., by about 0.5~dex) would be needed to fit
the data.
During the photoionization modeling presented above, we neglected
possible contributions from nearby stellar sources. However, we see
only small differences in the column densities of
low/intermediate-ionization phase gas (e.g., \ion{Si}{2}, \ion{Si}{3},
and \ion{Si}{4}) by a factor of 2 or 3, even after adding stellar
radiation from an {\rm O7} star with an effective temperature of
$T_{eff}$ = 38000~K (model\# C3 of \citealt{sch97}).  We adopt an
ionizing photon number density 10 times greater than that of the
extragalactic background radiation, as a maximum flux model. Thus,
contributions from nearby stellar sources, even if they exist, have no
significant qualitative effect on our conclusions.

We conclude that none of the assumption behind our photoionization
models have a qualitative impact on our conclusions.  In particular,
there is no effect that works against our conclusion of a very high
metallicity for System 3.

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

In this study, we applied photoionization models to constrain the
physical properties of three single-cloud weak \ion{Mg}{2} absorption
systems at $z$ $\sim$ 1.7.  Along with results presented in
\citet{lyn07}, these complete a sample of the seven $z$ $>$ 1.5
single-cloud weak \ion{Mg}{2} absorbers found in the VLT archive
\citep{nar07} for which metallicity constraints could be derived
because of simultaneous coverage of various metal lines and Ly$\alpha$.

We start this final section by summarizing the model constraints
derived for the three $z$ $\sim$ 1.7 single-cloud weak \ion{Mg}{2}
absorbers studied in this paper (\S~\ref{sec:summary}). We then
proceed to compare our results to the other four single-cloud weak
\ion{Mg}{2} systems at $z$ $\sim$ 1.7 from \citet{lyn07} in
\S~\ref{sec:comp_weakMgII}. In this section, we also compare to lower
redshift single-cloud weak \ion{Mg}{2} systems from the literature to
investigate evolutional trends. We address the high metallicity of
these systems in the context of DLAs (whose metallicities are usually
smaller than the solar value by a factor of 10 -- 30) in
\S~\ref{sec:comp_DLA}, sub-DLAs in \S~\ref{sec:comp_subDLA}, and
strong and multiple-cloud weak \ion{Mg}{2} systems (the strong ones
are usually associated with bright ($L$ $>$ 0.05$L^{*}$) galaxies
within 40$h^{-1}$ kpc; e.g., \citealt{ber91}) in
\S~\ref{sec:comp_strongMgII}, respectively.  We finally discuss a
possible correlation between metallicity and total hydrogen column
density in \S~\ref{sec:metal}.

\subsection{Summary of Results}
\label{sec:summary}

System~1 can be fit with a single low-ionization \ion{Mg}{2} cloud and
three higher ionization clouds (two of them offset in velocity from
the \ion{Mg}{2}) that give rise to \ion{C}{4} absorption.  The low
ionization cloud has $\log Z$ $>$ $-$0.8 and two of the high
ionization clouds must have metallicities at least this high as well.
Two offset low-metallicity clouds, producing only Ly$\alpha$\ absorption
appear to be clustered with this system.

System~2 is also fit with one \ion{Mg}{2} and three \ion{C}{4} clouds,
though all four clouds have similar, relatively high ionization
parameters ($\log U$ $\sim$ $-$1.9).  The metallicity of the
\ion{Mg}{2} cannot be well-constrained directly from comparison to the
Ly$\alpha$\ because the kinematical spread of the \ion{C}{4} lines give rise
to a large Ly$\alpha$\ equivalent width, however near solar metallicities
are possible, and even favored for a solar abundance pattern.  Also,
the blueward \ion{C}{4} cloud will over-produce absorption in Ly$\alpha$\
unless its metallicity is $\log Z$ $>$ $-$1.1.  A separate, offset
low-metallicity component was again needed to fit the red portion of
the Ly$\alpha$\ profile.

System~3 has the best metallicity constraint among the seven
single-cloud weak \ion{Mg}{2} absorbers at $z$ $>$ 1.5 because of the
absence of offset \ion{C}{4} components, thereby keeping the Ly$\alpha$\
profile very narrow.  Two phases are still required, with clouds of
different ionization parameters centered at the same velocity that
produce both the narrower low-ionization transitions, and the broader
high-ionization transitions.  Both the low- and the high-ionization
phases are constrained to have solar or super-solar metallicities.

The most striking thing about these model results, is that in all
three $z$ $\sim$ 1.7 systems, the metallicity of at least one phase of
gas is constrained to be greater than one tenth the solar value.  The
three systems all have at least two phases of gas: the low ionization
phase that arises in a layer of gas $\sim$1 -- 100~pc thick with a
density of 0.001 -- 0.1~{\ifmmode{{\rm cm}^{-3}}\else{cm$^{-3}$}\fi}, and the high ionization phase that
comes from a larger region ($\sim$ 0.1 -- 10~kpc) with a lower
density.  For these systems, even metallicities of the larger
high-ionization phase regions range from one tenth solar up to solar.
We will return to a discussion of the surprising issue of how regions
with $\log N_{\rm HI}$\ $<$ 15 (i.e., over a million times less than the
threshold for star formation) can be enriched in metals to the solar
value.

\subsection{Comparison to Other Single-Cloud Weak \ion{Mg}{2} Systems}
\label{sec:comp_weakMgII}

The properties of the three $z$ $\sim$ 1.7 single-cloud weak
\ion{Mg}{2} systems, derived from the observed profiles, are compared
to those of 23 other single-cloud weak \ion{Mg}{2} absorbers in
Figure~\ref{fig5}, both at similar and at lower redshifts.
Figure~\ref{fig5} shows that there is no significant evolution in the
observed \ion{Mg}{2}~$\lambda$2796 profiles.  Figure~\ref{fig6}
presents quantities derived from photoionization models, for our three
absorbers and for others taken from the literature, using the same
methods as we have used here \citep{rig02,cha03,din05,lyn07}.  The
specific constraints that are displayed in Figures~\ref{fig5} and
\ref{fig6} are also listed in Table~\ref{tab7}, with the relevant
references.

\citet{lyn07} derived similar constraints for the phase structure of
four different single-cloud weak \ion{Mg}{2} absorbers at $z$ $\sim$
1.7.  They also placed constraints on the metallicities of those
absorbers, of $\log Z$ $>$ $-$1.5, $>$ $-$1.5, $>$ $-$1.0, and $>$
$-$2.0 for the four different low-ionization phases.  They note that
the metallicities could be significantly higher than these lower
limits, however, these absorbers (by chance) tended to have larger
contributions to the Ly$\alpha$\ absorption from offset \ion{C}{4} clouds
and from Ly$\alpha$-only clouds.

Results from modeling lower redshift (0.4 $<$ $z$ $<$ 1.4)
single-cloud weak \ion{Mg}{2} systems also typically yield a two-phase
structure, with similar densities.  Metallicity constraints are again
often limited because of the difficulty in separating the
contributions of the two phases.  However, of the 11 cases of
single-cloud weak \ion{Mg}{2} absorbers for which some metallicity
constraint could be obtained, 2 cases require the metallicity of the
low-ionization phase to be be solar or even super-solar \citep{cha03}.
Furthermore, a total of 7/11 of the low redshift cases required a
metallicity of at least one tenth of the solar value.  We note that
the metallicities are often quite likely to be higher than these
strict lower limits because they are derived assuming that none of the
Ly$\alpha$\ absorption comes from the high-ionization clouds.

We conclude that some fraction of weak \ion{Mg}{2} systems have been
demonstrated to have solar or super-solar metallicity, at low redshift
(2 systems at {\ifmmode{z_{abs}}\else{$z_{abs}$}\fi}\ = 0.8181 and 0.9056; \citealt{cha03}) and even at
$z$ $\sim$ 1.7 (1 system at {\ifmmode{z_{abs}}\else{$z_{abs}$}\fi}\ = 1.7557; this paper), which
corresponds to an age of the Universe of 3.7~Gyr.  Values constrained
to be greater than 1/10th the solar value are even more common.
Furthermore, it is likely that in many cases for which lower limits on
metallicity have been derived the actual value is much higher.  The
metallicity constraint that we have derived for the {\ifmmode{z_{abs}}\else{$z_{abs}$}\fi}\ = 1.7557
absorber in this paper is high because this system is a rare case that
does not have offset \ion{C}{4} absorption that contributes to the
Ly$\alpha$\ absorption.  In other ways there is no reason to think it is
unusual.  The same applied for the low redshift weak \ion{Mg}{2}
systems with the highest metallicity constraints \citep{cha03}.  We
therefore know that single-cloud weak \ion{Mg}{2} absorbers at both
low and high redshift definitely have metallicities of at least 1/10th
solar, but that some and probably many have solar or even supersolar
metallicities.

Thus even before $z=1.5$ star formation must have polluted certain
environments with metal-rich gas.  Since \citet{lyn07} could not
constrain any $z$ $>$ 1.5 single-cloud weak \ion{Mg}{2} absorbers to
have close to solar metallicity, and since several lower redshift
solar metallicity cases were known, they tentatively suggested that
there might be a gradual build-up of metals in the single-cloud weak
\ion{Mg}{2} absorber population.  In this paper, we have shown there
is at least one strong counter-example, i.e. of $z$ $>$ 1.5 solar
metallicity cases, and therefore there appears to be no metallicity
evolution in the population.  This is not to say that the population
has a narrow range of physical properties: there is a large spread of
cloud densities and metallicities at all redshifts (see
Figure~\ref{fig6}).  However, it does suggest that common processes
are at work to create the metals, and that the metals are mixed into
fairly similar surrounding environments.

The existence of high metallicity (greater than solar) compact
intergalactic clouds ($\sim 100$~pc) has also been demonstrated by
\citet{sch07}, who derived robust lower limits on metallicities of
\ion{C}{4} absorbers using similar photoionization modeling
techniques.  \citet{sch07} find that the number density of this
population of high metallicity \ion{C}{4} absorbers is of the same
order as that of single-cloud weak \ion{Mg}{2} absorbers, and suggests
that the former might evolve from the latter as material expands.
More observations allowing direct connections between these
populations will be of great interest.  In fact, two of our three $z$
$>$ 1.5 weak \ion{Mg}{2} systems have \ion{C}{4} phases with
metallicities comparable to those derived by \citet{sch07} for the
high metallicity \ion{C}{4} cloud population.

\subsection{Comparison to Damped Ly$\alpha$\ Systems}
\label{sec:comp_DLA}

Damped Ly$\alpha$\ systems (DLAs), whose neutral hydrogen column densities
are greater than $\log N_{\rm HI}$\ $=$ 20.3, dominate the neutral gas in the
Universe at high redshift \citep[e.g.,][]{lan95,sto00,per03}. By
identifying metal absorption lines that correspond to high redshift
DLAs, it is possible to constrain the global metal abundances and
trace the evolution of metallicity in neutral gas from $z$ $=$ 5 to 0
\citep[e.g.,][]{pet94,pro00}.

The metallicity evolution, measured from DLAs, however, has two big
problems: (i) the global mean metallicity is less than 1/6th of the
solar abundance at $z$ $>$ 1 \citep[e.g.,][]{pet99,pro00}, and is
still less than the solar by a factor of $>$ 4 if the metallicity
redshift relationship is extrapolated down to $z$ $\sim$ 0 using a
least-squares linear fit \citep{pro03,kul05,per06b}
\footnote{At $z$ $<$ 1, there are only a few cases of metallicity
measurements in DLAs using high quality spectra, and in these cases
the metallicity is evaluated to be $\sim$ 1/10th of the solar value
\citep{var00,pet00}}, and (ii) the metallicity is also much smaller
than the theoretically expected value from the star forming activity
at higher redshift \citep[e.g.,][]{mad98,pet99,wol03}.  This
discrepancy (i.e., ``missing-metals problem'') remains an unresolved
problem.

There are at least two possible classes of ideas proposed for the
origin of this discrepancy: (i) a large amount of dust in metal-rich
DLAs obscures background quasars, which prevent us from detecting DLAs
with higher metallicities \citep[e.g.,][]{boi98,fal93,vla05}, and (ii)
the DLA region is sampling a lower metallicity than is found in the
star-forming parts of the galaxy \citep{ell05a,wol03}.

The dust obscuration scenario is supported by the small abundances of
depleted elements such as Cr and Fe, relative to undepleted elements
like Zn, in DLA systems \citep[e.g.,][]{pet97,kha04}. \citet{vla05}
estimated that $\sim$30~\%\ to 50\%\ DLAs are missed as a consequence
of obscuration.  In absorbers expected to have a lower dust abundance,
sub-DLA systems (super-LLSs) with $\log N_{\rm HI}$\ $=$ 19.0 -- 20.3, a higher
metallicity has been measured
\citep[e.g.,][]{per06a,pro06,pet00,kul07}. Although \citet{ell01} did
not see any significant differences in DLAs toward radio-selected
quasars, compared to those toward optically selected quasars at $z$
$>$ 2, the dust-depletion effect is expected to be more important at
lower redshift.

Alternatively, \citet{ell05a} and \citet{che05} suggest that there
could be an metallicity gradient as a function of a distance from the
galactic center, with DLAs produced at a larger impact parameter than
the stellar emission.  DLAs that are produced in the host galaxies of
gamma-ray bursts tend to have higher metallicities
\citep[e.g.,][]{djo04}, which could be more representative of typical
star-forming regions.  Also, small regions with DLA column densities
could be ejected from galaxies by superwinds or as a part of AGN
outflows \citep{mac99,ham97,gab06,bou06}.

Figure~\ref{fig7} summarizes the metallicities derived for DLAs at $z$
$=$ 0.4 to 4.8 from \citet{pro03} and \citet{kul05}, and compares to
the values for single-cloud weak \ion{Mg}{2} absorbers.  It is
striking that weak \ion{Mg}{2} absorbers, despite having drastically
smaller {\ifmmode{N_{\rm H\;I}}\else{$N_{\rm H\;I}$}\fi}\ values, have considerably higher metallicities than the
DLAs, both at low redshift and at $z$ $\sim$ 1.7.

Because they are not near regions with current star formation, the
dust-to-gas ratio in single-cloud weak \ion{Mg}{2} absorbers is
expected to be very small.  Even for DLAs, it has been demonstrated
that there is no bias against observing high metallicities due to dust
obscuration of their background quasars \cite{ell01,ell05b}.  Even if
the dust-to-gas ratio is similar to those in DLAs, their total dust
column density should be much smaller than those in DLAs, because of
their lower gas column densities.  Thus, they would not be subject to
the metallicity bias that was more likely to have affected DLAs.  Thus
the weak \ion{Mg}{2} absorbers provide an opportunity to see some high
metallicity regions of the universe at high redshift, though they are
not the same type of regions that would produce high metallicity DLA
absorption.

Relating to the idea that DLAs sample low metallicity parts of
galaxies, the weak \ion{Mg}{2} absorbers, despite their low column
densities, must somehow sample higher metallicity regions which are
not closely related to luminous galaxies.  This provides clear
evidence for metallicity inhomogeneities, which because of the size
constraints for single-cloud weak \ion{Mg}{2} absorbers, are on very
small scales (parsecs to hundreds of parsecs).

\subsection{Comparison to sub-DLA Absorbers}
\label{sec:comp_subDLA}

Sub-DLA absorbers, with $\log N_{\rm HI}$\ $=$ 19.0 -- 20.3, have been found to
have systematically higher metallicities than DLAs with the
{\ifmmode{N_{\rm H\;I}}\else{$N_{\rm H\;I}$}\fi}-weighted value larger by 0.5 -- 0.8 dex, at 0.6 $<$ $z$ $<$ 3.2
\citep{kul07}.  Several sub-DLA systems (super-LLSs) at 0.5 $<$ $z$
$<$ 2.0 have been found with metallicities that are solar or
super-solar \citep[e.g.,][]{per06a,pro06}.  Even small numbers of such
high metallicity systems contribute significantly to the total metal
mass density at high redshift \citep{pro06,kul07}.  The \citet{kul07}
mean values of metallicity for sub-DLAs, determined from [Zn/H], are
reproduced in our Figure~\ref{fig7} from their Table~1.  Although
there is a range of metallicities at each redshift, there is a clear
increase in the metallicity of sub-DLAs from $\langle z \rangle$ $=$
1.5 to $\langle z \rangle$ $=$ 0.9.  This is consistent with the
global star formation history in the universe over this period, and
the metallicity values are much more in line with expectations from
chemical evolution models than those for DLAs \citep{kul07,som01}.

Some single-cloud weak \ion{Mg}{2} absorbers, with $\log N_{\rm HI}$\ $\sim$ 15
-- 16, also have solar or super-solar metallicities, both at $\langle
z \rangle$ $=$ 0.9 and at $\langle z \rangle$ $=$ 1.5.  There is
considerable overlap between the metallicity constraints for sub-DLAs
and single-cloud weak \ion{Mg}{2} absorbers at low and intermediate
redshifts.  However, it is important to note that many of the weak
\ion{Mg}{2} metallicities are only lower limits, and it is quite
likely that the values are higher at least in some cases.  The reason
that we cannot derive higher metallicities for those systems could
simply be the separate \ion{C}{4} clouds that happen to overlap in
velocity.
It is also important to note that for $z$ $\sim$ $1.7$, the redshift
of the absorbers we have studied in this paper, there are not many
sub-DLA systems with super-solar metallicity, but all single-cloud
weak \ion{Mg}{2} absorbers potentially could have solar- or
super-solar metallicity because we can place only lower limits on
metallicity. Although two super-solar metallicity sub-DLAs were
investigated by \citet{pro06}, they were selected for study from a
much larger sample because of their extremely strong \ion{Zn}{2}
absorption. However, much larger data samples of sub-DLAs and weak
\ion{Mg}{2} systems will be required before concluding whether there
are any differences of metallicity between these two categories.

Another apparent difference between the metallicities of single-cloud
weak \ion{Mg}{2} absorbers and sub-DLAs, is that there is no clear
increase in the metallicity of the single-cloud weak \ion{Mg}{2}
absorbers over the range 0.4 $<$ $z$ $<$ 1.7.  This type of absorption
is apparently produced by the same types of regions present over this
redshift range, though they are known to be less common at low
redshifts and at $z$ $>$ 2 \citep{nar05,nar07,lyn06,lyn07}.

The sample size of weak \ion{Mg}{2} systems is now much smaller than
those for DLAs and Sub-DLAs, and only a few cases have been confirmed
to have super-solar metallicities. It is quite important to increase
the sample size of weak \ion{Mg}{2} absorbers, particularly to include
other systems with only one \ion{C}{4} component such that a strong
constraint on metallicity can be derived.  At $z < 1.5$ this will
involve high resolution observations with {\it HST}/COS or {\it
HST}/STIS of quasars for which high-resolution optical spectra are
also available.

\subsection{Comparison to Strong \ion{Mg}{2} and Multiple-Cloud Weak \ion{Mg}{2} Absorbers}
\label{sec:comp_strongMgII}

Intermediate in $\log N_{\rm HI}$\ between single-cloud weak \ion{Mg}{2}
absorption and sub-DLA's are Lyman limit and partial Lyman limit
systems, with $\log N_{\rm HI}$\ $=$ 16 -- 19.  These populations loosely
correspond to the bulk of the strong \ion{Mg}{2} absorber population,
with $W_r$(2796) $>$ 0.3~{\AA}, and to multiple-cloud weak \ion{Mg}{2}
absorbers.  Fewer of these systems have been studied in detail to
date, but it is still possible to derive good constraints on
metallicities in many cases, using photoionization modeling techniques
the same as those employed here.  The ten cases of strong and
multiple-cloud weak \ion{Mg}{2} absorbers at 0.7 $<$ $z$ $<$ 1.9
\citep{mas05,zon04,lyn07,pro99,din03}, shown on Figure~\ref{fig7},
have metallicities ranging from $-$1.5 $<$ $\log Z$ $<$ 0.6.  These
are consistent with the ranges for both sub-DLAs and single-cloud weak
\ion{Mg}{2} absorbers, but tend to be higher than those for DLAs.
Without larger samples, we cannot distinguish differences or
evolutionary trends, but again we note that many of these are
measurements of the metallicity, while for the single-cloud weak
\ion{Mg}{2} absorbers the low metallicity values are all lower limits.

\subsection{Metallicity vs. {\ifmmode{N_{\rm H\;I}}\else{$N_{\rm H\;I}$}\fi}}
\label{sec:metal}

In summary, Figure~\ref{fig7} shows that it is the lowest {\ifmmode{N_{\rm H\;I}}\else{$N_{\rm H\;I}$}\fi}\
systems, that appear to have the highest metallicities, both at $z$
$<$ 1 and at $z$ $\sim$ 1.7.  The data are even consistent with a
gradual increase in metallicity with decreasing {\ifmmode{N_{\rm H\;I}}\else{$N_{\rm H\;I}$}\fi}, though there is
a huge spread at any given value. Such a trend has already been
pointed out using compiled samples of DLAs and sub-DLAs
\citep{boi98,ake05,mei06,kha07}, extending down to sub-DLAs with
$\log N_{\rm HI}$\ $\sim$ 19 \citep{per03}. \citet{yor06} also suggested that a
similar trend continued to strong \ion{Mg}{2} systems with smaller
hydrogen column densities.

The individual systems plotted on Figure~\ref{fig7}, however, are only
the systems with detected \ion{Mg}{2} absorption.  In particular, for
single-cloud weak \ion{Mg}{2} absorbers we are sure to be missing
objects with metallicities significantly less than those of the
systems that we do detect.  In fact, those lower metallicity objects
are part of the much larger Ly$\alpha$\ forest population.
Figure~\ref{fig7} also shows the typical metallicities of Ly$\alpha$\ forest
clouds, both those with $\log N_{\rm HI}$\ $>$ 14.5, which have $\log Z$ \lower0.3em\hbox{$\,\buildrel <\over\sim\,$}\
$-$2.0 \citep{son96,cow95,tyt95}, and those with lower column
densities, which have a mean metallicity of $\log Z $ \lower0.3em\hbox{$\,\buildrel <\over\sim\,$}\ $-$3.2
\citep{cow98,lu98}.  These much lower values for the Ly$\alpha$\ forest
overlap with the DLA metallicities at $z$ $>$2.5, but are smaller than
those for low redshift DLAs.

The fact that the high metallicity population of single-cloud weak
\ion{Mg}{2} absorbers represents only a small fraction of all Ly$\alpha$\
forest clouds does not diminish their significance.  The single-cloud
weak \ion{Mg}{2} absorbers at $z$ $\sim$ 1 would account for 25 --
100\%\ of the Ly$\alpha$\ forest clouds with 15.8 $<$ $\log N_{\rm HI}$\ $<$ 16.8
\citep{rig02}.  More importantly, the cross-section of single-cloud
weak \ion{Mg}{2} absorbers in the plane of the sky is comparable to
that of galaxies at $z$ $\sim$ 1, considering the regions of galaxies
$\lower0.3em\hbox{$\,\buildrel <\over\sim\,$}$30 -- 40~kpc that produce Lyman limit absorption.  Thus there
are significant regions of the universe covered by these mysterious
near-solar or super-solar metallicity objects with small $\log N_{\rm HI}$.

As described above, one explanation proposed for the larger
metallicities for sub-DLAs than for DLAs is the dust obscuration of
quasars that have the highest metallicity DLAs in their foregrounds.
Although there may be some debate about whether this is a plausible
explanation of that trend, it is not likely to explain the continued
increase in metallicity toward the lowest $\log N_{\rm HI}$\ weak \ion{Mg}{2}
absorbers.  There is not likely to be a large selection effect due to
dust obscuration for even the sub-DLA systems.  It seems more
plausible that different types of absorbers are sampling different
types of regions in and around galaxies, which can have large
variations in metallicities even on small scales.  Different types of
galaxies, e.g. star-forming versus quiescent, will have different
fractions of area covered by regions of the different types. By
comparing to the mass-metallicity relation seen in star-forming
galaxies at local universe \citep{tre04} and higher redshift
\citep{sav05,erb06}, \citet{yor06} and \cite{kha07} proposed that the
DLAs are associated with low-mass ($<$ $10^9$ $M_{\odot}$) galaxies,
while sub-DLAs and weaker LLSs are probably the systems that arise in
massive spiral/elliptical galaxies.

The question remains: how do the single-cloud weak \ion{Mg}{2}
absorbers develop such high metallicities even at redshifts as high as
$z$ $\sim$ 1.7?  We know that the lines of sight that pass through
these objects have {\ifmmode{N_{\rm H\;I}}\else{$N_{\rm H\;I}$}\fi}\ five or six orders of magnitude below the
star formation threshold.  We also know that most single-cloud weak
\ion{Mg}{2} absorbers are not located very close to luminous galaxies,
but that they tend to be in the vicinity \citep{chu05,mil06}.
Outflowing gas from dwarf galaxies in the local universe
\citep[e.g.,][]{sto04} or from Lyman break galaxies at higher redshift
\citep[e.g.,][]{ade05} could partially contribute to the high
metallicities in single-cloud weak \ion{Mg}{2} systems.  \citet{sch07}
have argued that high metallicity intergalactic \ion{C}{4} clouds may
be ejected as small clumps in superwinds from star-forming galaxies.
However, at least some of the very high metallicity single-cloud weak
\ion{Mg}{2} absorbers show signs of \ion{Fe}{2} line detections
(e.g. our System~1 and several systems in \citealt{rig02}), although
sometimes it is detected with only a few $\sigma$ level. For those
systems, it has been stated that they are not $\alpha$-enhanced, and
that "in situ" star formation, including the less energetic Type~Ia
SNe that increase the iron to magnesium ratio, is responsible for
enriching the gas \citep{rig02}.
Here, readers should be aware that it is not the gas along the weak
\ion{Mg}{2} absorber line of sight that has formed stars, since it has
much too low a column density.  It must instead be a nearby molecular
cloud region that was actively star forming more than 1 billion years
before the time from which we observe the absorber.  That region would
itself produce DLA absorption if a line of sight passed through it.
It is known that many DLAs at $z$ $<$ 1 are associated only with dwarf
or low surface brightness galaxies \citep{rao03}, and that in at least
one case no galaxy is detected in deep narrow-band H$\alpha$ images
\citep{bou01}.

So then why do DLAs not have high metallicities as well?  And how do
we explain the general trend of increasing metallicity with decreasing
{\ifmmode{N_{\rm H\;I}}\else{$N_{\rm H\;I}$}\fi}?  The reason could be simple.  Large quantities of metals can be
produced in star forming regions in a variety of environments, inside
and outside of galaxies.  However, the metals will be dispersed into
widely different types of regions, depending upon the environment.
Within giant galaxies molecular cloud regions are typically surrounded
by a significant volume of high density, high column density gas.  So
the metallicity is diminished by a large dilution effect as the
ejected metals are spread through the galaxy region.  In an
intergalactic star-forming region, which is metal-polluted by {\it
in-situ} star formation or has been ejected from a nearby dwarf
galaxy, the clouds may not have a significant volume of high column
density gas around it, thus the dilution could be minimal.  This would
especially apply in a shallow potential well, where much of the
surrounding gas is ejected by the early Type~II SNe, and the remnant
low density gas remaining is greatly enhanced by the metals produced
by later Type~Ia's which encounter a relatively sparse surrounding
medium.  The same type of dilution effect could lead to DLAs having
lower metallicity than sub-DLAs.  \citet{pro06} suggests that the
super-solar metallicity sub-DLAs that they observed at $z$ $\sim$ 1.7,
extreme members of the sub-DLA population, are related to regions
affected by the feedback from actively star-forming regions.  The
metals would be spread to larger lower density sub-DLA regions by
energetic processes.

Obviously, the metallicity of a region is determined by a combination
of the amount of star formation activity in its vicinity and the
overall density of its gas.  In the case of the single-cloud weak
\ion{Mg}{2} absorbers, we conclude that the density effect dominates,
leading to unexpected patches of high metallicity outside of galaxies.
In the contrasting case of a DLA, there is just too much gas around it
to dilute the metals coming through it, and a low metallicity often
results.

\acknowledgments
This research was funded by the National Science Foundation (NSF)
under grant AST-04-07138 and by NASA under grant NAG5-6399. This work
was also partially supported by the Sumitomo foundation (070380). We
also acknowledge the ESO archive facility for providing the data. This
paper benefit from many constructive suggestions by an anonymous
referee.


Oxygen in the early Galaxy: OH lines as tracers of oxygen abundance in extremely metal-poor giant stars

Abstract: Oxygen is a powerful tracer element of Galactic chemical evolution. Unfortunately, only a few oxygen lines are available in the ultraviolet-infrared stellar spectra for the reliable determination of its abundance. Moreover, oxygen abundances obtained using different spectral lines often disagree significantly. In this contribution we therefore investigate whether the inadequate treatment of convection in 1D hydrostatic model atmospheres used in the abundance determinations may be responsible for this disagreement. For this purpose, we used VLT CRIRES spectra of three EMP giants, as well as 3D hydrodynamical CO$^5$BOLD and 1D hydrostatic LHD model atmospheres, to investigate the role of convection in the formation of infrared (IR) OH lines. Our results show that the presence of convection leads to significantly stronger IR OH lines. As a result, the difference in the oxygen abundance determined from IR OH lines with 3D hydrodynamical and classical 1D hydrostatic model atmospheres may reach -0.2 ... -0.3 dex. In case of the three EMP giants studied here, we obtain a good agrement between the 3D LTE oxygen abundances determined by us using vibrational-rotational IR OH lines in the spectral range of 1514-1626 nm, and oxygen abundances determined from forbidden [O I] 630 nm line in previous studies.

\section{Introduction}

Although oxygen is an important tracer of Galactic chemical evolution, only a few spectral lines are available for its diagnostics in extremely metal-poor (EMP) stars. In red giant stars, [O~I] 630~nm line can be used, however, in the EMP regime it is very weak and frequently blended with telluric oxygen lines ([O~I] and ${\rm O}_2$ band). Alternatively, oxygen abundances can be derived from ultraviolet (UV) and infrared (IR) OH lines but results obtained using molecular and atomic lines are frequently discrepant: for example, abundances obtained from IR~OH lines are typically $\sim0.3-0.4$~dex higher than those determined from the forbidden [O~I] 630~nm line.

One possible reason behind these discrepancies is that 1D hydrostatic model atmospheres that are used in the abundance analysis do not account for the horizontal fluctuations of thermodynamic quantities arising from convective motions in stellar atmospheres. Since formation of OH lines is very sensitive to the effects of convection \citep[see, e.g.,][]{CAT07,DKS13}, one may expect that application of 3D hydrodynamical model atmospheres in the abundance analysis may diminish, or even eliminate, the differences in oxygen abundances obtained from atomic and molecular lines.

In this contribution we therefore study whether a more realistic modeling of convection may bring the oxygen abundances obtained from molecular IR~OH lines into better agreement with those determined using forbidden [O~I] 630~nm line (note that [O~I] 630~nm line is insensitive to 3D hydrodynamical \textit{and} NLTE effects, see, e.g., \citealt[][]{A05}). To this end, we used spectra of three EMP giants obtained with the VLT~CRIRES spectrograph, as well as 3D hydrodynamical \COBOLD\ and 1D hydrostatic \LHD\ model atmospheres, to better understand the role of convection in the formation of IR~OH lines. All results obtained in this analysis will be presented in a forthcoming paper \citep[][]{DKB14}, while in this contribution we briefly summarize some of the most important findings.

\section{Methodology}

We used high resolution ($R = 50\,000$) and high signal-to-noise ratio ($S/N\approx400-600$) near-infrared ($H$-band) spectra of 3 EMP stars ($\FeH=-2.6$ to --3.1, see Table~1) obtained with the CRIRES spectrograph at the VLT~UT1 telescope (program ID 089.D-0079(A); a sample spectrum is shown in Fig.~\ref{fig:spectrum}). For each star, oxygen abundances were determined using the measured equivalent widths of vibrational-rotational IR~OH lines ($X^{-1}\Pi$) from the first-overtone sequence located in the spectral range of $1527-1626$~nm. Atmospheric parameters and iron abundances used in the abundance determination were taken from the literature.

Abundance analysis was done in 1D~LTE with the \ATLAS\ model atmospheres. The obtained 1D~LTE oxygen abundances were corrected for the 3D hydrodynamical effects, by applying 3D--1D abundance corrections to oxygen abundances determined from each individual IR~OH line. The 3D--1D abundance corrections (i.e., differences in the abundance of a given chemical element obtained from the same spectral line using 3D hydrodynamical and 1D hydrostatic model atmospheres) were computed with the 3D hydrodynamical \COBOLD\ \citep{FSL12} and 1D hydrostatic \LHD\ model atmosphere codes \citep{CLS08}. The two types of model atmospheres shared identical elemental abundances, atmospheric parameters, equation of state, and opacities. Abundance corrections were derived using a single 3D hydrodynamical model atmosphere computed using the following atmospheric parameters: $\teff=4600$~K, $\log g=1.6$, $\MoH=-2.5$. We note that differences between the atmospheric parameters of this 3D hydrodynamical model and those of the target EMP stars studied here did not exceed $\Delta \teff=150$~K, $\Delta \log g = 0.2$, and $\Delta \MoH = 0.6$. Such differences would lead to only minor deviations in the abundance corrections, the size of which (typically, a few hundredths of a dex) would render them unimportant in the context of the present study. The 3D hydrodynamical \COBOLD\ model atmosphere used in our work had a spatial resolution of $220\times220\times280$ elements, corresponding to the physical size of $3.85\times3.85\times2.21$~Gm in $x,y,z$ directions, respectively. The solution of radiative transfer equation was made using opacity binning technique \citep[][]{N82,L92}, by sorting opacities into 6 bins. The opacities were taken from the \MARCS\ model atmosphere package \citep[][]{GEK08}.

A sample VLT~CRIRES spectrum with the Gaussian fits to IR~OH lines is shown in Fig.~\ref{fig:spectrum}. The obtained 1D~LTE and 3D~LTE oxygen abundances, as well as abundances determined using forbidden [O~I] 630~nm line (literature data), are provided in Table~\ref{table:stars-abundances}.

\begin{table}[tb]
\caption{List of observed EMP stars and determined oxygen abundances.}
\label{table:stars-abundances}
\smallskip
\begin{center}
\scalebox{0.8}{
\setlength{\tabcolsep}{4pt}
\begin{tabular}{lcccccccc}
\hline
\noalign{\smallskip}
Star        &  $V$  &  $S/N$  & \teff, K  & \logg  & \FeH   & \multicolumn{2}{c}{Abundances from OH lines}           & Abundances from     \\
            &     &         &           &       & 1D LTE  & \multicolumn{2}{c}{(this work)}                          & [O I] 630 nm line   \\
            &     &         &           &       &         & $\OFe^{\rm{1D\,LTE}}$    &   $\OFe^{\rm{3D\,LTE}}$       & $\OFe^{\rm{1D\,LTE}}$ \\
\noalign{\smallskip}
\hline
\noalign{\smallskip}
HD 122563   & 6.2 & 560    & $4600^a$   & $1.6^a$  & $-2.8^b$  & $0.79 \pm 0.12$  & $0.56 \pm 0.13$ & $0.53^c$; $0.60^d$; $0.62^e$   \\
HD 186478   & 9.2 & 440    & $4700^b$   & $1.3^b$  & $-2.6^b$  & $1.02 \pm 0.10$  & $0.76 \pm 0.09$ & $0.75 \pm 0.11^e$              \\
BD -18:5550 & 9.4 & 600    & $4750^b$   & $1.4^b$  & $-3.1^b$  & $1.07 \pm 0.12$  & $0.84 \pm 0.15$ & $0.42 \pm 0.24^e$              \\
\noalign{\smallskip}
\hline
\end{tabular}
}
\end{center}
\begin{list}{}{}
\item[$^{\mathrm{a}}$] \citet{CTB12}; $^{\mathrm{b}}$ \citet{CDS04}; $^{\mathrm{c}}$ \citet{SKP91}; $^{\mathrm{d}}$ \citet{BMS03}; $^{\mathrm{e}}$ \citet{SCP05}.
\end{list}
\end{table}

\begin{figure}
\centering
\includegraphics[width=12cm]{fig01.eps}
\caption{Observed VLT~CRIRES spectrum of HD~122563 (black dots). Gaussian profiles fits to IR~OH lines used in the oxygen abundance determination are shown in red.\label{fig:spectrum}}
\end{figure}

\begin{figure}
\centering
\includegraphics[width=12cm]{fig02.eps}
\caption{\textbf{Top panel:} temperature profiles in the 3D hydrodynamical (gray scale is temperature probability density), average $\xtmean{\mbox{3D}}$ (red solid line), and 1D (red dashed line) model atmospheres, plotted versus the logarithm of Rosseland optical depth, $\log \tau_{\rm Ross}$. All models were computed using identical atmospheric parameters ($\teff=4600$~K, $\log g=1.6$, $\MoH=-2.5$), chemical composition, equation of state, and opacities. Horizontal bars indicate the range in optical depth where typical IR~OH line forms (i.e., its where equivalent width increases from 5\% to 95\%; parameters of this particular line were: wavelength $\lambda = 1605.5$~nm, excitation potential $\chi = 0.48$~eV): black and dashed blue bars correspond to the line forming regions in the full 3D and 1D model atmospheres (the equivalent widths of OH lines in both cases were 1.6~pm). \textbf{Bottom panel:} RMS value of horizontal temperature fluctuations in the 3D model (black line) and temperature difference between the average $\xtmean{\mbox{3D}}$ and 1D models (dashed blue line).\label{fig:T-profiles}}
\end{figure}

\section{Results and conclusions}

Results provided in Table~\ref{table:stars-abundances} show that oxygen abundances determined using 3D hydrodynamical model atmospheres are always smaller than those obtained with 1D hydrostatic models: 3D--1D differences for the three stars are in the range of $-0.23 \dots -0.26$~dex. These differences are a result of large horizontal fluctuations of thermodynamic and hydrodynamical quantities in 3D hydrodynamical model atmospheres of these stars that arise because of convective motions in their outer atmospheric layers. Horizontal temperature fluctuations are clearly seen in Fig.~\ref{fig:T-profiles}, where we show temperature profiles of the 3D hydrodynamical, average $\xtmean{\mbox{3D}}$, and 1D model atmospheres computed using identical atmospheric parameters, chemical composition, equation of state, and opacities. In fact, the presence of large temperature fluctuations (which increase towards the outer atmospheric layers characterized with lower $\log \tau_{\rm Ross}$ values) seems to be a general property of different types of stars and is seen routinely in their 3D hydrodynamical model structures \citep[see, e.g.,][]{LK12,KSL13,DKS13}. Large horizontal temperature fluctuations lead to larger OH line opacities in the low-temperature regions. This, in turn, results in stronger lines in 3D, negative abundance corrections, and thus, lower 3D~LTE abundances.

In fact, it is not only the horizontal temperature fluctuations that influence the formation and strengths of IR~OH lines. The differences between the average temperature profile of the 3D model and that in 1D hydrostatic model do play a role, too. Since the average $\xtmean{\mbox{3D}}$ model is slightly hotter in the optical depths where IR~OH lines form ($\log \tau_{\rm Ross}\approx-0.5\dots-4$; see Fig.~\ref{fig:T-profiles}), this results in somewhat weaker lines in $\xtmean{\mbox{3D}}$, and thus, small but positive abundance corrections. This slightly diminishes the effect of horizontal temperature fluctuations and leads to smaller (i.e., in the absolute sense) total 3D--1D abundance corrections.

At the same time, we obtain a good agreement between the 3D~LTE oxygen abundances and those obtained in 1D~LTE using forbidden [O~I] 630~nm line. The agreement is especially good in HD~122563 and HD~186478, where the difference between abundances obtained using the two indicators is smaller than $\sim0.02$~dex. Although differences are larger in BD~--18:5550, note that in this case 1D~LTE abundance measurement is less reliable, as indicated in \citet[][]{SCP05}.

The obtained results therefore demonstrate that the use of 3D hydrodynamical model atmospheres may allow to reconcile oxygen abundances obtained from IR~OH lines and [O~I] 630~nm line. Nevertheless, one should note that NLTE effects may play a role in the formation of OH lines at low metallicities. The impact here may be two-fold. First, formation and dissociation of OH molecules may be affected by non-equilibrium effects and thus non-equilibrium OH number densities may be different from those where equilibrium molecular formation is assumed. Second, IR~OH lines may form in the conditions that are far from LTE, thus NLTE line strengths may be different from those that would be expected in LTE. It is quite likely that both effects may be important in the atmospheres of EMP red giants, especially in the outer atmospheric layers were molecular lines form (note, however, that the forbidden [O~I] 630~nm line seems to be insensitive to 3D hydrodynamical and NLTE effects, see, e.g., \citealt{A05,DKS13}). In the future, these effects should be properly taken into account using full 3D~NLTE analysis methodology.

\acknowledgments{This work was supported by grants from the Research Council of Lithuania (MIP-065/2013) and the bilateral French-Lithuanian programme ``Gilibert'' (TAP~LZ~06/2013, Research Council of Lithuania). E.C. is grateful to the FONDATION MERAC for funding her fellowship. H.G.L. acknowledges financial support by the Sonderforschungsbereich SFB 881 ``The Milky Way System'' (subproject A4 and A5) of the German Research Foundation (DFG).}

\normalsize

Solar and Stellar Astrophysics (astro-ph.SR)