On the Distance and Age of the Pulsar Wind Nebula 3C58

Abstract: There is a growing community of astronomers presenting evidence that the pulsar wind nebula 3C58 is much older than the connection with the historical supernova of A.D 1181 would indicate. Most of the strong evidence against a young age for 3C58 relies heavily on the assumed distance of 3.2 kpc determined with HI absorption measurements. I have revisited this distance determination based on new HI data from the Canadian Galactic Plane Survey and added newly determined distances to objects in the neighbourhood, which are based on direct measurements by trigonometric parallax. This leads to a new more reliable distance estimate of 2 kpc for 3C58 and makes the connection between the pulsar wind nebula and the historical event from A.D. 1181 once again much more compelling.

\section{Introduction}

Of all the supernova remnants (SNRs) which are linked to historically observed
supernova events the connection between the pulsar wind nebula (PWN) 3C\,58 and
the supernova explosion observed A.D.~1181 by Chinese and Japanese astronomers
is probably the most disputed one. \citet{step71} claimed that there is a high
probability for a connection between the Guest Star from A.D.~1181 and the
supernova explosion that created 3C\,58. This has been revisited and confirmed
by \citet{step02}. The main arguments are the length of the visibility of the
1181 event, indicating a supernova explosion, and the lack of any other supernova remnant
candidates. For such a recent supernova there should
be a visible SNR. Although it is very difficult to argue with this argument,
strong evidence against such a young age for 3C\,58 has been mounting
\citep[for a list of the main arguments see e.g.][Table 3]{fese08}.

Most arguments against a young age for 3C\,58 rely heavily on the
assumed distance of 3.2~kpc \citep{robe93}. This distance was determined
kinematically from \ion{H}{i} absorption measurements by comparing the
resulting radial velocity with a flat rotation curve for the Outer Galaxy.
However, in particular for Perseus arm objects, this
leads to a significant overestimate of the distance. A
spiral shock in the Perseus arm is ``pushing'' objects towards us,
giving them a higher negative radial velocity, which makes them appear to
be farther away \citep{robe72}. Examples for this effect on
distances for SNRs and PWNe can be found
in \citet{koth02b,koth03} and \citet{fost04}. SNRs in particular rely on
kinematic distance estimates, since their distances cannot be easily determined
directly or by related stars as in the case of \ion{H}{ii} regions.

In this short article I present a new distance estimate for 3C\,58 and
based on this new distance I will discuss and re-evaluate evidence presented in
the literature against the historical connection of 3C\,58. I will show that the
new distance changes the
characteristics of this PWN quite dramatically, which leads to a higher
probability for its historical connection.

\section{Observations}

The \ion{H}{i} data I used for the determination of the new absorption profile towards
3C\,58 were obtained with the Synthesis Telescope at the Dominion Radio
Astrophysical Observatory \citep[DRAO,][]{DRAO} as part of the Canadian Galactic Plane
Survey \citep[CGPS,][]{tayl03}. Single antenna data are incorporated into the
interferometer maps to ensure accurate representation from the largest
structures down to the resolution limit of about 1\arcmin. The low spatial
frequency data was drawn from the Low-Resolution DRAO Survey of the CGPS region,
which was observed with the 26-m radio telescope at DRAO \citep{higg00}. The
resolution in the final maps varies slightly across the final images close to
$1\arcmin \times 1\arcmin {\rm cosec(DEC)}$. At the centre of 3C\,58 the
resolution of the \ion{H}{i} data is $58\arcsec \times 65\arcsec$ at an angle of
$75\deg$ for the major axis (counter-clockwise from the Galactic longitude
axis). The RMS noise
is about 3~K $T_B$ in each velocity channel of width 0.82446~km~s$^{-1}$ at
a velocity resolution of
1.3~km~s$^{-1}$.

\section{Results}

\subsection{\ion{H}{i} Absorption and Emission Associated with 3C\,58}

The latest distance
determination for 3C\,58 by \ion{H}{i} absorption measurements resulted in a
systemic velocity of
$\sim -38$~km~s$^{-1}$ and a Perseus arm location \citep{robe93}. This was translated
with a flat rotation model for our Galaxy
and the IAU supported values for the Sun's Galacto-centric distance of
$R_\odot = 8.5$~kpc and the Sun's circular motion around the Galactic centre of
$v_\odot = 220$~km\,s$^{-1}$ to a distance of 3.2~kpc.

\begin{figure}[!ht]
\centerline{
  \includegraphics[bb = 100 263 475 620,width=0.75\textwidth,clip]
  {kothes_roland_f1.eps}}
 \caption{{\bf Top}: \ion{H}{i} absorption profiles of 3C\,58. The relative absorption
     of the radio continuum signal at 1420~MHz is displayed as a function of radial
     velocity $v_{LSR}$. Left: Absorption profiles of 3 different positions are
     displayed as indicated in the bottom right panel. Right: Profile
     calculated using all pixels on 3C\,58 with $T_b \ge 100$\,K. Each pixel was
     weighted by its intensity.
     {\bf Bottom}: Left: \ion{H}{i} emission profiles used to determine the
     absorption profiles in the top left. These were averaged
     over the positions marked in the right panel. Right: Continuum image of
     3C\,58 taken from the CGPS. The "on" and "off" positions used to calculate the
     \ion{H}{i} absorption profiles are indicated by solid and dashed circles,
     respectively.}
   \label{abs}
\end{figure}

\begin{figure}[!ht]
\centerline{
 \includegraphics[bb = 35 129 553 475,width=0.85\textwidth,clip]
 {kothes_roland_f2.eps}}
\caption{\ion{H}{1} channel maps of 3C\,58 taken at the peaks of the
absorption spectra at +1 (top left) -7 (top centre), -34 (bottom left), and
-41~km\,s$^{-1}$ (bottom centre). The right panels show 10 channels averaged over
the interarm velocities between the Local and Perseus arm (top) and 10 channels
averaged over the velocities just outside the last absorption feature.
The continuum emission of 3C\,58 is indicated
by the white contours. The black dashed lines represent the \ion{H}{i} contours at the
levels indicated by the labels of the colour bars.}
\label{abs2}
\end{figure}

I used the CGPS data to determine a new \ion{H}{i} absorption profile towards 3C\,58
to study not only the integrated absorption spectrum, but also the changes over
the PWN. The resulting profiles are shown in Fig.~\ref{abs}. \ion{H}{i} channel maps
taken at the peaks of the absorption spectrum and channels averaged over the
interarm areas
are displayed in Fig.~\ref{abs2}. The integrated
absorption spectrum is literally identical to the one published by
\citet{robe93}. However, I found that the weak peak, seen
at around $-41$~km~s$^{-1}$ is actually a real absorption feature. It is not
seen over the entire PWN, but only on the right hand side at position 3
(see Figs.~\ref{abs} and \ref{abs2}).
This is likely the reason why this feature was not obvious in previously
published integrated absorption spectra.

\begin{figure}[!ht]
\centerline{
 \includegraphics[width=0.85\textwidth]
 {kothes_roland_f3.eps}}
\caption{\ion{H}{1} channel maps of the area around 3C\,58. Each channel is
0.82446~km\,s$^{-1}$ wide and the centre velocity is indicated. The location
and emission structure of 3C\,58 are marked by the white contours.}
\label{hiemi}
\end{figure}

\citet{wall94} discussed the possible association of 3C\,58 with a large
interstellar bubble. An inspection of the CGPS
\ion{H}{i} data set reveals that this double shell structure is actually a very smooth
object (see Fig.~\ref{hiemi}). The bubble is visible from about
$-33$ to $-38$\,km\,s$^{-1}$ and 3C\,58 is located in the top left area of it.
If 3C\,58 was located within
this bubble, the last absorption feature at about $-41$\,km\,s$^{-1}$ is likely
produced by the shell of the bubble expanding towards us.

\subsection{The Systemic Velocity of 3C\,58}

\citet{robe93} found a Perseus arm location at a systemic velocity of about
$v_{LSR} = -38$\,km\,s$^{-1}$, which they translated to a distance of 3.2\,kpc.
This systemic velocity was confirmed with the discovery of the large bubble
by \citet{wall94}. The CGPS \ion{H}{i} data (Figs.~\ref{abs} and \ref{abs2}) show an
additional absorption feature at about $-41$\,km\,s$^{-1}$ and a systemic
velocity of about $-36$\,km\,s$^{-1}$ for the bubble. There is no absorption
at velocities beyond $-41$\,km\,s$^{-1}$. In Fig.\ref{abs2} the map averaged
around $-25$\,km\,s$^{-1}$, representative of the interarm between the Local arm and
the Perseus arm shows significant absorption, however, the map averaged over velocities
beyond $-41$\,km\,s$^{-1}$ is free of absorption even though this velocity range
shows brighter emission than the interarm around $-25$\,km\,s$^{-1}$. This confirms the
Perseus arm location for 3C\,58. Presumably 3C\,58 is
located inside the bubble shown in Fig.~\ref{hiemi} and the last absorption feature,
which is only partly visible in the absorption of 3C\,58 is caused by a part of that bubble moving
towards us. That would give this bubble
an expansion velocity of about 5\,km\,s$^{-1}$. Using the equations from
\citet{mcgriff02a} we calculate a dynamic age of about $3\times 10^6$\,yr
assuming it is the result of a single supernova explosion and $6\times 10^6$\,yr for
a stellar wind bubble.

