Quasar Absorption Lines and the Intergalactic Medium

Abstract: The importance of HST for the study of quasar absorption lines and of the nature of the intergalactic medium is illustrated by reviewing selected results from past HST observations. Topics reviewed include the study of Ly-alpha absorbers at low redshift and the search for a diffuse IGM at high redshifts.

\section{Opening the Ultra-violet Window}

Soon after quasars were recognized as extragalactic sources (Schmidt
1963; Greenstein and Matthews 1963; Schmidt 1965) it was pointed out
that as their light travels to Earth any intervening matter will leave
its imprint on the spectra.  This is true whether the intervening
medium is full of diffuse hydrogen causing a uniform decrease of the
quasar continuum (Gunn \& Peterson 1965; Scheuer 1965) or discrete
clouds producing separate absorption lines (by hydrogen, possibly
associated with galaxies, Bahcall \& Salpeter 1965; by hydrogen and
other species in galactic halos, Bahcall \& Salpeter 1966, Bahcall \&
Spitzer 1969).  However, even before the discovery of quasars, Lyman
Spitzer (1956) had pointed out the importance of UV spectroscopy for
understanding the physical conditions of the gaseous content of the
Galaxy, and by implication the halos of other galaxies and the gaseous
content of the universe in general.  As he pointed out, the majority
of the strong resonance absorption lines occur in the rest frame UV
(see Figure 1).
\begin{figure}[htb]
\centering
\plotfiddle{jannuzib1.eps}{3.1truein}{0}{80}{80}{-270}{-180}
\caption{The majority of the strong resonance absorption lines have
rest wavelengths which occur in the UV.  In this figure are plotted
the redshift ranges over which selected lines are observable in the IR
(applies only to Mg II), optical (solid line), ultra-violet
(observable with HST; dashed line), or far-UV (not observable with
HST, but observable with HUT, FUSE, or other UV missions; thin dashed
line).  Note that Galactic Lyman-limit absorption precludes any
instrument from being able to observe distant He~I or He~II absorption
below respective redshifts of 0.56 and 2.0.}
\end{figure}

At high redshifts intervening absorption systems can be studied from
the ground (for an example of the current state of the art see Figure
2). However, the gaseous content of the nearby universe and the far-UV
lines (e.g. He~II) occurring at high redshift can only be observed in
the UV.  Although HST is not the first telescope with UV sensitive
spectrographs, it is the first to provide both the spectral resolution
and the sensitivity to allow the extensive observation of the quasar
absorption lines.

Here is a list of three of the many important studies that the
observation of quasar absorption lines at UV wavelengths makes
possible:
\noindent
1.) The evolution of the gaseous content of the universe can now be
traced by observing the changing number density per unit redshift of
Ly$-\alpha$ absorbers from the present (using HST, from $z_{\rm abs}
\approx 0.0$ to $z_{\rm abs} \approx 1.6$) back to when the universe
was 10\% of its current age (using groundbased telescopes like Keck to
observe the most distant quasars).  Such data are important for the
study of cosmology, star and galaxy formation, development of large
scale structures, and the composition of the ISM.  The dramatic
changes that occur from high to low redshifts are illustrated by
comparing Figures 2-4.  Quantifying and understanding these changes in detail
is the continued focus of the HST quasar absorption line key
project and  other research efforts as well.

\begin{figure}[htb]
\centering
\plotfiddle{jannuzib2.eps}{1.25truein}{0}{93}{93}{-290}{-390}
\caption{The large number of absorbers blueward of Ly$-\alpha$, ``the
forest'', is quite evident in this Keck~I and HIRES observation of
the $z=3.63$ quasar Q1222$+$2309 obtained by Womble et~al. (1996).}
\end{figure}
\begin{figure}[htb]
\centering
\plotfiddle{jannuzib3.eps}{1.27truein}{0}{93}{93}{-290}{-390}
\caption{The number of absorbers blueward of Ly$-\alpha$ is noticeably
less in this HST and FOS spectrum of PG~1634$+$706 (Bahcall
et~al. 1996).  Each Ly$-\alpha$ line is indicated with a vertical
tick mark above the location of the line.}
\end{figure}
\begin{figure}[hbt]
\centering
\plotfiddle{jannuzib4.eps}{1.25truein}{0}{93}{93}{-290}{-390}
\caption{By the low redshifts probed by this HST observation of the quasar
3C~273 (Bahcall et~al. 1991) the Ly$-\alpha$ forest seems to have been
completely removed.  However, the number of absorbers is still in
excess of a simple extrapolation of the evolution observed at
redshifts greater than 1.6 (Bahcall et~al. 1991; Morris et~al. 1991)
Each Ly$-\alpha$ line is indicated with a vertical tick mark above the
location of the line. The other absorption lines are caused by the ISM
of our own Galaxy.  Note the two Ly$-\alpha$ lines near zero redshift
at the approximate velocity of the Virgo cluster.}
\end{figure}

\noindent
2.) For low redshift absorbers it is now possible to study their
relation to individual galaxies, groups, or clusters -- i.e. test the
1969 proposition of Bahcall and Spitzer, ``that most of the absorption
lines observed in quasi-stellar sources with multiple absorption
redshifts are caused by gas in extended halos of normal galaxies.''

\noindent
3.) Although efforts to detect a diffuse intergalactic medium (IGM)
through continuous absorption by neutral hydrogen have failed, we can
now extend the search by looking for He~II absorption that may be
easier to detect since the lower level of the ionizing background at
short wavelengths means that the He~II fraction is larger than the
fraction of H~I.

In the following sections I will give examples of how HST has
contributed to each of these problems.  The examples are drawn from
the work of the HST quasar absorption line key project and a few other
groups to illustrate this progress, with apologies to the many other
researchers that have also contributed to the study of these and other
problems through the use of quasar absorption lines observed with HST.

\section{The HST Quasar Absorption Line Key Project}

\subsection{Design and Goals of the Survey}

The HST quasar absorption line survey was an HST key project for
cycles 1-3, with carryover observations extending into cycle 4.  Led
by John Bahcall, the survey had the ambitious goal of obtaining a
large and homogeneous catalogue of absorbers suitable for the study of
the nature of gaseous systems and their evolution (Bahcall et~al.
1993).  While the well known telescope problems in effect prior to the
servicing mission reduced the original scope of the survey, the key
project still successfully observed 89 quasars with the higher
resolution (R$=1300$) gratings of the Faint Object Spectrograph. A
small subset of the quasars were observed from 1150--3300 \AA, but the
majority were observed only between 2200--3300 \AA~ or 1600--3300 \AA,
depending on the redshift of the quasar.  Targets were selected to be
bright and have low Galactic extinction ($b>20$ degrees).  The
distribution of the targets is shown in Figure 5. Redshifts of the
observed quasars range between 0.25 and 2.0.  Details of the data
calibration and analysis can be found in Schneider et~al.~1993,
Jannuzi \& Hartig 1994, and Bahcall et~al.~1996.  In all of our
analysis (from line measurements to line identifications) we have
tried to remove subjective decision making from the process and
replace it with well tested algorithms implemented through computer
software.  This allows us to run the same software on simulated data
in order to improve our understanding of the limitations of both our
data and our analysis techniques.

\begin{figure}[htb]
\plotfiddle{jannuzib5.eps}{2.35truein}{0}{70}{70}{-220}{-190}
\caption{The distribution in galactic coordinates of the 89 quasars
observed as part of the HST quasar absorption line survey.}
\end{figure}

\subsection{Past Results}

The first nice surprise that HST presented to us was a larger number
of low redshift Ly$-\alpha$ absorption systems in the spectrum of
3C~273 than might have been expected from a simple extrapolation of
the evolution in the number density (per unit redshift) of such
systems observed at high redshift (Bahcall et~al.~1991; Morris et
al.~1991).  Early results have also been produced by HST on the nature
and evolution of metal line systems (e.g. Reimers \& Vogel 1993;
Bergeron et~al.~1994).  Our understanding of the evolution of
Lyman-limit systems from high redshifts down to $z=0.4$ has been
improved (Storrie-Lombardi et~al.~1994; Stengler-Larrea et~al.~1995)
and a first attempt has been made to measure the proximity effect in
the spectra of low redshift quasars (Kulkarni \& Fall 1993).

The key project catalogue of Ly$-\alpha$ absorbers makes it easier to
investigate the extent and nature of the relationship between
Ly$-\alpha$ absorbers and individual galaxies, groups, or clusters.
Many groups are actively working on this problem (e.g. incomplete or
single field surveys: Bahcall et~al.~1991, 1992; Morris et~al.~1993;
Spinrad et~al.~1993; to more extensive surveys in progress that have
presented partial results: Lanzetta et~al.~1995; Stocke et~al.~1995;
Le Brun et~al.~1995), but I have chosen to adapt a figure from Morris
et~al.'s (1993) study of the field of 3C~273 to illustrate both the
progress that has been made and how much more needs to be done
(Figure~6).  Despite a complete redshift survey of even the faintest
galaxies in the field, the study of the 3C~273 field gives a mixed
signal.  While some absorbers appear to be associated with the same
structures as the galaxies (as suggested by Lanzetta et~al., actually
part of the halos of the galaxies) other lines appear in voids (see
also Stocke et~al.) with no detected galaxy within 1 Mpc.  The Morris
et~al.~study is limited by the small number of Ly$-\alpha$ systems
along the line of sight toward 3C~273, resulting in a limited
comparison between the distribution of galaxies and absorbers.  In
fact no single line of sight provides enough Ly$-\alpha$ absorbers to
allow the accurate determination of the fraction of all absorbers
which are associated with galaxies or larger structures.  For some of
the other papers listed above the problem is similar or the galaxy
redshift surveys that they use are incomplete. Some of the other
surveys are also not able to address the relationship between the
absorbers and groups or clusters of galaxies because the galaxy survey
does not cover a large enough angular area to be able to identify a
cluster or group.  To determine accurately the fraction of Ly$-\alpha$
absorbers associated with galaxies and large scale structures requires
both the completion of the key project catalogue of absorbers and an
increase in the number of fields for which galaxy redshifts are
available (e.g. Sarajedini et~al.~1996).

\begin{figure}[htb]
\centering
\plotfiddle{jannuzib6.eps}{4.0truein}{0}{88}{88}{-258}{-40}
\caption{Displayed in these pie-diagrams are the positions of the
galaxies in the field of the quasar 3C~273 and the locations of
Ly$-\alpha$ lines detected in the FOS and GHRS ultra-violet spectra of
3C~273.  Angles have been exaggerated by a factor of 15 to improve the
clarity of the figure, but results in a distorted plot with spherical
structures appearing elongated transverse to the line of sight. Note
that while some of the Ly$-\alpha$ absorbers appear associated with
galaxies, several have no detected galaxy within 1 Mpc.  Adapted from
Morris et~al.~1993, see their paper for complete discussion.}
\end{figure}

