Extra-Mixing in Luminous Cool Red Giants. Hints from Evolved Stars with and without Li

Abstract: We present an analysis of Li abundances in low mass stars (LMS) during the Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stages, based on a new determination of their luminosities and evolutionary status. By applying recently suggested models for extra-mixing, induced by magnetic buoyancy, we show that both Li-rich and Li-poor stars can be accounted for. The simplest scenario implies the development of fast instabilities on the RGB, where Li is produced. When the fields increase in strength, buoyancy slows down and Li is destroyed. 3He is consumed, at variable rates. The process continues on the AGB, where however moderate mass circulation rates have little effect on Li due to the short time available. O-rich and C-rich stars show different histories of Li production/destruction, possibly indicative of different masses. More complex transport schemes are allowed by magnetic buoyancy, with larger effects on Li, but most normal LMS seem to show only the range of Li variation discussed here.

\section{Introduction}
Most K-type giants are depleted in lithium. For low-metallicity
objects an upper limit on the Li content was set by
\citet{gratton}, stating that field, metal-poor (hence low-mass)
stars on the upper RGB have $\log~ \epsilon$(Li) $\le$ 0

Mixing processes on the Main Sequence and up to the first
dredge-up imply that Li be gradually carried down from the
photosphere to regions of high temperature and destroyed, while
any amount of $^7$Be produced by H burning burns "on the flight"
and new Li is not produced. The observations of Li-poor red giants
are therefore well understood. It was further established three
decades ago \citep{sweigart} and supported by many observations
\citep[e.g.][]{gilroy,gilroybrown,kraft} that, after the first
dredge-up, low mass stars must experience other mixing episodes,
sufficiently slow to gradually carry to the surface a considerable
amount of $^{13}$C, so that the $^{12}$C/$^{13}$C ratio decreases
from 25-30 at the first dredge-up, down to 12-15 (in population I
stars) or to 4-8 (in population II red giants). The occurrence of
these mixing phenomena is made easy by the fact that, after
dredge-up, the H burning shell advances in a homogeneous region,
so that the natural barrier opposed to mass transport by a
chemical stratification is not present. The onset of this phase is
accompanied by a local decrease in the luminosity, followed by a
new rise, so that the same regions in the H-R diagram are crossed
repeatedly and the points representing coeval stars pile up,
producing a bump in the Luminosity Function. This is the so-called
"L-bump phase".

Any slow mixing process occurring at or after the L-bump would
further reduce not only Li, but also $^9$Be and $^3$He
\citep{bs99}, as confirmed by observations \citep{cast}. However,
a few red giants (about 2\%) show enhanced Li abundances, at
levels sometimes higher than in the present Interstellar Medium
($\log~\epsilon$(Li) $\simeq 3.3$). Various explanations have been
attempted, the most common assumption being that some form of fast
extra-mixing might transport $^{7}$Be from above the H-burning
shell to the envelope at a speed sufficient to overcome the rate
of p-captures; this is the so-called Cameron-Fowler mechanism
\citep{camfow}. Further production of $^7$Be along the path can
also occur \citep{bs99}. Very often, the cause of such
extra-mixing was looked after in the stellar rotation, inducing
shear instabilities, meridional circulation and diffusion
phenomena \citep{charb94,charb98,denvan,palacios03}. These
possibilities were frustrated by the understanding that the
stellar structure reacts to rotational distortions too quickly to
induce significant mixing on long time scales \citep{palacios06}.

Recently, two old ideas on mass transport have been revisited,
looking for alternative causes of deep circulation in stars. Both
attribute mass transport to the presence, in deep layers, of
materials lighter than in the envelope. The first mechanism ({\it
thermohaline} diffusion) attributes this to the activation of the
reaction $^3$He+$^3$He $->$ $^4$He + 2p \citep{eggleton}, reducing
the molecular weight. The second process is the buoyancy of
magnetized H-burning ashes \citep{busso07a,wasserburg,nordhaus},
based on the fact that magnetic bubbles are lighter than the
environment, due to the unbalance generated by magnetic pressure
($B^2/8\pi$). This requires a magnetic dynamo to occur in red
giants, as indeed observed \citep{and88}.

Triggering these mechanisms requires specific conditions: the
first one needs a suitable molecular weight inversion; the second
is possible if strong magnetic fields exist. The two processes
might be complementary and even occur together \citep{den09}.
However, while thermohaline mixing implies a diffusive, slow
transport, magnetic buoyancy can occur at different speeds,
depending on the amount of heat exchange with the environment
\citep{den09}. Therefore, a possibility of understanding which
mechanism is at play might be provided by observations of nuclei,
like Li, for which the production and destruction are sensitive to
the mixing velocity.

Such a test requires an accurate calibration of the luminosity of
Li-bearing stars, in order to attribute them to the proper
evolutionary stage. Therefore, we here apply bolometric
corrections for evolved stars, as derived in our recent work, to
infer the absolute HR diagram of red giants showing Li in their
spectrum. This helps in setting constraints on where extra-mixing
is active. In Section 2 we present a sample of cool red giants
showing highly dispersed Li abundances, ranging from very Li-rich
to very Li-poor, and we determine for them absolute magnitudes
with the help of our bolometric corrections and of recent distance
estimates. In Section 3 we compare the inferred magnitudes with
the Li content, as taken from the literature and we present
results of nucleosynthesis and mixing calculations, using the
transport velocities suggested by magnetic buoyancy. Finally, in
Section 4 preliminary conclusions are illustrated.

\section{Bolometric Corrections and Li-abundances in Red Giants}

Any growth in our understanding of AGB stars is today bound to
improvements on two crucial parameters: luminosities and mass loss
rates. For both, determinations are hampered by the large distance
uncertainties; both are also affected by the fact that most of the
radiation emitted by these sources can be observed only at
infrared wavelengths and therefore depends on space-borne
facilities, or maybe on the exploitation of special ground-based
locations, like Antarctica \citep{irait}.

In the last years we performed an extensive analysis of the
photometric properties of AGB stars, using available infrared data
from the ISO and MSX experiments and looking for relations between
their luminosities and their main chemical and physical
parameters. The aim of this work was to put constraints on basic
stellar parameters, thus offering calibrating tools to
evolutionary models. In this framework, large samples of C-rich,
S-type, M-type AGB stars were collected and studied. The first
results were published for carbon-rich and S-type stars in
\citet{guabus,busso07b,gua06}. In Figures 1 and 2 we present the
bolometric corrections we derived for samples of galactic C-rich
and O-rich (S-type) stars, respectively. The techniques adopted
and the source properties are described in \citet{gua06} and in
\citet{guabus}.

For the purposes of the present work we use the bolometric
corrections of Figures 1 and 2, applying them to a sample of AGB
stars showing highly dispersed Li abundances (from Li-rich to
Li-poor). For these sources we determine the apparent bolometric
magnitude and we use the recent revision of the Hipparcos
catalogue \citep{hipparcos}, to infer the absolute magnitude. The
AGB stars considered, their magnitudes and their Li abundances are
shown in Table 1. In Table 2 we present a second sample, made of
RGB stars. The bolometric magnitudes are taken from the
literature, with distances updated by \cite{hipparcos}. Some
explanations can be useful in presenting the sources in Tables 1
and 2.

\begin{figure}[t*]
\begin{center}
\includegraphics[scale=0.4, angle=0]{fig1.eps}
\caption{The bolometric correction for C stars, as a function of the
infrared color K-[12.5].}\label{fig1}
\end{center}
\end{figure}

\begin{figure}[t*]
\begin{center}
\includegraphics[scale=0.4, angle=0]{fig2.eps}
\caption{The bolometric correction for O-rich stars, derived from
MS- and S-type red giants, as a function of the infrared color
K-[12.5].}\label{fig2}
\end{center}
\end{figure}
A) Horizontal lines in the tables divide the stars of the sample
in a few sub-groups, according to their spectral classification;
from top to bottom we find:

\onecolumn
\begin{table}[t*]
\begin{center}
\caption{AGB sample}\label{table1}
\begin{tabular}{llcrrcl}
\hline Source Name & Sp. Type & D (kpc)  & m$_{bol}$ & M$_{bol}$ & $\delta$M$_{bol}$ & Log $\varepsilon$(Li)  \\
\hline\hline
WZ Cas  &   C9,2JLi(N1P)    &   0.88    &   3.16    &   $-$6.56   &   1.26    &   +5.80$^a$$^,$$^b$ \\
Y CVn   &   C5,4J(N3)   &   0.32    &   1.73    &   $-$5.80   &   0.49    &   $-$0.40$^a$$^,$$^b$    \\
RY Dra  &   C4,5J(N4P)  &   0.43    &   3.10    &   $-$5.07   &   0.80    &   $-$0.45$^a$$^,$$^b$   \\
T Lyr   &   C$-$J4:p C2 5 j3.5    &   0.72    &   2.86    &   $-$6.43   &   1.02    &   +2.20$^a$ \\
VX And  &   C$-$J4.5 C2 5 j5 MS5  &   0.39    &   3.88    &   $-$4.09   &   0.93    &   +2.60$^c$ \\
V614 Mon    &   C4,5J(R5)   &   0.47    &   4.43    &   $-$3.95   &   1.23    &   +1.30$^c$ \\
\hline
Z Psc   &   C7,2(N0)    &   0.38    &   3.53    &   $-$4.40   &   0.78    &   $-$1.30$^a$$^,$$^b$    \\
R Scl   &   C6,5EA(NP)  &   0.27    &   2.55    &   $-$4.61   &   0.68    &   $-$0.40$^a$$^,$$^b$    \\
U Cam   &   C3,9$-$C6,4E(N5)  &   0.97    &   3.25    &   $-$6.69   &   1.50    &   $-$0.75$^a$$^,$$^b$   \\
ST Cam  &   C5,4(N5)    &   0.99    &   3.21    &   $-$6.76   &   1.65    &   $-$0.50$^a$$^,$$^b$    \\
R Lep   &   C7,6E(N6E)  &   0.41    &   2.88    &   $-$5.20   &   1.17    &   $-$1.10$^a$$^,$$^b$    \\
W Ori   &   C5,4(N5)    &   0.38    &   2.19    &   $-$5.70   &   1.03    &   $-$1.05$^a$$^,$$^b$   \\
Y Tau   &   C6.5,4E(N3) &   0.36    &   3.07    &   $-$4.72   &   1.03    &   $-$0.80$^a$$^,$$^b$    \\
TU Gem  &   C6,4(N3)    &   0.52    &   3.55    &   $-$5.03   &   0.68    &   $-$0.55$^a$$^,$$^b$   \\
BL Ori  &   C6,3(NB,TC) &   0.60    &   3.37    &   $-$5.52   &   0.68    &   $-$1.10$^a$$^,$$^b$    \\
W CMa   &   C6,3(N) &   0.79    &   3.54    &   $-$5.95   &   0.68    &   $-$1.20$^a$$^,$$^b$    \\
X Cnc   &   C5,4(N3)    &   0.34    &   2.71    &   $-$4.94   &   0.83    &   $-$1.00$^a$$^,$$^b$  \\
Y Hya   &   C5,4(N3P)   &   0.39    &   3.19    &   $-$4.77   &   0.63    &   $-$0.90$^a$$^,$$^b$    \\
VY Uma  &   C6,3(N0)    &   0.38    &   3.15    &   $-$4.75   &   0.54    &   $-$0.50$^a$    \\
S Sct   &   C6,4(N3)    &   0.39    &   3.30    &   $-$4.66   &   0.73    &   $-$0.50$^a$$^,$$^b$    \\
V Aql   &   C5,4$-$C6,4(N6)   &   0.36    &   2.55    &   $-$5.23   &   0.79    &   $-$0.65$^a$$^,$$^b$   \\
UX Dra  &   C7,3(N0)    &   0.39    &   3.03    &   $-$4.90   &   0.49    &   +0.10$^a$$^,$$^b$     \\
AQ Sgr  &   C7,4(N3)    &   0.33    &   3.51    &   $-$4.11   &   0.74    &   $-$0.95$^a$$^,$$^b$   \\
V460 Cyg    &   C6,4(N1)    &   0.62    &   2.89    &   $-$6.08   &   0.77    &   $-$0.90$^a$$^,$$^b$    \\
TX Psc  &   C7,2(N0)(TC)    &   0.28    &   2.02    &   $-$5.21   &   0.48    &   $-$0.65$^a$$^,$$^b$   \\
CR Gem  &   C8,3E(N)    &   0.35    &   4.28 &   $-$3.44   &   1.78    &   +0.00$^a$   \\
RS Cyg  &   C8,2E(N0PE) &   0.47    &   3.99 &   $-$4.39   &   1.17    &   $-$0.50$^a$    \\
RV Mon  &   C4,4$-$C6,2(NB/R9)    &   0.52    &   4.56  &   $-$4.02   &   0.99    &   $-$1.00$^a$  \\
RY Mon  &   C5,5$-$C7,4:(N5/R)    &   0.43    &   3.66 &   $-$4.52   &   0.97    &   +0.50$^a$     \\
S Cep   &   C7,4E(N8E)  &   0.41    &   2.70  &   $-$5.35   &   0.82    &   +0.00$^a$   \\
TT Cyg  &   C5,4E(N3E)  &   0.56    &   4.57  &   $-$4.18   &   0.87    &   +0.70$^a$     \\
TU Tau  &   C$-$N4.5 C2 6 &   0.43    &   4.30 &   $-$3.89   &   1.35    &   $-$0.50$^a$    \\
V Cyg   &   C5,3E$-$C7,4E(NPE)    &   0.37    &   2.74 &   $-$5.08   &   1.51    &   $-$1.00$^a$  \\
\hline
Y Lyn   &   M6SIB$-$II    &   0.25    &   1.66    &   $-$5.33   &   0.77    &   $-$2.00$^d$  \\
RS Cnc  &   M6EIB$-$II(S) &   0.14    &   0.52    &   $-$5.21   &   0.41    &   $-$2.00$^d$  \\
V1981 Cyg   &   S4/1III &   0.30    &   3.43    &   $-$3.96   &   0.46    &   $-$2.00$^d$  \\
\hline
bet And &   M0+ IIIa    &   0.06    &   0.17    &   $-$3.74   &   0.32    &   $-$0.48$^e$   \\
CU Dra  &   M3III   &   0.11    &   2.09    &   $-$3.19   &   0.30    &   $-$1.50$^e$    \\
BY Boo  &   M4$-$4.5III   &   0.16    &   2.15    &   $-$3.84   &   0.34    &   $-$0.87$^e$   \\
OP Her  &   M5IIB$-$IIIA(S)   &   0.30    &   2.47    &   $-$4.92   &   0.42    &   $-$2.00$^d$  \\
\hline
V2652 Sgr   &  M9     &   $\simeq$ 8.0   &   9.98    &   $-$4.52   &  0.30     &   +2.00$^f$   \\
V3252 Sgr   &  M7     &  $\simeq$ 8.0   &    9.74   &   $-$4.76   &   0.30    &   +1.10$^f$ \\
V3537 Sgr   &  M9S     &  $\simeq$ 8.0   &   9.22    &   $-$5.28   &  0.30     &   +0.80$^f$ \\
V2017 Sgr   &  M7S     &  $\simeq$ 8.0   &   9.07    &   $-$5.43   &  0.30     &   +0.80$^f$ \\
\hline\hline
\end{tabular}
\medskip\\
For Log $\varepsilon$(Li): $^a$\citet{boffin}, $^b$\citet{denn},
$^c$\citet{abia}, $^d$\citet{vanture},
$^e$\citet{luck},$^f$\citet{utt},
$^g$\citet{lambert},$^h$\citet{mallik},
$^i$\citet{charb00}, $^j$\citet{brown}.\\
\end{center}
\end{table}

\begin{table}[t*]
\begin{center}
\caption{RGB sample}\label{table2}
\begin{tabular}{llcrrll}
\hline Source Name & Sp. Type & D (kpc)  & m$_{bol}$ & M$_{bol}$ & T$_{eff}$ (K) & Log $\varepsilon$(Li)  \\
\hline\hline
alpha Boo   &   K1.5III &   0.01    &   $-$0.82   &   $-$1.08   &   4290    &   $-$1.50$^g$    \\
beta UMi    &   K4III   &   0.04    &   1.02    &   $-$2.00   &   4000 &   $-$1.50$^g$    \\
alpha Ser   &   K2IIIb  &   0.02    &   1.92    &   0.14    &   4300 &   $-$0.28$^g$   \\
gamma Dra   &   K5III   &   0.05    &   1.18    &   $-$2.20   &   4000 &   $-$0.80$^g$$^,$$^h$    \\
\hline
HD 787  &   K4III   &   0.19    &   4.15 &   $-$2.25   &   3990$^1$    &   +1.80$^i$ \\
HD 30834    &   K3III   &   0.18    &   3.81 &   $-$2.52   &   3920$^1$    &   +1.80$^i$ \\
HD 39853    &   K5III   &   0.20    &   4.53 &   $-$2.01   &   3920$^1$    &   +2.80$^i$ \\
HD 121710   &   K3III   &   0.18    &   4.05 &   $-$2.20   &   4100$^2$    &   +1.50$^i$ \\
HD 146850   &   K3IIICNpv   &   0.33    &   4.98 &   $-$2.59   &   4200$^3$    &   +2.00$^i$   \\
HD 9746 &   K1III:  &   0.16    &   5.51 &   $-$0.47   &   4420$^2$    &   +3.75$^i$    \\
HD 148293   &   K2III   &   0.09    &   4.71 &   $-$0.06   &   4640$^2$     &   +2.00$^i$   \\
HD 183492   &   K0III   &   0.09    &   5.10  &   0.29    &   4720$^1$    &   +2.00$^i$   \\
HD 219025   &   K2IIIp  &   0.26    &   7.03 &   $-$0.06   &       &   +3.00$^i$   \\
\hline
HD 3817 &   G8III   &   0.11    &   4.85 &   $-$0.34   &   5020$^2$    &   $<$+0.70$^j$ \\
HD 3856 &   G9III$-$IV    &   0.16    &   4.90 &   $-$1.10   &   4750$^2$    &   +1.20$^j$ \\
HD 4627 &   G8III   &   0.18    &   5.28   &   $-$0.96   &   4610$^2$    &   $<$+0.00$^j$   \\
HD 19845    &   G9III   &   0.12    &   5.33 &   0.00    &   4830$^2$    &   $<$+0.30$^j$ \\
HD 78235    &   G8III   &   0.08    &   5.08 &   0.46    &   5000$^2$    &   +0.80$^j$ \\
HD 108225   &   G9III   &   0.08    &   4.66 &   0.14    &   4960$^2$    &   $<$+0.60$^j$ \\
HD 119126   &   G9III   &   0.10    &   5.14 &   0.07    &   4750$^2$    &   $<$+0.20$^j$ \\
HD 120048   &   G9III   &   0.13    &   5.49 &   $-$0.09   &   4890$^2$    &   +1.00$^j$  \\
HD 129336   &   G8III   &   0.12    &   5.09 &   $-$0.37   &   4940$^2$    &   $<$+0.50$^j$ \\
HD 134190   &   G7.5III &   0.08    &   4.75 &   0.26    &   4880$^2$    &   $<$+0.40$^j$ \\
HD 141680   &   G8III   &   0.08    &   4.59 &   $-$0.02   &   4750$^2$    &   $<$+0.20$^j$ \\
HD 152815   &   G8III   &   0.08    &   4.90 &   0.31    &   4900$^2$    &   $<$+0.40$^j$ \\
HD 158974   &   G8III   &   0.12    &   5.14 &   $-$0.28   &   4890$^2$    &   $<$+0.40$^j$ \\
HD 163532   &   G9III   &   0.12    &   4.40 &   $-$1.04   &   4710$^2$    &   $<$+0.20$^j$ \\
HD 171391   &   G8III   &   0.10    &   4.44 &   $-$0.58   &   4980$^2$    &   +1.20$^j$ \\
HD 176598   &   G8III   &   0.10    &   5.25 &   0.34    &   4925$^2$    &   +1.10$^j$ \\
HD 186675   &   G7III   &   0.09    &   4.49 &   $-$0.24   &   4910$^2$    &   $<$+0.30$^j$ \\
HD 192944   &   G8III   &   0.13    &   4.77 &   $-$0.87   &   4900$^2$    &   $<$+0.30$^j$ \\
HD 194013   &   G8III$-$IV    &   0.08    &   4.86 &   0.46    &   4850$^2$    &   $<$+0.40$^j$ \\
HD 201381   &   G8III   &   0.05    &   4.04 &   0.60    &   4960$^2$    &   $<$+0.50$^j$ \\
HD 211391   &   G8III   &   0.06    &   3.68 &   $-$0.11   &   4900$^2$    &   $<$+0.30$^j$ \\
HD 211554   &   G8III   &   0.19    &   4.87 &   $-$1.51   &   4920$^2$    &   +0.50$^j$ \\
HD 215030   &   G9III   &   0.09    &   5.41 &   0.60    &   4800$^2$    &   +0.60$^j$ \\
\hline\hline
\end{tabular}
\medskip\\
For T$_{eff}$: $^1$\citet{melo}, $^2$\citet{brown},
$^3$\citet{charb00}.\\ For Log $\varepsilon$(Li):
$^a$\citet{boffin}, $^b$\citet{denn}, $^c$\citet{abia},
$^d$\citet{vanture}, $^e$\citet{luck},$^f$\citet{utt},
$^g$\citet{lambert},$^h$\citet{mallik},
$^i$\citet{charb00}, $^j$\citet{brown}.\\
\end{center}
\end{table}