\section{Discussion}

\subsection{The Distance to 3C\,58}

The method by \citet{fost06} was used to determine a distance-velocity relation in
the direction of 3C\,58 (Fig.~\ref{dvplot}). For a
systemic velocity of $-36$\,km\,s$^{-1}$ there are two possible distances: in the
Perseus arm shock at about 2\,kpc and beyond Perseus arm at about 2.5 - 2.8\,kpc.
Interstellar material is compressed in the spiral shock, forms molecules, and then
stars. It takes a long times -- several $10^7$ years -- to migrate beyond the Perseus arm.
Therefore, considering the dynamical age of the bubble and the high probability that
the progenitor star of 3C\,58 formed in the spiral shock, the likelihood of 3C\,58 being
at the farther distance is very low.

There is a second independent method to determine a distance to 3C\,58, to relate
the PWN to a nearby object with a more reliable distance estimate.
The W\,3/4/5 \ion{H}{ii} region
complex and the related SNR HB\,3 are just a few degrees away from 3C\,58. Both, HB\,3
and W\,3/4/5, have very similar systemic velocities \citep[e.g.][]{rout91}. They are all
located in the Perseus arm and for W\,3 the distance was determined with trigonometric
parallax to related masers to be $1.95\pm0.04$\,kpc and $2.04\pm0.07$\,kpc
by \citet{xu06} and \citet{hach06}, respectively. This nicely agrees with the 2\,kpc
I estimated for 3C\,58 from \ion{H}{i} absorption.

\subsection{The Age of 3C\,58}

There are a number of arguments against a relation between the PWN 3C\,58 and
the historical supernova explosion from A.D. 1181 as outlined in
\citet[Table 3]{fese08}. I will discuss the two major groups of these arguments below.

\begin{figure}[!ht]
\centerline{
 \includegraphics[bb = 22 18 616 585,width=0.50\textwidth,clip]
 {kothes_roland_f4.eps}}
\caption{Plot of the distance-velocity relation in the direction of 3C\,58
determined with the method by\citet{fost06}.}
\label{dvplot}
\end{figure}

\subsubsection{Expansion Studies}

One important aspect of expansion studies that is not always fully appreciated is the
fact that a supernova explosion is a one-time event. For any expansion
study, observed features related to the explosion could have been decelerated but not
accelerated. On the other hand, this
may not be true for synchrotron emitting filaments related to a PWN, because those have
a continuous source of energy in the central pulsar. Hence, the interpretation of
expansion studies of synchrotron filaments is not straightforward.

A comparison of the two major expansion studies in optical
\citep[$\Rightarrow$ age $t \approx 3500$\,yr]{fese08} and radio
\citep[$\Rightarrow t \approx 7000$\,yr]{biet06}
already implies significant deceleration for the radio structures relative to the optical
filaments unless the energy released by the pulsar significantly accelerated all the
supernova ejecta. A comparison of the typical $10^{51}$\,erg released in a supernova
and the approximate energy released by the pulsar since birth of $\sim 10^{48}$\,erg
\citep{chev04} negates that possibility.

The optical filaments are created by material accelerated through the supernova
explosion. Because those filaments could have been decelerated but not accelerated
a simple averaging of the expansion velocity of individual filaments would not
necessarily lead to a good age estimate unless the scatter is entirely produced
by uncertainties in the observations or systematic errors, which is certainly not the case in the study
by \citet{fese08}. The fastest filament should be taken, since it presumably shows the
lowest deceleration. There is quite
a wide spread in velocity in the optical study by \citet{fese08}, which again
indicates that a lot of the emitting material must have been
decelerated. Therefore the assumption of insignificant deceleration for the
optical filaments is invalid as well.
This significantly  weakens the case for a large age based on the optical and
radio expansion studies.

\subsubsection{PWN Evolution and Energy}

\citet{chev04} discussed theoretical calculations of the evolutionary path of a PWN and
compared their results for swept up mass and internal energy with those from
observations. Both of these results rely heavily on the assumed distance.
\citet{bocc01} determined from their X-ray observations a swept-up mass
of $M_{sw} = 0.1 d_{3.2}^{2.5}$\,M$_\odot$, assuming a radius $R$ of $2.5\arcmin$, which results in 0.1\,M$_\odot$ for a
distance $d$ of 3.2\,kpc. The theoretical value from \citet{chev04} equates to:
$M_{sw} = \dot{E} R^{-2} t^3$ ($\dot{E}$: rotational energy loss rate of the pulsar)
resulting in 0.005\,M$_\odot$. Even with the relatively high
uncertainty of this kind of calculation the large discrepancy suggests a much larger
age for 3C\,58. To increase the theoretical value to the observed one the age of the
PWN has to be about $t \approx 2300$\,yr. The other argument is that the minimum energy
required to produce the synchrotron nebula is about $10^{48}$\,erg from equipartition
considerations. This value depends on the radio luminosity and thereby on
$d^2$. The total energy released from the pulsar into the nebula, however, which
can be approximated by $\dot{E} t$ is about $0.7\times 10^{48}$\,erg, which is much smaller.
This value is distance independent.

For a distance of 2\,kpc the minimum energy required to produce the synchrotron emission
decreases to $0.4\times 10^{48}$\,erg, which is now lower than the total energy released by
the pulsar, which remains unchanged. The values for the swept up masses equate to 0.030\,M$_\odot$ and
0.013\,M$_\odot$ for the X-ray observations and the theoretical calculations, respectively.
The two results for the swept-up mass still differ but not by much. To get those values to
be equal requires either an age of about 1100\,yr or a distance of 1.7\,kpc. For
a distance of 1.9\,kpc the required age would be reduced to 1000\,yr.

\section{Conclusion}

I derived a new more reliable distance of 2\,kpc to the PWN 3C\,58, by means of
\ion{H}{i} absorption in combination with a newly determined distance-velocity relation
and by relating this PWN to a nearby \ion{H}{ii} region SNR complex of well known
distance. This new distance changes many characteristics of this PWN quite
dramatically which once again makes the connection between 3C\,58 and the
historical supernova explosion from 1181 A.D. much more likely.

\acknowledgements

I would like to thank Tyler Foster for providing me with a distance-velocity
diagram in the direction of 3C\,58 using his modeling of Galactic hydrogen
distribution. The Dominion Radio Astrophysical Observatory is a National Facility
operated by the National Research Council. The Canadian Galactic
Plane Survey is a Canadian project with international partners, and is
supported by the Natural Sciences and Engineering Research Council
(NSERC).


Globular Clusters in Dwarf Elliptical Galaxies

Abstract: Globular clusters (GCs) provide one way of determining when the most important periods of star formation in dwarf elliptical (dE) galaxies occurred. Thus, they are also useful for comparing the star formation histories of dwarf and giant galaxies. We present results on the GC systems of 24 nearby dE galaxies imaged with the Hubble Space Telescope. We find that the specific GC frequency, S_N, increases with M_V. Also, S_N in dEs is more like the values found in giant ellipticals rather than in spirals, and the mean S_N for nucleated dEs is roughly a factor of 2 higher than for non-nucleated dEs. The colors suggest that most of the GCs formed before the bulk of the stars in the galaxies.

\section{Introduction}

Dwarf elliptical (dE) galaxies make up most of the objects in the
steeply rising, low-luminosity end of the galaxy luminosity function
in galaxy clusters \cite{t98}\cite{jpsp}.  Besides being the most
numerous type of galaxy in the universe, in hierarchical clustering
scenarios of galaxy formation progenitors of dEs are the building
blocks of massive galaxies.  Therefore, solutions to the questions of
when dEs formed and how they evolve have important cosmological
implications.  Unfortunately, there is still little consensus on these
issues or even about the connections between different types of dwarf
galaxies (see \cite{fb94} for a review).

The globular cluster systems of dE galaxies may provide an observational
probe of their most important episodes of star formation.  Globular
clusters (GCs) are thought to form under the extreme pressures produced by
very vigorous star formation \cite{ee97}\cite{hhm98}.  Since GCs have
instantaneous burst stellar populations, their integrated light contains
information about their ages and metallicities. Also, the specific globular
cluster frequency ($S_N = N_{c}10^{0.4(M_V+15)}$)
is dependent on Hubble type; it is a factor of 2--4
higher for giant ellipticals than for spirals and irregulars \cite{h91}.
This suggests a test for whether dEs have had star formation histories more
like giant ellipticals or to irregulars.  If they are related to Es, then
they should have $S_N \sim 2-5$, if they are more like late-type disks then
one would expect $S_N\sim1$.