The same key project spectra have also provided valuable information
on our own Galaxy's halo and ISM (Savage et~al.~1993), the emission
line properties and spectral energy distributions of quasars (Espey
et~al.~1994; Weymann et~al.~1996; Laor et~al.~1994, 1995; Sulentic
this conference), and warm x-ray absorbers (in the quasar 3C~351,
Mathur et~al.~1994).

\subsection{Some New and Future Results}

A continuing focus of the key project is to determine the nature and
evolution of the low redshift Ly$-\alpha$ absorbers.  The number
density of such systems as a function of redshift is summarized in
Figure~7 (see Bahcall et~al.~1996 for details).  At low redshift ($z<
1.3$), the key project data analyzed to date (about 10\% of the
expected final catalogue) is consistent with no evolution for $\gamma
= 0.58 \pm 0.50$ and $ dN/dz \propto (1+z)^{\gamma} $. This result is
derived from a maximum likelihood estimation for the observed lines in
those spectral regions where the 4.5 $\sigma $ detection limit is less
than 0.24~\AA.  We further find that the slope of the
observed low-redshift $dN/dz$ relation differs at the $2-4.5\sigma$
level of significance from the slope deduced from various ground-based
samples that refer to redshifts $z > 1.6$ (Lu et~al.~1991; Press et
al. 1993; Bechtold 1994).

\begin{figure}
\plotfiddle{jannuzib7.eps}{4.9truein}{0}{80}{80}{-220}{-98}
\caption{The number density of Ly$-\alpha$ absorption systems is shown as a
function of redshift.  At redshifts below 1.6 the data come from
Bahcall et~al.~1996, including approximately 10\% of the final HST key
project catalogue of Ly$-\alpha$ absorbers.  At higher redshifts two
samples (Lu et~al.~1991; Bechtold 1994) are plotted which give
similar but not identical results. The separate fits to the HST and Lu
et~al.~data are shown as the two solid lines, with slopes of
$\gamma=0.48\pm0.54$ and $2.68\pm0.27$ respectively. The dashed line
shows the best fit to a single power law for both the HST and the Lu
et~al.~data and has a slope of $1.58\pm0.13$. A KS test indicates that
this fit is only acceptable at the 2.7\% level. Fits using the
Bechtold (1994) data give similar results although a single power law
fit is not ruled out as strongly. The fits were done using a
maximum likelihood technique on the unbinned data. The binned data
have been placed on the figure for reference purposes only.}
\end{figure}

As the number of absorbers in the analyzed catalogue increases it
becomes possible to study the clustering properties of Ly$-\alpha$
absorbers.  While we have yet to detect any signal in the two-point
correlation function, we have found evidence that about half of the
extensive metal line systems seen at redshifts between 0.4 and 1.3 are
accompanied by highly-clustered clumps of Ly$-\alpha$ lines which are
physically associated with the metal-line systems (details in Bahcall
et~al.~1996).

Our understanding of both the redshift evolution of all absorption
systems and of their clustering properties will improve as we complete
the catalogue of absorption systems.  The last observation of the key
project was made in May of 1995.  At the time of this meeting all of
the quasar spectra have been reduced, lines measured, and the lines
are being identified. The key project results I have reviewed have
been based on only part of the total absorption line data set (see
Figure 8). While we have analyzed one sixth of the objects, the
remaining five sixths include most of the higher redshift objects and
four fifths of the observed redshift path length.  Expected
improvements upon completion of the catalogue include: 1.)  examining
the evolution of Ly$-\alpha$ systems not only as a function of
redshift, but also as a function of neutral column density and 2.)
confirming or refuting the preliminary evidence for clustering of
Ly$-\alpha$ absorbers around metal line systems.

\begin{figure}
\plottwo{jannuzib8.eps}{jannuzib9.eps}
\caption{Shown here are histograms depicting the number of quasars
observed as part of the absorption line key project which will
contribute information at each redshift.  The solid region indicates
the number of quasars whose data have been included in the published
papers and the outlined region indicates the complete sample currently
being analyzed. When the analysis is complete the size of the
absorption line catalogue will increase by a factor of 10, greatly
improving the determination of the evolution depicted in Figure 7.}
\end{figure}

\section{Has the IGM Finally Been Detected?}

We now leave the universe at low redshift behind and examine the
Herculean efforts that have been made to detect the diffuse
intergalactic medium with HST (a complete and detailed account of this
exciting area can be found in the contribution to these proceedings by
Dr. Jakobsen).  Excluding the detection of absorption assigned to weak
individual ``clouds'' (the low column density end of the Ly$-\alpha$
forest; the Bahcall-Salpeter effect), all efforts to detect absorption
by diffuse neutral hydrogen (the Gunn-Peterson effect) at ANY redshift
have failed.  The ionizing background radiation reduces the fraction
of H~I and He~I (Sargent et~al.~1980) and removes them as probes of
the diffuse IGM (note, He~I has been observed in high column density
systems, the first detection being made with HST by Reimers \& Vogel
1993).  The lower ionizing background at short wavelengths might leave
a higher fraction of He~II and provide a means of detecting the
diffuse IGM, but the short wavelength of He~II (304 \AA) means that it
can only be observed at high redshifts.  This means that ``clear''
quasars, as Jakobsen call them, must be found.  Such quasars must be
bright, have redshifts greater than 3 (to have He~II observable with
HST), and be free of significant absorption from intervening
``clouds'', particularly the high column density systems whose
Lyman-Limit absorption would preclude the observation of He~II.
Jakobsen et al. (1994) and Tytler et~al.~(personal communication) both
conducted searches during cycles 1-3 using HST and respectively the
FOC and FOS to find candidate ``clear'' quasars.  In cycle 4 both
groups succeeded in detecting He~II absorption in the spectra of
distant quasars.

While the details can be found in Dr. Jakobsen's contribution, here
are some bottom lines.  There are now three detections of absorption
due to He~II. Two quasars (Q0302$-$003 and PKS1935$-$692) exhibit
black, continuous absorption blueward of the expected wavelength of He~II
at the redshift of the observed quasars ($z > 3$ for the two observed
with HST, see figures in Jakobsen et~al.~1994 and Jakobsen's
contribution; Jakobsen and Tytler 1996). The lower limits on the
optical depth of He~II absorption are 1.7 in both cases. A third
quasar, HS1700$+$6414, was successfully observed by HUT (Davidsen
1995) and shows He~II absorption beginning at $z=2.7$, but the
absorption is not as strong as at the higher redshifts observed with
HST.  It appears that both the HST and HUT observations can be
interpreted as consistent with each other given the possible evolution
between redshifts of 3.2 to 2.7 (Jakobsen, this meeting).

One problem with any search for a diffuse component in the IGM is that
as we are able to detect and resolve lower column density systems we
remove absorption from the previously unresolved ``diffuse'' component
and move it into the ``cloud'' component.  Q0302$-$003 has been
observed with Keck and the HIRES echelle spectrograph and a population
of very low column density Ly$-\alpha$ clouds has been detected by
Songaila et~al.~(1995) and they report that the detected population is
extensive enough that it is possible to explain the observed He~II
absorption without invoking a diffuse IGM. Neither the HUT nor the HST
FOC and FOS observations have the spectral resolution necessary to
distinguish directly between the He~II ``forest'' and a more diffuse and
uniform absorption.  The issue is likely to remain unsettled until the
existing lines of sight (or additional new detections, hopefully with
brighter background quasars) are successfully observed at a high
enough spectral resolution that the He~II associated with the hydrogen
forest clouds can be resolved.

There is a second complication.  Cosmological simulations of the
universe at intermediate and high redshifts (e.g. Katz et~al.~1996)
indicate that we should now expect a complex distribution for the gas
in the IGM with filamentary structures covering a large range of
physical scales and conditions. There might not exist any component
that matches our expectations of a smooth or uniform ``diffuse IGM''.
It might be that the distinction between numerous, closely packed,
very diffuse (low column density) ``clouds'' and a more uniform
diffuse medium is purely a question of semantics, but the resolution
of this issue has implications for a variety of issues, including
understanding the physical conditions that exist during the formation
of galaxies (see Jakobsen, these proceedings for further discussion).

\section{End Matters}

