Blazars and Gamma Ray Bursts

Abstract: Blazars and Gamma Ray Bursts (GRBs) are the fastest objects known so far. The radiation we see from these sources originates in a jet of similar aperture angle, and we think it is the result of the conversion of some of the jet kinetic energy into random motion of the emitting particles. Mechanisms for producing, collimating and accelerating the jets in these sources are uncertain, and it is fruitful to compare the characteristics of both class of sources in search of enlightening similarities. I discuss some general characteristics of blazars and GRBs such as the power of their jets compared with what they can extract through accretion, and the dissipation mechanism operating in the jets of both classes of sources. In both classes, there is a well defined trend between the bolometric power and the frequency at which this power is mainly emitted, but blazars are "redder when brighter", while GRBs are "bluer when brighter". Finally, I discuss some recent exciting prospects to use blazars to put constraints on the cosmic IR-Optical-UV backgrounds, and to use GRBs as standard candles to measure the Universe.

\section{Introduction}

It makes sense to compare blazars and Gamma Ray Bursts (GRBs),
looking for similarities and differences.
The bulk velocity in these systems differs from the light speed
by less than 0.5 per cent in blazars, and probably less than
0.005 per cent in GRBs.
In blazars, jet must be accelerated  without ``wasting" much of their
power in radiation, while in GRBs the high initial internal
optical depth prevents us to
study in detail the starting phases of the jet acceleration process.
The central engine in both systems is not completely understood,
but the main ingredients should be the same: accretion onto a
black hole that rapidly spins.

Consider also that GRBs are {\it not} explosive events.
If their duration is directly linked to the activity of
the accretion, then they last at least for $10^6$ dynamical
times, and a factor 100 more if the "post--prompt" flares
observed by Swift are produced by the central engine.
Their duration is then analogous to $\sim$3000 years
in the life of an AGN with a $10^8$ solar mass black hole.

Both classes of objects emit most of their radiation in
the $\gamma$--ray band: the radiation processes for blazars
and for the afterglow of GRBs are the same:
synchrotron and inverse Compton process, while some uncertainty
remains for origin of the prompt emission of GRBs.

Both classes of objects seems to form a sequence, as a function
of their power, which is linked to the overall spectral energy
distribution (SED), but the two sequences are opposite to each other:
blazar are bluer when dimmer, while bluer GRB are brighter.

Since the Universe is transparent to hard X rays and more opaque
to high--energy (GeV and TeV) $\gamma$ rays (since $\gamma$ rays
can be absorbed in $\gamma$--$\gamma \to$ e$^\pm$ collisions
with IR, optical and UV background photons), we can use the fact that
blazars and GRBs are strong X--ray and $\gamma$--ray emitters to
study the far Universe.
This is particularly true for blazars in the forthcoming GLAST era,
when it will be possible to study the optical--UV background
in a way similar to what done today for the IR background
using TeV emitting BL Lacs.
Finally, and more importantly,
the GRB sequence is so tight that it can be used to standardize the
GRB energetics and to use GRB as standard candles up to redshifts
unreachable by type Ia supernovae.

\section{Bulk Lorentz factors}

\subsection{$\Gamma=10$ in blazars?}

The most direct estimates of the bulk Lorentz factor $\Gamma$
in blazars come from VLBI observations of superluminal speeds
$\beta_{\rm app}$
(see e.g. Jorstadt et al. 2001 and Marscher, these proceedings).
One derives (since $\Gamma\ge \beta_{\rm app}$) values
between 5 and 45 for FSRQ (powerful blazars with broad emission lines)
and much less for low power, TeV emitting BL Lacs
(see e.g. Edwards \& Piner 2002; Piner \& Edwards 2004; Giroletti et al., 2004).
In these very same sources, however, modeling of the SED
requires very large $\Gamma$--factors, especially when the TeV
flux is de--absorbed (see e.g. Krawczynski, Coppi \& Aharonian 2002;
Konopelko et al. 2003;  Katarzynski et al. 2006).
Furthermore, the intra--day variability in the radio band,
if intrinsic, implies $\Gamma\sim 100$ (see Wagner \& Witzel 1995
for a review).
Even if due to interstellar scintillation, as it seems,
the implied small angular size of the source still puts
severe limits on its bulk Lorentz factor.
Theoretical limits on the maximum possible bulk Lorentz factor
($\Gamma<100$) have been derived by Begelman, Rees \& Sikora (1994),
while Ghisellini \& Celotti (2002) argued that the jet minimizes
its total power for a specific value of its bulk Lorentz factor.
For a specific source, the jet bulk Lorentz factor can vary
(see Marscher, this proceedings), bringing some support to the
``internal shock scenario" (see below).
There is some evidence and theoretical support to the idea that the
jet may be structured, having a fast spine
surrounded by a slower layer (with $\Gamma$ greater than $\sim$3
to avoid to see the counterjet; Laing 1993; Chiaberge et al. 2001;
Ghisellini, Tavecchio \& Chiaberge 2005).

Summarizing, there are some evidences that the typical values of
$\Gamma$ for blazars are larger then what thought previously, and
that different blobs in the same jet can have different $\Gamma$.
Also, the values of $\Gamma$ in low power jets is still an issue,
due to the largely different values derived from spectral modeling
(on scales of 0.01--0.1 parsecs) and VLBI data ($\sim$ 1 pc).

\subsection{$\Gamma=100$ in GRBs?}

The most common argument used to demonstrate the need of large
bulk Lorentz factor in GRBs is the compactness argument.
This combines the most rapid variability timescales
(around 1--10 ms) with the observation, in a few bursts,
of GeV emission during the prompt phase (see e.g. the review by
Fishman \& Meegan 1995).
If taken together, these two facts imply that the source,
if it is not moving relativistically, would be opaque
to the $\gamma$--$\gamma \to$ e$^\pm$ process, and no
$\gamma$ rays would be seen.
The required $\Gamma$ to let the source be transparent
to the photon--photon process is $\Gamma>100$.
On the other hand the GeV emission may originate in the
afterglow, not during the prompt, and according to the
time analysis of Beloborodov, Stern \& Svensson (1998; 2000)
the power spectrum has a break for timescales shorter than 0.1--1
second, so the variability on 1--10 ms may not be indicative of
the entire size of the emitting region.
Therefore the limit $\Gamma>100$ becomes questionable,
and one of the key issues that AGILE and GLAST are called
to clarify is the presence of high energy emission
that can clearly be associated with the prompt phase of GRBs.
The other arguments to invoke a large $\Gamma$--factor
are theoretical and are somewhat model--dependent, being
linked to the amount of baryon loading and the requirement
that the internal shocks thought to dissipate the carried
kinetic energy are more efficient increasing the contrast
between the $\Gamma$--factors of the colliding shells.

\section{Power}

The power of the jets in blazars is not known precisely,
although several ways to estimate it have been proposed.
The limiting factor are the still unknown matter
content of jet, the still unknown energy contained
in the baryon component of the radio lobes, and
also the fact that the radiation produced by the jet
is strongly boosted by beaming and then depends
on the not precisely known bulk Lorentz factor and
viewing angles.

For GRBs, the energetics of the received radiation is
accurately measured, but we have still a large uncertainty
concerning the efficiency of emitting this radiation.
The current estimates are $E_{\rm iso}=10^{50}$ -- $10^{54}$
ergs for the total emitted energy assuming isotropic emission.
The collimation corrected values are a factor 10--1000 less.

\subsection{Blazars}

Historically, the first estimate of the jet power
was derived by estimating the energy content
of the radio lobes in FR II radiogalaxies,
assuming an equipartition magnetic field and
cold protons (therefore the minimum energy conditions).
If one also estimates the lifetimes of these structures,
then one has the average jet energy required
to energize the lobes.
In this way one derives $10^{43}$--$10^{44}$ erg s$^{-1}$ for FR I
radiogalaxies and $10^{46}$--$10^{47}$ erg s$^{-1}$ for FR II
radiogalaxies and radio--loud quasars.
Rawling \& Saunders (1991) noted that this power
linearly correlates with the power in the narrow lines,
which have a factor 100 less power.
Since narrow line luminosities are also a factor $\sim 100$
less than the accretion disk luminosities,
the jet power and the disk luminosities are of the same order.

One can also calculate the power carried by the jet by
inferring its density through modeling the observed SED
and requiring that the jet carries at least the particles
and the magnetic field necessary to make the radiation we see.
This has been done on the pc scale by Celotti \& Fabian (1993),
on sub--pc scale (the $\gamma$--ray emitting zone,
see e.g. Ghisellini \& Celotti 2002),
and on the hundreds of kpc scale (the X--ray jets seen by Chandra,
see e.g. Celotti, Ghisellini \& Chiaberge 2001; Tavecchio et al. 2000).
These studies suggest large values of the power transported
by the jet and require
the presence of a dynamically dominating proton component
(see also arguments by Sikora \& Madejski 2000).

The superb angular resolution of Chandra
made possible to use two new methods:
the first is to calculate the energy needed to produce
the large scale jet X--ray emission (thought to be produced
by inverse Compton scattering with the microwave photons,
Tavecchio et al. 2000; Celotti, Ghisellini \& Chiaberge 2001), and
the second is to calculate the energy contained in the
``X--ray bubbles" i.e. cavities though to be inflated by the jet
(see e.g. Allen et al. 2006).

The emerging picture is that the powerful jets (as the ones in FR II
radiogalaxies) remains relativistic
up to hundreds of kpc, and that therefore they dissipate only a
small fraction of their power, which remains intact to power the radio lobes.
In these sources the jet power is of the same order of the accretion
luminosity, or even slightly larger occasionally (i.e. during flares).
Low power jets (as the ones in FR I radiogalaxies), even if they are born
with larger bulk Lorentz factor, can decelerate between the zone
of production of most of their radiation and the VLBI zone.
In these systems there are no or very weak broad emission lines, and
no sign of the thermal accretion disk luminosities, which is smaller
than the jet power.

A still hotly debated issue is what is the main energy carriers of jets.
Proton--electron, e$^\pm$ pairs or Poynting flux?
Part of the answer may come from observations in the X--ray band,
where a feature due to bulk Comptonization is expected
(Begelman \& Sikora 1987;  Sikora, Begelman \& Rees 1994;
Sikora et al. 1997; Moderski et al. 2004).
The spectrum of this feature has been calculated recently
(Celotti, Ghisellini \& Fabian, 2006) and has two components,
both with a quasi--blackbody shape: the first, produced by
the bulk Compton process between the cold jet material and
the photons coming directly from the accretion disk, is predicted to
be at UV energies and it is thus unobservable, while the second,
occurring in the X--ray band, is produced by comptonizing the
broad line photons and it can be observed if
other contributing processes, namely the self--Compton flux,
do not hide it.
The detection of this component would not immediately prove that the
jet is matter dominated (it could be made by initially cold e$^\pm$
pairs), but we could determine both the total number of leptons
and the $\Gamma$--factor.
This information, together with other consideration
about the possible values of the magnetic field, would greatly
help in determining the real matter content.

\subsection{GRBs}