\twocolumn

\begin{itemize}
    \item In Table 1: 1) CJ stars, 2) the C(N) sample, 3) S-type sources, 4) M-type
giants and finally 5) a few O-rich AGB sources of the Galactic
bulge from \citet{utt}.
    \item In Table 2: 1) A sample of close-by, Li-poor K-type stars, 2) a similar
    one for Li-rich K-type sources, 3) a group of G-type, Li-depleted giants.
\end{itemize}

B) As mentioned, distances (in kpc) are in general from
\cite{hipparcos}. For three sources in Table 1, for which the
Hipparcos' parallaxes are not available, we use the data from
\citet{bergeat}. For the four stars of the bulge the estimates of
the distance and therefore also the absolute bolometric magnitudes
are taken from \citet{utt}.

C) The apparent bolometric magnitude of the sour\-ces in Table 1 was
obtained through the integration of the data from ISO, 2MASS, MSX
and IRAS-LRS at different wavelengths and by applying the bolometric
corrections suitable for their chemical type. Concerning
K-type stars, due to their known distance and
luminosity, both their apparent and absolute bolometric magnitude
are taken from the available literature.

D) In Table 1 a typical uncertainty of 0.25 mag can be attributed to
apparent magnitudes; the uncertainty on the distance is given by
\citet{hipparcos}.

E) For the two sub-groups at the bottom of Table 2 (second and third
one) we calculated the magnitudes from data available in the literature
\citep{charb00,takeda}, updated with new estimates of the distance.

F) The Li abundances were taken from various references: the full lists
are shown in the footnotes to the tables. In the case of
sources with multiple references we used an average of the data available.

G) The effective temperatures presented in Table 2 come
from the list of references indicated in the footnote to the same Table.

{\begin{figure}[t!]
   \centering
   {\includegraphics[width=\columnwidth]{fig3.eps}}
\caption{The sample of cool red giants showing Li in their spectra.
Arrows indicate the suggested evolutionary path induced by
extramixing. Dashed lines roughly limit the region occupied by C(N)
giants and their range of luminosities as inferred from the
empirical Luminosity Function of Figure 4. The dash-dotted line
indicates the maximum Li abundance observed in pop. II red
giants.}\label{fig3}%
    \end{figure}
}

\section{On the Luminosities of Red Giants Showing Li}

Figure 3 presents the selected sample of evolved stars showing Li
in their spectra, after correcting their luminosities with the
criteria described in the previous section. The plot shows several
groups of stars. As compared to the ISM Li abundances, many of the
stars in our sample are strongly depleted in Li
($\log~\varepsilon$(Li) $\le 1$). Two main outlying groups exist
(including mainly K giants) with Li abundances near
$\log~\varepsilon$(Li) $\simeq 2$. A few super-Li rich stars are
also present, with $\log~\varepsilon$(Li) $\ge 3$. We shall not
discuss  super-Li rich stars in this paper, as they are often
peculiar objects (like CJ giants). Their high Li abundance can be
produced either by a relatively fast extra-mixing in LMS or by Hot
Bottom Burning in more massive AGB stars.

The group of stars at the left of the plot, with
$\log~\varepsilon$(Li) $\simeq 2$, is made of K giants with
luminosities typical of the L-bump on the RGB. In agreement with
\citet{charb00} we consider them as being stars newly enriched in
Li. After the extensive depletion of Li in the Main Sequence and
in the early RGB phases, their envelopes must have seen some form
of rapid mixing implementing a Cameron-Fowler mechanism, thus
bringing to the surface fresh $^7$Be synthesized above the
H-burning shell.

The second moderately Li-rich group was interpreted by
\citet{charb00} as being formed by early-AGB stars, producing Li
(after some destruction in the late RGB phases). Inspection of
stellar models, however, reveals that the early-AGB stages
suitable to reproduce the luminosities and temperatures of these K
giants fall in a temporal phase where the H shell is extinguished,
so that no $^7$Be survives to be carried to the surface. We must
therefore tentatively conclude that a new production of Li on the
early-AGB at temperatures and luminosities compatible with
observations does not occur and therefore the observed Li must be
a relic of the previous production on the RGB. We actually
interpret this group of stars as being in the upper part of the
RGB itself (see next section).

On the right side of the plots we find AGB stars, either from the
galactic disk or from the galactic bulge: these last are O-rich,
while for the galactic disk most (although not all) Li
observations concern C-rich stars. As an eye-guide we have marked
by dashed lines the region occupied by C(N) giants, and in Figure
4 we also report their Luminosity Function, as determined by the
works mentioned in Section 2. Vertical dashed lines on the two
figures allow one to realize that Li-bearing C(N) stars belong to
the typical luminosity intervals of normal C stars. We believe
that also their Li abundances should represent typical trends.

\section{Modeling the Li Production and Destruction}

As discussed elsewhere \citep{palm2} mixing mechanisms induced by
magnetic buoyancy offer a scenario in which both fast and slow
matter circulation can occur. A simple example of fast transport
is by magnetic instabilities, when portions of a magnetic
$\Omega-$shaped loop detach and travel the radiative layers above
the H shell and the convective envelope at a velocity close to the
Alfv\'en speed. In this case Be is produced by normal H-burning
above the H-shell, refueled by fresh $^3$He carried down by the
mixing. When fields grow above a certain limit, buoyant
instabilities are released sporadically, transporting (without
local mixing) their composition to the surface. This form of Li
production has an intrinsic limit in the fact that no burning
occurs along the path, so that the maximum Li abundance reachable
is set by the equilibrium Be concentration near the H shell (this
corresponds to about $\log \varepsilon_{Max}(Li) = 2 - 2.5$).

Alternatively Li can be produced by a slower mixing process, where
$^3$He has time to burn along the path; the mechanism must however
be still fast enough to save some of the produced Be. In general,
producing or destroying Li depends on a delicate equilibrium of at
least four different parameters: i) the velocity of mixing; ii)
its maximum penetration (hence the maximum temperature reached,
$T_P$); iii) the initial Li content; and iv) the time available
for the process.

\begin{figure}[h*]
\begin{center}
\includegraphics[scale=0.45, angle=0]{fig4.eps}
\caption{The Luminosity Function of C(N) stars (see text).
Vertical dashed lines establish a rough fiducial interval for the
luminosities, as also indicated in Figure 3.}\label{fig4}

\end{center}
\end{figure}

We adopt the formalism by \citet{palm1}, derived from the idea
presented by \citet{busso07a} for magnetically-induced mixing.
With these tools we try here to account for the observed Li
abundances in evolved LMS.