The 1990's is the epoch of two revolutions in the the study of quasar
absorption lines.  Prior to HST and the Keck telescopes, quasar
absorption lines have been discussed and studied in distinct
subgroups, roughly separated by column density. At the extremes were
the Ly$-\alpha$ forest lines that were observed to be unclustered and
possibly composed of primordial material (based in part on the lack of
any detected metal line absorption) and the damped systems with their
high column densities and large gas masses identified as the
progenitors of spiral galaxies (e.g. Wolfe 1988).  Such divisions,
while still useful, are getting fuzzy as new results rapidly blur
distinctions.  Just one example (of many) is the detection of weak CIV
absorption associated with some fraction of low column density
Ly$-\alpha$ absorbers, systems that would have previously been
securely identified as part of the primordial ``forest clouds''
(e.g. Cowie et~al.~1995; Womble et~al.~1996). Such wonderful
observations require modification of the pre-HST-Keck picture of
absorption line systems.  How should we modify the old ``standard
picture''? I am not sure.  But I do think that a second revolution is
going to provide critical guidance in the development of the new more
complex and detailed models.  The second revolution is the progress
theorists have made in leaving behind spherical cloud, slab, or
mini-halo models and replacing them with the help of super computers
to generate full hydrodynamic and SPH simulations of the evolution of
the universe. Three groups are now able to not only generate
simulations of large scale structures, but also simulated quasar
absorption line spectra along numerous lines of sight through their
simulations that can be compared to real observations (see Zhang,
Y. et~al.~1995; Hernquist et~al.~1996 and Katz et~al.~1996; Cen
et~al. 1994, Miralda-Escud\'e et~al.~1996). The challenge ahead is to
extract the best set of observables from both the simulations and the
various data sets so that cosmological models might be discriminated
against. Furthermore, enough simulations (and observational data!)
need to be generated that the uniqueness of ``good fit'' models can be
tested.

In his introduction to the Hubble Deep Field project, Bob Williams
ably described how HST has opened up the distant universe to our view.
He speculated that one of the Hubble Space Telescope's lasting and
important legacies would be providing us our first ``clear'' images of
the early history of the universe.   In the future HST will also be
remembered for making possible unique studies of the more evolved and
nearby universe.  WFPC-2 is providing exquisite images of galactic
sources and nearby objects that reveal a wealth of previously
unobservable detail (see for examples the contributions to these
proceedings by Bally, Livio, Machetto, and O'Dell). But HST should
also be remembered for the unique information provided by its
spectrographs.  By making it possible to study quasar absorption lines
in the ultra-violet HST has already provided important data
about the gaseous content of the universe at both low and high
redshifts. This legacy will continue to grow as existing data is
further analyzed and when STIS makes its appearance on HST.

\acknowledgments
I thank Jill Bechtold, Simon Morris, Donna Womble and Wal Sargent, and
the HST quasar absorption line key project team for providing data
used in the figures and acknowledge valuable discussions with Jill
Bechtold, Peter Jakobsen, David Weinberg, and the entire HST quasar
absorption line survey key project team. Hans-Walter Rix and David
Weinberg provided useful comments on an early version of this paper. I
thank the meeting organizers, particularly the local organizers, for
managing to host such an enjoyable conference given the logistical
problems brought about by the Paris transit strike.


Observational View of Magnetic Fields in Active Galactic Nuclei Jets

Abstract: According to the currently favored picture, relativistic jets in active galactic nuclei (AGN) are launched in the vicinity of the black hole by magnetic fields extracting energy from the spinning black hole or the accretion disk. In the past decades, various models from shocks to magnetic reconnection have been proposed as the energy dissipation mechanism in the jets. This paper presents a short review on how linear polarization observations can be used to constrain the magnetic field structure in the jets of AGN, and how the observations can be used to constrain the various emission models.

\section{Introduction}
Active galactic nuclei (AGN) are among the most extreme objects in the
universe. They are the centers of distant galaxies hosting a
supermassive black hole, billions of times the mass of the Sun. About
10\% of AGN produce collimated, relativistic outflows or plasma jets
that shine brightly over the entire electromagnetic spectrum from
radio to gamma-ray energies. The formation and stability of these jets
are not yet fully understood. It is believed that energy is extracted
from the central black hole by large-scale magnetic fields
(\cite{blandford77}) while rotating inner regions of the accretion disk
form an outflowing wind that helps to collimate the jet
(\cite{blandford82}). In this picture, the jet is launched as a
magnetically dominated outflow with strong magnetic fields
accelerating the flow to relativistic velocities. General relativistic
magnetohydrodynamic simulations show that the efficiency of
the jet generation can be very high if the black hole is threaded by a
dynamically important magnetic field (\cite{tchekhovskoy11}).

While the simulations suggest that the jets are likely highly
magnetized near the black hole, their composition and magnetic field
structure in the parsec scales, where they produce most of the emission, is much less
clear. A link between the magnetic fields near the black hole and the
parsec-scale jets was established by \cite[Zamaninasab et al. (2014)]{zamaninasab14} who showed
that there is a correlation between the magnetic flux of the jet and
the accretion disk luminosity, as predicted by the magnetically
arrested disk (MAD) model (\cite{tchekhovskoy11}). This result
illustrates how studying the emission regions in the parsec-scale jets
may help us gain knowledge about the jet formation processes.

One of the main open questions is the nature of the emission region.
Traditionally, the flaring behavior of AGN jets has been well
explained by shock-in-jet models, especially in the radio bands
(e.g., \cite{hughes85, marscher85}). These models assume the jet to
be kinetically dominated at parsec scales. Shock models struggle to
explain the very fast variability detected at high-energy gamma-ray
bands, and models involving magnetic reconnection have been
invoked to explain the high-energy emission
(e.g., \cite{giannios09}). In this case the jet flow should be
magnetically dominated, although recent particle-in-cell simulations
show that even in reconnection models, the magnetic field and particle energies can be in
equipartition at the emission region (\cite{sironi15}).

Another major open question is the structure of the magnetic field at
the emission site. According to the simulations, the magnetic field is
helical close to the black hole but whether it preserves its order in
the parsec scales is debated upon. For example, in the model of
\cite[Marscher et al. (2008)]{marscher08}, the magnetic field is helical at the acceleration
and collimation zone of the jet but gets disrupted and turbulent
beyond a standing shock in the jet. Polarization observations of some
jets, however, indicate that the magnetic field could be helical even
on the parsec scales (e.g., \cite{asada02, gabuzda04,
  hovatta12}). Figure \ref{fig1} shows an example adapted from
\cite[Hovatta et al. (2012)]{hovatta12}.

\begin{figure}[b]
\begin{center}
 \includegraphics[width=\textwidth]{3C273_gradient_new.eps}
 \caption{Faraday rotation map of 3C 273 obtained with the VLBA at 8 -
 15 GHz. A rotation measure gradient is seen extending the full length
of the jet. This can be interpreted as a signature of a helical
magnetic field. Figure adapted from \cite[Hovatta et al. (2012)]{hovatta12}.}
   \label{fig1}
\end{center}
\end{figure}

In this contribution, I will discuss how linear polarization observations in the
radio and optical bands can be used to constrain jet emission
models and the magnetic field structure of the jets. This is not
intended as a comprehensive review of all possible models, but more to
give some examples on how observations are used to study magnetic
fields in AGN jets.

\section{Polarization observations}\label{obs}
The optical and radio emission of AGN jets is synchrotron emission,
which is intrinsically highly polarized. In an optically thin emission
region with a uniform magnetic field,
the polarization degree is up to 70\% (e.g., \cite{pacholczyk70}),
although such high polarization degree values are not typically seen
in AGN jets in the radio (e.g., \cite{aller03, lister05}) or optical
(e.g., \cite{angel80, pavlidou14}) bands. This has been taken as
evidence for disordered magnetic fields. The emission in AGN jets is
often described with the Stokes parameters \textit{I} (for total
intensity), \textit{Q} and \textit{U} (for linear polarization) and
\textit{V} (for circular polarization). I will only discuss linear
polarization here. Using the Stokes parameters, the polarization
degree and the electric vector position angle (EVPA) can be defined as
$m=(\sqrt{Q^2+U^2})/I$ and EVPA$=1/2\tan^{-1}(U/Q)$. In the simplest
  case, in an optically thin jet, the polarization position angle is
perpendicular to the magnetic field direction and one can use the EVPA observations
to infer the direction of the magnetic field.
However, one should note that in relativistic jets viewed at a small angle to the line of
sight of the observer, the situation is more complicated (\cite{lyutikov05}).

Another way to obtain information about magnetic fields in blazar jets
is via Faraday rotation observations. When synchrotron radiation
passes through magnetized plasma, Faraday rotation of the EVPAs
proportional to the line-of-sight magnetic field and electron density
occurs (e.g. \cite{burn66}). In the simplest case, the effect can
be described by a linear dependence between the observed electric
vector position angle (EVPA; $\chi_\mathrm{obs}$) and wavelength
squared ($\lambda^2$), given by
\begin{equation}\label{eq:RM}
\chi_\mathrm{obs} = \chi_0 + 0.81 \int n_e \mathbf{B} \cdot \mathbf{\mathrm{d}l} = \chi_0 + \mathrm{RM}\lambda^2,
\end{equation}
where $\chi_0$ is the intrinsic EVPA and RM is the rotation measure
(in rad/m$^2$), related to the electron density $n_e$ of the plasma (in cm$^{-3}$) and the magnetic
field component $\mathbf{B}$ (in $\mu$G) along the line of sight (in parsecs). The RM
can thus be estimated by observing the EVPA at different
frequencies. This will give us information on the line-of-sight
component of the magnetic field. For example, if the rotation measure is positive,
the magnetic field is coming towards the observer, and if negative it
is going away from the observer.
Therefore it was suggested that a gradient in a Faraday rotation
measure transverse to the jet direction could reveal helical magnetic field structures (\cite{blandford93}).

\section{Constraining emission models through polarimetry}
Polarization observations of AGN jets in the radio and optical bands
have been conducted since the 1970s (see e.g., \cite{saikia88}
for a review). At the same time, a large number of theoretical models
were developed to explain the total intensity
and polarization observations of AGN jets (e.g, \cite{jones77,
blandford79, laing80}). In the 1980s it was suggested
that most of the radio variability, especially in the cm-band, are due
to shocks in the relativistic jets where the magnetic field is
predominantly turbulent (\cite{hughes85, marscher85, jones85}).

The field of polarization modeling has received a new boost since an
exciting connection to high-energy gamma-ray emission was found
(\cite{marscher08}). A rotation in the optical polarization angle
was coincident with an ejection of a new very
long baseline interferometry (VLBI) component from the core, and a
very-high-energy (VHE; $>$ 100 GeV) detection of the source BL Lac by
the MAGIC telescope. A further connection between the optical
polarization and gamma-ray flares was suggested in two quasars 3C~279
(\cite{abdo10}) and PKS~1510-089 (\cite{marscher10, aleksic14}). In both of
these the optical polarization angle was seen to rotate by more than
180 degrees over a course of 20 to 50 days during which a sharp
gamma-ray flare was observed. Following all these observations, a
number of new models have been published in the last years
(e.g., \cite{marscher14, zhang14, hughes15, zhang15}). I will go through some of these models
with the emphasis on the observations that can be used to constrain
them.

\subsection{Turbulence in the jets}
Turbulent magnetic field produces stochastic variations to the
polarization angle and degree (\cite{jones85, marscher14}), such
as an EVPA rotation of any length, typically accompanied by a low
polarization degree. The challenge in constraining these types of
models is the stochastic nature of the variations so that it is not
possible to fit the individual light curves
(\cite{marscher14}). Instead, one should observe a large number of
objects and compare the statistical properties of the variations with
the turbulent models.

To study the connection between high-energy emission and optical
polarization in a statistical manner, the RoboPol program was initiated
in 2013 (\cite{pavlidou14}). The RoboPol instrument is mounted on the
1.3-m telescope in Skinakas Observatory in Crete, where the observing
season lasts from April until November. During its first three observing
seasons, RoboPol observed about 100 AGN twice per each week in search
of optical polarization angle rotations. The main goal was to study
the differences in gamma-ray detected and non-detected objects, and to
search for a connection between rotations and flaring behavior.

RoboPol has detected 40 optical polarization angle rotations during
the first three observing seasons, tripling the number of known rotations (\cite{blinov15, blinov16a,
blinov16b}). Comparing the rotations with random walk simulations showed
that while stochastic variability for any individual rotation could not
be ruled out, it was highly unlikely that all of them would be of
random walk origin. \cite[Kiehlmann et al. (2016)]{kiehlmann16} suggested that the smoothness
of the rotation could be used as an additional indicator when
comparing the rotations with random walk models. They also illustrated the
importance of good sampling in polarization observations by revealing
how additional data changed the rotation reported by
\cite[Abdo et al. (2010)]{abdo10} in the quasar 3C~279.

Figure \ref{fig2} shows an example of a EVPA rotation in the source
1ES~1727+502 studied in \cite[Hovatta et al. (2016)]{hovatta16}. Based
on random walk simulations, conducted in the same manner as in
\cite[Kiehlmann et al. (2016)]{kiehlmann16}, a stochastic origin for
the rotation cannot be excluded.

\begin{figure}[b]
\begin{center}
 \includegraphics[width=3.4in, angle=-90]{RBPLJ1728_rotation.eps}
 \caption{Polarization degree and EVPA of 1ES~1727+502 in 2014
   observed with RoboPol and the Nordic Optical Telescope (NOT). The
   EVPA shows a rotation of nearly 150 degrees in about 60 days. Figure adapted from \cite[Hovatta et al. (2016)]{hovatta16}.}
   \label{fig2}
\end{center}
\end{figure}

\subsection{Emission in a helical field}
A rotation in the polarization angle could also be due to an emission
feature tracing a magnetic field in the jet as suggested by
\cite[Marscher et al. (2008)]{marscher08} and \cite[Marscher et al. (2010)]{marscher10}. In this model, the magnetic
field in the jet is helical in the acceleration and collimation zone,
which is probed by the optical band, and a rotation is seen when the
emission feature moves along the magnetic field. The rotation may be
accompanied by flaring in other bands when the emission feature
reaches a standing shock in the jet.

Another alternative is a shock moving down a jet with a helical field,
in which case the EVPA rotation would be due to light travel time effects
when parts of the shock are seen at different times
(\cite{zhang14}). This model was used to successfully fit the EVPA rotation in the
quasar 3C~279 (\cite{zhang15}), although one should bear in mind the
caveat that with the additional data from \cite[Kiehlmann et al. (2016)]{kiehlmann16}, the
rotation was not as long as originally presented by
\cite[Abdo et al. (2010)]{abdo10}.

Whether the magnetic field in the jet is helical, especially in
regions beyond the acceleration and collimation zone, is still
unclear. As explained in Sect.~\ref{obs}, observations of a Faraday
rotation measure gradient across the jet could be an indication of
such a field. First signatures of a helical magnetic field were
observed in the quasar 3C 273 using multifrequency very long baseline
array (VLBA) observations (\cite{asada02, zavala05}).
Many more claims of such gradients have been made (e.g., \cite{gabuzda04}) but the issue has been controversial due to the
limited resolution across the jets in blazars (\cite{taylor10}).

In \cite[Hovatta et al. (2012)]{hovatta12}, we performed Monte Carlo simulations to quantify
the significance of the rotation measure gradients in the large sample
of parsec-scale jets in the MOJAVE (Monitoring of Jets in Active
galactic nuclei with VLBA Experiments) sample. Our observations
confirmed the gradient in 3C~273 (see Fig.~\ref{fig1}), and we also found significant
gradients in three other quasars. In \cite[Zamaninasab et al. (2013)]{zamaninasab13} we modelled
the gradient of the quasar 3C~454.3 using a magnetic field with both
helical and turbulent components. In the other two sources, the
gradient span only over a small portion of the jet and detailed
modeling was not possible.

Nowadays performing simulations to quantify the significance of the
gradients is a common practice (e.g., \cite{murphy13}) and the number
of significant gradients has steadily increased (e.g.,
\cite{gabuzda15}). However, rotation measure gradients can also arise from
changes in the density of the Faraday rotating material
(e.g. \cite{gomez11}) so that
detailed modeling of the gradients should be done in order to confirm
that they are indeed due to helical fields. It is also unclear whether
the helical field is within the jet or in a sheath layer around
the jet because Faraday rotation is a propagation effect originating
in the plasma outside the emission region.

\subsection{Shock in a jet}
As stated earlier, shock-in-jet models have been successful in
reproducing the observed flaring especially in the radio
bands. Typically, radiative transfer simulations are generated to
estimate the parameters of the shocks (e.g. \cite{hughes11, hughes15})
and these are then compared with observations at multiple bands
(\cite{aller14}). For example, as shown by \cite{hughes15},
the polarization degree and range of EVPA values can
be used to constrain the magnetic field geometry and jet orientation,
while total intensity behavior is indistinguishable in different models.

Interestingly, although the magnetic field is assumed to be predominantly
turbulent, in some cases an ordered field component (possibly helical)
could also be present (\cite{aller16}). This supports the
findings of e.g., \cite[Zamaninasab et al. (2013), Blinov et al. (2015)]{zamaninasab13, blinov15}, and
\cite[Kiehlmann et al. (2016)]{kiehlmann16}, who also suggest that there are both turbulent
and ordered magnetic field components / deterministic variations in the jets.

\subsection{Magnetic reconnection}
Magnetic reconnection (see \cite{kagan15} for a review) is a new hot topic in the field.
It was originally invoked to explain the fast high-energy variability
of blazars through the jet-in-jet model (\cite{giannios09}). In the
recent years, particle-in-cell simulations have evolved significantly,
and it is now possible to reliably simulate the complicated structure
of the reconnection layer (e.g., \cite{sironi15}). While the
reconnection models can already be compared to total intensity
variability of the sources (\cite{petropoulou16}), there are no
explicit models for the polarized emission.

\cite[Zhang et al. (2015)]{zhang15} stated that the EVPA swing in the 2009 flare of 3C~279
was reproduced by a model that ``favors magnetic energy dissipation
process during the flare'' because the magnetic field strength was
seen to decrease when the total intensity increased. However, this
model did not yet include detailed comparisons to a specific
reconnection model. Considering the growing interest in the magnetic
dissipation models, it is likely that in the next few years more
detailed models with observational predictions will become available.

\subsection{Statistical trends}
From the previous sections, it is clear that polarization observations
can be used to constrain various types of emission models and the
magnetic field structure in the jets. In order to generalize the
findings in individual sources into the AGN jet population, a statistical approach must be
used. This is what the RoboPol program aims for by observing a sample
of about 100 sources with high cadence. RoboPol has found, for
example, that any class of blazars from low to high synchrotron
peaking objects can show rotations (see
also \cite{jermak16, hovatta16}). However, there seems to be a
specific class of ``rotators'' that do so more often, and rotations
are more common in the low spectral peaking objects (\cite{blinov16b}).

The low synchrotron peaking objects also have higher optical polarization degree
than the high-peaking objects (\cite{angelakis16}), a trend that has
also been seen in the smaller sample studied in \cite[Jermak et al. (2016)]{jermak16} and
earlier in the radio band (\cite{lister11}). These general trends can
be explained with a simple, qualitative model where a shock moves down
a jet, which has both helical and turbulent magnetic field components (see
\cite{angelakis16} for details). Any model put forward to explain the
polarization behavior in individual flares, should also account for
these general trends.

\section{Future directions}
One challenge in connecting the observations of magnetic fields to the
theory of jet formation is that often the observations, especially in
the radio bands,  probe regions $10^3-10^5$ gravitational radii away from the black hole.
Optical observations do not suffer from this restriction, but for
example, Faraday rotation with $\lambda^2$ wavelength dependence, is
not typically seen at optical wavelengths. A solution can be provided by going
to millimeter-band observations, as demonstrated by
\cite[Mart\'i-Vidal et al. (2015)]{marti-vidal15}. They used ALMA polarization observations to
probe the Faraday rotation at the jet base of the lensed quasar PKS~1830-211, and found extremely high Faraday rotation of
$10^8$rad/m$^2$, which is two orders of magnitude higher than previous
observations in other sources (e.g., \cite{plambeck14}). They inferred
this as a signature of a very high
magnetic field at the base of the jet, in support of magnetically
launched jet models. Whether similar high RM values are a common
property of quasars remains to be seen with future ALMA
observations. Especially interesting will be future observations using
ALMA as part of the global VLBI array, which may be able to spatially resolve
the regions of high magnetic fields.

A major leap forward will come with the future X-ray polarization
missions because the amount of polarization for different high-energy emission
mechanisms in AGN jets, for example between leptonic and hadronic
models, is very different (\cite{zhang13}), and X-ray polarization
observations can be used to distinguish between the models.

\section{Summary}
Magnetic fields in AGN jets can be probed through polarization
observations, which are most easily done at radio and optical
bands. Modeling the polarization degree and position angle behavior
can be used to test different blazar emission models, while Faraday
rotation observations at radio, and now also in the millimeter bands,
can be used to probe the line-of-sight magnetic field component. A
statistical approach is favored when connecting the results to
different AGN jet populations.

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


WISE/2MASS-SDSS brown dwarfs candidates using Virtual Observatory tools

Abstract: Massive far-red and infrared imaging surveys in different bandpasses are the main contributors to the discovery of brown dwarfs. The Virtual Observatory represents an adequate framework to handle these vast datasets efficiently and filter them out according to specific requirements. A statistically significant number of BDs is mandatory for understanding their general properties better for identifing peculiar objects. WISE, an all-sky survey in the mid-infrared, provides an excellent opportunity to increase the number of BDs significantly, in particular those at the lower end of the temperature scale. We aim to demonstrate that VO tools are efficient in identifing and characterizing BDs by cross-correlating public catalogues released by large surveys. Using VO tools we performed a cross-match of the WISE Preliminary Release, the 2MASS Point Source and the SDSS Data Release 7 catalogues over the whole area of sky that they have in common ($\sim$ 4000 deg$^{2}$). Photometric and proper motion criteria were used to obtain a list of BD candidates. A temperature estimate is provided for each candidate based on their spectral energy distribution using VOSA, a VO tool for SED fitting. We derive the spectral types from the effective temperatures. Distances, calculated from the absolute magnitude- spectral type relation, place our candidates at 14-80 pc from the Sun, assuming that they are single. We have identified 31 BD candidates, 25 of which have already been reported in the literature. The remaining six candidates have been classified as L- (four) and T-type (two) objects. The high rate of recovery of known BDs ($\sim$ 90% of the T dwarfs catalogued in 2MASS) demonstrates the validity of our strategy to identify them with VO tools.

\section{Introduction}

Brown dwarfs (BDs), self-gravitating objects that form like stars but do not get enough mass to maintain a sufficiently high temperature and pressure in their cores for stable hydrogen
fusion, provide a natural link between very low-mass stars and gaseous giant planets. Theoretically proposed fifty years ago (\cite{1963ApJ...137.1121K}), it was not until
1995 when the first BDs were discovered (\cite{1995Natur.377..129R}; \cite{1995Natur.378..463N}). Since then, hundreds of BDs have been found mainly thanks to the advent of large-area optical and
near-infrared surveys such as the Two Micron All Sky Survey (2MASS, \cite{2006AJ....131.1163S}; \cite{2007AJ....134.1162L}), the DEep Near Infrared Survey
of the Southern Sky (DENIS, \cite{1997Msngr..87...27E}; \cite{1997A&A...327L..25D}; \cite{1999AJ....118.2466M}), the Sloan Digital Sky Survey
(SDSS, \cite{2000AJ....120.1579Y}; \cite{2006AJ....131.2722C}), the UKIRT Infrared Deep Survey (UKIDSS,  \cite{2007MNRAS.379.1599L};
\cite{2007MNRAS.379.1423L}; \cite{2010MNRAS.406.1885B}) and the Canada-France Brown Dwarf Survey (CFBDS, \cite{2008A&A...484..469D}; \cite{2011AJ....141..203A})\footnote{See http://dwarfarchives.org for an updated list of L and T dwarfs}.

The Wide-field Infrared Survey Explorer (WISE; \cite{2010AJ....140.1868W}) is a NASA mission that has mapped the sky at 3.4 (W1), 4.6 (W2), 12 (W3), and 22 (W4)
$\mu$m in 2010 with an angular resolution of 6.1", 6.4", 6.5", and 12.0", respectively. WISE achieved 5$\sigma$ point-source sensitivities better than
0.08, 0.11, 1 and 6 mJy in unconfused regions on the ecliptic in the four bands. The dataset obtained by the WISE imaging survey constitutes an excellent resource for finding new brown dwarfs, in particular the coldest members, because they emit a substantial part of the spectral flux in the mid-infrared. The recent discoveries of T dwarfs reported by different authors
(Burgasser et al. 2011; Mainzer et al. 2011; \cite{2011A&A...532L...5S}) confirm the high expectations placed on the mission. The WISE Preliminary Release\footnote{http://irsa.ipac.caltech.edu/Missions/wise.html} is available to the astronomical community since April 14, 2011 and includes photometric information for over 257 million objects observed during the first 105 days of the survey.

The Virtual Observatory\footnote{http://www.ivoa.net} (VO) is an international initiative designed to help the astronomical community in the exploration
of the digital, multi-wavelength universe that is resides in the
astronomical data archives. The VO is already an operational research infrastructure, as demonstrated by the growing number of VO-papers published in the last years ($>$ 50 since 2009)\footnote{http://www.euro-vo.org/pub/fc/papers.html}. In this work we have made use of VO tools to benefit from an easy data access and analysis.

We present here the identification of 31 BDs (25 known and 6 strong candidates not previously reported in the literature) identified in the sky area in common to the WISE Preliminary Release and the 2MASS Point Source and SDSS Data Release 7 catalogues. In Sect. 2 we describe the methodology devised to search for T dwarfs. In Sect. 3 we derive proper motions, effective temperatures ($T_\mathrm{eff}$), spectral types, and distances for our candidates. Finally, we summarize conclusions in Sect. 4.

\section{Candidate selection}

We built a VO-workflow with the STILTS\footnote{http://www.star.bris.ac.uk/$\sim$mbt/stilts/} scripting capabilities  with  criteria that must accomplish our preliminary list of candidates. The workflow consisted in the following steps:

\begin{itemize}

\item Cross-match. To avoid memory overflow problems associated to the data processing of large volumes of data, we divided the common WISE-2MASS-SDSS sky into overlapping circular regions of 30\arcmin. After different tests, we adopted a matching radius of 20\arcsec$ $ to ensure that objects with high proper motion were not left out while at the same time avoiding an unmanageable number of false positives. Only the closest 2MASS and SDSS counterparts to each WISE source were considered. \\

\item Filters on photometry. \\

\begin{itemize}

\item[$\bullet$] Xflg=0, Aflg=0, to avoid sources flagged in 2MASS as minor planets or contaminated by nearby extended sources.
\item[$\bullet$] 9 $<$ \textit{J}, to exclude of too bright sources that typically have associated large uncertainties.
\item[$\bullet$] Qflg(\textbf{\textit{J}}) $\neq$ "U", to discard sources with upper limits in the \textit{J} band.
\item[$\bullet$] \textit{g'} $>$ 22.2, \textit{r'} $>$ 22.2 (SDSS limiting magnitudes), because T dwarfs should not be detected in these bands.
\item[$\bullet$] \textit{i'} $>$ 21.3 (SDSS limiting magnitude) or (\textit{i'}-\textit{z'}) $\geq$ 3 (Hawley et al. 2002).
\item[$\bullet$] \textit{z'} $<$ 20.5 (SDSS  limiting magnitude)
\item [$\bullet$] (W1-W2)$>$0.5 , W2$<$15.5, (W2-W3)$<$ 2.5. The W1 and W2 WISE bands were specifically designed to distinguish T dwarfs from background
sources. W1 includes much of the CH$_{4}$ fundamental absorption band and W2 includes the pseudo-continuum peak that is present for all  objects cooler than 3000 K (\cite{2003ApJ...596..587B}; \cite{2004AJ....127.3516G}). \cite{2011ApJ...726...30M}  suggest a (W1-W2) $>$ 1.8 for cool T dwarfs. Because we search for T dwarfs at all spectral subtypes, we relaxed this criterion up to the stellar boundary, (W1-W2)$>$ 0.5. A side effect of this approach was the identification of four new L-type candidates. The (W2-W3)$<$ 2.5 criterion was selected to exclude extragalactic sources (\cite{2010AJ....140.1868W}). W2=15.5 indicates the limiting magnitude (5$\sigma$) in the W2 band.
\item[$\bullet$] (\textit{J}-W2)$>$1.8 and (\textit{z'}-\textit{J})$>$2.5 as suggested by \cite{2011ApJ...726...30M} and \cite{2002AJ....123.3409H}, respectively.\\
\end{itemize}

\item Filter on astrometry and proper motions (PMs). \\

\begin{itemize}
\item[$\bullet$] The angular separation between the WISE source and the 2MASS counterpart must be greater than 1\arcsec$ $ (WISE 2$\sigma$ astrometric accuracy). Assuming a minimum difference epoch between 2MASS and WISE of nine years (Feb 2001 - Jan 2010), this criterion implies that we are loosing objects with PMs lower than 0.11\arcsec/yr. On the other hand, the 20\arcsec $ $ matching radius will limit the detection of T dwarfs with PM higher than 1.53\arcsec/yr $ $given the maximum difference epoch between 2MASS and WISE (Jun 1997 - Apr 2010).
\item[$\bullet$]There must not be counterparts in the source table of the Supercosmos Science Archive (\cite{2001MNRAS.326.1279H}) at less than 1\arcsec$ $ from the WISE source.
\item[$\bullet$] The differences in the PMs derived from WISE-2MASS and WISE-SDSS must be less than 50\%.
\end{itemize}
\end{itemize}
This query returned 138 candidates. They were visually inspected using the scripting capabilities of Aladin\footnote{http://aladin.u-strasbg.fr/} (\cite{2000A&AS..143...33B}). Aladin is a VO-compliant software that allows users to visualize and analyze digitised astronomical images, and superimpose entries from astronomical catalogues or databases available from VO services. Sources from DENIS, 2MASS, SDSS, Supercosmos and UKIDSS as well as images from UKIDSS, 2MASS and SDSS were used in the analysis. One hundred and seven candidates were rejected due to several reasons (artifacts, presence of a nearby bright star, etc.), the most likely being the mismatch between WISE sources and 2MASS/SDSS counterparts: due to the different depth of the surveys it may happen that a 2MASS source is wrongly associated to a SDSS source without a real counterpart in 2MASS. Thirty one candidates passed the visual inspection. Twenty three sources were identified as known BDs (7 L- and 16 T-dwarfs) in a cross-match with SIMBAD\footnote{http://simbad.u-strasbg.fr/simbad/} and the DwarfArchives, providing an independent sanity check of our selection method. Two out of the remaining eight
candidates were published during the preparation of this letter: WISE J1625+1528 (\cite{2011arXiv1106.3105D}) and WISE J1627+3255
(\cite{2011AJ....142...57G}). Multiwavelength charts and the astrometric and photometric information of the six newly discovered BD candidates are given
in Fig.~\ref{secuencia} and Table~\ref{tabla}, respectively. Figure ~\ref{diagrama} compares in a color-color diagram our six candidates with a sample of
known L and T dwarfs (\cite{2011ApJ...726...30M}, \cite{2011A&A...532L...5S}, \cite{2011ApJ...735..116B}). All our candidates have  W1,W2 signal-to-noise ratios (SNR) $\geq$ 10 and can be considered as point
sources (w2rchi2 $\leq$ 3). They were also searched for common PM companions within a radius of 10\arcmin$ $ in the WISE and 2MASS catalogues but none was found.

Finally, we assessed the efficiency of our search by estimating the false negative rate (number of known T dwarfs that were not rediscovered in the analysis). Thirty six T dwarfs catalogued in the DwarfArchives were not found with our methodology. Most of them (33) were too faint to be detected in 2MASS and the remaining three did not meet some of the photometric criteria explained above.

\begin{figure*}[h!]
   \centering
   \includegraphics[width=19cm]{aladin_composicion.ps}
   \caption{Four-band charts showing the area around the six new BD candidates. The fields are 1.5\arcmin $\times$ 1.1\arcmin oriented with North up and East to the left. For each source, the SDSS(\textit{z'}), 2MASS(\textit{J}), WISE(W1) images are centered on the WISE(W2) coordinates. Fast motion between 2MASS and WISE epochs is clearly visible in WISE J0920+4538.}

   \label{secuencia}
\end{figure*}

\begin{table*}
\centering
\caption{Astrometry, photometry and physical parameters for the six new BDs candidates. See Sect. 3 for a detailed description on how effective temperatures, spectral types, absolute magnitudes and distances were obtained. \textit{J} and W2 absolute magnitudes were only computed for L and T dwarfs respectively. Also, the uncertainty in W3 is not provided for four of the candidates. According to the WISE source catalogue user's guide, this happens if the W3 profile-fit magnitude is a 95\% confidence upper limit or if the source is not measurable.}
\label{tabla}
\scriptsize

\hspace{-1cm}
\begin{tabular}{ccccccc}
\noalign{\smallskip}
\hline
\noalign{\smallskip}
Parameter		& WISE J0802+2527	& WISEJ0820+2632   &  WISE J0821+1443  	      & WISE J0830+4837   & WISE J0838+1511 & WISE J0920+4538 \\
\noalign{\smallskip}
\hline
\noalign{\smallskip}
RA(J2000)        	& 08 02 02.86 & 08 20 07.35 & 08 21 31.63 & 08 30 41.67 & 08 38 11.45 &  09 20 55.41  \\
Dec(J2000)       	&$+$25 27 18.5&$+$26 32 28.5&$+$14 43 19.4&$+$48 37 15.0&$+$15 11 15.1&$+$45 38 56.4    \\
$\mu_{\alpha}cos\delta$ ($\arcsec$/year)&   -0.011$\pm$0.026       &  0.087$\pm$0.042             &  -0.115$\pm$0.035    &  -0.125$\pm$0.026   &   -0.121$\pm$0.031     &  -0.075 $\pm$0.010      \\
$\mu_{\delta}$ ($\arcsec$/year)  	&  -0.125$\pm$0.055        &  -0.208$\pm$0.062             &  -0.297$\pm$0.054  &  -0.058$\pm$0.050    &   -0.032$\pm$0.051    &  -0.833$\pm$0.045       \\
\textit{z'}	         & 19.238$\pm$0.069 & 19.514$\pm$0.067 & 19.732$\pm$0.118 & 19.905$\pm$0.153 &19.907$\pm$0.142 &21.002$\pm$0.658 \\
Y (UKIDSS)               & 17.914$\pm$0.012 & 18.145$\pm$0.029 &                  &                  &                 &                 \\
\textit{J} (2MASS) 	 & 16.713$\pm$0.157 & 16.988$\pm$0.222 & 16.825$\pm$0.155 & 17.275$\pm$0.199 &16.645$\pm$0.159 &15.223$\pm$0.052 \\
\textit{J} (UKIDSS)      & 16.626$\pm$0.010 & 17.036$\pm$0.015 &                  &                  &                 &                  \\
\textit{H} (2MASS)       & 15.797$\pm$0.154 & 15.995$\pm$0.221 & 16.515$\pm$0.240 & 16.279$\pm$0.176 &16.207$\pm$0.172 &14.164$\pm$0.054 \\
\textit{H} (UKIDSS)      & 15.937$\pm$0.019 & 16.487$\pm$0.038 &                  &                  &                 &                  \\
\textit{K$_{s}$} (2MASS) & 15.397$\pm$0.151 & 15.881$\pm$0.251 &                  & 15.596$\pm$0.165 &                 &13.728$\pm$0.050 \\
\textit{K$_{s}$} (UKIDSS)& 15.248$\pm$0.016 & 15.883$\pm$0.028 &                  &                  &                 &                  \\
W1                       & 14.758$\pm$0.037 & 15.578$\pm$0.057 & 16.438$\pm$0.114 & 14.620$\pm$0.034 & 15.712$\pm$0.068&13.059$\pm$0.025 \\
W2                       & 14.256$\pm$0.060 & 15.009$\pm$0.099 & 14.283$\pm$0.064 & 14.000$\pm$0.047 & 14.568$\pm$0.087&12.391$\pm$0.026  \\
W3                       & 11.965$\pm$---   & 12.658$\pm$---   & 12.572$\pm$---   & 12.828$\pm$0.539 & 12.285$\pm$---  &11.063$\pm$0.115 \\
W1$-$W2 & 0.502 &0.569 & 2.155 & 0.62 & 1.144 & 0.668 \\
W2$-$W3 & 2.291 & 2.351& 1.711 & 1.172 & 2.283 & 1.328\\

$T_\mathrm{eff}$  (VOSA) (K)& 1800 & 2000 & 700 & 1700 & 900 & 1700 \\
SpT ($T_\mathrm{eff}$) & L3-L4 & L2-L3 & T8 & L4-L5 & T6-T7 & L4-L5 \\
$M_{J}$   & 13.08 & 12.48 &        & 13.31 &         & 13.31        \\
$M_{W2}$           &       &       & 13.47  &       &  13.07  &        \\
d (pc)                                 & 53 & 80 & 15  & 62 &  20  &  24\\
\noalign{\smallskip}
\hline
\noalign{\smallskip}
\noalign{\smallskip}
\hline
\noalign{\smallskip}
\end{tabular}

\end{table*}
\begin{figure}[t!]
   \centering
   \includegraphics[width=7cm]{DiagramaJHvsW1W2.ps}
   \caption{Colour-colour diagram of our six candidates (black full circles) compared with a sample of known L  and T dwarfs (triangles and circles,
respectively) observed with WISE (Mainzer et al. 2011; \cite{2011A&A...532L...5S}; Burgasser et al. 2011.)}
   \label{diagrama}
\end{figure}
\begin{figure*}[h!]
    \label{VOSA}
    \centering
    \includegraphics[width=11cm]{fit_vosa2.eps}
    \caption{VOSA SED fitting for our six BD candidates. Catalogue and synthetic photometric points are represented in red and blue, respectively. Points
not considered in the fitting are shown in black. Spectral bands are those listed in Table 1. The theoretical spectrum that best fits (and from which the synthetic photometry was computed by
convolving with the corresponding filter profiles) is also plotted in grey.}
 \end{figure*}

\section{Physical parameters}

\subsection{Proper motions}

The long baseline between WISE and 2MASS/SDSS, the excellent astrometry of these catalogues (less than 0.5\arcsec$ $ with respect to 2MASS in the case of WISE (\cite{2010AJ....140.1868W})), together with the fact that brown dwarfs are nearby objects and therefore typically show high proper motions facilitates the measurement of reliable angular separations between the WISE sources and their respective counterparts in the 2MASS and SDSS catalogues.

Proper motions were computed using a linear least-squares fit to the coordinates given in the WISE, 2MASS and SDSS catalogues and weighted by their astrometric errors. No correction for parallactic motion was considered. The small separation between the 2MASS and SDSS positions (less than 1$\arcsec$ for five of our candidate BDs) and the lack of 2MASS $K_{s}$ photometry for two of them are the main reasons for their omission in previous 2MASS/SDSS-based PM searches.

\subsection{Effective temperature and spectral types}
Effective temperatures for our candidate BDs were obtained from the spectral energy distribution (SED) $\chi^{2}$ fitting between the observed photometry and a suite of collections of theoretical models. The temperature determination was carried out using VOSA\footnote{http://svo.cab.inta-csic.es/theory/vosa/} (Virtual Observatory SED Analyzer, \cite{2008A&A...492..277B}). VOSA is a VO-tool designed to query several photometric catalogues that are accessible through VO services as well as VO-compliant theoretical models and perform a statistical test to determine which model reproduces the observed data best. $T_\mathrm{eff}$ is then estimated from the best fit. Two different model collections were used in our analysis: DUSTY (\cite{2001ApJ...556..357A}) and COND (\cite{2000ApJ...542..464C}), which have $T_\mathrm{eff}$ models ranging from 4000 to 100K and 3900 to 500K in 100K steps, respectively. Effective temperatures obtained with VOSA are shown in Fig. \ref{VOSA} and Table \ref{tabla}.

Spectral types were then estimated following the effective temperature--spectral type relation given in \cite{2005ARA&A..43..195K}. The 23 BDs found in our analysis  with spectral types previously reported in the literature were used as benchmarks to calibrate the method. Typical errors of $\sim$0.5 spectral types were found. The main contribution to the error budget comes from the degeneracy in the $T_\mathrm{eff}$-- spectral. type relation between L4 and T4.

\subsection{Distances}

Distances for the L dwarfs were calculated using the absolute \textit{J} magnitude (2MASS) -- spectral type relation given in
\cite{2003AJ....126.2421C}, whereas the distances to the T dwarfs were estimated using the absolute W2 magnitude (WISE) -- spectral type relation derived
by \cite{2011ApJ...735..116B}. The calculated values for the distances are given in Table \ref{tabla}. The uncertainties in the derived distances are dominated by the uncertainty in the spectral type because the errors in the 2MASS J photometry are typically only 0.1-0.2 mag, whereas the errors in the spectral type are 0.5 types. This leads to a 30\% uncertainty on distance, assuming all our candidates are single. An unresolved multiplicity would imply an underestimate of distance.

\section{Conclusions}

Taking advantage of VO tools, we reported the discovery of six BD candidates in the region of the sky in common to 2MASS Point Source, SDSS (Data Release 7) and the WISE Preliminary Release catalogues. The six candidates are clearly visible in the WISE W1 and W2 bands with SNR $\geq$ 10 and show physical parameters typical of L and T-type objects. The number of new candidates is remarkable, considering that 2MASS has been extensively searched for ultracool dwarfs.

These results clearly show how new surveys and the use of VO tools can help to mine older surveys: all but one of our BD candidates are very close to the 2MASS limiting magnitude, which puts them beyond the limits of previous, much shallower 2MASS-based BD searches (e.g. J$<$16, \cite{2007AJ....134.1162L}). Also remarkable is the fact that the methodology used in this paper is not limited to brown dwarfs but can be easily extrapolated to searches for other rare objects (e.g. high-z quasars).

Finally, our work clearly  demonstrates the suitability of exploiting WISE data following a VO methodology and increases the expectations of building an accurate census of substellar objects in the solar vicinity.

\begin{acknowledgements}
 This research has made use of the WISE data hosted at the NASA/ IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This publication has also made use of the SIMBAD, Vizier and Aladin services, operated at CDS, Strasbourg, France. Our research has benefitted from the M, L, and T dwarf compendium housed at DwarfArchives.org and maintained by Chris Gelino, Davy Kirkpatrick, and Adam Burgasser. Effective temperatures were estimated using VOSA, a VO-tool developed under the Spanish Virtual Observatory project supported from the Spanish MICINN through grant AyA2008-02156.
\end{acknowledgements}

Solar and Stellar Astrophysics (astro-ph.SR)


Bulge Formation

Abstract: The currently discussed theories of bulge formation are reviewed, including the primordial scenario, where bulges form rapidly and then accrete disks, the secular scenario, where bulges are formed by dynamical evolution of disks through bars and galaxy interactions, and some combinations of both, where formation of bulges and disks are more continuous and interleaved. The various scenarios make specific predictions about the relative masses, angular momenta, colours, metallicities of bulges relative to disks, and the bulge-to-disk ratio as a function of time. Dynamical processes relevant to the formation of bulges (bar instabilities, mergers) are described and tested against observed statistics. Current data suggest a dynamical feedback from gravitational instabilities in bulge and disk formation. It is very difficult to discriminate between the various scenarios from surveys at z=0 only, and observations at high redshift are presently the best hope for large progress.

\section{Introduction}

The problem of bulge formation is in rapid evolution:
not only many scenarios, dynamical processus, formation theories
have been proposed and studied, but there has been great
progress in observation of bulges: on ages, metallicities,
structures, etc... One can cite in particular:

\begin{itemize}

\item Detailed age study of individual globular clusters, with colour-magnitude
diagrams of individual stars, with the high spatial resolution of HST
(in the Milky Way bulge, and in Local Group galaxies)

\item More extinction-free studies of bulge structure through near-infrared
imaging, made recently possible at large scale (wide extragalactic
surveys, COBE for the Milky Way)

\item Morphological studies of galaxies at high redshift: this
has been extensively developped in this meeting, and is the
privileged tool to tackle galaxy evolution in situ.

\end{itemize}

There are very good recent reviews on the subject, namely in
proceedings of the STSci workshop held in 1998 (``When and how bulges
form ?'', ed. Carollo, Fergusson \& Wyse),  Renzini (1999), Silk \& Bouwens
(1999), Carlberg (1999) or the Annual Review on ``Galactic Bulges'' by
Wyse, Gilmore \& Franx (1997).

The main scenarios proposed for the fornation of bulges are:

\begin{itemize}

\item  Monolithic formation (Eggen, Lynden-Bell \& Sandage
 1962), or very early dissipative collapse at the beginning of galaxy
formation. This assumes that the gas experiences violent
3D star formation, so quickly that it had no time to
settle into a disk. This scenario was first proposed to explain
the almost spherical old metal-poor stellar halo
and the subsequent formation of disks with different thickness,
that were thought to be an age sequence following the
progressive gas settling. This is now proposed for
the central bulge, the disk being acquired later.

\item Secular dynamical evolution: the time-scale of these
processes are longer than the dynamical time (i.e. secular),
but however smaller than the Hubble time. They are due
to the dynamical interaction between the various components,
disk, bulge, halo. Gravitational instabilities, such as bars and
spirals, are able to transfer efficiently angular momentum, and
produce radial mass flows towards the center. Due to vertical
resonances, stars in the center  are elevated above the plane,
and contribute to bulge formation (Combes et al. 1990).

\item galaxy interactions, through merger and mass accretion.
It is well known that major mergers between spirals can
result in the formation of an elliptical galaxy (Toomre \&
Toomre 1972). Similarly, the accretion by a spiral galaxy of
a dwarf companion could contribute to the formation of
a spheroid at the center, assumimg that this minor merger
has not destroyed the disk.

\end{itemize}

In fact, all three of these main scenarii certainly occur, the
question is to estimate the relative role of each of them,
which is related to the time-scale of bulge formation.
All these processes may be included in the general frame
of hierachical galaxy formation. Disks are supposed to
form through gas cooling in a dark halo. Cooling runs
progressively from the center to the outer parts
(case of continuous gas infall), and produces an
inside-out formation of disks. Either this cooling is
first violent at the center, due to the short dynamical
time-scales, and the absence of a stabilising heavy stellar
object, and a spheroid can form first; or secular
dynamical evolution could afterwards transfer angular
momentum, and bring mass to the center.
 In all scenarii, the speed of evolution and rate of
star formation depends strongly on environment.
Galaxy interactions are both responsible of the merger
scenario, and also trigger or boost bar formation
and secular evolution.
Bulges at the center of galaxies accumulate all stars
formed, and in any case they are expected to be older than
the outer disk, and mainly with much more scattered
properties (ages, abundances).

\smallskip

We will first briefly review the constraints and clues
brought by observations related to bulge formation,
and then examine each scenario respectively.

\section{Clues from observations}

Observations of our own bulge is impeded by
dust extinction, confusion (crowding), contamination
by foreground stars, that make data on the Milky Way very uncertain.
There exist in the literature  a certain number of
prejudices, like the idea that ``bulges are old and metal-rich,
and small versions of ellipticals'', that are not today
completely confirmed: the reality is not so simple, as clearly
reviewed by Wyse et al (1997).

\subsection{ Metallicity}

In the Milky Way, the mean metallicity of the bulge
is the same as that of the disk in the solar neighborhood,
but in contrast to the disk, the bulge has a very wide scatter.
This characteristic allows the metallicity distribution of
the bulge to be explained by a closed box chemical model,
contrary to the disk: the latter has a very narrow distribution
of metallicities, that gives rise to the well-known G-dwarf problem
(its solution requires other processes, like gas infall, etc..).
In the Andromeda bulge also, super metal-rich globular
clusters are the exception (Jablonka et al 1998).

In the Milky Way bulge, there is no correlation between age, color, abundance
and kinematics, which could have given clues about the origins
(e.g. Rich 1997).
There does not seem to have been a starburst, since there
is no excess of $\alpha$ elements (Mc William \& Rich 1994).

\subsection {Age}

As for ages, however, we should not make the confusion:
{\it the bulge has not necessarily the age of its stars}.
In particular the bulge could have formed recently from
old stars. It is true namely if the bulge has formed from bars.
 Thanks to the clear overall structure of the MW given
by COBE-DIRBE, we know directly now that a bar is present,
 while it has long be assumed from gas kinematics only
(e.g. Peters 1975). The peanut/boxy bulge of the MW
has asymmetries due to a bar seen in perspective
(e.g. Blitz \& Spergel 1991, Zhao et al 1996). It is even difficult to distinguish
what is bar and what is bulge (Kuijken 1996). Since strong bars
are only a transient phase (see below) it is easy to extrapolate
how the stars presently in the bar will form the bulge.

\begin{figure}
\psfig{figure=bulge99_f1.ps,width=12cm,bbllx=1cm,bblly=5cm,bburx=19cm,bbury=20cm,angle=0}
\caption{The $B-K$ color in surface-brightness of the central disk, versus
that of the bulge, from de Jong (1996). The various symbols indicate morphological
types, as described at the bottom left. The dashed line means equality.   }
\label{fig1}
\end{figure}

If the bulges are thought to be older than disks, it might
be due to the implicit average over the whole disk.
In fact, when bulges and inner disks are compared, they
have the same integrated colors (Peletier \& Balcells 1996,
de Jong 1996, see fig \ref{fig1}). But there exist clear
   radial gradients of colors and metallicities.
Spiral galaxies become bluer with increasing radius.
Moreover, colors correlate well with surface brightness.
These colors and their gradients can best be explained
with history of star formation, and dust reddening is
not dominant, as shown by de Jong (1996):
outer parts are younger than central regions.
These observations support an inside-out galaxy formation.
Alternatively, dynamical secular evolution can
also account for these observations, through angular momentum transfer
(cf Tsujimoto et al 1995).

\subsection{ Are bulges similar or not to giant ellipticals ?}

Bulges are spheroids that follow the same fundamental plane
as elliptical galaxies. They also follow the same
luminosity-metallicity relation (Jablonka et al 1996).
But they are oblate rotaters (flattened by rotation),
while giant ellipticals are not rotating. However,
 low-luminosity ellipticals also rotate (Davies et al 1983).
There might be a continuity between bulges, lenticulars
and ellipticals. Most of the latters have been observed
with compact disks, e.g.  Rix \& White 1992, Rix et al. 1999).
 There is also continuity in the light profiles.
They can be fitted as an exponential of
$r^{1/n}$, while $n$ is a function of luminosity
($n=$4 for Ellipticals, $\infty$ for Bulges, Andredakis et al 1995).

\subsection{ High-redshift galaxies }

From the CFRS and LDSS, surveying galaxies at redshifts
below 1, there is very little evolution of the
luminosity function of red galaxies
(Lilly et al. 1995, 1998). Their
disk scale-length stays constant, up to $z=1$ (which
could mean either no-evolution, or stationary merging).
On the contrary, there is substantial evolution in the blue
objects. Does it mean that a
large fraction of bulges/spheroids are already formed?
(or mean age of stars $\propto$ age of Universe).
From HST images of the CFRS galaxies,
Shade et al (1995, 96) conclude that
red galaxies have large bulge-to-disks (B/D) ratios, and that
most blue galaxies are interacting (also forming bulges, nucleated
galaxies).

Steidel et al (1996) from their sample of z $\sim$ 3 galaxies
find in average objects smaller than today
(see also Bouwens et al 1998).
If violent starbursts are forming at high redshift, like the
ultra-luminous IR galaxies observed nearby, they should
dominate the sub-mm surveys. Already a large
fraction of the cosmic IR background has been resolved
into sources (Hughes et al. 1998). If these objects are
forming the spheroids (bulges and ellipticals), then their
epoch of formation is relatively recent $z < 2$ (Lilly et al. 1999).

From the high spatial resolution of HST, it becomes now
possible to study the morphology of high-redshift galaxies,
and address the evolution of the Hubble sequence.
At least, it is possible to determine the concentration
(C) and asymmetry (A) to classify galaxies (Abraham et al 1996, 99).
If there is a consensus among the various studies, it is
on the considerable increase of perturbed and interacting
galaxies ($\sim$ 40\% objects interacting).
But on the concentration, or bulge-to-disk ratios,
there is no unanimity. Abraham et al (1998) do not
find evolution in the B/D ratio, and
tend to favour monolithic bulge formation.
Marleau \& Simard (1998) find on the contrary many
more disk-like faint galaxies in the Hubble Deep Fields, i.e.
that there are less objects with high bulge-to-disk ratios
at high redshift (with respect to $z=0$).

\section{Monolithic dissipative formation}

 This kind of violent and intense star formation
is different to what is seen today in galactic disks: it
occurs in three-dimension (Elmegreen 1999).
In the deep potential well of a giant galaxy center,
almost entirely due to the dark halo, there cannot be
SF self-regulation by blow-out (supernovae, winds, pressure)
as in disks or dwarfs today (Meurer et al 1997).
Most of the time, in the ultra-luminous galaxies
of the nearby universe, the starburst occurs in nuclear rings or disks
(Downes \& Solomon 1998). This could be different in the
early universe, since there are no heavy disks formed.
The threshold for SF in 3D can be estimated from the virial theorem
(Spitzer 1942) as a critical central volumic gas density
$\rho_c = \rho_{vir} /(1+\beta/2)$ for gaseous potential
$\beta GM_{gas}/R$. The star formation rate SFR is then
$SFR \propto \rho_{vir}/t_{dyn} \propto  \rho_{vir}^{3/2}$
(an equivalent of Schmidt law). In localised dense regions,
the SFR could be very large, and an
extremely clumpy distribution of stars
is expected. Simulations illustrating this process
have been done by Noguchi (1998).

\section{Secular evolution}

A wide series of N-body simulations, supported by
observations, have established the various phases
of this secular dynamical evolution
(e.g. Sellwood \& Wilkinson 1993, Pfenniger 1993,
Martinet 1995, Buta \& Combes 1996). They can
be described as:

\begin{itemize}

\item Gravitational instabilities spontaneously
form bars, spirals, which transfer angular-momentum
through the disk.
Both stars and gas lose momentum to the wave.
Radial gas flows in particular produce large
central mass concentration.

\item Due to vertical resonances, the bar thickens
in the center (box or peanut-shape, Combes et al 1990),
and contributes to bulge formation (see fig \ref{fig2}).
 This scenario is compatible with the observed
correlation between scale-lengths of bulges and disks
(Courteau et al 1996).

\item Sufficient mass accumulation in the center
(1-5\% of the total disk mass) perturbs the rotation curve,
changes the orbit precession rates, and
destroy the bar, by creating perpendicular orbits
$x2$. The galaxy has now become a hot stable system,
with a central spheroid.

\item Through gas infall,
the disk can become unstable again  and reform a bar
(with a different pattern speed)

\end{itemize}

\begin{figure}
\psfig{figure=bulge99_f2.ps,width=17cm,bbllx=15mm,bblly=5mm,bburx=10cm,bbury=27cm,angle=-90}
\caption{Edge-on projections of a bar obtained in N-body simulations, for
three orientations of the bar with respect to the line of sight (90, 45 and 0$^\circ$).}
\label{fig2}
\end{figure}

\section{Accretion and mergers}

Ellipticals can form by the mergers of two spirals
(e.g. NGC 7252 prototype); but the most general case
is the formation by the merger of several smaller objects.
That ellipticals have accreted several smaller companions
during their life is supported by the frequency of shells
(as much as 50\% according to Schweizer \& Seitzer 1988).
Bulges could similarly be the result of minor mergers;
even today, there exist numerous companions around giant spirals.

In hierachical cosmological models, it is easy to
compute analytically the history of dark halo formation,
and their merging rate (Kauffmann et al 1993, Baugh et al 1996).
But the fate of the dissipative visible matter is still
not well-known. According to some star formation and feedback
parameters that fit the colour-magnitude relation of ellipticals
in clusters,  Kauffmann \& Charlot (1998) have shown that
massive galaxies must have formed continuously, and
must have assembled only recently,
in order to fit observed redshift distributions of K-band luminosity at $z=1$.

If the bulge comes from accretion, observations of our own MW
put constraints on the metallicity of the objects accreted
(they should be relatively high).
This put strong limits on the fraction of the bulge that has been accreted
recently (Unavane et al 1996).

\section{Combination of all scenarii}

In a hierachical scenario, there are many uncertainties
about the physical processes controlling  the
baryonic matter. Feedback must prevent efficient conversion of
gas into stars, to represent observations. There remains much freedom
to determine it more quantitatively, as well as the distribution of
angular momentum, etc...

The simplest recipe is to assume that the angular momentum of
the matter is not redistributed, and that the gas cools at the
radius where the cooling time becomes shorter than the dynamical time
(Mo et al 1998; van den Bosch 1999). This represents gas infall as
cooling flows in dark haloes, and the formation of disks as an
inside-out process. There are however
many hypotheses possible: for instance,
that the mass of the disk is a fixed fraction of the halo mass,
or that the angular momentum of the disk is a fixed fraction of the halo
angular momentum; also that the
disk is exponential in shape, and is stable (which gives a
condition on specific momentum). The efficiency
of galaxy formation ($\epsilon_{gf}$) is also a free parameter, as
a function of redshift; it could vary proportionally to
the Hubble constant, to ensure the Tully-Fisher relation
(van den Bosch 1999).
The specific angular momentum of the galaxies are supposed
to be acquired by tides (Fall \& Efstathiou 1980).

These scenarios, coupled with sufficient feedback (due to star-formation)
can explain high-surface brightness galaxies (HSB) and even LSB,
if a range of angular momentum is assumed (Dalcanton et al 1997).
 The following constraints can be satisfied:
1) the Tully-Fisher scaling relation;
2) the density-morphology relation (Dressler 1980);
3) the surface-density/size relation $\mu$ - R$_d$
(generalisation of Kormendy 1977).

A general feature found by van den Bosch (1999), whatever the cooling/feedback
hypotheses, is that the present observed disks lie very close to their
stability limit (with respect to strong gravitational instabilities), which suggests
a self-regulated formation, through dynamical processes:
the disk alone is unstable until a sufficient bulge stabilises it.
The gas infall could then alternatively favor bulge or disk formation
according to the bulge-to-disk mass ratio it encounters at a given
time. This could explain the observed coupling between
bulges and disks.

\begin{figure}
\psfig{figure=bulge99_f3.ps,width=12cm,bbllx=1cm,bblly=45mm,bburx=17cm,bbury=24cm,angle=0}
\caption{ Comparison of the bulge-to-total ratios (B/T) between observations
(histograms) and model predictions (lines), from Bouwens et al (1999). At $z=0$,
the data are from Peletier \& Balcells (1996, PB) and de Jong (1996, DJ).
At high $z$ from Schade et al (1996) and Lilly et al (1996). In all frames,
the models are: {\bf I} = secular evolution; {\bf II} = simultaneous formation;
{\bf III} = early bulge formation. }
\label{fig3}
\end{figure}

Bouwens et al. (1999) have recently confronted several bulge formation scenarii
to the observations. At $z=0$, it appears almost impossible to distinguish
between early and late bulge formation (cf fig \ref{fig3}). At high redshift,
it becomes possible, essentially because a large fraction of the star formation
occurs at the observed $z$. For instance, in fig \ref{fig3} at top-left, the model
of  simultaneous formation (of disk and bulge) departs from the two others,
since star-formation makes the bulge momentarily brighter.

Could the star formation history (SFH) constrain the period
of formation of bulges and ellipticals?
Since spheroids contain 30\% (Schechter \& Dressler 1987),
or 66\% (Fukugita et al 1998) of the stars in the Universe,
they cannot form too early, according to the Madau et al. (1996)
SFH curve. However, this constraint is only a lower-limit
on the formation time,
since the spheroids could be recently formed from old stars.

\section{Conclusions}

Bulges are very heterogeneous structures, with
large scatter in colors, metallicities and ages.
At large luminosities, bulges become similar to elliptical
galaxies, at low luminosities, they resemble more
spiral disks (or bars). There exist objects with a whole
range of properties ensuring continuity between these
two extremes, suggesting evolution or the combination
of several processes.

Bulge formation is certainly a combination between
three main scenarii:

a fraction of bulges could have formed early (at first collapse);
then secular dynamical evolution enrich them; in parallel,
according to environment, accretion and minor mergers
contribute to raise their mass.

Big spheroids (SOs, Es) can only form through major/minor mergers.
Present disks have been (re)-formed recently.
To be more quantitative, and precise the time-scale of bulge
formation, a privileged solution is to observe and classify
the morphology of galaxies as a function of redshift, together
with detailed balance of the comoving volumic density of
galaxies of each type. This is necessary to avoid confusion between
a no-evolution model, at a given type, or a stationary evolution flow
along the Hubble sequence.

\begin{moriondbib}
\bibitem{} Abraham R., Tanvir N.R., Santiago B.X. et al. 1996, \mnras {279}{L47}
\bibitem{} Abraham R., Ellis R.S., Fabian A.C., Tanvir N.R., Glazebrook K. 1999,
\mnras {303}{641}
\bibitem{} Andredakis Y.C., Peletier R.F., Balcells M., 1995,
\mnras {275} {874}
\bibitem{} Baugh C.M., Cole S., Frenk C.S., 1996, \mnras {283}{1361}
\bibitem{} Blitz L., Spergel D.N., 1991, \apj {379}{631}
\bibitem{} Bouwens R.J., Cayon L., Silk J., 1999, {\it preprint} (astro-ph/9812193)
\bibitem{} Bouwens R.J., Broadhurst T., Silk J. 1998, \apj {506} {557}
\bibitem{} Buta R., Combes F. 1996, {\it Fundamental of Cosmic Physics}
 {\bf 17}, 95-282
\bibitem{} Carlberg R.G., 1999,  in "When and
    How do Bulges Form and Evolve?", ed. by C.M. Carollo,
H.C. Ferguson \& R.F.G.  Wyse (Cambridge University Press),
(astro-ph/9903373)
\bibitem{} Combes F., Debbasch F., Friedli D., Pfenniger D., 1990,
\aa {233} {82}
\bibitem{} Courteau S., de Jong R.S., Broeils A.H., 1996, \apj {457}{73}
\bibitem{} Dalcanton J.J., Spergel D.N., Summers F.J., 1997, \apj {482}{659}
\bibitem{} Davies R.L., Efstathiou G., Fall S.M., Illingworth G., Schechter P.L.,
1983, \apj {266} {41}
\bibitem{} de Jong R.S., 1996, \aa {313} {377}
\bibitem{} Downes D., Solomon P.M., 1998, \apj {507}{615}
\bibitem{} Dressler A., 1980, \apj {236}{351}
\bibitem{} Eggen O.J., Lynden-Bell D., Sandage A., 1962, \apj {136} {748}
\bibitem{} Elmegreen B.G., 1999, \apj {517} {May 20}
\bibitem{} Fall S.M., Efstathiou G., 1980, \mnras {193}{189}
\bibitem{} Fukugita M., Hogan C.J., Peebles P.J.E., 1998, \apj {503}{518}
\bibitem{} Hughes D.H., Serjeant S., Dunlop J. et al.: 1998, \nat {394} {241}
\bibitem{} Jablonka P., Martin P., Arimoto N., 1996, \aj {112} {1415}
\bibitem{}  Jablonka P., Bica, E., Bonatto C. et al., 1998, \aa {335} {867}
\bibitem{} Kauffmann G., Charlot S., 1998, \mnras {297} {L23}
\bibitem{} Kauffmann G., White S.D.M., Guiderdoni B., 1993, \mnras {264}{201}
\bibitem{} Kormendy J., 1977, \apj {218} {333}
\bibitem{} Kuijken K. 1996, in {\it Unsolved Problems of the Milky Way},
IAU Symp. 169, ed. L. Blitz \& P. Teuben, p. 71 (Kluwer)
\bibitem{} Lilly S.J., Tresse L., Hammer F., Crampton D., LeFevre O., 1995,
\apj {455}{108}
\bibitem{} Lilly S.J., LeFevre O., Hammer F., Crampton D., 1996,
\apj {460}{L1}
\bibitem{} Lilly S.J., Shade D., Ellis R. et al., 1998,  \apj {500}{75}
\bibitem{} Lilly S.J., Eales S.A., Gear W.K. et al. 1999, in "When and
    How do Bulges Form and Evolve?", ed. by C.M. Carollo,
H.C. Ferguson \& R.F.G.  Wyse (Cambridge University Press),
(astro-ph/9903157)
\bibitem{} Madau P., Ferguson H.C., Dickinson M.E. et al. 1996,
\mnras {283}{1388}
\bibitem{} Marleau F.R., Simard L., 1998, \apj {507}{585}
\bibitem{} Martinet L., 1995,  {\it Fundamental of Cosmic Physics}
 {\bf 15}, 341-440
\bibitem{} McWilliam A., Rich R.M., 1994, \apjs {91}{749}
\bibitem{} Meurer G.R., Heckman T.M., Lehnert M.D., Leitherer C., Lowenthal J.,
1997, \aj {114} {54}
\bibitem{} Mo H.J., Mao S., White S.D.M., 1998, \mnras {295}{319}
\bibitem{} Noguchi M., 1999, {\it preprint}  (astro-ph/9806355)
\bibitem{} Peletier R.F., Balcells M., 1996, \aj {111} {2238}
\bibitem{} Peters W.L., 1975, \apj {196} {617}
\bibitem{} Pfenniger D., 1993, in {\it Physics of Nearby Galaxies:
Nature or Nurture}, ed. T.X. Thuan, C. Balkowski, J. Tran Thanh Van,
Ed. Frontieres, p. 519
\bibitem{} Renzini A., 1999,  in "When and
    How do Bulges Form and Evolve?", ed. by C.M. Carollo,
H.C. Ferguson \& R.F.G.  Wyse (Cambridge University Press),
(astro-ph/9902108)
\bibitem{} Rich M., 1997, in {\it ``The Central Regions of the Galaxy
and Galaxies''}, IAU Symp. 184, Kyoto (Kluwer Pub)
\bibitem{} Rix H-W., White S.D.M., 1992, \mnras {254} {389}
\bibitem{} Rix H-W., Carollo C.M., Freeman K., 1999, \apj {513} {L25}
\bibitem{} Schade D., Lilly S.J., Crampton D. et al., 1995, \apj {451}{L1}
\bibitem{} Schade D., Crampton D., Hammer F., LeFevre O., Lilly S.J., 1996,
\mnras {278}{95}
\bibitem{} Schechter P.L., Dressler A., 1987, \aj {94}{563}
\bibitem{} Schweizer F., Seitzer P.  1988, \apj {328}{88}
\bibitem{} Sellwood J.A., Wilkinson A., 1993, {\it Rep. Prog. Phys.}
{\bf 56}, 173
\bibitem{} Silk, J., Bouwens R.J., 1999, in {\it Proceedings of
``Galaxy Evolution: Connecting the distant Universe with
the local fossil record"}, Rencontres de Meudon, 21-25 sep 98,
ed. M. Spite, Kluwer (astro-ph/9812057)
\bibitem{} Spitzer L., 1942, \apj {95}{329}
\bibitem{} Steidel  C.C., Giavalisco M., Dickinson M., Adelberger K.L.
            1996, \aj  {112}  {352}
\bibitem{} Toomre A., Toomre J., 1972, \apj {178} {623}
\bibitem{} Tsujimoto T., Yoshii Y., Nomoto K., Shigeyama T.,
1995, \aa {302} {704}
\bibitem{} Unavane M., Wyse R.F.G., Gilmore G., 1996, \mnras {278}{727}
\bibitem{} van den Bosch F.C., 1999, \apj {507} {601}
\bibitem{} Wyse R.F.G., Gilmore G., Franx M., 1997,
{\it Ann. Rev. Astron. Astrophys.} {\bf 35},  637
\bibitem{} Zhao H., Rich R.M., Spergel D.N., 1996 \mnras {282} {175}
\end{moriondbib}
\vfill
\end{document}

Bulge Formation