High-Energy Neutrino Astronomy: where do we stand, where do we go?

Abstract: With the identification of a diffuse flux of astrophysical ("cosmic") neutrinos in the TeV-PeV energy range, IceCube has opened a new window to the Universe. However, the corresponding cosmic landscape is still uncharted: so far, the observed flux does not show any clear association with known source classes. In the present talk, I sketch the way from Baikal-NT200 to IceCube and summarize IceCube's recent astrophysics results. Finally, I describe the present projects to build even larger detectors: GVD in Lake Baikal, KM3NeT in the Mediterranean Sea and IceCube-Gen2 at the South Pole. These detectors will allow studying the high-energy neutrino sky in much more detail than the present arrays permit.

\section{Introduction}

The first conceptual ideas how to detect high energy neutrinos date back to the late fifties.
The long evolution towards detectors with a realistic discovery potential
started in the seventies, by the pioneering works
in the Pacific Ocean close to Hawaii (DUMAND).
The DUMAND 1978 design envisaged an array of about 20\,000 photomultipliers
spread over a 1.26 cubic kilometer volume of water.
This project was terminated in 1995, but the baton was
taken by the projects NT200 in Lake Baikal, AMANDA at the South Pole,
ANTARES in the Mediterranean Sea and, again at the South Pole, IceCube (see for detailed
information on the history and on corresponding references \cite{Spiering-History}).
But only now, half a century after the first concepts, a cubic kilometer detector
is in operation: IceCube at the South Pole.
With the discovery of a flux of high-energy neutrinos of astrophysical origin (``cosmic neutrinos'')
in 2013 \cite{Science-2013}, the IceCube Neutrino Observatory has opened a new window to the Universe of non-thermal cosmic processes.
A next generation of arrays is under construction or planned: KM3NeT in the
Mediterranean Sea \cite{KM3NeT}, the Gigaton Volume Detector GVD in Lake Baikal
\cite{Baikal},
and IceCube-Gen2 \cite{Gen2}.

The primary goal of these detectors is identifying the sources of high-energy cosmic rays.
In contrast to charged particles, neutrinos are not deflected in cosmic magnetic fields and keep their direction; in contrast to gamma rays they provide a direct, water-tight prove for the acceleration of hadrons in the emitting sources. This makes them unique tracers of sources of cosmic rays.
On the other hand, due to their small interaction cross section they are difficult to detect: The         ``neutrino effective area'' of the 1 km$^3$ IceCube detector (essentially the geometrical area multiplied with the interaction probability, the trigger efficiency and the transparency to neutrinos of the Earth)
is less than 1\,m$^2$ at 1\,TeV and of the order of 100\,m$^2$ at 100\,TeV \cite{point-sources}. It is therefore no surprise that it took several decades to detect cosmic neutrinos.

Neutrino telescopes are multi-purpose detectors. Apart from investigating cosmic neutrinos, they exploit atmospheric neutrinos to study neutrino oscillation, to search for sterile neutrinos or to test fundamental laws of physics. They are used to search for neutrinos from Dark Matter annihilations in the Sun or the Galactic halo, to search for exotic particles like magnetic monopoles, or to study muons from cosmic-ray induced air showers.

This paper focuses to the search for neutrinos from cosmic acceleration processes. I will further focus to the developments at Lake Baikal and at the South Pole where I have been involved myself over a long period.

\section{From Baikal NT200 to IceCube}

\subsection{NT200 and AMANDA}

First test deployments in Lake Baikal started in 1981.
The construction of the NT200 detector was started in 1993,
about 30 km South-West from the outflow of Lake Baikal into
the Angara river, at a distance of
3.6\,km to shore and at a depth of about 1.1\,km, and was completed in 1998.
NT200 was an array of 192 optical modules
(glass spheres containing large photomultipliers)
at eight strings, 72\,m in height and 43\,m in diameter
(see Fig.\,\ref{nt200}). Actually,
this is not much more than twice the size of Super-Kamiokande. First
upward going muons (i.e.\,neutrino events) were found already with the
3-string version from 1993, and then with the 4-string version from 1996:
the first ``underwater neutrinos'' ever! This was the first proof of principle
to detect neutrinos in open media and a breakthrough for the field.

\begin{figure}[ht]
\includegraphics[width=0.80\linewidth]{figures/nt200andgoldplated.pdf}
\caption{\small {\it Left}: The Baikal Neutrino Telescope NT200.
{\it Right}: One of the first upward moving
muons from a neutrino interaction recorded with the 4-string stage of the
detector in 1996 \cite{Baikal-atm-Balkanov-1999}. The Cherenkov light from the
muon is recorded by 19 channels.}
\label{nt200}
\end{figure}

In 1988, a new, spectacular idea appeared on stage: to use Antarctic ice instead
of water as target and as detector medium.
The project was named AMANDA (Antarctic Muon And Neutrino Detection Array).
AMANDA was deployed some hundred meters from the Amundsen-Scott station, first
at a too shallow depth (where bubbles disturb light propagation), and then
from 1996 -- 2000 at depths between 1500 and 2000\,m. It consisted of
677 optical modules at 19 strings.
AMANDA was switched off in April 2009, after more than 9 years of data
taking in its full configuration and with 6959 neutrino events collected. Naturally
this sample was dominated by atmospheric neutrinos.
No indication of point sources was found, and no excess of high-energy
events which might have pointed to an admixture of cosmic neutrinos.
Figure \ref{A-B-skyplot}
shows -- as a kind of sentimental reminiscence -- a combination of NT200 and AMANDA
data compiled in 2005 (including all NT200 data and 2 years AMANDA data).

\begin{figure}[ht]
\centering
\includegraphics[width=0.65\linewidth]{figures/Amanda+Baikal.png}
\caption{\small A skyplot in Galactic coordinates, combining neutrinos detected by NT200 (circles) and AMANDA (crosses).
Compiled in 2005 \cite{Igor-Christian-skyplot}.
}
\label{A-B-skyplot}
\end{figure}

AMANDA provided record limits on fluxes for cosmic neutrinos,
be it for diffuse fluxes (where the much smaller NT200 could compete for some years)
and for point sources -- steady as well as transient.
AMANDA also extended the measured spectrum of atmospheric neutrinos by nearly
two orders of magnitude, from a few TeV to 200\,TeV.

\subsection{Mediterranean Projects}

First site studies in the Mediterranean Sea were performed in 1989,
leading to the NESTOR project, with the goal to install towers
of hexagonal floors close to Pylos in Greece. The first floor was deployed only 15 years
later, in 2004.
But the project was further and further delayed and the community split into
a Greek project (NESTOR), an Italy-based project (NEMO) and a
project in France (ANTARES). Only the French site made it to
a working -- and actually excellently working! -- detector.

The construction of ANTARES started in 2002 with the deployment of a shore cable.
The detector in its final 12-string configuration was
installed in 2006--2008 and has been operational since then.
The strings have lateral distances of 60--70\,m, and each of them carries
25\,triplets of optical modules at depths of 2.1--2.4\,km.

ANTARES has demonstrated that a stable operation of a deep-sea detector is possible.
Similar in size to AMANDA,  it has collected more than 8000 upward-going
muon tracks over eight years of operation.
With its excellent view of the Galactic plane and
good angular resolution, the telescope could constrain the Galactic origin of the
cosmic neutrino flux reported by IceCube. ANTARES has explored he Southern sky and in particular
central regions of our Galaxy searching for point-like objects,
for extended regions of emission, and for signals from transient objects selected
through multi-messenger observations
\cite{Antares-results}.

\subsection{IceCube}

Like AMANDA, the IceCube Observatory \cite{IceCube} is located at the geographical South Pole. It consists of the main IceCube array with its subarray DeepCore and the surface array IceTop. IceCube comprises 5160 digital optical modules (DOMs) installed on 86 strings at ice
depths of 1450 to 2450\,m and covers 1 km$^3$ of ice.
A string carries 60~DOMs. DeepCore, a high-density sub-array of eight strings
at the center of IceCube, has smaller spacing and DOMs with more sensitive photomultipliers than
IceCube and sits in the midst of the clearest ice layers. This results in a
threshold of about 10\,GeV and opens a new venue for oscillation physics.
The threshold of the full IceCube detector is about 100\,GeV. In its final configuration,
IceCube takes data since spring 2011, with a duty cycle of more than 99\%. It collects about $10^5$ clean neutrino events per year, with nearly
99.9\% of them being of atmospheric origin.

\section{Where do we stand?}

\subsection{Diffuse Fluxes}