\begin{figure}[t*!!]
\begin{center}
\includegraphics[angle=0,width=7cm]{fig5.eps}
\caption{RGB stars showing Li and model results from the mixing
processes discussed in the text. Rhombs are GIII stars from
\citet{brown,takeda}. Filled squares are stars at the L-bump from
\citet{charb00}); crosses are K giants; filled squares are RGB
stars that remained Li-rich. The dotted line indicates Li
enrichment through fast magnetic instabilities. Solid lines show
the effects of slower transport at different mixing rates. This
might be induced by the buoyancy of large structures, exchanging
heat with the environment.}\label{fig5}
\end{center}
\end{figure}

\begin{figure}[t*!!]
\begin{center}
\includegraphics[angle=0,width=7cm]{fig6.eps}
\caption{AGB stars showing Li and model results from magnetic
mixing. Circles are C(N) stars (the empty ones with larger
luminosity uncertainties); asterisks show O-rich AGB objects in the
galactic bulge \citep{utt}: the most luminous of them also show
Tc in their spectra.}\label{fig6}
\end{center}
\end{figure}

We consider two different types of magnetic buoyancy, which can be
described as follows.

i) (\textbf{Model A}). Integrating the assumptions by
\citet{busso07a}, we consider a phase of magnetic field growth, at
the L-bump, where fields are assumed to be not strong enough to
promote the buoyancy of entire flux tubes; only local
instabilities can detach from regions near the H shell, traveling
at the Alfv\'en speed $B/\sqrt(4 \pi \rho)$, up to a few Km/sec.
In our exercise, the release of unstable magnetized bubbles is
fine-tuned to produce a rate of mass addition to the envelope of
10$^{-6}$ $M_{\odot}$/yr and the Alfv\'en velocity is computed adopting
field strengths from \citet{busso07a}. An equally fast downflow of
envelope material guarantees mass conservation, replenishing the
inventory of $^3$He in the radiative regions. No burning occurs in
the circulation, because the motion is very fast. The net effect,
for the production of Li, is the mixing of the "normal" $^7$Be
(produced near the H shell) into the envelope, where it will decay
to Li. Contemporarily, the $^3$He abundance at the surface
decreases.

ii) (\textbf{Model B}). As an opposite case, we consider that,
when the fields have grown sufficiently, larger structures, or
even entire flux tubes become buoyant, driving a circulation of
mass that slows down due to the gradual heat exchange between the
magnetized zones and the environment (heat exchanges grow with the
surface of the emerging structure). The speed is reduced down to
values of few cm/sec \citep{den09}. Again, $^3$He is destroyed at
various efficiencies, depending on the parameters mentioned above.
The new property of our model is that, if magnetic fields are at
the origin of the transport, then one has a justification for both
fast and slow regimes of circulation, and for both production or
destruction of Li; a single physical mechanism can explain very
different outcomes, depending on the parameters (and observations
will then provide the key information on the parameter values).
Models of extra-mixing so far presented lack this characteristic,
as diffusive processes induced by either rotational effects or
thermohaline mixing have a very low velocity.

Magnetic fields can even generate more complex schemes of mixing,
e.g. with upward and downward fluxes occupying different mass
fractions and traveling at different speeds. We shall afford a
more complete analysis, including those possibilities, in a
forthcoming paper. For the sake of illustration it can be now
sufficient to limit the discussion to the simple and rather
extreme cases of Models A and B.

As mentioned, we apply Model A at the luminosity bump, then we
switch to the slow circulation of Model B. This way of proceeding
is rather ad-hoc, although it might be qualitatively reasonable.
We do not know enough of magnetic fields in red giants to claim
that our hypotheses are really the correct ones, but they have at
least a basis in the observations, showing Li-rich stars at the
L-bump and (mostly) Li-poorer stars in more evolved stages (again
excluding super Li-rich objects). Ours is therefore an exercise,
aiming at showing how, with suitable assumptions, magnetic
buoyancy can in fact offer a framework where the puzzle of Li
production and destruction can find a solution.

With the above cautions in mind, Figures 5 and 6 present the
results of our effort. In Figure 5 we show, by a dotted line, the
effects of applying Model A over a composition typical of the
post-main sequence phase. Li is rather rapidly produced (in a few
million years) up to levels of $\log~\varepsilon$(Li) = 2 - 2.5
(depending on the depth of mixing and on the temperature
distribution). The example shown is for a 1.5 $M_{\odot}~$ model of half
solar metallicity but is rather typical of the intervals 0.01
$\le$ $Z$ $\le$ 0.02, 1.5 $\le$ $M$/$M_{\odot}~$ $\le$ 2. The solid lines
represent instead results of the application of Model B with
various (moderate) rates for mass circulation. As the Figure 5
shows, by varying this rate one can explain the whole range of Li
destruction displayed by observed red giants, and even an almost
complete Li survival, thus accounting for the few cool Li-rich
stars. This substantiates our suggestion that they are late RGB
objects where Li was not destroyed, not early-AGB stars where Li
was newly produced. The maximum temperature achieved by mixing
($T_P$) is chosen at the position where the equilibrium $^7$Be has
its maximum, near $\log T_{\rm H} - \log T_P$ = 0.3. One can
notice that, for this $T_P$ value, the destruction of Li becomes efficient
only for the higher $\dot M$ values shown. However, if $\dot M$
were further increased one might obtain an excessively fast
circulation, and Li would be produced. Li production occurs when
$\dot M$ is above a limiting value, somewhere between 10$^{-6}$
and 10$^{-5}$ M$_{\odot}$/yr, the actual value depending on $T_P$
and on the time available for the mixing to operate. The large
number of different possible outcomes thus forms a real zoo, in
which only observations can effectively constrain the parameters.

Finally, Figure 6 illustrates the effects of applying Model B,
with moderate mass circulation rates, to AGB stars of 1.5 - 2.0
$M_{\odot}$. No change would be seen by adopting the same rule for $T_P$
used for the RGB ($\log T_{\rm H} - \log T_P$ = 0.3) (the time
available is too short on the AGB). We have therefore pushed $T_P$
to its maximum possible value that still does not disturb the
stellar equilibrium ($\log T_{\rm H} - \log T_P$ = 0.1). Even so,
with moderate transport rates the changes are rather small and
would indicate that most of the Li production or destruction has
occurred previously, on the RGB. Stronger transport rates would
induce Li production also on the AGB, but they might also imply
strong difficulties in saving enough carbon to allow the stars to
become C-rich. These are general indications, not firm
conclusions, as we have seen that several parameters affect Li
abundances.

In conclusion our results can be summarized as follows:

i) Outlining the details of Li production and destruction requires
a careful analysis of the bolometric magnitudes of the observed
stars, as discussed in Section 2. ii) Magnetic buoyancy might
offer a framework for interpreting both moderately Li-rich and
Li-poor stars. iii) Extensive Li processing is generally
accompanied by $^3$He destruction, which fact should be relevant
for reconciling the inventory of this nucleus with Big Bang and
stellar nucleosynthesis \citep{bs99}. iv) When the $^3$He
destruction is very effective, and thermohaline mixing might be
therefore inhibited \citep{den09}, magnetic buoyancy might still
offer a viable mechanism for explaining extra-mixing effects. v)
O-rich AGB stars might experience extra-mixing processes
differently than C-rich ones, with the tendency to favour higher
Li abundances. We remember that also AGB luminosities
\citep{guabus} and N abundances \citep{sl90} suggest differences
between MS-S and C(N) stars, according to which only rarely the
former class of objects will evolve to the latter one. More often,
they will preserve a different surface composition and a different
C/O ratio because they have different initial masses.

In a different paper of this volume \citep{palm2} we examine
the consequences of the mixing schemes here discussed for the isotopic
ratios of CNO elements and for $^{26}$Al.

\section*{Acknowledgments}
This research was supported by the Italian Ministry of Research
under contract PRIN2006-022731. The idea of mixing by magnetic
buoyancy was developed in common with G.J. Wasserburg and K.M.
Nollett.

Solar and Stellar Astrophysics (astro-ph.SR)


Enhanced Magneto-optical Kerr Effect at Fe/Insulator Interfaces

Abstract: Using density functional theory calculations, we have found an enhanced magneto-optical Kerr effect in Fe/insulator interfaces. The results of our study indicate that interfacial Fe atoms in the Fe films have a low-dimensional nature, which causes the following two effects: (i) The diagonal component $\sigma_{xx}$ of the optical conductivity decreases dramatically because the hopping integral for electrons between Fe atoms is suppressed by the low dimensionality. (ii) The off-diagonal component $\sigma_{xy}$ of the optical conductivity does not change at low photon energies, but it is enhanced at photon energies around 2 eV, where we obtain enhanced orbital magnetic moments and spin-orbit correlations for the interfacial Fe atoms. A large Kerr angle develops in proportion to the ratio $\sigma_{xy}/\sigma_{xx}$. Our findings indicate an efficient way to enhance the effect of spin-orbit coupling at metal/insulator interfaces without using heavy elements.

\section{Introduction}
Magneto-optical Kerr effect (MOKE) is the phenomenon in which the light reflected from a magnetized material has a rotated plane of polarization.
Given its sensitivity, local probing nature, and experimental simplicity, this phenomenon has significantly impacted research on magnetic
materials~\cite{MOKE-Bader}. The MOKE originates from spin-orbit coupling in materials and has been extensively studied.
It has also been applied for magneto-optical data recording~\cite{MOKE-Data}. Recently, the effect has played a significant role in the rapidly developing field of spintronics. For example, the first experimental observation of the spin Hall effect~\cite{MOKE-SHE}, direct experimental observation of the skyrmion Hall effect~\cite{MOKE-SKX-Hall}, and investigations of spin-orbit torques in metallic and insulating magnetic bilayers~\cite{MOKE-SOT-1, MOKE-SOT-2, MOKE-SOT-3} were all performed using the MOKE. The effect is not limited to ferromagnetic materials, and a substantial MOKE has been proposed recently for some antiferromagnetic materials from theoretical calculations~\cite{MOKE-AF-1, MOKE-AF-2}.