Being transient events, for GRBs it is often more
meaningful to consider the total emitted energy instead
of their luminosities.
This can be calculated from the fluence (time integrated flux)
assuming isotropic emission, and the resulting $E_{\rm iso}$
values range between $10^{49}$ to $10^{54}$ erg (see Fig. \ref{fig1}).
The most energetic events emit approximately (if isotropic)
the equivalent of one solar mass entirely transformed in energy.
Since there will be some efficiency $\eta$ involved, the total
(i.e. kinetic) energy of the jet must be larger by the factor $1/\eta$.
These huge values must be corrected if the emission is collimated.
There are indeed good evidences that this is the case,
because the lightcurve of the afterglow often shows a break
(at 0.2--5 days) which, if achromatic, is what is expected when
the bulk Lorentz factor $\Gamma$ (which is decreasing during the afterglow)
becomes equal to $1/\theta_{\rm j}$, the inverse of the jet opening angle.
Before this time, the observer sees only a fraction of the entire
jet surface (corresponding to an angle $1/\Gamma< \theta_{\rm j})$
while after this time the available emitting area is the entire jet surface.
The jet break time therefore measures $\theta_{\rm j}$ (through some
assumptions about the efficiency and the circumburst density).
The isotropic energetics can then be corrected by the solid
angle $(1-\cos\theta_{\rm j}$), and the resulting
energetics $E_\gamma$ (see Fig. \ref{fig1}) are distributed in
a narrower range than  the isotropic ones, between $3\times 10^{49}$
and $10^{51}$ erg (Frail et al. 2001; Ghirlanda et al. 2004a).
Some independent constraints on the energetics
can be derived by the radio afterglow (see e.g. Frail, Waxman \& Kulkarni 2000
for GRB 970508), which can be thought as
the equivalent of the radio lobes in radiogalaxies, but unfortunately
these estimates could be made only for a few GRBs.