It has been predicted since long that the first evidence for extragalactic cosmic neutrinos would be provided by a diffuse flux rather than by single-source signals \cite{Lipari}.
The first tantalizing hint to cosmic neutrinos in IceCube came from two shower-like events with energies
$\approx 1$ PeV, discovered in 2012 and dubbed ``Ernie'' and ``Bert'' \cite{Ernie}.
A follow-up search of the same data (May 2010 to April 2012) with a lowered threshold (30 TeV)
provided 25 additional events. This analysis used only events starting in a fiducial volume of
about 0.4 km$^3$ (High Energy Starting Events, or ``HESE''), using the other 60\% of IceCube as veto against all sorts of background. Energy spectrum and zenith angle distribution
of the 27 events excluded an only-atmospheric origin with $4.1\sigma$ and suggested that
about 60\% were of cosmic origin,
at energies above 100 TeV even about 80\% \cite{Science-2013}. A four-year data set with
54 neutrinos provided another shower-like PeV event (deposited energy $\approx$ 2 PeV) and
confirmed a dominant cosmic contribution with nearly $6.5\sigma$. Very recently, the results from
a six-year sample have been presented \cite{ICRC-Kopper}, with 82 events above 30 TeV.
Figure 3
shows the  energies deposited by these events inside IceCube.

\begin{figure}[h]
\vspace{-0.1cm}
\centering
\includegraphics[width=0.5\textwidth]{figures/HESE-spectrum.pdf}
  \caption{\small
	Distribution of the energy deposited by 82 events from the six-year HESE analysis.
	Backgrounds of atmospheric origin
come from punch-through down-going muons and from atmospheric neutrinos.
While the flux of neutrinos from $\pi$ and K decays is well known (blue region), the neutrino flux
from charm decays in the atmosphere is uncertain and dominates
the uncertainty of all background sources (gray region with 1$\sigma$ uncertainties). The best-fit astrophysical spectra are shown as gray lines, for a single power-law spectrum as solid line, for a two power-law model as dashed line. See \cite{ICRC-Kopper} for details.
}
\label{l-HESE-spectrum}
\end{figure}

A $5.6\sigma$ excess of high-energy cosmic neutrinos is also seen in the spectrum of secondary muons generated by neutrinos that have traversed the Earth, with a zenith angle less than 5 degrees above the horizon
(``upward through-going muons''\cite{muon-7years}).
Figure 4
shows the median neutrino energy. It is calculated for each energy deposited by the muon in the detector, assuming the best-fit spectrum. The highest energy muon has deposited $2.6\pm 0.3$ PeV inside the instrumented volume, which corresponds to a most probable neutrino energy of about 9 PeV.

\begin{figure}[ht]
\vspace{-0.1cm}
\centering
\includegraphics[width=0.65\textwidth]{figures/muon-energy-spectrum.png}
  \caption{\small
Spectrum of the median neutrino energy derived from the energy deposit of through-going muons with zenith angles less than 5 degrees above horizon (8 years sample \cite{muon-8years}).
}
\label{l-muon-energy-spectrum}
\vspace{-0.1cm}
\end{figure}

While both analyses (HESE and through-going muons)  have reached a significance for a strong
non-atmospheric contribution of more than 5\,$\sigma$, the spectral indices
$\gamma$ of the astrophysical flux from both analyses disagree:  $\gamma = 2.92 \pm 0.33/0.29$ for the HESE events
(unbroken spectrum $E^{-\gamma}$) and $\gamma = 2.19 \pm 0.10$ \cite{muon-8years} for the throughgoing muons. Adding two more years to the HESE sample has resulted in an even softer energy spectrum since all events of the recent two years have energies below 200 TeV.  Fig.\,\ref{l-spectrum-combined}a
shows the two fits under the assumptions of a single-power law. The possibility that all but the three PeV HESE events emerge from pion/Kaon/charm decays in the atmosphere is excluded by the zenith angle distribution.

In \cite{combined} the flavor ratio of the astrophysical neutrino flux has been investigated. It is consistent with an observed flavor ratio $\nu_e : \nu_{\mu} : \nu_{\tau} = 1 : 1 : 1$ and also with source neutrino ratios 1:2:0 (pion decay) and 0:1:0 (pion decay with suppressed muon decay) while largely excluding 1:0:0 (neutrinos from neutron decay).

The hope to see any clustering of the HESE and muon-track events at highest energies has not fulfilled. An initial indication of clustering of HESE events close to the Galactic center has vanished with more statistics.
In addition, ANTARES has looked to a point source at the position of IceCube's initial excess and could
exclude that it is due to a point source, assuming that the extension of the source does not exceed
0.5 degrees
and that the spectrum follows an $E^{-2}$ shape \cite{Antares-check}.

A recent IceCube analysis has used 7 years of the medium-energy $\nu_{\mu}$ data (which are optimized to search for point sources, see next section) to set constraints on the diffuse emission of neutrinos from the Galactic plane  \cite{ICRC-Galacticplane}. The resulting limits are shown in
Fig.\,5b
and compared to the  flux of the HESE and highest-energy $\nu_{\mu}$ data. They exclude that more than 14\% of the observed diffuse astrophysical flux come from the Galactic plane. However, the limit is not far from model predictions (gray band). Joining IceCube and ANTARES data and exploiting cascade-like events in addition to the $\nu_{\mu}$ sample may drive the sensitivity into the region predicted by KRA models.

\vspace{-0.4cm}
\begin{figure}
\centering
\includegraphics[width=0.7\textwidth]{figures/spectrum-combined.pdf}
\vspace{0.3cm}
\includegraphics[width=0.7\textwidth]{figures/galplane-first.pdf}
\vspace{-0.5cm}
  \caption{\small
{\it Top}: Best-fit of the per-flavor neutrino fluxes as a function of energy. The black points with $1\sigma$ uncertainties are extracted from a combined likelihood fit of all background components together with an astrophysical flux component with an independent normalization in each band (assuming an $E^{-2}$ spectrum within each band and atmospheric neutrino and muon fluxes subtracted). The best-fit conventional flux and the upper limit for prompt neutrinos are shown separately, not taking into account the HESE self-veto which actually reduces their contribution.
The blue band shows the $1\sigma$ uncertainties of a single power-law fit to the HESE data. The pink band shows the fit for the muon neutrino data, again with $1\sigma$ uncertainties. Its length indicates the approximate range providing 90\% of the significance of this analysis \cite{ICRC-Kopper}.
{\it Bottom}:
Upper limits on the three flavor neutrino flux from the Galaxy with respect to KRA model predictions \cite{KRA} and the measured astrophysical flux \cite{ICRC-Galacticplane}. Dots, yellow and green bands have the same meaning as the bands in the top figure.
}
\label{l-spectrum-combined}
\end{figure}

\subsection{Search for steady point sources}

For the standard steady-source search, a sample of through-going muons with good angular resolution
(median error smaller $1\deg$) is selected. In the lower hemisphere, the Earth acts as filter against
muons generated in the atmosphere. In the upper hemisphere, a radical energy cut removes most of the
atmospheric muons which have a rather soft energy spectrum, but naturally also rejects all but the most energetic cosmic neutrinos. Therefore only hard-source spectra would result in a significant number of events from the upper hemisphere (for IceCube: South).

Figure \ref{l-skyplot7years} shows the all-sky plot of seven years
data, with 422 791  upward muons from neutrino interactions
and 289 078 downward muons, the latter almost all from atmospheric showers. The downward sample
contains also 961 tracks starting inside the detector, i.e. generated in neutrino interactions \cite{point-sources}.

\begin{figure}[h]
\centering
\includegraphics[width=0.85\textwidth]{figures/skyplot7years.png}
\vspace{-0.1cm}
  \caption{\small
	All-sky plot of seven years IceCube data in equatorial coordinates. Shown is the
	negative logarithm of the pre-trial p-value, assuming no clustering as null-hypothesis
	\cite{point-sources}.
	}
\label{l-skyplot7years}
\end{figure}

No significant excess is found, resulting in the flux constraints show in Figure \ref{l-pt-limits}. Apart from sensitivities and limits for selected sources, the discovery potential is shown, i.e. the flux that would lead to a $5\sigma$ discovery of a source in 50\% of the cases.

One can then compare these values to predictions for selected sources.
Fig.\,\ref{l-limitsblazars} compares our sensitivities and the obtained 90\% upper limits to predictions \cite{Petropoulou} for three blazars. The limits are within a factor 5 of the predictions, for Mkr 421 even slightly below predictions. Similar relations hold for the Crab nebula -- always optimistically assuming that the gamma flux observed from these sources is basically due to $\pi^0$ decay and not to inverse Compton scattering.

