Toward an Understanding of the Baldwin Effect

Abstract: After a short review of the quasar problem and what we can hope to learn from their emission line spectra, I discuss the current body of knowledge concerning quasar emission lines and their relationships to the local and ionizing continua. I propose a hypothesis that the Baldwin effect is due to a characteristic relationship between the continuum spectral energy distribution, the gas metallicity (Z), and the quasar luminosity. I suggest that such a relationship might arise naturally from a scenario involving massive galaxy formation and evolution driving the birth and evolution of the quasar central engines in terms of four fundamental parameters: M_bh, Mdot Mdot/M_bh, and Z.

\section{Introduction}

\subsection{Some history}

Quasars have long been mysterious points of light in our sky.
Originally turning up as star-like radio sources in the early radio
surveys of the late 1950s, their spectra remained a puzzle until
Maarten Schmidt's realization (Schmidt 1963) that their emission lines
are as those seen in emission line galaxies (Seyfert 1943), but greatly
redshifted. The ``discovery object'' was 3C~273 with a redshift of
0.158. This placed quasars billions of light years away and made them
the most distant and luminous entities in the universe. As more were
discovered, the redshift limit was pushed ever higher, meaning that we
were probing ever-earlier epochs in our Universe's history. Soon it was
realized that if we could establish quasars as``standard luminosity
candles'', we could determine their distances independent of their
redshifts and thus determine the expansion history of the Universe.
Thirty-five years later (my lifetime) we are still puzzling over the
story that quasars are telling us.

\subsection{The quasar paradigm}

While a comprehensive understanding of the quasar phenomenon remains
illusive, a standard model has emerged: a supermassive black hole in
the nucleus of a massive galaxy accretes matter at rate $\dot{M}$,
converting gravitational potential energy into light ($L = \eta \dot{M}
c^2$; $\eta \sim 0.1$) and mechanical energy (Blandford \& Rees 1978;
Rees 1984).  Evidence continues to trickle in that quasars represented
a dynamically active stage in the lives of massive galaxies (McLeod \&
Riecke 1995a, 1995b; Bahcall et al.\ 1997; McLure et al.\ 1998), with
the masses of their central black holes related to the host galaxy's
halo mass (Haehnelt, Natarajan, \& Rees 1998; Laor 1998, and these
proceedings).  The rapid appearance of quasars with redshifts $z > 2$
was likely associated with massive galaxy formation (Boyle \& Terlevich
1998; Haehnelt et al.\ 1998; Miller \& Percival 1998; Silk \& Rees
1998). Quasars are now extinct, though their low luminosity cousins,
the Seyfert galaxies, are still found in the here and now.  During
their reign, the quasar ``central engines'' poured out tremendous
numbers of UV/X-ray photons.  In the standard model the mass of the
central black hole and the mass accretion rate together determine the
quasar's luminosity and continuum spectral energy distribution (SED).

\subsection{Quasar emission lines and the stories they tell}

\subsubsection{``Sometimes you can't stick your head in the engine, so
you have to examine the exhaust.'' (Don Osterbrock c.\ 1987).}

The central engine of a quasar is manifested in the accretion of matter
onto a supermassive black hole. This process releases energy, much of
this in the form of light at ionizing energies (e.g., Mathews \&
Ferland 1987) which is generally unobservable due the opacity of the
intergalactic and interstellar media. Therefore, we cannot discern
directly how a large fraction of the energy in the central engine is
dissipated.  However, some of these continuum photons are intercepted
by gas of unknown origin moving rapidly in the black hole's
gravitational potential. These photons ionize and heat this non-LTE
gas, which responds by producing the quasar broad emission line
spectrum.  The recombination lines of Ly$\alpha$ $\lambda$1216 and
He~II $\lambda$1640 are a measure of the numbers of photons lying just
shortward of 912~\AA\/ and 228~\AA\/, respectively.  The ratio of the
intensities of the strongest line coolant, C~IV $\lambda$1549, to the
strongest recombination line, Ly$\alpha$, measures the cooling per
recombination, or the heating per photoionization, and thus reflects
the energy balance of the gas. Since the energy from the continuum
source absorbed by the gas must equal the energy emitted by the gas,
the emission lines serve as an indirect means to measure the ionizing
continuum. This and the likelihood that the emission line widths are
indicative of the central gravitational potential well mean that we can
hope to learn from the emission lines the nature of the central engines
of quasars.

\subsubsection{Origin \& Nature of the Emitting Gas.}

The emission line spectra of quasars must be leaving clues, albeit
apparently subtle ones, as to the origin and nature of the emitting
gas. For a discussion of this see my other contribution to these
conference proceedings.

\subsubsection{Chemical abundances.}

The strengths of the heavy element emission lines in even the highest
redshift quasars give testimony to the heavy element enrichment from
stellar evolution that had to occur within the first billion years of
cosmic history. If they can be understood, quasar emission lines will
elucidate the chemical evolution of galactic nuclei over the time span
of $\sim$~10\% -- 90\% the age of the Universe (Hamann \& Ferland 1993;
Hamann \& Ferland 1999; Hamann, these proceedings).

\subsubsection{Quasar evolution.}

Quasars are known to have evolved in number density and luminosity over
past 10 or so billion years. Their rapid rise at high redshifts
(presumably due to galaxy formation), their reign near $z \approx 2.5$,
and their demise at lower redshifts, as the source of gas ``ran dry''
and/or the accretion mechanism itself failed (Haiman \& Menou 1998), is
reasonably well established (Boyle, Shanks, \& Peterson 1988; Boyle
1991; Hewett, Foltz, \& Chaffee 1993; Goldschmidt \& Miller 1997; La
Franca \& Cristiani 1997; Graham, Clowes, \& Campusano 1998). Quasar
evolution may be tied to the dynamical and chemical evolution of
massive galaxies ($\S~1.2$). But we do not understand these things at
the fundamental level. Quasar emission lines, mirroring the action
occurring in galactic nuclei, are keys to unlocking these mysteries.

\subsubsection{Beacons of cosmic history?}
Since the 1970s we've known of the so-called ``Ly$\alpha$'' forest of
absorption lines in the spectra of quasars representing matter along
the line of sight to the quasar. The statistics in number and character
of the Ly$\alpha$ forest demonstrate an evolution with look-back time
in the gas that may represent galaxies, protogalaxies, or failed
galaxies. However, quasars themselves are strewn throughout most of
cosmic history and their intrinsic light must be telling us something
about that history. Their connection to galaxy evolution is an
example.  Another is if their light could be somehow established as a
standard luminosity candle, we might establish the expansion history
and fate of the Universe. Detailed analysis of quasar spectra became
possible in the late 1970s with the advent of sensitive electronic
spectrographs. It was at this time that a young astronomer by the name
of Jack Baldwin was compiling spectral characteristics of quasars, when
he came across something interesting.

\subsection{Observed line, continuum, \& luminosity correlations}

This is a brief summary of the situation. For a more complete review
the reader is directed to the conference proceedings contribution from
Osmer \& Shields.

\subsubsection{The Baldwin Effect.}

Uncharacteristic of Galactic emission line sources (e.g., novae, PNe,
H~II regions), the spectra of quasars are remarkably homogeneous. The
same emission lines sitting atop a ``blue bump'' continuum, with
similar intensity ratios, are observed again and again --- this, in the
face of 4--5 orders of magnitude in their luminosities and billions of
years of lookback time. However, upon comparison of the C~IV emission
line flux to that of the underlying continuum flux density
$f_{\lambda}$, Jack discovered that this ratio is in inverse proportion
to the monochromatic ultraviolet continuum luminosity $L_{uv}$
(Baldwin 1977). \begin{equation} \log W_{\lambda}(C~IV) \propto
\beta\log L_{uv}, ~\beta < 0 \end{equation} This effect is illustrated
dramatically, though qualitatively, in Figure~1. If calibrated, a
measurement of a quasar's C~IV equivalent width would indicate its
luminosity and thus its luminosity distance.

\begin{figure}
\plotfiddle{korista1fig1.eps}{3.0in}{270.}{70}{70}{-250}{400}
\caption{A qualitative illustration of the Baldwin effect. The solid
curve is the continuum normalized mean UV spectrum of the Seyfert~1
nucleus of NGC~5548 ($z = 0.0174$) from the 1993 monitoring campaign
(Korista et al.\ 1995). The dotted curve is a portion of the revised
version of the composite quasar spectrum from the Large Bright Quasar
Survey (Morris 1994), also continuum normalized. Effectively $\sim$~3
orders of magnitude in UV continuum luminosity separate these two
spectra. Note, too, that the narrow emission line spectrum is
completely missing from the quasar composite. The blue wing of the Sy~1
Ly$\alpha$ profile is contaminated by geocoronal Ly$\alpha$, and in the
quasar composite spectrum the region shortward of the peak of
Ly$\alpha$ is affected by ``Ly$\alpha$ forest'' absorption. This is a
log-log plot.}
\end{figure}

Since Baldwin's initial discovery, several studies have found similar
relationships for UV emission lines in addition to C~IV, and that each
emission line has its own logarithmic slope in the $W_{\lambda} -
L_{uv}$ relation, perhaps offering clues to its origin (Baldwin,
Wampler, \& Gaskell 1989; Kinney, Rivolo, \& Koratkar 1990; Zamorani et
al.\ 1992; Espey, Lanzetta, \& Turnshek 1993; Osmer, Porter, \& Green
1994; Laor et al.\ 1995; Zheng, Kriss, \& Davidson 1995; Turnshek
1997). Kinney et al.\ (1990) also showed that the Baldwin effect
extends down to Seyfert nuclei luminosities.  Unfortunately, many of
these same studies also demonstrated the existence of substantial
intrinsic scatter in the relations. Thus, the use of quasars as
``beacons of cosmic history'' via the Baldwin effect would have to
await some understanding of quasar spectra.

\subsubsection{Other correlations.}

A number of other important quasar spectral correlations have been
amassed since Baldwin's discovery in 1977 and these may provide clues
to the latter's origin.

Quasars are intrinsically luminous X-ray sources. A number of quasar
X-ray emission studies have found a relationship between the UV/X-ray
continuum flux ratio \begin{equation} \frac{f_{uv}}{f_{x}} =
\left(\frac{\nu_{uv}}{\nu_{x}}\right)^{-\alpha_{uvx}} \end{equation}
and the UV continuum luminosity $L_{uv}$, in that higher luminosity
quasars have characteristically steeper effective logarithmic UV to
X-ray spectral slopes $\alpha_{ox}$ (2500~\AA\/ -- 2 keV) and
$\alpha_{uvx}$ (1350~\AA\/ -- 1 keV) (Zamorani et al.\ 1981; Wilkes et
al.\ 1994; Wang, Lu, \& Zhau 1998). The Wilkes et al.\ (1994) relation
is \begin{equation} \alpha_{ox} \approx 1.53 +
0.11\log\left(\frac{L_{\nu}(2500)}{10^{30}}\right). \end{equation} This
relation extends down to the luminosities of Seyfert nuclei. It should
be noted that the significance of this relationship is a subject of
debate (LaFranca, Franceschini, \& Vio 1995; Avni, Worrall, \& Morgan
1995; Yuan, Siebert, \& Brinkmann 1998).

