Young Dwarf Galaxies and Cosmology

Abstract: We develop here the theme that extremely metal-deficient Blue Compact Dwarf Galaxies (BCDs), those with Z < Zsun/20, are young galaxies which did not start to form stars until ~ 100 Myr ago. They can be thus considered as primeval galaxies in our local volume of the universe, and are excellent laboratories for studying physical processes occuring at the time of galaxy formation in very metal-deficient environments. We use BCDs to derive a new value of the primordial helium abundance Yp = 0.244+/-0.001, higher than previous determinations. This corresponds to a baryon mass fraction Omega_b*h^2_{50} = 0.058+/-0.007. We discuss the problem of dust in a very low-metallicity environment and show how a large fraction (~ 75%) of the young stars are not visible, and that any derivation of the cosmic star formation rate based only on optical/UV fluxes would be an underestimate. We also show how the loss of Ly-alpha photons from starburst regions puts strong constraints on Ly-alpha searches of high-redshift galaxies.

\section{Introduction}

 The formation of galaxies is one of the most fundamental
 problems in astrophysics, and much effort has gone into the
 search for primeval galaxies (PG). A possible definition of a primeval
  galaxy is a young system undergoing
its first major burst of star formation. It is now widely
believed that the vast majority of galaxies underwent such a phase at
redshifts $\sim$ 2 or greater. In most galaxy formation scenarios,
young galaxies are predicted to show strong Ly$\alpha$ emission,
associated with the cooling of the primordial gas and the subsequent
formation of a large number of massive  ionizing stars
(Partridge \& Peebles 1967; Charlot \& Fall 1993).
 Yet, despite intensive searches, the predicted
widespread population of Ly$\alpha$ primeval galaxies has remained
elusive (Pritchett 1994).

      Several objects have been put forward as possible PG candidates,
 ranging from high-redshift radio galaxies to Ly$\alpha$
emitters found around quasars and damped Ly$\alpha$ systems,
mainly on the basis of very high luminosity and
star formation activity. However, most of these
candidate PGs already contain a substantial amount of heavy
elements, as evidenced by
the presence of strong P Cygni profiles and interstellar absorption in their
spectra (Steidel et al. 1996; Yee et al. 1996).  These
spectra are very similar to those of nearby starburst galaxies known to contain
old stellar populations (Leitherer et al. 1996).
Thus high-redshift galaxies discovered thus far  are not truly primeval.
Moreover, even if true PGs are discovered at high-redshift,
 it is difficult to study them
in detail because of their extreme
faintness and very compact angular size. We propose here to take a different
approach to the PG problem. Instead of searching for very high-redshift
galaxies in the process of forming, we look for nearby galaxies undergoing
their first burst of star formation, and hence satisfying the above definition
of a PG. The best candidates for such a search are blue compact dwarf
galaxies (BCD).

      BCDs are  low-luminosity
 extragalactic objects with $M_B$ $\geq$ --18 where intense star
formation is presently occuring, as evidenced by their blue $UBV$ colors, and
their optical spectra which show strong narrow emission lines superposed on a
stellar continuum which is rising toward the blue, similar to spectra of HII
regions. Star formation in BCDs cannot be continuous but must proceed by
bursts because of several observational constraints: 1) Gas is transformed
into stars at the rate of approximately 1 $M_\odot$ yr$^{-1}$, so that the
current burst cannot last more than about 10$^8$ yr before depleting the
neutral gas supply of $\sim$ 10$^8$ $M_\odot$;
2) Optical-infrared colors of BCDs give burst ages of about 10$^7$ yr;
and 3) Population synthesis of UV spectra of BCDs give invariably jumps
in the stellar luminosity function, indicative of starbursts (see Thuan 1991
for a review).

      Ever since their discovery, the question has arisen whether BCDs are
truly young systems where star formation is occuring for the first time, or old
galaxies with an old underlying stellar population on which the current
starburst is superposed (Searle, Sargent \& Bagnuolo 1973). Thuan (1983)
carried out a near-infrared $JHK$ survey of BCDs and concluded that all the
objects in his sample possessed an old underlying stellar population of K and
M giants. That result was not unambiguous as the $JHK$ observations were
centered on the star-forming regions and the near-infrared emission could be
contaminated by light from young supergiant stars. The advent of CCD
detectors allowed to look for the low-surface-brightness underlying component
directly. Loose \& Thuan (1985) undertook a CCD imaging survey of a large BCD
sample and found that nearly all galaxies ($\geq$ 95\%) in their sample show an
underlying extended low-surface-brightness component, on
which are superposed the high-surface-brightness star-forming regions.
Subsequent CCD surveys of BCDs have confirmed this initial result
(Papaderos et al. 1996, Telles \& Terlevich 1997).
Thus, most BCDs are not necessarily young galaxies.
However, there was a hint that extremely metal-deficient BCDs
do not contain an old stellar population and can be primordial.
{\sl Hubble Space Telescope (HST)}
imaging of I Zw 18, the most metal-deficient BCD known
 ($Z_\odot$/50, Searle \& Sargent 1972), to $V$ $\sim$ 26 by
 Hunter \& Thronson (1995) suggests that the stellar population is
dominated by young stars and that the colors of the underlying diffuse
component are consistent with those from a sea of unresolved B or early
 A stars, with no evidence for stars older than $\sim$10$^7$yr.

      For more than 20 years, I Zw 18 stood in a class by itself. The BCD
metallicity distribution ranges from $\sim$ $Z_\odot$/3 to $\sim$ $Z_\odot$/50,
peaking at $\sim$ $Z_\odot$/10, and dropping off sharply for $Z$
$\leq$ $Z_\odot$/10.
Intensive searches have been carried out to look for
low-metallicity BCDs but they have met until recently with limited success.
 Several years ago,
a new BCD sample has been assembled by Izotov et al. (1993) from
objective prism survey plates obtained with the 1m Schmidt telescope
at the Byurakan Observatory of the Armenian Academy
of Sciences during the Second Byurakan Survey (SBS). The most
interesting feature of the SBS is its metallicity
distribution (Izotov et al. 1992, Thuan et al. 1994): it contains significantly
 more low-metallicity BCDs than previous surveys. It has uncovered
about a dozen BCDs with $Z$ $\leq$ $Z_\odot$/15, more than doubling the
number of such known low-metallicity BCDs and filling in the metallicity gap
between I Zw 18 and previously known BCDs.

In section 2, we use the low-metallicity SBS sample together
with other metal-deficient BCDs to study heavy element abundance ratios in
very low-metallicity environments
and to argue that all galaxies with a metallicity
less than about 1/20 of solar metallicity are young,
i.e. they did not start to form stars
until about 100 Myr ago. In that sense very metal-deficient local
BCDs are truly PGs,
and their study can shed light on galaxy formation at high redshift.
Because metal-deficient BCDs are young, they constitute excellent laboratories
for determining the primordial Helium abundance. We use a large sample of
low-metallicity BCDs in section 3 to derive a new value
of the primordial helium abundance which is appreciably higher than the
 previously
accepted value, and a baryonic mass density more in agreement with
other observations.
We discuss next in detail two BCDs with $Z$ $\leq$ $Z_\odot$/20,
SBS 0335--052 (section 4) and SBS 1415+437 (section 5). We show that
photometric and spectroscopic data on these two galaxies
also point to a young age for the two BCDs, of less than 100 Myr.
In the concluding sections we show how BCDs can shed light on
cosmological issues such as dust in a very low-metallicity environment
(section 6), and the escape of Ly$\alpha$ photons from starburst galaxies
(section 7).

\section{Galaxies with Z $\leq$ Z$_\odot$/20 are younger than 100 Myr}

\subsection{$\alpha$-elements}

The study of the variations of one chemical element relative
to another is crucial for our understanding of the chemical evolution
of galaxies and for constraining models of stellar nucleosynthesis and
the shape of the initial mass function. In the case of BCDs,
it is  particularly important for understanding their
evolutionary status, whether they are young or old.
Izotov \& Thuan (1999) have obtained very high-quality ground-based
spectroscopic
observations of 54 supergiant H II regions in 50 low-metallicity blue
compact galaxies with oxygen abundances 12 + log O/H between 7.1 and  8.3
($Z_\odot$/50 $\leq$ $Z$ $\leq$ $Z_\odot$/4).  They
use the data to determine abundances for the
elements N, O, Ne, S, Ar and Fe. They also analyze {\sl Hubble Space
Telescope} ({\sl HST}) Faint Object Spectrograph archival spectra of
10 supergiant H II regions to derive C and Si abundances in a subsample of 7
BCDs.
 The best studied and most easily observed element in BCDs is oxygen.
Nucleosynthesis theory predicts it to be produced only by high-mass
($M$ $>$ 9 $M_\odot$) stars. We
shall use it as the reference chemical element and consider the behavior of
heavy element abundance ratios as a function of oxygen abundance.
Figures 1 d,e,f,g show the dependence of the
abundance ratios Ne/O, Si/O, S/O and Ar/O on oxygen abundance.
The elements neon, silicon, sulfur and argon are all products of
$\alpha$-processes during both hydrostatic and explosive nucleosynthesis
in the same massive stars which make oxygen. Therefore, the Ne/O, Si/O,
S/O and Ar/O ratios should be constant and show no dependence on the oxygen
abundance. As predicted by stellar nucleosynthesis theory,
 {\it none} of the above heavy element-to-oxygen abundance ratios
  depend on oxygen
abundance. The mean values of these element abundance ratios are directly
 related to the stellar yields and thus provide strong constraints on the
theory of massive stellar nucleosynthesis
 (Thuan et al. 1995, Izotov \& Thuan 1999).

\def1.0{0.5}
\begin{figure}[t]
\centering
 \vbox{
	\special{psfile="thuan-fig1l.ps" hscale=37 vscale=37
	       hoffset=0 voffset=-295 angle=0}}
	\special{psfile="thuan-fig1r.ps" hscale=37 vscale=37
	       hoffset=230 voffset=-295 angle=0}
\vspace{10.cm}
\caption[]{\label{Fig1}{\bf Left panel:} C/O, N/O and Fe/O abundance
ratios vs oxygen abundance for the sample of low-metallicity BCDs from
Izotov \& Thuan (1999). Solid lines for C/O are theoretical predictions
from high-mass stars (HMS) evolution (Woosley \& Weaver 1995); the dashed
line shows the  ratio predicted from the evolution of HMS and intermediate-mass
stars (IMS) (Renzini \& Voli 1981). The solid line for N/O is the mean
observed ratio adopted as the primary nitrogen-to-oxygen abundance ratio
produced by massive stars.
{\bf Right panel:} Ne/O, Si/O, S/O and Ar/O abundance ratios vs oxygen
abundance for the same sample of BCDs. Solid lines are ratios predicted
from HMS evolution.}
\end{figure}
\def1.0{1.0}