\begin{figure}[h]
\centering
\includegraphics[width=0.8\textwidth]{figures/limitsblazars.png}
	\vspace{-0.1cm}
	\caption{\small
	Discovery potential and sensitivity (red solid and dashed, respectively) versus declination, assuming an unbroken
	$E^{-2}$ neutrino spectrum. Upper limits of 32 pre-selected source candidates
	are given as red crosses, the blue line represents the upper limit for the most significant spots in each half of the sky (actual positions of the spots are given by blue stars). The gray line shows the results from ANTARES.
	See \cite{point-sources} for details.
	}
\label{l-pt-limits}
\end{figure}

\begin{figure}[h]
\vspace{-0.1cm}
\centering
\includegraphics[width=0.75\textwidth]{figures/limits-vs-declination.png}
\vspace{-1cm}
\caption{\small
Differential energy spectra versus neutrino energy for blazars of the BL Lac type compared to model predictions
\cite{Petropoulou}. Thick lines give the 90\% upper limits from IceCube, thin lines represent the model. The sensitivities of the Icecube search are shown as dashed line. 90\% upper limit and sensitivity are shown for the energy interval where 90\% the events originate that are most signal-like \cite{point-sources}.
}
\label{l-limitsblazars}
\end{figure}

From these figures one could conclude that an improved angular reconstruction and twice more data could bring us close to discovery. For blazars, however, this hope is downsized by various
blazar stacking analyses, none of them yielding an excess in the directions of blazars. The most recent one \cite{blazarstacking} indicates that only 4-6\% of the observed diffuse astrophysical muon neutrino flux could come from blazars.

\subsection{Search for transient sources}

To improve the signal-to-background ratio one can search for transient signals, preferentially
in coincidence with an observation in electromagnetic waves. Examples are flares of Actice Galactic Nuclei (AGN) or Gamma-Ray Bursts (GRB). GRBs are interesting objects since there
are models which assume that they are the dominant source of
the measured cosmic-ray flux at highest energies, either by neutron escape \cite{Ahlers}
or by escape of both neutrons and protons \cite{WB-GRB} from the relativistic fireball. Naturally models where protons a kept in the acceleration region and only neutrons escape and constitute the observed cosmic ray flux give a higher neutrino/cosmic ray ratio.

All three collaborations -- Baikal, ANTARES and IceCube -- have searched for neutrinos in local and spatial
coincidence with GRBs.  In particular
IceCube limits on neutrinos from GRBs have drastically improved over the recent years. A recent analysis has combined the searches for spatial and temporal coincidences
of upward and downward tracks and cascade-type events with 1172 GRBs. No  significant correlations between the gamma-ray signals and neutrinos have been observed. Figure \ref{l-GRB} shows exclusion contours
for double broken power-law spectra, with breaks from $E^{-1} \times \epsilon_b$ to $E^{-2}$ at
energy $\epsilon_b$, and from $E^{-2}$ to $E^{-4} \times (10\epsilon_b)^2$ at an
energy $10\epsilon_b$.

\begin{figure}[h]
\centering
\includegraphics[width=0.65\textwidth]{figures/GRB-limits-WB-1.pdf}
  \caption{\small
Excluded regions for 99\%, 90\% and 68\% confidence level of the generic
double broken power law neutrino spectrum as a function of the
first break energy $\epsilon_b$ and per-flavor quasi-diffuse flux normalization
derived from upward and downward muon tracks and all-sky cascades.
	}
\label{l-GRB}
\end{figure}

Both models, those with cosmic ray escape via neutrons and those which
allow additionally for cosmic ray escape via protons,
are excluded at over 90\% confidence level, with most of
the model assumption phase space excluded at over the
99\% confidence level, greatly constraining the hypothesis that GRBs
are significant producers of ultra-high energy cosmic rays in the prompt GRB
phase.

30 years after the discovery of supernova SN1987 it is worth to highlight that IceCube
runs a supernova trigger, with a duty time of more than 99\%. The trigger reacts to
a collective rise in photomultiplier counting rates on top of the dark-noise rate.
This rise would be due to the feeble signals
from $\nu_e$
reactions close to a photomultiplier. A 1987A-type supernova
at 30 kpc distance (edge of the Galaxy) would lead to an collective-rate enhancement with
a significance of about 20 standard deviations, and even at distance of the Large Magellanic Cloud
(50 kpc) the excess would reach $6-7\sigma$ \cite{SN}. IceCube is part of the SuperNova Early Warning
System (SNEWS \cite{SNEWS}), together with the underground neutrino detectors Borexino, Super-Kamiokande,
LVD and Kamland which, in the case of a significant coincidence from more than one of the detectors,
would alarm the astronomers community. However, no significant neutrino signal has been recorded yet,
neither with the analogous trigger of IceCube's predecessor AMANDA nor with that of IceCube.

\subsection{Real-time alert and follow-up programs}

With no steady sources of high-energy neutrinos observed so far, neutrinos produced
during transient astrophysical events are a viable alternative. High-energy neutrinos from the
prompt phase of GRBs or MeV neutrinos from a supernova collapse as discussed in the previous
section are just two examples. Coincident detections could enhance the significance of the
IceCube observation and, more generally, contribute to the mosaic of informations from different
messengers, providing a more complete picture of the source. Since IceCube and ANTARES have nearly
$4\pi$ acceptance (depending on energy), they could could trigger detections with pointing
devices like optical or gamma-ray telescopes, which otherwise would have been missed.

Both collaborations run a number of high-energy alert and follow-up programs \cite{alerts,alertsAnt}
which react to particular single events.
In the case of IceCube, neutrino alert candidates are identified in real-time at the South Pole. A brief message sent to
the North is automatically issued to the Gamma-Ray Coordinates Network (GCN \cite{AMON}) via the
Astrophysical Multimessenger Observatory Network (AMON \cite{AMON}). In parallel, quality checks are
applied and the directional and energy reconstruction refined. Results from that are completed
within a few hours and lead to an updated alert notification in the form of a CGN circular.
IceCube runs two of these alerts: a ``HESE Track alert'' which is issued if a track-like HESE event is recorded
(4.8 events expected per year, with 1.1 being of astrophysical origin) and
an ``EHE Track Alert''
which is based on a selection which originally targets cosmological neutrinos (10 PeV to 1 EeV) but
here is modified to be sensitive down to 500 TeV (about 5 alerts per
year).

Apart from these public alerts, IceCube also issues alerts
to optical, X-ray and gamma-ray observatories which are based on neutrino {\it multiplets}.
These alerts are based on individual agreements with these observatories.
The multiplets can be due to phenomena on the second-to-minute scale (high-energy neutrinos from relativistic jets in SN or GRB), or to phenomena of the hour-to-week scale
(like AGN flares).
None of the alerts yet has led to a significant correlation, although at least two cases have generated some
initial excitement.
The one \cite{doublet} was a neutrino doublet detected in March 2012 which triggered follow-up observations by the Palomar Transient Factory (PTF). PTF found a Type IIn supernova within an error radius of $0.54\deg$ of the direction of the doublet. A Pan-STARRS1 survey, however, showed that its explosion time was
at least 158 days before the neutrino alert,
so that a causal connection is unlikely.
The second case \cite{triplet} was the first triplet: three neutrinos arriving within 100 s of
one another at February 17, 2016. Follow-up observations by SWIFT's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime did not detect any likely electromagnetic counterpart. In a refined reconstruction, the directions of the events changed
slightly, so that the triplet turned to a double-doublet (error circle of the one event overlapping with those of the two others, but not all three with each other).
Still, these two cases impressively illustrate the potential of and challenges for future follow-up campaigns. Although no significant correlations have been detected so far, the IceCube/ANTARES alerts and the triggered electromagnetic-domain observations herald the era of multi-messenger observation. This remark also applies to the follow-up programs where IceCube scrutinizes its own data to search for correlations with signals from Gravitational Waves \cite{GW}.

\section{Where do we go?}

Four years after the detection of cosmic neutrinos, we have learned a lot about their
spectrum and flavor composition.
We have learned that blazar jets and GRBs can contribute only a small fraction to
the observed astrophysical neutrino flux. The spectral features of this flux
(single power law or two power law) open new questions about the contributing source classes.
No individual sources have been detected yet. The non-observation of neutrinos coinciding with GRBs strongly constrains models which attribute the highest-energy cosmic rays to GRBs. Neutrino events possibly related to supernova explosions have been observed, although with a non-negligible probability for a chance occurrence.
No neutrinos have been observed that could be attributed to the GZK effect \cite{BZ}, but the non-observation starts constraining evolution scenarios
for ultra-high energy cosmic rays sources (not addressed in this report).

