Building the Galactic halo from globular clusters: evidence from chemically unusual red giants

Abstract: We present a spectroscopic search for halo field stars that originally formed in globular clusters. Using moderate-resolution SDSS-III/SEGUE-2 spectra of 561 red giants with typical halo metallicities (-1.8 < [Fe/H] < -1.0), we identify 16 stars, 3% of the sample, with CN and CH bandstrength behavior indicating depleted carbon and enhanced nitrogen abundances relative to the rest of the data set. Since globular clusters are the only environment known in which stars form with this pattern of atypical light-element abundances, we claim that these stars are second-generation globular cluster stars that have been lost to the halo field via normal cluster mass-loss processes. Extrapolating from theoretical models of two-generation globular cluster formation, this result suggests that globular clusters contributed significant numbers of stars to the construction of the Galactic halo: we calculate that a minimum of 17% of the present-day mass of the stellar halo was originally formed in globular clusters. The ratio of CN-strong to CN-normal stars drops with Galactocentric distance, suggesting that the inner-halo population may be the primary repository of these stars.

\section{Introduction}
Light-element abundance inhomogeneities in globular clusters have been a
subject of ongoing study for over thirty years. Early observations of a
bimodal distribution of CN bandstrength in globular cluster red-giant
(RGB) stars, anticorrelated with CH bandstrength (e.g., Norris
\& Cottrell 1979\nocite{NC79}; Suntzeff 1981\nocite{S81}; Norris et
al. 1984\nocite{N84}), were interpreted as evidence of anticorrelated
ranges of carbon and nitrogen abundances in globular cluster stars
(e.g., Bell \& Dickens 1980\nocite{BD80}). Since that time,
observations have been extended to include more clusters, fainter
stars within individual clusters, and additional elemental abundances. We
now know that anticorrelated C-N, O-Na and Mg-Al ranges are present in
stars on the red-giant branches of all globular clusters that have been
thoroughly studied (e.g., Kayser et al. 2008\nocite{KHG08}; Carretta et
al. 2009\nocite{CBG09}; Smolinski et al. 2011b\nocite{SM11}), and
that the C-N anticorrelation is also found among main-sequence stars in
several clusters (e.g., Harbeck et al. 2003\nocite{HSG03}; Briley et al.
2004\nocite{BCS04}; Pancino et al. 2010\nocite{PR10}; Smolinski et
al. 2011b\nocite{SM11}). Light-element abundances in globular cluster
stars range from scaled-Solar, similar to halo field stars of the same
metallicity, to depleted in C, O and Mg and enhanced in N, Na and Al,
with between 30\% \citep{PR10} and 70\% \citep{CBG09} of stars in a
given cluster having an atypical abundance pattern.

The atypical abundance pattern resembles the result of
high-temperature hydrogen burning: the CNO cycle converts
carbon and oxygen into nitrogen, and the NeNa and MgAl cycles increase
the abundances of sodium and aluminium while depleting magnesium. The
presence of this abundance pattern in old, low-mass main-sequence stars, which
have neither the high-temperature fusion zones nor the capacity to
transport material between their cores and surfaces, requires that the
abundance variations be extrinsic. The fact that the abundance
differences between stars with typical and atypical abundance patterns
are equally large before and after first dredge-up, which occurs low on the RGB,
indicates that the abundance differences are present throughout the
stars, and are not merely surface pollution.

As a result, the prevailing model for the origin of primordial
light-element abundance variations in globular clusters is that the
stars with atypical light-element abundances are a nearly coeval
second generation formed from material processed by intermediate- or
high-mass stars in the first generation. There are several possible
sites for high-temperature fusion processing in the first
generation. AGB stars are a common suggestion (e.g., Parmentier et
al. 1999\nocite{PJM99}) because they have relatively slow winds, they
are a site of hot hydrogen burning (e.g., Karakas 2010\nocite{K10}), and
they evolve on timescales of $\sim 10^{8}$ years, quite fast compared to
the lifetime of a globular cluster. Other possible sources for
high-temperature-processed feedback material are rapidly rotating
massive stars \citep{DM07} and massive binary stars undergoing mass
transfer \citep{DMP09}, both of which could return more feedback mass in
a shorter amount of time than AGB stars, although the stars
themselves are less common.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ff1.eps}}
\caption[]{Surface gravity versus effective temperature for our
  initial (small gray points) and trimmed (open gray
  triangles) data sets. Solid black lines are 12 Gyr Dartmouth
  isochrones with (left to right) [Fe/H]$=-1.8$, $-1.4$, and $-1.0$.}
\end{figure}

However, all of these possible feedback sources suffer from what has
come to be known as the ``mass budget problem,'' the mismatch between
the current 1:1 ratio of first- to second-generation stars and the
concept that the second generation was built entirely from winds from
the first generation. Recent theoretical models for globular
cluster formation propose two solutions to the mass budget
problem: massive first generations and significant gas accretion. The
models of \citet{BC07} and \citet{DVD08} both require a first
generation with an initial mass $10$ to $20$ times its present mass,
and have second generations built entirely from AGB winds. In the
model of \citet{DVD08}, the second-generation star formation occurs
near the center of the cluster, and concludes when type Ia supernovae
begin. The Ia SNe drive an overall expansion of the cluster, causing
stars at large radii (primarily first-generation stars) to dissociate
from the cluster, thereby reducing the ratio of first- to second-generation
stars to present-day levels. AGB stars are also the source of chemical
differences between the first and second generation in the model of
\citet{CS11}, but their model involves a significant amount of gas
being accreted from the environment, reducing the mass required
for the first generation of stars. The authors investigate the
mass-accretion rate for proto-globular clusters in a variety of
cosmological environments, and conclude that clusters with initial
masses above $10^{4}{\rm M_{\odot}}$ would have been able to accrete
significant amounts of gas in the early Milky Way. They also calculate
that the minimum mass for sufficient gas accretion would be lower for
globular clusters forming as satellites of dwarf galaxies, a claim that could
be tested with observations of intermediate-age star clusters in the
Large Magellanic Cloud (LMC). Interestingly, the recent
high-resolution study of \citet{MCB11} finds no sign of light-element
anticorrelations among 14 red giant stars in the massive
intermediate-age LMC cluster NGC 1866. The authors suggest that the
minimum mass for self-enrichment in the LMC is on the order of
$10^{5}{\rm M_{\odot}}$. Although there are questions about the
likelihood that early globular clusters would accrete unprocessed gas
from their surroundings with precisely the right metallicity (e.g.,
Martell 2011\nocite{M11}), and the gas dynamics involved in mass loss
and dilution in globular clusters have been the subject of few
advanced numerical studies (e.g., Priestley et al. 2011\nocite{PR11};
D'Ercole et al. 2011\nocite{DDV11}), the \citet{CS11} model provides
an important discussion of the interplay between early cluster
self-enrichment and the galaxy-scale environment.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ff2.eps}}
\caption[]{Two sample spectra from our data set, demonstrating the
  effects of varying [C/Fe] and [N/Fe] at fixed stellar
  parameters. The upper spectrum, of SDSS J1500+1126, is CN-strong
  while the lower spectrum, of SDSS J1524+0550, is CN-normal. Shaded
  regions denote the bandpasses of $S(3839)$ and $S(CH)$: lines angled
  up to the right mark the sideband (more closely spaced lines) and
  science band (more broadly spaced lines) of $S(3839)$, and lines
  angled down to the right mark the sidebands (more closely spaced
  lines) and science band (more broadly spaced lines) of $S(CH)$.}
\end{figure}