\subsection{Iron, Carbon and Nitrogen}

  Because O, Ne, Si, S and Ar are made in the same high-mass stars, their
abundance ratios with respect to O are constant and not sensitive to
the age of the galaxy. By contrast, C, N and Fe can be produced by both high
and intermediate-mass (3 $M_\odot$ $\leq$ $M$ $\leq$ 9 $M_\odot$) stars,
and their abundance ratios with respect to O
give important information on the evolutionary status of BCDs.
The constancy of [O/Fe] for the BCDs  and its high
value compared to the Sun (Fig.1 c) suggests that all iron was produced by
massive stars, i.e. in SNe II only. Since the time delay between
iron production from SNe II and SNe Ia is about 1 -- 2 Gyr, it is likely
that BCDs with oxygen
abundance less than 12 + log O/H $\sim$ 8.2 ($Z_\odot$/5)
 are younger than 1 -- 2 Gyr.

 The behavior of the C/O and N/O ratios as a function of oxygen
abundance (Fig.1 a,b) puts more stringent constraints on the age of BCDs
(Thuan et al. 1995, Izotov \& Thuan 1999).
The behavior of the C/O and N/O ratios
is very different whether the BCG has 12 + log O/H smaller or greater than 7.6
($Z_\odot$/20).
The remarkably small scatter of the C/O and N/O abundance ratios in BCDs
with $Z$ $\leq$ $Z_\odot$/20 rules out any time-delay model
in which O is produced first by massive stars and C and N are produced later by
intermediate-mass stars, and supports a common origin of C,
N and O in the same first-generation massive stars. Thus, it is very
likely that the presently observed episode of star formation in BCDs with $Z$
 $\leq$ $Z_\odot$/20 is the first one in the history of the galaxy and the age of
the oldest stars in it do not exceed $\sim$ 40 Myr, the lifetime of a
9 $M_\odot$ star.
    The conclusion that BCGs with $Z$ $\leq$ $Z_\odot$/20 are young is
supported by the analysis of {\sl HST} WFPC2 images of some of these galaxies:
I Zw 18 ($Z_\odot$/50, Hunter \& Thronson 1995), SBS 0335--052 ($Z_\odot$/41,
 Thuan et al. 1997, see section 4),
SBS 1415+437 ($Z_\odot$/21, Thuan et al. 1999, see section 5),
T1214--277 and Tol 65 (respectively $Z_\odot$/21 and $Z_\odot$/22,
Izotov et al. 1999).

   The situation changes for BCDs with $Z$ $>$ $Z_\odot$/20. The scatter of the
C/O and N/O ratios increases significantly at a given O abundance, which was
interpreted by Thuan et al. (1995) and Izotov \& Thuan (1999) as due to the
additional production of primary N by
intermediate-mass stars, on top of the primary N production by high-mass stars.
Thus, since it takes at least 500 Myr (the lifetime of a 2 -- 3 $M_\odot$ star)
for C and N to be produced by
intermediate-mass stars, BCDs with $Z$ $>$ $Z_\odot$/20 must have had several
episodes of star formation before the present one and they must be at least
older than $\sim$ 100 Myr. This conclusion is in agreement with photometric
studies of these higher metallicity BCDs which, unlike their very
low-metallicity counterparts, have a red old instead of a blue young underlying
stellar component (Loose \& Thuan 1985, Papaderos et al. 1996).

\section{The primordial helium abundance}

We use the data of
Izotov et al. (1994, 1997a) and Izotov \& Thuan (1998b) to construct a
sample of 45 HII regions appropriate for the determination of  $Y_p$.
Our sample constitutes one of the largest and most homogeneous
(obtained, reduced and analyzed in the same way) data set now available for the
determination of  $Y_p$. This is the same data as used to analyze
heavy element abundances in section 2.
Linear regressions of the $Y$ -- O/H and $Y$ -- N/H relations, with
$Y$s determined from a self-consistent treatment of the five brightest optical
He I emission lines, gives $Y_p$ = 0.244$\pm$0.001
(Figures 2 a,b, Izotov \& Thuan 1998b).
This value agrees very well with
the mean $Y$ of the two most metal-deficient BCDs
known [I Zw 18 ($Y_p$ = 0.242$\pm$0.009, Izotov \& Thuan 1998a)
and SBS 0335--052 ($Y_p$ = 0.249$\pm$0.004, Izotov \& Thuan 1998b)], which
is $\bar{Y}$ = 0.245$\pm$0.006. Values as low as
$Y_p$ = 0.234 or $Y_p$ = 0.230, as those obtained by Pagel et al. (1992) and
Olive et al. (1997) are excluded. Part of the difference comes from the
fact that previous workers have neglected underlying He I stellar
absorption which would artificially lower $Y_p$ (Izotov \& Thuan 1988a).

\def1.0{0.5}
\begin{figure*}[t]
\centering
 \vbox{
	\special{psfile="thuan-fig2l.ps" hscale=37 vscale=37
	       hoffset=0 voffset=-295 angle=0}}
	\special{psfile="thuan-fig2r.ps" hscale=37 vscale=37
	       hoffset=230 voffset=-295 angle=0}
\vspace{10.cm}
\caption[]{\label{Fig2}{\bf Left panel:}
Linear regressions of (a) the helium mass fraction
$Y$ vs. oxygen abundance O/H and (b) the helium mass fraction $Y$ vs.
nitrogen abundance for our sample of
45 H II regions (Izotov \& Thuan 1998b). The $Y$s are derived self-consistently
by using the 5 brightest He I emission lines in the optical range.
Collisional and fluorescent enhancements, underlying He I stellar absorption and
Galactic Na I interstellar absorption are taken into account.
1$\sigma$ alternatives are shown by dashed lines.
{\bf Right panel:}
The abundances of a) $^4$He; b) D and c) $^7$Li as a
function of $\eta_{10}$$\equiv$10$^{10}$$\eta$, where $\eta$ is the
baryon-to-photon number ratio, as given by the standard
hot big bang nucleosynthesis model.}
\end{figure*}
\def1.0{1.0}

Our higher primordial helium mass fraction is consistent
with deuterium abundance measurements in high-redshift Ly$\alpha$ absorbing
systems by Tytler et al. (1996) as corrected by turbulence effects in the
clouds by Levshakov et al. 1998), and lithium abundance measurements in
low-metallicity halo stars (Bonifacio \& Molaro 1997). Figure 2c gives
a baryon-to-photon ratio $\eta$= 4.0$\pm$0.5, which corresponds to
a baryon mass fraction  $\Omega_b$$h^2_{50}$ = 0.058$\pm$0.007, in
good agreement with the value derived from X-ray observations of clusters of
galaxies (White \& Fabian 1995).
We derive a slope $dY/dZ$ = 2.3$\pm$1.0, considerably smaller than those derived
before.
 With this smaller slope and taking into account the errors,
 chemical evolution models with an outflow of well-mixed material can be
  built for star-forming dwarf
galaxies which satisfy all the observational constraints.

\section{The young dwarf galaxy SBS 0335--052}

\subsection{Age of the stellar component}

This BCD with $Z$ = $Z_\odot$/41 is the second most metal-deficient galaxy known
(after I Zw 18).
HST WFPC2 images of the BCD show that most of the star formation in SBS
0335--052 ($M_B$ = -- 16.7, $v$ = 4076 km s$^{-1}$)
occurs in 6 super-star clusters (SSCs) with --14.1 $\leq$ $M_V$
$\leq$ --11.9, within a region $\sim$ 520 pc in size (Figure 3,
Thuan et al. 1997).
Later processing by Papaderos et al. (1998) of the same HST images reveals
several fainter clusters (not SSCs). The SSCs are roughly aligned in the SE-NW
direction, and there is a systematic reddening of the $V-I$ color  of the SSCs
away from the brightest one, with a flattening of the color  of the clusters
beyond 520 pc (Figure 4).

\def1.0{0.5}
\begin{figure*}[t]
 \vbox{
	\special{psfile="thuan-fig3.ps" hscale=40 vscale=40
	       hoffset=0 voffset=-285 angle=0}}
	\special{psfile="thuan-fig4.ps" hscale=40 vscale=40
	       hoffset=250 voffset=0 angle=-90}
\vspace{8.5cm}
{\label{Fig3_4}{\bf Fig. 3 (left):}
{\sl HST} WFPC2 $V$ image  of SBS 0335--052 showing the high surface brightness
super-star clusters. At a distance of 54.3 Mpc, 1 arcsec corresponds to a linear size of
263 pc.
{\bf Fig. 4 (right):}
$(V-I)$ color vs. distance from the brightest super-star cluster at the SE tip
of SBS 0335--052.
The color gets redder with increasing distance.}
\end{figure*}
\def1.0{1.0}

\def1.0{0.5}
\begin{figure*}[t]
 \vbox{
	\special{psfile="thuan-fig5.ps" hscale=40 vscale=40
	       hoffset=0 voffset=-285 angle=0}}
	\special{psfile="thuan-fig6.ps" hscale=40 vscale=40
	       hoffset=250 voffset=0 angle=270}
\vspace{8.5cm}
{\label{Fig5_6}{\bf Fig.5 (left):}
Same $V$ image as in Fig.3, with the contrast adjusted to show the low surface
brightness underlying component of SBS 0335--052.  The supershell delineating the
large supernova cavity is clearly seen.

{\bf Fig.6 (right):} $(U-B)$ color profile (labeled E) of SBS 0335--052. The color is very blue and
the profile flat.
 }
\end{figure*}
\def1.0{1.0}

Thuan et al.
(1997) and Papaderos et al. (1998) attribute
most of the color variation of the SSCs
 to an age variation resulting from sequential propagating  star
formation. The $V-I$ colors are consistent with the picture that star formation
started at the location of the most distant cluster,
some 1.8 kpc away from the
location of the brightest and bluest cluster  at the South East
end of the galaxy, at about 100 Myr ago and  propagated
through the ISM to the latter, whose age is only $\sim$ 4 Myr,
with an average speed of
$\sim$ 18 km s$^{-1}$ (Thuan et al. 1997, Papaderos et al. 1998). Thus the
star-forming clusters have ages between 4 and 100 Myr. Does SBS 0335--052
possess an underlying older stellar population?
The extended underlying component
is shown in figure 5 where the contrast has been adjusted to display very
low surface-brightness features. The unusually blue colors of this underlying
 component