IceCube continues collecting data. A twofold statistics combined with improved directional precision,
also for cascade-like events,
and better understanding of systematics effects will considerably improve the understanding of what has
been observed so far and may even provide first detection of individual (point-like or extended) sources.
IceCube's capabilities, however, are limited by its size, the chance to detect neutrinos from the central part of our Galaxy are constrained by its location at the South Pole.

The next important steps are being done at the Northern hemisphere:  GVD in Lake Baikal and KM3NeT in the Mediterranean Sea.

\newpage

\subsection{Baikal-GVD}

Baikal-GVD (Gigaton Volume Detector) is configured in ``clusters'', where each cluster consists of eight strings, instrumented over a length of 520\,m with 36 optical modules. The OMs house 10'' Hamamatsu photomultipliers with a high-sensitive photocathode. After an extensive period of in-situ tests of single components and prototype strings, a first cluster with 24 OMs per string was deployed in Spring 2015. One year later, it was upgraded to a full cluster with 36 PMs per string, and in Spring 2017 a second cluster was added.
First preliminary results of the 2016 cluster have been presented at the ICRC 2017 \cite{Baikal-ICRC}.

Baikal-GVD will be built in two phases. Phase-1 will consist of eight clusters, each 120\,m in diameter,
with lateral distances of 300\,m (see Fig.\,\ref{GVD}). The effective volume for cascades in the 10-100 TeV energy range will be about 0.4\,km$^3$, the sensitivity to muons is negligible below 1\,TeV but rapidly raises in the multi-TeV range. Phase-1 is financed and planned to be completed in 2020/21. In a second phase, Baikal-GVD will be extended to 18 clusters and then surpass the cubic kilometer benchmark.

\begin{figure}[h]
\centering
\includegraphics[width=0.8\textwidth]{figures/GVD.pdf}
  \caption{\small
Schematic view of phase-1 of Baikal-GVD, consisting of 8 clusters, each with 120\,m diameter and 520\,m height.
A cluster consists of eight strings with 36 optical modules along each string.
	}
\label{GVD}
\end{figure}

\subsection{KM3NeT}

KM3NeT has two main, independent objectives: a) the discovery and subsequent observation of
high-energy cosmic neutrino sources and b) precise oscillation measurements and the determination of the mass hierarchy of neutrinos. For these purposes
the KM3NeT Collaboration plans to build an infrastructure distributed over three sites:
off-shore Toulon (France), Capo Passero (Sicily, Italy) and Pylos (Peloponnese, Greece).
In a configuration to be realized until 2020/22, KM3NeT will consist of three so-called building blocks
(``KM3NeT Phase-2'').
A building block comprises 115 strings, each string with 18 optical modules.Two
building blocks will be sparsely configured to fully explore the IceCube signal with a comparable instrumented
volume, different methodology, improved resolution and complementary field of view, including the Galactic
plane. These two blocks will be deployed at the Capo Passero site and are referred to as ARCA: Astroparticle Research with Cosmics in the Abyss. The third building block will be densely configured to precisely measure atmospheric neutrino oscillations.
This block, being deployed at the Toulon site, is referred to as ORCA: Oscillation Research with Cosmics in the Abyss (see Fig.\,\ref{KM3NeT}).

\begin{figure}[h]
\centering
\includegraphics[width=0.65\textwidth]{figures/KM3NeT.pdf}
  \caption{\small
The two incarnations of KM3NeT. The two ARCA blocks (bottom) have diameters of 1\,km and a height of about 600\,m and focus to high-energy neutrino astronomy. ORCA (top) is a shrinked version of ARCA with only 200\,m diameter and 100\,m height. Both ARCA and ORCA have 115 strings with 18 optical modules (OMs) per string. Top left, a drawing of an
OMs is shown. Each OM houses 31 small photomultipliers.
	}
\label{KM3NeT}
\end{figure}

A novel concept has been chosen for the KM3NeT optical module:
The 43\,cm glass spheres of the DOMs will be equipped with 31 PMTs of 7.5\,cm
diameter, with the
following advantages: a) The overall photocathode area exceeds that of a 25\,cm
PMT by more than a factor three; b) The individual readout of the PMTs results
in a very good separation between one- and two-photoelectron signals which is
essential for online data filtering; c) some directional information is
provided. This technical design has been validated with in situ prototypes. A
cross-sectional view of this DOM is shown at the top of Fig.\,\ref{KM3NeT}.

\vspace{0.3cm}
With a fully equipped ARCA, IceCube's cosmic neutrino flux could be detected
with high-significance within one year of operation. In practise the detector will
be deployed in stages allowing to reach the one-year sensitivity of two clusters much before
the second cluster is fully installed. Actually the same is true for Baikal-GVD, with a good chance that
GVD will find cosmic neutrinos before ARCA. ORCA could determine the neutrino mass hierarchy with at least 3\,$\sigma$ significance after three years of operation, i.e. as early as 2023.

\subsection{IceCube-Gen2}

The progress from IceCube over the next decade is limited by the modest
numbers of cosmic neutrinos measured, even in a cubic
kilometer array. In \cite{Gen2} a vision
for the next-generation IceCube neutrino observatory is presented.
At its heart is an expanded array of optical modules
with a volume of 7 to 10\,km$^3$.  This high-energy array will mainly address
the 100\,TeV to 100\,PeV scale. It has
the potential to deliver first
GZK neutrinos, of anti-electron neutrinos produced via
the Glashow resonance, and of PeV tau neutrinos,
where both particle showers associated with the production
and decay of the tau are observed (``double bang events'').

Another possible component of IceCube-Gen2 is the
PINGU sub-array. It targets -- similar to ORCA -- precision measurements
of the atmospheric oscillation parameters and the
determination of the neutrino mass hierarchy. The facility's
reach would further be enhanced by exploiting the
air-shower measurement and vetoing capabilities of an
extended surface array. Moreover, a radio array
(``ARA'', for Askarian Radio Array) will achieve improved
sensitivity to neutrinos in the $10^{16} - 10^{20}$ eV energy
range, including GZK neutrinos. Figure \ref{Gen2} sketches a possible design of IceCube-Gen2.

\newpage

\begin{figure}[h]
\centering
\includegraphics[width=0.8\textwidth]{figures/IceCube_Gen2.pdf}
  \caption{\small
Schematic view of IceCube Gen-2, comprising the existing IceCube array with
its densely equipped inner region DeepCore, the high-energy array of Gen2, the super-densely
equipped PINGU sub-detector, and an extended surface array. Not shown is the radio array ARA with
its size exceeding that of the basic surface array.
	}
\label{Gen2}
\end{figure}

For point sources, the high-energy array will have five times better
sensitivity than IceCube, and the rate for events at energies above a few hundred TeV will be ten times higher than for IceCube.

\subsection{Expectations}

Personal expectations do not necessarily agree with the optimistic time schedules of the experiments.
Still, one can certainly expect that at the early 2020s, the IceCube cosmic neutrino signal will be scrutinized by KM3NeT-ARCA and Baikal-GVD, with different experimental systematics and from the Northern hemisphere: what is the exact form of the spectrum, what is the flavor composition, what is the contribution from our Galaxy? The diffuse flux from the Galactic plane will be almost certainly discovered, at least if moderately conservative predictions are correct. The discovery of point or extended sources in our Galaxy is not guaranteed but seems likely. Some bet that the first discovery of individual sources will come from transient events correlated to electromagnetic or gravitational wave
observations. Actually, within the next few years, this seems to be the most promising chance for IceCube.

Another aspect of increased discovery potential are analyses combining data from different detectors. Within the Global Neutrino Network GNN \cite{GNN}, such analyses are being performed between IceCube and ANTARES \cite{common} and bear the chance to detect diffuse emission from the Galaxy rather soon. When Baikal-GVD and ARCA approach the cubic kilometer scale, these efforts will become even more important than now.

With the appearance of four or more ARCA blocks, with GVD Phase-2 and with IceCube Gen2 in the second
half or at the end of the 2020s,
one could have 3--5 km$^3$ instrumented volume in the North and 7--10 km$^3$ in the South.
I seems unlikely that this coordinated attack to the high-energy neutrino frontier would
fail to detect structures and individual sources. That, at the end, will allow charting the
neutrino landscape to which IceCube has enabled a first glance!

High Energy Astrophysical Phenomena (astro-ph.HE)


The Spatial Structure of the Evershed effect