Line ratios sensitive to the ionizing continuum SED, e.g.,
Ly$\alpha$/C~IV and Ly$\alpha$/O~VI, and the equivalent widths of these
three emission lines all correlate with $\alpha_{ox}$ and
$\alpha_{uvx}$. This is an expected result based upon energetics
arguments {\em if} $\alpha_{ox}$ and $\alpha_{uvx}$ are proxies for the
actual ionizing continuum SEDs of quasars (Netzer, Laor, \& Gondhalekar
1992; Zheng et al.\ 1995; Green 1996; Wang et al.\ 1998).  As explained
above these line ratios reflect the energy balance of the gas and thus
broadly measure the average ionizing photon energy of the incident
continuum. Relationships between $L_{uv}$ and these line ratios and the
line equivalent widths were also found in these studies.

Taken together, these correlations would indicate a general
relationship between the quasar's ionizing SED and its observed
ultraviolet luminosity, and this may play a major role in the origin of
the Baldwin effect.

\section{``Natural Selection'' and the emission line spectra of quasars}

\subsection{Emission line spectrum from a single cloud}

The emission lines of quasars arise mainly or exclusively via the
photoionization of gas, in the so-called broad line region, surrounding
a UV/X-ray ionizing continuum source.  Early successful photoionization
models of emission from single clouds are described in the works of
Davidson \& Netzer (1979) and Kwan \& Krolik (1981). The incident
continuum SED and gas chemical abundances are the basic ``raw
materials'' from which the emitted spectrum from a single cloud
arises.  For a given set of raw materials, the ionizing photon flux,
the gas density, and {\em especially} the ratio of these two parameters
--- proportional to the ionization parameter $U(H)$ --- set the
spectrum.  The gas column density and details concerning the radiation
transfer place the final touches on the emitted spectrum. The early
models found that there was a preferred value of $U(H) \approx 0.01$,
though a separate population of gas with a larger value must also exist
to account for the observed strength of the O~VI $\lambda$1034 line.

\subsection{An observational paradox and its resolution}

\begin{figure}
\plotfiddle{korista1fig2.eps}{6.0in}{0}{70}{70}{-205}{-30}
\caption{Contours of $\log W_{\lambda}$, referenced to the incident
continuum at 1215~\AA\/ and for full source coverage, for six emission
lines (or blends) are shown as a function of the hydrogen density and
flux of hydrogen ionizing photons. The total hydrogen column density is
$10^{23}~\rm{cm^{-2}}$. The $W_{\lambda}$ is in direct proportion to
the continuum reprocessing efficiency. The smallest (generally
outermost decade contour corresponds to 1~\AA\/, each solid line is 1
decade, and dotted lines represent 0.2 decade steps. The contours
generally decrease monotonically from the peak (solid triangle) to the
1~\AA\/ contour. The solid star is a reference point marking the old
``standard BLR'' parameters. ``Sum Ly$\alpha$~1216'' is the sum of
Ly$\alpha$~1216, He~II~1216, and O~V]~1218.}
\end{figure}

The work of many observational spectroscopists of the late 1970s and
early 1980s demonstrated the homogeneous nature of the emission line
spectra of quasars and their lower luminosity cousins, Seyfert nuclei.
Nature must somehow know to adjust the ionization parameter knob to the
nearly the same value through 5 decades in luminosity and 10 billion years
of lookback time.

From the mid 1980s to the mid 1990s, the monitoring of Seyfert nuclei
variability and emission line reverberation (Peterson 1993)
demonstrated that the broad line regions lie nearer to their continuum
sources than thought earlier, and contain gas with a range of particle
densities to accommodate a range in ionization through a region that
spans at least $\sim$~1.5 decades in radius. Rees, Netzer, \& Ferland
(1989) and Goad, O'Brien, \& Gondhalekar (1993) and others explored
models in which the gas density and column density were unique to and
allowed to vary with the distance from the ionizing photon source.
However, the means by which the quasar normalizes its gas properties or
why one power law of gas properties is chosen over another remained
unknown.  Somehow, over a huge range in ionizing continuum luminosity,
gas that is broadly distributed in space, particle density, and
ionization is able to produce the same sets of emission lines in
roughly similar ratios.  What's going on?

In the history of quasar photoionization models, the ``single cloud''
model was followed by the ``two cloud'' model, and then a radial power
law in cloud properties. Recently, Baldwin et al.\ (1995) extended the
``multi-cloud'' idea to its next logical step: ``Locally
Optimally-Emitting Clouds (LOC).'' What if nature provides the broad
line regions with a large pool of gas properties, such that at every
radius there exists a large range in particle and column densities?
Under these conditions the emitted spectrum is controlled by powerful
{\em selection effects} introduced by the atomic physics and basic
radiative transfer effects.  See Figure~2 and Korista, Baldwin, \&
Ferland (1997). Gas having properties that make it most efficient in
reprocessing the continuum into a particular line will emit most of
that line's luminosity. The gas need not fill the entire $n(H) -
\Phi\/(H)$ plane in Figure~2, but unless the gas in the BLR has a
remarkably restricted range of properties, we will in fact observe the
optimally emitting gas for most lines. And just as important, such an
ensemble spectrum reproduces a typical quasar spectrum.

The LOC model suggests that at work is a loose analog of Darwin's
natural selection and it offers natural explanations for: (1) the
observed homogeneity of the quasar spectra and the ``magic ionization
parameter'' problem, and (2) the ionization stratification inferred
from observations. Another important consequence of this idea is that
to a good approximation we can take a ``statistical mechanics''
approach to the modelling of quasar spectra, concerning ourselves with
the ensemble emission, rather than with the details of the properties
and arrangement of the individual gas clouds.  In this picture the
continuum SED and gas chemical abundances are the major drivers of the
emission line spectrum, not a magic ionization parameter.  In this case
we can hope to learn about these two fundamental parameters and the
possible stories they can tell concerning cosmic history and the nature
of quasars.

A word of caution should appear here. The emission line spectrum from
gas distributed in particle density and radius is in general much less
sensitive to changes in the gas abundances and continuum SED than that
emitted by a single ``cloud'' of fixed particle density and ionization
parameter.  This is because that as a cloud's thermodynamics adjusts to
changes in the SED or gas abundances, its reprocessing efficiency can
change immensely. In {\em any} distributed cloud model, there will be
other clouds with other properties at similar distances from the
ionizing source that will also react, and the integrated response is
dampened.  Under the distributed cloud hypothesis, the emitted spectrum
does remain sensitive to the continuum SED and gas abundances, and this
sensitivity will probably better represent nature.

\section{What are the parameters that govern the Baldwin Effect?}

I use the spectral synthesis code, Cloudy (90.04; Ferland et
al.\ 1998), and the emission from an ensemble of plane-parallel clouds
under the assumption of the LOC model to demonstrate that a
characteristic dependence of the continuum SED and gas abundances on
the quasar luminosity may provide an explanation for the Baldwin
effect. For simplicity and because of the spectrum's general
insensitivity to it, I held the column density fixed at $N(H) =
10^{23}~\rm{cm^{-2}}$. See Korista, Baldwin, \& Ferland (1998) for
details.

\subsection{The continuum SED}

For reasons of simplicity and because we do not know what ionizing
continuum SED is incident upon the broad line emitting gas, I will
demonstrate the effects of a changing continuum SED on the gas emission
using simple power laws, their spectral indeces $\alpha$ ranging from
$-2$ to $-1$ ($F_{\nu} \propto \nu^{\alpha}$).

Panel (b) of Figure~3 shows how the equivalent widths (here, always
measured with respect to the continuum at 1215~\AA\/) of 5 prominent
quasar emission lines vary with the power law index. The gas abundances
were fixed to roughly solar values in all cases. The emission line
equivalent widths ($W_{\lambda}$) all diminish with the softening in
the continuum. This is simply a result of conservation of energy ---
fewer high energy ionizing photons and the resulting lower gas
temperatures produce weaker lines relative to the {\em non-ionizing} UV
continuum at 1215~\AA\/. Panel (d) of Figure~3 shows that ratios
Ly$\alpha$/O~VI and Ly$\alpha$/C~IV are highly sensitive to the
continuum SED, but that N~V/He~II and especially N~V/C~IV are not.
Other SED-insensitive line ratios are $\lambda$1900~blend/C~IV,
$\lambda$1400~blend/C~IV, and N~III~990/C~III~977.

Figure~4 compares our simulations for fixed gas abundances to the
observations of Ly$\alpha$/O~VI (Zheng et al.\ 1995) and
Ly$\alpha$/C~IV (Wang et al.\ 1998). The observational points used
measurements of $\alpha_{ox}$ and $\alpha_{uvx}$, while of course the
simulations have assumed a simple power law with index $\alpha$. The
correspondence is good and demonstrates that the observed $\alpha_{ox}$
and $\alpha_{uvx}$ must stand as reasonably good proxies for the true
ionizing continuum SED.

Zheng \&  Malkan (1993) proposed that the Baldwin effect could be
explained via the combined effect of a strengthening non-ionizing UV
continuum and weakening ionizing continuum with luminosity.  That the
higher ionization lines in Figure~3b have steeper $W_{\lambda}$
dependencies on the spectral index is a reflection of these two
effects, and is generally consistent with the trends in the measured
slopes of the Baldwin effects of various emission lines (Espey et
al.\ 1993; Espey, these proceedings).

A SED--$L_{uv}$ relationship might be driving the Baldwin effect, but
it cannot be the whole story. When the $W_{\lambda}$(N~V) is measured
it has little or no dependence upon luminosity (Espey et al.\ 1993;
Osmer et al.\ 1994; Laor et al.\ 1995).  Figure~1 illustrates this
further. Panels (b) and (d) of Figure~3 predict that N~V/C~IV should be
luminosity independent, i.e., N~V should show a Baldwin effect as
strong as C~IV. Hamann \& Ferland (1993) showed that larger values of
N~V/C~IV and N~IV/He~II are systematically found at larger quasar
luminosities. This must mean that the continuum shape is not the only
parameter.

\begin{figure}
\plotfiddle{korista1fig3.eps}{4.0in}{270.}{80}{80}{-310}{410}
\caption{The equivalent widths of 5 prominent quasar emission lines and
their ratios as a function of metallicity ($Z$; panels [a], [c]) for a
fixed continuum SED ($\alpha = -1.4$), and as a function of the power
law index of the continuum SED ($\alpha$; panels [b], [d]) at a fixed
gas metallicity ($Z = 1$).}
\end{figure}

\begin{figure}
\plotfiddle{korista1fig4.eps}{4.0in}{270.}{80}{80}{-335}{450}
\caption{Observed emission line ratios vs.\ $\alpha_{ox}$ (solid dots)
and predicted emission line ratios vs.\ $\alpha$ (solid lines).}
\end{figure}

\subsection{Metallicity}

Here I will demonstrate the effects of a changing gas metallicity, $Z$,
for a fixed continuum SED ($\alpha = -1.4$). In a simple scheme to
change the metallicity I assumed that the metal abundances relative to
hydrogen scale from their solar values as $Z$. However, nitrogen is
expected to be selectively enhanced by ``secondary'' CNO processing in
stellar populations (at least when Z exceeds a few-tenths solar), and so
in these simulations $\rm{(N/H)} \propto Z^2$ while $\rm{(N/C)} \propto
Z$ and $\rm{(N/O)} \propto Z $ (Hamann \& Ferland 1993; Vila-Costas \&
Edmunds 1993; van Zee, Salzer, \& Haynes 1998). As such, nitrogen lines
are valuable tracers of the gas chemical enrichment (Ferland et
al.\ 1996).