For applications to magneto-optical devices, a substantial MOKE is required; this depends on the availability of materials with considerable spin-orbit coupling. Because it is a relativistic effect, spin-orbit coupling is small in many materials, but it is substantial in heavy elements such as Pt. To obtain large spin-orbit coupling, many different effects have recently been applied. These effects are as follows: spin-orbit splitting in the band structure due to the Rashba effect in systems without inversion symmetry~\cite{LSOC-Rashba-Fert, LSOC-Rashba-Zhang,LSOC-Rashba-Barnes}; charge-spin conversions due to spin-momentum locking of topological surface states~\cite{LSOC-TI-Lock-Fan, LSOC-TI-Lock-Shiomi, LSOC-TI-Lock-Ralph, LSOC-TI-Lock-Otani}; Dzyaloshinskii--Moriya interactions in bilayer systems with heavy elements~\cite{LSOC-Heavy-Fert,LSOC-Heavy-Bode, LSOC-Heavy-Heinze, LSOC-Heavy-Emori, LSOC-Heavy-Luchaire}; spin-orbit couplings due to impurities~\cite{LSOC-Imp-Guo, LSOC-Imp-AuFe, LSOC-Imp-AuPt, LSOC-Imp-CuBi-Niimi, LSOC-Imp-CuBi-Gu, LSOC-Imp-CuIr}.

There are a few ways to enhance the MOKE. One is to employ alloys~\cite{Kerr-Sugimoto,DFT-Kerr-Osterloh} or multilayers~\cite{DFT-Kerr-Guo,DFT-Kerr-Ebert} of transition metals such as Fe or Co and heavy elements such as Pt because heavy elements have large spin-orbit coupling.
Another way to enhance the MOKE is to utilize the plasma edge effect~\cite{DFT-Kerr-Feil-Haas}. By decreasing the diagonal component $\sigma_{xx}$ of the optical conductivity, the MOKE is enhanced because the Kerr rotation angle is proportional to $1/\sigma_{xx}$. Such a mechanism has been discussed to enhance the MOKE in Fe/Cu bilayers~\cite{Kerr-Katayama}. The photonic effect can also be used to enhance the MOKE because of multiple interference of light within the magnetic multilayers~\cite{Kerr-photonic-Inoue, Kerr-photonic-Fedyanin, Kerr-photonic-Koba}. Cavity enhancement of the MOKE has also been reported, in which a dielectric layer acts as a Fabry--Perot optical cavity~\cite{Kerr-cavity-2004, Kerr-cavity-2014, Kerr-cavity-2017}. Metamaterials, i.e., composite nanostructured materials, can also show enhancement of the MOKE~\cite{Kerr-meta-2007, Kerr-meta-2014, Kerr-meta-2017}.

\begin{figure}[tbp]
\includegraphics[width = 8.5 cm]{fig1.eps}
\caption{(a) The spin moment and (b) the orbital moment of a ten-monolayer supercell of Fe (Fe-10ML) as obtained by density functional theory calculations.
The 1st-ML and the 10th-ML are Fe/vacuum interfaces. We performed structural relaxation in the calculation.
For reference, the calculated value for bulk body-centered cubic iron (Fe-bcc) is shown by the dashed line. }
\label{F-spin-orbit-fe}
\end{figure}

In this paper, we propose an efficient way to enhance the MOKE using Fe/insulator interfaces.
In contrast to previous methods, this new method does not use heavy elements. Instead, it focuses on spin-orbit interactions, which can be enhanced by our method, while they have been ignored in previous applications of plasmonic or photonic effects.
We were motivated by the discovery of large magnetic anisotropies in ultrathin Fe/MgO films~\cite{Fe-MgO-Li, Fe-MgO-Boubeta,Fe-MgO-Beltran, Fe-MgO-Baumann, Fe-MgO-Bose}. The Fe/MgO interface is used in a wide range of devices, and enhanced spin-orbit interactions are expected because of the large magnetic anisotropy. In this paper, we use density functional theory calculations to show that atomically thin Fe layers as well as Fe/insulator multilayers comprising a few Fe layers and a few insulating layers of MgO or AlF$_3$ can produce a large MOKE and large spin-orbit interactions. Our results indicate an efficient way to obtain enhanced spin-orbit coupling effects at metal/insulator interfaces without using heavy elements.

\section{Density functional theory calculation results}

We first discuss the spin and orbital magnetic moments at a metallic Fe surface.
Figure \ref{F-spin-orbit-fe} shows the monolayer (ML)-resolved spin and orbital moments of a ten-monolayer Fe supercell (Fe-10ML) at the Fe (001) surface.
The 1st-ML and the 10th-ML are the two surface layers (Fe/vacuum interfaces), and we take the vacuum along the [001] direction to be 6ML thick. We performed density functional theory calculations using the \texttt{WIEN2K} package~\cite{DFT-Wien2k}.
In the calculations, the exchange-correlation interactions are described by the Perdew--Burke--Ernzerhof generalized gradient approximation~\cite {DFT-PBE}, and spin-orbit coupling is included using the second-variation method~\cite{DFT-SOC}. We used the atomic sphere radius parameter RMT = 2.24 for Fe, and employed the cutoff parameter RK$_{max}$ = 9. We used 12 $\times$ 12 $\times$ 1 k-point sampling for the Fe-10ML calculation. The magnetization direction is along the [001] direction. We performed structural relaxation for bulk body-centered cubit iron (Fe-bcc) and for Fe-10ML. In the Fe-10ML supercell, the lattice of MLs Nos. 1--4 was fixed to be equal to that of bulk Fe-bcc, the z component of the lattice of MLs Nos. 5--10 was relaxed, and the in-plane components of the lattice were fixed as those of bulk Fe-bcc.
We found the spin and orbital magnetic moments to be strongly enhanced in comparison to the values of bulk Fe-bcc. While this enhancement is limited mainly to the two surface layers, the spin and orbital moments of the internal layers decrease dramatically, thus approaching the value of bulk Fe, as shown in Fig. \ref{F-spin-orbit-fe}. This suggests that the surface Fe atoms are in a low-dimensional state intermediate between that of bulk Fe, which has small spin and orbital moments, and that of a single Fe atom state, for which the 3$d^6$ state has rather large spin and orbital moments.
Our result is in good agreement with previous calculations pertaining to the spin and orbital moments of 3$d$ transition metals in both surface and bulk cases~\cite{Bruno, Eriksson}.

Given that the MOKE originates from spin-orbit coupling, we investigated whether the enhanced spin and orbital moments at the Fe surface can induce a substantial Kerr effect. The Kerr angle is given as follows:
\begin{equation}
  \theta_{Kerr}\left(\omega\right) = -{\text Re}\frac{\epsilon_{xy}}{(\epsilon_{xx}-1)\sqrt{\epsilon_{xx}}},
   \label{E-kerr}
\end{equation}
where $\epsilon_{xx}$ and $\epsilon_{xy}$ are the diagonal and off-diagonal components of the dielectric tensor $\epsilon$, and $\omega$ is the photon energy, respectively. The dielectric tensor $\epsilon$ and the optical conductivity tensor $\sigma$ are related as follows:
\begin{equation}
\sigma\left(\omega\right) =\frac{\omega}{4\pi i}\left[\epsilon\left(\omega\right)-I\right],
\label{E-sigma}
\end{equation}
where $I$ is the unit tensor.
We performed density functional theory calculations with the \texttt{QUANTUM ESPRESSO} package~\cite{DFT-QE} along with maximally localized Wannier function calculations using the \texttt{wannier90} tool~\cite{DFT-W90} to obtain the optical conductivity tensor $\sigma$ and the Kerr angle $\theta_{Kerr}$ for the following three different cases: bulk Fe-bcc, Fe-10ML, and Fe-2ML.
In the calculations, the exchange-correlation interactions are described by the Perdew--Burke--Ernzerhof generalized gradient approximation, electron-ion interactions are represented by the Rabe--Rappe--Kaxiras--Joannopoulos ultrasoft pseudopotential, and spin-orbit coupling is included~\cite{QE-UPF}. The kinetic energy cutoff parameters for the wavefunctions (ecutwfc) and for the charge density and potential (ecutrho) are taken to be ecutwfc = 60 Ry and ecutrho = 600 Ry for Fe-bcc, ecutwfc = 60 Ry and ecutrho = 800 Ry for Fe-10ML, and ecutwfc = 60 Ry and ecutrho = 1200 Ry for Fe-2ML. In the self-consistent calculations, the k-point samplings were 16 $\times$ 16 $\times$ 16 for Fe-bcc, 12 $\times$ 12 $\times$ 1 for Fe-10ML, and 12 $\times$ 12 $\times$ 6 for Fe-2ML. To evaluate the optical conductivity tensor in the Wannier functional calculations, we used 25 $\times$ 25 $\times$ 25 k-point sampling. The magnetization direction is along the [001] direction. We used the structures previously obtained for bulk Fe-bcc and Fe-10ML. For the Fe-2ML calculation, we took the vacuum along the [001] direction to be 2ML thick, relaxed the z component of the lattice, and fixed the in-plane components of the lattice to be equal to those of bulk Fe-bcc.