Abstract: The line profile asymmetry defining the Evershed effect in sunspot penumbrae, disappears abruptly at the outer sunspot boundary over a horizontal distance of less than 500 km immediately at the outermost ends of the dark penumbral continuum structures ('filaments'). This is indicated from the non-magnetic lines Ni I 5435.9 (g=0.5) and Fe I 5434.5 (g=0), which loose their profile asymmetries at the same spatial location although formed at a vertical height distance of 300 km. The widely accepted 'canopy' picture conflicts with this finding which, in contrary, suggests that the gas motion associated with the Evershed effect disappears together with the continuum intensity drop at the outer sunspot (penumbral) border. Corresponding downflows must then be smaller than the spatial resolution presently achieved (i.e., <0.5 arcsec) and located at the immediate outer sunspot border. It is suggested that this location marks the sharp threshold of the equipartition between kinetic and magnetic energy density at the outer penumbral boundary.

\section{Introduction}

The characteristic signature of the Evershed effect is a line profile
asymmetry which often occurs as a pronounced line 'kink' (Maltby 1964)
and can not be explained by a velocity gradient with depth (Stellmacher
\& Wiehr 1980). Instead, one has to consider that the observed line profiles
represent superpositions of spatially unresolved penumbral structures: a
stronger 'main component' originating from the bright structures superposed
by a Doppler-shifted and weaker line 'satellite' from the dark structures
(Stellmacher \& Wiehr 1971). This is supported by the finding that pronounced
line 'kinks' only occur if the (dark) penumbral structures are oriented almost
along the line-of-sight, i.e. in center-side penumbrae of spots very close
to the limb (Wiehr 1995).

Observed differences between the Evershed effect
in center- and limb-side penumbrae suggest angles of inclination to the solar
surface between $6^o$ (Maltby 1964) and $15^o$ (Schr\"oter 1965). Such flat
structures cannot gradually 'disappear with height' as is to be expected for
a penumbral 'canopy' (cf. Solanki et al. 1992). In that picture the profile
asymmetry of Fe\,I\,5434\AA{} (formed about 500 km above the continuum level),
should 'disappear' over a horizontal distance of several arcsec beyond the
white-light penumbral border. This, however, disagrees with the results of
Wiehr \& Degenhardt (1992,1994).

On the other hand, the chromospheric H$\alpha$ superpenumbra extends far
beyond the edge of the continuum penumbra. But the opposite sign of the
chromospheric Evershed effect (e.g. Maltby 1975) indicates that it cannot
represent a continuation of the photospheric Evershed effect at higher
levels. Hence, a possible extension of structures, flow, and magnetic
field beyond the white-light penumbra needs further study. The present
paper investigates an upper limit for the abrupt disappearance and its
physical reason.

\section{Observations}

On October 22, 1994, the neighboring lines Fe\,I\,5434.5 (g=0) and
Ni\,I\,5435.9 (g=0.5) were observed with the Gregory telescope on Tenerife
(Wiehr 1986) in a sunspot at a heliocentric angle of $\vartheta = 65^o$
with a 1024 x 1024 pixel CCD at an exposure time of 0.1\,seconds. In order
to avoid the use of an image rotator, the time of observation was chosen
such that the Coud\'e rotation of the solar image yields an almost parallel
orientation of the slit to the limb nearest the spot. During the two hours
of observation the image rotation increased the angle between slit and limb
from $-15^o$ to $+15^o$ with respect to the exactly parallel orientation.

An accurate pointing of the 40\,$\mu$m slit at very small distances from
the outer penumbral edge (e.g. $\Delta<0.5$ arcsec) is difficult, since the
slit occurs as a 1/3 arcsec wide dark streak on the slit-jaw image just at
the location under study. However,  if the telescope's guiding automat
is switched off, a slow image drift of about 0.5 arcsec/min occurs, which
affects a kind of automatic scanning if a large number of CCD-spectra is
taken at a short time cadence. This yields spectra at a variety of slit
distances including such of few 100\,km inside or outside the penumbral
edge. Since the image rotation slowly varies the slit orientation on the
sun, a repetition of such a series of 'drifted spectra' gives slit
orientations of different azimuth angles with the spot border, thus
covering large parts of the limb- and of the center-side penumbral
boundaries. In the best spectra, a spatial resolution near 0.5\,arcsec
is achieved. Two examples are shown in figure\,1.

   \begin{figure}[ht]
   \hspace{-4mm}\includegraphics[width=9.9cm]{EvershedFig1.jpg}
   \caption{CCD spectra of Fe\,I\,5434.5 (left) and Ni\,I\,5435.9 (right)
a) immediately inside the center-side, (b) just outside the limb-side
penumbral border; the spatial extension of each spectrum covers 50\,arcsec.
The black dot marks a profile asymmetry {\it not} related to the Evershed
effect.}
   \label{Fig1}
    \end{figure}

   \begin{figure}[htb]
   \hspace{0mm}\includegraphics[width=9.2cm]{EvershedFig2.jpg}
   \caption{Fe\,I\,5434.5 and Ni\,I\,5435.9 with profile asymmetry
at a location of a continuum depression produced by the end of a dark
penumbral structure (solid line in the upper panel) in comparison with a
location immediately outside the penumbra (dotted line). Lower panel:
profile asymmetry only for Fe\,I\,5434.5 and not related with a continuum
depression (solid line), originating from a high velocity region (dot in
Fig.\,1b) outside the spot.}
   \label{Fig2}
    \end{figure}

\section{Results}

The numerous spectra obtained from slit positions in the immediate
vicinity of the penumbral border confirm (e.g. Wiehr \& Degenhardt 1994)
that line asymmetries occur exclusively at locations of pronounced
continuum depressions produced by dark penumbral structures (solid line
in Fig.\,2a). In the neighboring bright structures, both line profiles
are symmetric (dotted lines in Fig.\,2). {\it Immediately outside the
penumbra, where the slit just 'missed' the outer ends of the dark
structures, the line profile asymmetries disappear simultaneously for
both lines.}

The spatial distance between adjacent positions of the slit is partly
determined by slow image drift and mainly by the image motion. The latter
can be expected to be slightly smaller than the 0.5\,arcsec spatial
resolution in the spectra, which are additionally degraded by the
spectrograph optics. Assuming a spatial distance of the order of the
slit width of 1/3\,arcsec and correcting this value for the geometric
foreshortening at $\vartheta=65^o$, one obtains a {\it horizontal distance
of maximally 500\,km over which the line profile asymmetries disappear
simultaneously for the two lines, which formed at a vertical distance
of about 300\,km}.

Among the numerous spectra showing a strictly simultaneous disappearance
of the asymmetry of both lines, only one location is found (dot in Fig.\,1b)
which exhibits a line asymmetry of Fe\,I\,5434 but not of Ni\,I\,5436
(see Fig.\,2b). However, this structure shows no continuum depression
and is thus not a dark penumbral filament. It is located between a bright
continuum streak from a facula and another streak with opposite line
asymmetry (cf. Fig.\,1b) thus indicating a disturbed region at the penumbral
border (possibly an Ellerman bomb).

\section{Discussion}

The simultaneous disappearance of the profile asymmetry of both lines
together with the continuum depressions supports the idea that the profile
asymmetries are due to superposed line satellites from spatially unresolved
dark penumbral structures (Wiehr 1995). These line 'satellites' disappear
in photospheric layers together with the dark structures where they are
formed. If a 'canopy' existed, one would expect the asymmetry of Fe\,I\,5434,
(a line formed near the temperature minimum; see e.g. Degenhardt \& Wiehr
1994)) to disappear further away from the penumbral continuum border than
that of Ni\,I\,5436 (formed about 300\,km higher). In the 'canopy' scenario,
the simultaneous disappearance over a horizontal distance of less than 500\,km
would require a steep inclination of the flow (and thus of the magnetic
field) of at least $60^o$, which severely conflicts with the values observed
by Maltby (1964) and by Schr\"oter (1965).

The present observations, however, strongly suggest a nearly horizontal
orientation of the Evershed flow (and thus also of the magnetic field in the
dark penumbral structures) with an abrupt disappearance at the very border
of the white-light penumbra. A plausible explanation for such a sharp end is
the equipartition between the magnetic and kinetic energy densities: For a
photospheric density value of $\rho = 5\cdot 10^{-7}$ g/cm$^2$, a granular
velocity of 3\,km/s balances 740\,Gs, which is close to values observed at
the penumbral edge. If the magnetic field varies radially with $B\sim r^{-2}$
(cf. Wittmann 1974), the magnetic energy density $B^2/8\pi$ decreases with
$r^{-4}$. This results in a sharp threshold for an equipartition with the
kinetic energy density of the neighboring granules. Flow and field will
than (vertically) sink at the very penumbral border.[Such a 'disappearance'
of the flow may conflict with the siphon model (Meyer \& Schmidt 1968) but
hardly with the idea of 'turbulent pumping' (Parker 1974).]