(see the $U-B$ color profile labeled E in Figure 6) and its irregular, blotchy
and
filamentary structure suggest that a significant fraction of the light is of
 gaseous rather
than stellar origin. A supershell of 380 pc radius can be seen delineating
 a large
supernova cavity. However Izotov et al. (1997) found that the H$\beta$
equivalent
width in the underlying component is $\sim$ 3 times lower than the value
expected for pure gaseous emission, implying that two-thirds of the light
comes from an underlying stellar population. Papaderos et al. (1998)
have modeled the $UBVRI$ colors of this stellar component, after removal of
the ionized gas contamination. They found that the colors are consistent
with an underlying stellar population not older than 100 Myr.

\subsection {Age of the neutral component}

    Thus the stellar component in SBS 0355--052 is extremely young and
the BCD is likely undergoing star formation for the first time. If this is the case,
the neutral gas envelope surrounding the BCD must also be very metal-deficient.
A 21 cm VLA map of the BCD (Pustilnik et al. 1999)  has shown it to be embedded
in an extraordinarily large HI cloud seen nearly edge-on,
with dimensions some 64 by 24 kpc. This is
to be compared with the typical size of HI envelopes around BCDs which is more
like a few kiloparsecs in each dimension.  We can use the BCD as a background
light source shining through the HI envelope to probe the physical conditions of the
neutral gas. The Ly$\alpha$ line seen in absorption would give the column density
of atomic hydrogen, while the OI $\lambda$1302 line would give the column density
of the most abundant heavy element that remains neutral in the HI cloud. This
would allow us to set limits on the O/H abundance ratio in the neutral gas.
 Figure 7
shows the ultraviolet spectrum of SBS 0335--052 around the Ly$\alpha$ line obtained
by Thuan \& Izotov (1997)
with the Goddard High Resolution Spectrograph aboard the {\sl Hubble Space Telescope (HST)}.
A strong damped Ly$\alpha$ absorption line is seen along with several heavy element
interstellar absoption lines such as OI $\lambda$1302, SiII $\lambda$1304 and
SII $\lambda$1251, $\lambda$1254, and $\lambda$1259. The HI column
derived by fitting the Ly$\alpha$ absorption profile is $N$(H I) =
(7.0$\pm$0.5)$\times$10$^{21}$ cm$^{-2}$, the highest
derived thus far for a BCD, and $\sim$2 times larger than in
I Zw 18 (Kunth et al. 1994).  Comparison
with high-resolution quasar spectra implies that the O I
$\lambda$1302 line along with other heavy element interstellar absorption lines
such as Si II $\lambda$1304 and S II $\lambda$1251, $\lambda$1254,
$\lambda$1259 are not saturated, which allow us to derive abundances.
Assuming that these lines originate in the H I gas, we derive
extremely low abundances of
oxygen, silicon and sulfur, respectively 37000, 4000 and 116 times lower than
the solar values. The oxygen abundance is a whole 37 times lower than in the
neutral gas of I Zw 18. However, these highly discrepant deficiency factors
between different elements suggest that the absorption lines are
produced, \underline{not in the H I, but in the H II gas}.
Adopting that hypothesis, the derived abundance from the $UV$ absorption lines
are then consistent with that derived from the optical emission lines ($Z$
$\sim$ $Z_\odot$/41). The conclusion that the heavy element absorption lines
originate in the H II region is supported by the detection of several systems of
blueshifted S II $\lambda$1259, Si II $\lambda$1260, O I $\lambda$1302,
Si II $\lambda$1304, C II $\lambda$1335 absorption lines originating in
fast-moving clouds with velocities up to $\sim$ 1500 km s$^{-1}$, and also by
the presence of heavy element
absorption lines with excited lower levels. If this
conclusion holds, then the H I cloud in SBS 0335--052 is truly
\underline{primordial}, unpolluted by heavy elements (Thuan \& Izotov 1997).
In summary, all the known observational evidence suggests that SBS 0335--052
is truly a young galaxy.

\section{Another young dwarf galaxy: SBS 1415+437}

\subsection {Age from color-magnitude diagrams}

    Figure 8 shows  the {\sl HST} WFPC2 $I$ image of SBS 1415+437
   ($M_B$ = -- 14.0, $v$ = 607 km s$^{-1}$)  with the
contrast level adjusted so as to show the
low surface-brightness underlying extended component
(Thuan et al. 1999). The galaxy has an
elongated, comet-like shape with a bright H II region on its SW tip. Many
point sources  identified as luminous stars can be seen. To the SW
of the brightest H II region, two stellar clusters with resolved stars are
present. The luminous stars and the HII regions are aligned suggesting,
just as in SBS 0335--052, propagating star formation (from the NE to the SW).
The mode of star formation in the two BCDs is different however. While
SBS 0335--052 makes stars in luminous super-star clusters, star-formation in
SBS 1415+437  appears to be less extreme and is more similar to that
in I Zw 18 (Hunter \& Thronson 1995). The superior spatial resolution of
the {\sl HST}
allows to resolve individual stars and construct color-magnitude diagrams to
study the stellar populations in the BCD. To check the hypothesis of propagating
star formation  , we have derived stellar ages for 6 separate regions
in the BCD, labeled from I to VI as shown in Figure 8.

  The $(V-I)$ vs. $I$ color-magnitude diagram (CMD) for each region
  are shown in Figure 9 together
with stellar isochrones by Bertelli et al. (1994)
for a heavy element abundance equal to 1/20 the solar value.
 Each isochrone is marked
by the logarithm of the age in years. The region of
the  asymptotic giant branch (AGB) stars is
shown by a dashed line while the observational limits
are shown by dotted lines.
We adopt  $V$ = 27.5 mag and  $I$ = 27 mag as completeness
limits. Although there is evidence for spatial variations of the
internal dust extinction, as a first approximation, we correct the CMDs for
regions III -- VI by a constant extinction as derived by the Balmer decrement.
The CMDs for regions I and II have been not corrected for extinction for
lack of information.

Inspection of Figure 9 shows that there is a clear age gradient from region I to
region VI. Region I is the youngest (about 5 Myr), containing only main-sequence
stars. Region II is more evolved. It contains red
supergiants and has an age of $\geq$ 10 Myr .
Region III includes the brightest H II region
and shows a mixture of stellar populations.
The bulk of the stars in region III have ages ranging between $\leq$ 10 Myr
and 100 Myr.  The red stars with $(V-I)$  between 1.0 and 1.8
are likely to be red supergiants rather than AGB stars.
The stellar populations in regions IV - VI are
similar to those in region III except for the fact that very young
populations with age less 10 Myr are no more present.
We emphasize that the properties of the stellar populations
 in SBS 1415+437 are quite different from
those in other nearby low-metallicity dwarf galaxies with
HST color-magnitude diagrams. In those dwarfs, very red AGB
stars are present  indicating a larger age
(e.g. Dohm-Palmer et al. 1997; Schulte-Ladbeck et al. 1998).
Two general conclusions can be obtained from the above color-magnitude
analysis:
1) star formation in SBS 1415+437 is propagating from the NE to the SW;
2) there is no evidence for stars older than $\sim$ 100 Myr in the BCD.

\def1.0{0.5}
\begin{figure*}[t]
 \vbox{
	\special{psfile="thuan-fig7.ps" hscale=40 vscale=40
	       hoffset=-20 voffset=0 angle=-90}}
	\special{psfile="thuan-fig8.ps" hscale=40 vscale=40
	       hoffset=250 voffset=-270 angle=0}
\vspace{8.5cm}
{\label{Fig7_8}{\bf Fig.7 (left):}
HST GHRS spectrum of SBS 0335--052 around the Ly$\alpha$ line. A strong damped
Ly$\alpha$ absorption is seen along with several heavy element interstellar absorption lines.

{\bf Fig.8 (right):} {\sl HST} WFPC2 $I$ image of SBS 1415+437. The shape is cometary-like.
A bright supergiant HII region at the SW end, resolved stars in the body of the galaxy
and unresolved diffuse stellar emission can be seen. The different regions used for
stellar population color-magnitude diagram analysis are marked.

 }
\end{figure*}
\def1.0{1.0}

\def1.0{0.5}
\begin{figure*}[t]
 \vbox{
	\special{psfile="thuan-fig9.ps" hscale=37 vscale=37
	       hoffset=-15 voffset=-20 angle=-90}}
	\special{psfile="thuan-fig10.ps" hscale=40 vscale=40
	       hoffset=265 voffset=-290 angle=0}
\vspace{10.cm}
{\label{Fig9_10}{\bf Fig.9 (left):} Color-magnitude diagrams of stellar populations in the different
regions of
SBS 1415+437 as defined in Fig.8. There is a systematic age increase from region I
to region VI, implying propagating star formation from the NE to the SW.
Typical photometric errors are shown in Fig. 9a.

{\bf Fig.10 (right):} MMT spectrum of region V (see Fig.7) of SBS 1415+437. The upper panel shows the
spectrum uncorrected for extinction, while the lower panel shows the same spectrum
corrected for extinction. The best-fit model gives an age between 30 and 100 Myr.
 }
\end{figure*}
\def1.0{1.0}

\subsection{Age from spectral evolutionary synthesis models}

The second conclusion is supported by evolutionary synthesis modeling of the
spectrophotometric data of regions III to VI.
To estimate quantitalively the age of each region, we calculate a grid
of  spectral energy distributions (SED)
for stellar populations with ages varying between 10 Myr and 20 Gyr and heavy
element abundance $Z_\odot$/20,  using isochrones from
Bertelli et al. (1994) and the compilation of stellar atmosphere models from
Lejeune et al. (1998). A Salpeter IMF with slope --2.35, an upper
mass limit of 120 $M_\odot$ and a lower
mass limit of 0.6 $M_\odot$ were adopted. For stellar populations with
 age less than 10 Myr we use theoretical spectral energy distributions by
Schaerer \& Vacca (1998) for a heavy element abundance $Z_\odot$/20
and a Salpeter IMF. The stellar emission in SBS 1415+437
is contaminated by emission of ionized gas from supergiant H II regions.
Therefore, to study the stellar composition in the BCD,  it is necessary
to produce a synthetic SED which includes both stellar and ionized gaseous
emission. It is clear from Figure 10, the best fit to the
spectrum corrected for extinction (lower panel)
 gives an age between 30 and 100 Myr for region V (the
models are labelled by the logarithm of the age in years). A similar analysis for
the other regions give the same answer: none contains stellar populations older than
100 Myr. Thus SBS 1415+437, just like SBS 0335--052,  is also a young galaxy.

\section{Dust in an extremely metal-deficient environment}