In panel (a) of Figure~3 I show how the equivalent widths of the same
five emission lines vary with the metallicity. The strengths of
Ly$\alpha$ and especially He~II diminish with the increasing metal
opacity. In as far as hydrogen and helium are still the major opacity
sources, their ionized volumes diminish with increasing $Z$ due to the
increased numbers of electron donors. For $Z > Z_{\odot}$, the strengths
of C~IV and O~VI both diminish as the gas cools and the cooling shifts
to other lines.  However, since N~V is initially a weak line and
because the N/H abundance ratio is proportional to $Z^2$,
$W_{\lambda}$(N~V) $\propto Z$! Other initially weak lines are
generally correlated to a lesser extent with metallicity. However,
those lines that are formed near the front of the cloud (e.g., O~IV]
1402) remain roughly constant with $Z$. The temperatures here scale
strongly with $Z$ (i.e., cooler with increasing $Z$) and the He$^{++}$
zone where these lines form is diminshed in volume. In panel (c) of
Figure~3 I show the variations in the same five line ratios with
metallicity. For reasons given above, Ly$\alpha$/O~VI and
Ly$\alpha$/C~IV are relatively insensitive to the gas metallicity.
However, the N~V/He~II and N~V/C~IV ratios are roughly linear in $Z$.
We also find the following rough dependencies on metallicity:
$\lambda$1900/C~IV $\propto Z^{0.6}$, N~III~990/C~III~977 $\propto Z$,
and $\lambda$1400/C~IV $\propto Z^{0.3}$.

It is worth noting that the same line ratios that are insensitive to
the continuum SED are sensitive to $Z$, and those sensitive to the
continuum SED are insensitive to $Z$. While the relationship is not
always simple, this dichotomy of emission line ratio properties is not
a coincidence, and may prove to be a powerful tool in disentangling the
effects of these two fundamental parameters on the quasar spectrum.

\subsection{An explanation?}

Panels (a) and (b) of Figure~3 suggest a solution to the N~V
conundrum.  If, characteristically, the gas metallicity increases (with
N/H $\propto Z^2$) while the continuum SED softens with increasing
luminosity, the $W_{\lambda}$(N~V) will remain approximately constant
and N~V/C~IV and N~V/He~II will increase with luminosity, as observed.

The models do not tell us how to link the SED and $L_{uv}$ together, so
I adopted the empirical $\alpha_{ox} - L_{uv}$ relation found by Wilkes
et al.\ (1994; $H_o = 50~\rm{km~s^{-1}~Mpc^{-1}}$ and $q_o = 0$), and
adjusted slightly the logarithmic zero point to reference the continuum
at 1550~\AA\/. Guided by the N~V/He~II, N~V/C~IV, and $L_{uv}$
relations found by Hamann \& Ferland (1993) I guessed a $L_{uv} - Z$
relation: $Z \approx 1$ for the lowest luminosity quasars (i.e.,
Seyfert nuclei) and $Z \approx 10$ for the highest luminosity quasars.
In this way I adopted the following simple relation to represent the
hypothesis:  \begin{equation} \log Z =
0.183\log\left(\frac{L_{\nu}(1550)}{10^{30}}\right) + 0.477.
\end{equation} After coupling the adopted SED -- $L_{uv}$ and $Z$ --
$L_{uv}$ relations, the photoionization simulations then predict the
slopes, $\beta$, of the Baldwin effect for the different emission lines
and blends. We show these in Table~1, in descending order of the ions'
destruction ionization potentials.

\begin{table}
\caption{Simulated Baldwin relation slopes}
\begin{tabular}{lc}
\tableline
\multicolumn{1}{c}{Emission Line} & \multicolumn{1}{c}{$\beta$} \\
\tableline

O~VI $\lambda$1034 & -0.22 \\
N~V $\lambda$1240 & -0.047 \\
C~IV $\lambda$1549 & -0.20 \\
He~II $\lambda$1640 & -0.17 \\
$\lambda$1400 & -0.13 \\
$\lambda$1900 & -0.083 \\
Mg~II $\lambda$2800 & -0.092 \\
Ly$\alpha$~$\lambda$1216 & -0.10 \\
\tableline
\tableline
\end{tabular}
\end{table}

The predictions in Table~1 follow both the measured amplitudes and
trends of the Baldwin effect slopes reported in the literature
(Baldwin, Wampler, \& Gaskell 1989; Kinney, Rivolo, \& Koratkar 1990;
Zamorani et al.\ 1992; Espey, Lanzetta, \& Turnshek 1993; Osmer,
Porter, \& Green 1994; Laor et al.\ 1995; Zheng, Kriss, \& Davidson
1995; Turnshek 1997). Note, however, that it is difficult to compare
results from different works since the investigators assumed different
values of $q_o$ whose effect is non-uniform in the $W_{\lambda} -
L_{uv}$ plane. Variations in slope due to different assumptions of
$q_o$ are at the 0.04 decade level, comparable to the typically quoted
uncertainties in the measured slopes.

The interplay of the luminosity dependent SED -- $Z$ effects determines
the spectrum, and accounts for the differing Baldwin effect slopes
amongst the different emission lines. While it has been suggested that
the observed strengthening of the Ly$\alpha$/C~IV and
$\lambda$1900/C~IV ratios with luminosity implies a shift to lower
ionization parameters with increasing quasar luminosity (e.g.,
Mushotzky \& Ferland 1984), I would argue for characteristic shifts to
softer continuum SEDs and larger values of $Z$. The predicted changes
in Ly$\alpha$/C~IV with $L_{uv}$ are dominated by a shift in the
continuum SED, while changes in $\lambda$1900/C~IV are dominated by
changes in $Z$.

\section{Discussion}

While I have shown that a SED -- $Z$ -- $L_{uv}$ relationship holds
promise in explaining the origin of the Baldwin effect, the SED --
$L_{uv}$ and $Z$ -- $L_{uv}$ relations or at least their combined
effect must be firmly established and then calibrated, before we can
hope to understand what the Baldwin effect can tell us about cosmology
or cosmic history. I have presented several emission line ratios that
should provide the means to disentangle the effects of the continuum
SED and $Z$.

But why should we expect a SED -- $Z$ -- $L_{uv}$ relationship? We are
beginning to see quasars in the context of the environment of massive
galaxies ($\S~1$), rather than isolated, exotic beasts. The
co-evolution of massive galaxies and the quasars within them may offer
an answer to the above question, and thus may be the fundamental origin
of the Baldwin effect in quasar spectra. Below, I try to summarize this
emerging, yet, I caution, immature picture.

{\em 10 $>$ z $>$ 2.~} The rapid turn on of quasars was the result of
massive galaxy formation, during which the action of merging halos
(Press \& Schecter 1974) both created and drove gas into supermassive
black holes (Haehnelt, Natarajan, \& Rees 1998; Miller \& Percival
1998). Masses of the nuclear black holes, $M_{bh}$, correlated with
masses of forming galaxies or their dark halo masses (Haehnelt et
al.\ 1998).  More massive galaxies, with their deeper central
gravitational wells, were better able to retain gas deposited during
stellar evolution than were lower mass galaxies, so they had a more
rapid buildup of the heavy elements in their nuclear regions (Hamann \&
Ferland 1993; 1999).  The most massive black holes likely produced the
most luminous quasars, {\em initially}, since quasars turned on near or
above their Eddington limits:  $L/L_{Edd} \propto \dot{M}/M_{bh} \geq
1$. These most luminous quasars thus required the largest mass
accretion rates, $\dot{M}$, to sustain their luminosities.  While we
are far from understanding the mechanisms that generate the quasar
continua (see Antonucci 1988; 1998), the standard accretion models
predict softer UV bumps for larger $M_{bh}$ at a fixed $\dot{M}/M_{bh}
\propto L/L_{Edd}$ (Laor \& Netzer 1989; Sincell \& Krolik 1998), and
the current corona -- accretion disk models predict larger UV to X-ray
continuum flux ratios for larger values of $\dot{M}/M_{bh}$ (Czerny,
Witt, \& $\rm{\dot{Z}}$ycki 1996).  Therefore, more luminous quasars
might be expected to have typically softer continuum SEDs (e.g.,
steeper $\alpha_{ox}$).  With the largest required mass accretion
rates, the most massive and luminous quasars are expected to have the
shortest lifetimes (like massive stars), as the limited ``fuel'' supply
is exhausted.  The fuel exhaustion process is not well understood, but
there are some ideas. At sufficient luminosities relative to the disk
binding energy, a ``back reaction'' may occur, choking off the matter
accretion significantly (Haehnelt et al.\ 1998). Furthermore, the
self-similar infall solutions of Bertschinger (1985) and the
Press-Schecter theory of the collapse of dark matter halos (1974) both
predict a decline in the mass accretion rate with time. Even if
$\dot{M}$ were to remain constant, the ratio $\dot{M}/M_{bh}$ must
evolve to smaller values due to accretion alone. Thus the evolution of
an accreting black hole quite naturally leads to a drop in
$\dot{M}/M_{bh}$ with time (see also Haiman \& Menou 1998). Once
$\dot{M}/M_{bh}$ drops below a critical value, the accretion flow
radiation efficiency is expected to drop to values $\eta\/ \ll~0.1$,
and an advection dominated accretion flow (ADAF) may result (Narayan \&
Yi 1995).  In this mode, the UV to X-ray continuum luminosities are
expected to be more comparable. Thus, as $\dot{M}/M_{bh}$ in individual
quasars diminished with time, evolution to lower luminosities and
harder continuum SEDs (e.g., flatter $\alpha_{ox}$) likely resulted
(Caditz, Petrosian, \& Wandel 1991; Yi 1996; Haehnelt et al.\ 1998).

{\em z $<$ 2.~} By this time most of the massive galaxies had already
formed and so the major fire works were over --- most quasars had begun
their death march, dying of exhaustion. The peak in the {\em active}
black hole mass function likely slippled slowly to smaller values with
time (Siemiginowska \& Elvis 1997), and with it the corresponding
required $\dot{M}$, perhaps allowing for longer quasar lifetimes.  The
typical quasar SED would evolve as discussed above. The initial nuclear
gas metallicity would be characteristically smaller in the shallower
potential wells of the lower mass black holes. Further evolution to
smaller values in $Z$ might have occurred in quasars of any mass due to
dilution by lower $Z$ gas falling into the nuclear regions over the
life of the galaxy (Hamann \& Ferland 1993).  Later epoch galaxy
mergers could rekindle the nuclear fires and again enrich the gas,
temporarily extending the quasar epoch. But in a maturing and expanding
universe, galaxy merger rates slowed, and eventually all quasars died.
Remaining are the Seyfert nuclei with their typically smaller $M_{bh}$,
$\dot{M}$, and $\dot{M}/M_{bh}$, lying within the centers of more
modest galaxies, whose nuclear gas metallicities are within a factor of
2 of solar. Also present are the more massive, yet generally dormant,
black holes in the nuclei of present-day giant elliptical and other
galaxies (Margorrian et al.\ 1998).

The above picture is presently sketchy at best, and may be completely
wrong. Key uncertainties are: the lives of quasars and their energy
dissipation mechanism(s), galaxy formation and its relation to the
births of the quasars. If this picture has any validity it is easy to
see why the Baldwin effect has been such a hard nut to crack.
Different initial conditions and times of massive galaxy formation
resulted in separate quasar populations, characterized by $M_{bh}$,
$\dot{M}$, $\dot{M}/M_{bh}$, and $Z$, co-evolving in time. This
scenario suggests that these four parameters are fundamental to the
quasar in driving its continuum and emission line spectra.  Osmer,
Porter, \& Green (1994) demonstrated that the Baldwin effect is seen
globally across all redshifts {\em and} within narrow redshift
intervals. This would argue for the study of carefully selected samples
of quasars at many redshift or lookback time intervals.