These scenarios rely implicitly on the high mass and high density of
early globular clusters to retain and/or accrete gas well enough to
permit a second burst of star formation. It is therefore not surprising
that although stars with second-generation-like light-element abundances
are found in every globular cluster in the Milky Way, they are not found
at all in old open clusters. In studies of individual stars (e.g.,
Jacobson et al. 2008\nocite{JFP08}; Martell \& Smith 2009\nocite{MS09})
and mean abundances (e.g., de Silva et al. 2009\nocite{dSG09}) for old
open clusters, light-element abundance behavior is found to be
distinctly different from what is observed in globular clusters. They
are also quite {\it unlikely} to have formed in the halo field,
but they have recently been found there: \citet{MG10} (hereafter MG10)
searched the SEGUE survey (Sloan Extension for Galactic Understanding
and Exploration, Yanny et al. 2009\nocite{YR09}) for halo giants with
unusually strong UV/blue CN bands and identified $49$ (of roughly
$2000$) stars likely to have low carbon abundances and high nitrogen
abundances. Similarly, \citet{CBG10b} compiled a sample of 144 metal-poor disk,
halo and bulge stars from the literature and identified 2 of those as
Na-rich and likely to have originated in globular clusters. In both
cases, the authors interpret the field stars with second-generation
abundances as stars that formed in globular clusters and were later
transferred to the halo through mass-loss processes such as tidal
stripping by the Galaxy or two-body interactions within the cluster. It
is also claimed in MG10 that their figure of $2.5\%$ of halo stars
having second-generation abundances, seen through the lens of a globular
cluster formation model with strong early mass loss (such as, e.g.,
D'Ercole et al. 2008\nocite{DVD08}), implies that as much as $50\%$ of
halo field stars originally formed in globular clusters. \citet{CBG10b}
conclude that early mass loss from globular clusters is ``a major
building block of the halo'', and may also have contributed
significantly to the formation of the thick disk. This possibility has
been suggested theoretically as well: both \citet{BK08} and \citet{MK10}
find that early mass loss from globular clusters should contribute
significant numbers of stars to the halo field.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ff3.eps}}
\caption[]{CN bandstrength index $S(3839)$ versus absolute $r$
  magnitude for our full data set, divided into $0.1$-dex-wide bins in
[Fe/H]. Maximum [Fe/H] in each bin is given in the upper left corner
of each panel. Solid lines are best fits to the CN-normal stars in
each panel. The slope of the best-fit line and the separation between
CN-normal and CN-strong stars both decline with decreasing metallicity.}
\end{figure}

In this paper we use data from the SEGUE-2 moderate-resolution
spectroscopic survey, available as part of the eighth data release of
the Sloan Digital Sky Survey (Aihara et al. 2011\nocite{DR8};
Eisenstein et al. 2011\nocite{SDSSIII}; York et al. 2000\nocite{Y00}),
to search for CN-strong halo giants like those identified in
MG10. This will allow us to make two important extensions to that
work: first, by expanding the sample of stars surveyed, we can refine
our estimate of ${\rm f_{h}^{2G}}$, the present-day fraction of
CN-strong stars in the halo field; second, we calculate ${\rm f_
    {h}^{GC}}$, the fraction of halo stars with first- or
second-generation chemistry originating in globular clusters, as
predicted by several different two-generation cluster formation
models.

\section{The Data Set}

The SDSS and its extensions have acquired $ugriz$ photometry for
several hundred million stars. SEGUE, one of three sub-surveys that
together formed SDSS-II, extended the $ugriz$ imaging footprint of
SDSS-I \citep{FI96,GC98,Y00,PM03,GS06, SL02,DR1,DR2,DR3,DR4,DR5,DR6,DR7}
by approximately 3500 deg${\rm ^{2}}$, and also obtained $R \simeq 2000$
spectroscopy for approximately $240,000$ stars over a wavelength range
of $3800-9200\hbox{\AA}$. SEGUE included spectra for a collection of
Galactic globular and open clusters, which served as calibrators for the
${\rm T_{eff}}$, ${\rm log(g)}$, and [Fe/H] scales for all stars
observed by SDSS/SEGUE, as processed by the SEGUE Stellar Parameter
Pipeline (SSPP; Lee et al. 2008a; b\nocite{LBS08a}\nocite{LBS08b};
Allende Prieto et al. 2008\nocite{AP08}). The SSPP produces estimates of
${\rm T_{eff}}$, ${\rm log(g)}$, [Fe/H], and radial velocity, along with
the equivalent widths and/or line indices for $85$ atomic and molecular
absorption lines, by processing the calibrated spectra generated by the
standard SDSS spectroscopic reduction pipeline \citep{SL02}. See
\citet{LBS08a} for a detailed discussion of the approaches used by the
SSPP; \citet{SL11} provides details on the most recent updates to this
pipeline, along with additional validations. SEGUE-2 (C. Rockosi et al.,
in prep.) expanded the numbers of Galactic stars with available
low-resolution SDSS spectra by over $120,000$, concentrating in
particular on stars with greater distances in order to better sample the
outer-halo region of the Galaxy.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ff4.eps}}
\caption[]{CN-CH planes for our final data set, divided into $0.1$-dex-wide bins in [Fe/H]. The CH bandstrength index $S(CH)$ is plotted versus the CN bandstrength residual $\delta S(3839)$, and candidate globular cluster stars are drawn as open circles. As in Figure 3, the maximum metallicity in each bin is given in the upper left corner. As is explained in the text, we select as candidate globular cluster stars all those with relatively high $\delta S(3839)$ in the three highest-metallicity bins, and stars with both high $\delta S(3839)$ and $S(CH)$ below the mean $S(CH)$ of the CN-normal group in all other metallicity bins.}
\end{figure}

As in MG10, SEGUE-2 targets were selected from across a large range in
parameter space: we began by selecting all SDSS-III targets from
SEGUE-2 plates with [Fe/H]$\leq -1.0$, ${\rm log(g)}\leq 3.0$,
$(g-r)_{0}\geq 0.2$, and a mean signal-to-noise ratio (SNR) of at
least $20$ per pixel. Errors on the derived quantities and
quality-assurance flags were also required to fall within a certain
range: $\sigma_{\rm [Fe/H]}\leq 0.5$, $\sigma_{\rm log(g)}\leq 0.5$,
$\sigma_{\rm Teff}\leq 200$~K. This initial data set contained 2019
stars. We then identified four stars as carbon-rich using the C$_{2}$
indices defined in MG10, and removed them from the sample. RGB
stars were identified by dividing the data into $0.2$-dex-wide
metallicity bins and making a recursive 3-sigma selection about a
fiducial sequence in the $({\rm log(g)}$, ${\rm T_{eff}})$ plane. Since
our analysis is based on the CN and CH molecular absorption features,
which become quite weak at low metallicity, we also restrict our sample
to stars with [Fe/H] $\geq -1.8$ and a SNR per
pixel between $4000\hbox{\AA}$ and $4100\hbox{\AA}$ (SN$_{\rm blue}$)
above $15$. The RGB, metallicity, and SN$_{\rm blue}$ selections together
reduce our data set to $561$ stars.

\begin{table*}
\caption{Star identifiers, position, photometry and parameters for all
  program stars from this study and from MG10.}