\def1.0{0.5}
\begin{figure*}[t]
 \vbox{
	\special{psfile="thuan-fig11.ps" hscale=40 vscale=40
	       hoffset=0 voffset=-285 angle=0}}
	\special{psfile="thuan-fig12.ps" hscale=50 vscale=50
	       hoffset=200 voffset=-400 angle=0}
\vspace{8.5cm}
{\label{}{\bf Fig.11 (left):}
$(V-I)$ color map of SBS 0335--052. Blue is dark, red is light. The dust patches on
top of the super-star clusters are clearly seen.

{\bf Fig.12 (right):} The best-fit model (thick continuous line) for the mid-infrared spectral
energy distribution of SBS 0335--052
obtained by ISO (Thuan et al. 1999).  It is composed of a 260\,K blackbody
spectrum modified by an emissivity law proportional to $\nu^{1.5}$ and
extinguished by a screen of dust of optical thickness $A_{V}$ = 21 mag,
with an extinction curve similar to that observed
toward the Galactic center  (Lutz et al. 1996).}
\end{figure*}
\def1.0{1.0}

SBS 0335--052, although it is extremely metal-deficient, clearly contains dust.
It is seen as
white patches in the $V-I$ color map of the BCD (Figure 11), and is spatially
mixed with the SSCs.
Thuan et al. (1999) have obtained ISO mid-infrared observations of the BCD.
  With
$L_{12{\mu}m}/L_{B}$ of 2.15, the galaxy is unexpectedly bright in the
mid-infrared for such a low-metallicity object. The mid-infrared
spectrum (Figure 12) is very unusual when compared to spectra of other
 star-forming
galaxies: 1) there is no emission from the so-called Unidentified
Infrared Bands, which we interpret
as an effect of destruction of carbon-based dust by the very high UV
energy density in SBS 0335--052 ; and 2) there are
no evident fine structure ionic lines, even though neon and sulfur lines
are usually quite bright in starburst galaxies. The absence of these lines can
be explained by a very strong continuum which decreases the equivalent
widths of the lines, making them difficult to detect.
  The spectral energy distribution from 5 to
17 microns can be fitted with a grey-body spectrum modified by extinction.
Silicate grains can account for the unusual shape of the MIR spectrum. The
required extinction law is similar to that observed toward the
Galactic Center and the optical depth is $A_{V}\,\sim$\,19-21 mag. Such a
large optical depth implies that a large fraction (as much as 75\%) of the
current star-formation activity in SBS 0335--052 is
hidden by silicate dust whose mass is in the range 3$\times$10$^3$ --
5$\times$10$^5$ $M_\odot$.

Thus nearly primordial environments can contain a significant amount of dust.
If the hidden star formation in SBS 0335--052 is typical of young galaxies at
high redshifts, then the cosmic star formation rate derived from UV/optical
fluxes would be underestimated.
This is in fact the result found by Flores et al. (1999) who carried out
deep 15 micron ISO surveys of distant galaxies with a
median redshift of $\sim$ 0.76. Those authors found that the
cosmic SFR derived from FIR luminosities, from 43 to 123 microns,
assuming that the FIR luminosity comes mostly from dust heating by young stars,
 is 2-3 times higher  than the SFR
estimated previously from optical/UV fluxes (Madau et al. 1996).

\section{The escape of Ly$\alpha$ photons}

   The study of nearby young galaxies can also shed light on a long standing
problem concerning the Lyman-alpha emission of primeval galaxies. In any
galaxy formation scenario, young galaxies are predicted to show strong
Ly$\alpha$ emission (rest frame equivalent width of $\sim$ 100\AA) associated
with the formation of a large number of massive ionizing stars (e.g. Charlot \&
Fall 1993). Yet, despite intensive searches, the predicted population of
Ly$\alpha$ primeval galaxies remained elusive (Pritchet 1994). The absence of
Ly$\alpha$ emission in high-redshift galaxies is reminiscent of the behavior
of Ly$\alpha$ emission in nearby
starburst and BCD galaxies, which is either absent
or greatly diminished.
The Ly$\alpha$/H$\beta$ line intensity ratio in those
galaxies with detected Ly$\alpha$ emission does not exceed 10, significantly
lower than the theoretical recombination ratio of 33. Some galaxies show strong
Ly$\alpha$ absorption rather than emission.
The favored explanation for such a reduction in Ly$\alpha$
emission from recombination values is redistribution of Ly$\alpha$
photons by multiple scattering in the H I envelope or
absorption of these
Ly$\alpha$ photons by dust in the star-forming region.
The latter mechanism would imply increasing Ly$\alpha$/H$\beta$ line intensity
ratios with decreasing metallicities, since presumably low-metallicity
objects contain less dust, and hence suffer less destruction of Ly$\alpha$
photons (Terlevich et al. 1993).
{ \sl HST} observations of the two most metal-deficient
BCDs known, I Zw 18 (Kunth et al. 1994) and SBS 0335--052 (Thuan et al.
1997, section 4) , show Ly$\alpha$ not in emission but absorption
(Figure 7), which goes against a
Ly$\alpha$ strength-metallicity anticorrelation.
Lequeux et al. (1995), in their study of the BCD Haro 2 ($Z_\odot$/3) which
shows Ly$\alpha$ in emission, have
argued that Ly$\alpha$ photons can escape when the neutral material where
the absorption occurs is outflowing with a velocity of $\leq$
200 km s$^{-1}$ with respect to the star-forming region.
The Ly$\alpha$ emission is redshifted with
respect to both the H II region and the expanding absorbing shell, the motion
of which is probably powered by stellar winds and supernovae.
This explanation does not apply to the BCD
T1214--277 ($Z_\odot$/23) which is the
lowest metallicity galaxy known with detected Ly$\alpha$ emission (Thuan
\& Izotov 1997).
Contrary to Haro 2,
 Ly$\alpha$ emission is not redshifted with respect to the H II gas velocity.
Thus the escape of Ly$\alpha$ photons in T1214--277 is not a consequence of the
motions of the neutral H I envelope.

As for  SBS 0335--052, dust extinction may play some role as it is
directly seen (Figure 11). While there is
evidence for fast gas motions in SBS 0335--052 with velocities up to
$\sim$ 1500 km s$^{-1}$, the H I gaseous envelope appears to be static with
respect to the H II region, as the 21 cm and emission-line velocities are in
good agreement. Thus, with its extremely large H I column density, the
redistribution of Ly$\alpha$ photons in SBS 0335--052 by multiple scattering
over the large volume of the H I cloud probably plays also an important role
in diminishing the intensity of the Ly$\alpha$ line. The orientation of the
HI cloud may also play a role. In the case of SBS 0335--052 , the HI envelope
 is reasonably flattened (Pustilnik et al. 1999) suggesting it is seen nearly
edge-on. In that case, Ly$\alpha$
photons escape more easily along directions perpendicular to the line of sight
than along it. Moreover, as discussed by Giavalisco et al. (1996), the escape of
 Ly$\alpha$
photons may be controlled not only by the geometry but also by the
porosity of the neutral gas.

In summary, there is no unique mechanism which controls the appearance
of Ly$\alpha$ emission in nearby young dwarf galaxies. Dust extinction may
play a role, but the velocity structure of the H I gas, the orientation of
the H I cloud and its porosity may also be important factors.
The fact that some  BCDs do not show Ly$\alpha$ in emission implies that
Ly$\alpha$ searches for high-redshift galaxies will always be incomplete.

\acknowledgements{We thank the partial financial support of NSF grant
AST-9616863 and NASA grant JPL961535, and Polichronis Papaderos, Simon
Pustilnik and Marc Sauvage for useful discussions.}

\begin{bloisbib}
\bibitem{} Bertelli, G., Bressan, A., Chiosi, C., Fagotto, F., \&
Nasi, E. 1994, A\&AS, 106, 275
\bibitem{} Bonifacio P., Molaro P., 1997, MNRAS, 285, 847
\bibitem{} Charlot, S. \& Fall, S.M. 1993, ApJ, 415, 580
\bibitem{} Dohm-Palmer, R. C., Skillman, E. D., Saha, A., Tolstoy, E.,
Mateo, M., Gallagher, J., Hoessel, J., Chiosi, C., \& Dufour, R. J. 1997,
AJ, 114, 2527
\bibitem{} Flores, H., Hammer, F., Thuan, T.X., Cesarsky, C., Desert, F.-X.,
Omont, A., Lilly, S.J., Eales, S., Crampton, D. \& Le Fevre, O. 1999, ApJ, in press
\bibitem{} Giavalisco, M., Koratkar, A., \& Calzetti, D. 1996, ApJ,
461, 831
\bibitem{} Hunter, D. A., \& Thronson, H. A. 1995, ApJ, 452, 238
\bibitem{} Izotov, Y.I. \& Thuan, T.X. 1998a, ApJ, 497, 227
\bibitem{} Izotov, Y.I. \& Thuan, T.X. 1998b, ApJ, 500, 188
\bibitem{} Izotov, Y.I. \& Thuan, T.X. 1999, ApJ, February 1
\bibitem{} Izotov, Y.I., Thuan, T.X., \& Lipovetsky, V.A. 1994, ApJ, 435, 647
\bibitem{} Izotov, Y.I., Thuan, T.X., \& Lipovetsky, V.A. 1997a, ApJS, 108, 1
\bibitem{} Izotov, Y.I., Thuan, T.X., \& Papaderos, P. 1999, in preparation
\bibitem{} Izotov,Y.I., Lipovetsky, V.A., Guseva, N.G., Stepanian, J.A., Erastova, L.K.,
\& Kniazev, A.Y. 1992, in The Feedback of Chemical Evolution on the Stellar Content of
Galaxies, ed. D. Alloin \& G. Stasinska, (Paris: Observatoire de Paris), 127
\bibitem{} Izotov, Y.I., Guseva, N.G., Lipovetsky, V.A., Kniazev, A.Y., Neizvestny, S.I.,
\& Stepanian, J.A. 1993, Astron. Astrophys. Trans., 3, 197
\bibitem{} Izotov, Y. I., Lipovetsky, V. A., Chaffee, F. H., Foltz, C. B.,
Guseva, N. G., \& Kniazev, A. Y. 1997b, ApJ, 476, 698
\bibitem{} Kunth, D., Lequeux, J., Sargent, W. L. W., Viallefond, F.
1994, A\&A, 282, 709
\bibitem{} Leitherer, C., Vacca, W.D., Conti, P.S., Filippenko, A.V.,
 Robert, C. \& Sargent, W.L.W. 1996, ApJ, 465, 717