To test this hypothesis and study the properties of the GCs and nuclei of
dE galaxies, we have begun a {\it Hubble Space Telescope (HST)} snapshot
survey of dEs in nearby groups and clusters.  The original sample was split
evenly between nucleated (dE,N) and non-nucleated (dE,noN) galaxies.
In {\it HST} Cycle 6 we
obtained observations of 24 galaxies in the Virgo and Fornax Clusters and
the Leo Group that are suitable for measuring $S_N$.  We are able to measure
($V\!-\!I$) colors of GCs down to $V\approx25$ \cite{de6}.

\section{Results}

The main results of this survey at present are summarized below.

\begin{itemize}

\item {\bf Nuclei:} The high resolution of {\it HST} allows us to
refine the classification of the dEs into nucleated and non-nucleated
categories.  For example, FCC~324 was classified as a non-nucleated dS0
from ground-based photographic material, but now clearly shows a
nucleus.  VCC~9 had been classified as nucleated but the ``nucleus''
is offset from the centroid of the background light by about 160~pc.
Therefore, we classify it as non-nucleated.  Three other galaxies do
not meet our criteria for being nucleated but their more luminous GCs
are within $\sim60$~pc of the galaxies' centroids.  These latter four
galaxies may have {\it offset nuclei} and they are shown in
Figure~1.

\begin{figure}
\epsscale{0.58}
\plotone{fig1.ps}
\caption{{\it HST} snapshots of four dE galaxies with possible offset
nuclei.  Each image shows a physical region about 3~kpc across.}
\end{figure}

\item {\bf Luminosity Function:} The luminosity function (LF) of the
GC candidates from the entire sample is consistent with the standard
Gaussian-shaped GCLF with a peak at $M_V^0=-7.4$ and a width
$\sigma_V\approx1$~mag.  Also, with a mean color of $\langle
V\!-\!I\rangle\approx0.95$, most of the GCs appear to be
old and metal-poor ($[{\rm Fe/H}]\approx-1.5$) like GCs in the Galaxy
and in nearby giant ellipticals.

\item {\bf Metallicities:} If differences in $(V\!-\!I)$ are due to
metallicity, then in general we find
\begin{displaymath}
[{\rm Fe/H}]_{GC} <  [{\rm Fe/H}]_{nuc} <  [{\rm Fe/H}]_{gal}.
\end{displaymath}
The GCs are between 0.0 and 0.4 dex more {\it metal-poor} than the
stars in bodies of the galaxies.  The GCs in a couple dEs have large
color spreads that imply metallicity spreads of up to 1 dex.  The most
luminous nuclei are $\sim0.3$ dex more {\it metal-rich} than the GCs
in these galaxies.

\item {\bf Trends in $S_N$:} There is a trend of increasing $S_N$ with
$M_V$ (Figure~2). Also, there is evidence for a difference between
dE,N and dE,noN galaxies.  The dE,Ns define the upper envelope of this
trend (solid curve).  The trend for the dE,noN is much shallower, it
is basically flat (dashed curve).  The galaxies with offset nuclei
shown in Figure~1 seem to have $S_N$ more like the nucleated galaxies.

\begin{figure}
\plotfiddle{fig2.eps}{8.7cm}{-90}{52}{52}{-210}{285}
\caption{$\log(S_N)$ vs. $M_V$ for all dE galaxies with measured $S_N$.
For galaxies with $S_N=0$, $\log(S_N)$ has been set to 0.0.
The circular symbols are from the {\it HST} snapshot sample \cite{de6}.
Triangles are from the ground-based work of Durrell et al. \cite{d96}
and squares are galaxies in the Local Group (LG).  The triangle labeled
``LG dSph'' represents the $\sim12$ LG dwarf spheroidals (non-nucleated)
that have $S_N=0$. The lines are weighted least-squares fits.}
\end{figure}

\item {\bf Averages for $S_N$:} The differences in $S_N$ suggested by
Figure~2 are also reflected in the average values of $S_N$ for the two types
are given below.  These averages include dE galaxies in the Local Group and
Durrell et al.'s \cite{d96} ground-based sample with $M_V<-13.5$.  The
uncertainties are standard deviations of the mean.

\begin{center}
\begin{tabular}{lc}
Type & $\langle S_N\rangle$ \\ \hline
dE,N & $6.4\pm1.2$ \\
dE,noN & $3.3\pm0.5$ \\ \hline
\end{tabular}
\end{center}

\end{itemize}

\section{Discussion}

The basic result of this work is that
the mean values of $S_N$ in dEs is in the range 3--6, very similar to
values for $S_N$ in giant ellipticals in rich clusters \cite{h91}.  Thus it
would seem that most dEs, and especially the dE,Ns, have experienced similar
star formation histories as giant Es.  The colors of the GCs suggest that
most of the
cluster formation occurred very early out of relatively metal-poor gas.
Subsequent star formation from more enriched gas produced most of the
stars in the bodies of the galaxies.

The differences in $S_N$ between dE,Ns and dE,noNs may be a result of
environment.  Nucleated dEs are preferentially found in the centers of
rich clusters and have spatial distributions like giant Es \cite{fs89}.
Thus, they would have formed in a high pressure, turbulent environment
that would be conducive to GC formation \cite{ee97}.  Gravitational torques
from interactions could drive gas towards
the centers and help with the formation of nuclei.
On the other hand, dE,noNs are more uniformly distributed, like the
spirals and irregulars \cite{fs89}.  Thus, they may have formed in
lower-pressure environments less suited to the formation of GCs or nuclei.
Numerical simulations would be useful for testing this scenario.
Finally, some dE,noNs have $S_N=0$, making themselves candidates for
being stripped irregulars.

The trend of $S_N$ with $M_V$ is more difficult to explain.  One possibility
is that dynamical friction is causing many of the more massive GCs to
conglomerate into the nucleus.  The dynamical friction timescales are short
enough for this to occur.  A second
option is that the most luminous dEs are stripped spirals; the bulges
of spirals do have similarities to dE,Ns \cite{c97}.

This work has shown the usefulness of using GCs to compare the star
formation histories of dwarf galaxies and giant galaxies.  However,
several issues still need to be resolved.  First, it is not clear whether the
trends in $S_N$ continue to fainter magnitudes or whether there
is just an increase in scatter.  Therefore, in {\it HST} Cycle~7 we
will image up to 30 more faint dEs.  Also, the mean value and scatter
in $S_N$ for dwarf irregulars (dIs) is very uncertain since very few
dIs have $S_N$ measurements.  Therefore, we will also image up to 30
dIs in Cycle~7.  Hopefully, this will provide a clearer picture of
the connection between different types of dwarf galaxies.

\acknowledgements

B.W.M. would like to thank the Leids Kerkhoven Bosscha Fonds and the
Netherlands Foundation for Research in Astronomy (ASTRON) for
financial support to attend this meeting.

\begin{moriondbib}
\bibitem{c97} Carollo, C. M., Stiavelli, M., de~Zeeuw P. T., \& Mack J.
1997, \aj {114} {2366}
\bibitem{d96} Durrell, P., Harris, W. E., Geisler, D., \& Pudritz R. 1996,
\aj {112} {972}
\bibitem{ee97} Elmegreen, B. G., \& Efremov, Y. N. 1997, \apj {480} {235}
\bibitem{fb94} Ferguson, H. C., \& Binggeli, B. 1994, A\&AR {6} {67}
\bibitem{fs89} Ferguson, H. C., \& Sandage, A. 1989, \apj {346} {L53}
\bibitem{h91} Harris, W. E. 1991, ARA\&A {29} {543}
\bibitem{hhm98} Harris, W. E., Harris G. L. S., \& McLaughlin, D. E. 1998,
{\it Astron. J.} in press
\bibitem{jpsp} Jones, J. B., et al. 1998, this volume (and astro-ph/9805287)
\bibitem{de6} Miller, B. W., Lotz, J., Ferguson, H. C., Stiavelli, M., \&
Whitmore, B. C. 1998, in preparation
\bibitem{t98} Trentham, N. 1998, this volume (and astro-ph/9804013)
\end{moriondbib}
\vfill
\end{document}


1/8 Doping Anomalies and Oxygen Vacancies in Underdoped Superconducting Cuprates