\centering
\begin{tabular}{c c c c c c c c c c c c}
\hline \hline
SDSS ID & Plate & MJD & FiberID & $\alpha$ & $\delta$ & $g_{0}$ &
$(g-r)_{0}$ & [Fe/H] & ${\rm T}_{eff}$ & log(g) & Survey\\
\hline
SDSS J1605+1638 & 3289 & 54910 & 132 & 16:05:52.72 & 16:38:59.21 & 16.04 & 0.663 & -1.07 & 4965 & 2.11 & SEGUE-2\\
SDSS J0842+2216 & 3373 & 54940 & 51 & 08:42:16.23 & 22:16:58.27 & 16.14 & 0.715 & -1.05 & 4913 & 2.06 & SEGUE-2\\
SDSS J1528+0759 & 3308 & 54919 & 529 & 15:28:51.34 & 07:59:08.62 & 17.11 & 0.696 & -1.04 & 4950 & 2.29 & SEGUE-2\\
SDSS J1557+3629 & 3478 & 55008 & 229 & 15:57:48.95 & 36:29:33.04 & 16.77 & 0.735 & -1.03 & 4910 & 2.29 & SEGUE-2\\
SDSS J1011+1810 & 3178 & 54848 & 269 & 10:11:13.50 & 18:10:32.73 & 17.22 & 0.689 & -1.00 & 4934 & 2.37 & SEGUE-2\\
SDSS J2243+1322 & 3128 & 54776 & 104 & 22:43:31.49 & 13:22:33.68 & 17.72 & 0.689 & -1.09 & 4877 & 2.55 & SEGUE-2\\
SDSS J1016+1737 & 3178 & 54848 & 55 & 10:16:42.76 & 17:37:23.51 & 17.44 & 0.619 & -1.05 & 5066 & 2.50 & SEGUE-2\\
SDSS J1206+3021 & 3181 & 54860 & 73 & 12:06:49.60 & 30:21:49.39 & 18.06 & 0.739 & -1.05 & 4983 & 2.53 & SEGUE-2\\
SDSS J1210+1042 & 3214 & 54866 & 14 & 12:10:11.17 & 10:42:28.64 & 16.80 & 0.652 & -1.07 & 5025 & 2.56 & SEGUE-2\\
SDSS J1617+0502 & 3298 & 54924 & 48 & 16:17:42.83 & 05:02:32.70 & 17.54 & 0.663 & -1.06 & 4991 & 2.56 & SEGUE-2\\
\hline
\end{tabular}
\end{table*}

Some stars from SEGUE-1 were intentionally re-observed in SEGUE-2,
as a way to check for consistency in data reduction and analysis between
the two surveys. However, none of the stars in this study were also in
the MG10 data set. Table 1 lists SDSS identifiers, the plug plate
number, MJD, and fiber number of the SEGUE observation (useful for
selecting individual stars from the SDSS Catalog Archive
Server\footnote{http://skyservice.pha.jhu.edu/CasJobs/login.aspx}),
right ascension, declination, photometry, SSPP-derived stellar
parameters, and survey name (SEGUE-1 or SEGUE-2) for each star in the
final trimmed data set and the final data set from MG10. Only the
first ten rows of Table 1 are given in the print version of this
paper, to indicate form and content; a full version is available at the
CDS \footnote{via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5)
  or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/}. Figure 1 shows surface gravity
versus effective temperature for our initial (small gray points) and
trimmed (open gray triangles) data sets, using the ${\rm T_{eff}}$ and
${\rm log(g)}$ values derived from the spectra by the SSPP. The
overplotted solid lines are 12 Gyr Dartmouth isochrones \citep{DC08}
with metallicities of [Fe/H]$=-1.8$, $-1.4$, and $-1.0$, spanning the
metallicity range of our data set.

\section{Analysis}

Because we are searching for old, low-mass red giant stars with
globular-cluster-like abundance patterns, our analysis follows a
process typical in the globular cluster literature. Specifically, we are
using index-based techniques developed for globular cluster studies as
indicators of carbon and nitrogen abundance, and we infer the
presence of the full second-generation light-element abundance pattern
from those indices. Overall, we followed the procedure from
MG10 fairly closely in analyzing our data, in order for the results to
be as compatible as possible.

\subsection{Bandstrength indices}

We use the bandstrength indices $S(3839)$ \citep{N81} and $S(CH)$
\citep{MSB08}, which were initially developed for studies of
light-element abundance variations in globular clusters, to identify
stars with strong CN bands and weak CH bands relative to the majority of
the halo field. Indices measure the magnitude difference between the
integrated flux in the feature (the ``science band'') and one or two
nearby, independent regions of the spectrum (the ``sidebands''), in the
sense that more absorption in the feature results in a larger index
value. Since molecular abundance is strongly controlled by the abundance
of the minority species, the CH band is indicative of [C/Fe] while the
CN band traces [N/Fe]. We are using carbon and nitrogen abundances as an
indicator of the ``second-generation'' abundance pattern of depleted
carbon, oxygen, and magnesium and enhanced nitrogen, sodium, and aluminum
because that full abundance pattern is consistently found together in
Galactic globular clusters.

Figure 2 shows typical spectra from our data set, chosen to illustrate
the range of variation in CN and CH bandstrength. The upper spectrum
(of SDSS J1500+1126) is of a CN-strong star with ${\rm T_{eff}}=5250$,
[Fe/H]$=-1.70$, $S(3839)=0.005$ and $S(CH)=0.70$, and the lower
spectrum (of SDSS J1524+0550) has similar stellar parameters but
relatively weaker CN and stronger CH bands (${\rm T_{eff}}=5210$,
[Fe/H]$=-1.74$, $S(3839)=-0.189$ and $S(CH)=0.79$). The spectra are
both normalized near $4500\hbox{\AA}$ and offset vertically. The
shaded areas with lines tipping up to the right mark the sideband
(more closely spaced lines) and science band (more broadly spaced lines) of
$S(3839)$, and the shaded areas with lines angled down to the right
mark the sidebands (more closely spaced lines) and science band (more
broadly spaced lines) of $S(CH)$. The bandpasses of both indices are
given in Table 2. Errors on the bandstrengths were estimated as in
MG10, using a Monte Carlo technique to randomly sample the SEGUE noise
vectors and add that noise to the spectra, then measure the standard
deviation of the bandstrength over 100 realizations. Typical values
for $\sigma_{S(3839)}$ and $\sigma_{S(CH)}$ are $0.015$ and $0.007$,
respectively.

\subsection{Identification of CN-Strong candidates}

We identify candidate second-generation halo stars in our data set in
two stages, first by selecting those stars with unusually high $S(3839)$
relative to other stars of similar metallicity and evolutionary stage,
and then by taking just the subset of those CN-strong stars with
appropriate $S(CH)$. Absolute magnitudes for our stars were calculated
using a spectroscopic parallax method similar to that used in MG10. We
used 12-Gyr Dartmouth isochrones \citep{DC08}, interpolated to the
metallicities of our target stars, then use dereddened $(g-r)_{0}$
colors to predict absolute ${\rm M}_{r}$ magnitudes.

Figure 3 shows $S(3839)$ versus ${\rm M}_{r}$ for our 561 program stars,
divided into $0.1$-dex-wide metallicity bins. The solid lines
overplotted are linear fits to the CN-normal stars, and the maximum
metallicity in each bin is given in the upper right corner of each
panel. The distinctive characteristics of each panel are the dominance
of the CN-normal group, the slope of the best-fit line, and the
separation of the CN-normal and CN-strong groups. In globular
clusters, the fraction of CN-strong stars is roughly $50\%$ \citep{K94}, but the
slope of the CN bandstrength-luminosity relation and the separation in CN
bandstrength between the two groups behaves similarly to what we
observe in field stars. As in MG10, the slope and the separation both
decline with decreasing metallicity, and the two CN bandstrength
groups become difficult to distinguish at the lowest metallicites.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ff5.eps}}
\caption[]{Heliocentric distance as a function of absolute $r$ magnitude for our final data set (upper panel) and for the MG10 data set (lower panel). CN-normal stars are plotted as small dots, while CN-strong stars are shown as filled stars. Only a randomly selected $25\%$ of the CN-normal MG10 stars is shown to reduce crowding. In both panels, the lower boundary is set by the bright limit of the SEGUE survey ($g=14$), and the upper boundary is set by our requirement that ${\rm SN_{blue}}>15$.}
\end{figure}