A 'canopy' formed by such an extension of the bright penumbral structures
does not contradict R\"uedi et a1. (1995) who deduced significantly
steeper inclinations from the He\,I\,line at 10830\AA{} than expected
from a neighboring Si\,I\,line. The latter, however, is formed in deep
layers which are dominated by flat dark structures. The helium line,
instead, originates from high layers, which are dominated by the steep
bright structures.

\begin{acknowledgements}
The author acknowledges numerous discussions with C.\,R.
de\,Boer, G.\,Stellmacher, H.\,Balthasar, M.\,Sch\"ussler and a critical
review of the manuscript by C.\,Zwaan.
\end{acknowledgements}

\eject

Solar and Stellar Astrophysics (astro-ph.SR)


Does Every Quasar Harbor A Blazar?

Abstract: Assuming there is a blazar type continuum in every radio-loud quasar, we find that the free-free heating due to the beamed infrared continuum can greatly enhance collisionally excited lines, and thus explain the stronger CIV $\lambda$1549 line emission observed in radio loud quasars. We further predict that the CIV line should show variability {\it not} associated with observed continuum or Ly$\alpha$ variability.

\section{Introduction}

The dichotomy between radio-loud and radio-quiet QSOs
remains one of the greatest
puzzles in studies of AGN. While their optical-UV spectra appear
similar,
quasars (we use the term ``quasar'' for radio loud QSO)
have CIV $\lambda$1549 line emission $\sim$50\% stronger
compared with radio-quiet QSOs (Francis, Hooper, \& Impey 1993).
Marziani et al. (1996) find that the luminosities of CIV and H$\beta$ are
strongly correlated in radio-quiet QSOs, but only weakly
correlated in quasars. Some observed line profile differences could
also be the result of an additional
beamed optical-UV synchrotron continuum in quasars,
but not in radio-quiet QSOs (Wills et al. 1993; Francis et al. 1993).

Unified schemes for quasars hold that
there is a beamed synchrotron continuum in
every quasar (e.g., Urry \& Padovani 1995).
If the beam is toward the observer, the quasar is
seen as a blazar. Blazars are known to be highly variable on time
scales from days to years. For example, BL Lacertae
brightened by 3 magnitudes in R band during the
outburst in the summer of 1997 (Bloom et al. 1997).
The outburst lasted for a few months,
and during the outburst it shows intraday and even intrahour variation
of $\sim$1 magnitude.

Bregman et al. (1986) found that the Ly$\alpha$ line flux varied in
proportion to the continuum flux density in the blazar 3C446.
Stephens \& Miller (1984) found that the CIII] $\lambda$1909 line flux varied
with the continuum in 3C446 and suggested that
the beamed continuum might affect the emission line spectra.
Koratkar et al. (1998) found no variation in the flux of Ly$\alpha$ emission line
in HST spectra of the blazar 3C 279 while the
continuum varied by a factor
of 50. Cohen et al. (1987) found possible Mg II line variability in
the blazar AO 0235$+$164.
Corbett et al. (1996) found the H$\alpha$ line in BL Lac to be
variable, but not as a result of the beamed continuum.
The emission line variability in blazars is not
easy to measure because of the large contribution from the
variable continuum.
Most previous studies of AGN spectral variability
have focused on objects known to have large continuum variability
and concentrated on the emission line variations following the
{\it observed} continuum variations (e.g., Maoz et al. 1994).

In this Letter we investigate the interactions between  the
beamed  synchrotron radiation and the BELR. Although the beam may illuminate
only a small portion of the Broad Emission Line Region (BELR),
it can provide a large flux of infrared
photons that heat the gas by the
free-free mechanism (Ferland et al. 1992), greatly enhancing the CIV
$\lambda$1549 line emission. This enhanced emission can significantly
contribute to the
line flux from the entire BELR, explaining the stronger line emission
seen in quasars.
The CIV line emission
is thus expected to vary in response to the rapid and
large amplitude variations of the
blazar continuum, although we may not see the continuum variation if
it is beamed away from us.

\section{Additional Beamed Continua in Quasars}

Theoretically, a simple blazar synchrotron continuum is
beamed with a half opening angle
${\theta}={\sin}^{-1}(1/\gamma)$
for a Lorentz factor $\gamma$.
The fraction of the sky area illuminated by the beam is $1/2{\gamma}^2$,
or ${\sim}2\%$ for ${\gamma}=5$, and ${\sim}5\%$ for ${\gamma}=3$.
The real beaming angle could be larger (Lind \& Blandford 1985),
and can be estimated from the observed fraction of blazars among quasars
as a whole (e.g., Urry \& Padovani 1995).

If the beamed synchrotron continuum originates closer to
the central black hole than
the BELR (${\sim}0.1$ pc in luminous quasars,
see e.g., Netzer 1990) then it may affect
the emission line spectrum. In NGC\,1275 the VLBI jet
probably originates within 0.1 pc at $43$ GHz (Krichbaum et al. 1992).
The higher frequency synchrotron
radiation of interest (${\sim}10^{12}$ Hz) is likely to
originate closer to the black hole.
Some theoretical models predict infrared synchrotron emission
originating ${\sim}0.1$ pc from the black hole
and that lower frequency emission originates farther away
from the center (see e.g., Marscher 1996).

The beamed synchrotron continuum provides a high flux
of IR photons because of
its steep spectrum, which extends to mm wavelengths.
In contrast, the normal QSO
continuum provides only weak IR flux, with much of
the observed IR emission
almost certainly arising from dust beyond the
BELR (Barvainis, Antonucci, \&
Coleman 1992).  Moreover, if the observed IR continuum
in radio-quiet QSOs were
within BELR distances it would significantly overpredict
the observed CIV strength.

\section{Models}

Ferland et al. (1992) showed that a high flux of IR photons
incident on BELR clouds could enhance the flux of the CIV and OVI
emission lines by free-free
heating of the BELR.
As discussed above, a strong IR central continuum is
unlikely in radio-quiet
QSOs.  However, for BELR gas exposed exposed to the
beamed synchrotron radiation,
free-free heating is the most important
heating process.

The heating rate is proportional to (e.g., Netzer 1990)
\begin{equation}
 \displaystyle
  F_{\nu}e^{-\tau_{\nu}(ff)}T^{-1/2}_eN_eN_i{\nu}^{-3}(1-e^{-h\nu/kT_e}).
\end{equation}
Collisional excitation results in an important emission line cooling
process in optically thick BELR gas, and
the rate is proportional to
\begin{equation}
 \displaystyle
  N_e\Omega_{ij}e^{-E_{ij}/kT_e}T^{-1/2}_e,
\end{equation}
where $\Omega_{ij}$ is the collision strength. For CIV $\lambda$1549,
$E_{ij}\sim 8$eV, and the cooling rate increases with temperature.
At ${\sim}16$ eV, this line intensity stops
increasing with $T_e$ because the collision cross sections become smaller
for higher velocity electrons.
Roughly speaking, the IR heating leads to a higher electron temperature,
which results in stronger CIV $\lambda$1549.

We use a fixed $L_{\rm agn} = 10^{46}$ erg s$^{-1}$  for an
isotropic AGN continuum and use the spectral energy distribution (SED)
of Mathews \& Ferland (1987) except that the submillimeter break
is at 10\micron.  The beamed continuum is conservatively assumed to have
 $f_{\nu}~{\propto}~{\nu}^{-1}$ from
$10^{-3.5}$ to $10^{-2}$ Ryd, and
$f_{\nu}~{\propto}~{\nu}^{-2}$ from $10^{-2}$ to $10^6$ Ryd.
At lower energies, $f_{\nu}~{\propto}~{\nu}^{2.5}$.
We assume that  the beamed continuum has constant
intensity within the small opening angle $2\theta$, and
use $L_{\rm beam}$ to represent the apparent luminosity as seen by
an observer within the beam.
Adding $L_{\rm beam}$ essentially
does not change the ionization parameter because of the steepness of
the spectrum.