As shown in Figs. \ref{F-kerr-surf}(a)--(c), for the bulk Fe-bcc case, our results are in good agreement with experimental measurements of the diagonal component $\sigma_{xx}$ of the optical conductivity~\cite{Exp-Fe-xx}, the off-diagonal component $\sigma_{xy}$ of the optical conductivity~\cite{Exp-Fe-xy}, and the Kerr angle $\theta_{Kerr}$~\cite{Exp-Fe-Kerr}. Moreover, our results are consistent with the previous calculations of $\sigma$ and $\theta_{Kerr}$ for the case of bulk Fe-bcc~\cite {DFT-Kerr-Oppeneer, DFT-Kerr-Guo, DFT-Kerr-Ebert, DFT-Kerr-Mainkar, DFT-Kerr-Delin, DFT-Kerr-Yao, DFT-Kerr-Yates, DFT-Kerr-Sangalli}. Interestingly, compared to the bulk Fe-bcc case the calculated results for $\sigma_{xx}$, $\sigma_{xy}$, and $\theta_{Kerr}$ for the Fe-10ML case are closer to the experimental values, as shown in Figs. \ref{F-kerr-surf}(a)--(c). We note that intraband transitions are not included in our calculations. The intraband contribution, as represented by the Drude formula, has usually been added to obtain good comparisons between the numerical and experimental results for the MOKE in metallic Fe, Co, and Ni using the Drude parameters extracted from experimental data~\cite{DFT-Kerr-Delin}. Interestingly, the good agreement between our calculations for the Fe-10ML case and the experiments provides an alternative way to understand the Kerr experiment in metallic Fe when the Drude contribution is not included. In addition, the Fe-10ML case is still dominated by Fe atoms with bulk-like properties, as shown in Fig. \ref{F-spin-orbit-fe}; thus there is no enhancement of $\theta_{Kerr}$ in the Fe-10ML case.

\begin{figure}[tbp]
\includegraphics[width = 8.5 cm]{fig2.eps}
\caption{(a) The real part of the diagonal component Re($\sigma_{xx}$) of the optical conductivity; (b) the imaginary part of the off-diagonal component Im($\sigma_{xy}$) of the optical conductivity multiplied by the photon energy $\omega$; (c) the Kerr angle $\theta_{Kerr}$, each as a function of the photon energy $\omega$. Density functional theory calculations for the following three cases are shown: bulk Fe-bcc, Fe-10ML, and Fe-2ML.
The experimental values of Re($\sigma_{xx}$) \cite{Exp-Fe-xx}, $\omega$Im($\sigma_{xy}$) \cite{Exp-Fe-xy}, and $\theta_{Kerr}$ \cite{Exp-Fe-Kerr} for bulk Fe-bcc are shown for comparison.}
\label{F-kerr-surf}
\end{figure}

These findings for the Fe-10ML case led us to consider the Fe-2ML case, which consists entirely of surface Fe atoms, i.e., the two monolayers are the two surfaces (Fe/vacuum interfaces). As shown in Fig. \ref{F-kerr-surf}(c), quite a large Kerr angle $\theta_{Kerr}$ is obtained for the Fe-2ML case, particularly for photon energies $\omega$ below 2 eV. Two factors rooted in the Fe surface produce this large Kerr angle.
First, the diagonal component $\sigma_{xx}$ of the optical conductivity decreases dramatically, as shown in Fig. \ref{F-kerr-surf}(a),
because the hopping integral for electrons between Fe atoms is suppressed owing to the low dimensionality of the atomically thin Fe layers.
Second, the off-diagonal component $\sigma_{xy}$ of the optical conductivity, in contrast to $\sigma_{xx}$, does not change much for photon energies $\omega<$  1eV, and it even increases above the value of $\sigma_{xy}$ for bulk Fe at photon energies close to $\omega \sim$ 2eV, as shown in Fig. \ref{F-kerr-surf}(b).
According to Eqs. (\ref{E-kerr}) and (\ref{E-sigma}), the large values of $\sigma_{xy}/\sigma_{xx}$ and $\theta_{Kerr}$ are obtained because of decreased $\sigma_{xx}$ and large $\sigma_{xy}$.

\begin{figure}[tbp]
\includegraphics[width = 8.5 cm]{fig3.eps}
\caption{(a) The real part of the diagonal component Re($\sigma_{xx}$) of the optical
conductivity; (b) the imaginary part of the off-diagonal
component Im($\sigma_{xy}$) of the optical conductivity multiplied by the photon energy $\omega$;
(c) the Kerr angle $\theta_{Kerr}$, each as a function of
the photon energy $\omega$, for Fe/MgO and Fe/AlF$_3$ multilayers, where the
Fe layer is two monolayers thick. The calculated results for bulk Fe-bcc are shown for reference.}
\label{F-kerr-2ml}
\end{figure}

Motivated by the enhanced Kerr effect we found for the Fe-2ML case, we continued with studies of other Fe/insulator interfaces. To maintain the large Kerr effect we obtained for an Fe layer comprising solely interfacial Fe atoms, we fixed the Fe layer to be two monolayers thick.
For the Fe/insulator multilayers, we added 2ML of MgO or AlF$_3$ along the [001] direction as the insulator layers. We again relaxed the z component of the lattice, and we fixed the in-plane components of the lattice to be equal to those of bulk Fe-bcc. As shown in Fig. \ref{F-kerr-2ml}(c), we obtained a large Kerr angle, particularly for photon energies $\omega$ below 2 eV owing to the same mechanism as in the previous Fe-2ML case, i.e., the diagonal component $\sigma_{xx}$ of the optical conductivity decreases dramatically, as shown in Fig. \ref{F-kerr-2ml}(a), while the off-diagonal component $\sigma_{xy}$ of the optical conductivity does not change for photon energies lower than 0.5 eV, and it even increases above the value for bulk Fe at photon energies around 2 eV, as shown in Fig. \ref{F-kerr-2ml}(b). These calculations indicate that the enhancement of the Kerr effect in Fe/insulator multilayers is rather robust against the choice of insulator.

\begin{figure}[tbp]
\includegraphics[width = 8.5 cm]{fig4.eps}
\caption{(a) Spin moment $S^z$, (b) orbital moment $L^z$, and (c) spin-orbit correlation $S^zL^z$ for Fe-2ML, Fe/MgO and Fe/AlF$_3$ multilayers
calculated using density functional theory. The structures are the same as those in Figs. \ref{F-kerr-surf} and \ref{F-kerr-2ml}. The calculated results for bulk Fe-bcc are shown for reference. }
\label{F-spin-orbit-2ml}
\end{figure}

Since $\sigma_{xy}$ originates from spin-orbit coupling in the systems, we studied the spin and orbital magnetic moments and their correlations for the Fe-2ML, Fe/MgO, and Fe/AlF$_3$ multilayers, for which the structures are the same as those mentioned above.
As shown in Fig. \ref{F-spin-orbit-2ml}(a), the spin moments $S^z$ of Fe-2ML and of Fe/AlF$_3$ increase to $\sim$ 3 $\mu_B$, while those of Fe/MgO remain close to 2.2 $\mu_B$, which is the calculated value for bulk Fe. As shown in Fig. \ref{F-spin-orbit-2ml}(b), the orbital moment $L^z$ increases in all three cases. Consequently, we define an effective spin-orbit correlation, as shown in Fig. \ref{F-spin-orbit-2ml}(c), which is the product of $S^z$ in Fig. \ref{F-spin-orbit-2ml}(a) and $L^z$ in Fig. \ref{F-spin-orbit-2ml}(b). The spin-orbit correlation $S^zL^z$ increases for all three cases.

\begin{figure}[tbp]
\includegraphics[width = 8.5 cm]{fig5.eps}
\caption{ Density of states (DOS) for bulk Fe-bcc and for Fe-2ML calculated using density functional theory. The Fermi level is $E_F$ = 0 eV.}
\label{F-dos-fe}
\end{figure}

To understand the reduction of $\sigma_{xx}$ in low dimensions, we studied the density of states (DOS) for bulk Fe-bcc and for Fe-2ML, for which the structures are the same as described above. As shown in Fig. \ref{F-dos-fe}, compared to the DOS of Fe-bcc, the DOS of Fe-2ML for spin up (majority spin) below the Fermi level ($E_F$ = 0 eV) becomes narrow and shifts down about 1.5 eV, and the DOS of Fe-2ML for spin down (minority spin) below the Fermi level becomes narrow and shifts up about 1.5 eV. The effect of the reduced dimensionality is a narrowing of the bands and an increase in the exchange splitting. This in turn decreases the diagonal part $\sigma_{xx}$ of the optical conductivity.

\section{Discussion}
We briefly compare our results with those of the previous studies of the enhancement of the Kerr angle $\theta_{Kerr}$.
It has been suggested that by considering alloys~\cite{Kerr-Sugimoto,DFT-Kerr-Osterloh} or multilayers~\cite{DFT-Kerr-Guo,DFT-Kerr-Ebert} of transition metals such as Fe or Co and heavy element such as Pt, $\theta_{Kerr}$ can be enhanced because heavy elements have large spin-orbit coupling. By contrast, we have found a way to enhance $\theta_{Kerr}$ and spin-orbit correlations {\it without using heavy elements} because avoiding the use of heavy elements is ideal for applications. Another mechanism for enhancing $\theta_{Kerr}$ is the plasma edge effect~\cite{DFT-Kerr-Feil-Haas}, which requires a dielectric constant satisfying the relation Re(${\epsilon_{xx}}$) = 1. Since $\theta_{Kerr}$ is proportional to $\epsilon_{xy}/(\epsilon_{xx}-1)$ as defined in Eq.(\ref{E-kerr}), the condition Re($\epsilon_{xx}$) = 1 induces a large $\theta_{Kerr}$ if Im($\epsilon_{xx}$) is small. Considering the relationship between optical conductivity and dielectric constant, $\sigma_{xx}=\frac{\omega}{4\pi i}(\epsilon_{xx}-1)$ given in Eq. (\ref{E-sigma}), the condition Re($\epsilon_{xx}$) = 1 and the condition Im($\sigma_{xx}$) = 0 are the same. Our study proposes a new method for enhancing $\theta_{Kerr}$ by considering Fe/insulator interfaces. Our study simultaneously considers both $\sigma_{xx}$ and $\sigma_{xy}$ on an equal footing, while the effect of $\sigma_{xy}$ has been ignored in studies that focused on the the plasma edge effect~\cite{DFT-Kerr-Feil-Haas}.