Further sources of intrinsic scatter in any continuum -- emission line
correlations can arise from geometry -- orientation effects due to
possible anisotropies in the line and continuum emission (Netzer 1987;
Netzer et al.\ 1992), as well as short time scale continuum -- emission
line reverberation (Kinney et al.\ 1990; Pogge \& Peterson 1992;
Shields, Ferland, \& Peterson 1995). The latter should be less
important for the higher luminosity quasars.  These and other nuisances
should ``average out'' in large samples of objects --- and perhaps the
luminosity distances can be recovered through a process analogous to
``main sequence fitting'' of stellar clusters, as suggested by Jack
Baldwin at this meeting. For a thorough discussion of these issues, see
the conference proceedings contribution of Osmer \& Shields.

While the jury is still out as to whether the Baldwin effect in quasar
spectra can ever be used as a standard luminosity candle, it may be
telling us about the workings of the quasar central engine, the gas
enrichment in the nuclei of massive galaxies, and possibly about the
evolution of massive galaxies themselves.

\acknowledgments

Jack Baldwin and Gary Ferland are my collaborators in the work
presented here. I also acknowledge helpful discussions with Fred Hamann
and Joe Shields. I am grateful to the Baldwin family for their gracious
hospitality during my two week stay in La Serena.


The supernovae associated with gamma-ray bursts

Abstract: The connection between long GRBs and supernovae is now well established. I briefly review the evidence in favor of this connection and summarise where we are observationally. I also use a few events to exemplify what should be done and what type of data are needed. I also look at what we can learn from looking at SNe not associated with GRBs and see how GRBs fit into the broad picture of stellar explosions.

\section{Introduction}

The connection between gamma-ray bursts (GRBs) and supernovae (SNe)
has existed for longer than the ``official'' history of GRBs.  Colgate
\cite{Colgate68} predicted that exploding stars should emit $\gamma$
rays. Even though it turned out that this model for GRB prompt
emission is not correct, it meant that a connection between SNe and
GRBs was on everybody's mind from the beginning. Indeed, Klebesadel,
Strong \& Olson \cite{KSO73} searched for spatial and temporal
coincidence between explosive events (novae and supernovae) and
GRBs. Although they didn't find any coinciding event, this was not
enough to discourage further similar searches.  It is only in the
afterglow era that major progress could finally be made and that a
clear link between long GRBs and SNe would be firmly established.

This review is not an attempt to be exhaustive so there is no
discussion of all known or candidate SNe associated with GRB. I will
rather try to show where we are now as a community and to place these
SNe in a general context. Several recent reviews took a close look at
the evidence for each SN/GRB and should be consulted for specific
information about any particular object (see e.g. \cite{WB06, HB12}).

\section{From suspicion to proof}