Solar abundances are assumed for the BELR gas. The calculations
are done using CLOUDY (Ferland 1996).
We first try a simple BELR model ``A'' resembling a model given by
Ferland \& Persson (1989).
They use a fixed hydrogen density $n_H=10^{9.5}$ cm$^{-3}$,
and a fixed ionization
parameter $U = 10^{-0.5}$. Because we are using
$L_{\rm agn}=10^{46}$ erg s$^{-1}$
for the input continuum,  $U = 10^{-0.5}$ constrains the cloud distance to
$r=10^{17.67}$cm. The column density for the cloud is set to be
$10^{25.5}$ cm$^{-2}$.

In the more realistic model ``B'' we follow the model of Goad, O'Brien, \&
Gondhalekar (1993),
in which discrete clouds are distributed from $10^{17}$ cm to $10^{19}$ cm
from the continuum source. The density starts from $10^{13}$ cm$^{-3}$,
and has a
power law distribution $n~{\propto}~r^{-2}$. The column density is
set to be $N(r)~{\propto}~r^{-4/3}$, starting from $10^{23}$ cm$^{-2}$ at
the inner radius. The differential covering factor varies as
${\propto}r^{-1/6}$,
to fulfill the mass conservation
of the clouds. The calculated emission line luminosities are integrated
from inner to outer radius.

\section{Results}

In Tables 1 and 2 we list the calculated emission line luminosities for
models A and B.
For the H$\beta$ line we give the logarithm of luminosity in units of
erg s$^{-1}$, and the other line
luminosities are given as the ratio to H$\beta$.

Column 2 is the result for  $L_{\rm agn}=10^{46}$ erg s$^{-1}$ with no
$L_{\rm beam}$;
column 3 is for $L_{\rm agn}=10^{46}$ erg s$^{-1}$ mixed with
$L_{\rm beam}=10^{46}$ erg s$^{-1}$, and  $L_{\rm beam}$
affects only 2\% of the sky area. A small covering factor (${\sim}0.1$) for the
clouds is assumed everywhere in the sky, including in the beam.
Columns 4 and 5 correspond to
$L_{\rm beam}=10^{46.5}$ and $10^{47}$ erg s$^{-1}$.
The observed line ratios in column 6 are taken from Netzer (1990).

From Tables 1 and 2 it can be seen  that among strong lines,  the
CIV $\lambda$1549 is enhanced by ${\sim}$50\% when
$L_{\rm beam}=10^{47}$ erg s$^{-1}$,
while hydrogen lines show little change.
The CIII $\lambda$977 is also strongly dependent on the
IR heating. These two lines actually set up an upper limit
of $L_{\rm beam}~{\sim}~10^{47}$
erg s$^{-1}$, otherwise they would be too strong. However, we note
that CIII $\lambda$977 is difficult to observe because of
intervening absorption, especially at higher redshifts.
Model A predicts a slight enhancement of OVI $\lambda$1035 and
NV $\lambda$1240. We do not consider their enhancement in model
B, because the strengths of these two lines are not well predicted.

In Fig. 1 we plot for model A the abundance ratio C$^{3+}$/C,
and the electron temperature versus the depth into the cloud.
As expected, the C$^{3+}$ density is not strongly affected by the beamed
continuum, and CIV $\lambda$1549 is an effective coolant.

In Fig. 2, we give cloud surface temperature and
various line luminosities versus radius in model B.
The free-free heating is more important at small radii, since it is
proportional to the photon flux as shown in eq. (1), while photoelectric
heating by ionizing photons is determined by the hydrogen number
density and the recombination coefficient and does not depend
explicitly on photon flux (Ferland et al. 1992).

Some line profiles with and without a beamed continuum are compared
in Fig. 3. We show the profiles here primarily for the purpose of showing
the enhancement of integrated emission lines flux due to the
beamed continuum, though affecting only a small portion of the
BELR clouds. The BELR cloud distribution and velocity field
are probably anisotropic with velocities having a strong
disk component in addition to a random component (e.g., Wills \& Browne
1986). For simplicity,
here we consider only spherically distributed clouds with a
Keplerian velocity field, and the velocity direction is assumed to be
random. The maximum cloud velocity is $10^4$ km s$^{-1}$ at $r=10^{17}$ cm.
The flat tops of the profiles come from the lowest gas velocity,
$10^3$ km s$^{-1}$ at $r=10^{19}$ cm.

\section{Discussion}

The IR-optical continuum of the
synchrotron radiation is steep, hence
the fraction of hydrogen ionizing photons it can
offer is small compared with the normal AGN continuum.
In addition, the beam illuminates only a small fraction of
the BELR clouds, while the emission lines we see are integrated
over all the clouds. Consequently, the
effect on  the H$\alpha$ and Ly$\alpha$ lines is negligible.
For collisionally excited lines such as CIV $\lambda$1549,
the line emissivity can be enhanced by over two orders
of magnitude in the beam, and the line intensity integrated
over the whole BELR consequently increases by a factor of $\sim$50\%.

If there is a beamed blazar type continuum in every radio loud
quasar, we should be able to see CIV line variability
{\it not} associated with observed continuum
variability or Ly$\alpha$ variability. This phenomenon
could well have escaped
detection in previous observations because high redshift
quasars have not been extensively monitored. Also, this phenomenon
may be seen only when the blazar continuum is in the outburst stage.
We are planning observations to search for CIV $\lambda$1549
line variability among $z~{\sim}~2.5$ quasars,
which are thought to be less variable compared with low luminosity
AGNs. If the CIV line variability at the level $\sim$50\% is observed,
while the lines for comparison (Ly$\alpha$, CIII]
$\lambda$1909 or MgII $\lambda$2798)
do not change significantly, it
will be  strong support for our models and for quasar unified schemes.
In addition, the disk-wind
model (e.g., Emmering, Blandford, \& Shlosman 1992) for the BELR
will be ruled out, since it does not include any
gas in the polar regions.

We also expect quasars to have statistically stronger CIV $\lambda$1549.
The strength of CIV $\lambda$1549
depends strongly on the blazar continuum, and should
not be correlated well with H$\beta$ line strength as shown by
Marziani et al. (1996).
The intensity variation
may be seen in the line core or line wings, and thus offer
information about the gas kinematics in the polar regions, e.g., whether
the gas clouds have lower velocity or have large radial velocities
driven by the jet.

In reality, the line emission depends on the details of BELR
geometry.  If the line is strongly beamed backward
(e.g., Ferland et al. 1992), we may not be able to see
the enhanced CIV emission when the angle between the line of sight
and the beam is large. However,
the excited lines are at least more isotropic than the beamed
continuum, and we do expect to see the line variability without
seeing the beam.

If the phenomenon suggested by us is not detected, this
will still give strong constraints on models. Three possibilities are left.
First, the IR synchrotron radiation originates
outside the BELR. Second, the observed blazar continuum variability
may be caused by a change of the jet direction rather than a change
in strength of the beam (e.g., Abraham \& Carrara 1998).
Third, there is no
broad line emitting gas in the polar regions. The last possibility
would favor the disk-wind model and may rule out the bloated star model
(e.g., Alexander \& Netzer 1994).

Finally we note that if the BELR is made of small clouds
confined by a hot medium (Krolik, McKee, \& Tarter 1981), then the jet
might destroy the clouds. However, if the BELR is made of
bloated stars,
the jet may enhance the mass loss rate of the stars and increase the
amount of broad line emitting gas in the beam.

\acknowledgements

We thank Greg Shields for helpful discussions, and Gary Ferland
for advice and making his code CLOUDY available. We also thank the referee,
Ski Antonucci, for valuable comments.
This work is supported by NASA LTSA grant NAG5--3431.

Does Every Quasar Harbor A Blazar?


A Generic Information Extraction System for String Constraints

Abstract: String constraint solving, and the underlying theory of word equations, are highly interesting research topics both for practitioners and theoreticians working in the wide area of satisfiability modulo theories. As string constraint solving algorithms, a.k.a. string solvers, gained a more prominent role in the formal analysis of string-heavy programs, especially in connection to symbolic code execution and security protocol verification, we can witness an ever-growing number of benchmarks collecting string solving instances from real-world applications as well as an ever-growing need for more efficient and reliable solvers, especially for the aforementioned real-world instances. Thus, it seems that the string solving area (and the developers, theoreticians, and end-users active in it) could greatly benefit from a better understanding and processing of the existing string solving benchmarks. In this context, we propose SMTQUERY: an SMT-LIB benchmark analysis tool for string constraints. SMTQUERY is implemented in Python 3, and offers a collection of analysis and information extraction tools for a comprehensive data base of string benchmarks (presented in SMT-LIB format), based on an SQL-centred language called QLANG.

\subsection{Further Examples}

Next, we show some examples of benchmark analysis, realizable by our~tool.\looseness=-1