\section{Conclusions}
In summary, using density functional theory calculations, we have found an enhanced MOKE in Fe/insulator multilayers.
The Kerr effect is rather robust against the choice of insulator (vacuum, MgO, or AlF$_3$ in this study), but it requires the Fe layer to be atomically thin and to be dominated by interfacial Fe atoms. Our calculations suggest that the enhanced Kerr effect originates from the low-dimensional nature of the interfacial Fe atoms, which causes the following two effects:
(i) The diagonal component $\sigma_{xx}$ of the optical conductivity decreases dramatically because the hopping integral for electrons between the Fe atoms is suppressed by the low dimensionality of the atomically thin Fe films.
The off-diagonal component $\sigma_{xy}$ of the optical conductivity does not change at low photon energies, and it is enhanced at photon energies around 2 eV, where enhanced orbital magnetic moments and spin-orbit correlations are obtained for the interfacial Fe atoms. The Kerr angle $\theta_{Kerr}$ is proportional to the ratio $\sigma_{xy}/\sigma_{xx}$ by definition, and a large $\theta_{Kerr}$ develops at Fe/insulator interfaces with decreasing $\sigma_{xx}$ and large $\sigma_{xy}$.
Our results show an efficient way to design novel devices, such as ultrathin ferromagnetic films~\cite{Okada} and ferromagnetic nanogranular films~\cite{Kobayashi-1,Kobayashi-2} by manipulating the substantial spin-orbit coupling effect at metal/insulator interfaces without using heavy elements.

\section*{Acknowledgment}
The authors acknowledge N. Kobayashi, K. Ikeda, and H. Masumoto for many valuable discussions about experiments on ferromagnetic nanogranular films.


The Far-infrared Continuum of Quasars

Abstract: ISO provides a key new far-infrared window through which to observe the multi-wavelength spectral energy distributions (SEDs) of quasars and active galactic nuclei (AGN). It allows us, for the first time, to observe a substantial fraction of the quasar population in the far-IR, and to obtain simultaneous, multi-wavelength observations from 5--200 microns. With these data we can study the behavior of the IR continuum in comparison with expectations from competing thermal and non-thermal models. A key to determining which mechanism dominates, is the measurement of the peak wavelength of the emission and the shape of the far-IR--mm turnover. Turnovers which are steeper than frequency^2.5 indicate thermal dust emission in the far-IR.Preliminary results from our ISO data show broad, fairly smooth, IR continuum emission with far-IR turnovers generally too steep to be explained by non-thermal synchrotron emission. Assuming thermal emission throughout leads to a wide inferred temperature range of 50-1000 K. The hotter material, often called the AGN component, probably originates in dust close to and heated by the central source, e.g. the ubiquitous molecular torus. The cooler emission is too strong to be due purely to cool, host galaxy dust, and so indicates either the presence of a starburst in addition to the AGN or AGN-heated dust covering a wider range of temperatures than present in the standard, optically thick torus models.

\section{INTRODUCTION AND ANALYSIS}
The energy output of quasars and AGN is comparable
thoughout the far-IR to X-ray spectral regions, so studies of their
emission mechanisms require observations throughout this range.
Prior to ISO, the far-IR continuum was observed by IRAS, providing our
first look at quasars in this wavelength range and showing them to be strong
emitters from 12--100\,$\mu$m. The poor spatial resolution and limited
wavelength range of IRAS restricted observations to the IR-brightest subsets
of quasars and AGN and left open many questions as to the range of behavior
present and the emission mechanisms responsible in the various types of
AGN. With ISO we are now able to address some of these limitations by
studying a large number of different kinds of AGN and quasars.

We observed 72 quasars and AGN covering a full range of
redshift and luminosity and many different types of SEDs.
The observations were made with ISOPHOT (\cite{phot96}) and
cover 5--200\,$\mu$m in 8 bands. The status of this
project, including a description of the sample and global comparisons of
our results with those of IRAS, is summarised elsewhere (Hooper {\it et~al.} ,
this volume). Here we discuss our preliminary results in
the context of the scientific questions we wish to address, with
particular emphasis on the high-redshift objects in the sample.

The data were processed with PIA (\cite{Gabriel97})\footnote{PIA is a joint
development by the ESA Astrophysics Division and the ISOPHOT Consortium}
through the AAP level, from which custom
scripts (originating with M. Haas \& S. Mueller at MPIA) were used to
determine source fluxes and uncertainties.  Default vignetting values from
PIA were applied to the chopped data, and the flux scale was set by FCS
observations.  Sky noise estimates based on \cite{herb98}
were included in the error budget.
For further details of the data reduction please refer to Hooper {\it et~al.}\ (this
volume).

\section{SPECTRAL ENERGY DISTRIBUTIONS (SEDs).}
Figure~\ref{fig:seds} shows typical examples of the SEDs of a radio-loud
(upper, RLQ) and radio-quiet (lower, RQQ) quasar (\cite{ewm94}).
Throughout the far-IR--ultraviolet (UV) spectral region, the SEDs are
remarkably similar, showing the big blue bump in the optical-UV, the
broad IR bump, and an inflection between the two centered at $\sim 1$\,$\mu$m.
The differences between the two radio classes are evident in the radio
region where, by definition, the RLQs are several orders of
magnitude more luminous, and in the soft X-ray region where the RLQs are
a factor of $\sim 3$ brighter and have a harder spectrum. This latter
difference is generally ascribed to an additional synchrotron
self-Compton component linked to the radio emission (see~\cite{Wilkes99}
for a review of quasar SEDs).

\begin{figure}[h]
    \leavevmode
  \centerline{\epsfig{file=bwilkes1.ps,height=3.8truein,angle=-90}}
  \caption{\em Rest-frame radio--X-ray spectral energy distributions (SEDs) for
low-redshift radio-loud (upper, RLQ) and radio-quiet (lower, RQQ)
quasars (Fig.1 from \protect\cite{ewm94}). These are typical
examples of SEDs in the
current literature. The display of $\nu$L$_{\nu}$ vs $\nu$ highlights the
energy output and the structure of the SED. }
  \label{fig:seds}
\end{figure}

The radio--IR SEDs of core-dominated RLQs are generally smooth
while in both RQQs and lobe-dominated RLQs
the far-IR cutoff is steep and does not extrapolate into the radio.
This difference has resulted in the general belief that the IR emission in
core-dominated RLQs is dominated by beamed, non-thermal and
in RQQs by thermal emission mechanisms (see~\cite{Wilkes99} for a
review). The median RLQ SED in
Figures~\ref{fig:1351},\ref{fig:1202},\ref{fig:1946}
was generated from a combination of lobe- and
core-dominated RLQs. In a unification scenario
(see~\cite{Barthel99} for a review) the beamed IR component
is expected to be weaker in lobe-dominated sources where the beam is in
the plane of the sky. Thus it may not be surprising that the RLQ SED,
formed by combining
SEDs from the two radio subclasses, is similar to that of
the RQQs. A detailed comparison of the SEDs of core- and
lobe-dominated RLQs is required to confirm this scenario.
Such a transition from non-thermal to thermal emission with decreasing
core-dominance is demonstrated
in the preliminary ISOPHOT results of the European Core Program on Quasars
(\cite{Haas98}; see also this volume).

The main questions we wish to address with our ISO sample are (i) to better
understand the interplay between
thermal and non-thermal emission in the various types of source (ii)
to investigate the
relative contributions of cool (galactic/starburst) and warm (AGN-related)
dust (iii) to characterize and measure the continua in order to
constrain models (iv) to study the evolution of the IR SEDs.

\begin{figure}[h]
    \leavevmode
\centering
\vskip -0.2in
  \centerline{\epsfig{file=bwilkes2u.eps,height=3.5truein}}
\vskip -1.2in
  \centerline{\epsfig{file=bwilkes2l.eps,height=3.5truein}}
\vskip -1.0in
  \caption{\em The rest-frame SED of the low-redshift (z=0.088), RQQ
PG1351+640 from \protect\cite{ewm94} with the ISO data points (open triangles)
added. {\bf Upper}: the radio--UV SED showing a
steep ($\alpha_{fir} > 2.9$) far-IR cut-off;
{\bf Lower}: the IR--optical SED
with representative grey body (assuming $\beta = 2$)
curves superposed.
}
  \label{fig:1351}
\end{figure}

\subsection{LOW REDSHIFT QUASAR: PG1351+640}
\label{sec:1351}