We measure the quantity $\delta S(3839)$, the difference between an
individual star's CN bandstrength and the best-fit line at the same
absolute magnitude, in the plane of Figure 3. The eight panels of
Figure 4 correspond to the panels of Figure 3, and show $\delta
S(3839)$ versus $S(CH)$ for each metallicity bin, with candidate
second-generation field stars plotted as open circles and all other stars
plotted as filled circles. Because the $S(CH)$ index has a reduced
sensitivity to carbon abundance at the upper end of our metallicity
range, we accept all CN-strong stars with metallicities above
[Fe/H]$=-1.3$ as candidates, and we select only those CN-strong stars
with $S(CH)$ below the mean of the CN-normal group for stars with
metallicity below [Fe/H]$=-1.3$. These are the same criteria as were
used in MG10, and they return $16$ candidate globular cluster stars,
$3\%$ of our sample. This is very similar to the MG10 result based on
SEGUE-I stars.

\subsection{Stellar distances}

We calculate heliocentric distances D to our target stars from the
distance modulus $(r_{0}-{\rm M}_{r})$. Errors in distance were
calculated by randomly sampling the error on $(g-r)_{0}$, then
recalculating the distance. The standard deviation in $100$ of these
realizations was then adopted as $\sigma_{\rm D}$. One of the goals of
the SEGUE-2 survey was to observe stars at greater distances than were
observed in SEGUE-1, and while our data set covers a heliocentric
distance range very similar to that of MG10, the average star at any
fixed ${\rm M}_{r}$ is more distant for the current data set than for
the MG10 data set. The upper panel of Figure 5 shows our calculated
distances versus absolute ${\rm M}_{r}$ magnitudes for the SEGUE-2 data
set, with CN-normal stars drawn as small dots and CN-strong stars drawn
as filled stars. The lower panel shows the analogous data for the MG10
data set, with only a randomly selected $25\%$ of the MG10 CN-normal
stars plotted for visual clarity. In both panels, the stars fill a band
that is restricted at small distances by the SEGUE bright limit of
$g=14$ and at large distances by our requirement that SN$_{\rm blue}
\geq 15.0$.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ff6.eps}}
\caption[]{CN bandstrength index $\delta S(3839)$ versus
  Galactocentric distance ${\rm r_{GC}}$ for our final data set (upper
  panel) and the MG10 sample (lower panel). Symbols are the same as in
  Figure 5. There is a sharp drop in the number and frequency of
  CN-strong stars at roughly $20$ kpc. }
\end{figure}

Heliocentric distances are converted geometrically to Galactocentric
distances ${\rm r_{GC}}$, using the relation ${\rm
r_{GC}^{2}=D^{2}+(8\hbox{kpc})^{2}-2D\cos(\gamma)}$, in which $\gamma$
is the angle on the sky between the Galactic center and the star in
question. Table 3 lists plate, MJD, and FiberID identifiers, along
with our measurements of $S(3839)$ and $S(CH)$, our determination of a
star's CN class (normal or strong), absolute ${\rm M}_{r}$ magnitude,
heliocentric distance, Galactocentric distance, and survey name (SEGUE-1
or SEGUE-2), for the same stars as in Table 1. As with Table 1, only
the first ten rows are given in the print version of this paper, and a
full version of the table is available from CDS. Figure 6 shows
$\delta S(3839)$ versus ${\rm r_{GC}}$ for our final data set (upper
panel) and the MG10 sample (lower panel), using the same symbols as in
Figure 5. The frequency of CN-strong stars appears to drop at roughly
$20$ kpc, with only one CN-strong star beyond ${\rm r_{GC}}=30$ kpc.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ff7.eps}}
\caption[]{The ratio of CN-strong to CN-normal stars, ${\rm f}_{s}$, as a function of Galactocentric distance for the present data set and the MG10 data combined. Vertical error bars represent the Poisson errors on ${\rm f}_{s}$ in each distance bin. The drop in the frequency of CN-strong stars seen in Figure 5 at $20$ kpc is also seen here.}
\end{figure}

\begin{table}
\caption{Bandstrength index definitions}
\centering
\begin{tabular}{l c c c}
\hline \hline
Index & Blue sideband & Science band & Red sideband\\
\hline
$S(3839)$ & - & $3846$-$3883\hbox{\AA}$ & $3883$-$3916\hbox{\AA}$\\
$S(CH)$ & $4050$-$4100\hbox{\AA}$ & $4280$-$4320\hbox{\AA}$ & $4330$-$4350\hbox{\AA}$\\
\hline
\end{tabular}
\end{table}

To further investigate this result, we visualize the relationship
between the frequency of CN-strong stars and Galactocentric distance
in two ways: Figure 7 shows ${\rm f}_{s}$, the ratio of CN-strong to
CN-normal stars, versus Galactocentric distance for the present data
set combined with MG10. Vertical error bars represent Poisson errors in
${\rm f}_{s}$. The fraction of CN-strong stars is lower at larger ${\rm
r}_{GC}$, with a transition at $\simeq 20$ kpc: relatively small-${\rm
r}_{GC}$ stars are consistent with a value for ${\rm f}_{s}$ around
$0.04$ and relatively large-${\rm r}_{GC}$ stars are consistent with a
value for ${\rm f}_{s}$ around $0.01$. Figure 8 shows cumulative
distribution functions of $\delta S(3839)$ for stars with ${\rm r_{GC}}
< 20$ kpc (solid line) and ${\rm r_{GC}} \geq 20$ kpc (dotted line). The
cumulative distribution curve for stars with smaller Galactocentric
distances climbs continuously for $\delta S(3839) > 0.05$, while the
curve for more distant stars is flatter in that range, indicating that
there are relatively more CN-strong stars at smaller Galactocentric
distances. A one-sided Kolmogorov-Smirnov test returns a probability of
$0.15$ that the $\delta S(3839)$ distributions in these two distance
ranges were drawn from the same parent population. This result is
only a very marginal rejection of the single population hypothesis;
however, it clearly points to a reasonable possibility that two parent
populations with differing CN-strong fractions exist.

\begin{table*}
\caption{Measured and derived parameters for all program stars from
  this study and from MG10.}
