\section{Introduction}

The abundances of the rare light elements, Li, Be, and B, have revealed an
amazing amount of information about an array of astrophysical issues.  These
include cosmology and nucleosynthesis in the Big Bang, Galactic chemical
evolution, cosmic ray theory, and stellar interiors and stellar evolution.  It
is important to add more observations of Be to the large number of Li
observations in order to make new progress in these research fields.
Beryllium is the least complicated of the three light elements because it has
only one known source of production and only one stable isotope, $^9$Be.  It
is also less susceptible to depletion than Li in stellar interiors.

Observations of Primas et al.~(2000a, 2000b) of Be in very metal poor stars
led them to suggest the possibility of a Be plateau at very low [Fe/H],
possibly similar to the Spite lithium plateau (Spite \& Spite 1982).  The Li
plateau in older stars is attributed to nucleosynthesis of Li during the Big
Bang.  If such a plateau exists, could it be similar to the Li plateau?
Whereas the Li plateau is produced by nuclear reactions during the Big Bang,
very little Be is produced in standard Big Bang nucleosynthesis (BBN):
N(Be)/N(H) = 10$^{-17}$ (Thomas et al.~1994).  Models that include
inhomogeneities in the early universe could produce a different mixture of
light elements from the standard models (e.g. Malaney \& Mathews 1992).
Inhomogeneous models can create far more Be than standard models and the
predictions of Orito et al.~(1997) are 10$^{-15}$ to 10$^{-14}$; the implied
plateau corresponds to 5 x 10$^{-14}$, which is tantalizingly close to the
observed values.

Another possible cause of a Be ``plateau'' could be the creation of Be through
spallation reactions in contained superbubbles created by multiple supernovae
(Parizot 2000) or near hypernovae (Fields et al.~2002, Nakamura et al.~2006).
Enrichment of the ISM in Be by such processes early in the history of the
galaxy could result in a detectable plateau.