\begin{figure}
\vskip -1 true cm
\begin{tabular}{ll}
\hskip -1 true cm
\includegraphics[width=0.55\textwidth]{sequence.ps}
&
\hskip -0.5 true cm
\includegraphics[width=0.55\textwidth]{ghirlanda.ps}
\end{tabular}
\caption{
{\bf Left:} The blazar sequence.
The average SED of blazars belonging to different
bins of radio luminosity.
As can be seen, the radio luminosity tracks well
the bolometric one.
More powerful blazars are redder, in the sense that the
peak frequencies of both the synchrotron and the
inverse Compton components have smaller values than in
low power blazars. From Fossati et al., (1998) and
Donato et al. (2001).
{\bf Right:} The GRB sequences.
The peak energy of the spectrum of the prompt emission
(time integrated) of GRBs correlates both with the
isotropic bolometric energetics (on the right
of the figure, the so called ``Amati" relation",
$E_{\rm peak}\propto E_{\rm iso}^{1/2}$)
and with the collimation corrected energy (on the
left, the so called ``Ghirlanda" relation, as
derived assuming that the circumburst density
has a wind ($\propto r^{-2}$) profile.
In this case the correlation is linear.
Note that ``bluer" objects are brighter.
Adapted from Ghirlanda et al. 2004 and Nava et al. 2006.}
\label{fig1}
\end{figure}

\section{Families}

\subsection{The blazar sequence}

Fossati et al. (1998), considering $\sim$100 blazars belonging to different
complete (flux limited) samples, constructed the average
SED for blazars in different bins of radio luminosities,
finding a well defined trend, illustrated in Fig. \ref{fig1}.
The SED is always characterized by two broad peaks,
which can be interpreted as due to the synchrotron and the
inverse Compton processes.
The peak frequency of both peaks shifts to smaller values
(the blazar become redder) when increasing the bolometric luminosity.
At the same time, the high energy peak becomes more prominent.
Note also the ``valley" between the two broad peaks, located
in the X--ray band.
This minimum of the SED strongly constrains the location
of the emission region, at least in powerful blazars with
broad emission lines, which reveals the presence of a
(normally) radiating accretion disk.
The argument is as follows: if the emitting region of the jet
is too close to the accretion disk, then the produced $\gamma$ rays
can interact with the X rays produced by the hot corona of the disk,
producing e$^\pm$ pairs.
These pairs are born relativistic, and rapidly cool by inverse Compton
scattering the UV radiation of the disk.
A significant fraction of this radiation would then fall in the ``valley",
contrary to what is observed.
As an order of magnitude, one can derive that the jet dissipation region
must be at least at $\sim100 R_{\rm s}$ (Schwarzschild radii) from the disk.
Note also that this argument apply only if the jet emission region
produces $\gamma$ rays: if for some reason they are not produced
(in a significant quantity), then the argument does not apply, and
we can have dissipation also close to the disk
(Katarzynski \& Ghisellini 2006).

The blazar sequence has been interpreted by Ghisellini et al. (1998)
as due to the different radiative cooling suffered by the
emitting electrons in sources of different powers.
In summary, powerful jets always live in systems where also
the broad emission lines are present and strong.
The line photons, used as seeds for the inverse Compton scattering,
efficiently cool the relativistic electrons, producing a break
(a cooling break) in the electron distribution at low energies.
Low power BL Lacs live in systems of no or weak emission lines:
only electrons of very high energy cool fast, and this produces
a break in the electron distribution at high energies.
The resulting SED is therefore bluer, with less power in the
high energy peak.

\subsection{The GRB sequences}

In very recent years, several correlations have been discovered,
linking the spectral properties of GRBs and their total energetics
or their peak luminosities (for a review see Ghirlanda, Ghisellini \& Firmani 2006).
All of them involve the energy at which the time integrated spectrum
peaks [in a $E F(E)$ diagram].
The fact that GRBs are ``bluer when brighter"
was indeed the surprising discovery after the finding
(by Frail et al. 2001 among others), that the collimation corrected
energetics are more narrowly distributed than the isotropic ones.
In summary, the recently found correlations are
(the exponents are approximate):
\begin{itemize}
\item
$E_{\rm iso}\propto E_{\rm p}^2$: the so called ``Amati" relation
(Amati et al. 2001; Amati 2006);
\item
$L_{\rm p, iso}\propto E_{\rm p}^2$; the so called ``Yonetoku" relation,
(Yonetoku et al. 2004), where $L_{\rm p, iso}$ is the peak luminosity
of the prompt emission;
\item
$E_{\gamma}\propto E_{\rm p}$: the so called ``Ghirlanda" relation,
where $E_\gamma =(1-\cos\theta_{\rm j})E_{\rm iso}$ is the collimation
corrected energetics.
It is linear if the circumburst density has a $r^{-2}$
(wind--like) profile (Nava et al. 2006), and it is of the form
$E_\gamma\propto E_{\rm p}^{3/2}$ if the circumburst medium is
homogeneous (Ghirlanda, Ghisellini \& Lazzati 2004);
\item
$L_{\gamma}\propto E_{\rm p}^{1.2}$: where $L_\gamma$ is
the collimation corrected peak luminosity
(Ghirlanda, Ghisellini \& Firmani 2006);
\item
$E_{\rm iso}\propto E_{\rm p}^2 t_{\rm j}^{-1}$, where $t_{\rm j}$ is the time
when the afterglow light curve shows a steepening (Liang \& Zhang 2005);
\item
$L_{\rm p, iso}\propto E_{\rm p}^{3/2} T_{0.45}^{-1/2}$:
the so--called ``Firmani" relation,
where $T_{0.45}$ is the time of enhanced burst emission, also used
for characterizing the variability behavior (Firmani et al. 2006).
\end{itemize}
In Fig. \ref{fig1} (right panel) we show the Amati relation and,
on the same plot, the Ghirlanda relation in the ``wind case".
The derivation of the jet angle involves some assumptions, such as the
value of $\eta$ (assumed to be equal for all bursts) and the
density profile of the circumburst medium (homogeneous or wind--like).
Note that in the wind--like case, the $E_\gamma$--$E_{\rm p}$ correlation
becomes linear (this is the correlation shown in Fig. \ref{fig1};
from Nava et al. 2006).
A linear correlation means that the number of photons in the peak
is the same for all bursts, and that the correlation is linear also in the
comoving frame of the bursts, since we transform both $E_\gamma$ and
$E_{\rm p}$ by the same factor $\Gamma$.
The equally tight correlation found recently by Firmani et al. (2006),
has the virtue of being model--independent and assumption free.
The tightest correlations (with reduced $\chi^2\sim 1$) are the Ghirlanda,
the Liang \& Zhang and the Firmani ones.
The latter has also the virtue to relate quantities belonging to the
prompt emission only, while the other two need to measure $t_{\rm j}$
by measuring the lightcurve of the afterglow.

There have been a few attempts to explain these correlations,
especially the Amati and the Ghirlanda ones.

Eichler \& Levinson (2005) considered them to be due to a viewing
angle effect, and consider a ``ring" geometry for the emission site.
In their view, all the bursts are intrinsically the same:
they have the same energetics and intrinsic $E_{\rm p}$.

Rees \& Meszaros (2005) pointed out that a tight correlation
between the energetics and the spectral peak suggest black body
emission, since in this case the number of free parameters is minimum.
They then suggest that the Amati relation can be obtained
if there is a correlation between the amount of dissipated power and the
location at which dissipation occurs.
Recently, Thompson (2006) and Thompson, Meszaros \& Rees (2006)
have suggested that the dissipation may instead occur at
more or less the same location, thought as the radius of the
progenitor star, and therefore of the order of $10^{10}$ cm
from the black hole.
The other ingredient is that, at this radius, the bulk Lorentz factor
should be of the order of $1/\theta_{\rm j}$.
This is justified by causality arguments, even if, in Thompson (2006)
the possibility to have a faster on--axis spine is mentioned.
This interpretation of the Amati correlation requires that most
of the observed energetics and the spectral peak corresponds
to a blackbody.
According to this idea, therefore, the observed spectra should
have a relatively strong blackbody component, especially visible
when analyzing time resolved spectra (the temperature may evolve in
time and the blackbody not be well visible in the time integrated spectrum).
Attempts to analyze the spectra this way have been done
(Ghirlanda et al. 20003; Ryde et al., 2005; Bosnjak et al. 2006)
and the decomposition of the time resolved spectrum as a sum
of a power law and a blackbody indeed yield acceptable results
when considering the BATSE energy range only (i.e. between
30 keV and $\sim$1 MeV).
However, the index of the fitted power law is in this case
systematically softer than the low energy index ($\alpha$)
derived fitting a Band function, and the extrapolation to lower energies
becomes in conflict with the data of the Wide Field Camera
of $Beppo$SAX, for those bursts observed both by BATSE and $Beppo$SAX
(Ghirlanda 2006; Bosnjak et al. 2006 in prep.).

\section{Dissipation}

In 1978, Rees proposed that the jet of M87 could be powered
by collisions among different part of the jet itself,
moving at different speeds.
When colliding, they would produce shocks, giving rise to the
non--thermal radiation we see.
Although born in the AGN field, this idea of {\it internal shocks}
grew up more robustly in the GRB field, to become the
``paradigm" to explain their prompt emission
(see e.g. Rees \& M\'esz\'aros 1994).
Faster and later shells can then catch up slower earlier ones,
dissipating part of their bulk kinetic energy into radiation.
However, all shells are and remain relativistic: after the collision
the merged shells move with a bulk Lorentz factor which is
intermediate between the two initial ones.
It is then clear that this mechanism has a limited efficiency, because
only a small fraction of the bulk energy can be converted into radiation
(unless the contrast between the two initial Lorentz factor is huge,
see Beloborodov 2000; Takagi \& Kobayashi 2005).

This weak point for GRBs becomes instead a point in favor for
blazars: in fact, for them, we {\it require} a small efficiency,
for letting most of the bulk energy go un--dissipated to the
outer radio--lobes.
The other main point in favor for internal shocks in blazars is that
it naturally explains why the main dissipation event, in the
relativistic jets, should occur at a few hundreds $R_{\rm s}$:
if two shells are initially separated by $R_0 \sim$ a few $R_{\rm s}$
and are moving --say-- with $\Gamma=10$ and $\Gamma=20$, they will
collide at $R_{coll}\sim R_0 \Gamma^2$, i.e. exactly where it is required.

A detailed study via numerical simulations of the predictions of
this model in blazars has been carried out by Spada et al. (2001)
for powerful objects (like 3C 279) and by Guetta et al. (2004)
for less powerful objects, like Mkn 421 and BL Lac itself,
with encouraging results.
More recently, however, new ideas came out for the TeV BL Lac case,
involving jet deceleration (Georganopoulos \& Kazanas 2003) or
interactions between a fast jet spine and a slower layer (Ghisellini,
Tavecchio \& Chiaberge 2005).

As discussed above, $R_{coll}\sim 10^2 R_{\rm s}$ when the typical
$\Gamma$--factor of the shells is around 10.
What happens if the source produces, for some time, shells with
a smaller $\Gamma$?
This possibility has been explored by
Katarzynski \& Ghisellini (2006), who assumed a jet working
at a constant average power, with shells of equal kinetic
energy $\Gamma Mc^2$ ($M$ is the mass of the shell).
Collisions of shells with low $\Gamma$ (and therefore large $M$)
occur closer to the jet apex than collisions between shells with
larger $\Gamma$.
This means that they dissipate when they are denser and more compact:
the magnetic field is larger, and the external radiation produced
by e.g. the broad line region, as seen
in the comoving frame, is less boosted.
All these facts means that the energy is radiated more through the
synchrotron process rather than by the Inverse Compton mechanism.
In turn, this means that fewer $\gamma$ rays are produced, and
even if they are partly transformed into pairs, they might not
``fill the valley" between the observed two broad peaks.
In this case we should observe a large (even dramatic) variability
in the IR--UV synchrotron and in the
X--ray self--Compton part of the spectrum, and much less
in the MeV--GeV band produced by scattering broad line photons.
This is an ``economic" way to explain, for instance, the huge
flare observed in the blazar 3C 454.3, in which the optical and the
X--ray flux increased by a factor 100 with respect to its ``normal"
level.
AGILE and GLAST will easily confirm or falsify these ideas.

\section{The maximum jet power}

If relativistic jets are powered (at least initially)
by a Poynting flux, we can derive a simple expression for
the maximum value of their power, under some reasonable assumptions.
The Blandford \& Znajek (1977) power can be written as:
\begin{equation}
L_{\rm BZ}\, \sim \, 6\times 10^{20}
\left({ a\over m}\right)^2 \left( {M_{\rm BH}\over M_\odot} \right)^2 B^2
\,\,\,\, {\rm erg~s^{-1}}
\end{equation}
where $(a/m)$ is the specific black
hole angular momentum ($\sim 1$ for maximally rotating black holes),
and the magnetic field $B$ is in Gauss.
Assume that the value of magnetic energy density
$U_B\equiv B^2/(8\pi)$ close to the hole is a fraction $\epsilon_B$
of the available gravitational energy:
\begin{equation}
U_B \, = \, \epsilon_B {GM_{\rm BH} \rho \over R} \,=\,
\epsilon_B {R_s \over R} {\rho c^2 \over 2}
\end{equation}
The density $\rho$ is linked to the accretion rate $\dot M$ through
\begin{equation}
\dot M  \, = \, 2\pi RH \rho \beta_{\rm R}c
\end{equation}
where $\beta_{\rm R} c$ is the radial infalling velocity.
The mass accretion rate $\dot M$ is linked to the observed luminosity
produced by the disk
\begin{equation}
L_{\rm disk} \, =\, \eta \dot M c^2
\end{equation}
The BZ jet power can then be written as:
\begin{equation}
L_{\rm BZ} \, \sim\, \left({ a\over m}\right)^2 { R_{\rm s}^3 \over R^2 H}
{\epsilon_B  \over \eta}
{L_{\rm disk} \over \beta_{\rm R} }
\end{equation}
The maximum jet power is obtained setting $R\sim H\sim R_{\rm s}$ and
$a/m$, $\epsilon_B$, $\beta_{\rm R}$ equal to unity.
In this case
\begin{equation}
L_{\rm BZ, max} \, \sim\, {L_{\rm disk} \over \eta } \, =\, \dot M c^2
\end{equation}
This is in qualitative agreement with what can be estimated in blazars
and microquasars, and also in GRBs.
High power blazars (and therefore FR II radiogalaxies) can have
jets with a similar power than what emitted by the accretion disks,
especially during (non--thermal) flares, when $P_{\rm jet}$ can even
be slightly larger.
For low power BL Lacs (and therefore FR I radiogalaxies) the luminosity
emitted by the accretion disk is largely sub-Eddington, and the efficiency
$\eta$ can be much smaller than the canonical 0.1 value.
This means that these systems can be largely jet dominated.
Of course the most extreme objects in this respect are GRBs,
in which most of the energy originated in the accretion process
is advected into the hole, due to the huge densities and optical depths.
According to the simple estimate given above, the most efficient
``engines" in the universe are jets, not accretion disks.

\section{Cosmological implications}

Hard X--ray photons are the ones suffering less from any absorption
process in their travel from the sources to Earth, while higher
energy photons can interact with IR--optical and UV photons and
be absorbed through the $\gamma$--$\gamma \to$ e$^\pm$ process.
The fact that blazars and GRBs are strong emitters at both hard X--ray
and $\gamma$--ray frequencies allow us to use them as torch-lights to
study the amount of IR--UV cosmic backgrounds.
Furthermore, we have seen that the GRBs sequences (the Ghirlanda
and the Firmani relations) can be used to standardize their energetics,
enabling us to use them as standard candles.

\begin{figure}
\vskip -0.5 true cm
\begin{tabular}{ll}
\hskip -0.7 true cm
\includegraphics[width=0.55\textwidth]{1101.ps}
&
\includegraphics[width=0.45\textwidth]{0906.ps}
\end{tabular}
\caption{
{\bf Left:}
The TeV spectrum of 1ES 1101--232 (red points)
de--absorbed assuming different IR backgrounds
(blue points), decreasing from left to right.
The de--absorbed spectrum is very hard for the
left and central panel spectra.
The minimum possible IR background value
corresponds to a an intrinsic photon spectrum
$\propto \nu^{-1.5}$ which is
considered by Aharonian et al. (2006) has
the most likely.
This conclusion was however questioned by
Katarzynski et al. (2006), who pointed out
that the limiting slope for the inverse Compton
photon spectrum, when the electron distribution has
a low energy cut--off, is $\propto \nu^{-2/3}$,
making stronger IR backgrounds still possible.
(Figure from Aharonian et al. 2006.)
{\bf Right:}
The SED of Q0906+6930. The solid line is a model.
Also shown is the limiting sensitivity of GLAST (red line),
to show that GLAST can well determine the $\gamma$--ray
spectrum  and any possible cut--off due to photon--photon
collisions. From Romani (2006).
}
\label{fig2}
\end{figure}

\subsection{Blazars}

Cherenkov telescopes have detected a dozen blazars in the TeV band,
and EGRET detected $\sim$60 blazars in the GeV band.
Blazars, as a class, are as strong $\gamma$--ray emitters.
For the TeV emitting BL Lacs one of the key issues debated now is
the determination of the poorly constrained cosmic IR background:
we can use the absorption of their TeV photons as a tool to determine
the level of the IR backgrounds.
This however requires to assume an intrinsic spectrum,
which is model--dependent.
On the other hand, when we will have a sufficient number of sources
distributed in redshift, it will be possible to determine the
absorption cut--off energy in more model--independent way.
In Fig. \ref{fig2} we report the observed and the de--absorbed
TeV spectrum of the blazar 1ES 1101--232 (Aharonian et al. 2006).
Different IR backgrounds yield very different de--absorbed spectrum
and Aharonian et al. (2006) argued that to have a ``reasonable"
intrinsic spectrum (the one with the photon spectral index $\Gamma=1.5$)
it is necessary to assume the minimum possible IR background.
This conclusion was however questioned by
Katarzynski et al. (2006), who pointed out
that the hardest possible intrinsic spectrum
is the one having $\Gamma=-2/3$, corresponding to an
underlying electron distribution with a low energy cut--off.
In this case less extreme IR backgrounds are still possible.

In the AGILE and especially GLAST era the same
study done now for the IR cosmic backgrounds can be done
also for the cosmic optical and UV backgrounds,
which are responsible for the absorption of GeV
photons (see e.g. Chen, Reyes \& Ritz 2004).
This means also to detect sources at relatively large
redshifts ($z>2$, because the Universe starts to
be opaque to tens of GeV photons for these distances),
and to be able to detect a cut--off in their high
energy spectrum.
As a ``preview" of what GLAST can do for
cosmology using blazars, we show in Fig. \ref{fig2}
the SED of  Q0906+6930, the most distant blazars.
This source was recently associated by Romani et al. (2004)
to a previously unidentified EGRET source,
and it lies at the redshift $z=5.48$.
GLAST can likely find hundreds of blazars at redshifts
large enough to study the intervening absorption of
their GeV photons, once we can disentangle the absorption
made by the locally (i.e. accretion disk and broad line region)
produced optical--UV photons (Tavecchio et al. in preparation).

\subsection{GRBs}

The Ghirlanda, the Laing \& Zhang and the Firmani relations are
tight enough to enable their use to standardize the energetics
of GRBs, and therefore we can use them as standard candles.
The obvious advantage is that GRBs are brighter than SN Ia,
therefore we can explore and measure the Universe
beyond $z\sim 1.7$, the current redshift limit of SN Ia.
Furthermore, GRBs are free from extinction uncertainties
which may affect SN Ia.
There is however a problem: the correlations are found
using a particular cosmological problem, and we pretend
to use them to find the cosmological parameters.
This {\it circularity problem} has been dealt with in
several ways.

\begin{figure}
\vskip -1.5 true cm
\begin{tabular}{ll}
\hskip -1.5 true cm
\includegraphics[width=0.7\textwidth]{firmani_hd.ps}
&
\hskip -2.5 true cm \includegraphics[width=0.725\textwidth]{firmani_oo.ps}
\end{tabular}
\vskip -1 true cm
\caption{
{\bf Left:} The Hubble diagram including SN Ia of the Legacy Survey
(red points)
and the GRBs analyzed in Firmani et al. (2006, blue points).
{\bf Right:}  The constraints on $\Omega_{\rm M}$--$\Omega_{\Lambda}$
set by considering SN Ia alone (green contours) and SN and GRBs
together (red shaded area).
It can be appreciated how much GRBs help in constraining
the cosmological parameters, and also the consistency with the
so--called "concordance cosmological model" (black cross).
(Both figures from Firmani et al. 2006b)
}
\label{fig3}
\end{figure}

\begin{enumerate}

\item
If we had a robust theoretical interpretation of the
(Ghirlanda, Firmani or Laing \& Zhang)
correlation, then we could trust the theoretically derived slope and
normalization to use it for any choice of $\Omega_{\rm M}$
and $\Omega_{\Lambda}$
(see Dai et al. 2004 for what happens assuming
a fixed Ghirlanda correlation of the form $E_\gamma \propto E_{\rm p}^{3/2}$).

\item
Ghirlanda et al. (2004b) proposed to use the scatter around the
fitting line of the $E_{\rm peak}$--$E_\gamma$ correlation as
an indication of the best cosmology.
In other words, they find a given Ghirlanda correlation using
a pair of \mbox{$\Omega_{\rm M}$}, \mbox{$\Omega_{\Lambda}$}, and measure the scatter of the points
around this correlation through a $\chi^2$ statistics.
Then they change \mbox{$\Omega_{\rm M}$}, \mbox{$\Omega_{\Lambda}$}, find another correlation,
another scatter and another $\chi^2$ which can be compared
with the previous one.
Iterating for a grid of \mbox{$\Omega_{\rm M}$}, \mbox{$\Omega_{\Lambda}$} ~sets, they can assign to any point
in the \mbox{$\Omega_{\rm M}$}, \mbox{$\Omega_{\Lambda}$} ~plane a value of $\chi^2$, and therefore draw
the contours level of probability.

\item
Firmani et al. (2005) proposed a more advanced method,
based on a Bayesian approach, which is currently the best
to cure the circularity problem.
The basic idea is to take into account the information that
the correlation is unique, even if we do not know its slope and normalization
(see Firmani et al. 2005 and Ghisellini et al. 2005 for the description
of the method).

\item
If one had many GRBs at low redshifts (i.e. $z<0.1$) the slopes
and normalizations of the correlations would be largely
cosmology--independent, but there are not many nearby bursts.
On the other hand, the same goal can be achieved
using bursts in a narrow redshift bin:
as suggested by Ghirlanda et al. (2006) and Liang \& Zhang (2006)
a dozen of objects in a $\Delta z \sim 0.2$ around $z=1$ are sufficient.

\end{enumerate}

The other important point to be stressed is that GRBs should be considered
as {\it complementary}, not {\it alternative} to SN Ia.
Broadly speaking, one can think to the two classes of objects as standard candles
at large, located at low (SN Ia) and large (GRBs) redshifts.
The fact that also SN Ia need to be ``standardized" (through the
``stretching" relation: more powerful SN Ia have lightcurves that decay
more slowly, see e.g. Phyllips et al. 1993) even requires a very similar
treatment for SN Ia and GRBs (even more so for the SN Ia of the
Supernova Legacy Survey; Astier et al. 2006).

Applying the Bayesian approach to the Firmani relation, Firmani et al. (2006b)
already derived a very significant reduction in the uncertainties
in the \mbox{$\Omega_{\rm M}$}, \mbox{$\Omega_{\Lambda}$}~ plane, when considering GRBs and SN Ia together.
This can be seen in Fig. \ref{fig3}: on the left there is the
classical Hubble diagram for SN Ia and GRBs, and on the right the
\mbox{$\Omega_{\rm M}$}, \mbox{$\Omega_{\Lambda}$}~ constraints when considering GRB only (large elliptical
contours) and GRB+SN (smaller shaded area).
One can appreciate how the addition of GRBs can improve the
constraints, and also the fact that the joint fit is perfectly
consistent with the so--called ``concordance cosmological model"
(i.e. \mbox{$\Omega_{\rm M}$}$\sim0.3$ and \mbox{$\Omega_{\Lambda}$}$\sim 0.7$).

\section{Conclusions}

The values of the bulk Lorentz factor $\Gamma$ estimated for blazars,
have been increased during the years (except the $\Gamma$=100 peak reached
when intraday variability was thought to be completely intrinsic).
Values of 40--50 are commonly proposed in the literature,
at least to explain $\gamma$--ray flares and especially for low power,
TeV emitting, BL Lacs, and are also occasionally derived through
the superluminal speeds of high frequency VLBI knots.
On the other hand the ``old" limits in GRBs (i.e. $\Gamma>100$)
coming from the compactness arguments rely on associating
the high energy $\gamma$--ray emission to the prompt.
This is to be confirmed, because, instead, the high energy emission could
belong to the afterglow.
For this issue AGILE and GLAST will be crucial.

Dissipation of the kinetic power of the jets into radiation
could happen at a variable distance from the black hole.
Internal shocks are natural candidates as a mechanism able
to transport bulk kinetic energy and to dissipate at the right place,
namely at 100--300 $R_{\rm s}$ in powerful blazars.

The sequences of blazars and GRBs are opposite: redder when brighter
for blazars, bluer when brighter for GRBs.
For blazars, the sequence has been explained by the increasing
dominance of the radiative cooling over the particle "heating"
as the luminosity increases.
For GRBs, it is less clear.

The power of these systems can exceed the luminosity of the
accretion disk, and be of the order of $L_{\rm jet} \sim \dot M c^2$.

There are good prospects for the use of blazars and GRBs for cosmology.
For blazars, they are linked to the existing TeV Cherenkov observations
and the coming GeV--band AGILE and GLAST satellites.
Indeed, the use of distant blazars to measure the optical and UV
backgrounds can turn out to be one of the key programs for GLAST.
Also for GRBs GLAST can be crucial, due to its BATSE--like
instruments, able to measure the peak energy in a reliable way
for a sufficient number of GRBs, whose location could be obtained
by a simultaneous detection by SWIFT or by ground based
robotic telescope images.

\vskip 0.3 true cm
{\bf Acknowledgements  ---} I thank Annalisa Celotti, Giancarlo Ghirlanda
and Fabrizio Tavecchio for years of fruitful collaboration.

Blazars and Gamma Ray Bursts


Spatially inhomogeneous discrete states of a superconductor upon injection of nonequilibrium quasiparticles from a point contact with a normal metal

Abstract: The current-voltage characteristics (IVC) of $S-c-N$ point contacts of superconductors with a small coherence length ${{\xi}_{0}}$ reveal steps with discrete values of the differential resistance. This peculiarity is associated with a transition of the contact region of the superconductor to a spatially inhomogeneous state under the influence of the current injection of nonequilibrium quasiparticles penetrating the superconductor to a depth ${{l}_{E}}$ . The role of the relaxation of the disbalance between the occupancies of electron- and hole-like branches of the quasiparticle spectrum is manifested in the displacement of the position of the singularities on the IVC towards higher energies upon an increase in the magnetic field and/or temperature. This effect was observed in superconductors with different ratios of the contact diameter $d$ and ${{\xi}_{0}}$ or ${{l}_{E}}$ in the series $Ta\ (d\ll {{\xi}_{0}},\ {{l}_{E}})\to NbS{{e}_{2}}$, $N{{b}_{3}}Sn\ (d\gtrsim {{\xi}_{0}},\ {{l}_{E}})\to YB{{a}_{1.25}}S{{r}_{0.75}}C{{u}_{3}}{{O}_{7-\delta}}(d\gg {{\xi}_{0}},\ {{l}_{E}})$. Apparently, a jumplike displacement of the boundary between regions with suppressed and equilibrium values of the energy gap near the contact is responsible for the oscillations observed on the IVC of point contacts between single crystals of high- temperature superconductors (HTS) and a normal metal. The resistance-periodic step-like structure of the IVC allows us to estimate the penetration depth ${{l}_{E}}$ of the electric field in $NbS{{e}_{2}}$ and $YBaSrCuO$.

\section{INTRODUCTION}
The state of a superconductor in the contact region of $S-c-N$ or $S-c-S$ structures ($c$ stands for a constriction or microbridge) in the current regime has been attracting the attention of researchers for a long time. The simplest model \cite{Iwanyshyn} taking into account the violation of superconducting state by a critical current density or Joule heating assumes a local equilibrium between quasiparticle excitations, the condensate, and phonons. A similar model was worked out by the authors of Refs. \cite{Hahn1,Hahn2,Hahn3} for explaining the $dV/dI\left( V \right)$  oscillations of the $Ta-Ag$  point contacts. However, the investigations of metals by the point-contact spectroscopy methods showed that in most cases there is no local equilibrium between electrons and phonons in the current concentration region \cite{Yanson1}. Moreover, there is no equilibrium in a superconductor between quasiparticle excitations and the condensate, which is manifested in the disbalance between the population densities of the electron- and hole-like branches in the spectrum of quasiparticle excitations. It has been shown by Yanson et al. \cite{Yanson2,Yanson3} that the jumplike decrease in the excess current with increasing bias voltage $V$  for $S-c-N$ point contacts $Ta-Ta$, $Ta-Cu$, and $Ta-Au$ in which the spectral regime of current passage is observed in the entire interval of applied bias voltages, i.e., the inequality $d<{{l}_{\varepsilon }}$ , ${{\left( {{l}_{i}}\ {{l}_{\varepsilon }} \right)}^{{1}/{2}\;}}$  ($d$  is the contact diameter and ${{l}_{\varepsilon }}$ , ${{l}_{i}}$  are the inelastic and elastic electron mean free paths) associated with the attainment of the critical concentration by the nonequilibrium quasiparticles is observed and is essentially of nonthermal type. In highly pure $Ta-Cu$ contacts of size $d\ll {{\xi }_{0}}$  (${{\xi }_{0}}$  is the coherence length of the superconductor), this decrease is of the order of several percent and is accompanied by an increase in the differential resistance of the same order of magnitude.

In the present work, a steplike structure of the IVC is observed in $S-c-N$ type point contacts between the single crystals of $NbS{{e}_{2}}$, $YB{{a}_{1.25}}S{{r}_{0.75}}C{{u}_{3}}{{O}_{7-\delta }}$ and noble metals $Cu$ , $Ag$ . This structure is reminiscent of the IVC structure associated with the formation of the phase slip centers in thin superconducting filaments. Such a structure appears when ${{l}_{E}}$ is of the order of, or smaller than, $d$ . It is shown that the formation of this structure has a nonthermal origin and is associated with an essentially nonequilibrium state of the superconductor in the contact region. In the same way as the emergence of a nonequilibrium singularity on the IVC of tantalum point-contacts upon a decrease in temperature is associated usually with the characteristic energies $h{{\omega }_{ph}}$  of the phonon spectrum \cite{Yanson3}, the jumps in the excess current and differential resistance on the IVC appearing upon a temperature decrease of superconducting point contacts with a small ${{\xi }_{0}}$  are also fixed near certain characteristic energies $e{{V}_{ph}}$  corresponding to the characteristic energies of the phonon spectrum \cite{Bobrov} or other excitations with which the conduction electrons effectively interact \cite{Yanson4}. On the other hand, the multiplicity of the jumps in the differential resistance may lead to oscillating $dV/dI\left( V \right)$  dependences which are frequently observed in experiments \cite{Kirtley}.
\subsection{DISCUSSION OF EXPERIMENTAL RESULTS}
\textbf{1.} $\textit{\textbf{NbS}}{{\textit{\textbf{e}}}_{\textit{\textbf{2}}}}$. Figure \ref{Fig1} shows the IVC of an $S-c-N$ point contact between a freshly cleaved single crystal of $NbS{{e}_{2}}$ and a copper needle oriented at right angles to the basal plane (see the upper inset to Fig.\ref{Fig1}).
\begin{figure}[]
\includegraphics[width=8cm,angle=0]{Fig1.pdf}
\caption[]{Current-voltage characteristics of the $NbSe_2-Cu$ junction having a resistance $170\ \Omega $  at $T=1.68\ K$ . The dot-and-dash line shows the predicted IVC in the normal state, the dashed line corresponds to discrete values of differential resistance at the steps 60, 117, 150, and $170\ \Omega $. The upper inset shows the geometry of the experiment, and the lower inset shows the initial segement of the IVC (the dashed lines are the tangents to various segments of the IVC.}
\label{Fig1}
\end{figure}
\begin{figure}[]
\includegraphics[width=8cm,angle=0]{Fig2.pdf}
\caption[]{The $S-c-N$ model of a point contact in the form of a circular hole of diameter $d$ in an impenetrable partition ($d>{{l}_{E}}>{{\xi }_{0}}$). The dashed curves are the boundaries between different nonequilibrium state regions of the superconductor, which are separated from the contact center by integral multiples of ${{l}_{E}}$.}
\label{Fig2}
\end{figure}
 For $eV>{{\Delta }_{\bot }}$ (${{\Delta }_{\bot }}\lesssim 0.4-0.6\ meV$   is the energy gap along the $c$-axis), the electric field penetrates the superconductor to a depth ${{l}_{E}}\approx {{l}_{\varepsilon }}$  (${{l}_{\varepsilon }}$  is the mean inelastic relaxation time for quasiparticles with energies $0<\varepsilon <eV$), and the differential resistance is apparently determined by the resistance of the ellipsoid with $z\approx {{l}_{E}}\lesssim d$ (Fig.\ref{Fig2}). This resistance is much higher than the zero-bias resistance (${{R}_{0}}=30\ \Omega $) (see the lower inset to Fig.\ref{Fig1}) associated with the normal face. In the model considered by us, it is significant that the contact diameter is small in comparison with the inelastic relaxation length of quasiparticles in the normal faces ($N$) for all bias voltages $eV$. For noble metals, this is obviously valid for $d\lesssim {{10}^{-5\ }}cm$. Hence, electrons (holes) are incident from the normal metal on the $N-S$  interface with an excess energy relative to the electrochemical pair potential in the superconductor (right up to $eV$).

Subsequent regions of IVC in Fig.\ref{Fig1} are reminiscent of peculiarities caused due to the formation of phase slip centers in a thin superconducting filament. The differential resistance ${{R}_{D}}$  changes abruptly on quasilinear segments of IVC: ${{R}_{1}}=62$, ${{R}_{2}}=117$, ${{R}_{3}}=150$, ${{R}_{4}}=174\ \Omega $. A consistent decrease in the value of the increment $\Delta {{R}_{D}}$ in successive steps is in qualitative agreement with the model with a jumpwise displacement of the ${S}'S$  boundary into the bulk of the superconductor (see Fig.\ref{Fig2}). Here, ${S}'$ is the superconducting region adjoining the constriction with a large nonequilibrium concentration of quasiparticle excitations into which the electric field penetrates, and $S$  is an equilibrium superconductor- corresponding to the state away from the constriction. In this model, we assume that the penetration depth ${{\lambda }_{L}}$  magnetic field into the superconductor is larger than $d$. It is obvious that even for voltages that exceed the region where the steps are observed, the excess current on the IVC continues to be finite and nearly constant. This means that the order parameter (energy gap) in the contact region ${S}'$ also has a finite value whose order of magnitude is equal to its unperturbed value.
\begin{figure}[]
\includegraphics[width=8cm,angle=0]{Fig3.pdf}
\caption[]{Stepwise variations in the IVC of an $NbSe_2-Cu$  contact at temperatures $T=1.85$  (1) and $3.22\ K$  (2). The dot-and-dash line is parallel to IVC for large values of $V$ $\left( {{R}_{N}}=60\ \Omega  \right)$.}
\label{Fig3}
\end{figure}

Figure \ref{Fig3} shows the IVC of another junction, viz., $NbSe_2-Cu$, at different temperatures. As the temperature \textbf{decreases}, the "critical" current ${{I}_{c1}}$  and power ${{P}_{c1}}={{I}_{c1}}{{V}_{c1}}$ corresponding to a transition from the first to the second step decrease instead of increasing, as would be the case if the IVC jumps were due to the trivial destruction of superconductivity as a result of heating or the attainment of the critical current density along certain current paths. The observed temperature dependence of the IVC points towards a significant role of nonequilibrium quasiparticles. Among other things, the concentration of such particles is determined by the relaxation time which increases with decreasing temperature. Consequently, the "critical" concentration of nonequilibrium quasiparticles, which is one of the reasons behind the phase lamination in a nonequilibrium superconductor at low temperatures, is attained at a lower injection rate (lower ${{I}_{c1}}$ and ${{P}_{c1}}$).

It is interesting to note that the normal increment (${{R}_{2}}-{{R}_{1}}$) in the differential resistance at the quasilinear region corresponding to the second step is double the increment (${{R}_{1}}-{{R}_{0}}$). Since the resistance ${{R}_{n}}$  for small values of $2n{{l}_{E}}/d$  is proportional to $n$  (see \eqref{eq__1}), it can be surmised that the resulting structure will be identical to that of the phase slip centers. Indeed, for $eV\gg \Delta $, the main part of the current in the superconductor near the $N-S$  boundary is transported by the quasiparticle excitations which recombine to form pairs at a distance ${{l}_{E}}$ from the constriction. If the current density $j(\mathbf{r})$  at this distance is less than ${{j}_{c}}$, such a state is stable. Otherwise, a narrow layer (phase slip plane or PSP) is formed at a distance $2{{l}_{E}}$  from the constriction. This layer has a suppressed order parameter $\Delta $  and generates new quasiparticles capable of transporting the "supercritical" current further into the bulk of the superconductor to a depth $3{{l}_{E}}$. Such a process is repeated with increasing current until the spatial region covered by it exceeds the current spread region of the order $d$. It is obvious that the sharply steplike nature of the IVC may be observed only under the condition ${{l}_{E}}<d$, and the number of observed jumps is of the order $d/2{{l}_{E}}$. The contribution of the superconducting face to the contact resistance is given by
\begin{equation} \label{eq__1}
{{R}_{n}}=\left( \rho /\pi d \right)\operatorname{tanh}^{-1}2n{{l}_{E}}/d=2\rho n{{l}_{E}}/{{d}^{2}},\ \ \ \left( n{{l}_{E}}\ll d \right)
\end{equation}
($n$  is an odd integer).
It follows hence that the resistance of the entire contact for large values of $V$ and at the first step are given, respectively, by
\begin{equation} \label{eq__2}
{{R}_{\infty }}=\rho /2d+{{R}_{0}};\ \ {{R}_{1}}=2\rho {{l}_{E}}/{{d}^{2}}+{{R}_{0}}
\end{equation}
while the ratio $d/{{l}_{E}}$  is defined as
\begin{equation} \label{eq__3}
d/{{l}_{E}}=\frac{4}{\pi }\left( \frac{{{R}_{\infty }}-{{R}_{0}}}{{{R}_{1}}-{{R}_{0}}} \right)
\end{equation}
The contact diameter can be estimated from ${{R}_{\infty }}$ \eqref{eq__2}. Unfortunately, the exact value of the resistivity $\rho$ in the region adjoining the contact is now known. It may be much larger than ${{\rho }_{\bot }}\sim {{10}^{-3}}\ \Omega \cdot cm$ measured for a bulk sample along the $c$-axis. Of the two different $NbSe_2-Cu$ contacts whose characteristics are presented in Figs. \ref{Fig1} and \ref{Fig3}, we choose the latter for numerical estimates since it has a lower resistance, and hence a value of $\rho$  closer to the bulk value, for nearly the same number of steps on the IVC (and for the same value of the ratio $d/{{l}_{E}}$). Substituting ${{R}_{\infty }}=62\ \Omega $  and ${{R}_{0}}=10\ \Omega $ into \eqref{eq__2}, we obtain $d\gtrsim 960\ \mathrm{\AA}$, ${{l}_{E}}=145\ \mathrm{\AA}$. The coherence length in $NbS{{e}_{2}}$ is quite small: ${{\xi }_{\parallel }}=77\ \mathrm{\AA}$, ${{\xi }_{\bot }}=23\ \mathrm{\AA}$. Hence the formation of narrow PSP layers with a suppressed gap (Fig. \ref{Fig2} is quite possible.

\textbf{2.} $\textit{\textbf{YBa}}_{\textit{\textbf{1.25}}}\textit{\textbf{Sr}}_{\textit{\textbf{0.75}}}\textit{\textbf{Cu}}_{\textit{\textbf{3}}}\textit{\textbf{O}}_{7-\delta}$. Obviously, the number of observed steps must increase as the inequality $d>{{l}_{E}}>{{\xi }_{0}}$, becomes more stringent. However, in this case $d$ must not exceed the inelastic mean free path of phonons, since the thermal regime will set in otherwise. In this connection, it is interesting to use high-temperature superconductors which have record-low values of ${{\xi }_{0}}$  and $l_{\varepsilon}$. Figure \ref{Fig4}a
\begin{figure}[]
\includegraphics[width=8cm,angle=0]{Fig4.pdf}
\caption[]{(a)\ -\ stepwise structure of IVC for $YBaSrCuO-Ag$ point contacts at $T=4.2\ K$  and ${{T}_{N}}=336\ \Omega $:
1 - $I(eV)$, \\2 - ${{V}_{1}}(eV)\sim dV/dI$;\\
(b),(c)\ -\ IVC and their first derivatives ${{V}_{1}}(eV)$ $\left( {{R}_{N}}=83\ \Omega  \right)$(b)   and $68\ \Omega $(c)).}
\label{Fig4}
\end{figure}
shows the IVC of a point contact formed by $YB{{a}_{1.25}}S{{r}_{0.75}}C{{u}_{3}}{{O}_{7-\delta }}$ and silver as a result of first touches between the electrodes in liquid helium. The contact geometry is shown in the inset. It can be seen that the IVC comprises a large number of (about twenty) steps with a gradually increasing slope $dV/dI$. For subsequent touches between the electrodes, the stepwise structure of the IVC is indistinguishable but its existence can be clearly seen in the first derivative plots (Fig. \ref{Fig4}b). The number of steps decreases in comparison with the preceding contact and their arrangement reveals a regularity manifested in a smooth decrease in the distance between spikes in ${{V}_{1}}(eV)\sim dV/dI$ with increasing step number. The envelope ${{V}_{1}}(eV)$, shown by a dashed line, has minima at $\pm 42\ meV$  which can be interpreted as manifestations of the energy gap $\Delta$ due to a small "tunneling" component of the current emerging as a result of the reflection of a part of the electrons at the physical boundary between the superconductor and the normal metal.
The nonlinearity of IVC for $V=0$  vanishes at $T=70\ K$  (this temperature corresponds to the end of the resistive transition of the bulk sample into the superconducting state). Taking this temperature as $T_c$, we obtain $2\Delta /k{{T}_{c}}\simeq 14$, which is anomalously large (as compared to the value given by the BCS theory) and is nearly equal to the frequently encountered values in the tunnel- or point-contact studies of high-temperature superconductors.

In the same series of experiments, subsequent mutual displacements of electrodes lead to the disappearance of the above-mentioned structure of $dV/dI$  caused, in our opinion, by the discrete nature of the displacement of the ${S}'S$  boundary, and also to a considerable decrease in the gap value $\Delta $  to about $25\ meV$, which points towards a partial degradation of superconducting properties near the contact (Fig. \ref{Fig4}c).

Thus, it can be concluded that the regular stepwise nature of IVC of point contacts involving superconductors with short ${{\xi }_{0}}$  and ${{l}_{E}}$ is observed only for an unperturbed superconductor with a perfect lattice. The disorder, which is created during the formation of the point contact and is manifested, among other things, in a reduced (or averaged over different directions) gap, prevents the formation of a regular system of phase boundaries $S'S$  in the contact region. This explains the absence of a regular structure of $dV/dI(V)$  in the studies involving HTS ceramic materials.

Figure \ref{Fig5}
\begin{figure}[]
\includegraphics[width=8cm,angle=0]{Fig5.pdf}
\caption[]{Periodic structure of the $dV/dI$ characteristic of a $YBaSrCuO-Ag$ contact (with $\left( {R_N}\simeq105\ \Omega  \right)$ at $T=10$ (1) and $20\ K$ (2).}
\label{Fig5}
\end{figure}
shows the ${{V}_{1}}\left( eV \right)$ structure for another contact formed by the same single crystal, recorded at two different temperatures. As the temperature increases, a small displacement of some peaks towards higher voltages can be observed. This displacement is analogous to the one discussed above for $NbSe_2$ (Fig. \ref{Fig3}). Such an effect is observed more clearly upon the application of an external magnetic field (Fig. \ref{Fig6}).
\begin{figure}[]
\includegraphics[width=8cm,angle=0]{Fig6.pdf}
\caption[]{Effect of magnetic field on the peak-like peculiarities of the IVC derivative of a $YBaSrCuO-Ag$ contact at $T=4.2\ K$, $R_N=30\ \Omega$: $H=0$  (1) and $H=40\ kOe$ (2).}
\label{Fig6}
\end{figure}
The characteristics shown in Fig. \ref{Fig6} correspond to a polycrystalline sample obtained by cooling the melt. This may be the reason behind the irregular structure of ${{V}_{1}}\left( eV \right)$. It can be seen that the magnetic field considerably decreases the intensity of peaks. The ${{V}_{1}}\left( eV \right)$ peaks are broadened and displaced towards \textbf{higher} energies.

Hence, it can be concluded on the basis of the temperature and field dependences presented above that the sharp peaks in the ${{V}_{1}}\left( eV \right)$ dependence are not associated with the usual mechanism of the violation of current paths in inhomogeneous HTS samples due to heating or a local increase in the critical current density. Using formulas \eqref{eq__1}-\eqref{eq__3}, we can estimate the order of magnitude of the effective diffusion depth ${{l}_{E}}$ quasiparticle excitations in an HTS. Using the value ${{\rho }_{\parallel }}=200\cdot {{10}^{-6}}\ \Omega \cdot cm$ for ${{R}_{\infty }}=100\ \Omega $, we obtain the current diameter $d\approx {{10}^{-6}}\ cm$. Proceeding from Fig. \ref{Fig4}b as a typical case, the number of layers of thickness ${{l}_{E}}$ in the current concentration region is estimated at about 10. Consequently, ${{l}_{E}}\approx d/10\sim {{10}^{-7}}cm$, which is smaller than the corresponding value for $NbSe_2$ by an order of magnitude.

Let us clarify the meaning of the effective diffusion depth ${{l}_{E}}$ of quasiparticles. Quasiparticles with an excess energy relative to the electrochemical pair potential (right up to $eV$) are injected, into the HTS layer adjoining the $N-S$ boundary. However, if $eV\gg \Delta $, the quasiparticles rapidly relax to states above the energy gap, and hence the value of  $l_E$ is mainly determined by the quasiparticles having an excitation energy of the order of $\Delta $. For such quasiparticles, it is reasonable to use the values $\rho (T)$ at $T=T_c$, as was indeed done by us in this paper.

For the sake of comparison, it is worthwhile to consider the results of point-contact investigations in $Ta$, \cite{Yanson3} whose EPI function is well known. As a rule, the $S-c-N$ contacts with noble metals reveal only one nonequilibrium singularity. It was mentioned above the this singularity has a rather weak intensity and can therefore be associated with the formation of an ${S}'S$ boundary at the periphery of the current spread region. In many cases, the nucleation of this singularity was ascribed to the characteristic phonon energies on the $eV$ axis.
It is also worth noting that an analysis of the point contacts formed by the $Nb_3Sn$ single crystals and $Cu$ also reveals a stepwise structure of the IVC whose main features are analogous to those for the $NbSe_2$ contacts described above.


Ultra-compact (X-ray) binaries

Abstract: A short review of ultra-compact binaries, focused on ultra-compact X-ray binaries, is followed by a discussion of recent results of our VLT campaign to obtain optical spectra of (candidate) ultra-compact X-ray binaries. We find evidence for carbon/oxygen as well as helium/nitrogen discs and no evidence for (traces) of hydrogen. This suggests that the donors in the observed systems are white dwarfs. However, we also find large differences between the two C/O discs of which we have good spectra, which highlights the need for a better understanding of the optical spectra.

\section{Introduction: ultra-compact binaries}

Ultra-compact binaries are double stars with periods less than about
one hour. In such short-period binaries the components must be so
close together that ordinary, hydrogen rich, stars do not fit in
\citep[e.g.][]{nrj86}. Both stars in ultra-compact binaries thus must
be compact stars: the cores of evolved giants, white dwarfs, neutron
stars or black holes.

Of special interest are ultra-compact binaries in which mass is
transferred from one component to the other. Observationally two types
are known: ones with an white dwarf accretor (called AM CVn systems,
after the first detected system) and ones with a neutron star
accretor, called ultra-compact X-ray binaries (UCXBs).  There are
currently 16 AM CVn stars known, plus two candidate systems. For more
information we refer to a recent review \citep{nel05}. Four new
systems were reported in the last year, found in the SDSS
spectroscopic database \citep{ahh+05}. The number of known systems is
increasing rapidly: 9 of the known systems and one of the two
candidates were discovered in the last five years. The two candidate
systems V407 Vul and RX J0806.3+1527, are both found as X-ray sources,
showing a 100 per cent modulation on 9.5 and 5.4 mins
respectively. One interpretation is that these periods are orbital
periods and thus they could be part of the AM CVn family. However,
their characteristics, in particular the systematic shortening of
their periods \citep{str04a,hrw+03,rhw+05,str05}, are difficult to
reconcile with them being AM CVn stars, but see
\citet{mn05}. Alternative models that have been proposed are a face-on
longer-period binary in which we only see the spin period of an
accreting white dwarf \citep{nhw02} and a pair of detached white
dwarfs generating X-ray emission through unipolar induction
\citep{wcr02}. The latter has recently attracted new attention, with
claims that the original model cannot explain the observations
\citep{mn05} as well as suggestions that in a somewhat different form
the model can at least explain RX J0806.3+1527 \citep{dis05}.

(Ultra-compact) X-ray binaries are rare objects. The known UCXBs are much
further away than the AM CVns stars which hinders optical studies. Therefore
it is only with the advent of 8m class telescopes that high-quality optical
spectra have become available, and consequently much of our current knowledge
is based on X-ray observations and photometry at different wavelengths. There
are 12 systems with known or suggested orbital periods. Four more systems are
identified based on similarities with known systems in either their X-ray
spectra \citep{jpc00} or as a result of their optical faintness in outburst,
as expected for the small accretion discs in these systems \citep{vm94}. Six
systems are found in a globular cluster \citep[see][for reviews of globular
cluster X-ray sources]{kdi+03,vl04}.

\subsection{Why study ultra-compact binaries}

There are still many open questions in
the study of ultra-compact binaries which we would like to answer for the
following reasons:
\begin{enumerate}
\item Ultra-compact binaries are strong
gravitational-wave sources in the low-frequency regime where the ESA/NASA LISA
mission will be sensitive.
\item They are unique astrophysical laboratories because of the
combination of compact objects, short period variability and peculiar
chemical composition.
\item They are important test
of binary evolution theory, in particular the common-envelope phase.
\end{enumerate}

\begin{table*}
\caption{Overview of the UCXB (candidates). Data from \citet{vm94},
\citet{dma00}, \citet{lvv01}, \citet{kdi+03}, \citet{pou05} unless
otherwise indicated. }
\begin{tabular}{lllllllll}\hline
Name & Period & P/T & B & d & $M_V$ & \hspace*{-0.5cm}log $\frac{L_X}{L_{\rm Edd}}$ & remark & Ref. \\
 & (min) & & & (kpc) & & & & \\
\hline
4U 1820-30    & 11   & P & y & 7.6   & 3.7   & -0.7  & NGC 6624    & \\
4U 1543-624   & 18(?)& P & n & ?     & ?     & ?     & Ne-rich(?)  & 1 \\
4U 1850-087   & 21(?)& P & y & 8.2   & 5.1   & -2.4  & NGC 6712    & 2\\
M 15 X-2      & 23   & P & y & 10.3  & ?     & -2.2  & M 15        & 3, 4\\
XTE J1807-294 & 40   & T & n & ?     & ?     & ?     & msPSR         & \\
4U 1626-67    & 42   & P & n & $<4^a$& $>5.5$&$<-2.3$& \multicolumn{2}{l}{Ne-rich, slow PSR}\\
XTE J1751-305 & 42   & T & n & ?     & ?     & ?     & msPSR         &\\
XTE J0929-314 & 44   & T & n & ?     & ?     & ?     & msPSR         &\\
NGC 6652 B    & 44?  & P?& y?& 9.6   & 4.7   & -4.7  & MS colours  & 5\\
4U 1916-05    & 50   & P & y & 9     & 5     & -2    & Dipper      & \\
4U 0614+09    & 50?  & P & y & $<3$  & $>5.4$& $<-2$ & Ne-rich(?)  & 6, 7\\
XB 1832-330   &55?$^b$& P & y?& 9.6   & 3.7   & -2.3  & NGC 6652(A) & 5\\ \hline
XB 0512-401   & ?    & P & y & 12.1  &5.6$^c$& -2.0  & NGC 1851    & \\
2S 0918-549   & ?    & P & y & 4.2   & 6.9   & -2.5  & Ne-rich(?)  & 6\\
4U 1822-00    & ?    & P & n &$<20^d$&$>5.5$ & $<-1$ & Ne-rich(?)  & \\
XB 1905+000   & ?    & ? & y & 8     & 4.9   & -2.1  &             & 1\\ \hline
\end{tabular}
\textbf{Notes:} $a$ Assuming a distance from the plane $z < 1$kpc;
$b$ Or its aliases 2.2 and 4.4 hr;
$c$ $M_B$; $d$ Assuming a reasonable maximum distance. For $d = 8$kpc, $M_V$ = 7.5 and log $L_X/L_{\rm Edd} = -1.7$\\
\textbf{References:} 1 \citet{wc04} and refs therein ; 2 \citet{hcn+96}; 3 \citet{dkz05}; 4 \citet{wa01}; 5 \citet{heg01}; 6 \citet{njm+04} and refs therein; 7 \citet{ojd+05}
\label{tab:candidates}
\end{table*}

\section{A short review on ultra-compact X-ray binaries}

\subsection{Observational properties}

An overview of the observational properties of the known and candidate
UCXBs is given in Table~\ref{tab:candidates}, ordered by (suggested)
orbital period (if known). The table lists the name, the orbital
period, the type of source (persistent versus transient), the
occurrence of type I X-ray bursts, the distance, absolute V-band
magnitude, X-ray luminosity in terms of the Eddington limit (taken as
$2 \times 10^{38}$ erg/s), remarks (msPSR = millisecond pulsar, for
discussion of the meaning of Ne-rich, see Sect.~\ref{X-ray}) and
references.

\subsubsection{X-ray properties}\label{X-ray}

UCXBs have so far been studied mainly in X-rays.  We can only
summarise the results briefly. There are two types of UCXBs:
persistent and transient systems. According to \citet{db03} systems
with periods longer than about 30 minutes should be transients
\emph{if the current mass-transfer rate is similar to its long-time
average}.  The observations show that there are at least 2 and
possibly 5 persistent systems with periods longer than 30
min. However, especially in comparing persistent and transient sources
there are many biases. Interestingly, the three transients all have
periods very close together between 40 and 44 min and all three of
them harbour accretion powered millisecond X--ray pulsars.

\textbf{(millisecond) pulsations}\\ An intriguing fact is that three
of the six known accreting millisecond pulsars reside in UCXBs, e.g.
\citet{mss+02,rss02}, while only 16 UCXBs are known. 4U 1626-67 is a
peculiar system because it harbours a young neutron star with a 7s
spin period.

\textbf{X-ray spectra}\\ The X-ray spectra lacked resolution until
recent observations obtained with the gratings on board the {\it
Chandra} satellite.  The X--ray spectra are difficult to interpret due
to the lack of clear emission features (except in the 7s accreting
pulsar 4U 1626-67 which shows strong O and Ne emission lines in its
X-ray spectrum \citep{sch+01}) and contributions of interstellar
absorption. A number of UCXBs have been (first) identified based on
their X-ray spectra, when \citet{jpc00} noted a similarity between the
spectrum of 4U 1850-087 and three more LMXBs.  The common feature in
the X-ray spectra was suggested to be due to an enhancement of neon in
these systems, however, this is uncertain. These observations were
interpreted as evidence for the donors being ONe, or more likely CO
white dwarfs in which Ne has sunk to the core in an earlier
evolutionary phase and is now exposed, after accretion has peeled--off
the outer layers \citep{ynh02}. Recently it has been noted that if it
is really the O/Ne ratio that is anomalous, and not the Ne abundance
itself, this points towards a He white dwarf donor, in which the O
abundance is strongly reduced due to CNO processing
\citep{icv05}. Interestingly, one system that shows a similar feature
in its X-ray spectrum \citep[4U~1556-605][]{ffm+03}, shows strong
hydrogen and helium emission in its optical spectrum, suggesting that
it is not an UCXB (Nelemans \& Jonker, in prep).

\textbf{Type I X-ray bursts}\\ For a number of UCXBs type I X-ray
bursts have been found, especially for the ones in globular clusters,
but that may be a selection effect. The peculiar chemical composition
of the accreting material will influence the burst properties and in
principle can give an extra tool to study the chemical composition
\citep[e.g.][]{cb01}. Recently a number of peculiar bursts have been
observed in 2S~0918-549 and 4U~0614+09 \citep{icv05,kuu05} which
\citet{icv05} suggest are due to thick layers of He that are
burned. It is clear that there is still uncertainty how we should
interpret the different pieces of information, because in some cases
the bursts and the optical spectra suggest conflicting compositions as
we will discuss below!

\subsubsection{Optical and UV properties}

\begin{figure*}
  \includegraphics[angle=-90,width=\textwidth,clip]{Mvs}
  \caption{Absolute magnitudes as function of period (left) and
  $\Sigma$ (right) for a number of fiducial LMXBs and UCXBs with their
  mass transfer rate at the long-term average as expected from binary
  evolution calculations (solid dots). The observed (candidate) UCXBs
  are indicated with the plusses (period and magnitude reliable), open
  circles (periods and magnitudes not very reliable) and dashed line
  (periods unknown).  The straight line in the right panel gives the
  \citet{vm94} fit to the observations used by them.}
\label{fig:Mvs}
\end{figure*}

\textbf{Absolute magnitudes}\\ \citet{vm94} investigated the absolute
magnitudes of LMXBs, based on the assumption that most of the optical
flux is generated by X-ray irradiation of the disc. This means that
the two dominant factors determining the absolute magnitude of LMXBs
are the X-ray luminosity and the size of the disc. They derived a
parameter $\Sigma \propto L_X^{1/2} P_{\rm orb}^{2/3}$ which should
correlate with the absolute magnitude. They indeed found a strong
correlation.  Because of their small disc size UCXBs are expected to
be relatively dim. This has served as an indication for systems to be
classified as UCXB candidates. Based on the irradiation model of
\citet{ak93} we calculated the absolute magnitudes of a number of
fiducial LMXBs which we compare to the observed systems
(Fig.~\ref{fig:Mvs}).  Indeed most sources seem to be consistent with
being UCXBs, although NGC 6652 B (see Sect.~\ref{GC}) is such a faint
X-ray source that the irradiated disc model is probably not applicable
for this source and the system more likely is a quiescent LMXB
\citep[see][]{heg01}.

\textbf{Photometry}\\ The faintness of the optical counterparts has inhibited
high-quality spectra for a long time, but photometric studies have revealed
many of the know orbital periods. Recent developments in photometry are the
used of ULTRACAM to study 4U~0614+09, yielding a suggested orbital periods of
50 min \citep{ojd+05}. Recent photometry of 4U~1543-624 yielded a possible
orbital period of 18 min \citep{wc04}.

\begin{figure*}
  \includegraphics[angle=-90,width=\textwidth]{0614_ALL}
  \caption{VLT spectra of 4U 0614+09, showing lines from a
  carbon-oxygen accretion disc. From \citet{njm+04}}
\label{fig:2003}
\end{figure*}

\textbf{Spectra}\\ We recently started a systematic spectroscopic
study of (candidate) UCXBs. The main aim was to confirm/reject
candidates and to study the chemical composition of the donors stars
in these systems, which hold clues to the formation of UCXBs (see
Sect.~\ref{formation}). The first results are published in
\citet{njm+04} and are summarised in Fig.~\ref{fig:2003}. We
identified the features in the spectrum of 4U 0614+09 as relatively
low ionisation states of carbon and oxygen. This clearly identifies
this system as an UCXB and suggests the donor in this system is a
carbon-oxygen white dwarf. The similarity of the spectrum of 4U
1543-624 suggests it is a similar system, while for 2S 0918-549 the
spectrum didn't have a high enough S/N ratio to draw firm conclusions,
but it is also is consistent with being a similar system (and clearly
does not show the characteristic strong hydrogen emission lines of
low-mass X-ray binaries). We therefore concluded that all these
systems are UCXBs. In Sect.~\ref{results} we present new results.

\subsubsection{Globular cluster sources}\label{GC}

A special mention should be made about UCXBs in globular clusters
(GCs). Five of the 12 systems with known orbital periods reside in
GCs, an even larger enhancement compared to the disc than for
longer-period X-ray binaries. For reviews of X-ray sources in globular
clusters see \citet{ver04,vl04}. The shortest period system
4U~1820-30, which with an orbital period of 11 min for a long time has
been the shortest period binary known and is the most well-known
source of this class. The positional accuracy and better sensitivity
of {\it Chandra} and {\it XMM} are clearly crucial for better
understanding these systems. For instance in M 15 the bright X-ray
source has been resolved in two \citep{wa01}, the brightest one
associated with a faint FUV source that shows a clear modulation on 23
min \citep{dkz05}.  In NGC 6652 the suggested counterpart of the
brightest source turns out to be the counterpart of a (much fainter)
X-ray source \citep[source B][]{heg01}. It has a suggested periodicity
of 44 min \citep{dma00} but its colours are consistent with a
main-sequence star \citep{heg01}. The counterpart of the bright source
(source A), has a possible period of 55 min \citep[but with aliases at
2.2 and 4.4 hr][]{heg01}.

\subsection{Formation and evolution}\label{formation}

For a detailed discussion of the formation of AM~CVn systems and UCXBs
in the field we refer to \citet{npv+00,prp02} and references
therein. In short there are three routes, differentiated by the nature
of the donor star: (i) a white dwarf donor when a detached binary with
a white dwarf and either a neutron star/black hole or white dwarf
comes into contact due to angular momentum losses via gravitational-wave radiation; (ii) a (semi-degenerate) helium star donor that
evolved from a helium core burning star that filled its Roche lobe to
a white dwarf or neutron star/black hole and (iii) the core of a star
that filled its Roche lobe to a white dwarf or neutron star/black hole
at the end of the main sequence and thus has a helium rich core. In
globular clusters the formation of X-ray binaries is probably
dominated by dynamical interactions \citep[see][]{ver04}.

These different formation scenarios result in principle in different
chemical compositions, but there is some overlap. The three formation
scenarios will yield the following chemical composition:
\begin{enumerate}
\item White dwarf donors: He + CNO processed (i.e. mainly N) or C/O
\item He star donor: He + He burning products
\item Evolved secondaries: He + CNO processed + some H
\end{enumerate}

\subsection{Open questions}

The rapid increase in the number of (candidate) UCXBs in recent years
shows that this is an area where new technological developments play
an important role. However, with the increase in the sample size a
number of new questions have risen, while many of the old questions
still remain. A subjective list of the current most important open
questions is
\begin{itemize}
\item What are the chemical compositions of the donor stars and what
  are the best measurements to determine this? In particular the
  interpretation of the optical spectra (see below) and the type I
  X-ray bursts deserve attention.
\item What are the orbital periods of the remaining systems? Only when
  we have found these periods we can really start to look at the
  period distribution. The peculiar fact that most systems had an
  orbital period of about 40 min is already being diluted by new
  period determinations at 18 and 23 min. However, still the
  clustering of the three transients at very similar periods remains.
\item What is the determining factor in the difference between
  transient and persistent sources?
\item Why are so many of the accreting millisecond pulsars found in
UCXBs? Is this due to the low average mass accretion rate that is
expected in the longest-period UCXBs?
\end{itemize}

\section{VLT spectroscopy of ultra-compact X-ray binaries}\label{results}

In order to tackle some of the open questions mentioned in the last
section, we continued our spectroscopic study of (candidate) UCXBs.

\subsection{Observations}

We used the FORS2 spectrograph on the 8.2m Very Large Telescope (VLT)
of the European Southern Observatory at Paranal to obtain optical
spectra of our candidate UCXBs. The observations were taken in the
spring of 2003 and 2004.

\subsection{Results}

\begin{figure*}
\begin{center}
  \includegraphics[angle=-90,width=0.9\textwidth]{1916_models}
  \caption{VLT spectra of 4U 1916-05 (bottom). Above we show an LTE
  model with only helium and nitrogen (black line) and a similar model
  with all metals at solar abundance (grey line).}
\label{fig:1916}
\end{center}
\end{figure*}

\textbf{4U 1916-05}\\ In 2004 we obtained a spectrum of the known UCXB
4U 1916-05, which has a 50 min orbital period. In Fig.~\ref{fig:1916}
we show the spectrum together with an LTE model of a disc of pure
helium with (overabundant) nitrogen (black line) and a similar model
with metals at solar abundance (grey line). We conclude that this is
the first clear detection of a helium donor in an optical
spectrum. However, the spectrum does not look at all like we had
anticipated (i.e. strong helium emission lines, like the spectra of AM
CVn stars ES Cet and GP Com). This means that really high S/N spectra
are needed to determine the chemical composition of the donors and
that in our earlier results, the spectrum of 2S 0918-549 clearly is
not good enough to distinguish between helium and C/O and might be
more likely a helium donor because of the properties of its X-ray
bursts \citep{icv05}.

\begin{figure*}
\begin{center}
  \includegraphics[angle=-90,width=0.8\textwidth]{0929}
  \caption{Spectrum of XTE J0929-314 taken during outburst. Above and
  below an LTE model for a CO (top) and a He/N (bottom) disc are
  shown.}
\label{fig:0929}
\end{center}
\end{figure*}

\begin{figure}
\begin{center}
  \includegraphics[angle=-90,width=\columnwidth,clip]{0929_Halpha}
  \caption{Region around H$\alpha$ of the outburst spectrum of XTE
  J0929-314 on the basis of which \citet{ccg+02} claim H$\alpha$
  emission which would contradict the interpretation of the donors as
  white dwarfs. The 3 individual spectra are shown at the bottom, the
  average as the top spectrum.}
\label{fig:0929_Ha}
\end{center}
\end{figure}

\textbf{XTE J0929-314}\\ XTE J0929-314 is one of the three transient UCXBs
and was discovered in 2002 \citep{rss02}. In outburst a spectrum was
taken with the ESO 3.6m telescope in which emission at H$\alpha$ and
around 4650\AA~was reported \citep{ccg+02}. In particular the
H$\alpha$ emission is important, since this would clearly point to a
donor with still some hydrogen left. However, the evidence for
H$\alpha$ emission is rather weak as is shown in
Fig.~\ref{fig:0929_Ha}, where we display the spectra. The formal
significance of the H$\alpha$ emission is 3$\sigma$, but one has to
keep in mind that there are many emission lines in the night sky
around H$\alpha$ that could lead to systematic errors. The emission
around 4650\AA~is clearly significant, as seen in Fig.~\ref{fig:0929},
where we also plot simple C/O and He LTE model. The spectrum is not
good enough to determine the composition, although the lack of
emission around 4500\AA~suggest a C/O donor.

\begin{figure*}
 \begin{center}
  \includegraphics[angle=-90,width=0.8\textwidth]{1626_0614}

  \caption{Spectra of 4U 1626-67 and 4U 0614+09. The spectra are quite
  similar above 4500\AA, but below they are very different and
  interestingly the spectrum of 4U 1626-67 follows the LTE model that
  we used for the line identification of 4U 0614+09 much better than
  4U 0614+09 itself!}
\label{fig:1626_0614}
 \end{center}
\end{figure*}

\textbf{4U 1626-67 and 4U 0614+09}\\ We also obtained a spectrum of 4U
1626-67. In \citet{nj04} we already showed a comparison between the
spectrum of 4U 1626-67 and our earlier 4U~0614+09 spectrum which
looked very much like each other, in agreement with the interpretation
of the donors of both these systems being C/O rich. Indeed 4U 1626-67
follows the simple C/O LTE model that we used for 4U 0614+09
reasonably well even at bluer wavelengths than out 4U 0614+09 spectrum
reached. However, a second spectrum of 4U 0614+09 (obtained with the
VLT by Klaus Werner, Werner et al.~in prep.) with the same setup and
wavelength coverage as our 4U 1626-67 spectrum clearly does
\emph{not} follow the simple LTE model! (see
Fig.~\ref{fig:1626_0614}). It is clear that the last word on the
interpretation of these spectra is not yet said.

Ultra-compact (X-ray) binaries


Hybrid Session-based News Recommendation using Recurrent Neural Networks

Abstract: We describe a hybrid meta-architecture -- the CHAMELEON -- for session-based news recommendation that is able to leverage a variety of information types using Recurrent Neural Networks. We evaluated our approach on two public datasets, using a temporal evaluation protocol that simulates the dynamics of a news portal in a realistic way. Our results confirm the benefits of modeling the sequence of session clicks with RNNs and leveraging side information about users and articles, resulting in significantly higher recommendation accuracy and catalog coverage than other session-based algorithms.

\section{Introduction}

Recommender systems help users to deal with information overload by providing tailored item suggestions to them. One of the earliest application domains is the recommendation of online \textit{news} \cite{karimi2018news}. News recommendation is sometimes considered as being particularly difficult, as it has a number of distinctive characteristics \cite{Zheng:2018:DDR:3178876.3185994}. Among other challenges, news recommenders have to deal with a constant stream of news articles being published, which at the same time can become outdated very quickly. Another challenge is that the system often cannot rely on long-term user preference profiles. Typically, most users are not logged in and their short-term reading interests must be estimated from only a few logged interactions, leading to a \textit{session-based recommendation problem} \cite{QuadranaetalCSUR2018}.

In recent years, we observed an increased interest in the problem of session-based recommendation, where the task is to recommend relevant items given an ongoing user session. Recurrent Neural Networks (RNN) represent a natural choice for sequence prediction tasks, as they can learn models from sequential data. \textit{GRU4Rec} \cite{hidasi2016} was one of the first neural session-based recommendation techniques, and a number of other approaches were proposed in recent years that rely on deep learning architectures, as in \cite{Liu2018stamp, Li2017narm}.

However, as shown in \cite{jannach2017recurrent,ludewig2018evaluation,LudewigMauro2019}, neural approaches that only rely on logged item interactions have certain limitations and they can, depending on the experimental setting, be outperformed by simpler approaches based, e.g., nearest-neighbor techniques.
Differently from previous works, we therefore leverage multiple types of side information with RNNs, including textual article embeddings, as well as the context of users and articles. Furthermore, we propose a meta-architecture to address the aforementioned challenges of recommending in the news domain.

\section{Technical Contribution}
Our approach is based upon CHAMELEON \cite{moreira2018chameleon, moreira2018news, moreira2019inra, moreira2019contextual}, which is a Deep Learning Meta-Architecture for News Recommendation. It supports session-based news recommendation scenarios, modeling the sequence of user clicks using Recurrent Neural Networks. The resulting system is a hybrid recommender system, which addresses the permanent user and item cold-start problem in the news domain by leveraging the textual content of news articles, the article context (e.g., recent popularity and recency) and the user context (e.g., time, location, device, previous session clicks).

Figure~\ref{figure:chameleon_instantiation} shows our instantiation of the \emph{CHAMELEON} framework with its two main modules:  the \emph{ACR} module on the left creates distributed representations of articles' textual content. The \emph{NAR} module on the right is responsible to generate next-click predictions. The \emph{NAR} module is trained on a ranking loss function based on similarities, which is designed to recommend fresh articles without retraining. As proposed for the \textit{DSSM} loss function \cite{huang2013learning}, it is trained to maximize the likelihood of correctly predicting the next clicked article given a user session.

\begin{figure}[h!t]
	\centering\includegraphics[width=1.0\linewidth]{figures/chameleon_instantiation_journal.jpg}	\caption{An architecture instantiation of \textit{CHAMELEON}}
	\label{figure:chameleon_instantiation}
\end{figure}

\section{Evaluation Protocol}

The evaluation was performed as follows: (1) Recommenders are continuously trained on users' sessions ordered by time and grouped by hours. Each five hours, the recommenders are evaluated on sessions from the next hour; (2) For each session in the evaluation set, we incrementally revealed one click after the other to the recommender; and (3) For each click to be predicted, we created a set containing 50 negative samples articles (not clicked by the user in her session) and compute top-N metrics about accuracy, item coverage, and novelty.

\begin{figure}[h!t]
	\centering\includegraphics[width=0.99\linewidth]{figures/chameleon_eval_protocol.jpg}	\caption{Illustration of the evaluation protocol. After training for 5 hours, we evaluate using the sessions of the next hour.}
	\label{figure:eval_protocol}
\end{figure}

Experiments were performed with two public real-world datasets from the \emph{G1} \cite{moreira2018news} and \emph{Adressa} \cite{gulla2017adressa} news portals, described in Table~\ref{tab:datasets}.

We evaluated the following recommendation quality factors for the top-N ranked items: accuracy -- Hit Rate (\emph{HR@n}) and Mean Reciprocal Rank (MRR@n); item coverage -- (\emph{COV}) (i.e., the number of distinct articles that appeared in any top-N list divided by the number of recommendable articles); and novelty -- \emph{ESI-R}, which is based on item popularity, returning higher values when recommending long-tail items.

As baseline algorithms for session-based recommendation, we have used: two neural approaches (\emph{GRU4Rec} \cite{hidasi2016} and \emph{SR-GNN} \cite{wu2019session}); association rules-based methods (\emph{Co-Occurrence (CO)} and \emph{Sequential Rules (SR)} \cite{ludewig2018evaluation}); neighborhood-based methods (\emph{Item-kNN} \cite{hidasi2016} and \emph{Vector Multiplication Session-Based kNN (V-SkNN)} \cite{jannach2017recurrent}); and two other classical methods (\emph{Recently Popular (RP)} \cite{ludmann2017recommending} and \emph{Content-Based (CB)}.

\begin{table}[h!t]
\centering
\footnotesize
\caption{Statistics of the datasets used for the experiments.}
\label{tab:datasets}
\vspace{10pt}
\begin{tabular}{p{3.6cm}rr}
\hline
 & \textit{Globo.com (G1)}
 & \textit{Adressa} \\ \hline
Language  & Portuguese & Norwegian  \\
Period (days)  & 16 & 16 \\
\# users   & 322,897 & 314,661  \\
\# sessions & 1,048,594 & 982,210 \\
\# clicks  & 2,988,181 & 2,648,999 \\
\# articles   & 46,033 & 13,820  \\

Avg. Sessions length (clicks)  & 2.84  & 2.70  \\
\hline
\end{tabular}

\end{table}

\section{Results}

The evaluation results are presented in Table~\ref{tab:metrics_results}, as originally reported in \cite{moreira2019contextual}. The best results for a metric are printed in bold face and marked with * if they are significantly different \footnote{As errors around the reported averages were normally distributed, we used paired Student's t-tests with Bonferroni correction at $p<0.001$ for significance tests.} from all other algorithms.

\begin{table}[h!t]
\centering
\caption{Evaluation of recommendation quality factors}
\label{tab:metrics_results}
\footnotesize
\begin{tabular}{p{1.8cm}p{1.0cm}p{1.0cm}p{1.1cm}p{1.4cm}}
\hline
 \textit{Recommender}
 & \textit{HR@10}
 & \textit{MRR@10}
 & \textit{COV@10}
 & \textit{ESI-R@10}
 \\
\hline
\multicolumn{5}{p{4cm}}{\textbf{G1 dataset}}\\ \hline
\textit{CHAMELEON} & \textbf{0.6738}* & \textbf{0.3458}* & 0.6373 & 6.4177   \\
\textit{SR} & 0.5900 & 0.2889 & 0.2763 & 5.9747 \\
\textit{Item-kNN} & 0.5707 & 0.2801 & 0.3913 & 6.5909   \\
\textit{CO} & 0.5689 & 0.2626 & 0.2499 & 5.5728 \\
\textit{V-SkNN} & 0.5467 & 0.2494 & 0.1355 & 5.1760  \\
\textit{SR-GNN} & 0.5144 & 0.2467 & 0.3196 & 5.4280  \\
\textit{GRU4Rec} & 0.4669 & 0.2092 & 0.6333 & 5.2332   \\
\textit{RP} & 0.4577 & 0.1993 & 0.0218 & 4.4904  \\
\textit{CB} & 0.3643 & 0.1676 & \textbf{0.6774} & \textbf{8.1531}*  \\
\hline
\multicolumn{5}{p{4cm}}{\textbf{Adressa dataset}}\\ \hline
\textit{CHAMELEON} & \textbf{0.7018}* & \textbf{0.3421}* & 0.7926 & 5.3410  \\
\textit{SR} & 0.6288 & 0.3022 & 0.4604 & 5.4443  \\
\textit{Item-kNN} & 0.6179 & 0.2819 & 0.5314 & 5.4675  \\
\textit{CO}  & 0.6131 & 0.2768 & 0.4220 & 5.0789  \\
\textit{V-SkNN} & 0.6140 & 0.2723 & 0.1997 & 4.6018 \\
\textit{SR-GNN} & 0.6122 & 0.2991 & 0.5197 & 5.1013  \\
\textit{GRU4Rec} & 0.4958 & 0.2200 & 0.5143 & 5.0571   \\
\textit{RP} & 0.5648 & 0.2481 & 0.0542 & 4.1465  \\
\textit{CB}  & 0.3307 & 0.1253 & \textbf{0.8875}* & \textbf{7.6715}* \\
\hline
\end{tabular}
\end{table}

\section{Conclusion}

CHAMELEON was specifically designed to address news recommendation challenges such as (a) the short lifetime of the recommendable items and (b) the lack of longer-term preference profiles of the users.
In the extensive experiments performed, CHAMELEON was able to provide recommendations with much higher accuracy than all other evaluated algorithms, and it led to the second best item coverage.

Machine Learning (cs.LG); Information Retrieval (cs.IR); Machine Learning (stat.ML)