\centering
\begin{tabular}{c c c c c c c c c c}
\hline \hline
Plate & MJD & FiberID & $S(3839)$ & $S(CH)$ & CN class & ${\rm M}_{r}$
& D (kpc) & R$_{GC}$ (kpc) & Survey\\
\hline
3289 & 54910 & 132 & -0.026 & 0.716 & Normal & 0.31 & 10.3 & 17.7 & SEGUE-2\\
3373 & 54940 & 51 & 0.004 & 0.760 & Normal & -0.12 & 12.9 & 17.2 & SEGUE-2\\
3308 & 54919 & 529 & -0.008 & 0.783 & Normal & 0.06 & 18.6 & 25.2 & SEGUE-2\\
3478 & 55008 & 229 & 0.014 & 0.754 & Normal & -0.23 & 18.0 & 25.5 & SEGUE-2\\
3178 & 54848 & 269 & 0.035 & 0.763 & Normal & 0.22 & 18.2 & 22.6 & SEGUE-2\\
3128 & 54776 & 104 & 0.025 & 0.710 & Normal & 0.01 & 25.3 & 30.3 & SEGUE-2\\
3178 & 54848 & 55 & -0.153 & 0.789 & Normal & 0.80 & 15.9 & 20.5 & SEGUE-2\\
3181 & 54860 & 73 & -0.106 & 0.819 & Normal & -0.32 & 33.8 & 38.2 & SEGUE-2\\
3214 & 54866 & 14 & -0.083 & 0.798 & Normal & 0.40 & 14.1 & 19.3 & SEGUE-2\\
3298 & 54924 & 48 & -0.105 & 0.784 & Normal & 0.33 & 20.4 & 27.6 & SEGUE-2\\
\hline
\end{tabular}
\end{table*}

We interpret this drop-off in the frequency of CN-strong stars with Galactocentric distance as a
possible sign of a transition from the inner-halo population to the
outer-halo population (as discussed by, e.g., Carollo et al. 2007;
2010\nocite{CB07}\nocite{CB10}). Note that this also corresponds to
the transition zone reported by \citet{dJY10}, based on fits to Hess
diagrams from the SEGUE vertical photometry stripes. Recent
theoretical studies of galaxy formation (e.g., Oser et
al. 2010\nocite{OO10}; Font et al. 2011\nocite{FM11}) predict that
stars formed {\it in situ} are the dominant population at small
Galactocentric radii, while stars accreted during minor mergers are
the majority at larger Galactocentric distances. Studies of abundances
in individual stars in dwarf galaxies do not often report stars with
light-element abundances similar to second-generation globular cluster
stars. In \citet{SVT03} there is only one star out of the 15 surveyed
(Fornax 21) with even mildly elevated [Na/Fe], \citet{SA04} found
fairly low [Na/Fe] abundances in the three stars they observed in Ursa
Minor, and \citet{GS05} find the same result for nine RGB stars in
Sculptor. \citet{KG08} also report one star in Carina (of eight with
measured sodium abundance) with moderate [Na/Fe], and \citet{LH10}
report one star in Fornax (of 70 with measured sodium abundances) with
moderate [Na/Fe], but neither study finds any stars with the level of
sodium abundance enhancement typically found in second-generation
globular cluster stars (e.g., Carretta et al. 2009\nocite{CBG09}).

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{ff8.eps}}
\caption[]{Cumulative distribution functions for $\delta S(3839)$ for
  stars with ${\rm r_{GC}} < 20$ kpc (solid line) and ${\rm
    r_{GC}}\geq 20$ kpc. The fraction of stars rises more quickly at
  large values of $\delta S(3839)$ for stars with smaller
  Galactocentric radii, indicating that CN-strong stars are more
  common in the inner halo.}
\end{figure}

A schematic picture of the
Galaxy in which the inner halo has significant contributions from
globular clusters, while the outer halo is dominated by stars that
formed in dwarf galaxies, fits reasonably well with our observation that
a higher fraction of stars in the inner halo have light-element
abundances similar to those of second-generation globular cluster stars.
Spectroscopic surveys currently underway (such as APOGEE; Allende Prieto
et al. 2008\nocite{APM08}), surveys presently being planned as
companions to the Gaia mission (e.g., Cacciari et al. 2009\nocite{C09}),
and studies based on nearby stars with accurate proper motions and halo
kinematics (such as, e.g., Carollo et al. 2010\nocite{CB10}) will
be important ways to explore this question of inner/outer halo
origins further: all of them expand the effective search space,
either by surveying to larger distances or by focusing on stars with
large apogalactic distances.

\section{Interpretation}

The homogeneity of metallicity in present-day globular clusters seems to
favor massive models of globular cluster formation, while the observed
range in oxygen and sodium abundances (see, e.g., Carretta et al.
2009\nocite{CBG09}) indicates that dilution of feedback material and
incomplete mixing of the gas that makes up the second stellar generation
play a central role in setting final light-element abundances for
individual stars. In order to evaluate the veracity of various
globular cluster formation models, it is necessary to explore their
implications for the formation of the Galactic stellar halo. Massive
models, which require a large fraction of first-generation stars to be
lost to the halo field at early times, will produce a large enhancement
between $f_{h}^{2G}$ and $f_{h}^{GC}$.

The recent paper of \citet{SC11} considers this question from a
theoretical standpoint, using the globular cluster self-enrichment model
of \citet{DM07} to calculate the mass in first-generation stars needed
to produce a significant second generation. In the \citet{DM07} model,
rapidly rotating massive stars provide chemical feedback between the
first and second generations, and \citet{SC11} explore the effects of
varying the slope of the high-mass end of the stellar initial mass
function (IMF) on the resulting second-generation population. Requiring
that the present-day ratio of first- to second-generation stars be 1:2,
the authors calculate that initial globular cluster masses must have
been 8 to 10 times their present-day values for a Salpeter IMF. They
then compare the mass in first-generation stars that was lost from
globular clusters at early times to the mass of the present-day Galactic
stellar halo. Since roughly $35\%$ of the mass in a Salpeter IMF is in
low-mass, long-lived stars, and the present-day Galactic globular
cluster system comprises $2\%$ of the mass of the halo, they estimate
that escaped first-generation globular cluster stars make up $5-8\%$ of
the mass of the halo.

\citet{SC11} then consider the origin of the second-generation stars found
in the halo by MG10 and by \citet{CBG10b}, specifically whether they
can reasonably be understood as having been lost to the halo
during the strong early mass loss that removed roughly $90\%$ of the
first-generation stars or whether they were transferred during later
episodes of globular cluster tidal disruption. Folding in some
assumptions involving the initial cluster mass function (ICMF), they conclude
that all second-generation stars in the halo were lost at early
times. This raises the fraction of cluster mass that was lost in this
early phase, which necessarily increases the initial cluster mass to
15 to 25 times the present-day mass and also raises the fraction of
halo stars that originally formed in globular clusters to as much as
$20\%$.

However, given the observed examples of globular cluster dissolution
in progress, e.g., Palomar 5 (Rockosi et al. 2002\nocite{ROG02};
Odenkirchen et al. 2003\nocite{OGD03}), NGC 5466 \citep{BEI06} and the
GD-1 stellar stream \citep{GD06}, as well as the relative emptiness of
the more inhospitable regions of the \citet{GO97} ``vital diagram'', we
expect cluster dissolution to be the dominant source of halo field stars
with second-generation abundance patterns. Theoretical studies of the
ICMF, early star cluster mass loss, and the evolution of the cluster
mass function (e.g., Baumgardt et al. 2008\nocite{BK08}; Parmentier et
al. 2009\nocite{PG09}; Decressin et al. 2010\nocite{DB10}) all suggest
that far more star clusters were initially formed than have survived to
the present day, especially at the lower end of the ICMF. The study of
\citet{JG10}, which used matched-filter techniques to search for faint
tidal tails around Galactic globular clusters in the SDSS imaging
footprint, found stars with cluster-like photometry outside the tidal
radii of a number of clusters without pronounced tidal tails, indicating
that stars continue to escape from globular clusters even in the absence
of dramatic tidal features. In the limit where cluster dissolution is
the sole source of second-generation field stars, we can estimate the
number of present-day globular clusters that would need to be disrupted
to provide the observed $2.5\%$ of halo stars with second-generation
chemistry, ${\rm N_{d}}$. For a typical present-day 1:1 ratio of first-
to second-generation stars, and assuming a stellar halo mass of
$10^{9}{\rm M_{\odot}}$ (e.g., Freeman \& Bland-Hawthorn
2002\nocite{FB02}) and a typical present-day globular cluster mass of $5
\times 10^{5} {\rm M_{\odot}}$, ${\rm N_{d}} \simeq \frac{2 \times {\rm
f_{h}^{2G}} \times {\rm M_{halo}}}{\langle {\rm M_{gc}} \rangle} \simeq
100$. Since the current number of of globular clusters is on the order
of $150$, this means that the initial population of globular clusters
was significantly larger, but this is not an outrageous claim:
\citet{Mv05} calculate that the ratio between the initial and
present-day number of globular clusters in the Galaxy is $3/2$, whereas
we derive $5/3$.