The lowest metal stars ([Fe/H] $<$ $-$3.0) were formed early in the evolution
of the Galaxy.  By tracing the abundance of Be in these old stars,
insight can be gained into the chemical history of the Galaxy.  Beryllium is
formed mainly through spallation reactions, which was proposed by Reeves,
Fowler \& Hoyle (1970) and further described by Meneguzzi, Audouze \& Reeves
(1971).  The basic process is that high energy ($\sim$150 MeV) protons and
neutrons bombard interstellar nuclei of C, N and O creating lighter isotopes.
It has also been suggested that the ``bullets'' and the ``targets'' might be
reversed near supernovae where C, N, and O nuclei could be accelerated into
the local interstellar gas including protons and neutrons (see, for example,
Duncan et al.~1997, 1998, Lemoine, Vangioni-Flam \& Cass\'e 1998). 

Recent studies of Be in low metal stars include Boesgaard, Levesque \& Rich
(2008), Tan, Shi \& Zhao (2009) and Smiljanic et al.~(2009).  All of these
studies examine the Galactic evolution of Be through the relationships between
the abundances of Be with Fe and Be with O and/or alpha-elements.  The large
study of Smiljanic et al.~(2009) includes 90 stars and the determination of Be
abundances of 73 halo and thick disk stars.  They relied on Hipparcos
parallaxes to find log $g$, and adopted published values for temperature,
[Fe/H], [$\alpha$/Fe], and microturbulent velocity.  Their stars range in
[Fe/H] from $-$0.5 to $-$2.5.

In this work we have obtained high signal-to-noise (S/N), high resolution
spectra with the upgraded HIRES on the Keck I telescope of 24 stars covering
values of [Fe/H] of $-$2.32 to $-$3.45 to examine Be abundances at as low
metallicities as possible.  We supplement this sample with reanalyzed
Keck/HIRES data from Boesgaard et al.~(1999a) and Boesgaard (2007) and
Subaru/HDs data from Boesgaard \& Novicki (2006).  Our total sample is 49
dwarfs and subgiants ranging in [Fe/H] from $-$0.5 to $-$3.5, with 34 stars
below $-$2.0.  This is the most metal-poor set of stars which have been
analyzed for Be abundances.  Another difference from previous studies is that
we have determined parameters spectroscopically.


\section{Observations and Data Reduction}
The Be II resonance lines, located in the UV spectral region at 3130.421 and
3131.065 \AA, were used to determine the Be abundances.  High S/N and
high-resolution spectra are required to observe the Be doublet, which is weak
in low metallicity stars.  The upgraded High Resolution Echelle Spectrograph
(HIRES) instrument on the Keck I telescope achieves the necessary requirements
(Vogt et al.~1994).  The spectral resolution for our setup was 48,000, and the
light collecting area of the 10 m Keck telescope enables us to obtain S/N near
100 per pixel for our metal-poor sample.  The pixel size on the upgraded CCD
is 15 $\mu$.  Such stars are uncommon and thus most are faint, especially in
the UV wavelength region of the Be II lines.  The observations reported here
were carried out over 11 nights using the blue cross-disperser.  The spectral
range is approximately 3000 to 6000 \AA.  The new HIRES CCD is composed of
three chips, one each for UV/blue (on which the Be doublet is located), green
and red wavelengths.  The quantum efficiency near 3100 \AA\ is about 93\%,
making exposures of stars down to V = 12 possible.

Stars were selected based on their metallicity and observability.  Stars with
[Fe/H] values of less than $-$2.4 dex were selected from a literature search.
To obtain the best spectra, stars were observed near the meridian, at the
lowest possible airmass and atmospheric dispersion.  Our aim was to observe
stars with a declination range of $-20 < \delta < +60$; the latitude of Mauna
Kea is +19$^{\deg}$ 45'.  In order to obtain spectra of multiple objects with
the required S/N, a cutoff at V of 12.5 was chosen.  Table 1 shows
observations of the stars analyzed with V magnitudes, literature values of
[Fe/H], the UT dates of the observations, the total exposure times and S/N per
pixel near the Be doublet; in the cases where there are exposures on multiple
nights the exposure times and S/N values are the totals.  The median S/N = 104
and the mean is 106 per pixel.

HIRES data reduction was carried out using the IDL HIREDUX pipeline and
standard IRAF routines.  The pipeline was used to bias-subtract,
flatten, extract and wavelength calibrate the spectrum from each individual
exposure.  IRAF was used to co-add multiple exposures of the same star and to
fit continua to each combined, calibrated spectrum.  G 64-37 was observed with
HDS on Subaru, which has a similar but smaller wavelength range, with a
maximum of $\sim$ 4500 \AA; see Boesgaard \& Novicki (2006).

In addition to the new data obtained with the upgraded HIRES we have included
data from Boesgaard et al.~(1999a), as reanalyzed here, data from the Subaru 8
m telescope and HDS by Boesgaard \& Novicki (2006) and data from Keck I with
HIRES for three standard stars by Boesgaard (2007).

\section{Data Analysis}

The data were analyzed using IRAF and MOOG, a stellar synthesis and analysis
program\footnote{http://verdi.as.utexas.edu/moog.html} (Sneden 1973).  We
used the 2002 version, which includes the UV opacity edges of the major atomic
species. Equivalent widths of spectral lines were measured using IRAF
routines.  MOOG was used for two purposes: stellar parameter determination and
abundance determinations.  MOOG's ``abfind'' driver was used with equivalent
widths of Fe I, Fe II, Ti I, and Ti II lines to measure the metallicity,
[Fe/H], and determine $T_{\rm eff}$ and the surface gravity, log $g$, of the
stars in our sample.  We used the ``synth'' driver to create synthetic spectra
in order to measure Be and O abundances.

\subsection{Stellar Parameters}
Stellar parameters are needed to generate the stellar models used for spectral
synthesis and abundance determination.  We have chosen to determine our
stellar parameters spectroscopically.  We are dealing with spectral lines and
need to know the temperature in the region of the atmosphere where the lines
are formed.  These temperatures seem more relevant for abundance work than the
continuum-based color temperatures.  It is generally the case that
spectroscopic temperatures are lower than color temperatures because the
continuum is formed deeper in the stellar atmosphere where the temperatures
are higher.

Several Fe I and Fe II and Ti I and Ti II lines fall within the range of the
HIRES CCD and can be used to determine [Fe/H], $T_{\rm eff}$ and log $g$,
using an iterative method similar to that of Stephens (1999).
Microturbulence, $\xi$, has a negligible effect on any of the measurements in
this project, so a standard value of 1.5 km/s was used in all models (Magain
1984).  We measured equivalent widths of 30-60 Fe I and 10 Fe II lines on the
green and red chips of the HIRES CCD as well as 3 Fe I lines from the UV chip
near the Be doublet.  The Fe I lines covered a range in excitation potential
from 0.86 to 3.98 eV and so could be used to derive the temperature.  In
addition we measured equivalent widths of 10-15 Ti I lines and 8-12 Ti II
lines.  A list of the lines measured, along with the excitation potential and
the $gf$ values, are given in Table 2.  The values for $gf$ come from the
compilation presented in Stephens (1999) and Stephens \& Boesgaard (2002)
which gives the appropriate source references.

We have been careful to use only weak lines which would be on the linear
portion of the curve of growth.  This corresponds to log (W/$\lambda$) $<$
$-$4.82 which is $\sim$75 m\AA\ at 5000 \AA.  For some stars we applied a even
stricter limit of log (W/$\lambda$) $<$ $-$5.15 which is $\sim$35 m\AA\ at
5000 \AA.  The measured equivalent widths are given in the tables in the
Appendix.  For two stars, G 64-12 and G 64-37, we adopted the parameters of
Stephens \& Boesgaard (2002) which were also determined spectroscopically;
they used the same method as we did here, but they had a larger wavelength
coverage toward the red and thus more useable lines.

Using an initial estimate for the stellar parameters and the measured
equivalent widths, MOOG's ``abfind'' driver was then used to calculate an Fe
abundance for each line, as well as an average [Fe/H] from the equivalent
widths.  MOOG also calculates the slope of each line's calculated abundance
versus its excitation potential (EP).  If the temperature is correct, there
should be no trend of abundance with EP.

Atomic abundances derived from two ions of the same species should give
similar results.  If the abundances do not agree, the value of log $g$ can be
adjusted and the abundances recalculated with MOOG.  For the purposes of this
project, MOOG was run iteratively by first calculating the abundance from Ti I
lines using the temperature found via the Fe I lines, then finding the
abundance from Ti II and comparing the two sets of abundances.  The value for
log $g$ is adjusted and the process can then be run iteratively between Fe I
and Ti I plus Ti II abundances until [Fe/H], $T_{\rm eff}$ and log $g$ are
calculated.  (We have used the ionization balance between Ti I and Ti II,
rather than Fe I and Fe II out of concern about NLTE effects on Fe I lines in
metal-poor stars, e.g.~Th\'evenin \& Idiart 1999.)  The stellar parameters
used to calculate abundances are shown in the second, third and fourth columns
of Table 3.  Typical errors are $T_{\rm eff}$ $\pm$80 K, log $g$ $\pm$ 0.2
dex, [Fe/H] $\pm$0.15 dex.  Although our random errors on the temperatures
determined spectroscopically are typically $\pm$80 K, we point out that there
may be large (100-250 K) systematic differences in temperature determinations
by other methods such as UBV and other photometric indices (and the potential
uncertainties in the reddening corrections), the infrared flux method,
Balmer-line profiles.

Although we have relied on the Magain (1984) value of 1.5 km s$^{-1}$ for the
microturbulence, we have checked the effect of using 1.0 km s$^{-1}$ on the
stellar parameters.  On average with $\xi$ = 1.0 km s$^{-1}$, the value of
$T_{\rm eff}$ is lower by 66 K, log $g$ is lower by 0.15 and [Fe/H] is lower
by 0.03 dex.  All of these are within the range of the uncertainities in our
parameter determinations.

Ten of our stars have been observed for Li by Hosford et al.~(2009).  They use
the excitation balance of Fe I to constrain their temperature scale and they
constrain log $g$ using theoretical isochrones.  They present two sets of
parameters for most of the stars appropriate for a main-sequence star and a
subgiant-branch star.  Our temperatures for the 10 stars in common are
typically lower than theirs by $-$163 K on average.  Our log $g$ values are
usually more in alignment with their subgiant branch gravities and on average
are lower by $-$0.05, well within the uncertainty of the determination.  While
the Li abundances are sensitive to temperature, but not to log $g$, the Be
abundances are quite insensitive to temperature, but are sensitive to log $g$
(see $\S$3.2 below.)

\subsection{Abundance Determination and Errors}

The difficulty in determining A(Be) in metal poor stars can be seen in Figure
1, which shows a comparison of a very metal poor star, BD +3 740 at [Fe/H] =
$-$2.95 and one of much higher metallicity, HD 194598 at [Fe/H] = $-$1.23,
from Boesgaard (2007).  Both Be II lines are considerably weaker in the low
metallicity star, BD + 3 740, and high signal-to-noise ratios are needed in
spectra from such stars.  The blending lines are also weaker, however, and the
continuum placement is easier in the metal poor stars.

Once stellar parameters were calculated for each star, they were used to
interpolate models from the Kurucz (1993) grid.  The models, coupled with a
line list (adapted from Boesgaard et al.~1999a) generate a synthetic spectrum,
which can then be compared to the data.  The data are shifted in wavelength
space to match the model and gaussian smoothing is applied to the model
comparable to the shape of the line.  Residuals of the synthetic spectrum
versus the data can be viewed in MOOG in order to determine appropriate values
for the wavelength shift and gaussian smoothing.  The synthesis of the Be
region of two stars is shown in Figure 2.  Both Be and O abundances are
optimized for the best fit in Figure 2.

The Be abundances determined this way are listed in Table 3 for each star.

Inasmuch as Li is more fragile with respect to nuclear destruction than Be, we
looked at Li abundance determinations to ascertain whether there might be Be
deficiencies.  All of our stars except one (LP 752-17) have been observed for
the abundance of Li, and all have normal Li-plateau abundances.  Based on the
criterion that normal Li-plateau values imply normal Be abundances, we
conclude that none of these 23 stars has undergone Be depletion or dilution.
Eleven of our stars appear in the Li compilation of Ryan et al.~(1996) and
nine are in the study of Ryan et al.~(1999).  The mean A(Li) = 2.17 $\pm$0.12.

Although O is one of the most abundant elements in the universe, the
determination of the O content in stars is not straight-forward.  The three
spectral features most commonly used are the OH lines in the UV, the O I
triplet near 7774 \AA\, and the [O I] line at 6300 \AA.  All three features
present drawbacks, see e.g. Boesgaard (2001), Asplund \& Garc\'\i a P\'erez
(2001).  In this work we are constrained to use the OH features in the UV
which are near the Be II resonance lines.  Our spectra do not extend past 6000
\AA.  Furthermore, the [O I] line is very weak in metal-poor stars, and in
dwarfs and subdwarfs in particular; the O abundances found from [O I] for red
giants may reflect altered O abundances due to nuclear processing of O along
with mixing to the stellar surface.

There are electronic transition lines of OH that fall on the UV chip which we
used to calculate the O abundance of each star using the ``synth'' driver as
well.  Some of these lines fall within the Be synthesis region, just redward
of the Be line at 3130 \AA; two others region are further toward the red at
3139 \AA\ and 3140 \AA.  The line list for the 3139 region covers 1.5 \AA\ and
contains seven OH lines, the strongest of which is at 3139.17 \AA\ with a log
$gf$ of $-$1.71, as well as atomic lines.  The line list in the 3140 region
covers 1.1 \AA\ and five OH transitions plus atomic lines.  Generally, the O
abundance agreed with the initial abundance calculated in the Be region.  The
average values are listed in Table 2, where the 3130 \AA\ line has 2 times the
weight of the other two lines.  The synthesis fits for the two other OH line
regions that we used are shown in Figure 3 for HD 140283 which has [Fe/H] =
$-$2.56.

Errors in the abundance determinations are a result of error in the stellar
parameters as well as the accuracy of the synthesis fits.  Errors in the
abundances were determined by adjusting the stellar parameters [Fe/H], $T_{\rm
eff}$ and log $g$ to see how it affected the abundance calculation.  The
largest effect on the Be abundances of any of the stellar parameters is due to
uncertainties in log $g$.  For instance, a change in $T_{\rm eff}$ in $\pm$80
K results in a change in A(Be) of $\pm$0.02 and a change in [Fe/H] of $\pm$0.2
dex results in a change in A(Be) of $\pm$0.01.  However, a change in log $g$
of $\pm$0.2 results in a change of about $\pm$0.10 dex in A(Be).  These errors
are added in quadrature, which results in a typical error in A(Be) of about
$\pm$0.12 dex.  The individual errors are given in Table 3.

The errors on O abundances also result from the uncertainties in the stellar
parameters and were determined in a similar manner.  For the OH features,
however, the error in $T_{\rm eff}$ has the largest effect on the determined
abundance, followed by log $g$.  A change in $T_{\rm eff}$ of $\pm$80 K
changes [O/H] by $\pm$0.15 to 0.17 dex, while a change in log $g$ of +0.2 only
results in a change of [O/H] by +0.05 to +0.07 dex.  These uncertainties were
added in quadrature and appear in Table 3.

Although the ``missing'' UV opacity has been blamed for inaccuracies in the Be
abundances (Balachandran \& Bell 1998), in these stars of such low
metallicity, there is no problem of continuum placement and no issue of
missing opacity.  Similarly, the O abundances from the UV lines of OH are
basically unaffected by opacity or continuum placement issues.  We point out
again that the 2002 version of MOOG includes the UV opacity edges of the
dominant elements.  In considering this issue Smiljanic et al.~(2009)
calculated the effect of an increase of 0.20 in [Fe/H] (a factor of 1.6 as
suggested by Balachandran \& Bell 1998) on Be for a star with [Fe/H] = $-$0.5;
they found that the difference in A(Be) was a negligible effect of only 0.022
dex.

Another source of uncertainty in the O abundances found from the OH features
is in our use of standard LTE 1-D hydrostatic model atmospheres.  Asplund \&
Garc\'\i a P\'erez (2001) use 3-D hydrodynamical model atmospheres to assess
the effect of their more sophisticated models on the O abundance.  They find
that the corrections to [O/Fe] increase with decreasing [Fe/H] and also
increase with decreasing gravity ranging from 0.3 to 0.4 dex at [Fe/H] =
$-$2.0 to 0.5 to 0.6 dex at [Fe/H] = $-$3.0.  They suggest that some of the
calculated decrease in O may be mitigated by a full calculation of non-LTE
effects.

We note that the absolute solar O abundance does not affect our [O/H] values
as the MOOG input value is just subtracted out.  The same is true for the Be
abundance; it does not matter whether we use the meteoritic Be abundance or
some solar abundance.

\subsection{New Be and O abundances from Previous Observations}

We have reanalyzed the Keck/HIRES observations in Boesgaard et al.~(1999a,
1999b) to use the newer version of MOOG from 2002.  These observations were
made with the original HIRES and due to the limited wavelength coverage of
that CCD, we could not determine the stellar parameters spectroscopically.
Instead we used the the parameters from the 1999 papers with temperatures on
the Carney (1983) scale.  That temperature scale is consistent with the
spectroscopic temperatures derived for the stars in this paper.  For three
stars we have lowered the temperature and achieve better fits.  The new Be
abundances were found by the same method used here - spectrum synthesis - and
with the same line list.  We have rederived the O abundances but use only
three OH features (3130, 3139, and 3140 \AA), as was done here for the newly
observed stars.  The new abundance results and the 1$\sigma$ error estimates
are presented here in Table 4 along with the stellar parameters and the
1$\sigma$ errors on those.

The Be abundances do not differ much from those presented in Boesgaard et
al.~(1999a): on average they differ by $-$0.005 $\pm$0.077 with a range in
$\Delta$A(Be) of $-$0.14 to +0.15 (in the sense of new minus old).  The O
abundances differ systematically when we use just the three clean OH features.
The difference in $\Delta$[O/H] ranges from $-$0.55 to +0.12 and the average
is $-$0.22 $\pm$0.16 (again in the sense of new minus old).  The differences
result from our use of the three best OH features and ignoring the other weak
and blended OH lines used in the earlier study.

\subsection{Be and New O Abundances for Published Metal-Poor Stars}

We have incorporated in this study the results on Be in metal-poor stars
published by Boesgaard \& Novicki (2006) (four stars), and the three Li-Be
normal stars from Boesgaard (2007).  These stars, their parameters and
abundances are given in Table 4.  We have determined O abundances in these
stars in the same manner as above, from 3 OH features.  Those results are
presented in Table 5.

\subsection{Be Abundance Comparisons}

There are two stars in common in this study with Boesgaard et al.~(1999a): BD
+3$^{\arcdeg}$ 740 and BD $-$13$^{\arcdeg}$ 3442.  For BD +3$^{\arcdeg}$ 740
the parameters determined spectroscopically and photometrically agree well.
The temperatures differ by $-$80 K, log $g$ by 0.21, [Fe/H] by $-$0.06, [O/H]
by $-$0.10, and A(Be) by $-$0.03 (in the sense of this study minus the
reanalyzed (Table 3) values from the data from Boesgaard et al.~1999a).  In
the case of BD $-$13$^{\arcdeg}$ 3442 the spectroscopic gravity determined
here is considerably higher (4.11 versus 3.50) and thus A(Be) is lower by 0.6
dex.  The low value for log $g$ from Boesgaard et al.~(1999a) is very
uncertain ($\pm$0.44 in Table 3) so we adopt the results from this study.  In
addition, the S/N of the spectrum used here is 183 compared to the earlier one
at 129.

There are two stars analyzed in this project that were previously considered
by Primas et al.~(2000a, 2000b), LP 815-43 and G 64-12.  The values of A(Be)
calculated from our data are slightly different from theirs: for LP 815-43 we
derived A(Be)= $-$0.95, while Primas et al.~(2000a) found A(Be)= $-$1.09 and
for G 64-12 we found A(Be)= $-$1.43, whereas Primas et al.~(2000b) derived
A(Be)= $-$1.15.  Primas et al.~(2000a, 2000b) adopted different stellar
parameters for their stellar models, which would lead to different abundances.
The parameters used in this project were determined spectroscopically and fit
the data well.  When the stellar parameters from Primas et al.~(2000a) were
applied to our data for LP 815-43, A(Be) was found to be about $-$1.10, but
other spectroscopic features in the Be region did not fit the data as well as
with their parameters as with ours.

Although there are no stars from Smiljanic et al.~(2009) in common with our
new Keck data set in Table 3, there are 11 in common with our revised
Boesgaard et al.~(1999a) values in Table 4.  Even though they used log $g$
values derived from parallax and ours are from ionization balance, with the
exception of one star (HD 219617) the agreement in log $g$ is good with a mean
difference of +0.01 $\pm$0.20 (in the sense of this study minus Smiljanic et
al.).  The differences in A(Be) range from $-$0.17 to +0.19 dex (in the sense
of this study minus Smiljanic et al.).  The mean difference is +0.02 $\pm$0.13
dex, indicating agreement within the individual errors.  Two of the seven
stars in Table 4 were also studied by Smiljanic et al.  For HD 194598 our
value for A(Be) is +0.05 dex higher and for HD 195633 our value is +0.18 dex
higher.

There are four stars in common from Tables 4 and 5 with Tan et al.~(2009): HD
76932, HD 132475, HD 140283, and HD 195633.  The mean difference in our values
minus theirs for A(Be) is $-$0.03 $\pm$0.16 dex with the largest difference
($-$0.24 dex) being that for HD 140283 caused by stellar parameter
differences.  

Interestingly, there are nine stars in common to the studies of Tan et
al.~(2009) and Smiljanic et al.~(2009) from the same ESO/UVES data set.  Their
agreement for A(Be) is good with a mean difference (Smiljanic minus Tan) of
+0.07 $\pm$0.09 dex.

\section{Results and Discussion}
In this paper we discuss the abundance relationships between Be and Fe,
between Be and O, and between Fe and O.  We have acquired additional
observations of Be in 34 more stars with [Fe/H] between $-$0.7 and $-$2.5.  A
paper by Boesgaard et al.~(in preparation) will address the issue of the
possibility of a spread in Be at a given [Fe/H] and [O/H] in the larger sample
of 82 stars.  That larger sample will also examine the potential connection
between A(Be) and [O/Fe] and the use of Be as a chronometer.  Preliminary
results from a smaller sample on these matters have appeared in a conference
proceeding by Boesgaard et al.~(2008).

\subsection{Be and Fe Relationships}
The Be abundances from Tables 3, 4, and 5 are shown in Figure 4 plotted
against the values of [Fe/H].  The abundances show a trend of increasing A(Be)
with increasing [Fe/H].  

The straight line between A(Be) and [Fe/H] is

A(Be) = 0.92 ($\pm$0.04)[Fe/H] + 1.41 ($\pm$0.09).

One of our stars in Table 3, LP 831-70, has only an upper limit on A(Be); it
is plotted as an inverted triangle in the figures and it has not been included
in the linear fits.  The slope of the least squares fit is 0.92 $\pm$0.04.
This slope is similar to that found by Boesgaard et al.~(1999a) of 0.96
$\pm$0.04.  It is in agreement with the value of Molaro et al.~(1997) of 1.07
$\pm$0.07, but less steep than the recent values of 1.16 to 1.27 $\pm$0.07
found for the samples in Smiljanic et al.~(2009).  Tan et al.~(2009) augmented
their sample with the results for A(Be) and [Fe/H] from Boesgaard et
al.~(1999a) in their Figure 7; considering their sample of 21 stars (omitting
four ``abnormal'' Be values) they derive a slope of 1.10 $\pm$0.07.

In Figure 5 we show two fits: one with a slope of 0.56 $\pm$0.07 for the stars
with [Fe/H] $<$$-$2.2 and one of 1.07 $\pm$0.06 for stars with [Fe/H]
$>$$-$2.2.  The slope for the higher [Fe/H] stars is in good agreement with
those reported in the earlier studies, most of which contain only stars with
[Fe/H] $>$ $-$2.5.  The shallower slope for the metal-poorest stars may
indicate a slower, but gradual increase in Be as Fe increases in halo stars.
On the other hand, the fit with one line is statistically no better or worse
than the one with two lines.  We have also tried making the break at [Fe/H] =
$-$2.7; the fit is less good and the slopes are nearly the same at 0.57
$\pm$0.12 and 1.05 $\pm$0.05.  The two slopes intersect at $-$2.3, not $-$2.7.

The possible nucleosynthesis of Be in the Big Bang has been discussed by
several authors (e.g.~Boyd et al.~1989, Malaney and Fowler 1989, Alibes et
al.~2002, Jedamzik and Rehm 2006).  Only in some inhomogeneous BB cosmologies
is there enough Be produced to observe using current instruments.  Our results
do not seem to indicate a Be plateau at very low metallicities, as has been
suggested by Primas et al.~(2000a).  Analyzing stars of even lower
metallicities might provide further information, but observations become very
time consuming, due to the faintness of the extremely low metal stars and the
need for high S/N in order to observe the increasingly weak Be II lines.  Our
results in Figures 4 and 5 show no evidence for a plateau of constant A(Be),
but a slower increase in Be with Fe at very low metallicities.

A more sensitive way to search for evidence of a plateau is to examine the Be
abundance normalized to Fe, shown in Figure 6.  The four stars with [Fe/H] $<$
$-$3.0 do show enhanced Be relative to Fe.  The best fit relationship is

[Be/Fe] = $-$0.08 ($\pm$0.04) [Fe/H] $-$0.01 ($\pm$0.09).

The dashed line in the figure corresponds to the mean value of [Be/Fe] for the
48 stars which is +0.17 $\pm$0.21.  This probable error bar is shown at the
left in Figure 6 along the mean value line for [Be/Fe].  Although there is
considerable scatter in [Be/Fe] (the range is $-$0.21 to +0.65), the stars
with the highest values of [Be/Fe] are generally the lowest metallicity stars.
But, as in Figures 4 and 5, there is not a plateau with a constant A(Be) or
[Be/Fe], but rather a smooth trend.  The removal of the two Be-rich stars, HD
94028 and HD 132475 near [Fe/H] of $-$1.5, reduces the mean only to +0.16 but
shows a clearer trend of high [Be/Fe] at low [Fe/H].  (We note that Smiljanic
et al.~2009 found a lower log $g$ for HD 94028 and thus a lower Be abundance,
putting it in alignment with the other stars.)


\subsection{Be and O Relationships}

The nucleosynthesis of Be is directly related to O as the spallation of
abundant atoms like C, N, and O produce smaller atoms like Li, Be, and B.  The
spallation can occur 1) as high energy cosmic rays bombard CNO atoms in the
ambient interstellar gas or 2) during supernovae explosions where atoms of CNO
and protons are accelerated to high energies in the vicinity.  These two
mechanisms predict different relationships between Be and O.  

In the more traditional GCR spallation the slope between A(Be) and [O/H] would
be 2 because the number of O atoms is proportional to the cumulative number of
SN IIa ($N$) while the energetic protons are proportional to the instantaneous
number of SN II ($dN$).  The integral of $N dN$ is $kN^2$.

In the immediate vicinity of SN II the number of Be atoms would be
proportional to the number of O atoms; this would result in a slope of 1
between A(Be) and [O/H].

In Figure 7 we show the least-squares fit between A(Be) and [O/H].  The
equation for the straight line is

A(Be) = 1.21 ($\pm$0.08) + 1.32 ($\pm$0.14).

There are reasons to suspect that Galactic evolution will play a role in the
rate of the production of SN II.  The early generations of stars will be
massive stars which will produce supernova and begin to enrich the next
generations of stars.  These early SN will expel C and O which has been
produced by nuclear fusion in their interiors prior to the supernova phase.
The dominant mechanism creating Be atoms would thus be \#2 above: acceleration
of CNO from the SN explosion.  Eventually, after multiple generations of stars
have formed, the dominant mechanism would be \#1 above: spallation by high
energy cosmic rays in the interstellar gas.

In Figure 7 one can discern a possible change in slope so we have tried
fitting two lines to the data.  This result is shown in Figure 8.  The oldest
stars with the least O, [O/H] $<$ $-$1.4, have a shallower slope of 0.74
$\pm$0.11, while the stars with [O/H] $>$ $-$1.8 have a steeper slope of 1.59
$\pm$0.15.  The transition occurs near [O/H] of $-$1.6.

For the stars with the lowest [O/H] the relationship is

A(Be) = 0.74 ($\pm$0.11) [O/H] + 0.40 ($\pm$0.22).

For the stars with the higher values of [O/H] the relationship is

A(Be) = 1.59 ($\pm$0.15) [O/H] + 1.75 ($\pm$0.20).

A changing slope is consistent with the change expected as the dominant
production method for the spallation of Be changes.  Other effects enter into
the determination of the relationship between Be and O, such as mass outflow
which may be lost during the process of star formation.  Such effects would
modify the predicted slopes of 1 and 2 to lower values.  While we do not plan
to discuss the models of Galactic chemical evolution, we point to these
observations as a potential constraint on such models.

With respect to the [O/H] abundances in Figures 7 and 8, we are concerned
about the reduction in the O abundances through the use of 3-D hydrodynamic
models as discussed by Asplund \& Garc\'\i a P\'erez (2001).  The decrease in
O becomes larger as the metallicity decreases.  The stars with the lowest
values of A(Be) would have decreased [O/H] values.  The low Be points in those
two figures would move to the left along the x-axis, thus favoring the 2-slope
solution.

\subsection{Fe and O Relationships}

Our data enable us to reenter the controversy regarding the form of the
relation between O and Fe.  Our sample contains 49 stars which are all dwarfs
and subgiants - no giants; the range in log $g$ is 3.20 to 4.52.  The
temperatures range from turn-off stars at 6400 K to dwarf stars at 5500 K.
Our star sample covers 3 orders of magnitude in [Fe/H] from $-$0.5 to $-$3.5
with the bulk of the stars, 34, being below $-$2.0 and six of them $\leq$
$-$3.0.  This is a large uniform collection of similarly-analyzed stars which
pushes the boundary to lower metallicities in an unevolved sample than was
done, for example by Nissen et al.~(2002), Fulbright \& Johnson (2003),
Garc\'\i a P\'erez et al.~(2006).  Although the Cayrel et al.~(2004) study
does go to low values of [Fe/H], their stars are low-gravity evolved stars
with log $g$ ranging from 0.7 to 2.7; they find a mean [O/Fe] of 0.7 with a
dispersion of 0.17 dex and 0.47 after an uncertain correction for stellar
surface inhomogeneities.

In Figure 9 we show the tight, linear relationship between [O/H] and [Fe/H]
with a well-fitting slope of 0.70 $\pm$0.03.  The lowest metallicity stars
merge well into the stars with higher values of [Fe/H].

The relationship is

[O/H] = 0.692 ($\pm$0.032) [Fe/H] $-$ 0.076 ($\pm$0.074)

Figure 10 shows the usual relationship where the element is normalized to Fe,
in this case [O/Fe].  This can be fit well with a linear relationship where
[O/Fe] is $\sim$1 at [Fe/H] = $-$3.5 and progresses towards 0 at [Fe/H] = 0.
The equation is

[O/Fe] = $-$0.30 ($\pm$0.03) [Fe/H] $-$0.05 ($\pm$0.07)

The continued upward trend toward lower metallicities in this relationship is
similar to that found by Israelian et al.~(1998, 2001), Boesgaard et
al.~(1999b), Nissen et al.~(2002) although none of those studies went to as
low values of [Fe/H] as in this work.  Fulbright \& Johnson (2002) show a very
similar trend with their ad hoc temperature which forces the parameters to
yield agreement between the permitted and forbidden O features (see their
Figure 13).  

\subsection{Implications}

The linearity of the relation between [O/H] and [Fe/H] in Figure 9 is
striking.  And it implies that the 2 slope fit between A(Be) and [O/H] (Figure
8) probably requires the 2 slope fit between A(Be) and [Fe/H] (Figure 5).  In
turn this indicates that the upturn in [Be/Fe] in Figure 6 at the lowest
metallicities is significant.  It is not evidence for a ``plateau'' in the Be
abundance, but rather a change in the chemical history of the formation of Be.
As discussed in $\S$4.2 the spallation reactions occurring in the vicinity of
SN II would dominate in the early episodes of Galactic evolution and there
would be a flatter slope between Be and O.  While still continuing to produce
Be by this means, the standard GCR spallation would be more dominant in later
eras, thus a steeper slope between Be and O.  Again, we note that other
stellar formation and evolution effects enter into the predicted relation and
reduce both slopes.

\section{Summary and Conclusion}
We have used the upgraded HIRES on the Keck I 10m telescope to obtain high
resolution, high signal-to-noise spectra of 24 very metal deficient stars from
[Fe/H] = $-$3.45 to $-$2.32.  Of the 24 stars 19 have [Fe/H] $<$$-$2.5 and the
other 5 are near $-$2.3.  (Even the large Be sample of Smiljanic et al.~(2009)
contains only three stars in our metallicity range, and those are at the high
end of our range, with [Fe/H] = $-$2.30, $-$2.35, and $-$2.48 and the latter
has S/N of only 45.)  Our data set was supplemented by spectra of 25 stars
from Boesgaard et al.~(1999a), Boesgaard \& Novicki (2006) and Boesgaard
(2007).  For all these stars Be abundances determined from the Be II doublet at
3130 and 3131 \AA\ and O abundances were found from 3 features of OH in the UV
spectral region.  For the newly observed stars stellar parameters were found
spectroscopically using equivalent widths of Fe I for $T_{\rm eff}$, and Ti I
and Ti II for log g.  An iterative procedure was followed for the parameter
determination and the Fe abundance.  The stellar synthesis program MOOG was
used in both {\it abfind} and {\it synth} modes to find the elemental
abundances of Be, O, and Fe.

The full sample of 49 stars reveals a tight correlation between [O/H] and
[Fe/H] with a slope of 0.69 $\pm$0.03 over 3 orders of magnitude in [Fe/H]
from $-$0.5 to $-$3.5 and over 2 orders of magnitude in [O/H] from $-$0.4 to
$-$2.5.  This is shown in Figure 9.  One of the implications of this linear
fit is that the stellar parameters and the abundance determinations give
consistent results for the two sets of stars: the 24 very metal-poor stars
(with spectroscopically determined parameters) and the reanalysis of the data
from the three earlier papers (with photometrically-derived temperatures).

For one star we could only determine an upper limit on the Be abundance so
our sample size is 48 for the Be data.  Of the 48 stars 34 are below $-$2.0 in
[Fe/H] and five of those are $\leq$$-$3.0.  This sample by far the largest
sample of Be abundances found for low metallicity stars.  In the sample of
Smiljanic et al.~(2009) of 73 Be detections, only four have [Fe/H] $<$ $-$2.0
and four more at $-$2.0.

We derive a linear relationship between A(Be) and [Fe/H] which has a slope of
0.92 $\pm$0.04 and a linear relationship between A(Be) and [O/H] with a slope
of 1.21 $\pm$0.08.  Both plots can also be well fit by two linear
relationships.  The tightness of the relationship between the two abscissae,
[Fe/H] and [O/H], lends support for the two-slope fits with A(Be).  The slope
change between A(Be) and [O/H] from 0.74 $\pm$0.11 at low [O/H] to 1.59 $\pm$
0.15 occurs at [O/H] = $-$1.6.  See Figure 8.  The corresponding change point
occurs at [Fe/H] = $-$2.25.  See Figure 5.  This ``kink'' occurs near $-$0.9
for A(Be).

The change in slope between A(Be) and [O/H] would be a natural consequence of
the change in the dominant mechanism of the formation of Be.  In the early
eras of the Galaxy, when massive stars become supernovae, the Be would be
formed predominantly in the vicinity of the SN II where CNO atoms are
accelerated to high energies and break up as they bombard protons, for
example, and form lighter atoms such as Be.  The slope would be a maximum of
1.0 as we find 0.74 for the stars with the lowest values of [O/H].  Later the
spallation reactions could occur as high energy cosmic rays bombard CNO atoms
in the interstellar gas.  In that case the slope would be a maximum of 2.0,
where we find 1.59.  The change between the two mechanism's dominance occurs
as [O/H] increases from $-$1.8 to $-$1.4.

As we have found, the relationship between [O/H] and [Fe/H] is linear; the
change in the slope between A(Be) and [Fe/H] also results from the change in
the dominant mechanism that forms Be.  Whereas the formation of Be and O are
related, the same cannot be said of Be and Fe.  The Fe atoms are a product of
stellar nucleosynthesis via the e-process and not directly connected to the
formation of Be.  The slope change is due to Be production, not Fe.

One of our goals was to search for a Be plateau at low metallicities,
analogous to the Li plateau.  We searched for this in our plot of [Be/Fe] vs
[Fe/H] in Figure 6.  We did not find evidence for a constant Be abundance at
the lowest values of [Fe/H], but rather an indication of a gradual increase in
[Be/Fe] toward the lowest [Fe/H] values.  As Figure 6 shows, the four stars
with the lowest [Fe/H] have enhanced [Be/Fe] equivalent to the one sigma
error.  These data do not indicate that Be has been formed in an
inhomogeneous-style Big Bang.

We have a large sample of low metal stars in which O and Fe abundances have
been determined consistently so we can use the connection between [O/Fe] and
[Fe/H] to learn about the early overproduction of O (formed in massive stars)
relative to Fe (formed in intermediate-mass stars).  We looked at the
enhancement in [O/Fe] as a function of [Fe/H] in Figure 10.  We find that the
O enhancement decrease steadily from [O/Fe] near $-$1.0 at [Fe/H] near $-$3.5
to 0.0 near [Fe/H] = 0.0.  The slope of the fit between the two parameters is
$-$0.30 $\pm$0.03.  The [O/Fe] enhancement diminishes steadily rather than
showing a constant enhancement below [Fe/H] of $-$2.0 of 0.7 $\pm$0.17 as
found by Cayrel et al.~(2004).

\acknowledgements We are grateful to the various Support Astronomers and
Observing Assistants who helped us over the 11 observing runs for this
project.  We acknowledge support from NSF through grant AST 05-05899 to A.M.B.

\clearpage