Figure~\ref{fig:1351} shows the rest-frame SED of a low-redshift (z=0.088),
RQQ PG1351+640. The ISO data points are shown as open triangles and the
remainder of the multi-wavelength data are from~\cite{ewm94}.
This SED is typical of the low-redshift RQQs in our sample and
compares well with the low-redshift median.
The current limit on the
far-IR turnover is $\alpha_{fir} > 2.9$, too steep to be due to
non-thermal synchrotron emission, for which $\alpha_{fir} < 2.5$ is required
and $\alpha < 1$ is more typical (\cite{Gear94}). We thus conclude that
the far-IR emission in this quasars is dominated by thermal emission from dust.
\begin{table*}[t]
\centering
\caption{The Broad-band IR luminosities of the High redshift Quasars.}
\label{tab:lum}
\footnotesize
\begin{tabular}{lllllc}
\hline
Name & Redshift & L(7--25\,$\mu$m)$^1$ & L(2--25\,$\mu$m)$^1$ & L(2--200\,$\mu$m)$^1$ &
\# ISO det'ns \\
\hline
1107+48 & 3.010 & 1.7E+47 & 6.9E+47 & - & 2 \\
BR 1202$-$0727$^2$ & 4.690 & 1.7E+48 & 8.5E+48 & 8.9E+48 & 3 \\
1407+265 & 0.944 & 7.5E+46 & 1.2E+47 & - & 2\\
1413+1143 & 2.551 & 8.8E+47 & 1.6E+48 & 2.0E+48 & 7\\
1718+481 & 1.084 & 2.7E+46 & 8.2E+46 & - & 3\\
HS 1946+7658 & 3.030 & 4.2E+47 & 1.1E+48 & 2.1E+48 & 4 \\
3C422 & 0.942 & 7.6E+46 & 1.4E+47 & - & 2 \\
\hline
1351+640 & 0.088 & 2.5E+45 & 3.1E+45 & 5.2E+45 & 6 \\
\hline
\end{tabular}
\begin{minipage}{5.5truein}
1: Luminosity over indicated bandpass, assuming $H_o$=50 km s$^{-1}$ Mpc$^{-1}$
and $q_o$=0 (erg s$^{-1}$) \\
2: Probably a gravitationally lensed source (\cite{Omont96})
\end{minipage}
\end{table*}

The IR SED is smooth and indicates a wide and probably continuous
range of temperatures: 40--500K. In a purely thermal scenario this
smooth continuum can be explained by
a combination of cool host galaxy dust, warm ($\sim 100$K) starburst-related
dust (e.g.\ \cite{RR95}) and hot (T$_{eff} \sim 800-1300$K)
AGN-related dust (e.g.\ from a torus, \cite{PK93}).
Alternatively the emission could originate in
AGN-related dust which is sufficiently optically thick and/or geometrically
arranged (e.g.\ the warped disk model of \cite{Sand89} or a patchy torus) to
produce comparable emission over the wide range of observed
temperatures. This latter scenario is supported by the results of van
Bemmel \& Barthel (this volume) which demonstrate that RLQs have
significantly brighter far-IR emission than radio galaxies in samples
which were matched
so that the only difference should be one of orientation.

We also note that, with these data alone, we cannot rule out a significant
contribution from a non-thermal component in the mid- and near-IR.

\section{HIGH REDSHIFT QUASARS}

\begin{figure}[h]
    \leavevmode
\vskip -0.2in
  \centerline{\epsfig{file=bwilkes3u.eps,height=3.5truein}}
\vskip -1.2in
  \centerline{\epsfig{file=bwilkes3l.eps,height=3.5truein}}
\vskip -1.1in
  \caption{\em The rest-frame, radio--UV SED of the high-redshift, RQQ
BR 1202$-$0727 including the ISO data points (open triangles). The
median SED for low redshift RQQs (\protect\cite{ewm94}), normalised to
the optical SED, is superposed
for comparison. {\bf Upper}:
grey-body curves (assuming $\beta$=2)
at two extreme temperatures are superposed
{\bf Lower}: the SED of an ULIRG is superposed,
normalised to the highest H-band luminosity host galaxy observed in low redshift
sources (\protect\cite{McLeod97}).}
  \label{fig:1202}
\end{figure}

Our sample includes 9 quasars with redshift $> 0.9$. The 7 for which we have
reliable ISO detections are listed in Table~\ref{tab:lum} along with their
redshifts, preliminary estimates of several broad-band IR luminosities
and the number of ISO detections used in this determination.
The luminosities were determined via linear interpolation between the
data points and are not tabulated where there were insufficient
data to provide useful constraints. The IR-band luminosities,
determined assuming $H_o$=50 km s$^{-1}$ Mpc$^{-1}$
and $q_o$=0, are high,
ranging from 10$^{46-49}$ erg s$^{-1}$ ($\sim 10^{13-16} {\rm L_{\odot}}$).
For comparison, the low-redshift RQQ,
PG1351+640 (Table~\ref{tab:lum}) has luminosities $\sim 10^{45}$ erg s$^{-1}$
($\sim 10^{12} {\rm L_{\odot}}$).

The rest-frame, radio--UV SED of BR 1202$-$0727, the object with the
highest luminosity, L(2--200\,$\mu$m)
$\sim 10^{49}$ erg s$^{-1}$, is displayed in
Figure~\ref{fig:1202}.
Our ISO fluxes and upper limits are shown as open triangles
and our ground-based near-IR photometry,
obtained at the Steward Observatory 61$"$ telescope in March 1996, as open
squares. Multiwavelength data are from~\cite{Isaak94} and
\cite{McMahon94}. The median, low-redshift, RQQ SED (\cite{ewm94}), normalised
at $\sim 2$\,$\mu$m, is superposed for comparison.
Relative to the rest-frame optical emission, the IR emission in this
object is $\sim 2$ orders of magnitude stronger than is typical.
This is a complex source. It is a luminous (dust
$\sim 10^{11}{\rm M_{\odot}}$, \cite{Ohta96}), double (\cite{Omont96}),
CO source with a Ly$\alpha$ emission line companion (\cite{Hu96}).
Given its unprecendentedly strong mm and IR luminosities and double
nature it seems highly likely that the source is graviationally lensed
and so amplified by $\sim \times 10$. However lensing
cannot explain the unusual SED of BR 1202-0727, since partial lensing would
boost the smaller, optical/UV source rather than the
probably extended IR emitting region.

Two extreme temperature grey-body curves (assuming an emissivity index,
$\beta =2$) are also superposed in Figure~\ref{fig:1202} (upper)
for illustration. Again a wide range of temperature is present, $50-1000$K
although the upper limits in the mid-IR prevent any constraint on the
emission at intermediate temperatures in this source.
In Figure~\ref{fig:1202} (lower) an SED representing a ULIRG with the
maximum H-band host galaxy luminosity based on low redshift quasars (\cite{McLeod97})
is superposed. This comparison
illustrates that the far-IR emission is comparable or in excess of that
expected from a ULIRG and far exceeds that expected from the host galaxy
(see also \cite{Wilkes97}).

The rest-frame, radio--UV SED of HS 1946+7658 is displayed in
Figure~\ref{fig:1946}.
Our ISO fluxes and upper limits are shown as open triangles and the
remainder of the multiwavelength data were taken
from~\cite{Kuhn94} and~\cite{Kuhn96}.
We also superpose grey-body curves (assuming dust with emissivity index,
$\beta =2$) at two extreme temperatures to illustrate the range of
temperature present in the data. The SED of HS 1946+7658
is remarkably similar to the low redshift, RQQ median, extending the
conclusion of~\cite{Kuhn94} based on the optical SED for this source.
As with the low redshift RQQ
PG1351+640, the IR continuum appears smooth.
In a pure thermal scenario, this suggests the presence of dust covering a
continuous range from $\sim <80-500$K rather than
two distinct low and high temperature components.
The low temperature limit is determined by the longest
wavelength at which we have data rather than a real cut-off in the
source.

\begin{figure}[h]
    \leavevmode
\vskip -0.2in
  \centerline{\epsfig{file=bwilkes4.eps,height=3.5truein}}
\vskip -1.1in
  \caption{\em The rest-frame, radio--UV SED of the high-redshift, RQQ
HS 1946+7658 including the ISO data points shown as open triangles.
The median SED for low redshift RQQs (\protect\cite{ewm94}) is superposed
for comparison. Grey-body curves (assuming $\beta$=2)
at two temperatures are superposed to illustrate the
range of temperature present in a pure thermal emission scenario.}
  \label{fig:1946}
\end{figure}

In summary, the broad-band IR luminosities, ranging
from 10$^{46-48}$ erg s$^{-1}$ ($\sim 10^{13-15} {\rm L_{\odot}}$, assuming
BR 1202-0727 is lensed),
are $\sim 2$ orders of magnitude
higher than those of Ultraluminous IR galaxies (ULIRGs) whose SEDs peak
in the IR with luminosities $\sim 10^{11-13} {\rm L_{\odot}}$ (e.g.\
\cite{Sand96}), even allowing for the different assumed Hubble
constant (\cite{Sand96} use $H_o = 75$ km s$^{-1}$ Mpc$^{-1}$).
In the far-IR these luminosities are too high to be due to the
host galaxy.
As for the low redshift sources (\ref{sec:1351}), the inferred temperature
ranges are generally too broad to be fit by standard, AGN
molecular torus models (\cite{PK93}), which can reproduce the
higher temperature emission. Thus more complex models or a strong starburst
component are required to provide the lower temperature emission.

\section*{ACKNOWLEDGMENTS}
We wish to thank the staff at the ISOPHOT data center in Heidelberg, in particular
Martin Haas, for their help and hospitality while working on these data.
We also thank the ISO staff at IPAC for their continual support and help
both during visits and remotely. This work is supported by NASA grant
NAGW-3134.

The Far-infrared Continuum of Quasars