Brown Dwarfs as Galactic Chronometers

Abstract: Brown dwarfs are natural clocks, cooling and dimming over time due to insufficient core fusion. They are also numerous and present in nearly all Galactic environments, making them potentially useful chronometers for a variety of Galactic studies. For this potential to be realized, however, precise and accurate ages for individual sources are required, a prospect made difficult by the complex atmospheres and spectra of low-temperature brown dwarfs; degeneracy between mass, age and luminosity; and the lack of useful age trends in magnetic activity and rotation. In this contribution, I review five ways in which ages for brown dwarfs are uniquely determined, discuss their applicability and limitations, and give current empirical precisions.

\section{Introduction}

Brown dwarfs are very low-mass stars
whose masses (M $\lesssim$ 0.075~{M$_{\odot}$}) are insufficient to sustain
the core hydrogen fusion reactions that balance radiative energy losses \citep{1963ApJ...137.1121K,1963PThPh..30..460H}.  Supported from further gravitational contraction
by electron degeneracy pressure, evolved brown dwarfs continually cool
and dim over time as they radiate away their initial
contraction energy, ultimately
achieving photospheric conditions that can be similar to those of giant planets.
The first examples of brown dwarfs were identified as recently as
1995 \citep{1995Natur.378..463N,1995Natur.377..129R}.
Today, there are hundreds known in nearly
all Galactic environments, identified largely
in wide-field, red and near-infrared imaging surveys such as 2MASS, DENIS, SDSS and UKIDSS.
The known population of brown dwarfs encompasses
the late-type M ({T$_{eff}$} $\approx$ 2500--3500~K), L ({T$_{eff}$} $\approx$ 1400--2500~K) and T spectral classes ({T$_{eff}$} $\approx$ 600--1400~K; e.g., \citealt{2004AJ....127.2948V}),
while efforts are currently underway to find even colder
members of the putative Y dwarf class (see review by \citealt{2005ARA&A..43..195K}).