Whatever the details, it seems clear that first-generation stars lost
from globular clusters at early times are an important contribution to
the construction of the Galactic halo, if two-generation self-enrichment
scenarios for globular cluster formation are indeed the correct model.
The models of \citet{DM07}, \citet{DVD08}, \citet{CBG10b}, and
\citet{VM10} all predict that $90\%$ of the first generation ought to
be lost from globular clusters at early times, while \citet{C11}
estimates that as much as $95\%$ of the first generation ought to be
lost. The mass of first-generation stars lost during this phase by the
${\rm N_{GC}}$ globular clusters in the present-day system can be
expressed as ${\rm M_{1G}^{lost}}={\rm M_{GC~system}}\times \frac{1}{2}
\times (\frac{1}{1-{\rm f_{lost}}}-1)$. Considering the early
contributions of first-generation stars from globular clusters that have
survived to the present day as well as the clusters that have completely
dissolved, we estimate that ${\rm f_{h}^{GC}}$ has a minimum value of
$\frac{\rm M_{1G}^{lost}}{\rm M_{halo}} + (\frac{\rm N_{d}}{\rm N_{GC}}
\times \frac{\rm M_{1G}^{lost}}{\rm M_{halo}}) + {\rm M_{GC~system}} =
2\% \times (\frac{9}{2} + \frac{12}{3}) = 17\%$ of halo field stars,
with both first- and second-generation abundance patterns, originally
formed within globular clusters, plus an additional unknown contribution
of first-generation stars from clusters too low-mass to self-enrich that
have dispersed by the present day.

\section{Future challenges}

Based on the current two-generation models for globular cluster
formation, the early phase of globular cluster formation must have
contributed a significant number of first-generation stars to the
Galactic halo. From the data presented herein, by MG10, and by
\citet{CBG10b}, it is also clear that second-generation globular
cluster stars are a component of the halo field. Further
investigations of the process of globular cluster formation will help
to clarify the role that globular clusters played in the early assembly
of the Galactic halo. There are two specific directions that new
observations could take that would be particularly insightful:
investigating the relationship between the minimum mass for cluster
self-enrichment and the larger galactic environment, and searching for
star clusters still in the process of self-enrichment.

These studies will necessarily require observations of star clusters
outside the Milky Way, as a way to compare cluster formation in
different large-scale environments. Observing very young star clusters
may offer the opportunity to see the process of self-enrichment in
progress. \citet{CBG10b} suggest that clusters in the process of violent
relaxation (at a few $\times 10^{7}$ years old) should have a central
compact cluster, an extended, non-bound halo and outflowing gas. It is
suggested in \citet{C11} that the extended clusters discovered in M31 by
\citet{HT05} are in this state (see also Huxor et al. 2011\nocite{HF11}).
If this is the case, the stars at large radius in the extended clusters
should have almost entirely first-generation chemistry, and strong
radial gradients in light-element abundances, most readily observable in
CN and CH molecular features, should be present.

To explore the question of the minimum cluster mass for self-enrichment,
it would be useful to consider star clusters in lower-mass galaxies.
\citet{CS11} claim that the weaker tidal field of the Large Magellanic
Cloud (LMC) ought to permit the formation of a second stellar generation
at lower cluster mass than in the Milky Way. Spectroscopic observations
of red giant-branch stars in intermediate-age populous clusters in the
LMC would be particularly interesting because recent photometric work
(e.g., Mackey et al. 2008\nocite{MBN08}; Milone et al.
2009\nocite{MBP09}; Goudfrooij et al. 2011\nocite{GP11}) has uncovered
broadened or split main-sequence turnoffs in a large fraction of them.
Since they are much younger than Galactic globular clusters (typically
$1$-$4$ Gyr), age differences on the order of a few hundred million
years are visible at the turnoff, and this age difference is
suggestively similar to the age difference expected between first- and
second-generation stars in Galactic globular clusters. Obtaining spectra
of turnoff stars at the distance of the LMC would be extremely
difficult -- indeed, high-resolution spectroscopic studies of RGB stars in
these clusters are limited to fairly small samples (e.g., Mucciarelli et
al. 2008\nocite{MCO08}). However, a lower-resolution study using a
multiobject spectrograph could collect CN and CH data analogous to the
data presented in this paper, and would allow a search for star-to-star
variations in carbon and nitrogen abundance for a larger data set per
cluster. A modified version of the SSPP has been developed for
application to spectra with resolving power as low as $R \sim 1000$),
and is already being used with a variety of low-resolution data from
non-SDSS sources (e.g., see Li et al. 2010\nocite{LC10}; Humphreys et
al 2011\nocite{HB11}). These data would also permit an investigation
of theoretical claims (e.g., Conroy \& Spergel 2011\nocite{CS11}) that
lower-mass galaxies permit the formation of multiple stellar generations in
lower-mass star clusters, by comparing the CN-CH behavior of
intermediate-age LMC clusters across a range of masses.

\newpage

Astrophysics of Galaxies (astro-ph.GA); Solar and Stellar Astrophysics (astro-ph.SR)


Superconductivity up to 30 K in the vicinity of quantum critical point in BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$

Abstract: We report bulk superconductivity induced by an isovalent doping of phosphorus in BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$. The P-for-As substitution results in shrinkage of lattice, especially for the FeAs block layers. The resistivity anomaly associated with the spin-density-wave (SDW) transition in the undoped compound is gradually suppressed by the P doping. Superconductivity with the maximum $T_c$ of 30 K emerges at $x$=0.32, coinciding with a magnetic quantum critical point (QCP) which is evidenced by the disappearance of SDW order and the linear temperature-dependent resistivity in the normal state. The $T_c$ values were found to decrease with further P doping, and no superconductivity was observed down to 2 K for $x\geq$ 0.77. The appearance of superconductivity in the vicinity of QCP hints to the superconductivity mechanism in iron-based arsenides.

\section{Introduction}
The discovery of high-temperature superconductivity in iron
arsenides\cite{Kamihara08,chen-nature} was an exciting breakthrough
in the year 2008,\cite{science} which has led to a number of new
superconductors\cite{122,111,11} all containing antifluorite-like
FeAs or FeSe layers. In most cases, superconductivity was induced by
doping charge carriers into a stoichiometric parent compound which
shows a collinear antiferromagnetic spin-density-wave ground
state\cite{daiPC}. The chemical doping was initially carried out
only outside the FeAs layers,\cite{Kamihara08,122,Th,Sr} with the
electronic phase diagram\cite{CXH,Zhao,Gooch} resembling that of
cuprate superconductors. Later, however, it was found that the even
doping at the Fe site (with Co\cite{Co} or Ni\cite{Ni}) could induce
superconductivity. Besides, superconductivity was able to be
stabilized in the parent compounds simply by applying hydrostatic
pressures.\cite{P-Ca122} The latter two results contrast sharply
with those in cuprates.