Abstract: Measurements of thermal conductivity ($\kappa$) versus temperature and doping for several cuprate superconductors are discussed. Suppressed values of the normal-state $\kappa$ and the slope change in $\kappa$ at $T_c$, observed near 1/8 doping in both YBa$_2$Cu$_3$O$_{6+x}$ and Hg cuprates, are attributed to local lattice distortions and the suppression of the superconducting condensate, respectively. Both phenomena are proposed to arise from small domains of localized planar holes, presumably a manifestation of phase separation. It is suggested that the phase behavior of $\kappa$ reflects stripe dynamics and the 1/8 doping anomalies stripe pinning by oxygen-vacancy clusters.

\section*{Introduction}
There is mounting experimental evidence that novel charge- and spin-ordered
phases are generic to underdoped cuprates\cite{Cuprates}. The doped holes in
Nd-doped La$_{1-x}$Sr$_x$CuO$_4$ segregate into periodic stripes that separate
antiferromagnetically ordered, hole-poor domains.  Lattice distortions pin
these stripes, and this pinning is most effective for planar hole
concentrations near $p$=1/8 where the stripe modulation wavelength is
commensurate with the lattice.  In the absence of pinning, stripe modulations
are presumed to be fluctuating and/or disordered, and this is the emerging
picture for La$_{1-x}$Sr$_x$CuO$_4$ (La-214) and YBa$_2$Cu$_3$O$_{6+x}$
(Y-123)\cite{DynamicStripes}.  We have recently demonstrated, through
measurements of thermal conductivity ($\kappa$)\cite{CohnHg}, that both Y-123
and HgBa$_2$Ca$_{m-1}$Cu$_m$O$_{2m+2+\delta}$ [Hg-12($m$-1)$m$, $m$=1, 2, 3]
exhibit doping anomalies near $p$=1/8 that can be attributed to the presence
of localized charge and associated lattice distortions.  Here we discuss the
oxygen doping behavior of the anomalies in Hg cuprates and present new results
on vacancy-free Ca-doped YBa$_2$Cu$_4$O$_8$ (Y-124) that suggest elastic
properties may reflect stripe dynamics and stripe fragments may pin near
oxygen-vacancy clusters.
\begin{figure}[t!]
\centerline{\epsfig{file=hts99_f1.psc,height=7.52in,width=5.8125in}}%
\vskip -3.4in
\caption{(a) Thermal resistivity at $T$=100K
relative to that at $p_{opt}=0.16$ for Y-123 polycrystals and the
ab-plane of single crystals \protect\cite{Popoviciu} [$W_{opt}=0.21$ mK/W
(0.08mK/W) for the polycrystal (crystals)]. The hole concentration
per planar Cu atom, $p$, was determined from thermopower
measurements and the empirical relations with bond valence sums
59 established by Tallon {\it et al.}\protect\cite{PfromTEP}.  Also plotted
($\times$'s) is the relative intensity of anomalous $^{63}$Cu NQR
signals\protect\cite{HamScal}. The dashed curve is a guide to the eye. (b)
The normalized slope change in $\kappa(T)$ at $T_c$ {\it vs}
doping for each of the Y-123 specimens in (a), for the a-axis of
single-crystal Y-124\protect\cite{Cohn124xtals}, and for polycrystal
Y$_{1-x}$Ca$_x$-124\protect\cite{CohnYCa124}. The Y$_{1-x}$Ca$_x$-124 data
are multiplied by 2.3 so that they coincide with the Y-124 crystal
data at x=0. Solid (open) symbols are referred to the left (right)
ordinate. Also shown ($\times$'s) are the normalized electronic
specific heat jump\protect\cite{Loram}, $\Delta\gamma/\gamma$, and (solid
curve) the $\mu$SR depolarization rate\protect\cite{muSRY123} (in
$\mu$s$^{-1}$), divided by 1.4 and 2.7, respectively, and
referred to the right ordinate. (c) Thermal resistivity for
Hg-1201 and Hg-1212 polycrystals relative to that of
Hg-1223\protect\cite{CohnHg}. Dashed curves are guides. (d) The normalized
slope change in $\kappa(T)$ at $T_c$ {\it vs} doping for Hg
cuprates.  The solid line is $1.71-250(p-0.157)^2$.  Dashed curves
are guides.}
\label{fig1}
\end{figure}

The $p$=1/8 features (Fig.~\ref{fig1}) are evident in the doping behavior of
the normal-state thermal resistivity, $W=1/\kappa$, and the normalized change
in temperature derivative of $\kappa$ that occurs at $T_c$, $\Gamma\equiv
-d(\kappa^s/\kappa^n)/dt|_{t=1}$ [$t=T/T_c$ and $\kappa^s (\kappa^n)$ is the
thermal conductivity in the superconducting (normal) state].  For Y-123, there
is a compelling correlation of $W(p)$ and $\Gamma(p)$ with the doping behavior
of anomalous $^{63}$Cu NQR spectral weight\cite{HamScal}, attributed to
localized holes, and the electronic specific heat jump\cite{Loram},
$\Delta\gamma/\gamma$, respectively [crosses in Fig.'s 1 (a) and (b)]. Thus
$W$ probes lattice distortions associated with localized holes, and $\Gamma$
the change in low-energy spectral weight induced by superconductivity.

Since both $\Gamma$ and $\Delta\gamma/\gamma$ provide bulk measures of the
superfluid volume, it is significant that the muon spin rotation ($\mu$SR)
depolarization rate [$\sigma_0$ in Fig. 1 (b)], proportional to the superfluid
density, exhibits no anomalous behavior near 1/8 doping.  The $\mu$SR signal
originates in regions of the specimen where there is a flux lattice. The
apparent discrepancy between $\sigma_0$ and $\Gamma$ or $\Delta\gamma/\gamma$
is resolved if the material is inhomogeneous, composed of non-superconducting
clusters embedded in a superconducting network.  The suppression of $\Gamma$
and $\Delta\gamma/\gamma$ below the scaled $\sigma_0$ curve in Fig. 1 (a) are
then measures of the non-superconducting volume fraction. Taken together, the
$W$ and $\Gamma$ data imply that the non-superconducting regions are comprised
of localized holes and associated lattice distortions, i.e. polarons.  It is
plausible that these hole-localized regions are stripe domains akin to those
inferred from neutron scattering\cite{Cuprates}.  That they produce lattice
thermal resistance implies they are static on the timescale of the average
phonon lifetime, estimated as a few ps in the ab-plane of Y-123\cite{CohnHg}.
Given that $T_c$ is not substantially suppressed near $p$=1/8 implies that
these regions do not percolate, and thus the picture is one of small clusters
of localized holes and their associated lattice distortions, perhaps no larger
than the stripe unit cell ($2a\times8a$ where $a$ is the lattice
constant\cite{Cuprates}), separated by a distance comparable to the phonon
mean free path (about 100${\rm\AA}\sim 25a$ at 100K).

For the Hg materials the $p$=1/8 enhancement of $W$ and suppression of
$\Gamma$ is most prominent in single-layer Hg-1201, less so in double-layer
Hg-1212, and absent or negligible in three-layer Hg-1223.  This trend follows
that of the oxygen vacancy concentration: a single HgO$_{\delta}$ layer per
unit cell contributes charge to $m$ planes in Hg-12($m$-1)$m$ so that the
oxygen vacancy concentration, $1-\delta$, increases with decreasing $m$
\cite{Occupancy}.  The absence of suppression in $\Gamma$ near 1/8 doping for
Hg-1223 [Fig.~1~(d)] suggests that this material has sufficiently few
localized-hole domains that their effects in $W$ and $\Gamma$ are
unobservable.  Thus we employ the Hg-1223 $W(p)$ data as a reference and plot
the differences for the other two compounds in Fig. 1 (c). Comparing
Fig.'s~1~(c) and (d) we see that for both Hg-1201 and Hg-1212 $W$ is enhanced
and $\Gamma$ is suppressed relative to values for Hg-1223 in common ranges of
$p$, with maximal differences near $p$=1/8.

\begin{figure}[t!]
\vskip -.5in \centerline{\epsfig{file=hts99_f2.psc,height=7.2in}}
\vskip -3in
\caption{(a) Thermal conductivity of Hg-1223 polycrystal versus doping.
The dotted and dashed curves are proposed electronic (right ordinate) and
lattice thermal conductivities, respectively.  The solid line is a guide to
the eye. (b) sound velocity data for single-crystal La-214 from
Ref.~\protect\cite{LSCOelastic}.} \label{fig2}
\end{figure}
\section*{Stripes Alter the phonon dispersion?}

That 1/8 doping anomalies are observed in both Y-123 and Hg cuprates suggests
they are a generic feature of underdoped CuO$_2$ planes.  In this context the
doping behavior of $\kappa$ for Hg-1223 [Fig. 2 (a)] is of interest since it
presumably approximates the phase behavior in the absence of localized holes.
We expect the electronic contribution, $\kappa_e$, to increase smoothly with
increasing $p$. The upturn in $\kappa$ above optimal doping ($p$=0.16) is
presumably attributable to this rising $\kappa_e$.  The most reliable estimate
of $\kappa_e$ in the cuprates comes from thermal Hall conductivity
measurements\cite{Krishana} on Y-123 which imply $\kappa_e/\kappa\simeq 0.1$
for in-plane heat conduction near optimal doping.  The data in Fig. 1 (a), (b)
suggest that this ratio is roughly the same in polycrystals, thus motivating
the dotted and dashed curves in Fig. 2 (a) as an educated guess for the
electronic and lattice terms, respectively, in Hg-1223. The lattice
conductivity, $\kappa_L=\kappa-\kappa_e$, predominates in the underdoped
regime and is peaked near $p$=1/8.

There is experimental support for the proposal that this peak in $\kappa$ is
associated with doping-dependent changes in the phonon dispersion which are
generic to underdoped cuprates.  Figure 2~(b) shows the behavior of the
normal-state, transverse shear elastic constant (proportional to the square of
the sound velocity) for single-crystal La-214\cite{LSCOelastic}, the only
material for which doping-dependent measurements for all the main symmetry
directions have been reported to our knowledge. A substantial hardening of the
lattice in the underdoped regime, with a maximum near $p=x\simeq$1/8, was
observed for all symmetries, indicating that the changes with doping are
systemic.

It is possible that the phase behavior of the elastic constants for La-214 and
$\kappa_L$ for Hg-1223 reflect a renormalization of the lattice dispersion due
to changes in stripe dynamics with doping. At $p$=1/8 the stripes are
maximally commensurate with the lattice, and their fluctuations should be
minimal. Enhanced fluctuations are to be expected at higher doping, due to the
destabilizing role of repulsive interactions in stripes that neutron
scattering results (Yamada {\it et al.}\cite{DynamicStripes}) suggest, are
charge compressed.  For $p<$1/8, the cluster spin-glass state observed by
$\mu$SR studies\cite{Niedermayer} is characterized by increasing magnetic
disorder with decreasing $p$ in the range $0.08 < p <0.12$, possibly
associated with increasing disorder\cite{Emery} in the stripe period.  It is
plausible that disorder in the stripe system at both higher and lower doping
about $p$=1/8 induces a softening of the lattice that is reflected in Fig.~2.
Theoretical investigations of elastic coupling to the stripe system would
certainly be of interest.

\begin{figure}[b!]
\centerline{\epsfig{file=hts99_f3.psc,height=6.47in,width=5in}}
\vskip -3in \caption{(a) The enhancement in thermal resistivity
and suppression in $\Gamma$ for the Hg compounds at $p$=1/8,
relative to values for Hg-1223, plotted versus the oscillator
strength ratio of Raman-scattering defect modes
(Ref.~\protect\cite{HgRaman}).} \label{fig3}
\end{figure}
\section*{Oxygen Vacancies and 1/8 Anomalies}

Returning now to the enhancement of $W$ and suppression of $\Gamma$ near
$p$=1/8, our measurements suggest these phenomena reflect stripe pinning (or
reduced stripe fluctuations) associated with oxygen vacancy clusters.  Raman
studies of defect modes in the Hg materials\cite{HgRaman} provide a measure of
the density of vacancy clusters. The 590 cm$^{-1}$ Raman mode in Hg cuprates
is attributed to c-axis vibrations of apical oxygen in the presence of oxygen
vacancies on each of the four nearest-neighbor sites in the HgO$_{\delta}$
layers.  The amount by which the thermal resistivity and $\Gamma$ values of
Hg-1201 and Hg-1212 differ from those of Hg-1223 at $p$=1/8 both correlate
well with the integrated oscillator strength of this vibrational mode
normalized by that of the 570 cm$^{-1}$ mode (Fig.~\ref{fig3}), the latter
attributed to apical vibrations in the presence of fewer than four vacancies,
and common in the spectra\cite{HgRaman} of all three Hg materials.  We infer
that oxygen-vacancy clusters, of size four or more in the Hg cuprates, can pin
a stripe fragment.

Supporting the role of oxygen vacancies in the observed 1/8 features for Y-123
are our recent measurements on oxygen-stoichiometric
Y$_{1-x}$Ca$_x$-124\cite{CohnYCa124}. The $\Gamma$ values for $x$=0, 0.10,
0.15 are plotted in Fig.~1~(b).  There is no evidence for suppression of
$\Gamma$ near 1/8 doping, and the data follow the scaled $\mu$SR curve quite
well.  Within the context of the interpretation we have outlined, we conclude
that Ca substitution for Y (at the level of $\leq$ 15\%) is not as effective
in stripe pinning as are oxygen vacancy clusters.  It remains to be determined
by what mechanism these clusters induce pinning.

\section*{acknowledgements}
This work was supported by NSF Grant No. DMR-9631236.


Small-Scale Mixing of the Elements

Abstract: The progressive mixing and contamination of interstellar gas by supernovae and other processes following the passage of a spiral density wave is reviewed, with an emphasis on the Solar neighborhood. Regions of star formation should begin their lives with an inhomogeneous mixture of abundances as a result of their chaotic and large-scale formation processes. These inhomogeneities should continue to increase during several generations of star formation until the gas enters the interarm region. Then cloud dispersal by large interarm tidal forces, high rates of shear, and internal star formation should all lead to cloud-to-cloud mixing, while increased ionization, heating, and evaporation should lead to mixing on the atomic level. Gradients in the dispersion of elemental abundances are expected for galaxy disks.

\section{Introduction: The Star-Formation Process}

Much of the interstellar gas is hierarchically structured, with
relative velocities that increase approximately as the square root of
the separation.  These and other power law correlations between cloud
properties imply that the structure is scale-free, in which case the
title of this paper, "small-scale mixing" can actually refer to a wide
range of scales.  We restrict the discussion, though, to contamination
and mixing processes that operate on scales smaller than or
equal to one scale height.  This includes much of the gas dynamics
prior to star formation, essentially all of it during star formation,
and some of it after star formation.  We also consider only gas
processes, and not stellar migration or other stellar mixing
processes, even though these stellar processes may dominate the
observed star-to-star abundance dispersion at any particular galactic
radius.  A third restriction is to galaxies like the Milky Way, and in
some discussions, to the Solar neighborhood of the Milky Way, because
this is the type of region where we know the most about abundance
dispersions.

\section{Mixing of the ISM Before Star Formation Begins}

\begin{figure}[t]
\plotfiddle{elmegreen_fig1.ps}{70truemm}{0}{40}{40}{-130}{-50}
\caption{An image of M51 made from
the ratio of 15$\mu$ emission to B band emission, showing both
the warm emitting dust and the cold absorbing dust as bright features.
The gas concentrations that lead to star formation in the spiral arms
are easily seen at this resolution of 15'.  Image from Block et al.
(1997).  }
\label{fig:m51}
\end{figure}

Most star formation seems to begin on a large scale, comparable to the
thickness of the Galaxy disk.  In galaxies like ours, which contain
strong spiral density waves, it begins in the spiral arms in giant gas
clouds that typically have a mass of $\sim10^7$ M$_\odot$.  These are
observed in atomic gas in many galaxies, such as the Milky Way
(Grabelsky et al.  1987; Elmegreen \& Elmegreen 1987), M100 (Knapen et
al.  1993), and M31 (Lada et al.  1988), and they are observed in
molecular gas in other galaxies, such as M51 (Rand \& Kulkarni 1990).

Figure \ref{fig:m51} shows the gas distribution for M51 made from a
combination of dust absorption in the B band and dust emission at 15
$\mu$ from ISO (Block et al.  1997).  The large knots in the spiral
arms, separated by several kiloparsecs, are the main regions of star
formation, as shown by H$\alpha$ and CO.  Several giant shells can
also be seen in the northern and southern arms in the region of
corotation.

Spiral arm gas concentrations like those shown in figure \ref{fig:m51}
presumably form by the gravitational collapse of
density-wave-compressed gas, involving collapse motions along the
direction of the arm.  The existing theory of this process gives the
correct cloud spacing of $\sim3 $ spiral arm widths, or about
2.5 kpc in galaxies like ours, and it gives the correct
cloud mass and a short enough time scale for these clouds to form
while the gas is still in the arm (Elmegreen 1994).  The formation of
such giant clouds is not limited to density wave arms, however, since
equally large clouds appear in galaxies with only weak arms
(Thornley \& Mundy 1997a,b)
and in irregular galaxies (Cohen et al. 1988; Hunter 1997).  The
peculiar conditions in the crest of a strong spiral arm, namely, the
high gas density, the low rate of shear, and the low galactic tidal
force, make this place highly favored over other sites in galaxies
that have spiral waves.

The onset of star formation in spiral arms implies that star-forming
clouds begin their lives in violent conditions, possibly inside shock
fronts, over a galactic radial range of $\sim1$ kpc, considering a
typical arm pitch angle of $\sim15^\circ$.  This means that the
pre-star-forming gas will be highly turbulent, and it will have a
range of metallicities from the background galactic gradient that may
be $\pm0.05$ dex.

The onset of star formation inside such a cloud is relatively quick.
It seems to being after a time comparable to the turbulent crossing
time, regardless of scale.  This time is is $\sim0.7S(pc)^{0.5}$ My
for scale size $S$, as determined from the molecular cloud
size-linewidth correlation, shown in figure \ref{fig:sl}.  At the top
of the figure is the linewidth versus size, and at the bottom is the
ratio of size to linewidth, versus size.  The {\it duration} of star
formation in a region is also proportional to the local crossing time,
perhaps equal to several times this for a wide range of scales
(Elmegreen \& Efremov 1996).  That is, small regions or small clumps
form their stars quickly and disperse, while large regions or large
clouds take a long time.  A schematic diagram of this trend is shown
in Figure \ref{fig:hierarchies}.  The correlation between crossing
time and size, shown at the bottom of figure \ref{fig:sl}, is
similar to the correlation between duration of star formation and
size, shown in figure \ref{fig:hierarchies}.

\begin{figure}[t]
\plotfiddle{elmegreen_fig2.ps}{80truemm}{0}{100}{100}{-300}{-260}
\caption{Size and linewidth data from GMC surveys. Symbols are for GMC
surveys = {\it dots}: Solomon et al. (1987), {\it open circles}: Dame et al. (1986);
quiescent clouds = {\it filled triangles}: Falgarone et al. (1992), {\it open
triangles}: Williams et al. (1994; the Maddalena-Thaddeus cloud),
{\it inverted open triangles}: Lemme et al. (1995; L1498),
{\it inverted filled triangles}: Loren (1989; Ophiuchus),
OB associations = {\it crosses}: Williams et al. (1994; Rosette),
{\it plus signs}: Stutzki \& G\"usten (1990; M17).}
\label{fig:sl}
\end{figure}

A quick onset to star formation in a cloud born under violent and
turbulent conditions implies that any non-uniformity of abundances
that may enter the cloud, from an initial non-uniformity in the
pre-cloud gas, for example, or from a galactic gradient, will not
smooth out much before star formation actually begins.	That is, {\it
there is not much time for mixing before star formation.} Indeed, the
dominant mode of star formation in galaxy disks involves giant cloud
complexes that are likely to be initially inhomogeneous in elemental
distribution at a level of perhaps $\pm0.05$ dex.  The stars that form
in these clouds should be equally inhomogeneous.

\begin{figure}[t]
\plotfiddle{elmegreen_fig3.ps}{70truemm}{0}{100}{100}{-280}{-310}
\caption{Schematic diagram showing how the duration of star formation
in regions of various sizes increases with the square root of the
size.}
\label{fig:hierarchies}
\end{figure}

\section{Mixing of Supernova Debris During Star Formation}

Hierarchical structure in star-forming regions implies that star
formation is concentrated into a small fraction of the volume, and the
ensuing supernova contamination is localized.  For fractal clouds, the
volume filling factor of dense gas is \begin{equation}f_{\rm
dense\;\;gas}={\rm peak\;density\;contrast}^{(D/3)-1} \sim0.1
\end{equation} for fractal dimension $D=2.3$ (Elmegreen \& Falgarone
1996) and peak density contrast of $\sim10^4$.	Thus, in any region,
star formation occurs in only 10\% of the volume.

The typical size of a supernova remnant in a region of star formation
is only $\sim10-50$ pc.  From Cioffi, McKee \& Bertschinger (1988) the
final size of a remnant following the pressure-driven snowplow phase
is
\begin{equation}
R_{final}=4.93R_{PDS}\left({{E_{51}^{1/14}n_0^{1/7}}\over{a_{10}}}\right)^{3/7}
\end{equation}
where
\begin{equation}
R_{PDS}(pc)=14.0\left({{E_{51}^{2/7}}\over{n_0^{3/7}}}\right).
\end{equation}
Here $E_{51}$ is the supernova energy, $n_0$ is the preshock density,
and $a_{10}$ is the preshock velocity dispersion in units of $10$ km
s$^{-1}$.  Thus
\begin{equation}
R_{final}(pc)=69E_{51}^{0.32}n_0^{-0.37}a_{10}^{-0.43}.
\end{equation}
But $a({km\;s^{-1}})\sim0.7S(pc)^{1/2}$ from the molecular cloud
correlations, and $n_0\sim5000S(pc)^{-1}$ for pressure $\rho
a^2\sim4\times10^5$ k$_B$ in a star-forming region.  Then we get
\begin{equation}
R_{final}(pc)\sim10E_{51}^{0.32}S(pc)^{0.15},
\end{equation}
which is relatively insensitive to scale.  If the supernova occurs in
a region that is cleared of gas by previous HII regions and wind
pressures, then it can become larger than this.

The size of a region of star formation in which OB stars form is
typically larger than the size of a supernova remnant in such a
region.  There {\it is} a characteristic (or minimum) size for such OB
star-forming regions even though there is no characteristic size for
star-forming regions in general.  This is because it takes a minimum
number of $\sim10^4$ stars before an O-type star is likely to form at
all under a random sampling of the initial stellar mass function.
Thus low mass clouds tend to form only low mass stars (because these
stars are common) and high mass clouds tend to form OB stars in
addition to low mass stars.  The star-forming regions in which
supernovae are likely to occur span a distance scale greater than
$\sim10$ pc, possibly even greater than $\sim 100$ pc, as suggested
also by the plotted points for OB subgroups and associations in figure
\ref{fig:sl}.  These scales are comparable to or larger than the
associated supernova remnant sizes, so {\it supernova contamination in
a region of star formation should be spotty and localized}.  The
supernova debris will therefore mix with only part of the star-forming
gas, and may even get highly concentrated in several small regions if
the explosion happens to occur at close range to the cloud.  Note that
$R_{final}$ in the above equation is small, so according to this, a
lot of contamination will be at a range less than $\sim10$ pc.

Debris from multiple supernovae in an OB association should
contaminate a few dense clouds and leave the rest in a shadow or far
away.  Thus {\it the clumpy, hierarchical nature of star-forming
regions should lead to abundance inhomogeneities during the star
formation/supernova epoch}, which is a $3-20$ My year period following
the onset of star formation in a region where OB stars form.  Studies
of supernova remnants in regions of star formation were made by Junkes
et al.	(1992) and others.

\section{Penetration of Supernova Debris into a Shocked Cloud}

In a simple bow shock, there is a contact discontinuity between the
dense cloud and the SN gas, and this surface has no direct mixing in
equilibrium conditions.  Kelvin-Helmholtz and Rayleigh-Taylor
instabilities on this surface will mix the gas, however, and these
instabilities should be responsible for contamination of the cloud
by supernova material.	In numerical simulations of this process, such
mixing is invariable {\it behind the cloud}, and not in the cloud core
(e.g., Xu \& Stone 1995).

It is unknown if SN debris also mixes with the gas in a compressed
cloud core.  Stone \& Norman (1992) simulate a shock/cloud collision
with $\gamma=5/3$ and suggest that a Richtmyer-Meshkov instability
might lead to mixing at the front surface and inside the core.	Foster
\& Boss (1996, 1997) consider a shock/cloud collision with $\gamma=1$
and get a very thin region between the cloud and the shock in the SN
ejecta.  This thin region allows Rayleigh-Taylor fingers at the head
of the cloud to penetrate into the core and aid mixing with dense,
pre-star-formation material.

Supernova interactions with fractal or hierarchically structured clouds
have not yet been modeled.

\section{Mixing of SN-contaminated Clouds after SF}

The size-linewidth-time relation in a hierarchical cloud implies that
many small regions of SF come and go before the large region is
finished.  This is also evident from the schematic diagram in figure
\ref{fig:hierarchies}.	One implication of this is that the large
star-forming clouds, which cannot be easily destroyed by internal star
formation, probably move around and get recycled as even larger scales
in the hierarchy continue to form stars.  In this case, GMC's can have
several epochs of supernova contamination, and each one can increase
the dispersion of the final stellar abundances in a stochastic fashion
before all star formation stops in the region.

An example is offered by the formation in the Solar neighborhood.
Most of it began 30-50 My ago when the neighborhood was shocked by the
local spiral arm, which is now the Carina arm.	This assumes a pattern
speed for the Galactic spiral system of 13.5 km s$^{-1}$ kpc$^{-1}$
from Yuan (1969).  The first generation of star formation
probably made the
Cas-Tau association (Blaauw 1984), in addition to some
clusters, such as $\zeta$ Sculptoris and
the Pleiades.
This first generation also could have caused the expansion of
Lindblad's ring, which now seems to contain the Orion, Perseus, and
Sco-Cen OB associations along the periphery (see review in P\"oppel
1997).	The Cas-Tau association at the center of this ring now has
very little gas, so presumably its residual gas and other spiral arm
gas got collected into
Lindblad's ring, and eventually into the local associations.  Thus the
current supernova contamination of Orion and other local clouds is
probably the second time this has happened since the spiral wave
passed by.

The tilt of Gould's Belt has been attributed to the impact of an
extragalactic cloud (Franco et al. 1988; Comer\'on \& Torra 1993),
but there are other explanations
such as a Parker instability (Gomez de Castro \&
Pudritz 1992; Shibata \& Matsumoto 1991), and
general perturbations from the LMC (Edelsohn \& Elmegreen 1997).
Maybe the spiral arm itself is responsible for the tilt.  The shock in
a spiral arm contains more than enough speed to elevate the
gas to one scale height; any asymmetry in the spiral potential
relative to the gas midplane could conceivable lead to a tilted gas
distribution when giant clouds form.

The primary disadvantages of the cloud impact model are (1) that the
local star formation seemed to begin just when the local region was at
or near the crest of a spiral arm, and such timing for an impact would
be unlikely, and (2) that corrugations like this are common in the
spiral arms of the inner galaxy (Pandey, Bhatt, \& Mhara 1988; Sanders,
Solomon,
\& Scoville 1984) and other galaxies (Florido et al. 1991),
and it is unlikely all were
formed by cloud impacts.  Nevertheless, the impact model
is intriguing from the point of view of abundance anomalies
in Orion and other local regions, because an impacting high velocity
cloud could have contaminated the disk with low metallicity material
(Cunha \& Lambert 1992; Edvardsson et al. 1995).

\section{Mixing of SN-contaminated Cloud Debris after SF}

Ultimately the star-forming gas gets converted to low density by
ionization, heating, and evaporation.  These processes cause {\it
atomic level} mixing.

The large-scale star-forming regions are also stretched by increased
shear and tidal forces when they enter the interarm.  These processes
cause {\it cloud-to-cloud} mixing rather than atom-to-atom mixing, but
they also expose the dense gas to stray stellar enegy which can lead
to more ionization, heating, and evaporation.  It follows that {\it
passing from an arm to the interarm region should help homogenize
the gas at both the cloud and the atomic level.}

In the inner part of the Milky Way, the next spiral density wave arm
arrives rather quickly, in only 50 My years or less.  This is before
cloud destruction is complete so there is a significant amount of
dense molecular gas and star formation still in the interarm.
For example, figure \ref{fig:m51} shows a lot more gas and dense
knotty emission in the interarm region at small radii than large
radii.

Based on such analyses, one might expect cloud self-destruction to be
more complete in the interarm regions in the outer part of the Milky
Way than the inner part.  This implies that the interarm
homogenization of supernova enriched gas should change with radius,
and so should the relative dispersion in abundances.

Post-star-formation mixing is also possible in chimney vents and giant
bubbles that send enriched material into the halo where it can
disperse over a radial range of perhaps a kiloparsec or more.  In the
hot gas, the supernova debris can mix atom-by-atom, so this form of
dispersal is at the atomic level.

\begin{figure}[t]
\plotfiddle{elmegreen_fig4.ps}{60truemm}{0}{100}{100}{-160}{-130}
\caption{Schematic diagram of the change with spiral arm phase of the
elemental abundance and the dispersion in the elemental abundance,
based on the qualitative model presented in the text.
Features related to star formation (SF), spiral arms and interarms,
spiral density wave (SDW) shocks, and supernovae (SN) are shown.
}
\label{fig:zscheme} \end{figure}

\section{Conclusions}

1. Pre-star formation clouds should be inhomogeneous as a result of
their chaotic and large-scale formation processes.\\
2.  Inhomogeneities should increase during star formation as a result
of the clumpy, hierarchical structure of young stars.  The penetration
depth of SN debris into triggered clouds is unknown, however.\\
3.  Short-time cloud recycling should continue to increase the
metallicity and the metallicity dispersion until all star formation
stops.\\
4.  Ionization, heating, evaporation, and hot chimneys should
homogenize the gas at the atomic level, while general galactic shear
should mix at the cloud-to-cloud level.\\
5.  Passing to an interarm region should disperse a cloud complex and
aid the homogenization process by allowing deeper penetration of
background radiation.

A schematic diagram of the progressive change in abundance and
abundance dispersion with phase in a spiral density wave is shown in
figure \ref{fig:zscheme}.  The contaminations from Type I and Type II
supernova debris are tracked separately, because only the Type II
debris (e.g., Oxygen) correlates with star formation.  The Type I
debris (e.g.  Iron) is added more continuously to the gas as the comoving
old stars evolve into supernovae.  The
vertical scales on the left and right are pure conjecture, but they
sensibly follow from the discussion given above.

Small-Scale Mixing of the Elements


Raman spectroscopy on etched graphene nanoribbons

Abstract: We investigate etched single-layer graphene nanoribbons with different widths ranging from 30 to 130 nm by confocal Raman spectroscopy. We show that the D-line intensity only depends on the edge-region of the nanoribbon and that consequently the fabrication process does not introduce bulk defects. In contrast, the G- and the 2D-lines scale linearly with the irradiated area and therefore with the width of the ribbons. We further give indications that the D- to G-line ratio can be used to gain information about the crystallographic orientation of the underlying graphene. Finally, we perform polarization angle dependent measurements to analyze the nanoribbon edge-regions.

\section*{Introduction}

Graphene nanoribbons~\cite{han07,lin08,wan08,mol09b,han09,mol10,gal10} attract increasing attention due to the possibility of building graphene-based nanoelectronics as for example field-effect transistors~\cite{wan08,zha08} or quantum dot devices~\cite{sta08a,pon08}. In contrast to two-dimensional gapless bulk graphene~\cite{gei07}, it has been shown that confinement~\cite{fer07a,yan07a}, disorder~\cite{sol07} and edge effects~\cite{yan07a} introduce a transport gap in graphene nanoribbons.
The fabrication technique may influence the transport properties of the nanoribbons in terms of added bulk and/or edge disorder. Disorder is expected to strongly influence the scaling behavior of the energy gap as a function of the nanoribbon width and the local doping profile~\cite{mol09b,han09,mol10,gal10}. In addition, very little is known about the edge structure and theoretical investigations of the vibrational properties of nanoribbons have only been started very recently~\cite{gil09}.
Raman spectroscopy on carbon (nano)materials~\cite{mal09,fer07} has been recognized as a powerful technique not only for probing selected phonons, but also for identifying the number of graphene layers~\cite{fer06,dav07a}, for determining local doping levels~\cite{sta07aa}, for studying electron-phonon coupling~\cite{pis07} and thus for the electronic properties themselves.

In this letter we report Raman spectroscopy experiments on etched graphene nanoribbons with different widths $w$ ranging from 30 to 130~nm (schematic in Fig.~1a). We show that the characteristic signatures of single-layer graphene (SLG) in the Raman spectra are still well preserved, that the absolute G- and 2D-line intensities scale with the nanoribbon width, whereas the D-line intensity does not. Consequently, the D-line intensity depends only on the edge-region of the nanoribbon including the edge  roughness which can be further analyzed by performing polarization dependent measurements.

\section*{Fabrication}

The nanoribbon fabrication is based on the mechanical exfoliation of natural graphite~\cite{nov04}, electron beam lithography, reactive ion etching and metal evaporation. For details see Ref.~\cite{mol09b}. The graphene flakes have been identified to consist of a  single-layer by measuring the Raman full width at half maximum (FWHM) of the 2D-line prior to processing~\cite{fer06,dav07a}. In Fig.~1b we show a scanning force microscope (SFM) image of six etched graphene nanoribbons with different widths. In total, we have studied three such nanoribbon arrays fabricated from three different single-layer flakes resulting in a total of more than 20 individual nanoribbons.

\section*{Experimental Setup}
All Raman spectra were acquired using a green laser (532~nm, $\hbar \omega_L=$ 2.33 eV). Employing a long working distance focusing lens (numerical aperture~=~0.80), we obtain a spot size with a diameter $d~\approx$ 400~nm. The incident light was -- unless stated differently -- linearly polarized parallel to the macroscopic edge of the nanoribbons ($\theta=0^\circ$) and detection was always insensitive to polarization. The laser power was set to $\approx$ 2~mW in order to exclude heating effects and all measurements were conducted at room temperature.

\begin{figure}[bt]
	\centering
	\includegraphics[width=0.5\linewidth]{Figure1
	\caption[FIG1]
(color online)
(a) Schematic illustration of etched nanoribbons with widths $w$. The yellow circle represents the Raman laser spot of the linearly polarized light with angle $\theta$.
(b) SFM image of six etched graphene nanoribbons (vertical white bars) with corresponding gold contacts (yellow) which were used to align and identify the nanoribbons during the Raman measurements. The nanoribbon width $w$ decreases from left to right from 110~nm down to 80~nm and was determined by the SFM cross-section shown in the lower part of the image (white line). The laser spot-size is represented by the white circle.
(c) Raman line scan along path $C$ [see panel (b)]. For each data point, a Raman spectrum was recorded and the intensity of the D-line was calculated by summing up the detector counts (from 1306-1375 cm$^{-1}$).
(d) Same as above, but for the G-line (1558-1609 cm$^{-1}$) and
(e) the 2D-line (2618-2734 cm$^{-1}$).
(f) Maximum intensity of the G- and the 2D-line plotted versus the nanoribbon width
}
	\label{fig:LineScan}
\end{figure}

\begin{figure}[t]
	\centering
	\includegraphics[width=0.5\linewidth]{Figure2}
	\caption[FIG2]
(color online)
A selection of Raman spectra from nanoribbons of different widths. For each spectrum, the silicon background signal was removed.
In the inset of panel (e), magnified versions of G- and 2D-lines are shown. The blue dots represent measured data points whereas the black lines represent best fits.
	}
	\label{fig:SpectraOverview}
\end{figure}

\section*{Results}

A selection of Raman spectra of nanoribbons with different widths $w$ is shown in Fig.~2. In all measurements, the characteristic signatures of graphene Raman spectra can be well identified: The defect induced D-line ($\approx$1340~cm$^{-1}$), the G-line ($\approx$1580-1590~cm$^{-1}$), the 2D-line ($\approx$2680~cm$^{-1}$) and for the thinner ribbons an additional defect induced D'-line ($\approx$1620~cm$^{-1}$). For a detailed review on the microscopic nature of all these lines, see Refs.~\cite{tho00,rei04,fer06,dav07a,fer07,pis07,mal09}.

The FWHM of the 2D-line with values below 40~cm$^{-1}$ (see right labels in Figs.~2a-d) is a strong signature for the single-layer nature of the graphene nanoribbons~\cite{fer06,dav07a}. In order to enhance the signal for thin nanoribbons, we recorded the spectra of the 30~nm nanoribbon on a sub-array of several closely spaced nanoribbons with equal width (approx. 6 ribbons irradiated by the laser). This sample averaging may also explain the rather large but still single-layer FWHM of the 2D-line of nearly 50~cm$^{-1}$. In contrast to defectless bulk graphene where significant D- and D'-lines are rarely present, we observe strong D- and D'-lines in the recorded spectra as the edges act as defects and allow elastic inter-valley scattering of electrons~\cite{can04,can04a,cas09}.

As demonstrated in Fig.~1c, the D-line intensity does not depend on the width of the nanoribbons. The peak width in Figs.~1c-e arises from strong spatial oversampling as the laser spot size ($\approx$ 400~nm) is significantly larger than the step size ($\approx$ 70~nm). As the laser spot size is also significantly larger than the nanoribbon width ($d\gg w$), the total irradiated edge length is approximately $2\times d$ (left and right edge) -- independent of the nanoribbon width. As long as there is no bulk disorder present, the intensity of the D-line is not expected to depend on the nanoribbon width. This is observed in Fig.~1c and leads to the important conclusion that the reactive ion etching process used to pattern nanoribbons does not introduce a detectable amount of bulk defects into our graphene nanostructures.

In contrast to the D-line, we observe a width dependence of the G- and 2D-line intensities as shown in Figs.~1d,e. The intensity of the G-line is a function of the amount of irradiated sp$^2$ bound carbon atoms~\cite{mal09} and therefore expected to directly depend on the graphene nanoribbon area below the laser spot ($\approx d\times w$). The same proportionality also holds for the 2D-line~\cite{tui70,cas01,Laz08}. This expected linear dependence of the G- and 2D-line heights is confirmed in Fig.~1f (see lines fitted to circles and triangles). It is crucial to note that Fig.~1f does not allow to determine the width of a possibly existing Raman-inactive edge-region~\cite{han07}. This is due to offsets in Figs.~1c-e (background of the Si substrate) and due to the inherent noise in the measurement, which masks the small line heights resulting from the short integration time. From the fact that we can measure 30~nm wide nanoribbons (see Fig. 2e), we conclude that such a Raman-inactive edge-region  must either be significantly smaller than 15~nm or scale with the nanoribbon width.

\begin{figure}[t]
	\centerin
	\includegraphics[width=0.5\linewidth]{Figure3}
	\caption[FIG3]{
(color online)
(a)
Intensity ratio $I(D)/I(G)$ versus $w^{-1}$ for different nanoribbons located on the three different arrays ($\#1-\#3$).
(b)
Polarization angle dependence of I(D)/I(G) for nanoribbons of different width. The inner gray circle is located at I(D)/I(G)=0 and the outer at I(D)/I(G) = 3.25. The data points are fitted with $I(D)/I(G)(\theta)=a[b+(1-b)\cos^2\theta]$.
(c)
Ribbon width dependence of the fit parameters $a$ and $b$.
	}
	\label{fig:Orientation}
\end{figure}

In order to compare the three arrays ($\#1-\#3$), we show in Fig.~3a the intensity ratio $I(D)/I(G)$ versus the inverse ribbon width. Each array is made from a single SLG flake and all nanoribbons on a specific array were etched in the same processing step and are oriented parallel to each other. As a measure for the intensity, the peak area was used~\cite{com0f}. Based on the previous argument, the relation $I(D)/I(G) \propto w^{-1}$  is expected and observed. It is important to note that this relation only holds in the limit $d\gg w$ (otherwise there is no edge located under the laser spot). Therefore, a linear fit was performed for each array and the intersection of the fit with the abscissa  was pinned to $w^{-1}=1/d$.

Interestingly, the fits for the different arrays exhibit different slopes. There are two possible explanations: (i) The edge roughness differs for the different arrays: Rougher edges result in more defects illuminated by the laser spot and therefore a higher D-line intensity. As all the arrays were fabricated using the same process, this is improbable. (ii) The slopes can also be attributed to different crystallographic orientations of the graphene relative to the edge. This is because perfect zig-zag edges do not activate the D-line due to momentum conservation~\cite{can04,can04a}, whereas armchair segments do. While the plasma etched edges are certainly rough, the ratio of zig-zag to armchair segments should nevertheless depend on the overall crystallographic orientation relative to the nanoribbon orientation.

In Fig.~3b, we show the $I(D)/I(G)$ ratio for different nanoribbons of one array as a function of the incident photon polarization angle  $\theta$. For each ribbon and each angle, a spectrum was recorded and both the G- and the D-line were fitted each with a single Lorentzian. Geometric considerations suggesting mirror planes at 0$^{\circ}$ and 90$^{\circ}$ are confirmed. For Raman measurements in graphene, only phonon wave vectors perpendicular to the incident light polarization contribute to the signal~\cite{gru03,can04a}. As phonon wave vectors along the principal nanoribbon axis (i.e. $\theta = $90$^{\circ}$ and 270$^{\circ}$) see little of the edges, the D-line is expected to be more suppressed than at all other angles.

Following Casiraghi et al.~\cite{cas09} describing similar measurements on graphene edges and taking into account that $I(G)$ does not depend on $\theta$~\cite{com0c}, the following expression can be derived: $I(D)/I(G)(\theta)=a[b+(1-b)\cos^2\theta]$, where $a=I(D)_{max}/I(G) \propto 1/w$ and $b=I(D)_{min}/I(D)_{max}$. This relation was used to fit the data and indeed a strong width dependence of $a\propto 1/w$ and an overall constant of $b \approx$ 0.55 only depending on the structure of the edges are found (Fig.~3.c). According to Ref.~\cite{cas09}, a lower bound of the edge disorder correlation length can be estimated by $\xi=2 v_F/(\omega_L b)\approx$ 1~nm, which is in reasonable agreement with current fabrication limitations~\cite{com0a}.

\section*{Conclusion}
In conclusion, we have demonstrated that the D- and D'-line depend only on the nanoribbon edge-region  whereas the G- and the 2D-line scale with the illuminated area. We have shown that our fabrication process does not introduce bulk defects and that the $I(D)/I(G)$ ratio can give indications about the crystallographic orientation of graphene. These insights may help in designing further experiments and designing future graphene nanoelectronic devices.

\section*{Acknowledgment}
The authors thank C. Casiraghi, A. C. Ferrari, M. Haluska, A. Jorio and L. Wirtz for helpful discussions and C. Hierold for providing access to the Raman spectrometer. Support by SNSF and NCCR nanoscience is gratefully acknowledged.

Mesoscale and Nanoscale Physics (cond-mat.mes-hall)