GRB980425 was quickly associated with a supernova, SN~1998bw. The
spectra of the SN clearly showed that it was a type Ic -- meaning that
the spectrum showed no sign of hydrogen nor helium. The association
between the GRB and the SN meant that the burst was extremely faint
however ($E_{iso} < 10^{48}$ erg).  SN~1998bw itself was also
peculiar, as it was a very bright SN, showed very broad lines (BL) and
had a large expansion velocity (e.g. \cite{Galama98}.  Because of this
highly unusual nature -- very weak and very nearby while the SN was
very bright, there was still a question mark over this event and some
fraction of the community was still waiting for a GRB at a genuinely
``cosmological'' distance.

This happened with GRB030329. The burst was close enough that a SN
could be seen, yet it was distant and bright enough that it could be
classified as a genuine GRB. Unambiguous SN features became visible in
afterglow spectra a few days after the burst (\cite{Stanek03,
  Hjorth03}). The SN was shown to be very similar to SN~1998bw
\cite{Matheson03} and was also classified as a SN~Ic-BL.  This event
removed any doubt as to the association between SNe and GRBs.

Since then, several other GRBs have had indisputable spectral evidence
supporting the presence of a SN (see Table~\ref{tab:sne}). Several
good quality spectra (between 2 and many) have been obtained also for
GRB031203 (SN~2003lw \cite{Malesani04}), XRF060218 (SN~2006aj
\cite{Pian06}), GRB100316D (SN~2010bh, \cite{Chornock10, Bufano12}),
GRB120422A (SN~2012bz, \cite{Malesani12, Sanchez12})\footnote{This
  last event is still a developing story as of this writing. While the
  presence of a SN is clear, the SN itself is not yet fully
  characterised.}  These SNe share the characteristics of SN~1998bw,
in the sense that they are of type Ic and they have broad lines,
indicative of high expansion velocity.

There were attempts to obtain a spectrum of the SN for a number of
other GRBs but the evidence is not as solid as one would wish.  The
reader is referred to \cite{HB12} for a detailed discussion of each
event and associated evidence.

\subsection{Light curves bumps}

Short of obtaining a spectrum around the time of the SN peak
luminosity, a well-sampled light curve in several filters can also
show evidence for a SN. This will take the form of a ``slowing `down''
of the light curve or even a rebrightening. A fairly large number of
bursts have displayed such a behaviour and in most cases, the bump can
be reasonably well fit with a SN light curve (see
Table~\ref{tab:sne}).  These SNe are represented as ``modified
SN~1998bw'', in the sense that the adopted model light curve is
usually that of SN~1998bw, stretched or compressed in time and shifted
in brightness (and properly accounting for the different luminosity
distance and \textit{k}-correction).

The brightness of these bumps is usually within a factor of a few of
SN~1998bw (e.g. \cite{Zeh04, Ferrero06, Cano11a}) albeit slightly
fainter on average, which is not surprising given that SN~1998bw was a
very bright event.
At this point, one can say that in nearly every case where we should
have been able to see a bump, a bump has been seen and it corresponds,
in brightness and timing, to a SN.

\begin{table}
\begin{tabular}{l}
\hline
\hline
\textbf{a)} GRBs with highest quality spectroscopic evidence:\\
GRB980425 $\equiv$ SN~1998bw, GRB030329 $\equiv$ SN~2003dh, GRB031203
$\equiv$ SN~2003lw, \\
GRB060218 $\equiv$ SN~2006aj, GRB100316D $\equiv$ SN~2010bh, GRB091127
$\equiv$ SN~2009nz, \\
GRB120422A $\equiv$ SN~2012bz \\
\hline
\textbf{b)} Other GRBs with SN evidence: \\
GRB970228, GRB980326, GRB990712, GRB991208, GRB000911, GRB011121,\\
GRB020305, GRB020405, GRB020410, GRB020903, GRB021211, GRB030723,\\
GRB040624, GRB041006, GRB050416A, GRB050525A, GRB050824, GRB060729,\\
GRB070419A, GRB080319B, GRB081007, GRB090618, GRB100418A, GRB101219B,\\
GRB101225A\\
\hline
\hline
\end{tabular}
\caption{A list of GRBs where some evidence for a SN has been found
  (see Table~9.1 in \cite{HB12} for references related to any
  particular GRB). The first category includes those events with a
  very clear spectroscopic signature; the second category includes
  event where the evidence ranges from good to poor.}
\label{tab:sne}
\end{table}

\subsection{Properties of SNe of GRBs}

We want to understand the explosions so that we can ultimately
constrain the progenitors of GRBs. A good multi-colour light curve,
coupled with a velocity measurement at maximum light, can be modelled
analytically \cite{Arnett82} to obtain explosion parameters: kinetic
energy $E_k$, mass of the ejecta $M_{ej}$, and mass of $^{56}$Ni
synthesised. This has been done for a few events (see \cite{Cano11b,
  Olivares12, Bufano12} for the case of XRF100316D). This approach may
be subject to systematic errors however and modelling based on light
curve and time-series spectroscopy yields more robust results
(e.g. \cite{Iwamoto98, Mazzali03, Maeda07}).  This provides a better
understanding of the ejecta's structure and yields more reliable
results than simple light curve analysis because it avoids several
simplifications.
The results of such modelling give a range of values for each parameter:
the kinetic energy is $E_k = 2-60 \times 10^{51}$ erg, the ejected
mass ranges from 2 $M_\odot$ to 13 $M_\odot$ and the mass of $^{56}$Ni
synthesised in the explosion is between 0.1 $M_\odot$ and 0.7
$M_\odot$.

\section{Cautionary tales}

In some cases, the observational data may not allow us to detect a SN
signature.  One such example is GRB070419A. While there was an initial
claim for a SN \cite{Hill09} based on a photometric bump, a full
analysis of all the light curve data showed that the late time
behaviour could be fully explained by a fairly shallow decay of the
afterglow \cite{Melandri09}. In order to determine the SN brightness,
we need to know very well what the afterglow is doing. This means that
a well-sampled GRB light curve is necessary for tens of days; not only
for a few hours.  Observing the source past the SN maximum is needed
as well; if any ``rebrightening'' is seen, the fast fading of the SN
will leave a different signature in the light curve than the power
law decay of the afterglow.

Another example is given in Fig.~\ref{fig:lc0329}. The left panel is a
light curve displaying a more or less power law decay. From this light
curve, it is impossible to tell whether there is a SN hidden in the
data or not.
The right panel of the figure shows the full light curve, which is for
GRB030329. Now, obviously a SN was there to be found in the data.
Assume for a moment that this burst had been at $z\sim 0.5$, hence
much fainter, we would have had much less data, the light curve would
have looked like the left panel of Fig.~\ref{fig:lc0329} and it would
have been very hard to find a SN in the light curve. The point here is
that even if we do not see an obvious SN bump, there may well be one
but we need good enough data to separate the relative contributions of
the afterglow and SN. The case of GRB060729 is particularly
instructive: there is no SN bump to be seen simply because the
afterglow was very bright but a careful analysis revealed a bright SN
in the light curve \cite{Cano11a}.

\begin{figure}
\includegraphics[angle=270, scale=0.55]{dual_lc.eps}
\caption{\emph{Left:} An $R$-band light curve for a seemingly
  anonymous GRB. The afterglow decays by several magnitudes over $\sim
  2.5$ decades in time. At first sight, there is no SN bump.
  \emph{Right:} The whole $R$-band light curve of GRB030329 which is
  the event on the left. There is no doubt that there was a SN
  associated with this event.}
\label{fig:lc0329}
\end{figure}

Another example is provided by GRB101225. There are two ``competing''
interpretations regarding the nature of this event.  One is a
cosmological GRB ($z\simeq 0.33$) with a SN bump \cite{Thoene11}, the
other is a minor planetary body tidally disrupted by a (Galactic)
neutron star \cite{Campana11}.
The point here is not to decide which interpretation is best but to
emphasise again that \emph{the right type of data at the right time}
will go a long way towards removing any ambiguities in the analysis.
This burst serves as a stark reminder that nothing replaces
spectroscopy.

\section{GRBs or SNe we do not see}

The central point of this section is about negatives: not finding what
we expect or finding what we do not expect.

\subsection{Observing a GRB and finding no SN}

There are times when, try as you may, there is no SN to be found. Two
events in particular, GRB060505 and GRB060614, are now etched in the
collective consciousness.  Both were at a redshift low enough that a
SN even substantially fainter than SN~1998bw would have been detected
easily. The limits on the brightness of any SN are actually very
strong: in both cases any SN would have to have been at least 100
times fainter than SN~1998bw at peak (e.g. \cite{Fynbo06, Gal-Yam06,
  DellaValle06}).  Several other events also showed no sign of a SN,
although in some cases the lack of redshift makes the evidence weaker
(\cite{Soderberg05, Levan05}).

While there may still be some lingering questions regarding the nature
of these bursts, they may be in a different class than most long GRBs.
A possibility is that they have slightly different
progenitors. Indeed, it may well be possible to make a GRB (i.e. a
relativistic jet) without making a bright SN (e.g. \cite{Woosley93,
  Tominaga07, Moriya10}). A bright burst of $\gamma$ rays may be the
only electromagnetic signature of such events.  Obviously we need to
find more of these ``failed'' SNe. This is why it is important to
monitor low redshift GRBs for several weeks. The absence of a SN may
tell us even more than an actual detection.

\subsection{Observing a SN and serendipitously finding a GRB}

Given that GRBs emit their radiation in a narrow cone, we miss most of
them (i.e. those pointing away from Earth). The emission becomes
visible however at late times in the radio when the afterglow is
emitting essentially isotropically.  Because we can separate normal
SNe from late-time jet emission via their luminosities (jets are much
more luminous than SNe, e.g. \cite{Soderberg10} and references
therein), this offers a way of finding GRBs whose high-energy emission
was directed away from us.

Observations of SNe at radio wavelengths have recently uncovered such
a case.  The radio luminosity of SN~2009bb places it squarely among
the jet-associated SNe \cite{Soderberg10}, yet no GRB was
observed. Another possible case is SN~2007gr \cite{Paragi10}.
These recent observations of nearby SNe show that the number of jets
found via their SN emission may be comparable to the number of SNe
found via their early gamma-ray emission.  If the two techniques have
comparable yields, they should both be pursued as each complements the
other.

\section{The SNe of GRBs in a general context}

From their relative rates, we know that GRBs represent a very small
fraction of all SNe. From their observed properties however, they do
stretch the parameter space of stellar explosions.
We know that local ``normal'' SNe and GRBs mark the explosions of
massive stars. What sets GRBs apart from other SNe is the large amount
of energy that goes into relativistic ejecta.  What we ultimately want
to know is why a small fraction of SNe manage to do that while most of
them do not. Comparing the respective observational properties of these
two types of explosions is a fruitful way of understanding their
differences.

One similarity is the spectroscopic type, Ic, which tells us that GRB
progenitors are stripped-envelope stars.  They also share the same
``locations'' on the respective host galaxies (\cite{Fruchter06,
  Kelly08}): local SNe~Ic and GRBs are much more concentrated on the
light of their hosts than other types of SNe.  The SNe of GRBs show,
on average, larger ejection velocities than local SNe~Ic, they have
very broad lines. They also have larger energies and peak brightnesses
than local SNe~Ic (see Fig.~\ref{fig:peak_mag}).

\begin{figure}
\begin{center}
  \includegraphics[width=.6\textwidth]{peak_mag.eps}
\end{center}
\vspace*{-10mm}
\caption{Cumulative distribution of peak magnitudes for local SNe Ibc
  (red dashed line) and the SNe associated with GRBs (blue solid
  line). On average, the SNe of GRBs are brighter than local SNe Ibc
  (based on data in \cite{Cano11a}). This should be interpreted with
  some caution however as observations of a SN/GRB at an intermediate
  redshift (say $z\simeq 0.3-0.5$) is easier when the SN is bright. }
\label{fig:peak_mag}
\end{figure}

The properties of GRB host galaxies are instructive as well; they are
different from those of SNe~Ic. GRBs are found in relatively low
luminosity irregular galaxies \cite{Fruchter06}. GRB are also found in
more metal-poor environments than SNe~Ic (e.g. \cite{Kelly11,
  Modjaz11}, see also \cite{Savaglio11} however). In this context, it
is interesting to note that local \emph{broad-line} SNe~Ic tend to
prefer low metallicity environments as well (\cite{Modjaz11}).

\section{Conclusion}

A bump on a light curve in itself is moderately interesting; it
becomes much more valuable when combined with time-series spectroscopy
because this allows us to characterise the progenitor star.
Photometric monitoring low-$z$ GRB is warranted however because this
is the only way to find other SN-less GRB (like GRB060614). It is
important to know what fraction of bursts can be formed via this
channel.  This will eventually help us understand whether these
objects are a separate class altogether or part of a continuum
encompassing very faint and very bright SNe.

It can be frustrating to wait for nature to be so kind as to offer us
a low redshift GRB where the SN can be studied in detail. Observing
nearby SNe can help us tremendously to understand GRBs in general.
Radio monitoring of SNe has been proved to be a good way to find
relativistic ejecta in the local universe. On the other hand, ongoing
optical surveys of the sky (PTF, Pan-STARRS, SkyMapper, La Silla
QUEST, CRTS, etc.) already find hundreds of SNe of all types every year.
This will lead to a better understanding of the evolution and
explosion of massive stars that will improve our understanding of
GRBs.
This is like making a jigsaw puzzle with a 1000 pieces with one blank
piece among them. Obviously we do not know where it fits in the global
picture but once the other 999 pieces are in place, we will will know
exactly where the blank one fits and, looking at the general picture,
we will have a very good idea as to what this piece should look like.

High Energy Astrophysical Phenomena (astro-ph.HE); Cosmology and Nongalactic Astrophysics (astro-ph.CO)


Colors of Inner Disk Classical Kuiper Belt Objects

Abstract: We present new optical  broadband colors, obtained with the Keck 1 and Vatican Advanced Technology  telescopes, for six objects in the inner classical Kuiper Belt. Objects in the inner  classical  Kuiper Belt are of interest as they may represent  the  surviving  members of the  primordial  Kuiper Belt that formed interior  to the  current  position of the 3:2  resonance  with  Neptune,  the current  position  of the  plutinos,  or,  alternatively,  they may be objects formed at a  different  heliocentric  distance  that were then  moved to their present   locations.  The  six  new  colors,  combined  with  four  previously published,  show that the 10 inner  belt  objects  with  known  colors  form a neutral  clump and a reddish  clump in B-R  color.  Nonparametric  statistical tests show no significant difference between the B-R color distribution of the inner disk objects compared to the color  distributions of Centaurs, plutinos, or scattered disk objects.  However, the B-R color  distribution  of the inner classical Kuiper belt objects does differ  significantly from the distribution of colors in the cold (low  inclination)  main classical Kuiper belt. The cold main classical  objects are predominately  red, while the inner classical belt objects are a mixture of neutral and red. The color  difference may reveal the existence of a gradient in the composition and /or surface  processing history in the primordial Kuiper Belt, or indicate that the inner disk objects are not dynamically analogous to the cold main classical belt objects.

\section{Introduction}

Ongoing discovery and orbital and physical characterization of
Kuiper belt objects (KBOs) and related outer solar system
bodies promise new insight into the formation and evolution of
our solar system.
However, the evolving dynamical picture shows that the orbits of many,
if not most,
objects
have been substantially changed over time.
Thus it may be difficult to study any  differences between
objects that formed at different heliocentric distances, as the
current heliocentric distances of individual objects may differ radically
from the
heliocentric distances at the time of their formation.

\citet{Gomes03} hypothesized a dynamical mechanism to move the
dynamically hot (high
values of orbital
inclination (\iital) and eccentricity (\eital))
KBOs from a massive protoplanetary disk
truncated at
30$-$35 AU from the Sun
outward to the present Kuiper Belt region.
This mechanism would not work for the cold KBO population
(low values of orbital
inclination (\iital) and eccentricity (\eital))
\citep{SSBNMorbiAJ}.
In the Gomes model, the cold KBO population would either have formed
{\it in situ} or would have been moved by a different
mechanism.
In the context of the Nice model of the outer solar system,
a mechanism has been proposed that could move even cold KBO
objects from a disk truncated at $\sim$ 35 AU to their present
distances \citep{LevisonNice08}.
In this model, the dynamically cold KBO objects would have formed in the
outer
regions of the massive protoplanetary disk, perhaps in the region
around 30 AU.

The vast majority of the known
dynamically cold classical belt objects have semimajor axes
(\aital) between 42 and 48 AU.
There are also low \eital, low \iital
classical objects with \aital $<$ 39.4 AU.
This semimajor axis is less than
the current position of the Neptune 3:2 resonance,
the semimajor axis at which Pluto and
the plutino KBOs are currently found.
The population of these objects, known as inner
classical
belt KBOs (ICKBOs), with \aital between about 36 and 39.4 AU, is much
lower than
that of the main classical belt objects ( \aital between 42 and 48 AU).
Because of the relative dearth of known classical belt objects interior to
the
Neptune 3:2 resonance,
no published physical studies have specifically targeted these
objects.
This class of objects extends
the \aital range of objects  compared to the
\aital range spanned by the main belt alone.
Thus it is
potentially useful in the
study of any gradients of chemical or surface processing
history with formation distance for minor bodies in the outer solar
system.
Here we present new optical colors of a sample of inner classical
KBOs.
We compare the color distribution of these objects with
that of other classes of outer solar system objects.

\section{The Inner Classical Kuiper Belt}

A dynamical classification scheme for minor bodies in the outer solar
system
is presented by \citet{GlaMarVan08} (hereafter GMV).
They define  inner classical belt objects as those that have
semimajor axes (\aital)\ less than  39.4
AU, are not Centaurs, are not scattered disk objects, and are
not in resonant orbits.
As explained in GMV, the stable inner classical belt is not disconnected
from
the main classical belt, which contains
objects on stable, non-resonant orbits with
\aital = 42$-$48 AU.

GMV only classify objects
that met certain criteria for having sufficient
astrometry as of May 2006.
They list 17 as inner classical Kuiper belt objects
(hereafter ICKBO) objects and over 250  main classical Kuiper
belt objects (hereafter MCKBOs).
As more recent astrometry is available, we searched  the
continuously updated listing
\footnote{http://www.boulder.swri.edu/$\sim$buie/kbo/kbofollowup.html}
of
the Deep Ecliptic Survey (DES) team
\citep{Elliot05} for possible ICKBOs in
addition to those identified in GMV.
We searched for objects
with \aital $<$ 39.4 AU which the DES team definitively classifies as
``classical" objects based on their detailed orbital integrations.
Two additional members of the inner classical belt
(144897 and 2003 QA92) were found
in this way.
Neither of these objects
are classified in GMV.

Of the 17 ICKBOs identified in GMV, five are classified by the
DES team as ``scattered near" objects, rather than classical
objects.
All five of these objects have orbital inclinations over 10\degr.
Thus, there is some disagreement over the classification of some
objects in the \aital $<$ 39.4 AU region.
However, all but one of the objects for which
colors are available have orbital inclinations less than 10\degr,
and
so are probably members of the inner classical belt.
If we use the GMV classification, augmented by the 2 additional objects,
the inner classical belt has 19 known objects.
\citet{CFEPS09} estimates that the inner disk KBOs may
have a population 10 to 20 times smaller than the main belt.

Evidence that the main classical Kuiper Belt is composed of a ``cold"
and ``hot" population has been presented by \citet{Brown01}.
There is some evidence of differences in the physical properties
between
the
cold and hot populations (see \citet{Pei08} and references therein).
However the dividing line in inclination between the cold and hot
populations is not sharp, and indeed a simple dividing line may not
be useful, as the two populations may overlap.
\citet{Brown01} models the two populations as separate Gaussians in
inclination, with the cold population having \sigmac \ =
2.2\degr \ and the hot  having \sigmah \ around 17\degr.
\citet{Gulbis06} use a inclination of about 5\degr \ to separate
cold (core) and hot (halo) classical objects.
However,
\citet{Pei08} find that the colors of classical main belt objects
are uniformly red up to an inclination of 12\degr- that is, they do
not see a break in colors at 5\degr.

Our initial goal was to measure the color of ICKBOs with the lowest
values of \iital and \eital, as these could be an extension of
the main cold classical belt towards the Sun.
If this were the case, the inner objects would extend the semimajor axis
range of the classical belt and provide a larger range of semimajor axis
in which to look for correlations between semimajor axis and
physical properties than the \aital range provided by the main cold
classical belt KBOs alone.
However, new results on the structure of the Kuiper Belt obtained
from surveys analyzed with observational biases taken into account
\citep{CFEPS09}
indicate the possibility that the inner belt objects are not analogous
to the cold main belt
but are perhaps more analogous
to the hot classical belt objects.
Such hot objects originated at  significantly different
heliocentric distances compared to their present locations.
\citet{CFEPS09} argue, from the number of expected and observed
inner disk objects at low \iital,
that the inner disk is most
likely devoid
of a cold component, but this conclusion is uncertain due to
the small number of inner disk objects found so far in their
survey.

In contradiction to the \citet{CFEPS09} results,
\citet{LykMuk07}
specifically posit a {\it cold} inner disk.
These authors  state that cold
classical KBOs are located in the inner disk region- 37 AU $<$ a $<$ 40AU (q $<$ 37 AU)-
as well as, of course, in the main belt region
(42 AU $<$ a $<$ 47.5 AU).

Thus, to summarize, the inner classical disk
may be
an extension of
the cold main classical belt, or it may be a hotter component that will
not
tell us directly about any  changes in physical properties
as a function of formation radius in the cold population.
To the extent that optical colors can be related to formation
radii, colors of inner objects might give a hint as to the
origin of this population.
More definite dynamical information on the origin of the inner objects
must await
a larger sample of objects.
Also, careful modeling of the observational biases must be
taken into account \citep{CFEPS09}
to specify whether the ICKBOs are part of a dynamically hot or cold
population.

\section{Observations}

BVR imaging observations of five ICKBOs were obtained with
the Low Resolution Imaging Spectrometer (LRIS) \citep{Oke95}
on the Keck 1 telescope on Mauna Kea on
20 September 2006.
We operated LRIS in its dual-channel imaging mode that allows
simultaneous imaging in a blue and a red channel.
We took 300 second images through a B (438 nm)
filter on the blue side, and alternated between a V (547 nm)
and an R (642 nm) filter on the red side, with a dichroic element
(D460) separating incoming light into blue and red channels.

There were some very intermittent light clouds during our
observations.
Repeatability of magnitudes in the B images indicate that the
clouds attenuated
light by no more
than a few tenths of a magnitude at worst.
Clouds are almost true neutral absorbers in the optical
region \citep{Serkowski70, Honeycutt71}, with B$-$V and
V$-$R colors made bluer by no more than 0.01 mag for
1 mag of extinction.
Thanks to the dual-channel imaging operation mode, our colors should not
be affected to any
significant extent by the presence of these thin clouds.
Observations of Landolt standard stars \citep{Landolt92}
were used to derive
color transformation equations to the
Johnson- Kron- Cousins system, as in our previous
observations with this instrument
\citep{TegRom03}.
Nine high quality standard stars from
three Landolt fields
(SA 110,
PG 2213$-$006 and SA 95)
were used to derive transformation equations between
instrumental and system colors.
These transformation equations were indistinguishable,
within errors, to those shown in \citet{TegRom03}.
The solar color in this system is B$-$R = 0.99 \citep{Pei08}.

New colors of 119951 reported here were obtained
with  B (450 nm), V (550 nm), and R (650 nm)
glass filters in front of a 2048 $\times$ 2048 pixel
CCD camera at the f/9 aplanatic Gregorian focus of
the 1.8-m Vatican Advanced Technology Telescope
(VATT; the Alice P. Lennon telescope and Thomas J.
Bannan facility) on Mt. Graham, Arizona.
We binned the 15 micron pixels 2 $\times$ 2,
yielding 1024 $\times$ 1024 pixels, covering 6.4 $\times$ 6.4
arcmin of the sky at 0.375 arcsec per pixel.
Observations of 119951 were obtained on 7 Jun 2005,
during a 5 night run of photometric weather.
Multiple observations of
sixteen stars in three Landolt fields (PG 1633+099,
PG 1323$-$086  and  SA 110) were used to derive transformation
equations.
The scatter of the individual colors of the standard stars
with respect to the equation was typically less than
0.01 mag.

Results of the new observations are reported in Table 1.
The first column gives the name of the object, if one has been
assigned, otherwise the provisional designation is listed.
The second column gives the permanent number of the object,
if one has been assigned.
Other columns give the optical BVR colors and the error
in the color ({$\sigma$} is the dispersion in the individual color
measurements and {\it n} the number of independent measurements).
V magnitudes are quoted to only a tenth of a magnitude,
as thin clouds may have
affected these measurements.
Also, the time span over which each objects was observed was between
1 and 2 hours, not long enough to get a
proper rotational magnitude sampling for objects with substantial
lightcurve variations.
The last three columns give the semimajor axis (\aital),
eccentricity (\eital) and inclination (\iital).

We also searched for previously published colors
of ICKBOs, primarily using the MBOSS online database
\footnote{http://www.sc.eso.org/$\sim$ohainaut/MBOSS/ }
\citep{HaiDel02}.
Four additional objects were found, and for convenience their
BVR colors are listed in Table 1.

Although the sample of 10 colors seems small, it is about half
of the total known population of ICKBOs.
The remaining objects are small and mostly very faint,
and will require significant time on very large telescopes to obtain
good colors.
A histogram of the B$-$R colors of the objects in Table 1 is shown in
Figure 1.

\section{Colors of Inner Disk Objects}

The histogram of B$-$R colors for the ICKBOs shown in Figure 1
certainly suggest the possibility of a bimodal
color distribution.
However,
the ``dip test" \citep{Har85,Pei03} does not indicate a
statistically
significant bimodal signal for the present sample (~50\% rejection of
unimodal hypothesis).
The ``dip test" may, in a sense, be too restrictive
when the question
is whether or not
two classes of
objects are present.
It {\it is} possible to combine two distinct classes of objects
and get a combined distribution with a {\it single} mode \citep{Zhu07}.
The definitive discussion of unimodality or bimodality of
the inner disk object colors must wait until a
larger sample of colors is available.

We use two different nonparametric or distribution-free statistical
tests to compare the B$-$R color distributions
of ICKBOs with samples of other classes of outer solar system objects.
The Wilcoxon rank-sum test,
also known as the Mann- Whitney  test,
calculates the probability that two populations
have the same mean, so could arise from the same parent population
\citep{AldRoe77,Langley71}.
This test  compares the sum of the rankings of the
values of the individual samples
in a pooled, ordered union of the two samples.
The test is useful for comparing samples which are not
distributed normally, as the test requires no assumption of
normal or Gaussian distributed populations.
We also use the more familiar Kolmogorov- Smirnov (or K-S) test
\citep{PressNumRec}.
We use
implementations of these tests found in the
statistical package R \footnote{http://www.R-project.org}.

We first compare the B$-$R colors  of the 10 ICKBOs with B$-$R colors
to samples of 26 Centaurs and 17 scattered disk
object (SDOs)
from \citet{SSBNcen}.
Next we compare the ICKBO sample to the sample
of 41 plutinos with B$-$R colors \citep{RomTeg07}.
Histograms of the B$-$R colors of the samples of
Centaurs, SDOs and plutinos used are shown in Figure 1.
The Wilcoxon test shows that the probability that the distributions
of 10 ICKBOs and the 26 Centaurs have the same mean color is 90\%.
For the ICKBOs and the plutinos the probability is 36\%, and
for the ICKBOs and the SDOs the probability is 80\%.
Thus, we find no evidence of any statistically significant difference
between the average color of
ICKBO and Centaur or SDO B$-$R color distributions.
While there is a greater difference between the ICKBO
and the plutino B$-$R color distributions, it is only at the one sigma
level and so  is not statistically significant.

If the inner objects are an extension of the cold main belt population,
then any difference in color distribution could indicate a radial
gradient in physical properties.
As discussed earlier, an alternative view is that the inner disk
objects are not part of the cold classical disk.
If this is so, the
objects we have studied are just the low \iital and \eital
members of a broader inclination distribution akin to the hot main belt
objects.
Main belt KBOs, particularly those that have low \eital and \iital,
are very predominately red \citep{TegRom00,TruBro02,Dor02},
with  B$-$R colors typically around 1.7.
From the results of \citet{Pei08},
hereafter PLJ,
we first use B$-$R colors of objects up to \iital of 12\degr \
as the cold main belt population distribution to compare to the
ICKBO colors.

Our first comparison is between the B$-$R colors of the 10 objects in Table 1
and
the 46 objects from PLJ with \iital up to 12\degr.
The Wilcoxon test shows that the probability that the distributions
of B$-$R colors of 10 ICKBOs and the 46 objects from PLJ
have the same average color  is 0.1\%.
The
Kolmogorov-Smirnov two-sample test
applied to the same samples also
indicates  the probability that the distributions
of B$-$R colors of 10 ICKBOs and the 46 objects from PLJ
have the same average color is 0.1\%.

Next we restrict the samples to \iital $<$ 5\degr, as several groups
propose this as a boundary between cold and hot classical populations
\citep{Gulbis06,Noll2pop08}.
This leaves  7 objects from Table 1 and 34 objects from PLJ.
The Wilcoxon test shows that the probability that the distributions
of the 7 ICKBOs and the 34 objects from PLJ
have the same average color is 1\%.
The
Kolmogorov-Smirnov two-sample test
applied to the same samples
indicates  probability that the distributions
of 7 ICKBOs and the 34 objects from PLJ
have the same average color is 2\%.
The decrease in sample size of the \iital $<$ 5\degr \ samples
presumably accounts for the slight increase in probability for the
smaller samples.

\section{Discussion}

When we began this particular project,
we thought that the low \eital and \iital inner belt objects were a
natural sunward extension of the cold classical population found
between 42 $<$ a $<$ 48 AU at low \eital and \iital . Our original aim
was simply to
measure colors for a sample of these objects to compare the colors with
the cold main belt objects.
As discussed in Section II,
the dynamical status of the low
\eital , \iital inner disk objects is a major question. The importance
of
these
objects is only now being recognized; for example, the Nice model of
\citet{LevisonNice08} said nothing specifically about the inner
classical
disk.

Much as we would like to have a definite answer to the precise
relationship between the inner-classical KBOs and the cold-classical
KBOs, neither our color data nor the dynamical models can resolve this
issue at present. One of the main issues is the dynamical origin of the
inner classical objects, as discussed in section II above. So long as
this issue remains unsettled, we are left with admittedly frustrating
result of two different possible ``conclusions" for how our colors may
relate these objects to the other populations. Though of course more
color data are always welcome, and could perhaps provide a stronger
statement about the connection, or lack of connection, with the plutino
population, that will have to wait for further discoveries of brighter
inner classical objects amenable to measuring colors. In the meanwhile,
we await further developments from the dynamicists to propose a possible
resolution to this dilemma.

\section{Conclusions}

The B$-$R color distribution of a sample of 10 inner disk Kuiper Belt
objects is shown to be not inconsistent with the color distribution
of samples
of plutinos, Centaurs and scattered disk objects.
The current inner disk sample has both red and neutral colored objects,
reminiscent of the Centaurs,
which have a bimodal B$-$R color distribution.
However,
we cannot claim that the ICKBOs have a statistically significant bimodal
signal with the present data.

The average B$-$R color of the ICKBOs are inconsistent, at the 98\% level
or higher, depending on specific sample and statistical test used,
with the colors of a sample of cold (low inclination) classical
KBOs with semimajor axes between 42 and 48 AU.
Possible conclusions of this are:

1) The inner disk objects we observed, even though they are of
low \eital and \iital, are members of a population analogous to the
{\it hot}
classical KBOs.
The hot classical KBOs do not show a predominately red optical
color distribution, as do the cold classical objects.

2) If the inner disk objects are in fact members of a cold inner disk
population which is a sunward continuation of the cold classical
objects between 42 and 48 AU, then there is a radial color
gradient in the colors of this cold disk population of KBOs.

\acknowledgements

We thank the NASA Planetary Astronomy program for financial support of
this research and the NASA Keck and Vatican
Observatory telescope allocation committees for consistent allocation
of telescope time.
We thank Brett Gladman and J.J. Kavelaars for discussions of the inner
classical
belt which helped motivate this work.

\vfil\eject

Earth and Planetary Astrophysics (astro-ph.EP)


Tuning the superconducting and magnetic properties in Fe_ySe_0.25Te_0.75 by varying the Fe-content

Abstract: The superconducting and magnetic properties of Fe$_{y}$Se$_{0.25}$Te$_{0.75}$ single crystals ($0.9\leq y \leq1.1$) were studied by means of x-ray diffraction, SQUID magnetometry, muon spin rotation, and elastic neutron diffraction. The samples with $y<1$ exhibit coexistence of bulk superconductivity and incommensurate magnetism. The magnetic order remains incommensurate for $y\geq 1$, but with increasing Fe content superconductivity is suppressed and the magnetic correlation length increases. The results show that the superconducting and the magnetic properties of the Fe$_{y}$Se$_{1-x}$Te$_{x}$ can be tuned not only by varying the Se/Te ratio but also by changing the Fe content.

\section{Introduction}
The iron chalcogenide family of superconductors Fe$_{y}Ch$ ($Ch=$~Se/Te) was discovered in 2008\cite{Hsu08} shortly after the first report of high temperature superconductivity (HTS) in the iron pnictides.\cite{Kamihara08} This family stands out because of its simple crystal structure relative to other Fe-based superconductors.\cite{Rotter08,Wang08,Ogino09,Zhu09} In common with the other Fe-based HTS the parent phase, FeTe, exhibits antiferromagnetic order,\cite{Fruchart75,Bao09,Li09} and superconductivity appears only upon substitution of Te with Se or S.\cite{Yeh08,Mizugushi09} The superconducting transition temperature $T_{\rm c}$ is lower than in most of the other Fe-based superconductors, reaching a value of $\simeq 14\,$K in Fe$_{y}$Se$_{1-x}$Te$_{x}$ at optimal Se-Te ratio and $\simeq 36\,$K at high pressures.\cite{Hsu08,Khasanov09_FeSeTe,Margadonna09,Sales09} Despite its relatively low $T_{\rm c}$, the binary Fe$Ch$ system is attractive for fundamental investigations of the interplay between magnetism and superconductivity because of (i) its simple crystallographic structure, (ii) the relative ease with which single crystals can be grown, and (iii) the similarity of the Fermi-surface topology with that of other Fe-based superconductors.\cite{Subedi08}

Recently, it has been reported that the superconducting and magnetic properties of Fe$_{y}$Se$_{x}$Te$_{1-x}$ not only depend on the Se-Te ratio but also strongly on the Fe content.\cite{Viennois10,Zhang09,Liu09,Wen09} Here we report a systematic investigation of the magnetic and superconducting properties, and their interplay, of Fe$_{y}$Se$_{0.25}$Te$_{0.75}$ with different nominal Fe content in the range $0.9\leq y \leq1.1$. Samples with low Fe content ($y<1$) are found to be bulk superconductors with coexistent magnetic order that sets in at a temperature below $T_{\rm c}$. Stoichiometric samples ($y\sim1$) show filamentary superconductivity and magnetic order. Fe-rich samples ($y>1$) are almost purely magnetic with only traces of superconductivity. Interestingly, the magnetic order was observed to be incommensurate\cite{Khasanov09_FeSeTe,Babkevich10} throughout the entire range of nominal Fe content investigated, although the correlation length increases with increasing Fe content.

The Fe$_y$Se$_x$Te$_{1-x}$ samples were prepared within a wide range of nominal Fe content from $y=0.9$ to 1.1 ($y=0.90$, 0.95, 0.98, 1.00, 1.01, 1.02, 1.03, 1.07, 1.10), with a fixed Se to Te ratio of $x=0.25$. The single crystals were prepared in the form of rods with masses of $m\sim4-5\,$g by a modified Bridgman method. For a detailed description of the procedure see Ref. \onlinecite{Bendele10_FeSeTe}.

To establish the crystal structure and the stoichiometry, the samples were investigated by single-crystal X-ray diffraction (XRD) at room temperature. Data reduction and numerical absorption correction were performed using the Bruker AXS Inc software package.\cite{XRD1} The crystal structure was solved by direct method and refined on $F^2$, employing the programs SHELXS-97 and SHELXL-97.\cite{XRD2} All crystals reveal a tetragonal lattice (space group $P4/nmm$) with the lattice parameters $a$ and $c$ presented in Table \ref{table_xrd}. The refined Se/Te ratio Se$_{\rm occ}$ is within standard deviation ($\pm 3$\,\%) close to the nominal content. The three samples with the nominal Fe content of 0.95, 0.98, and 1.03 reveal 100\,\% ($\pm 2$\,\%) Fe occupation Fe$_{\rm occ}$ at the $2b$ site [($1/4,3/4,1/2$) for the space-group $P4/nmm$, origin choice 2]. The sample with $y=1.07$ shows an occupation of 100\,\% Fe at the $2b$ position and some sharp maximum (2.7\,\AA\ from the Se/Te atom) on a difference Fourier map $F_0 - F_{\rm c}$ indicating that the remaining Fe (8\,\%) partially occupies the interstitial $2c$ site ($3/4,3/4,0.2105(2)$).

\begin{table}[t!]\centering\caption{Summary of structural parameters obtained from single-crystal XRD of selected compositions of Fe$_y$Se$_{0.25}$Te$_{0.75}$. $h_{\rm Se/Te}$ denotes the height of $Ch$ above the Fe plane.}\hspace{20mm}
\begin{tabular}{l|c|c|c|c}\toprule
 &$y=0.95$ & $y=0.98$ & $y=1.03$ & $y=1.07$\\\hline
$a$-axis (\AA)	& 3.8125(5)	& 3.8096(3)	& 3.8090(4)	& 3.8104(3)	\\
$c$-axis (\AA)	& 6.1699(13)	& 6.1524(6)	& 6.1562(8)	& 6.1717(10)	\\
$h_{\rm Se/Te}$(\AA)& 1.717(2)	& 1.709(1)	& 1.715(1)	& 1.733(1)	\\
Fe$_{\rm occ}$	& 1.00(2)	& 1.00(2)	& 1.00(2)	& $1.00+0.08$\footnote{
1.00 at the $2b$ and 0.08 at the $2c$ site}	\\
Se$_{\rm occ}$	& 0.25(3)	& 0.26(3)	& 0.24(3)	& 0.23(3)	\\
\toprule
\end{tabular}
\label{table_xrd}
\end{table}
The superconducting properties of Fe$_{y}$Se$_{0.25}$Te$_{0.75}$ were studied on plate-like samples with a typical mass $\sim 50\,$mg that were always mounted with the flat surface ($ab$-plane) parallel to the magnetic field to minimize the demagnetization effect on the magnetic moment. The Zero-field cooled (ZFC) susceptibility measurements were performed with a Quantum Design 7\,Tesla Magnetic Property Measurement System (MPMS-XL7) SQUID Magnetometer in a magnetic field of $\mu_0 H=0.3\,$mT using the Reciprocating Sample Option (RSO). The data are shown in Fig.~\ref{fig_magn_musr}(a) for samples with representative doping. The measurements indicate that Fe$_{y}$Se$_{0.25}$Te$_{0.75}$ starting from a nominal Fe content of $y=0.90$ up to $y=0.98$ exhibits bulk superconductivity, since $\chi_{\rm DC}(2\, \rm{K})\simeq -1$ close to ideal diamagnetism. However, only the samples with the lowest Fe content ($y=0.90$ and $0.95$) show a saturation of the magnetic moment to $\chi\simeq-1$ expected for a Meissner state. Already for $y=0.98$ the transition is broad and tends to saturate only below $2\,$K. The onset of the superconducting transition decreases with increasing nominal Fe content: $T_{\rm c}^{\rm onset}\simeq11.5\,$K, $\simeq10.6\,$K and $\simeq8\,$K for $y=0.90$, 0.95 and 0.98. The samples with a nominal Fe content higher than $y\simeq1$ show only traces of superconductivity with a superconducting volume fraction of less than $\approx 30\,$\% at low temperatures for $y\sim 1.00$ and less than $\approx 5\,$\% for $y\geq 1.03$. The onset of the superconducting transition of the samples showing traces of supercondutivity is always at $T_{\rm c}^{\rm onset}\approx 9$~K. In order to distinguish samples with bulk superconductivity, $T_{\rm c}$ was defined as the midpoint of the superconducting transition (namely $\chi_{\rm DC}=-0.5$). Thus $T_{\rm c}\simeq9.7$, 8.5, and 6.5\,K for the compositions $y=0.90$, 0.95, and 0.98.

The magnetic response of the samples was investigated by Zero-field (ZF), transverse field (TF) and longitudinal field (LF) muon spin rotation ($\mu$SR) experiments carried out at the $\pi$M3 beam line at S$\mu$S at the Paul Scherrer Institute (PSI), Switzerland. In TF geometry the muons stopping in magnetic parts of the samples lose their polarization relatively fast, because the field at the muon stopping site is a superposition of the internal field and the applied external field of $11.8\,$mT. The internal field distribution is examined in ZF measurements, whereas LF experiments provide information whether the internal field is static or dynamic (fluctuating).

The ZF time spectra (not shown) of the bulk superconducting samples ($y=0.90$ and 0.95) show no difference between $T=20$\,K and $7\,$K. However, at lower temperatures an additional fast drop of the muon spin polarization $P(t)$ develops, and the $\mu$SR time spectra at $1.8\,$K and $20\,$K do not coincide any more. This suggests that magnetic ordering observed by $\mu$SR develops in the sample at temperatures below $T_c$. The data were described using the function:
\begin{equation}
 P^{\rm ZF}(t)=P^{\rm ZF}_{\rm fast}(0)e^{-\Lambda^{\rm ZF}_{\rm fast}t}+P^{\rm ZF}_{\rm slow}(0)e^{-\Lambda^{\rm ZF}_{\rm slow}t}
 \label{eq_ZF}
\end{equation}
Here, $P^{\rm ZF}_{\rm fast (slow)}$ and $\Lambda^{\rm ZF}_{\rm fast (slow)}$ are the initial ZF muon spin polarization and the exponential depolarization rate of the fast (slow) relaxing component, respectively. Upon increasing the nominal Fe content to $y\leq 0.98$ the change in the relaxation occurs at temperatures above $T_{\rm c}$.
At higher temperatures ($T\sim 130\,$K) an additional change of the ZF depolarization was observed, whose origin requires further investigation.

\begin{figure}[t!]
\centering
\vspace{-0cm}
\includegraphics[width=1\linewidth]{fig_magn_musr_v3.eps}
\vspace{-1cm}
\caption{(color online) (a) Temperature dependence of the volume susceptibility $\chi_{\rm DC}$ of representative compositions ($y=0.95$, 1.00, and 1.10) of single-crystal Fe$_y$Se$_{0.25}$Te$_{0.75}$. The onset of the superconducting transition $T_{\rm c}^{\rm onset}$ and the midpoint corresponding to $\chi_{\rm DC}=-0.5$ are indicated. (b) Temperature dependence of the muon-spin polarization of the slow relaxing component ($P_{\rm slow}^{TF}$). The magnetic transition $T_N$ is determined from a fit to a Fermi-type function. }
\label{fig_magn_musr}
\end{figure}

The magnetic ordering temperature was investigated by means of TF $\mu$SR [see Fig.~1(b)]. The TF time spectra (not shown) can be divided into a fast drop of the muon spin polarization within the first 100\,ns and a slow relaxing part below the temperature where magnetism starts to develop. Accordingly, the signal was divided into two parts:
\begin{eqnarray}
 \nonumber P^{\rm TF}(t)=P_{\rm fast}^{\rm TF}(0)\exp[-\Lambda^{\rm TF}_{\rm fast} t]\cos(\gamma_\mu B t+\phi)\\
         + P_{\rm slow}^{\rm TF}(0)\exp[-\Lambda^{\rm TF}_{\rm slow} t]\cos(\gamma_\mu B t+\phi)
 \label{eq_TF},
\end{eqnarray}
where $\gamma_\mu/2\pi=135.5\,$MHz/T is the muon gyromagnetic ratio, $\phi$ is the initial phase of the muon ensemble, and $\Lambda_{\rm fast (slow)}$ the exponential relaxation rate. The fast relaxing component $P_{\rm fast}$, attributed to the development of magnetism, represents the magnetic volume fraction of the sample and increases with decreasing temperature [Fig.~1(b)]. In the bulk superconducting samples ($y=0.90$ and 0.95) it occupies at the lowest investigated temperature ($T=1.6\,$K) more than $60\,$\% of the signal, indicating that $\geq 60\,$\% of the sample is magnetically ordered. Furthermore, the relaxation of the slow relaxing part increases just below $T_{\rm c}$ indicating the formation of a vortex lattice in the superconducting state. Upon increasing the nominal Fe content to $y\leq 0.98$ the samples were found to be $100\,$\% magnetic at $1.6\,$K as the muon spin polarization drops to zero within the first $100\,$ns.

LF measurements reveal that the magnetic order is static in the bulk superconducting samples since the muon spin polarization recovers almost $100\,$\% at $B^{\rm LF}=0.64\,$T and the muon spins decouple from the static internal fields $B_{\rm int}$ at $B^{\rm LF}\simeq 10\cdot B_{\rm int}$.\cite{Schenk86} Thus the static internal field in the superconducting samples is $B_{\rm int}\gtrsim 0.1\,$T at the muon stopping site.

The ordering temperature $T_{\rm N}$ in Fig.~1(b) was determined by fitting a Fermi-type function $f(T)=[1+\exp( \beta(T_{\rm N}-T))]^{-1}$ ($\beta^{-1}$ is the with of transition) to the data [solid lines in Fig.~1(b)].\cite{Khasanov08_isotope} It develops with increasing nominal Fe content from $T_{\rm N}=1.7\,$K at the lowest Fe content ($y=0.90$) to 2.7\,K, 5\,K, 10\,K, 12.5\,K, and 17\,K for $y=0.95$, 0.98, 1.00, 1.01, and 1.02 respectively. Above $y=1.03$ it seems to saturate at $T_{\rm N}=27$\,K (see \onlinecite{Khasanov09_FeSeTe}), 28\,K and 29\,K for $y=1.03$, 1.07, and 1.10.

\begin{figure}[t!]
\centering
\vspace{-0cm}
\includegraphics[width=1\linewidth]{fig_neutron_v4.eps}
\caption{(color online) Neutron diffraction from Fe$_{y}$Se$_{0.25}$Te$_{0.75}$ samples for $y =  0.95$, 1.00 and 1.10. (a) Intensity profiles along the $(h,0,0.5)$ direction at $T = 2$\,K  after subtraction of a background signal measured at $T\approx$50\,K. The solid line shows a fit to a Lorentzian function convoluted with the resolution function of the instrument. (b) Temperature dependence of the intensity of the magnetic reflection. In both panels the $y=1$ and $y=1.1$ data are displaced vertically for clarity.}
\label{fig_neutron}
\end{figure}

To investigate the change in the magnetic correlations with Fe content, we performed neutron diffraction measurements on single crystals of Fe$_{y}$Se$_{0.25}$Te$_{0.75}$ for $y = 0.95$, 1.00, and 1.10 at 1.5\,K. The experiments were performed on the triple-axis spectrometer TASP\cite{TASP} at the SINQ spallation source (PSI).\cite{SINQ} The instrument was configured for a neutron wave vector of 2.66\,\AA$^{-1}$ with no collimation. Each sample was aligned on nuclear Bragg reflections to an accuracy of better than 0.008\,r.l.u. at 1.5\,K.
Figure~\ref{fig_neutron}(a) shows elastic-scans in the $(h,0,l)$ scattering plane. A background recorded at around 50\,K (above the magnetic ordering temperature) was subtracted in order to isolate the magnetic contribution to the scattering at 1.5\,K. Diffuse incommensurate magnetic peaks centered at ${\bf q} = (1/2 - \delta,0,1/2)$, with $\delta \approx 0.03$, are observed in all three samples. These results together with our previous neutron diffraction studies show that samples in the entire range $0.95<y<1.10$ have incommensurate magnetic order.\cite{Khasanov09_FeSeTe,Babkevich10} We note that our results hint at a possible reduction in incommensurability $\delta$ for the $y=1.00$ sample. However, due to the broad nature of the peak this shift may not be related to the sample. As can be seen from Fig.~\ref{fig_neutron}(a), the magnetic peaks along $(h,0,0.5)$ appear to become broader with reduced Fe-content with correlation lengths along $a$ deduced to be 7.1(5)\,\AA, 8.4(6)\,\AA, and 13.8(8)\,\AA\ for $y = 0.95$, 1.00, and 1.10 respectively. This would suggest that the magnetic correlations become more short-ranged on lowering the Fe content.

\begin{figure}[b!]
\centering
\vspace{-0cm}
\includegraphics[width=1\linewidth]{fig_phase_diagram_v10.eps}
\caption{(color online) (a) Variation of $T_{\rm c}$, $T_{\rm N}$ and the Se/Te height $h_{\rm Se/Te}$ in Fe$_y$Se$_{0.25}$Te$_{0.75}$ as a function of $y$. (b) Three dimensional phase diagram of $T_{\rm c}$ and $T_{\rm N}$ of Fe$_y$Se$_{x}$Te$_{1-x}$ as a function of $x$ and $y$. M and SC denote the magnetic and superconducting phases, respectively. The circles are from this work, and the squares are data taken from Khasanov {\it et al.}\cite{Khasanov09_FeSeTe}}
\label{fig_Phase_Diagram}
\end{figure}

The temperature at which magnetic order sets in appears also to be dependent on $y$ in Fe$_{y}$Se$_{0.25}$Te$_{0.75}$ as shown in Fig.~\ref{fig_neutron}(b). In the $y=1.10$ sample magnetic order is found to develop below $\approx 50$\,K. However, for $y=1.00$ and $y=0.95$ the transition temperature is reduced to around 30\,K. These results are in good agreement with the change in $T_{\rm N}$ with doping observed by $\mu$SR, but due to the difference in fluctuation rates sampled by neutrons and muons, the temperatures at which spin freezing occurs are not the same.

The results of the magnetization and $\mu$SR measurements of Fe$_{y}$Se$_{0.25}$Te$_{0.75}$ are summarized in a phase diagram in Fig.~\ref{fig_Phase_Diagram}(a).
In the region of $y <1$ a coexistence of superconductivity and magnetism is observed. A recent nuclear magnetic resonance (NMR) study on BaFe$_{2-x}$Co$_x$As$_2$ showed the appearance of magnetic order on all Fe sites and ruled out nanoscale segregation in this material.\cite{Laplace09} More reasonable is a coexistence of the two order parameters on atomic scale, as was already suggested earlier \cite{Khasanov09_FeSeTe,Bendele10,Laplace09} and in agreement with models recently proposed.\cite{Vorontsov09,Cvetkovic09}

It has been predicted theoretically that the height of the chalcogenide above the Fe-plane $h_{\rm Se/Te}$ can influence the magnetic ordering in Fe$_y$Se$_{x}$Te$_{1-x}$.\cite{Moon10} For that reason $h_{\rm Se/Te}$ derived from XRD measurements is plotted in addition to $T_{\rm c}$ and $T_{\rm N}$ in Fig.~\ref{fig_Phase_Diagram}(a). The crossover from the purely magnetic to the coexistence region occurs at $h_{\rm Se/Te}\simeq1.72\,$\AA. This result is in good agreement with Density-Functional Theory (DFT) calculations, which predict a change from double- to single-stripe AFM ordering at $h_{\rm Se/Te}\sim 1.71-1.72\,$\AA.\cite{Moon10} However, the theoretical calculations are based on Se substitution only and do not take into account any excess Fe. Furthermore, $h_{\rm Se/Te}$ might not be the only factor since bulk superconductivity appears only for $y\leq0.98$ whereas the proposed critical height $h_{\rm Se/Te}=1.72$\,\AA\ is reached already at $y=1.03$.

The amount of excess Fe seems to play a major role in this system as only for $y\leq1$ does bulk superconductivity occur. Here the magnetic correlations in the system become more short-ranged and lead to less well correlated magnetic order as compared with $y>1$. Thus, only upon reducing the magnetic correlations by lowering the amount of Fe in the system does it become superconducting. Nevertheless, it seems unlikely that the excess Fe acts as isolated magnetic moments that destroy superconductivity. It might on the other hand act as a magnetic electron donor\cite{Zhang09} that suppresses superconductivity and induces weakly localized electronic states.\cite{Wen09,Liu09}

The tentative three-dimensional phase diagram of the transition temperatures $T_{\rm c}$ and $T_{\rm N}$ of Fe$_y$Se$_{x}$Te$_{1-x}$ for $0\leq x\leq 0.5$ and $0.9\leq y\leq 1.1$ is shown in Fig. \ref{fig_Phase_Diagram}(b). Fe$_y$Te is always antiferromagnetically ordered.\cite{Bao09,Li09} Upon substituting Te by Se the order becomes weaker while superconductivity is enhanced, and finally the system becomes bulk superconducting. This behavior can be tuned not only by the substitution of Se but also by adjusting the Fe content. The superconductivity is suggested to be of multi-band nature, where different doping channels might be involved.\cite{Bendele10_FeSeTe}

To conclude, we have found that the phase diagram of Fe$_y$Se$_{0.25}$Te$_{0.75}$ in the range $0.9\leq y \leq 1.1$ exhibits a strong dependence of its superconducting and magnetic phases on $y$. In the low Fe content region $y\leq 1$ bulk superconductivity and incommensurate magnetism coexist. With increasing $y$ the magnetic order becomes correlated over a longer range and superconductivity vanishes.  This work emphasizes that not only the Se/Te ratio, but also the Fe content is important in controlling the magnetic and superconducting properties of the iron chalcogenides.

This work was performed at the S$\mu$S and SINQ, Paul Scherrer Institute (PSI, Switzerland). It was partially supported by the Swiss National Foundation and the NCCR program MaNEP.