Because brown dwarfs cool over time, their spectral properties are inherently time-dependent,
making them potentially useful chronometers for Galactic studies (much like white dwarfs;
see contributions by M.\ Salaris, J.\ Kalirai, S.\ Catal{\'a}n and H.\ Richer).
However, the primary observables of a brown dwarf---temperature, luminosity and spectral type---depend on both mass and age (and weakly on metallicity).  This degeneracy complicates characterizations of individual sources and mixed
populations.  Unfortunately, traditional stellar age-dating methods do not appear to be applicable
for brown dwarfs.
Magnetic activity metrics, such as the frequency and strength
of optical or X-ray non-thermal emission, appear to be largely age-invariant
(e.g., \citealt{2006A&A...448..293S}) and quiescent emission drops off
precipitously in the early L dwarfs as cool photospheres are decoupled from
field lines
\footnote{Interestingly, radio emission does not drop off for cooler brown
dwarfs, and may even
increase relative to bolometric flux, although there are currently few detections (e.g., \citealt{2006ApJ...648..629B}).}
(e.g., \citealt{2002ApJ...571..469M}; see contribution by A.\ West).  Long-term angular momentum loss in brown
dwarfs is far more muted than in stars, and there is no clear rotation-activity
relation for L dwarfs (e.g., \citealt{2008ApJ...684.1390R}).
Exploitation of the cooling properties of brown dwarfs is therefore a favorable approach
for determining their ages.

In this contribution, I review five methods currently employed to age-date brown dwarfs and summarize
their applicability, inherent limitations and current (typical) precisions.
Table~1 provides a summary of the methods discussed in detail below.

\begin{table}
  \begin{center}
  \caption{Age-dating Methods for Brown Dwarfs.}
  \label{tab1}
 {\scriptsize
  \begin{tabular}{|l|l|l|l|c|l|}
 \hline
{\bf Technique} &
{\bf Applicable for} &
{\bf Pros} &
{\bf Cons} &
{\bf Precision} &
{\bf Examples} \\
 & & & & &
{\bf in the} \\
 & & & & &
{\bf Literature} \\
 \hline
Cluster & nearby clusters; &  precise ages based on  & generally restricted to & $\sim$10\% clusters; & 1,2,3,4 \\
members \&  & companions to & stellar/cluster work; & young clusters; wide & $\sim$50--100\% for & \\
companions  & age-dated stars & calibration for other & (resolvable) \& close & companions & \\
  & & techniques and & (RV) companions rare; & & \\
  & & evolutionary models &  must verify coevality/ & & \\
  & & &  association; atmospheres & & \\
  & & &  variable for $t \lesssim$ 10~Myr & & \\
 \hline
6708~{\AA} Li~I & $t$ $\sim$ 10--200~Myr; &  consistent predictions   & requires high sensitivity, & $\sim$10\% for  & 5,6,7 \\
absorption  & resolved binaries; & from different models; & high resolution spectra; & young clusters; & \\
  & individual field  & straightforward test; & low brightness region; & upper/lower & \\
  & sources with & largely insensitive to  & not useful for T dwarfs & limits for all & \\
  & {T$_{eff}$} $>$ 1500~K & atmospheric properties & or for $t \lesssim$ 10~Myr; & others & \\
  & and $t \lesssim$ 2~Gyr  &  & relies on accurate & & \\
  & (limits only)  & & evolutionary models & & \\
 \hline
Mass & astrometric/RV & precise masses    & suitable systems are & $\sim$10--20\% & 8,9,10,11 \\
standards  & binaries & yield precise ages; & rare; long-term & & \\
  & &  weakly sensitive  & follow-up required; & & \\
  &  & to atmospheric & relies on accurate & & \\
  &  & properties & evolutionary models & & \\
 \hline
Surface & any source with & applicable to & low precision; other & $\sim$50--100\% & 12,13,14 \\
gravities  & a well-measured & individual sources;    & factors (e.g., metallicity) & & \\
 & spectrum & particularly useful & complicate analysis; & & \\
  & & for T dwarfs & relies on accurate & & \\
  &  & & evolutionary and & & \\
  &  & & atmospheric models & & \\
 \hline
Kinematics & well-defined & useful statistical test & very low precision; & $\sim$100--300\% & 15,16,17 \\
 & groups or & for various subclasses;   & large groups required & & \\
 & populations & insensitive to   & to beat statistical noise; & & \\
  &  & evolutionary or &  susceptible to & \\
  &  &  atmospheric models & selection biases & \\
 \hline
  \end{tabular}
  }
 \end{center}
\vspace{1mm}
 \scriptsize{
 {\it References:}
 (1) \citet{1998A&A...336..490B}; (2) \citet{2003ApJ...593.1093L}; (3) \citet{2001ApJ...556..373G};
 (4) \citet{2007ApJ...658..617B}; (5) \citet{1998ApJ...499L.199S}; (6) \citet{2005MNRAS.358...13J};
 (7) \citet{2005ApJ...634..616L}; (8) \citet{2004ApJ...615..958Z};
 (9) \citet{2006Natur.440..311S}; (10) \citet{2008arXiv0807.0238L}; (11) \citet{2008arXiv0807.2450D};
 (12) \citet{2004ApJ...609..854M}; (13) \citet{2006ApJ...639.1095B}; (14) \citet{2007ApJ...656.1136S};
 (15) \citet{2007AJ....133.2258S}; (16) \citet{2007ApJ...666.1205Z};
 (17) \citet{2008arXiv0809.3008F}.
 }
\end{table}

\section{Cluster Members and Companions}

The most straightforward way to age-date an individual brown dwarf is to borrow
from its environment, a tactic that is suitable for members of coeval
clusters/associations and companions to age-dated stars.
Brown dwarfs are well-sampled
down to and below the deuterium fusing mass limit (M $\lesssim$ 0.013~{M$_{\odot}$})
in the youngest nearby clusters ($t$ $\lesssim$ 5~Myr), as their
luminosities are greater at early ages. Brown dwarfs have
also been identified in somewhat older (10--50~Myr)``loose associations'' in the vicinity
of the Sun ($d$ $\lesssim$ 50--100~pc; e.g., \citealt{2008arXiv0808.3153K}).
For older and more distant clusters ($d$ $\gtrsim$ 1~kpc),
decreasing surface temperatures and compact radii exacerbate the sensitivity
issues that plague low-mass stellar studies (see contribution by G.\ Piotto).  There are as yet
no known brown dwarfs in globular clusters, despite detection of the end
of the stellar main sequence in systems such as NGC~6397 \citep{2008AJ....135.2141R}.

For brown dwarfs in young clusters,
numerous studies have examined age-related trends in
colors (e.g., \citealt{2008MNRAS.385.1771J}),
spectral characteristics (e.g., \citealt{2007ApJ...657..511A}),
accretion timescales (e.g., \citealt{2005ApJ...626..498M}) and
circum(sub)stellar disk evolution (e.g., \citealt{2007ApJ...660.1517S}).  These
have been coarsely quantified, and appear to be most useful at very
young ages ($t \lesssim$ 10~Myr).
Surface properties and luminosities are highly variable at these ages
due to sensitivity to formation conditions
(e.g., \citealt{2002A&A...382..563B}), ongoing accretion,
complex magnetic effects (e.g., \citealt{2007ApJ...671L.149R})
and age spreads within a cluster (see contribution by R.\ Jeffries).
Hence, while brown dwarfs in clusters with ages
spanning $\sim$1--650~Myr are now well-documented, with age uncertainties
as good as 10\% (for the LDB technique; see $\S$~3), useful predictive trends are probably
limited to ages of $\gtrsim$10~Myr.

Known brown dwarf companions to main sequence stars now number a few dozen,
spanning a wide range of separation, age and composition.  Many of these systems are widely-separated so that their brown dwarf
companions can be directly studied.
Age uncertainties for companion brown dwarfs depend on stellar dating methods which are generally
more uncertain (50--100\%; e.g., \citealt{2008arXiv0807.0238L}) than cluster ages.  Searches for substellar companions to more precisely age-dated
white dwarfs (e.g. \citealt{2008MNRAS.388..838D,2008ApJ...681.1470F}) and subgiants \citep{2006MNRAS.368.1281P} have so far met with limited success.
Nevertheless, brown dwarf companions to age-dated stars serve as important benchmarks for calibrating
other age-dating methods at late ages ($\gtrsim$500~Myr) and are fundamental
for testing evolutionary models (see contribution by T.\ Dupuy).

\section{6708~{\AA} {Li~I} Absorption}

Lithium is fused at a lower temperature
than hydrogen (2.5$\times$10$^6$~K), resulting in a somewhat lower fusing mass limit
(0.065~{M$_{\odot}$}; \citealt{1997ApJ...482..442B}).
Because the interiors of low mass stars and brown dwarfs are fully convective at early ages, an
object with a mass above this limit will fully deplete its initial reservoir of lithium within $\sim$200~Myr.
Hence, any system older than this
which exhibits {Li~I} absorption has a mass less than 0.065~{M$_{\odot}$}
and is therefore a brown dwarf (e.g., \citealt{1992ApJ...389L..83R}).  With a mass limit, one can use evolutionary models to determine an age limit.

In the age range $\sim$10--200~Myr, the degree of lithium depletion in low-mass stars
and brown dwarfs
is itself mass-dependent, occurring earlier in more massive stars which first achieve
the necessary core temperatures.  Hence, over this range the age of an individual source can be
precisely constrained if its mass is known.  A more practical approach, however, is to
ascertain the age of a group of coeval low-mass objects based on which sources do or do not
exhibit {Li~I} absorption; this is the lithium depletion boundary (LDB) technique.  Different
evolutionary models yield remarkably similar predictions for the location of the LDB over a broad
range of ages \citep{2004ApJ...604..272B}, and the boundary itself is readily identifiable in color-magnitude
diagrams.  As such, this technique has been used to age-date several nearby young clusters and
associations (e.g., \citealt{1998ApJ...499L.199S,1999ApJ...522L..53B,2005MNRAS.358...13J,2008arXiv0808.3584M}).  LDB studies have also provided independent
confirmation of other cluster-dating methods such as isochrone fitting \citep{2005MNRAS.358...13J}.
A variant of the LDB technique for coeval binary systems has been proposed by \citet{2005ApJ...634..616L}, in which a system that exhibits {Li~I} absorption in the secondary but not
in the primary can have both lower and upper age limits assigned to it (note that the presence/absence of {Li~I} in both components simply sets a single upper/lower age limit).  This technique requires resolved optical spectroscopy of both components and can be pursued only for
a few (rare) wide low-mass pairs (e.g., Burgasser et al., in prep.)
or using high spatial-resolution spectroscopy (e.g. \citealt{2006A&A...456..253M}).
No single brown dwarf pair straddling the LDB has yet been identified.

Despite its utility, the detection of LI~I absorption in brown dwarf spectra has
limitations.  The 6708~{\AA} line lies in an relatively faint spectral region for cool L-type dwarfs, so sensitive spectral observations on large telescopes are typically required to
detect (or convincingly rule out) this feature.  For optically-brighter
M-type brown dwarfs,
high-resolution observations are generally required to distinguish {Li~I} absorption from overlapping
molecular absorption features.  Young brown dwarfs ($t \lesssim$ 50~Myr)
with low surface gravities show weakened alkali line absorption (see $\S$~5),
including {Li~I}, making it again necessary to obtain sensitive, high-resolution
observations \citep{2008arXiv0808.3153K}.  For brown dwarfs cooler than $\sim$1500~K
(i.e., the T dwarfs), lithium is chemically
depleted in the photosphere through its conversion to LiCl, LiF or LiOH \citep{1999ApJ...519..793L}.  As such, practical age constraints using {Li~I}
can only be made for systems younger than $\sim$2~Gyr.

\section{Mass standards}

One way of breaking the mass/age/luminosity degeneracy for an individual brown dwarf is to explicitly
measure its mass.  This is feasible for sufficiently tight brown dwarf binaries for which
radial velocity (RV) and/or astrometric orbits can be measured.
Of the $\sim$100 very low mass (M$_1$,M$_2$ $\leq$ 0.1~{M$_{\odot}$}) binary systems
now known, only a handful have sufficiently short periods that
large portions of their RV orbits (e.g., \citealt{2007ApJ...666L.113J,2008ApJ...678L.125B}), astrometric orbits (e.g., \citealt{2001ApJ...560..390L,2004A&A...423..341B, 2008arXiv0807.0238L,2008arXiv0807.2450D}), or both \citep{2004ApJ...615..958Z,2006Natur.440..311S} have been measured.  With measurable total system masses or mass functions, individual masses can be
estimated from relative photometry or directly
determined from recoil motion in both components (e.g., \citealt{2006Natur.440..311S}).  The individual masses and component
luminosities can then be compared to evolutionary models to determine ages.

\citet{2008arXiv0807.0238L} have suggested that such ``mass standards'' provide more precise constraints on the
physical properties (including ages) of brown dwarfs as compared to ``age standards'', namely companions to
main sequence stars.  Orbital masses can currently be constrained to
roughly 5--10\% precision, translating into 10--20\% uncertainties in ages based
on evolutionary tracks (versus 50--100\% for main-sequence stars).
More importantly, brown dwarf binaries with
mass measurements and independent
age determinations---i.e., companions to age-dated stars and cluster members---can provide
specific tests of the evolutionary models themselves.  Further details are provided
in the contribution by T.\ Dupuy.

\section{Surface gravity}

Only 10--20\% of brown dwarfs are found to be multiple (e.g., \citealt{2006ApJS..166..585B}) and few of these are suitable for orbital mass measurements.  A proxy for mass
is surface gravity, which can be determined directly from a brown dwarf's spectrum.
For {T$_{eff}$} $\lesssim$ 2500~K and ages greater than $\sim$50~Myr, evolutionary
models predict that brown dwarf surface gravities ($g \propto M/R^2$) are roughly proportional to mass due to
near-constant radii (roughly equal to Jupiter's radius).
Surface gravity is also proportional to photospheric pressure ($P_{ph} \propto g/{\kappa}_R$, where $\kappa_R$ is the Rosseland mean opacity),  which in turn influences the chemistry, line broadening and (in some cases) opacities of absorbing species in the photosphere.
Hence, ``gravity-sensitive'' features in a brown dwarf's spectrum can be used to
infer its mass and, through evolutionary models, its age.

Examples of gravity-sensitive features include the optical and near-infrared VO bands and alkali lines in late-type M and L dwarfs, all of which evolve considerably between field dwarfs ({$\log{g}$} $\approx$ 5 cgs),
young cluster dwarfs ({$\log{g}$} $\approx$ 3--4) and giant stars ({$\log{g}$} $\approx$ 0; e.g., \citealt{1999ApJ...525..466L}).
Enhanced VO absorption and weakened alkali line absorption is a characteristic trait
of young brown dwarf spectra
(e.g., \citealt{2003ApJ...593.1074G,2007ApJ...657..511A}).
Alkali features in particular are useful for cooler brown dwarfs
as VO condenses out of the photosphere.
Quantitative analyses of these features have been used to distinguish ``young'' ($\gtrsim$10~Myr) from ``very young'' sources thus far ($\lesssim$5~Myr; e.g., \citealt{2007AJ....133..439C}).
More robust metrics await larger and more fully-characterized samples.

Another important surface gravity diagnostic is collision-induced H$_2$
absorption, a smooth opacity source spanning a
broad swath of the infrared (e.g., \citealt{1969ApJ...156..989L,1997A&A...324..185B}).
H$_2$ absorption is weakened in the low-pressure atmospheres of
young cluster brown dwarfs, resulting in reddened near-infrared spectral energy
distributions and colors; in particular, a characteristic, triangular-shaped H-band (1.7~{$\mu$m})
flux peak (e.g., \citealt{2001MNRAS.326..695L,2006ApJ...639.1120K}).
The proximity of many young and reddened brown dwarfs ($<$ 100~pc) rules out ISM absorption as a primary source for this reddening.   \citet{2008MNRAS.385.1771J} have exploited this trend by using a
proper-motion selected sample of nearby young cluster candidate members to infer an
age/color/luminosity relation for brown dwarfs younger than 0.7~Gyr, with a stated accuracy
of $\pm$0.2 dex in log(age), or about 60\% fractional uncertainty.
Kinematically older low-mass stars and brown dwarfs in the Galactic disk
(e.g., \citealt{2008arXiv0809.3008F}) and halo populations (e.g., \citealt{2003ApJ...592.1186B})
exhibit unusually blue near-infrared colors due to enhanced H$_2$ absorption.
However, differences in metallicities and condensate cloud properties can
muddle surface gravity determinations in these sources by modulating the
photospheric pressure through opacity effects (changing $\kappa_R$; e.g.,
\citealt{2000ApJ...535..965L,2008ApJ...686..528L}).

The use of H$_2$ absorption as a surface gravity indicator is particularly useful
for T dwarfs, as H$_2$ dominates the $K$-band (2.1~{$\mu$m}) opacity
and significantly influences near-infrared colors (e.g., \citealt{2002ApJ...564..421B,2004AJ....127.3553K}).
Several groups now employ this feature to estimate the atmospheric properties of individual T dwarfs
(e.g.,  \citealt{2006ApJ...639.1095B,2007MNRAS.381.1400W,2008MNRAS.tmp.1183B}), typically through the use of
spectral indices that separately sample surface gravity (e.g., the $K$-band) and temperature variations (e.g., {H$_2$O} or {CH$_4$} bands).  These indices are compared
to atmospheric models calibrated by one or more benchmarks (e.g., a companion to a precisely age-dated star), and evolutionary models are used to determine individual masses and ages.
Typical uncertainties of {$\log{g}$} $\sim$0.3~dex translate into 50-100\% uncertainties in age, comparable to uncertainties for main
sequence stars.
Again, variations in metallicity can mimic variations in surface gravity, although a third
diagnostic such as luminosity can break this degeneracy (e.g., \citealt{2007ApJ...658..617B}).
As atmosphere models improve in fidelity, parameters are increasingly inferred
from direct fits to spectral data, with comparable uncertainties
(e.g., \citealt{2007ApJ...656.1136S,2008ApJ...678.1372C}).

\section{Kinematics}

When a sufficiently large enough sample of brown dwarfs is assembled, one can
apply standard kinematic analyses, building from the assumption
that gravitational perturbations lead to increased velocity
dispersions with time (e.g., \citealt{1953ApJ...118..106S}; see
contribution by B.\ Nordstr\"om). Velocity dispersions are typically characterized
by a time-dependent power-law form, i.e., $\sigma \propto (1+t/\tau)^{\alpha}$ (e.g., \citealt{1977A&A....60..263W,2002MNRAS.337..731H}). Other statistics, such as Galactic scale height, can also be tied to age through kinematic simulations (e.g., \citealt{2008AJ....135..785W}) to calibrate
secondary age diagnostics such as magnetic activity (see contribution by A.\ West).

Samples of very low mass stars and
brown dwarfs have only recently become large enough that kinematic studies are feasible.
The largest samples (over 800 sources) have been based on proper motion measurements (e.g., \citealt{2007AJ....133.2258S,2008MNRAS.390.1517C,2008arXiv0809.3008F}).  For field dwarfs,
dispersion in tangential velocities for both magnitude- and volume-limited samples indicate a mean age
in the range 2--8~Gyr, largely invariant with spectral type.  This is consistent with the ages of
more massive field stars and population synthesis models
(e.g., \citealt{2004ApJS..155..191B}).  However, when field samples are broken down by color \citep{2008arXiv0809.3008F} or presence of magnetic activity \citep{2007AJ....133.2258S}, distinct age groupings are inferred, indicating that both very old (i.e., thick disk or halo) and very young (i.e., thin disk or young association) brown dwarf populations coexist in the immediate vicinity of the Sun.  Indeed, one of the major results from these studies is the identification of widely-dispersed brown dwarf members of nearby, young moving groups such as the Hyades (e.g., \citealt{2007MNRAS.378L..24B}).

With only two dimensions of motion measured, proper motion samples may produce biased dispersion measurements depending on the area of sky covered by a sample.  Full 3D velocities require
RV measurements which are more expensive and have thus far been obtained only for a small fraction of the known brown dwarf population
(e.g., \citealt{2006AJ....132..663B,2007ApJ...666.1198B}).
A seminal study by \citet{2007ApJ...666.1205Z} of
21 nearby, late-type dwarfs with parallax, proper motion and RV measurements found
considerably smaller 3D velocity dispersions for L and T dwarfs than GKM stars, suggesting that local brown dwarfs
are young ($t \sim$ 0.5--4~Gyr).  The discrepancy between this result
and the proper motion studies may be attributable to small number statistics and/or contamination
by young moving groups; $\sim$40\% of the brown dwarfs in the \citet{2007ApJ...666.1205Z} sample appear
to be associated with the Hyades.  Resolving this discrepancy requires
larger RV samples, which has the side benefit of potentially uncovering RV variables
that can be used as mass standards (e.g., \citealt{2006AJ....132..663B}).

\section{Improvements and Future Work}

With several methods for age-dating brown dwarfs over a broad range of phase space now
available,
opportunities to use these objects as chronometers for various
Galactic studies look to be increasingly promising; e.g., age-dating
planetary systems, examining cluster age spreads,
testing Galactic disk dynamical heating models, and direct measures of the
substellar mass function and birthrate in the field and other populations.  However, there are areas where improvements
in uncertainties are needed and basic assumptions tested.  Surface gravity determinations
in particular require better constraints, since these
enable age-dating of individual sources.  In the short term, improvements
in spectral models can help in this endeavor; however, a sufficiently sampled grid of
benchmark sources
may obviate the need for models entirely.
Benchmarks should increasingly arise from mass standards, for which
age constraints are more precise; these additionally provide
necessary tests of evolutionary models upon which most of the age-dating
techniques hinge.    Improved angular resolution and sensitivity
with $JWST$ and the next generation of large ($>25$m) telescopes
will increase resolved binary sample sizes by expanding the volume in which they
can be found and monitored.  These facilities will also aid searches for
substellar cluster members in old open field and globular clusters and, perhaps more
importantly, mass standards in these clusters to facilitate more stringent tests of evolutionary models
(T.\ Dupuy, priv.\ comm.).
Finally, a larger, more complete sample of brown dwarfs with precise RV measurements
will both improve our statistical constraints on the age of the local brown dwarf
population (and subpopulations) while additionally aiding in the search for mass standards.

\bigskip

\noindent The author thanks T.\ Dupuy, S.\ Leggett, M.\ Liu, \& A.\ West for helpful comments.

Brown Dwarfs as Galactic Chronometers


Neutron Diffraction Study of anomalous negative thermal expansion in non-superconducting PrFe1-xRuxAsO

Abstract: Neutron powder diffraction has been used to investigate the structural and magnetic behavior of the isoelectronically doped Fe pnictide material PrFe1-xRuxAsO. Substitution of Ru for Fe suppresses the structural and magnetic phase transitions that occur in the undoped compound PrFeAsO. Contrary to the behavior usually observed in 1111 pnictide materials, the suppression of both the structural and magnetic transitions does not result in the emergence of superconductivity or any other new ground state. Interestingly, PrFeAsO itself shows an unusual negative thermal expansion (NTE) along the c-axis, from 60K down to at least 4K; this does not occur in superconducting samples such as those formed by doping with fluorine on the oxygen site. We find that NTE is present for all concentrations of PrFe1-xRuxAsO with x ranging from 0.05 to 0.75. These results suggest that the absence of superconductivity in these materials could be related to the presence of NTE.

\section{\label{sec:level1}Introduction}

Since the discovery of Fe-based superconductors in 2008\cite{kamihara08, lynn09} there has been enormous interest in obtaining a detailed understanding of the properties of closely related materials.  Studies of substitutionally altered or doped materials have proved particularly valuable and have led to the discovery of several superconducting families.  Here we report new results for the isoelectronically doped 1111 pnictide material PrFe$_{1-x}$Ru$_x$AsO.

The 1111 FeAs based materials (typified by LaFeAsO) share several common characteristics with most of the other Fe-based superconductors. As the temperature is lowered the undoped parent material does not reach a superconducting ground state, but undergoes a structural transition, followed or accompanied by an antiferromagnetic spin density wave (SDW) transition of the moments associated with the Fe atoms and sometimes antiferromagnetic (AFM) ordering of the rare earth moments\cite{lynn09}.  For superconductivity to emerge in the material, all of these transitions must be eliminated or suppressed down to a sufficiently low temperature\cite{canfield09, sanna10}.  This can be done by doping with the appropriate element, e.g. F on the O-site in CeFeAsO\cite{zhao08}, Co on the Fe-site in LaFeAsO\cite{sefat08, wang09}, etc.  In certain members of the Fe oxypnictides, superconductivity can also be induced by simply applying pressure externally, e.g. LaFeAsO under 120kbar for a T$_C$ of 21K\cite{okada08} and BaFe$_2$As$_2$ under 40kbar for a T$_C$ of 29K\cite{alireza09}.  The application of pressure does not always produce this result, for example CeFeAsO achieved no superconductivity under pressures up to 500kbar\cite{zocco11}.

There are situations where AFM ordering and superconductivity can coexist, for example in 122 materials such as in BaFe$_{2-x}$Co$_x$As$_2$\cite{nandi10}.  Although this occurs in several 122 materials \cite{lumsden10} the general presumption that superconductivity and AFM ordering are competing ground states remains valid\cite{kim11, tanabe11}.  Notably coexistence is much less common in 1111 materials such as LaFeAsO, where long range AFM ordering is usually completely destroyed before superconductivity can emerge, although there are some exceptions such as SmFeAsO$_x$F$_{1-x}$\cite{drew09}. The suppression of the structural and magnetic transitions can also lead to new non-superconducting ground states, for example ferromagnetism\cite{luo10}.  It has been concluded that suppressing the structural and magnetic transitions is usually necessary, but not sufficient, for inducing superconductivity\cite{canfield09}.

\begin{figure}
\includegraphics[width=65mm]{figure1}
\caption{The PrFe$_{1-x}$Ru$_x$AsO room temperature tetragonal unit cell showing the rare-earth oxygen (top, bottom) and Fe-As (middle) layers, the Fe-As pyramid and the Fe-As-Fe angle $\phi$.  Note that in the orthorhombic structure there are two inequivalent values for $\phi$.}
\end{figure}

The unit cell of the 1111 parent compound PrFeAsO is illustrated in Figure 1.  Most of the materials with the highest reported superconducting T$_C$'s are doped on either site of the rare-earth-oxygen layer\cite{lynn09, tanabe11}.  Superconducting T$_C$'s are generally lower for materials doped on the Fe-As layer, possibly due to undesirable lattice distortions and disorder in the conducting Fe-As layer\cite{sharma10, delacruz10}.  To better understand this, isoelectronic substitution on the Fe site has proved very instructive since unlike carrier doping it enables the study of materials without directly altering the Fermi level\cite{mcguire09a}.  Moreover, isoelectronic doping with Ru can lead to superconductivity in certain Fe-oxypnictides, such as the 122 materials, BaFe$_{2-x}$Ru$_x$As$_2$\cite{sharma10} and SrFe$_{2-x}$Ru$_x$As$_2$\cite{qi09}.  In contrast, however, Ru/Fe substitution in 1111 materials has not been observed to lead to superconductivity even after the complete suppression of the structural and magnetic transitions: the studied systems include La(Fe,Ru)AsO\cite{sanna10}, Pr(Fe,Ru)AsO\cite{mcguire09a}, Sm(Fe,Ru)AsO and Gd(Fe,Ru)AsO\cite{pal10}.

McGuire \textit{et al.} showed that Ru doping in PrFeAsO does not induce superconductivity down to at least T = 2K\cite{mcguire09a}.  Even when the structural and magnetic transitions were suppressed completely there was no indication of any new ground state\cite{mcguire09a}.  Bulk and transport measurements on PrFe$_{1-x}$Ru$_x$AsO showed hints of anomalies associated with a possible phase transition for concentrations up to x = 0.5 or 0.6\cite{mcguire09a}.  This is a rather high level when compared with carrier doped Fe-pnictides, where the transitions are suppressed with less than 10\% doping in most cases.  Thus, Ru-doping provides a wide stoichiometric window for studying the progressive suppression of the structural and magnetic transitions.   Ru-doping studies have been also carried out on F-doped superconducting samples, such as LaFe$_{1-x}$Ru$_x$AsO$_{0.89}$F$_{0.11}$\cite{satomi10} and SmFe$_{1-x}$Ru$_x$AsO$_{0.85}$F$_{0.15}$\cite{tropeano10}.  The general results show that the superconductivity in these samples is not sensitive to isoelectronic doping except at very high levels of the order of x = 0.3 and above.

A phenomenon of particular interest in the  parent compound PrFeAsO is the observation of negative thermal expansion (NTE) along the c-axis, where the sample gradually expands along the c direction as it is cooled below 60K\cite{kimber08}.  Notably such NTE is absent in superconducting F-doped PrFeAsO\cite{kimber08}.  Recent results have shown that NTE also occurs in NpFeAsO were it is associated with magnetic ordering of the Np ions\cite{klimczuk12}.  A different NTE behavior has been observed in the superconducting samples of Co doped BaFe$_2$As$_2$ where a sudden onset of NTE along the c axis appears at T$_C$.  This feature is absent in the non-superconducting parent BaFe$_2$As$_2$\cite{daluz09}.  It seems worthwhile to explore whether or not there is a systematic correlation between the NTE and the presence or absence of superconductivity. Despite voluminous reports of the dependence of lattice constants on stoichiometry in Fe-As based materials there has been surprisingly little published data showing the temperature dependence of the lattice constants, especially the c-axis.

In this paper we report new results on Ru doped PrFe$_{1-x}$Ru$_x$AsO, , including detailed studies of the temperature dependence of lattice parameters.  As mentioned above undoped PrFeAsO shows NTE.  As discussed below, we observe that Ru-doping in PrFeAsO suppresses the structural and magnetic phase transitions without leading to superconductivity; however NTE along the c axis persists for all doping levels up to at least 75\% Ru on the Fe site. This paper is organized as follows: Sample synthesis and preliminary characterization is described in section II, the results of neutron scattering experiments are described and discussed in section III, with conclusions in section IV.

\section{\label{sec:level2}Sample synthesis and characterization}

\begin{figure}
\includegraphics[width=70mm]{figure2}
\caption{(a): Heat capacity vs T for PrFe$_{1-x}$Ru$_x$AsO, the peak at 14K indicates the AFM ordering of Pr in PrFeAsO. No peak is observed for x $\ge$ 0.1;
(b):$\chi$$^{-1}$ vs T, the inset shows a closer look at the lower temperature range for $\chi$ vs T.}
\end{figure}

Polycrystalline samples of PrFe$_{1-x}$Ru$_x$AsO with x = 0.05, 0.1, 0.33, 0.4, 0.5, 0.6 and 0.75 were synthesized using methods reported previously by McGuire \textit{et al.}\cite{mcguire09a}.  In each case, a stoichiometric mixture of powders of PrAs, Fe$_2$O$_3$, RuO$_2$, Fe and Ru was crushed and mixed thoroughly in a glove box, pressed into pellets and received two firings at 1200-1250$^o$C.  Transport property data and results of electronic density of state calculations for the same series of materials have been reported previously\cite{mcguire09a}.
The undoped parent compound PrFeAsO undergoes a tetragonal-orthorhombic structural transition at T $\sim$ 150 K, followed by a SDW transition at T $\sim$ 140 K associated with the AFM ordering of Fe, and then finally AFM ordering of the Pr moments at 14 K\cite{mcguire09b}.  Previously published transport measurements for PrFe$_{1-x}$Ru$_x$AsO have shown evidence for the suppression of long range order of the Pr moments for x $\ge$ 0.1 and structural/SDW transitions for x ${ >  0.5}$\cite{mcguire09a}.

Prior to the neutron powder diffraction experiments, we have carried out heat capacity and magnetization measurements.  As seen in Fig. 2(a) and (b), the AFM ordering of Pr at 14K in the x=0 material is indicated by a peak in the heat capacity and a kink in the magnetization.  These features disappear for x $\ge$ 0.1.

\section{\label{sec:level3}Experimental results and discussion}

\begin{figure}
\includegraphics[width=90mm]{figure3}
\caption{(Upper): Fully refined neutron powder diffraction pattern collected at HB2A of HFIR for PrFe$_{1-x}$Ru$_x$AsO at 4K; (Lower): 220$_N$ peak for x = 0.1, 0.33 and 0.75 at T = 200K (red circles) and 4K (black squares).  Peak splitting is clearly visible only for x=0.1.}
\end{figure}

\begin{figure}
\includegraphics[width=90mm]{figure4}
\caption{The ratio $\chi^2_{tetra}$ given by the Rietveld refinement with the P4/nmm (tetragonal, space group 129) structure, to $\chi^2_{ortho}$ obtained from refining to the Cmma (orthorhombic, space group 67) structure. Any value larger than 1.0 indicates the Cmma model provides a better fit than P4/nmm.  The represents data collected from both HB-2A and POWGEN.}
\end{figure}

\begin{figure}
\includegraphics[width=80mm]{figure5a}
\includegraphics[width=83mm]{figure5b}
\caption{(a): The intensity of the 101$_M$ peak as a function of temperature.  The AFM ordering of both Pr and Fe are visible in x = 0.05 but not x = 0.1.  The intensity of the 101$_M$ peak reaches a domelike maximum around 60K for x = 0.05.  Solid dots represent single data points collected from counting at the center peak position; hollow dots represent normalized data from fits of entire magnetic peaks.  The lines and shadings are guides to the eye;
(b): Intensity of the 100$_M$ peak as a function of temperature, showing the size difference of the Pr moment between x = 0.05 and 0.1.}
\end{figure}
\begin{figure}
\includegraphics[width=85mm]{figure6}
\caption{Neutron powder diffraction data collected from POWGEN refined to the orthorhombic Cmma structure for x = 0.33 at T = 12K.}
\end{figure}

\begin{table*}
\caption{Structural parameters and $\chi$$^2$ values from Rietveld refinement of NPD data from HB2A (for Ru\% = 10 only) \\and POWGEN, at T = 12K.
(Space group: Cmma (\#67) for 10\% and 33\%, P4/nmm (\#129) otherwise.)}

\begin{ruledtabular}
\begin{tabular}{ccccccc}

\textrm
Ru\%  & a/b (\AA)				&c (\AA)		&z$_{As}$		&z$_{Pr}$		&$\phi$	&$\chi$ $^2$\\
\colrule
10	&3.96865(8)/3.98509(8)		&8.5622(2)   &0.6550(5)  	&0.1402(7) 	&71.89(9)/72.23(9)	&3.34	\\
33	&4.00063(3)/4.00733(3)		&8.5154(1)	&0.6570(3)	&0.1389(5)	&71.95(7)/72.09(7)	&3.64	\\
40	&4.01162(4)				&8.4911(2)	&0.6582(5)	&0.1378(6)	&71.97(7)	&3.51	\\
50	&4.02491(6)				&8.4613(2)	&0.6626(2)	&0.1328(5)	&71.43(7)	&2.82	\\
60	&4.03351(7)				&8.4364(3)	&0.6603(6)	&0.1333(9)	&71.93(6)	&4.04	\\
75	&4.05088(4)				&8.3980(1)	&0.6587(4)	&0.1366(6)	&72.41(7)	&3.26	\\

\end{tabular}
\end{ruledtabular}
\end{table*}
	We performed neutron powder diffraction on PrFe$_{1-x}$Ru$_x$AsO using the neutron powder diffractometer HB-2A (for x = 0.1, 0.33 and 0.75) and the fixed incident-energy triple-axis HB-1A (for x = 0.05 and 0.1) of HFIR at ORNL to study the structural and magnetic transitions respectively.  Additional measurements were performed at the high-resolution powder diffractometer POWGEN of SNS at ORNL to study the effect of Ru doping on the anomalous NTE and also to have a closer look at the suppression of the structural transition (for x = 0.33, 0.4, 0.5, 0.6 and 0.75).

	The detailed structures were fitted via Rietveld refinement using the FULLPROF program\cite{caravajal93}.  Fig. 3 (upper) presents the refined neutron powder diffraction data collected at HB-2A for the x = 0.1 sample at the base temperature of 4K. 	Fig. 3 (lower) shows a closer look at the evolution of the 220 nuclear peak for x = 0.1, 0.33 and 0.75. The structural transition is evident in the x = 0.1 sample as the peak splitting in 220$_N$ is easily discernible; this is much less obvious for samples with x $\ge$ 0.33.  Fig. 4 shows the temperature dependence of the relative goodness of fits obtained by fitting the full powder patterns to each of the two different known structural models (tetragonal P4/nmm vs. orthorhombic Cmma.)   It is apparent that the orthorhombic model fits better than the tetragonal model for x = 0.1 at T $\le$ 120K and for x = 0.33 at T $\le$ 75K.  The inferred structural phase transition temperature in the x = 0.33 sample is consistent with the anomaly near 80K observed in the resistivity data\cite{mcguire09a}.  The diffraction measurements for higher Ru concentrations provide no evidence for any structural transition, suggesting that it is suppressed for Ru doping between 33-40\%.

	Interestingly for other isoelectronically doped 1111 materials the structural phase transition has been observed to persist to dopings well above the 10\% level;  for example in CeFeAs$_{1-x}$P$_x$O  both the transition to an orthorhombic structure and antiferromagnetic order are present up to x = 0.4\cite{luo10}.   In contrast, a much smaller level of carrier doping is needed to affect the transitions. In RFeAsO$_{1-x}$F$_x$ all structural and magnetic transitions were suppressed with less than 10\% doping\cite{kamihara08, rotundu09, ren08a, ren08b, ren08c}.  In fact the doping level resulting in the highest observed superconducting transition temperature occurs at x = 0.1 for R = La\cite{kamihara08}, 0.11 for Pr\cite{rotundu09, ren08a}, 0.11 for Nd\cite{ren08b} and 0.1 for Sm\cite{ren08c}, with superconductivity usually disappearing around x = 0.2\cite{ren08a, ren08b, ren08c}.

	  We now turn to a detailed examination of the magnetic transition that is observed in PrFe$_{1-x}$Ru$_x$AsO with 5\% Ru doping.  In Fig. 5 we present temperature dependent measurements of the intensities of the 101$_M$ and 100$_M$ magnetic peaks measured using HB-1A for doping levels x = 0.05 and 0.1.  The 101$_M$ peak shows contributions from both Pr and Fe moments while the 100$_M$ peak is sensitive only to the Pr moments.  It is apparent in Fig. 5(a) and (b) that magnetic ordering of both Pr and Fe exists at x = 0.05 but is suppressed at x = 0.1.  Fe orders antiferromagnetically in the x = 0.05 sample at 120K which may be compared to 137K\cite{mcguire09b} observed in the undoped PrFeAsO; no AFM order is visible in x = 0.1.    Evidently in PrFe$_{1-x}$Ru$_x$AsO, the magnetic transition is more sensitive to Ru doping than is the structural transition.  In any case neither transition is observed for x $\ge$ 0.4.

	The temperature dependence of the magnetic intensity of the 101$_M$ peak is particularly interesting. As seen in Fig 5(a), in the x = 0.05 sample the intensity of the 101$_M$ peak reaches a dome-like maximum at T $\sim$ 60 K, followed by a modest decrease as the temperature is lowered, increasing again below T $\sim$ 14 K, the AFM ordering temperature for Pr.    This feature has also been observed in undoped PrFeAsO at T $\sim$ 60K\cite{kimber08}.  Critical fluctuations in the Pr magnetic subsystem may contribute to this phenomenon\cite{kimber08}.

	To examine the temperature dependence of the crystal structure and lattice constants and in particular the behavior of the NTE we performed higher resolution neutron powder diffraction at POWGEN for PrFe$_{1-x}$Ru$_x$AsO with x = 0.33, 0.4, 0.5, 0.6 and 0.75.

	Fig. 6 illustrates a typical refinement using the powder pattern for x = 0.33.  	Table I shows a summary of the lattice and crystallographic parameters obtained at 12 K for several different concentrations of Ru. The angle $\phi$ connecting near-neighbor Fe atoms via an As atom (see Fig. 1) is also included.  To within the resolution of the experiment for each concentration the crystallographic parameters $z_{As}$, $z_{Pr}$, and the angle $\phi$ have only a very small temperature dependence below 200 K. It should be noted however that the precision with which these parameters are determined is an order of magnitude coarser than the precision with which the lattice constants are determined.  This is expected since the lattice constants are determined precisely by the peak positions alone whereas the other parameters are refined by fitting peak intensities.  On the other hand there is a clearly observable dependence of the lattice constants on the temperature over the entire measured range.
\begin{figure}
\includegraphics[width=80mm]{figure7a}
\includegraphics[width=80mm]{figure7b}
\caption{(Upper): Doping dependence of a and c at 12 K; (Lower): Temperature dependence of a- and c-axis plotted together, showing the contrasting thermal behavior between the two.  The red line is the Gruneisen model prediction for the c axis based on fits to the data for T $\ge$ 40K as described in the text.}
\end{figure}
\begin{figure}

\includegraphics[width=80mm]{figure8}
\caption{Unit cell volume as a function of temperature.  The quantity $\Delta$vol$/vol_{140K}$ is defined as (vol(T)-vol(T=140K))/vol(T=140K). The dotted red line is Gruneisen prediction for $x = 0.33$ fitted using the data for T $\ge$ 40K.}
\end{figure}

	Fig. 7 (upper) shows the contrast in the behavior of the a- and c-axes' doping dependence at 12K.  The trends are consistent with the room temperature x-ray data reported by McGuire \textit{et al.}\cite{mcguire09a}, which they noted may be related to both the larger radius of Ru compared to Fe as well as constraints on Ru-As and Fe-As bond distances.  Fig. 7 (lower) shows the temperature dependence of the lattice parameters for the concentration x = 0.4. The dashed red line is the extrapolated behavior for the c-axis based on a Gruneisen fit to the volume defined by c$^3$ for the data at 40K and above.  This is a representative concentration as the others show similar behavior. At low temperatures (i.e. for T $\le$ 50K), the lattice parameters deviate from the Gruneisen model fit in opposite directions, somewhat conserving the overall cell volume.  This is illustrated in Fig. 8, which shows the temperature dependence of the unit cell volume. A Gruneisen model fits the entire cell volume better than it fits either the a- or c-axis alone.  This behavior is not uncommon among materials that are anisotropic, where expansion along one direction is often accompanied by contraction along another in order to preserve the overall volume; graphite is one well known example\cite{bailey70}.  Despite this compensation a modest NTE of the unit cell volume is still visible in all samples for T $\le$ 20K.

\begin{figure}
\includegraphics[width=80mm]{figure9a}
\includegraphics[width=80mm]{figure9b}
\includegraphics[width=80mm]{figure9c}
\caption{ Detailed temperature dependence of of the lattice constants for several values of Ru doping. (Upper): The quantity $\Delta$c$/c_{140K}$ is defined as (c(T)-c(T=140K))/c(T=140K).  NTE is clearly visible in all samples for T $\le$ 50K.  The solid line is a representative Gruneisen model fit for x = 0.33 also shown in Fig. 8. (Middle): A similar plot, showing the relative change in a- and b-axis with respect to lattice parameter a measured at 140K, which is above the structural transition for x = 0.1 and 0.33;(Lower): A similar plot for a-axis for x = 0.4, 0.5 and 0.75, showing the lack of NTE at low temperatures.}
\end{figure}

	NTE in the c-axis is a prominent feature observed previously in the undoped parent compound: the c-axis reaches a minimum at T = 50K, and then expands as the system is further cooled\cite{kimber08}.  Kimber \textit{et al.} ascribed the presence of c-axis NTE in the non-superconducting parent compound to the effects of spin-lattice coupling\cite{kimber08}.  This was based on their observation that both magnetic order and NTE are absent in the superconducting PrFeAsO$_{1-x}$F$_{x}$, and the proximity of the temperature dependence of the NTE and that of the dome-like feature observed in the intensity of magnetic peaks associated with the ordering of the Fe moments\cite{kimber08}.    In this context, the observation of NTE in PrFe$_{1-x}$Ru$_x$AsO with $x > 0.1$ is surprising, since long range magnetic order is absent.  Fig. 9 shows the temperature dependence of the fractional change in the refined lattice parameters a, b and c for several concentrations.  As seen clearly in Fig. 9 (upper), the anomalous NTE along the c-axis is observed in all samples up to at least x = 0.75.  Figure 10 shows a phase diagram of PrFe$_{1-x}$Ru$_x$AsO with the region of c-axis NTE indicated.

In contrast the thermal behavior of the a-axis is almost linear at low temperatures, with the fractional change almost independent of the Ru concentration.   The existence of c-axis NTE in concentrations without magnetic order suggests that the origin of this behavior might involve more than simply spin-lattice coupling.  Although there is little published data on the c-axis temperature dependence of superconducting 1111 materials, the NTE is absent among those for which data is accessible, including PrFeAsO$_{0.85}$F$_{0.15}$\cite{kimber08}, and LaFeAsO$_{1-x}$F$_{x}$\cite{kondrat09, mcguire08}.  To date we know of no example in the literature of a superconducting 1111 compound that does show NTE.

	It is interesting to compare the low temperature thermal variation of the a- and c-axes with the composition dependence shown in Fig. 7.  In both cases the parameters vary in the opposite direction.  The reason for the opposite temperature dependence of the a- and c-axes remains an open question.  The combination of NTE in the c axis and near conservation of volume may lead to subtle but important changes in the effective magnetic and electronic interactions of the Fe-As layer both internally and with the rest of the sample.   Therefore, it can possibly be inferred that the NTE is linked to the lack of superconductivity (or any new ground state) in the Ru doped samples.

\begin{figure}
\includegraphics[width=85mm]{figure10}
\caption{Phase diagram of PrFe$_{1-x}$Ru$_x$AsO constructed with data reported in this article, showing temperatures for the tetragonal-orthorhombic transition (T$_{Structural}$), the magnetic Fe ordering (T$_{SDW}$), and the magnetic Pr ordering (T$_{Pr}$), in the measured Ru concentrations.  The region of c-axis NTE presents in all Ru concentrations below T = 50K.}
\end{figure}
\section{\label{sec:level4}Conclusion}
	We have found that magnetic order on both Fe and Pr sites in Pr-1111 is suppressed with 10\% of isoelectronic doping by Ru/Fe substitution, and the tetragonal-orthorhombic structural transition with 40\%.  We have studied the effect of Ru doping on the anomalous NTE along the c-axis and found that the c-axis NTE is present in all of our samples.  For all concentrations the lattice parameter along the c-axis reaches a minimum around 50K.  The NTE is extremely resilient to Ru doping and is measurable up to at least x = 0.75.  The c-axis NTE extends across the entire composition range studied: in compositions which remain tetragonal to the lowest temperatures, those which exhibit tetragonal-orthorhombic distortions with no long-range magnetic ordering, and those that undergo both structural and magnetic phase transitions.
 	Other anomalous behaviors observed include the near linear temperature expansion of a-axis at low temperature, which also has little discernible effect from doping.  The contrasting thermal behavior between a- and c-axis somewhat conserves the unit cell volume.  The temperature dependence of the unit cell volume agrees better with the Gruneisen prediction than the a- or c-axis alone.
	It would be interesting to investigate whether the apparent relationship of NTE to the lack of superconductivity is indeed a general rule for 1111 materials.  Such knowledge should provide important insights into the underlying mechanism for superconductivity in Fe based materials.

\section{\label{sec:level4}Acknowledgments}

This research was supported by the US Department of Energy, Office of Basic Energy Sciences.  Experiments were performed at SNS and HFIR, sponsored by the Scientific User Facilities Division.  V. O. G., A. H., and S. E. N. were supported by the Scientific User Facilities Division.  M. A. M. and D. M. were supported by U.S. Department of Energy, Office of Basic Energy Sciences, Materials Sciences and Engineering Division.  Y.Y. was supported by the Office of Basic Sciences, US Department of Energy, through the EPSCoR, Grant No. DE-FG02-08ER46528.


Radio Quiet AGN

Abstract: Active Galactic Nuclei are powered by accretion onto massive black holes. Although radio-quiet objects are not as spectacular sources of very high energy photons as radio-loud ones this class of objects also represents a challenge for modeling high energy processes close to a black hole. Both a hot optically thin plasma and a cooler optically thick accretion disk are usually thought to be present in the vicinity of a black hole although the details of the accretion flow are still under discussion. The role of the disk seems to decrease with a drop in the Eddington ratio: in sources like quasars and Narrow Line Seyfert 1 galaxies disk flow dominates while in Seyfert galaxies the disk retreats, and in sources like LINERS or Sgr A* a disk is most likely absent. Shocks and reconnections are possibly taking place in an inner hot flow and in the magnetic corona above the cold disk. Uncollimated outflow is also present and it may carry significant fraction of available mass and energy.

\section{Introduction}
Active Galactic Nuclei (hereafter AGN) are the most powerful persistent and compact
sources of radiation
in the Universe. As such, they are clearly of interest from the point of view of high
energy processes. AGN activity is caused by accretion of material onto the massive black hole residing at
the center of a host galaxy. Therefore, in those objects we can observe the behavior of the matter in a strong gravity field.

Although the basic engine is always the same ---  extraction of gravitational energy of the infalling material ---  AGN are quite an inhomogeneous class of objects. The most
important differentiating property is perhaps the radio loudness. A small fraction of
those sources are strong radio emitters and display spectacular jets and/or radio extended radio lobes.  The majority, however, forms a population of radio
quiet sources. The transition between the two populations is rather smooth, but
nevertheless it is convenient to introduce two separate classes of objects, with
a focus on strong relativistic jets as the main property of radio-loud sources.
The customary border is set at the 5 GHz radio to B band flux ratio
$ \log F_{5GHz}/\log F_B = 10$. The properties of the radio-loud objects are discussed by
R. Moderski and A. Celotti (these proceedings). Here I will concentrate on the
radio-quiet population and their relevance to high energy astrophysics.

\section{Basic facts}

AGN span a very broad range of luminosities. The brightest AGN (found among quasars)
have bolometric luminosities almost up to $\sim 10^{48}$ erg s$^{-1}$ \cite{woo}. They
are mostly found among the high redshift objects. Nowadays, numerous quasars
are found up to $z \sim 6$ \cite{highz} as a result of massive surveys (e.g.
SDSS \cite{sdss}).
Nearby Seyfert galaxies are a few orders of magnitude fainter. The lower limit for a
nuclear activity is
unspecified since the determination of the weak nuclear activity is
observationally difficult. However, it is now widely believed that all regular
galaxies with bulges contain a supermassive black hole and must show some level of
activity. In this sense we can also count SGR A* as an example of weak activity
(see S. Nayakshin, these proceedings),
with its occasional flares  reaching up to the level of $10^{35}$ erg s$^{-1}$.

Black hole mass estimates indicate much narrower range,
$\sim 10^6 -  10^{10} M_{\odot}$ \cite{wandel},
than the luminosity range
 which means that the Eddington ratio differs considerably between the sources.

\subsection{Classification}

Radio-quiet AGN are divided into several classes, with two parameters
most plausibly underlying this classification scheme: inclination angle
and the Eddington ratio.

Generally, all AGN are
divided into type 1 and type 2 objects. Type 1 objects have very broad emission lines while
type 2 objects show only narrow lines in their optical/UV spectra. It is generally accepted
that class 2 is simply an obscured
version of type 1 objects so we have no direct view of the nucleus in those sources.
This obscuration is due to the material predominantly located in the equatorial plane,
in a form of a dusty/molecular torus.
It can be noted, however, that there may be intrinsic type 2 sources among low
Eddington ratio AGN \cite{tran}. It only means that classification based on the width
of the hydrogen lines may not always represent well the characteristics of the
central engine. In further text we will discuss only type 1 AGN.

Another classification, mostly historical, divides type 1 radio quiet AGN into quasars,
Seyfert galaxies, Narrow Line Seyfert 1 galaxies and LINERS.
The transition between these classes (connected mostly with
the luminosity and the Eddington ratio) is smooth. Moreover, some quasars are also
classified as NLS1 galaxies,
and some LINERS perhaps are mostly starburst so the issue is confusing but we will
preserve these subgroups to indicate the luminosity class of the discussed sources. Also
for bright quasars the transition to the corresponding 'Narrow Line Type 1' class happens at
much higher width ($\sim 4000$ km s$^{-1}$ \cite{bachev}) than for Seyfert
galaxies ($\sim 2000$ km s$^{-1}$). This can be easily understood if the transition
happens at a fixed Eddington ratio and the difference
in black hole mass is taken into account (see formulae in \cite{BLR}).

\subsection{Broad band spectra}

The radiation spectrum of radio quiet AGN is very broad and span from radio to
gamma band. The spectrum is best studied for high luminosity sources. A schematic
view, representative for sources with high Eddington ratio (quasars and Narrow
Line Seyfert 1 galaxies; $L/L_{Edd} \sim 1$) is shown in Fig.~\ref{fig:spectrum} \cite{elvis94}.
The spectrum usually peaks in the far UV/soft X-ray range, not accessible to the
observations due to the Galactic extinction. However, the observational gap can be
partially filled  by combining low redshift and high redshift sources into a single
composite spectrum \cite{laor}. The extension of the spectra beyond
100 keV is not well studied and we return to this point in Sect.~\ref{sect:gamma_spectra}.

\begin{figure}
\psfig{figure=spectrum.ps,height=2.4 in}
\caption{Schematic view of the radiation spectrum of a bright radio quiet (continuous line) and radio-loud (dashed line) AGN.
\label{fig:spectrum}}
\end{figure}

The spectra of AGN with lower Eddington ratio like Seyfert 1 galaxies
($L/L_{Edd} \sim 0.01$) show less
pronounced far UV peak and relatively stronger X-ray emission, and the X-ray spectrum
is harder. The black holes accreting at the lowest rate ($L/L_{Edd} \sim 10^{-4}$)
like those in LINERS or LLAGN do not show any strong UV
emission while they are still relatively bright in X-rays. The broad band spectra of
such objects are difficult to determine since the optical spectrum is in this case
dominated by the host galaxy. Their X-ray spectra are as hard as in Seyfert 1 galaxies.

Practically all radio-quiet AGN are rather variable \cite{markowitz,optical}. The shortest
variability timescales (hundreds of seconds) are seen in X-ray band, in Seyfert galaxies.
Optical/UV emission varies in timescales of days in those sources while in quasars
monitoring over years is needed to see significant changes in optical band \cite{giveon}.

\section{Radio-quiet AGN as sources of high energy photons and particles}

Radio-quiet AGN are not as obviously attractive from the point of view of the high
energy physics as are radio-loud AGN. However, they also represent certain challenge
for such studies. In order to show that we will discuss now the
gamma-ray emission from these sources and the most plausible environment of its production.

\subsection{Extension of the AGN spectra into high energies}
\label{sect:gamma_spectra}

The direct signature of the importance of high energy processes in a given object is the presence of the gamma-ray emission.
Unfortunately,  spectral measurements going beyond 100 keV were performed only for a few sources,
and the results were uncertain. Fits to the OSSE composite
spectrum \cite{Zdzia-OSSE} gave the high
energy cut-off $120^{+220}_{-60}$ keV for Seyfert 1 and $130^{+220}_{-50}$ keV for Seyfert 2
galaxies (at such high energies the obscuration by the torus is relatively unimportant).
{\it Beppo-SAX} data, however, did not give such a strong constraints
(cut-off energy $231^{+ \infty}_{-168}$ keV for Seyfert 1 galaxies, no cut-off
for Seyfert 2 galaxies \cite{deluit}). Comparison of the model to the observed X-ray background
suggest that if the typical photon index is $\sim 1.9$ the spectrum extends up to $\sim 300$
keV \cite{background} while observations of high redshift quasars give results from 100 keV to 500 keV
and more, depending on the object and the adopted geometry \cite{sobo1,sobo2}.

Overall similarity of AGN and galactic black holes may suggest that sources with
rather hard spectra (photon index 1.9 or less) are thermal (electron temperature of order of 100 keV). Such a spectrum is well explained as a result
of the Comptonization of the soft photons by a predominantly thermal plasma with the electron
temperature around $10^9$ K. Ions may have much higher temperature, or may not be fully
thermalized, but we have no direct methods of estimating their properties. Efficient
thermalization is supposed to be achieved through synchrotron self-absorption \cite{ghisellini}.
Sources with steep spectra (photon index 2.0 and more) are likely to be significantly non-thermal, with a population of electrons having a power law distribution of energies, as expected in case of effective acceleration and ineffective thermalization.

Direct observation constraints should come from Astro-E2 and GLAST.

\subsection{accretion flow geometry}

\begin{figure}
\parbox{\textwidth}{
\parbox{0.5\textwidth}{
\epsfxsize=0.45\textwidth
\epsfbox[18 200 600 720]{fig2b.eps}}
\hfil
\parbox{0.5\textwidth}{
\epsfxsize=0.45\textwidth
\vskip -1.7 true cm
\epsfbox[18 200 400 0]{Geom2.ps}}
}
\caption{Left: Optical/UV composite spectrum of very bright quasars
(Francis  et al. 1991, continuous line), of moderately bright quasars
(Laor et al. 1997, long-dash line), together with an exemplary spectrum of
a simple Shakura-Sunyaev disk (short-dash line). Right: Schematic view
of the accretion flow in a bright AGN.
\label{fig:geom}}
\end{figure}

The broad band spectrum clearly shows that the accretion flow is
(at least) a two-phase medium. The profound optical/UV/soft X-ray bump (Big
Blue Bump) so characteristic for high $L/L_{Edd}$ sources is well
modeled as a thermal emission of a Keplerian, optically thick and geometrically
thin disk (see Fig.~\ref{fig:geom}). The hard X-ray emission must come from a
hot optically thin plasma in the vicinity of the disk. The IR emission is due
to reprocessing of a part of the radiation by circumnuclear dust, usually referred to
as a dusty-molecular torus. The presence of this torus is responsible for shielding
the central parts for highly inclined observers, as is the case for type 2 objects.
The torus is most probably clumpy, and quite possibly it is rather a kind of outflowing
dusty wind instead of a structure in the hydrostatic equilibrium.

The location and geometry of the hot plasma is uncertain. A plausible possibility is shown
in Fig.~\ref{fig:geom}. The disk is covered by magnetic loops emerging from its interior
\cite{galeev}. The
mechanism behind this is the magneto-rotational instability (MRI) operating in the disk interior.
This instability is responsible for the disk viscosity, and the corresponding
Shakura-Sunyaev parameter $\alpha$ is $\sim 0.01$, according to numerical simulations.
Occasionally, large loops emerge high above the disk surface \cite{MilStone00},
and the magnetic field
reconnection results in formation of a flare (hard X-ray flash). Therefore, a corona
similar to the solar corona forms above the disk. In the innermost part the
cold disk can be evaporated due to the interaction with such an active corona \cite{aga2000}.
A flow towards the black hole proceeds through a
hot plasma phase \cite{SLE,ichimaru,narayan}, and a (perhaps significant) fraction of the plasma can be expelled out
in the process. In strongly radio-loud objects this outflow takes a form of a
relativistic well-collimated jet but in radio-quiet objects the outflow is slower
and uncollimated. The mechanism of the outflow is unknown.

Within this picture, we would expect the formation of non-thermal plasma in the coronal loops above the disk while in the hot phase the electrons would be predominantly thermal. The relative role of the
two media is determined by the disk trucation radius which in turn is determined by $L/L_{Edd}$. At
$L/L_{Edd} \sim 0.1 - 0.5 $ or more the disk is supposed to extend down to the marginally stable orbit.

The geometry shown in Fig.~\ref{fig:geom} is not a unique possibility. Other accretion flow models include (i) the lamp-post model in which the disk always extends up to the marginally stable orbit and the high energy emission comes from a shock formed in the outflowing plasma and localized on the symmetry axis \cite{guy,miniutti} (ii) disk always extending down to the marginally stable orbit with magnetic flares above it \cite{galeev} (iii) accretion in the form of clumps of cold material embedded in a hot plasma \cite{collin96}. In those models the spectral differences between
high Eddington and low Eddington sources are less naturally explained.

\section{X-ray spectroscopy as a tool to study the plasma motion}

Given the lack of sufficient spatial resolution the key to the flow geometry
lies in X-ray spectroscopy. First extremely useful results came from Ginga \cite{ginga12} and ASCA \cite{tanaka95} data, now Chandra and XMM-Newton offer still higher quality data. A number of blueshifted absorption lines was identified which allowed to measure the outflow velocity of the material, and measured a number of emission line profiles, including the famous iron K$\alpha$ line.

\subsection{outflow}

Chandra and XMM-Newton observations confirmed the earlier findings that a partially ionized warm
absorber exists in many Seyfert 1 and NLS1 galaxies \cite{Blustin}. Outflow velocities range from
hundreds km s$^{-1}$ in Seyfert 1 \cite{kaspi,kaastra} galaxies to a fraction of light speed in
NLS1 \cite{chartas,pounds1211}. The distance of this outflowing material from the black hole is
difficult to assign. Simple arguments based on the assumption of a roughly Keplerian velocity
suggest that the slow outflow originates somewhere between the Broad Line Region and a Narrow
Line Region, i.e. at a few parsecs from the nucleus, and the fast outflow originates closer in,
at distances of order of hundreds of Schwarzschild radii. The amount of outflowing material may be
quite large - fast outflows may actually carry out considerable fraction of the inflowing
material \cite{kingpounds}. The estimates are difficult since a fraction of the material may be
completely ionized and leave no direct signature in the soft X-ray spectrum. However, this
material will scatter a fraction of the nuclear emission and may lead to modification of the
optical spectrum of an AGN through the irradiation of the outer disk \cite{uni,loska}.

\subsection{quasi-Keplerian motion}

\begin{figure}
\vskip -7.5 true cm
\parbox{\textwidth}{
\parbox{0.30\textwidth}{
\epsfxsize=0.25\textwidth
\vskip -1 true cm
\epsfbox[0 0 400 900]{plama.eps}}
\hfil
\parbox{0.30\textwidth}{
\epsfxsize=0.25\textwidth
\vskip -3 true cm
\epsfbox[0 0 300 900]{pierscienie.eps}}
\hfil
\parbox{0.30\textwidth}{
\epsfxsize=0.25\textwidth
\vskip 2 true cm
\hskip -2 truecm
\epsfbox[0 0 300 900]{disk_spots.eps}}
}
\parbox{\textwidth}{
\epsfxsize=0.35\textwidth
\center
\vskip -3 true cm
\epsfbox{rms_krakow.eps}}
\caption{Model of the variability in MCG-6-15-30 in X-ray band. Each magnetic flare
creates a hot spot at the disk surface (left), the energy-dependent single spot emission
is integrated assuming certain inclination angle of an observer, a time-dependent spot
distribution across the disk surface is generated (right), and the time-dependent
spectrum is calculated taking into account the relativistic effects. Finally, fractional
variability amplitude (upper curve) and point-to-point fractional variability amplitude
are obtained and compared to the observed one.
\label{fig:flary}}
\end{figure}

Chandra and XMM-Newton confirmed that in some objects the very broad emission
lines (iron line and soft X-ray lines) are seen, with the shape well
represented as the effect of the relativistic smearing due to the Keplerian
motion of the emitting material located very close to the black hole \cite{reviron}.
Some observations indicate that the inner radius of the disk increases when the source
becomes fainter \cite{matt}.

However, the observed variability of the iron line is not well understood and remains a
major puzzle \cite{longinotti,bianchi}. Models require some fine tuning in order to reproduce the
observed trends.

Models which explore full observational information have even more difficulties. We attempted
to reproduce both the fractional variability amplitude and point-to-point fractional variability
amplitude for MCG-6-15-30. The model assumed that magnetic flares are randomly distributed
above the disk surface. The flare flux, flare duration and the probability of a flare at a given
radius were assumed to have a power law dependence on the radius. Black hole was assumed to
rotate, disk was assumed to extend down to the marginally stable orbit appropriate for an adopted
Kerr parameter. Each flare created a hot spot on the disk surface by irradiation and the
reflection (including iron line and other soft lines) was calculated using the code
{\sc titan/noar} \cite{titan} and taking into account that the incident flux depends on the
distance from the flare center. Flare was assumed to be short-lasting so the vertical structure
of the disk was taken from an unilluminated model and the disk expansion was neglected. Predicted
spectrum was integrated in specific time bins, as in the
data, all relativistic effects were included using the code {\sc ky} \cite{ky}.

Within this scenario, we could not find a parameter range which would reproduce the shape of both
fractional variability functions. Either still the model is too simple (we did not consider
the dependence of the shape of the reprocessed radiation on the radius) or an important element
is still missing. There is some evidence that variations in the outflowing warm absorber can
contribute significantly. It was even suggested that warm absorber is mostly responsible for the
observed shape of the spectrum in soft X-ray band \cite{gierdone}, and the broad iron line is
partially an artefact \cite{pounds05}.

\section{Shocks and magnetic field reconnections as a heating mechanism of the hot plasma}

Although the geometry of the hot material is still under discussion, some mechanisms of
plasma heating are needed to explain the observed X-ray emission.

Models that can apply to magnetic flares above the disk surface were discussed in numerous
papers \cite{roma,liu}.

In ADAF type solutions most of the gravitational energy is used to heat ions, and subsequently
directed towards electrons through Coulomb interaction \cite{SLE,macha}. However, a fraction
of energy will be inevitably used to heat electrons directly \cite{gena1} thus limiting the low
radiative efficiency of the flow.

Several other aspects were also discussed, like disk-corona coupling \cite{kun}, disk evaporation
and the ion irradiation of the disk \cite{spruit} but the picture
is far from being complete.

Observationally, very interesting results were obtained in the
context of a galactic source (microquasar) GRS~1915+105 \cite{febel} and they may apply to
AGN when timescales are corrected by the black hole mass ratio. Methods, like Fourier-resolved
spectroscopy, successfully applied to galactic sources \cite{zycki2003}, may be also helpful although the
AGN data at present are hardly of the appropriate quality.

Therefore, observational determination of the gamma-ray emission from radio-quiet objects

\section*{Acknowledgments}
We thank Aneta Siemiginowska and Piotr \. Zycki for very helpful discussions.
Part of this work was supported by the grant PBZ-KBN-054/P03/2001
of the Polish
State Committee for Scientific Research, and by the Laboratoire Europe\' en Associ\' e Astrophysique
Pologne-France.

\section*{References}

Radio Quiet AGN