\bibitem{} Lejeune, T., Cuisinier, F., \& Buser, R. 1998, A\&AS, 130, 65
\bibitem{} Lequeux, J., Kunth, D., Mas-Hesse, J. M., \& Sargent, W. L. W.
1995, A\&A, 301, 18
\bibitem{} Levshakov, S.A., Kegel, W.H., \& Takahara,F. 1998, ApJ, 499, L1
\bibitem{} Loose, H.-H., \& Thuan, T.X. 1985, in Star-forming Dwarf Galaxies and
Related Objects, ed. D. Kunth, T.X. Thuan, \& J.T.T. Van (Gif-sur-Yvette: Editions
Frontieres), 73
\bibitem{} Lutz, D., et al. 1996, A\&A, 315, L269
\bibitem{} Madau, P., Ferguson, H., Dickinson, M., Giavalisco, M., Steidel, C.,
\& Fruchter, A. 1996, MNRAS, 283, 1388
\bibitem{} Olive, K.A., Skillman, E.D., \& Steigman, G. 1997, ApJ, 483, 788
\bibitem{} Pagel, B.E.J., Simonson, E.A., Terlevich, R.J. \& Edmunds, M.G. 1992,
MNRAS, 255, 325
\bibitem{} Papaderos, P., Loose, H.-H., Thuan, T. X., \& Fricke, K. J.
1996, A\&AS, 120, 207
\bibitem{} Papaderos, P., Izotov, Y.I., Fricke, K.J., Thuan, T.X. \& Guseva, N.G.
1998, A\&A, 338, 43
\bibitem{} Partridge, R.B., \& Peebles, P.J.E. 1967, ApJ, 147, 868
\bibitem{} Pritchet, C.J. 1994, PASP, 106, 1052
\bibitem{} Pustilnik, S. A., Thuan, T. X., Brinks, E., Lipovetsky, V. A., \&
Izotov, Y. I. 1999, in preparation
\bibitem{} Renzini A., \& Voli M., 1981, A\&A, 94, 175
\bibitem{} Schaerer, D., \& Vacca, W. D. 1998, ApJ, 497, 618
\bibitem{} Searle, L., \& Sargent, W. L. W. 1972, ApJ, 173, 25
\bibitem{} Searle, L., Sargent, W.L.W. \& Bagnuolo, W.G. 1973, ApJ, 179, 427
\bibitem{} Schulte-Ladbeck, R. E., Crone, M. M., \& Hopp, U. 1998,
ApJ, 493, L23
\bibitem{} Steidel, C.C., Giavalisco, M., Pettini, M., Dickinson, M.
\& Adelberger, K.L. 1996, ApJ, 462, L17
\bibitem{} Telles, E., \& Terlevich, R. 1997, MNRAS, 286, 183
\bibitem{} Terlevich, E., D\'{i}az, A. I., Terlevich, R., \&
Garc\'{i}a Vargas, M. L. 1993, MNRAS, 260, 3
\bibitem{} Thuan, T.X. 1983, ApJ, 268, 667
\bibitem{} Thuan, T.X. 1991, in Massive Stars in Starbursts, ed. C. Leitherer , N.R. Walborn,
T.M. Heckman, \& C.A. Norman (Cambridge: Cambridge Univ. Press), 183
\bibitem{} Thuan, T.X., \& Izotov, Y.I. 1997, ApJ, 489, 623
\bibitem{} Thuan T. X., Izotov Y. I., \& Lipovetsky V. A., 1995, ApJ, 445, 108
\bibitem{} Thuan, T. X., Izotov, Y. I., \& Lipovetsky, V. A. 1997, ApJ, 477,
 661
\bibitem{} Thuan, T. X., Sauvage, M, \& Madden, S. 1999, ApJ, May 10
\bibitem{} Thuan, T.X., Izotov, Y.I., \& Foltz, C.B. 1999, ApJ, in press
\bibitem{} Thuan, T.X., Izotov, Y.I., Lipovetsky, V.A. \& Pustilnik, S.A. 1994,
in Dwarf Galaxies, ed. G. Meylan \& P. Prugniel (Garching: European Southern Observatory),
421
\bibitem{} Tytler, D., Fan, X.-M., \& Burles, S. 1996, Nature, 381, 207
\bibitem{} White, D. \& Fabian, A. 1995, MNRAS, 273, 72
\bibitem{} Woosley S. E., \& Weaver T. A., 1995, ApJS, 101, 181
\bibitem{} Yee, H.K.C., Ellingston, E., Bechtold, J. Carlberg, R.G., \& Cuillandre, J.-C.
1996, AJ, 111, 1783
\end{bloisbib}
\vfill
\end{document}

Young Dwarf Galaxies and Cosmology


Doping evolution and polar surface reconstruction of the infinite-layer cuprate Sr$_{1-x}$La$_{x}$CuO$_{2}$

Abstract: We use angle-resolved photoemission spectroscopy to study the doping evolution of infinite-layer Sr$_{1-x}$La$_{x}$CuO$_{2}$ thin films grown by molecular-beam epitaxy. At low doping, the material exhibits a dispersive lower Hubbard band typical of the superconducting cuprate parent compounds. As carriers are added to the system, a continuous evolution from charge-transfer insulator to superconductor is observed, with the initial lower Hubbard band pinned well below the Fermi level and the development of a coherent low-energy band with electron doping. This two-component spectral function emphasizes the important role that strong local correlations play even at relatively high doping levels. Electron diffraction probes reveal a ${p(2\times2)}$ surface reconstruction of the material at low doping levels. Using a number of simple assumptions, we develop a model of this reconstruction based on the polar nature of the infinite-layer structure. Finally, we provide evidence for a thickness-controlled transition in ultrathin films of SrCuO$_2$ grown on nonpolar SrTiO$_3$, highlighting the diverse structural changes that can occur in polar complex oxide thin films.

\section{Introduction}

The generic aspects of the hole-doped side of the cuprate phase diagram have been firmly established by investigating a multitude of hole-doped material families. In contrast, there exist only two families from which most of our understanding of electron doping in the cuprates is derived: \textit{Re}$_{2-x}$Ce$_{x}$CuO$_{4}$, where \textit{Re} is a rare earth element, and ``infinite-layer'' Sr$_{1-x}$La$_x$CuO$_2$ \cite{armitage2010}. One conspicuous difference between electron and hole doping is the strength of $(\pi,\pi)$ antiferromagnetic order, which can persist up to an electron doping of $x = 0.14$ \cite{luke1990} and may even coexist with superconductivity in some materials \cite{yamada2003,alff2003,harter2012}. In all of the hole-doped cuprates, on the other hand, antiferromagnetism is rapidly suppressed by $x \approx 0.03$ and does not coexist with superconductivity. These facts highlight a clear asymmetry in the doping phase diagram of the cuprates, with important ramifications for theories of high-temperature superconductivity. Sr$_{1-x}$La$_x$CuO$_2$ is an ideal material for studying electron doping in the cuprates because of its uncomplicated structure, which consists of perfectly square and flat CuO$_2$ planes separated by simple Sr$_{1-x}$La$_x$ charge reservoir layers.

Bulk growth of Sr$_{1-x}$La$_x$CuO$_2$ is limited to polycrystalline samples due to the required high growth pressures. High quality thin films, however, can be grown by taking advantage of epitaxial stabilization \cite{tsukahara1991,norton1994,karimoto2001,karimoto2004,leca2006}. Indeed, within the last decade, complex oxide thin films have attracted increasing attention for the diverse electronic systems that they host \cite{heber2009,mannhart2010}. Polar structures with alternating charged atomic planes are particularly prevalent among the array of oxide structures commonly studied \cite{ohtomo2004,nakagawa2006,reyren2007}.  As thin films of these materials are grown, this alternating polarity leads to a thermodynamically unstable electric potential divergence---a so-called ``polar catastrophe''---which is prevented by a structural or electronic reconstruction. Undoped SrCuO$_2$, consisting of alternating Sr$^{2+}$ and (CuO$_2$)$^{2-}$ planes and no mobile charges to screen the electrostatic potential buildup, is such an example.

In this paper we use \textit{in situ} angle-resolved photoemission spectroscopy (ARPES) and electron diffraction to study both the doping evolution as well as the polar reconstruction in epitaxially-stabilized Sr$_{1-x}$La$_{x}$CuO$_{2}$ thin films grown by molecular-beam epitaxy. In Sec.\ \ref{sectionMethods}, we describe the film growth and experimental details. In Sec.\ \ref{sectionUndoped}, we present measurements of Sr$_{1-x}$La$_x$CuO$_2$ at low doping levels, showing a dispersive lower Hubbard band (LHB) characteristic of other parent cuprates. Section \ref{sectionDoping} demonstrates that with increased electron doping, a continuous evolution from insulator to superconductor occurs as spectral weight fills in the charge-transfer gap. In Sec.\ \ref{sectionReconstruction}, we describe electron diffraction probes that show evidence of a surface reconstruction consistent with the mitigation of a polar catastrophe in Sr$_{1-x}$La$_{x}$CuO$_{2}$. We also present evidence supporting a recent theoretical prediction of a thickness-controlled structural transition in ultrathin films of SrCuO$_2$ grown on nonpolar SrTiO$_3$. Finally, Sec.\ \ref{sectionConclusion} offers conclusions and implications of our work.