Very recently, we observed superconductivity at 26 K in
EuFe$_2$As$_2$ by the partial substitution of arsenic with
phosphorus.\cite{ren} This isovalent substitution does not introduce
additional electrons or holes into the system, thus it can be
roughly understood in terms of chemical pressures. Owing to the
ferromagnetic ordering of Eu$^{2+}$ moments below 20 K coexisted
with superconductivity, the Meissner effect was not obvious.
Subsequently, we found that the P/As substitution induced
superconductivity in a prototype "1111" material LaFeAsO, but the
superconducting critical temperature $T_c$ is only $\sim$10
K,\cite{Wang2009} which is substantially lower than the $T_c$ value
under physical pressures\cite{Okada}. Here we employ the same doping
strategy on another prototype "122" parent compound
BaFe$_{2}$As$_{2}$\cite{Rotter}. Bulk superconductivity up to 30 K
with strong Meissner diamagnetism has been observed for $x$=0.32.
Our systematic measurements of resistivity and magnetization reveal
an electronic phase diagram characterized by a dome-like
superconducting region located near a magnetic quantum critical
point (QCP). The present work is consistent with a recent
theoretical prediction which shows a unique type of QCP resulting
from a competition between electronic localization and
itinerancy.\cite{DaiJH}

\section{Experimental}

Polycrystalline samples of BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$ (with
the nominal P content $x$=0, 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45,
0.5, 0.6, 0.7, 0.8 and 0.9 and 1) were synthesized by solid state
reaction with the high-purity elements Ba, Fe, As and P. Mixture of
the elements in the stoichiometric ratio was pressed into pellets in
an argon-filled glove-box. The pellets were sealed in evacuated
quartz tubes and annealed at 600 K for 12 h, then at 1173 K for 40
h. The resultant was reground, pelletized, and annealed in evacuated
quartz tubes at 1273 K for 24 h. Finally the samples were cooled
down to room temperature by shutting off the furnace.

Powder x-ray diffraction (XRD) was performed at room temperature
using a D/Max-rA diffractometer with Cu-K$_{\alpha}$ radiation and a
graphite monochromator. The structural refinement was performed
using the programme RIETAN 2000.\cite{Izumi} The electrical
resistivity was measured using a standard four-probe method. The
measurements of dc magnetization were performed down to 1.8 K on a
Quantum Design Magnetic Property Measurement System (MPMS-5). Both
the zero-field-cooling (ZFC) and field-cooling (FC) protocols were
employed under the magnetic field of 10 Oe.

\section{Results and discussion}

BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$ system crystallizes in
ThCr$_2$Si$_2$-type structure with the space group I4/\emph{mmm}.
Except the main phase, minor secondary phase was identified as
Fe$_{2}$P. Figure 1(a) shows the XRD pattern and its Rietveld
refinement profile of a typical BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$
sample with $x$=0.35 at room temperature. The two-phase Rietveld
refinement is successful, as indicated by the fairly good
reliability factors $R_{wp}$=8.5\% and $S$=1.8. It was found that
the fitted P content is a little lower than the nominal one,
consistent with the existence of impurity Fe$_{2}$P (about 10\% in
molar ratio). The refined lattice parameters plotted as functions of
$x$ are shown in Fig. 1(b). The cell volume decreases monotonically
with the P doping. We notice that, for the samples with no
significant Fe$_{2}$P impurity (\emph{e.g.}, $x$=0.45 and 0.9), the
data points of the cell volume fall exactly on the straight line, in
accordance with the Vegard's law. Therefore, the actual P content
can be estimated from the measured cell volume, which gives the
corrected P content $x'$=0.02, 0.13, 0.2, 0.24, 0.32, 0.36, 0.47,
0.56, 0.63 and 0.77 for the nominal $x$ values of 0.1, 0.2, 0.25,
0.3, 0.35, 0.4, 0.5, 0.6, 0.7 and 0.8, respectively. This correction
is basically (within the uncertainty of $\pm$0.02) consistent with
the Rietveld refinement results.