\smallskip
\emph{\S\ A variable dependency analysis.} To speed up the solving process for a particular string constraint, one might be interested in splitting a formula into multiple, independent sub-formulae. A relatively na\"ive way of splitting a formula is to determine whether the parts of the input formula in which each variable occurs, and see how they overlap. We can use \textsc{SMTQuery}{} to visualize the interactions between variables within a formula, using the AST data structure. Whenever an AST is built starting from an SMT-LIB file, we store the variable occurrences for each node, which is actually the only information we need to build a simple variable-dependency graph. So, we have to define the generation of the corresponding plot by implementing the aforementioned class interface in \texttt{smtquery.extract}, call it \texttt{VarDepPlot}, and register it, after implementing the required logic in the \texttt{PullExtractor}. Now we can simply ask a query, e.g. \qline{Extract VarDepPlot From * Where hasWEQ}, creating a variable-dependency plot for each registered instance which contains at least a single word equation.

\begin{figure}[t]
  \begin{center}
    \begin{tikzpicture}[->,>=stealth',shorten >=1pt,auto,node distance=0.75cm,semithick]
      \node[]   (0)  at (0,0)               {\includegraphics[width=\textwidth]{vardep}};
    \end{tikzpicture}
  \end{center}
  \vspace*{-1cm}
\caption{Example variable dependency plot.}
\vspace*{-0.5cm}
\label{fig:varplot}
\end{figure}

Posing the query, e.g., have the instance of Listing~\ref{lst:smtexxample2} leading to the plot of Figure~\ref{fig:varplot} where variables $G,H$ and $J$ are top nodes and the corresponding assertions are to bottom nodes.

\begin{lstlisting}[caption={Instance for variable dependency plot},label=lst:smtexxample2,captionpos=b]
(set-logic QF_S)
(declare-fun H () String)
(declare-fun G () String)
(declare-fun J () String)
(assert (= (str.++  "aacfb" G "abdeddaaa")  (str.++  "aacfbdffebaaaaac" H "aaa") ))
(assert (= (str.++  "efbaeecedaaecfceffaffaedfcebcf" J "aeaadcbe")  (str.++  "e" J "aeecedaaecfceffaffaedfcebcf" J "aeaadcbe") ))
(assert (= (str.++  "bdceafbededddcfcacffdeaefcfa" J "dbabcdebee")  (str.++  "bdceafbededddcfcacffdeaefcfa" J "dbabcdebee") ))
(check-sat)
\end{lstlisting}

The edges in the figure indicate the presence of a variable, allowing us to split the instance accordingly.

\emph{\S\ Analyzing the performance of a string solver.} To review the performance of a particular solver, one is usually interested in getting a comparison w.r.t. other solvers. \textsc{SMTQuery}{} allows exporting a summary table and the commonly used cactus plot to compare string solvers on benchmarks. For instance, we want to see whether \textsc{cvc5}{} is performing well enough on the \textsc{Woorpje}{} benchmark set. First, we need to trigger \textsc{cvc5}{} on the respective benchmark by executing, e.g., \qline{Select Name From woorpje where isSAT(CVC5)}. Then, we can obtain a summary table by posing the query \qline{Extract ResultsTable From woorpje} and a cactus plot by simply changing the extractor asking the query \qline{Extract CactusPlot From woorpje}. The results are visualized in Figure~\ref{fig:cactusS}.

\begin{figure}[t]
  \begin{center}
    \begin{tikzpicture}[->,>=stealth',shorten >=1pt,auto,node distance=6cm,semithick]
      \node[]   (0)  at (0,0)               {\includegraphics[width=6cm]{cactus_small}};
      \node[right of=0,text width=5cm]   (0) {\begin{lstlisting}[label=lst:cactus,captionpos=b,basicstyle=\ttfamily\tiny,numbers=none]
                    Z3Str3    Z3Seq     CVC4
----------------  --------  -------  -------
SAT               512       604      594
UNSAT             170       164      164
Unknown            13         0        0
Timeout           114        41       51
Crash               0         0        0
Time w/o Timeout   97.5154  123.816  134.242
Total Time        669.859   329.158  389.631
\end{lstlisting}};
    \end{tikzpicture}
  \end{center}
  \vspace*{-1cm}
\caption{Cactus plot and terminal results for an example query.}
\label{fig:cactusS}
\vspace*{-0.5cm}
\end{figure}

\emph{\S\ Modifying Instances.} Analyzing the real-world benchmarks w.r.t. to a particular type of constraints can be achieved by simply neglecting all others. For example, to analyze the structure of the occurring word equations, one may simply pose the query \qline{Extract SMTLib From * Where hasWEQ Apply Restrict2WEQ} to obtain cleaned SMT-LIB files. Naturally, also in this case, we can use any extraction function implemented to acquire the data which we are interested in.

\emph{\S\ Finding and analyzing sub-theories.} In \cite{words21} we have analyzed a large set of benchmarks w.r.t. regular expression membership queries. This kind of queries plays a central role in verifying security policies, by allowing to restrict the set of possible input strings by a regular language \cite{cav21}. The inspection of the respective benchmarks was performed using a sequence of different handcrafted scripts, restricted to a particular use case. \textsc{SMTQuery}{} provides the means to easily extract this data by simply defining predicates analyzing the regular languages (i.e., regular expressions) occurring in the benchmarks. For example, to gather all instances solely containing regular membership constraints asking whether a string without variables or a single variable is a member of a regular expression without complement or intersection is achieved by posing the query \qline{Select Name From * Where isSimpleRegex}. The key difference is that the definition of the particular predicate is much simpler, due to the extendable structure of \textsc{SMTQuery}{}. As such, we can now simply combine the acquired information with newly developed predicates. Since this analysis lead to a well performing algorithm, presented in \cite{cav21}, we are optimistic that our tool can be used to extract such relevant data, ultimately leading to better techniques in the area of solving string constraints.

\emph{\S\ Analyzing the structure of instances.} \textsc{SMTQuery}{} also offers the possibility of a more in-depth analysis of the (syntactic-)structure of the instance. For instance, knowing that all string variables occurring in a formula are subject to constant-length upper bounds allows us to rephrase the problem as a constraint satisfiability problem over finite domains, and ultimately may lead to faster solutions for it. To extract a list of instances having only length-upper-bounded variables, we can pose the query \qline{Select Name From * Where isUpperBounded}. The predicate analyzes the syntax of the constraints and extracts relevant information.
Another interesting aspect, which can potentially lead to a better choice of an algorithm for solving a particular instance, is the analysis of its combinatorial structure. For example, if we know that each variable is occurring at most twice inside a formula, or that each word equation is of the form $\mathsf{x} \doteq \alpha$ where $\mathsf{x}$ is a variable not occurring in $\alpha$ nor anywhere else in the formula, we can use customized solving techniques, and solve the instances more efficiently. Obtaining such information is done by using the predicate \qline{isQuadratic}, resp.  \qline{isPatternMatching}.

\begin{comment}
\emph{\S\ Efficiency.} All operations implemented so far in \textsc{SMTQuery}{} have a reasonable run-time. Clearly, it is inherent that some operations require a rather long execution time: firstly, we analyze a huge set of instances, and, secondly, the analysis we apply might involve complicated predicates, which are provably computationally hard. Our approach to speeding-up this process is to allow the incremental inclusion in the stored-data of the results of various queries, on which one can build efficiently more complex queries. Summing up, our goal was to build a tool allowing information extraction from benchmarks of string constraints, using a state of the art home-computer, without having to go deep in implementation details.\looseness=-1

Let us give some precise examples. For instance, we have recorded runtimes (on a standard PC) of 19807 instances of the Kaluza-set, and obtained the following results: database initialization (21s), obtaining results from \textsc{cvc5}{}, \textsc{Z3str3}, \textsc{Z3seq} (612s), using several apply-functions without caching (198s) and with caching (32s), extracting plots (200s-300s). The runtimes of the apply-functions are similar due to the underlying data-structure, which efficiently stores/provides useful information.

More details and examples regarding runtimes are given on the website of BASC {\url{https://b4sc.github.io/}}.
\end{comment}
\section{Introduction}
\input{new/intro}

\section{Potential Applications of \textsc{SMTQuery}{} and User Base}
\input{new/users}

\section{Architecture of \textsc{SMTQuery}{}}
\input{new/database}

\section{Use cases and examples}
\input{new/examples}

\section{Conclusion and Future Work:}
\input{new/conclusion}

\newpage

\bibliographystyle{splncs04}

Logic in Computer Science (cs.LO); Databases (cs.DB); Formal Languages and Automata Theory (cs.FL)