\begin{figure*}
\includegraphics{Figure1.pdf}
\caption{\label{figure1}{The low-energy electronic structure of Sr$_{0.99}$La$_{0.01}$CuO$_{2}$.  (a) Energy distribution curves (EDCs), offset for clarity, along a diagonal cut from $(0,0)$ to $(\pi,\pi)$ through the LHB, which is visible as a bump at the foot of the valence band. The bold red line shows the background EDC subtracted from the data in the remaining panels in order to enhance the LHB feature.  (b) Momentum space map of spectral weight at a binding energy of 0.5 eV showing the LHB at $(\pi/2,\pi/2)$ and equivalent points.  (c) EDC at $(\pi/2,\pi/2)$ after background subtraction. The peak has a Franck-Condon lineshape (demonstrated schematically in the lower left inset) and can be fit to a Gaussian (red curve). The ``foot'' in the low-binding-energy region deviates slightly from the Gaussian and most likely reflects low-energy levels occupied by the small amount of dopants added to the sample.  (d) Dispersion of the LHB as a function of momentum along Cut I from $(0,0)$ to $(\pi,\pi)$ (red points) and along perpendicular Cut II from $(\pi,0)$ to $(2\pi,-\pi)$ (blue points), as determined by Gaussian fitting. The smooth curves show the dispersion predicted by the $t$-$t'$-$t''$-$J$ model with $J = 150$ meV.  (e) Experimental spectrum along Cut I after background subtraction. The white dots reproduce the dispersion shown in panel (d).}}
\end{figure*}

\section{Experimental details}
\label{sectionMethods}

Sr$_{1-x}$La$_{x}$CuO$_{2}$ thin films ($x$ = 0, 0.01, 0.05, and 0.10) with a thickness of 60 unit cells (20 nm) and terminated with CuO$_{2}$ were deposited using a Veeco GEN10 dual-chamber oxide molecular-beam epitaxy system. Films were grown epitaxially on (110) GdScO$_{3}$ substrates, which have a distorted perovskite structure with a pseudocubic lattice constant of 3.968 \AA \cite{uecker2008}. Shuttered layer-by-layer deposition was performed in a background of 80\% pure distilled O$_{3}$ at a pressure of $1\times10^{-6}$ Torr and with a substrate growth temperature of 510 $^\circ$C. Depositions were monitored using reflection high-energy electron diffraction (RHEED). After growth, samples were reduced by vacuum annealing at 520 $^\circ$C for 30 minutes in order to eliminate excess oxygen atoms. The films were then cooled to 200 $^\circ$C before immediate transfer under ultra-high vacuum into the ARPES chamber. Samples with $x = 0.10$ were superconducting, exhibiting bulk resistance transitions in the range 25 $\pm$ 5 K. Further details of the film growth can be found in Ref.\ \citenum{maritato2013}.

ARPES measurements were performed with a VG Scienta R4000 electron spectrometer and He-I$\alpha$ photons (21.2 eV) at a base pressure of $7\times10^{-11}$ Torr and with an instrumental resolution better than $\Delta E$ = 20 meV and $\Delta k$ = 0.03 \AA$^{-1}$. The sample temperature was held at 200, 30, and 10 K for $x$ = 0.01, 0.05, and 0.10, respectively. The Fermi level $E_{\mathrm{F}}$ was determined by measuring polycrystalline gold in electrical contact with the sample. Experimental results were confirmed by studying multiple samples. After ARPES measurements, samples were characterized by \textit{in situ} low-energy electron diffraction (LEED) to examine surface structure and quality.  X-ray photoelectron spectroscopy was used to verify the stoichiometry of the films, measuring a difference of less than 0.01 between the measured and nominal lanthanum doping level $x$.

\begin{figure*}
\includegraphics{Figure2.pdf}
\caption{\label{figure2}{Momentum space evolution of spectral intensity with doping. Constant energy spectral maps for (a) $x = 0.01$, (b) $x = 0.05$, and (c) $x = 0.10$, integrated within $\pm$50 meV of the specified binding energy. At the Fermi level, the insulating $x = 0.01$ sample shows no spectral weight, while the $x = 0.05$ sample shows an accumulation of weight at $(\pi,0)$. By $x = 0.10$, the electron pocket at $(\pi,0)$ is well established and additional spectral weight is apparent at $(\pi/2,\pi/2)$. This weight is due to the finite integration region of the map rather than a true band crossing at the Fermi level \cite{harter2012}. At higher binding energies, there is clear evidence for a coexisting LHB with intensity near $(\pi/2,\pi/2)$ for all three doping levels. The circled regions at the top of panel (b) show shadow band reflections from spectral weight at $(\pi,0)$ due to a $p(2\times2)$ surface reconstruction in this sample, as discussed in Sec.\ \ref{sectionReconstruction}.}}
\end{figure*}

\section{Parent electronic structure}
\label{sectionUndoped}

The undoped parent compounds of the cuprates are charge-transfer insulators in which strong local Coulomb interactions dominate over a conventional band structure picture. Instead, the low-energy electronic structure is composed of a LHB, which typically has a bandwidth of $\sim$0.3 eV and a maximum at $(\pi/2,\pi/2)$ \cite{dagotto1994}. In Fig.\ \ref{figure1}, we show ARPES data for Sr$_{0.99}$La$_{0.01}$CuO$_{2}$, where $x = 0.01$ was intentionally added to prevent electrostatic charging of the sample (observed in stoichiometric SrCuO$_{2}$ films), pinning the chemical potential near the bottom of the upper Hubbard band. A small shoulder in the tail of the valence band is clearly present. After subtraction of a background EDC obtained by averaging the valence band tail near $(0,0)$, we observe a dispersive peak with a broad lineshape, characteristic of the LHB in other parent cuprates \cite{shen2004,armitage2002}. The spectral shape of the LHB is due to Franck-Condon broadening in which the coupling to a bosonic mode causes the spectral function to split into a set of discrete peaks. Each peak represents a resonance with a different boson occupation number, with the true quasiparticle pole residing in the low-binding-energy tail of intensity \cite{shen2004}. At $(\pi/2,\pi/2)$, the lineshape can be well-fit to a simple three-parameter Gaussian function, shown in Fig.\ \ref{figure1}(c), with the intensity maximum in the LHB at a binding energy of 0.81 eV and a full-width at half-maximum (FWHM) of 0.41 eV. This compares well with the LHB observed in other cuprate parent compounds, such as Ca$_2$CuO$_2$Cl$_2$ \cite{shen2004} and Nd$_2$CuO$_4$ \cite{armitage2002}, where the FWHM of the LHB is measured to be 0.34 eV and 0.36 eV, respectively.

As Fig.\ \ref{figure1}(d) shows, the LHB is dispersive, exhibiting a symmetric energy maximum at $(\pi/2,\pi/2)$. The $t$-$t'$-$t''$-$J$ model is often used to describe the motion of a doped carrier in a two-dimensional antiferromagnetic insulator, with the $t'$ and $t''$ parameters necessary to reproduce the clearly observed dispersion in the $(\pi,0)$ to $(0,\pi)$ direction. By fixing $J = 150$ meV and allowing $t'$ and $t''$ to vary in order to match the experimentally measured curvature along the $(0,0)$ to $(\pi,\pi)$ and transverse directions, we obtain a good fit to the data with $t' = -53$ meV and $t'' = 66$ meV. This is similar to the values $J = 140$ meV, $t' = -38$ meV, and $t'' = 22$ meV from Ref.\ \citenum{tohyama2000}, obtained by fitting to self-consistent Born approximation calculations of Sr$_2$CuO$_2$Cl$_2$, another undoped cuprate. The similarity suggests a universality of the electronic structure of the cuprates in the limit of low doping. While weakly interacting formalisms such as density functional theory predict a metallic state, the $t$-$t'$-$t''$-$J$ model in conjunction with Franck-Condon broadening can describe the observed electronic structure, accounting for the bandwidth change from $8t \sim 3$ eV to $2J \sim 0.3$ eV and the dispersion symmetry around $(\pi/2,\pi/2)$.

\section{Doping evolution of low-energy states}
\label{sectionDoping}

\begin{figure*}
\includegraphics{Figure3.pdf}
\caption{\label{figure3}{Evolution of electronic structure with doping. Spectra along $(0,0)$ to $(\pi,\pi)$ after background subtraction (a) for $x = 0.01$, (b) for $x = 0.05$, and (c) for $x = 0.10$. As the doping level increases, the spectral weight of the LHB shifts to lower binding energy and gradually fills in the charge-transfer gap. At $x = 0.10$, a coherent band near the Fermi level is visible.   (d) Momentum distribution curve (MDC)-derived dispersion for $x =$ 0.01, 0.05, and 0.10 (also shown as white lines in the preceding panels). The LHB maximum appears to shift away from $(\pi/2,\pi/2)$ with increased doping. The ``boomerang'' phenomenon is clearly visible for $x = 0.10$.  (e) Schematic diagram showing the qualitative form of the spectral function for $x = 0.10$. Spectral weight fills in the charge-transfer gap, forming a coherent band on top of the remnant LHB. The ``boomerang'' phenomenon in MDC-derived dispersions, shown by the black dashed line, is an artifact arising from the presence of two bands.}}
\end{figure*}

\begin{figure}
\includegraphics{Figure4.pdf}
\caption{\label{figure4}{Energy distribution curves.  (a) Doping dependence of EDCs at $(\pi/2,\pi/2)$ for $x =$ 0.01, 0.05, and 0.10. As the doping level increases, states are filled in near the Fermi level.  (b) An illustration of the two-component spectral function for $x = 0.10$. The EDC is made up of a coherent low-energy band and an incoherent high-energy remnant LHB.}}
\end{figure}

Upon addition of electrons into the CuO$_2$ planes via the substitution of trivalent lanthanum for divalent strontium in the intervening layers, the upper Hubbard band in Sr$_{1-x}$La$_{x}$CuO$_{2}$ evolves into a metallic Fermi surface with spectral weight first appearing at $(\pi,0)$ and equivalent points, as evidenced by the Fermi surface map of the $x = 0.05$ sample in Fig.\ \ref{figure2}(b). This behavior is identical to that seen in the electron-doped \textit{Re}$_{2-x}$Ce$_{x}$CuO$_{4}$ family \cite{armitage2002}. At $x = 0.10$, the Fermi surface consists solely of an electron pocket centered at $(\pi,0)$, although spectral weight is still apparent at $(\pi/2,\pi/2)$. As Ref.\ \citenum{harter2012} describes, this unusual Fermi surface is due to strong $(\pi,\pi)$ antiferromagnetism, which gaps the quasiparticles near $(\pi/2,\pi/2)$. Recent studies of the \textit{Re}$_{2-x}$Ce$_{x}$CuO$_{4}$ family have shown a clear link between oxygen content and the gapping of portions of the Fermi surface by antiferromagnetism \cite{richard2007,richard2008,song2012}. Chemical and structural similarities between the two families strongly suggests that the energy pseudogap near $(\pi/2,\pi/2)$ observed for Sr$_{0.90}$La$_{0.10}$CuO$_{2}$ may also be due to incomplete oxygen reduction. A further discussion of the oxygen reduction step can be found in Sec.\ \ref{sectionReconstruction}. Figure \ref{figure2} also shows maps at higher binding energies, where all three doping levels show evidence of a remnant LHB which takes the form of diffuse regions of spectral weight centered near $(\pi/2,\pi/2)$ and coexistent with the coherent bands dispersing through the Fermi level.

As Fig.\ \ref{figure3} shows, along $(0,0)$ to $(\pi,\pi)$ we observe that rather than closing abruptly, the charge-transfer gap gradually fills with spectral weight upon electron doping from 1\% to 10\%. This is accompanied by a qualitative change in the nature of the low-energy excitations of the system: the localized states of the LHB give way to a coherent itinerant band dispersing through the Fermi level. Figures \ref{figure3}(c,d) highlight a notable feature in the dispersion derived from a MDC analysis: at approximately 0.2 eV, the dispersion appears to ``boomerang'' backwards. This phenomenon is observed at multiple equivalent points in momentum space and is inconsistent with single-band physics. Instead, the effect likely arises from a two-component spectral function illustrated in Fig.\ \ref{figure4}: a coherent low-energy band forming the Fermi surface and a large contribution of incoherent spectral weight at higher energies derived from a remnant LHB that survives even at $x = 0.10$. This behavior differs from the hole-doped cuprates or the \textit{Re}$_{2-x}$Ce$_{x}$CuO$_{4}$ family, where so-called ``waterfalls'' are observed at higher binding energies \cite{graf2007,inosov2007,ikeda2009,moritz2009}. Recent theoretical calculations that take into account strong electron correlations have predicted such a coexistence of a low-energy band and an incoherent high-energy branch at 10\% electron doping \cite{senechal2005,weber2010}. Figure \ref{figure3}(d) also shows that the position of the LHB maximum appears to shift from $(\pi/2,\pi/2)$ towards $(0,0)$ with doping at a rate of approximately $8.5\times10^{-3}$ $(\pi/a)$/\%, which is quantitatively similar to the behavior observed for hole-doped Ca$_{2-x}$Na$_{x}$CuO$_{2}$Cl$_{2}$, where the LHB shifts with doping at a rate of $7.5\times10^{-3}$ $(\pi/a)$/\% \cite{shen2004}. Interestingly, the maximum shifts in the same direction for both electron and hole doping, counter to what would be expected from band structure calculations. The prominence of the remnant LHB in the experimental data highlights the important role that strong local electron correlations play in the electronic structure of Sr$_{1-x}$La$_{x}$CuO$_{2}$ even at relatively high doping levels.

We observe an interesting feature in the Fermi surface map of the $x = 0.05$ sample, as circled in Fig.\ \ref{figure2}(b): weak but finite spectral weight at $(0,0)$ and $(\pi,\pi)$. Bands at these locations in momentum space are not expected by tight-binding, density functional theory, or the $t$-$t'$-$t''$-$J$ model. Instead, it appears that the observed intensity is the result of a $p(2\times2)$ surface reconstruction of the sample, which causes shadows of the real spectral weight at $(\pi,0)$ to be reflected onto these locations in momentum space. At higher binding energies, evidence of this reconstruction is absent because reflections from intense regions of the remnant LHB fall onto each other. This surface reconstruction is described further in the next section.

\section{Polar surface reconstruction}
\label{sectionReconstruction}

\begin{figure}
\includegraphics{Figure5.pdf}
\caption{\label{figure5}{Structural change induced by the oxygen reduction step.  (a) RHEED image along the [100]$_p$ azimuth before the vacuum annealing step for an $x = 0.10$ film. White arrows highlight extra diffraction streaks present in all as-grown films.  (b) RHEED image after oxygen reduction for the same film. The extra RHEED streaks vanish during the annealing step.  (c) LEED image of an unannealed $x = 0.10$ film taken with 100 eV electrons. Circles show where Bragg peaks are located for films with the proper infinite-layer structure; some Bragg peaks are missing while a number of incommensurate peaks are visible.}}
\end{figure}

The metastable infinite-layer structure, which can only be grown in polycrystalline form under high pressure or as an epitaxially-stabilized thin film, lies near a manifold of other structural phases, such as the edge-sharing chain-type structure \cite{gambardella1992}. One consequence is an elevated sensitivity to oxygen stoichiometry.  For example, all as-grown Sr$_{1-x}$La$_{x}$CuO$_{2}$ films require a vacuum annealing step in order to eliminate excess oxygen and form the infinite-layer structure. This feature is shared by both families of electron-doped cuprates and is generally believed to be related to the absence of apical oxygen atoms in their respective crystal structures. As Fig.\ \ref{figure5} shows, RHEED patterns before and after the oxygen reduction step show a marked structural change, and LEED performed on an unreduced film shows a number of extra diffraction peaks incommensurate with the tetragonal infinite-layer reciprocal lattice.

\begin{figure}
\includegraphics{Figure6.pdf}
\caption{\label{figure6}{Evidence of surface reconstruction by electron diffraction.  RHEED images along the [100]$_p$ azimuth after the vacuum annealing step for Sr$_{1-x}$La$_{x}$CuO$_{2}$ films with (a) $x = 0.10$ and (b) $x = 0$. The latter image shows extra diffraction streaks consistent with a doubled lattice constant.  LEED images taken with 100 eV electrons for films with (c) $x = 0.10$ and (d) $x = 0.05$. The latter image again shows very clear evidence of a $p(2\times2)$ surface reconstruction causing a doubling of the unit cell in both the $a$ and $b$ directions.}}
\end{figure}

Even within the correct bulk structural phase, many nominally tetragonal transition metal oxides, such as Sr$_2$RuO$_4$ \cite{damascelli2000} and SrTiO$_3$ \cite{erdman2002}, are known to support surface reconstructions because of their complex surface chemistry. Electron diffraction probes are sensitive to such reconstructions. As demonstrated in Fig.\ \ref{figure6}, both RHEED, performed after growth at high temperature, and LEED, performed at low temperature, indicate that some Sr$_{1-x}$La$_{x}$CuO$_{2}$ samples with low doping levels show a $p(2\times2)$ surface reconstruction not observed by bulk x-ray diffraction. Within the ionic limit at low doping, the infinite-layer structure is intrinsically polar, alternating between charged Sr$^{2+}$ and (CuO$_2$)$^{2-}$ layers. This introduces a thermodynamic instability towards reconstruction as films of the material are grown. We argue below that the observed $p(2\times2)$ reconstruction is most likely a result of the polar surface of Sr$_{1-x}$La$_{x}$CuO$_{2}$, and that a polar catastrophe is avoided in this material by the formation of ordered oxygen vacancies on the topmost CuO$_2$ plane.

\begin{figure}
\includegraphics{Figure7.pdf}
\caption{\label{figure7}{Proposed model of surface reconstruction.  (a) One orientation of the surface structure composed of oxygen vacancies and consistent with a $p(2\times2)$ reconstruction. Three other orientations related by in-plane 90$^\circ$ rotations are also permitted. Upper left and lower right shaded squares shows the original unit cell and the doubled unit cell, respectively.  (b) Layer-by-layer view showing the proposed oxygen vacancy reconstruction of the terminal CuO$_2$ plane (highlighted in yellow). The vacancies result in a net $-1$ charge per unit cell on the topmost plane and a transfer of $-1$ charge per unit cell to the bottom of the film, avoiding the electric potential divergence associated with a polar catastrophe.}}
\end{figure}

The measured Sr$_{1-x}$La$_{x}$CuO$_{2}$ films are terminated with a CuO$_2$ layer, and the extra RHEED streaks associated with the reconstruction, shown in Fig.\ \ref{figure6}(b), become more intense during the vacuum annealing step that occurs at the end of film growth. This strongly suggests that the reconstruction is related to the removal of oxygen from the topmost CuO$_2$ plane. Indeed, half an oxygen vacancy per unit cell will change the net charge of the terminal CuO$_2$ atomic layer from $-2$ to $-1$, resulting in a cancellation of the divergent electric potential associated with alternating charged layers \cite{nakagawa2006}. Under this constraint, only one structure is consistent with a $p(2\times2)$ symmetry. One orientation of this structure is displayed in Fig.\ \ref{figure7}, with three others related by in-plane 90$^\circ$ rotations. The measured widths of both RHEED streaks and LEED spots indicate that such oxygen vacancy ordering persists over an in-plane length scale of at least $\sim$10 unit cells. We emphasize that the clear $(\pm1/2,\pm1/2)$ diffraction peaks in LEED measurements definitively rule out the possibility that the surface consists of domains of $p(2\times1)$ and $p(1\times2)$ reconstructions, a fact greatly limiting the set of possible oxygen vacancy structures to consider. If a single domain of our proposed model existed on the film surface, an anisotropic LEED structure factor would be observed. Because we instead observe four-fold symmetric diffraction patterns, the films likely contain domains of all four rotational orientations.

Higher doping levels lead to better metallic screening of the polar electric potential divergence, and the band gap of a material can play a large role in its polar reconstruction \cite{gu2009}. It is therefore natural to expect the tendency towards a polar surface reconstruction to be diminished at higher doping levels.  Indeed, we note that for $x \leq 0.05$, about half of films showed evidence of a reconstruction either by RHEED or LEED, while for $x \approx 0.10$, only one out of eight films showed the phenomenon. The observed film-to-film variability at fixed doping may be related to the fact that the formation of long-range structural order at the surface of the film, necessary for observing the reconstruction with diffraction probes, is likely sensitive to temperature, oxygen partial pressure, and other growth parameters.

\begin{figure}
\includegraphics{Figure8.pdf}
\caption{\label{figure8}{Fermi surface map for a Sr$_{0.95}$La$_{0.05}$CuO$_{2}$ film taken at 20 K showing spectral weight within $E_F \pm 50$ meV and normalized to a featureless background at high binding energy. The Fermi surface is composed of small electron pockets centered at $(\pi,0)$ and equivalent points. Weak spectral weight visible at $(0,0)$ and $(\pi,\pi)$ is a result of shadow band reflections (shown by yellow arrows) due to a $p(2\times2)$ surface reconstruction in this film.}}
\end{figure}

Despite a dramatic reconstruction of the terminal CuO$_2$ surface in Sr$_{1-x}$La$_{x}$CuO$_{2}$, ARPES measurements are mostly consistent with a pristine unreconstructed material, with the Fermi surface enclosing the correct Luttinger volume. Figure \ref{figure8} shows a Fermi surface map for an $x = 0.05$ sample, identical to the top panel of Fig.\ \ref{figure2}(b), where there is some evidence of a $p(2\times2)$ reconstruction in the form of weak shadow bands at $(0,0)$ and $(\pi,\pi)$. With this exception, the ARPES data can be analyzed without considering the reconstruction. If the model presented above is correct, the large number of oxygen vacancies on the film surface will naturally alter the valence states of the topmost copper and oxygen atoms, moving their energies away from the Fermi level. Low-energy photoemission on such films will therefore effectively probe the first buried CuO$_2$ plane in the material. This will slightly suppress the photoemission intensity near the Fermi level because photoelectrons must travel through the top atomic layer before leaving the sample, explaining the high relative background intensity that is observed in the photoemission data of all Sr$_{1-x}$La$_{x}$CuO$_{2}$ films studied.

\begin{figure}
\includegraphics{Figure9.pdf}
\caption{\label{figure9}{Oscillations of the intensity of the [10] diffraction rod (RHEED oscillations) at the start of film growth.  (a) Thickness-controlled transition for a SrCuO$_2$ film grown on nonpolar (001) SrTiO$_3$. Deposition of the first four unit cells results in a dramatically varying RHEED intensity (highlighted in gray). Stable oscillations are observed only after depositing the fourth unit cell.  (b) Stable oscillations for a film grown on (110) GdScO$_3$. Diagrams at right illustrate individual atomic layers for a two unit cell thick film, with arrows representing the growth direction.}}
\end{figure}

Recently, Ref.\ \citenum{zhong2012} predicted that ultrathin films of polar SrCuO$_2$ grown on nonpolar SrTiO$_3$ substrates would exhibit a thickness-controlled transition from a chain-type structure for $\le4$ unit cells to the infinite-layer structure for $\ge5$ unit cells. As Fig.\ \ref{figure9} shows, we observe evidence for such a transition in RHEED oscillations during the growth of the first few unit cells of SrCuO$_2$ on (001) SrTiO$_3$. The first four unit cells consistently show a different pattern of oscillations, and only after deposition of the fourth unit cell do the typical oscillations begin. Interestingly, when films are grown on polar (110) GdScO$_3$ substrates, this behavior is suppressed and the first four deposition periods show oscillations qualitatively similar to those at later times. This atomic reconstruction, occurring during the formation of the first few unit cells of SrCuO$_2$, is likely related to electrostatic instabilities of both the film and the substrate. This may explain why film growth on a polar substrate such as GdScO$_3$ results in qualitatively different RHEED oscillations even during deposition of the first atomic layers.

\section{Conclusion}
\label{sectionConclusion}

In conclusion, we have used \textit{in situ} ARPES in conjunction with molecular-beam epitaxy to study the doping evolution of thin films of the infinite-layer electron doped cuprate Sr$_{1-x}$La$_{x}$CuO$_{2}$. At low doping, a dispersive LHB characteristic of cuprate parent compounds is observed. With the addition of electron carriers, a transition from charge-transfer insulator to metallic superconductor is observed as spectral weight, first appearing at $(\pi,0)$, gradually fills in the charge-transfer gap. Most notably, we observe the coexistence of a remnant LHB with the coherent low-energy states, even for $x = 0.10$. Electron diffraction was used to study the $p(2\times2)$ polar surface reconstruction observed in low-doped films. This reconstruction can be explained with a simple model by considering the polar nature of Sr$_{1-x}$La$_{x}$CuO$_{2}$. Finally, we have presented evidence supporting the theoretical prediction of a thickness-controlled transition in ultrathin films of Sr$_{1-x}$La$_{x}$CuO$_{2}$ grown on nonpolar substrates. Our work shows that strong local correlations, forming a remnant LHB, remain important in the cuprates even at high electron-doping levels. Furthermore, it appears that in films of Sr$_{1-x}$La$_{x}$CuO$_{2}$ in the ionic limit at low doping, a polar catastrophe is avoided via a structural rather than an electronic reconstruction. Our observations highlight the need for deeper investigation into the stability and structural changes that occur in polar complex oxide thin films and their surfaces.

\section*{ACKNOWLEDGMENTS}

We are grateful to Haofei Wei for help with sample characterization. This work was supported by the Air Force Office of Scientific Research (Grant No.\ FA9550-11-1-0033) and the National Science Foundation through the MRSEC program (Grant No.\ DMR-1120296). L.M.\ was supported by the Army Research Office (Grant No.\ W911NF-09-1-0415), D. E. S. acknowledges support from the National Science Foundation under Grant No.\ DGE-0707428 and through the IGERT program under Grant No.\ DGE-0654193, and E. J. M. acknowledges support from an NSERC PGS.

Superconductivity (cond-mat.supr-con); Materials Science (cond-mat.mtrl-sci); Strongly Correlated Electrons (cond-mat.str-el)


A new diluted magnetic semiconductor: The half-metallic ferromagnet CoTi(1-x)FexSb

Abstract: Half-Heusler compounds with 18 valence electrons are semi-conducting. It will be shown that doping with electrons results in half-metallic ferromagnets, similar to the case of diluted semi-conductors. CoTiSb is known to be a semi-conducting Half-Heusler compound. Doping by Fe is expected to result in ferromagnetic order. It was found that Ti can be replaced by up to about 10% Fe while its crystal structure still remains C1b, which was proved by X-ray powder diffraction. SQUID magnetometry revealed a magnetic moment of 0.32 mB per unit cell at 5K.

\section{Introduction}
\label{IN} Within the last decade a new idea has revolutionised
electronic devices: To use the spin of electrons in addition to its
charge. This field of application is called spintronics \cite{Coe01}.
Half-metallic ferromagnets (HMF) seem to be a suitable class of
materials which meet all requirements of spintronics. One reason is
their exceptional electronic structure: They are metals in one spin
direction and semiconductors in the other. The Half-Heusler compound
NiMnSb was the first material being predicted by electronic structure
calculations to be a HMF \cite{GME83}.

Many attempts have been made to prepare semi-conducting compounds which
also have ferromagnetic properties. Mn-doped GaAs was considered to be
a suitable compound, but its Curie temperature is only about 150~K
\cite{EWC02} and thus still far away from being suitable for
application in electronic devices (see \cite{Mac05} for a recent
review). Heusler compounds exhibit a large variety of different
electronic and magnetic properties ranging from semi-conducting to
ferromagnetic. This work focuses on the search for suitable compounds
that bridge both semi-conducting and ferromagnetic properties. In
particular, the research concentrates on finding half-metallic
ferromagnetic Half-Heusler alloys with transition temperatures well
above room temperature.

\section{Results and discussion}

CoTi$_{1-x}$Fe$_x$Sb samples were prepared by arc melting of
stoichiometric amounts of the constituents in an argon atmosphere at
10$^{-4}$~mbar. The polycrystalline ingots were then annealed in an
evacuated quartz tube at 1000~K for 7~days.

\subsection{Electronic properties}

The electronic structure was calculated for pure and Fe doped CoTiSb in
order to examine their electronic and magnetic structure. Self
consistent calculations were performed by means of the
full-relativistic, spin polarised Korringa-Kohn-Rostocker (KKR) method
in combination with coherent potential approximation (CPA)
\cite{Ebe99}. As predicted by Tobola {\it et al} \cite{Tob00}, CoTiSb
turned out to be a semiconductor with a gap at the Fermi energy for
both spin directions. Partial replacement of titanium by iron (10\%)
(assuming that Fe is statistically distributed on the Ti positions),
has the result that the semiconductor is converted into a HMF: The DOS
at the Fermi energy is clearly different from zero for only one spin
direction, while remaining zero for the other. Measurements of the
conductivity revealed the metallic character of the Fe-doped samples in
contrast to the semi-conducting behaviour of the pure sample.

\subsection{Structural properties}

  \begin{figure}[t]
\includegraphics[scale =1.0]{balke_JMMM_cm-f1.eps}
\centering
\caption{X-ray diffraction pattern of CoTiSb and CoTi$_{0.9}$Fe$_{0.1}$Sb.
        The spectra were excited by Cu K$_\alpha$ radiation}
    \label{fig-1}
  \end{figure}

It was carefully proved that Fe occupies Ti positions and not an other
vacant site. For that purpose, the X-ray powder diffraction pattern of
the iron-substituted compounds were compared to that from pure CoTiSb.
It was observed that no additional diffraction reflexes appear up to
10\% Fe-doping, thus confirming that the structure remains the same as
of pure CoTiSb, that is C1$_b$. The diffraction pattern of the pure and
a 10\% Fe doped sample are compared in Fig.\ref{fig-1}. The lattice
parameter $a = 0.5883(4)$nm for CoTiSb changes slightly (0.1\%) if 10\%
of Ti is replaced by Fe. The R values for the best fit in the Rietveld
refinement are R$_i$ = 6\% and R$_p$ = 11\% confirming the high degree
of site order of the sample.

\subsection{Magnetic properties}

The magnetic properties were investigated by means of magnetometry
(SQUID). To determine the saturation magnetisation hysteresis loops of
the CoTi$_{1-x}$Fe$_x$Sb samples for $x=0.01$ , 0.02, 0.05 and 0.1 were
measured. Already the sample with only 1\% Fe exhibit ferromagnetism at
room temperature. The total spin magnetic moment of the samples can be
predicted from the generalised Slater-Pauling rule \cite{Kue84} to be:
$m=N_V-18$, where $N_V$ is the number of valence electrons per unit
cell.

CoTiSb carries no magnetic moment according to the Slater-Pauling rule
because it has overall 18 valence electrons. Doping Fe on the Ti
positions should lead to a magnetic moment. In the case of
CoTi$_{0.9}$Fe$_{0.1}$Sb one can predict from the Slater-Pauling rule a
saturation magnetisation of 0.4~$\mu_{B}$ per formula unit due to the
18.4 valence electrons of the compound. In Fig.\ref{fig-2}. the
hysteresis loops of CoTi$_{0.9}$Fe$_{0.1}$Sb at 5K, 300K and 768K are
shown. Ferromagnetic behaviour is still observed at 768K which points
to a very high Curie temperature. The saturation magnetisation at 768K
has still half of the value at 5K.
The calculated magnetic moments for the Fe-doped samples fit as well as
the experimental values to the values one expected from
the Slater-Pauling rule (see inset Fig.2).

For more details and other
properties of the series of CoTi$_{1-x}$Fe$_x$Sb with the Fe
concentration ranging from $x=0$ to 0.1 see Ref.\cite{KBF06}.

  \begin{figure}[t]
\includegraphics[scale =1.0]{balke_JMMM_cm-f2.eps}
\centering
\caption{Magnetic properties of CoTi$_{0.9}$Fe$_{0.1}$Sb. \newline
         The magnetisation as a function of the applied magnetic field
         at 5K, 300K and 768K is shown.
         Inset: The saturation magnetisation as a function von the Fe content is shown.
         Compared are the experimental results with the calculated values.}
    \label{fig-2}
  \end{figure}

\section{Summary}

In summary, it was shown that Fe substituted \\
CoTi$_{1-x}$Fe$_x$Sb is
a new diluted semiconductor with a high Curie temperature of above
770K. Self consistent band structure calculations predict the alloy to
exhibit half-metallic ferromagnetism. The future idea is to design
devices which consist of alternating semi-conducting CoTiSb and
half-metallic ferromagnetic CoTi$_{1-x}$Fe$_x$Sb layers. In that case,
one may be able to design a device with the semi-conducting and
ferromagnetic parts being of almost the same material.