\begin{figure}
\center
\includegraphics[width=8cm]{fig1.eps}
\caption{(a) X-ray powder diffraction pattern and its Rietveld
refinement profile for the BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$
($x$=0.35) sample. The secondary minor phase is identified as
Fe$_{2}$P. (b) Relative lattice constants (compared with the undoped
compound) as functions of nominal P content $x$. The solid straight
line obeys the Vegard's law.}
\end{figure}

Table 1 compares the crystal structures between the undoped and
P-doped ($x'$=0.32) BaFe$_{2}$As$_{2}$. A remarkable difference is
seen from the position of As/P. The P doping makes the pnictogen
atoms move toward the iron sheets, which leads to a dramatic
decrease of FeAs-layer thickness. This explains why the decrease in
$c$-axis is more pronounced (compared with the $a$-axis), as shown
in Fig. 1(b). In fact, the spacing of FeAs layers decreases little
with P doping. This indicates that P doping impacts mainly on the
FeAs block layers. As a result, the bond angle of As-Fe-As along
[110] directions increases substantially.

\begin{table}
\caption{\label{tab:table1}Comparison of room-temperature crystal
structures between the undoped and P-doped ($x'$=0.32)
BaFe$_{2}$As$_{2}$ with the space group I$4/mmm$. The atomic
coordinates are as follows: Ba (0, 0, 0); Fe (0.5, 0, 0.25); As/P
(0, 0, $z$).}

\center
\begin{tabular}{lcr}
\hline
Compounds&BaFe$_{2}$As$_{2}$&BaFe$_{2}$(As$_{0.68}$P$_{0.32}$)$_{2}$\\
\hline
$a$ ({\AA}) & 3.9612(3) &3.9231(3)\\
$c$ ({\AA}) & 13.006(1) &12.805(1)\\
$V$ ({\AA}$^3$) & 204.08(2) & 197.08(2)\\
$z$ of As/P & 0.3546(2) & 0.3530(2)\\
FeAs-layer thickness ({\AA}) & 2.721(2) &2.638(2)\\
FeAs-layer spacing ({\AA}) & 3.782(2) &3.765(2)\\
As-Fe-As angle ($^{\circ}$) & 111.0(1) &112.2(1)\\
\hline
\end{tabular}

\end{table}

Figure 2 shows the temperature dependence of resistivity for
BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$. For the undoped
BaFe$_{2}$As$_{2}$, $\rho$ drops rapidly below 140 K, consistent
with previous report\cite{Rotter}. The resistivity anomaly was
revealed to be associated with a simultaneous structural and
spin-density-wave (SDW) transition.\cite{Bao} Upon doping with P,
the transition is suppressed, exhibiting the decrease of transition
temperature ($T_{SDW}$) and the weakening of the resistivity drop.
For example, $T_{SDW}$ is suppressed to $\sim$ 90 K for the sample
of $x'$=0.13, where only a resistivity kink appears. When the P
content increases to 20\%, a sudden decrease in resistivity is seen
below 10 K in addition to an observable kink at $\sim$ 70 K. The
resistivity drop at lower temperature is due to a superconducting
transition (confirmed by the magnetic measurements below). In the
intermediate doping range (0.24$\leq x\leq$0.63), superconductivity
with higher $T_c$s was observed. The maximum transition temperature
is 30 K for $x'$=0.32. The resistive transition tail is probably due
to the influence of grain boundaries and/or minor ferromagnetic
Fe$_2$P impurity, rather than sample inhomogeneity, since nearly
perfect diamagnetism was measured (see below). For the higher P
doping ($x\geq$0.77), no superconductivity was observed down to 3 K.
One notes that the residual resistance ratio (RRR), defined by the
ratio of the resistance at room temperature to the residual
resistance at very low temperature, is 20, 35 and 75 for $x'$=0.77,
0.9 and 1.0, respectively. The high RRR values suggest high quality
of the samples.

\begin{figure}
\center
\includegraphics[width=8cm]{fig2.eps}
\caption{Temperature dependence of resistivity for
BaFe$_{2}$(As$_{1-x'}$P$_{x'}$)$_{2}$ samples. The data are
normalized for the convenience of comparison. All the curves are
labeled with the corrected P content. The insets display expanded
plot for showing the fit of the normal-state resistivity with the
formula $\rho(T)=\rho_{0}+\alpha T^{n}$. }
\end{figure}

To confirm superconductivity, the temperature dependence of dc
magnetic susceptibility for BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$ was
measured, as shown in Fig. 3. The sample of $x'$=0.2 shows strong
diamagnetic transition below 7 K, consistent with the drop in
resistivity below 10 K shown in Fig. 2. In spite of the large
superconducting fraction, however, it is difficult to ascertain
whether SDW order coexists microscopically with superconductivity or
not. For $x$=0.32, a stronger diamagnetic signal is observed below
30 K. The magnetic shielding is even beyond the ideal value of 100\%
apparently, since the demagnetization effect is not included here.
The Meissner volume fraction achieves $\sim$ 30\%, demonstrating the
bulk nature of superconductivity. Other superconducting samples also
displays similar bulk property.

\begin{figure}
\center
\includegraphics[width=8cm]{fig3.eps}
\caption{Superconducting diamagnetic transitions in
BaFe$_{2}$(As$_{1-x'}$P$_{x'}$)$_{2}$. Note that the magnetic
background coming from the ferromagnetic impurity of Fe$_2$P has
been removed. No correction of demagnetization has been made.}
\end{figure}

As is known, the normal-state resistivity contains valuable
information on electronic scattering and even superconducting
mechanisms. Dense polycrystalline samples may show intrinsic
transport property because the grain boundaries scatter electrons
like impurities. Thus the following analysis is valid even though
our data were obtained from polycrystalline samples. Assuming that
the intrinsic scattering produces a power-law resistivity at low
temperatures, the total resistivity can be expressed as
$\rho(T)=\rho_{0}+\alpha T^{n}$, where $\rho_0$ represents extrinsic
resistivity coming from the impurity/defect scattering. By fitting
the $\rho(T)$ data in the range of $T_c< T\leq$50 K, the exponent
$n$ was obtained. The inset of Fig. 4 shows that most samples have
the $n$-value close to 2.0, suggesting Fermi-liquid state due to
electron-electron interactions. Interestingly, $n$ is close to unity
(i.e., linear $T$-dependent normal-state resistivity) for the
samples with $x'$=0.3$\sim$0.45. The linear low-temperature
resistivity belongs to a non-Fermi liquid behaviour, as usually
observed in the regime of quantum criticality. Remember that the SDW
order is suppressed with the P doping, we argue that there exists a
QCP at $x'\sim$1/3 in the BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$
system. This result is quite consistent with a recent theoretical
prediction which shows a unique type of QCP resulting from a
competition between electronic localization and
itinerancy\cite{DaiJH}. We aware that similar QCP was also proposed
in (K,Sr)Fe$_{2}$As$_{2}$\cite{Gooch} system.

The electronic phase diagram is concluded in Fig. 4. The undoped
BaFe$_2$As$_2$ has the ground state of SDW order. With the P doping,
$T_{SDW}$ decreases gradually, forming a phase region of SDW. The
suppression of SDW can be easily explained in terms of $J_{1}-J_{2}$
model\cite{Si,Yildirim}, because $J_{2}$ is expected to decrease
upon P doping. On the other hand, the scenario of Fermi surface
nesting\cite{Dong} is not quite straightforward to meet the
experimental observation. In the vicinity of above-mentioned QCP, an
asymmetric superconducting (SC) dome emerges. The optimal doping
level (with the maximum $T_c$) coincides with the QCP, where a
non-Fermi-liquid behviour shows up. This strongly suggests that spin
fluctuations play an important role for the superconductivity. In
the overlapped area of SDW and SC, it is not clear whether SDW
coexists microscopically with SC at this moment. The upper-right
metallic region displays Fermi-liquid-like state, except that
non-Fermi-liquid behaviour appears near the QCP. This crossover in
electronic excitation states resembles those of
(K,Sr)Fe$_{2}$As$_{2}$\cite{Gooch} and
Ba(Fe$_{1-x}$Co$_{x}$)$_{2}$As$_{2}$\cite{PD122Co} systems.

\begin{figure}
\center
\includegraphics[width=8cm]{fig4.eps}
\caption{Electronic phase diagram in
BaFe$_{2}$(As$_{1-x'}$P$_{x'}$)$_{2}$. SDW and SC denote
spin-density-wave and superconductivity, respectively. The inset
shows the variations of the exponent $n$ in the formula
$\rho(T)=\rho_{0}+\alpha T^{n}$, obtained by the fitting of
normal-state resistivity at low temperatures.}
\end{figure}

It was reported that BaFe$_2$As$_2$ became superconducting up to 29
K under the high pressure of $\sim$4 GPa.\cite{P-SrBa122} The $T_c$
value is very close to the present one in
BaFe$_{2}$(As$_{1-x}$P$_{x}$)$_{2}$. This implies that chemical
pressure generated by P doping acts in a similar way to the physical
pressure. Meanwhile, the disorder induced by the P doping hardly
influences superconductivity. The relatively low $T_c^{max}$ value
(compared with 38 K in Ba$_{0.6}$K$_{0.4}$Fe$_2$As$_2$\cite{122})
can be well understood by an empirical structural rule for $T_c$
variations.\cite{Zhao} A very recent work\cite{Kimber} also reveals
similarities between structural distortions under pressure and
chemical doping in BaFe$_2$As$_2$. Therefore, it is the structural
details of FeAs layers (corresponding to specific electronic
structures), rather than charge doping level, that controls the
ground states in iron arsenides.

To summarize, BaFe$_2$As$_2$ has been tuned superconducting through
an isovalent doping of phosphorus at the arsenic site, for the first
time. Unlike previous doping strategy with changing the
concentration of charge carrier, the present P doping does not alter
the number of valence electrons. Structural measurements indicates
that the chemical pressure generated by the P doping impacts mainly
on the FeAs layers. The P-for-As substitution suppresses the SDW
ordering, and favors superconductivity, similarly to phenomena of
the carrier doping\cite{122,CXH,Zhao}. A magnetic QCP is evidenced
simultaneously by the disappearance of SDW ordering and the
non-Fermi-liquid behaviour. The maximum $T_c$ of 30 K is observed
where the QCP is located. Theoretical studies\cite{Moriya} indicate
that superconductivity around magnetic QCP is associated with the
spin fluctuation mechanism. Thus the present result suggests that
spin fluctuations play crucial role for the superconductivity in
iron-based arsenides.

\section*{Acknowledgments}
This work is supported by the NSF of China (No. 10674119), National
Basic Research Program of China (No. 2007CB925001) and the PCSIRT of
the Ministry of Education of China (IRT0754).

\section*{References}

Superconductivity (cond-mat.supr-con); Strongly Correlated Electrons (cond-mat.str-el)