Pressure-driven phase transitions in bulk HfS$_2$

Abstract: The effect of hydrostatic pressure up to 27 GPa on the Raman scattering (RS) in bulk HfS$_2$ is investigated. There are two transformations of RS spectra, which take place during compression at pressure between 5.7 GPa and 9.8 GPa as well as between 12.8 GPa and 15.2 GPa. Seven vibrational modes can be observed after the transformation, as compared to four modes before the transformation. The observed change suggests structural change in the material of yet unknown nature. The frequencies of the RS modes observed above the transformation change linearly with pressure and corresponding pressure coefficients have been determined. The other transition manifests itself as a change in the RS lineshape. While a series of well-defined RS modes are observed under pressure below the transition, broad spectral bands can be seen at higher pressure. The overall lineshape of the spectra resembles that of disordered materials. The lineshape does not change during decompression, which suggests permanent nature of the high-pressure transition.

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

Layered transition metal dichalcogenides (TMDs) and in particular their few-layer structures have been drawing the attention of researchers for more than a decade now.
While about 60 layered TMDs are recognized, until now, the attention of researchers is focused on Mo- and W-based compounds. Although basic properties of several TMDs in their bulk form are known for years \cite{wilson1969}, there is still room for studies of $e.g.$ their properties under hydrostatic pressure.

In this work we address the effect of hydrostatic pressure on hafnium disulphide ($\text{HfS}_{2}$ ), a member of group IVB TMD.
$\text{HfS}_{2}$  has been recently shown to have very effective electrical response \cite{xu2015ultrasensitive, kanazawa2016transistor}, which justifies a need to uncover basic properties of the material.
The evolution of lattice dynamics of $\text{HfS}_{2}$  with hydrostatic pressure is studied.
We follow RS spectra as a function of pressure to investigate possible structural transformations of the material as a first-order phase transition has been recently reported in Ref.\cite{ibanez2018high}.

\section{Samples and experimental setups \label{sec:Methods}}

\begin{figure}[t]
	\includegraphics[width=.5\textwidth]{figure_0.pdf}
\caption{Raman scattering spectrum of bulk $\text{HfS}_{2}$  at ambient pressure.} \label{fig:0}
\end{figure}

\begin{figure*}[t!]
	\includegraphics[width=\textwidth]{fig1_2_combined.pdf}
\caption{Raman scattering spectra of  bulk $\text{HfS}_{2}$   as a function of hydrostatic pressure for compression (a) and decompression (b) cycle} \label{fig:1}
\end{figure*}

Single crystals of $\text{HfS}_{2}$  have been synthesized in two-zone chemical vapour transport (CVT) furnace.
First, precursor materials were sealed in quartz tubes 25 cm long and of diameter of 8 cm.
Iodine was used as a transport agent.
The reaction and growth temperatures were set to 1050 K and 950 K, respectively.
The growth was performed continuously for 120 hours.
With higher stability of $\text{HfS}_{2}$  against iodine, CVT growth first starts at elevated temperatures.
In the vapour phase HfI and HfI2 can also dominate the direction of reaction.
After 120 hours the reaction was stopped automatically and the reactor was cooled down to room temperature in 5 hours.
High quality $\text{HfS}_{2}$  single crystals of 0.8-1 cm$^2$ size were grown in the low-temperature part of the reactor.
The pattern of powder X-Ray diffraction (XRD) of the investigated material was found to be in perfect agreement with that of the of the octahedral 1T phase of $\text{HfS}_{2}$  \cite{lukovsky1973IR}.
A very minor amount of an unidentified phase was also detected in the crystal by XRD technique \cite{ibanezPrivate}.

The high-pressure RS measurements were conducted at room temperature in a nonmagnetic diamond anvil cell (DAC) made of BeCu alloy.
Crystallographic oil was adopted as the pressure transmitting medium. The fluorescence line of ruby was used to calibrate pressure.
RS spectra were excited using \mbox{$\lambda$=561~nm} (2.21~eV) radiation from a diode-pumped laser. The excitation light was focused by means of a 50x long-working distance objective with a 0.42 numerical aperture (NA) producing a spot of about 1~$\upmu$m diameter. The signal was collected via the same microscope objective, sent through 0.75~m monochromator and then detected by using a liquid nitrogen-cooled charge-coupled device (CCD) camera.
The spectral resolution of the setup in the studied energy range of the RS experiment was equal to 0.7~cm$^{-1}$.

\section{Results \label{sec:Results}}

There are six normal vibrational modes at the $\Gamma$ point of 1T-$\text{HfS}_{2}$  \cite{lukovsky1973IR}: \\
\vspace{-35pt}
\begin{center}
\begin{equation*}
     \Gamma=A_{1g}+E_g+2A_{2u}+2E_u
\end{equation*}
\end{center}

\noindent The A$_{1g}$ and E$_g$ modes are Raman-active, the infrared-active modes A$_{2u}$ and E$_u$ are split into LO and TO branches (for the atom displacement pattern see $e.g.$ \cite{neal2021Raman}), and other A$_{2u}$ and E$_u$ modes are acoustic\cite{lukovsky1973IR} .

The measured RS spectrum of bulk $\text{HfS}_{2}$  at ambient pressure consists of four modes (see Fig.~\ref{fig:0}), the frequencies of which
equal approx. to: 141~cm$^{-1}$, 266~cm$^{-1}$, 325~cm$^{-1}$, and 344~cm$^{-1}$.
This is consistent with previous studies \cite{ibanez2018high,roubi1988resonance,cingolani1987raman}. There is an agreement in literature on the attribution of peaks at 266 cm$^{-1}$ and 344 cm$^{-1}$ to in-plane E$_g$ and out-of-plane A$_{1g}$ Raman-active modes, respectively. The peak at 325 cm$^{-1}$ is attributed to infrared-active E$_u$ (LO) \cite{cingolani1987raman,neal2021Raman} or A$_{2u}$(LO) \cite{roubi1988resonance} alike mode.
Despite the formal Raman-inactivity, both LO and TO components of infrared-active modes may be observed in RS spectrum due to long range Coulomb forces resulting from charges localised on atoms \cite{roubi1988resonance}.
A resonant nature of RS involving such phonons was previously reported in Ref. \citenum{iwasaki1982Raman}.
Temperature-dependent RS measurements \cite{golasa2017resonance} or measurements with modulation of excitation light \cite{Molas2017} might help to verify the attribution of the mode.
For the sake of this work, we assume the latter attribution to A$_{2u}$(LO) \cite{ibanez2018high}.

The assignment of the weak feature at 141~cm$^{-1}$, previously reported in Ref.~\citenum{ibanez2018high}, is also not clear.
Its attribution to a second-order difference processes was suggested, as such peaks can be observed in RS spectra of some TMDs \cite{chen1974, golasa2013}.
The E$_{u}$(TO) was also suggested as the signature of the mode \cite{roubi1988resonance, neal2021Raman}.
The peak in this work will be referred to as $\omega_{1}$.

The effect of hydrostatic pressure on RS spectra of bulk $\text{HfS}_{2}$  during compression and decompression can be appreciated in Fig.~\ref{fig:1}(a) and (b), respectively.
Let us first focus on the results obtained during compression process.
The lineshape of the RS spectra does not significantly change up to 5.7 GPa.
The observed mode experience blueshifts with increasing of the hydrostatic pressure. Moreover, a substantial enhancement of the A$_{1g}$ peak can be noticed in the spectrum at 5.7 GPa.

\begin{figure}[h]
	\centering
	\includegraphics[width=.45\textwidth]{figure_3.pdf}
\caption{The evolution of the Raman scattering peak energies in bulk $\text{HfS}_{2}$  with hydrostatic pressure. Straight lines correspond to linear fits to experimental data.}\label{fig:3}
\end{figure}

It was found that the energies of the observed modes increase linearly with increasing pressure, which can be appreciated in Fig.~\ref{fig:3}.
Assuming the linear dependence of the energies on pressure:\\
\vspace{-35pt}
\begin{center}
\begin{equation*}
     E(P)=E_0+\alpha P
\end{equation*}
\end{center}
\vspace{-5pt}
\noindent one can fit the experimental data with the results shown in Table \ref{tab:table}.

The obtained pressure coefficients are in general agreement with those recently reported by Ib{\'a}{\~n}ez and coworkers \cite{ibanez2018high}.
The initial values $E_0$ and pressure coefficients $\alpha$ for the observed modes span from 1.59~cm$^{-1}$/GPa to 4.36~cm$^{-1}$/GPa for A$_{1g}$ mode.

\begin{table}[t]
\caption{\label{tab:table} Results of linear fits to the pressure evolution of the Raman scattering peaks in bulk $\text{HfS}_{2}$  }
\begin{ruledtabular}
\begin{tabular}{c c c}
 Peak & $E_0$ &  $\alpha$ \\
  & [cm$^{-1}$] & [cm$^{-1}$/GPa] \\ \hline
  $\omega_1$ & 141.4 & 1.59 \\
 E$_g$ & 266.1 & 2.14 \\
 A$_{2u}$(LO) & 325.5 & 2.63\\
 A$_{1g}$ & 344.0 & 4.36 \\\hline
 I & 95.1 & -0.13 \\
 II & 101.8 & 0.69 \\
 III & 142.0 & 0.75 \\
 IV & 151.6 & 0.77 \\
 V & 277.3 & 1.68 \\
 VI & 315.3 & 1.69 \\
 VII & 334.9 & 2.61 \\
\end{tabular}
\end{ruledtabular}
\end{table}

\begin{figure}[b]
	\includegraphics[width=.5\textwidth]{figure_4.pdf}
\caption{Raman scattering spectrum of bulk $\text{HfS}_{2}$  in hydrostatic pressure of 9.8~GPa.} \label{fig:4}
\end{figure}

Substantial change in the line-shape of the RS spectrum can be observed at pressure of 9.8 GPa, which corresponding spectrum is presented in Fig.~\ref{fig:4}.
It can be appreciated in the Figure that the~A$_{1g}$ mode substantially weakens.
Simultaneously E$_g$, A$_{2u}$(LO), and $\omega_{1}$ disappear from the spectra and a set of new lines emerge in their stead. The~peaks are denoted with I-VII.
Similar lineshape of the spectra can be observed under hydrostatic pressure up to 12.8~GPa with the~A$_{1g}$ mode completely vanishing above 10~GPa.
The structural change of $\text{HfS}_{2}$ , which affects the RS spectra, was previously suggested in  Ref.~\citenum{ibanez2018high} with a new peak on the higher-frequency side of the A$_{1g}$ mode.
In our experiment, the transformation manifests itself in the emergence of peak VII at the lower-frequency side of the A$_{1g}$, as well as in changes to the low-energy part of the spectrum, in which two pairs of lines: I-II and III-IV can be observed.
The pressure coefficients $\alpha$ of the I-VII modes do not correspond in to those of modes detected at lower pressure.
In particular, the pressure coefficient $\alpha$ (2.61~cm$^{-1}$/GPa) for mode VII is substantially lower than the pressure coefficient for A$_{1g}$ mode although both linear energy evolutions have rather similar $E_0$ initial values  (334.9~cm$^{-1}$ for VII mode vs 344.0~cm$^{-1}$ for A$_{1g}$).
Similar effect can observed for E$_g$ and V mode with $E_0$ and $\alpha$ equal respectively to 266.1~cm$^{-1}$ and vs 277.3~cm$^{-1}$ and 2.14~cm$^{-1}$/GPa vs 1.68~cm$^{-1}$/GPa (see Table \ref{tab:table}).
Modes II-IV share relatively small pressure coefficient $\alpha$ equal to approx. 0.7~cm$^{-1}$/GPa) which is significantly smaller than the lowest coefficient for modes below the pressure-induced transformation.
Finally, the pressure coefficient for the I mode is negative. The properties substantially differ from those observed at lower pressure one can conclude that the phase transition of yet unknown origin takes place between 5.7 GPa and 9.8 GPa.

Next transformation of the RS spectrum can be observed above 15.2 GPa (see Fig.~\ref{fig:1}(a)).
Instead of seven well-defined vibrational modes present in the RS spectrum at lower pressure, broad spectral bands can be detected at 15.2 GPa or higher pressures.
The overall lineshape of the RS spectra at highest hydrostatic pressure resembles that of disordered materials \cite{golasa2015disorder, minguzzi2015disorder}.
The disorder and/or defects affecting crystalline structure can localize phonons allowing for momentum non-conserving processes.
As a result, the inelastic scattering is possible with phonons from whole Brilloiun zone.
Therefore the RS spectrum of disordered materials reflect the total density of phonon modes.
Adopting such an approach one can conclude that the low-energy band (below 150~cm$^{-1}$) corresponds to acoustic phonons and the high-energy band (at approx. 300~cm$^{-1}$) is due to optical vibrations.
The emergence of the disorder-like RS under highest hydrostatic pressure may be related to the contribution of an undefined phase observed in X-ray diffraction measurements as it was mentioned earlier.
As can be appreciated in the Fig.~\ref{fig:1}(b), the lineshape of RS spectra does not change during decompression process.
It suggests permanent nature of the high-pressure transition, most likely into an amorphous phase.

\section{Conclusions \label{sec:Conclusions}}

We have studied Raman scattering in bulk $\text{HfS}_{2}$  as a function of hydrostatic pressure up to 27~GPa. There are two transformation of Raman scattering spectrum observed during compression between 5.7~GPa and 9.8~GPa as well as between 12.8~GPa and 15.2~GPa.
The spectrum after the lower-pressure transformation consists of seven well-defined vibrational modes as compared to four modes observed below the transformation.
The modification in the spectrum suggests structural change in the material.
Frequencies of the modes observed above the transformation change linearly with pressure and the corresponding pressure coefficients have been determined.
Eelatively small (II-IV) or negative pressure coefficients (I) of the vibrational modes in the low-energy part of the spectra have been noticed.
The high-pressure transformation of the Raman scattering spectrum results in the change of its lineshape.
While a series of well-defined vibrational modes are observed under pressure below the transition, broad spectral bands can be detected at higher pressure.
The lineshape of the spectra resembles that of disordered materials, which suggests amorphisation process involved.
The lineshape does not change during decompression confirming permanent nature of the high-pressure transition.

\section*{Acknowledgment}
The work has been supported by the National Science Centre, Poland (grants no.  2017/27/B/ST3/00205 and 2018/31/B/ST3/02111). Z.M. and W.Z. acknowledge support from the National Natural Science Foundation of China (Grants No. 61627813 and 61971024), the International Collaboration Project (No. B16001), the Beihang Hefei Innovation Research Institute (project no. BHKX-19-02).

\bibliographystyle{apsrev4-1}


Boron in the very metal-poor star BD-13 3442

Abstract: The Goddard High Resolution Spectrograph (GHRS) of the Hubble Space Telescope (HST) has been used to observe the boron 2500 A region of BD-13 3442. At a metallicity of [Fe/H]=-3.00, this is the most metal-poor star ever observed for B. Nearly 26 hours of exposure time resulted in a detection. Spectrum synthesis using the latest Kurucz model atmospheres yields an LTE boron abundance of log e(B)=0.01 +/- 0.20. This value is consistent with the linear relation of slope ~1.0 between log e(B[LTE]) and [Fe/H] found for 10 halo and disk stars by Duncan et al. 1997. Using the NLTE correction of Kiselman and Carlsson (1996), the NLTE boron abundance is log e(B)=0.93 +/- 0.20. This is also consistent with the NLTE relation determined by Duncan et al. (1997) where the slope of log e(B[NLTE]) vs. [Fe/H] is ~0.7. These data support a model in which most production of B and Be comes from the spallation of energetic C and O nuclei onto protons and He nuclei, probably in the vicinity of massive supernovae in star-forming regions, rather than the spallation of cosmic ray protons and alpha particles onto CNO nuclei in the general interstellar medium.

\section{Introduction}

The light elements lithium, beryllium, and boron are of great interest
out of proportion to their very low abundances, having implications
in Big Bang Nucleosynthesis and stellar structure, as well as in
constraints on models of galactic chemical evolution.

The ``canonical'' theory of the origin of the elements Li, Be, and B
was first presented by Reeves, Fowler, \& Hoyle (1970) and further
developed by Meneguzzi,
Audouze, \& Reeves (1971), and then Reeves, Audouze, Fowler, \&
Schramm (1973).
In this model, most light element formation can be accounted for by
galactic cosmic rays (GCR) impinging on the interstellar medium (ISM),
assuming a constant flux of GCRs through the life of the
Galaxy and making reasonable assumptions about CR confinement by the
Galactic magnetic field.  Meneguzzi et~al.\ (1971) also
introduced the idea of a large (up to three orders of magnitude)
increase in the low energy (5-40 MeV nucleon$^{-1}$) CR flux; since CRs
in this energy range are mostly shielded from the Solar System by the
solar wind, they are not detectable.  This additional CR flux increased
the production of all light elements and matched the isotopic ratios
and total abundances to the accuracy known at the time.

Reeves \& Meyer (1978) added the additional constraint that models
should match not only present-day abundances but their evolution
throughout the life of the Galaxy.  Their conclusions were similar to
MAR, except that they had to introduce infall of light-element-free
matter into the Galactic disk to match the evolution with time.  In
retrospect, it can be seen that the data they were fitting were sparse
and not very precise.  With the launch of the the {\it Hubble Space
Telescope} ({\it HST}) and the availability of uv-sensitive CCD
detectors (B is usually observed at $\lambda$2500 and Be at
$\lambda$3130), data are now much more numerous and accurate, and
abundances can be traced from the epoch of formation of the Galactic
halo until the present day.

In the past several years, the evolution of Li, Be, and B has been used
as a test of different models of the chemical and dynamical evolution
of the Galaxy.  For example, in the models of Vangioni-Flam et~al.\
(1990), Ryan et~al.\ (1992),
and Prantzos et~al.\ (PCV; 1993), light element production
depends on the intensity and shape of the GCR spectrum, which in turn
depends on the supernova (SN) and massive star formation rates.  It
also depends on the rise of the (progenitor) CNO abundances and the
decline of the gas mass fraction, which is affected by rates of infall
of fresh (unprocessed) material and outflow, e.g. by SN heating.  Other
things being equal, at early times when target CNO abundances were low,
light element production would be much lower for a given CR flux than
presently, when the ISM abundances are higher.  PCV found that even
with these numerous adjustable parameters, no time-independent CR
spectrum can reproduce the evolution of light element abundances.  By
assuming a very particular form of time variation of the CR flux
(greatly enhanced at early epochs), they
were able to (barely) fit the evolution of the abundances.

The present investigation supports a different solution to the problem
of the origin of the light elements.  Duncan et~al.\ (1997) present B
abundances in a large number of stars ranging in metallicity from
$\sim$solar
to [Fe/H]~$\sim -2.8$.  They find that (LTE) B follows metals in direct
proportion from the earliest times (very metal-poor stars) to the
present, with little if any change of slope between halo and disk
metallicities.  A straightforward interpretation of this is that the
rate of production of B and Be does {\it not} depend on the CNO
abundances in the ISM, and that the production site is
associated with the production site for metals.  This would be true if
the spallation process most important for light element production is
not primarily protons and $\alpha$ particles colliding with CNO nuclei
in the ISM but rather C and O nuclei colliding with
ambient protons and $\alpha$ particles, probably in regions of massive star
formation (cf. Vangioni-Flam et~al.\ 1996 and Ramaty et~al.\ 1997).

This paper focuses specifically on the B measurement in \bd-13; it is
consistent with (and was used to help determine) the
relationships between LTE and NLTE B and [Fe/H] seen in Duncan et~al.\ (1997).
The data suggest that B production at the lowest metallicities occurs
in the same way as today, and thus support the new
description of galactic light element production.

It is possible that recent GRO satellite observations of gamma rays
from the Orion Nebula (Bloemen et~al.\ 1994; Bykov \& Bloemen 1994)
provide direct evidence of C and O spallation occurring today, even
though not all instruments on GRO detected evidence of such spallation
(Murphy et~al.\ 1996).  The light element data alone, however, are the
strongest evidence in support of a new model of their production.

\section{Observations and data reduction}

The Goddard High-Resolution Spectrograph (GHRS) of the {\it Hubble
Space Telescope} ({\it HST}) was used with the G270M grating to obtain
spectra of resolution 26\,000 (0.025 \AA\ pixel$^{-1}$) in the
$\lambda$2500 B\,{\sc i}\ region of \bd-13.  Data was collected for a total
of 25.39 hours exposure time in six separate visits of the satellite
over several months from December 1994 to May 1995.  The S/N of the
final spectrum was 100 per pixel (200 per resolution element), limited
only by photon statistics.  The spectrum obtained appears in Fig.~1.

These data were combined with that from
10 other stars, including three previously analyzed
by Duncan, Lambert, \& Lemke (1992) and seven new stars
([Fe/H] = $-$0.45 to $-$2.75), in a separate
investigation (Duncan et~al.\ 1997).  Special care was taken
to analyze each of the stars in a similar manner, and the stars from
Duncan, Lambert, \& Lemke (1992) were re-analyzed to assure
consistency.  See Duncan et~al.\ (1997) for more details.

Data reduction followed the standard {\it HST} procedure, using the
IRAF {\sf stsdas} package, combining the quarter-stepped,  FP-SPLIT
data within each visit
with the tasks {\sf poffsets} and {\sf specalign}, which perform a
generalized least-squares solution for the photocathode granularity and
the true spectrum.
Since {\sf specalign} equally weights the input
spectra, exposures from separate visits were shifted and combined using
the separate tasks {\sf shiftlines} and {\sf scombine}, weighting the
data by the number of counts in the raw spectrum from each visit.  By
comparing subsets of the data, it was confirmed that the photon
statistics limit was reached and that the line in question was present
in both data sets if the data was arbitrarily divided in half, and not
substantially altered by any exposure or by texture on the
photocathode.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{figs/fig1.eps}}
\caption{Spectrum of \bd-13. The instrumental resolution is 0.096 \AA,
and the observational data shown are oversampled by a factor of about 2.}
\end{figure}

\section{Abundance analysis}

Spectrum synthesis using the latest Kurucz model atmospheres was used
to determine log $ \epsilon $(B), with particular attention paid to the errors.
An accurate determination of the uncertainties is required in the
subsequent discussion of the results, to permit judgement of the
goodness-of-fit of models of Galactic chemical evolution and to
discriminate among the different proposed mechanisms of light element
formation, particularly because this star is the most metal-poor yet
investigated for boron.

\subsection{Spectrum synthesis}

Spectrum synthesis was done using the SYNTHE program distributed by
Kurucz (1993) on CD-ROM.  This program assumes local thermodynamic
equilibrium (LTE) in determining level populations and calculating the
emergent spectrum.  Scripts written by Steve Allen (UC Santa Cruz) were
used to run the program on Unix SparcStations. The most recent grid of
model atmospheres, also released by Kurucz on CD-ROM, includes the
blanketing of almost 60 million lines, both atomic and molecular.  As
in our parallel study of the Galactic evolution of boron, (Duncan
et~al.\ 1997), we used the linelist of the Duncan et~al.\ (1997b) study of
B in the Hyades giants, which consists almost entirely of
laboratory-measured lines.  It has been carefully tested in order to
fit both the Hyades giants and various metal-poor stars in the B
$\lambda$2500 region.  This list was constructed by starting with
Kurucz's LOWLINES list (CD-ROM \#~1), adjusting (reducing) some
{\it gf}\ values where necessary, and adding some predicted Fe\,{\sc i}\ lines in
order to achieve a good match to the absorption features near the
B\,{\sc i}\ lines.  The {\it gf}\ value for the Co\,{\sc i}\ 2496.708~\AA\ line which
blends with the short wavelength component of the B\,{\sc i}\ doublet is the
same as that of the higher resolution study of Edvardsson
et~al.\ (1994).  In our linelist, the {\it gf}\ values and wavelengths of the
B\,{\sc i}\ resonant doublet were selected from O'Brian \& Lawler (1992).
They differ by only 0.03~dex and $<$0.01~\AA\ from the values given by
Kurucz.  Our abundances are derived from fitting the shorter
wavelength, less blended component of the B doublet, with the longer
wavelength component examined as a consistency check.

As in Duncan et~al.\ (1997), the Kurucz default abundance of
log $ \epsilon $(B)=2.57 is used as our solar B value.  This decision
affects the choice of the zero point of the abundance scale and does
not affect the interpretation of trends in the B data discussed in the
present paper.

The Kurucz model grid is given in steps of 250~K in ${T_{\rm eff}}$, 0.50 in
$\log g$, and 0.50~dex in [Fe/H] (at least for halo
stars).  We determined the
physical parameters and their uncertainties for the program star from
an extensive search of the literature and performed calculations at the
nearest grid points.  Within the range of values determined from the
literature, models from the Kurucz grid which best fit our entire uv spectrum
were chosen.
Experiments with linear interpolation between grid models
were also performed.

The Kurucz models we used assume all metals scale together, whereas it
is known that $\alpha$-elements are enhanced at low metallicities.  In
Duncan et~al.\ (1997), testing showed that the effect of
$\alpha$-enhancement on the structure of the model atmosphere produces
very small effects on the calculated B abundances.  In particular, the
three stars originally studied by Duncan, Lambert, \& Lemke (1992) were
analyzed by them using atmospheres with $\alpha$-enhanced composition
from Gustafsson et~al.\ (1975).  Analyzing them
with the methods of the present paper produced differences in B
abundance within 0.04~dex, much lower than the other errors
associated with the abundance determination.  The direct effect of an
enhanced Co abundance could be larger due to the particular line
Co\,{\sc i}\ $2496.708$ \AA\ which is blended with the shorter wavelength B
component.  This is discussed further in the section on uncertainties.

The atmospheric models we used were computed with a microturbulent
velocity of 2~km~s$^{-1}$. We did not find it necessary to change to the
commonly accepted value of 1.5~km~s$^{-1}$ for halo stars (Magain 1989),
because analysis of several curves of growth showed that the dependence
of B abundances on this stellar parameter is negligible,
though microturbulence would affect the blending Co line, which we
discuss more fully below.

Fig.~2 presents spectrum synthesis fits for \bd-13; the fits
presented are calculated to have log $ \epsilon $(B)=+0.3, +0.1, and $-$0.4,
the latter indicates what the abundance would be if this star
had the solar B to metals ratio.  Table~1 presents
the stellar parameters and B abundance for \bd-13.

Table~1 also presents parameters from the three most metal-poor stars
of our recent study
(Duncan et~al.\ 1997):
HD~140283, which for several years was the most metal-poor star
with detected B (Duncan et~al.\ 1992), and BD~$+$3~740.
The metallicities determined by us for these stars are consistent
with the ranges found by other workers (see compilation by Cayrel de
Strobel et~al.\ 1997).
Fig.~3 compares the observed spectra of \bd-13\ with that of
HD~140283 and BD~$+$3~740.
It is clear that \bd-13\ is the most metal-poor of the three stars.

\begin{table}
\caption{Stellar parameters and B abundances}
\begin{flushleft}
\begin{tabular}{cccccc} \hline
star & ${T_{\rm eff}}$ & $\log g$ & [Fe/H] & log$\epsilon(B)_{LTE}$ &
	log$\epsilon(B)_{NLTE}$ \\
\hline
\bd-13 & 6250  & 3.75 & $-$3.00 & 0.01 & 0.93 \\
BD~+3~740 & 6125 & 3.5 & $-$2.75 & 0.22 & 1.04 \\
HD~140283 & 5640 & 3.5 & $-$2.60 & $-$0.10 & 0.35 \\
\hline
\end{tabular}
\end{flushleft}
\end{table}

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{figs/fig2a.eps}}
\resizebox{\hsize}{!}{\includegraphics{figs/fig2b.eps}}
\caption{{\bf a} Spectrum synthesis fits for \bd-13. {\bf b} Enlarged
view of fits near the B lines.}
\end{figure}

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{figs/fig3.eps}}
\caption{Comparison of \bd-13, HD~140283, and BD~$+$3~740.}
\end{figure}

\subsection{Uncertainties}

Considerable care was taken in determining the error of the abundance
determination. Sources of error we considered included uncertainties in
stellar ${T_{\rm eff}}$, $\log g$, and [Fe/H] (all of which change the atmosphere
used in the analysis), continuum placement, photon statistics in the
points defining the line itself, and the effect of metal lines which
blend with the B features.

In such a metal-poor star the continuum is relatively easy to define
but the B line is weak.  Table~2 presents, in order, effects due to
uncertainties in stellar ${T_{\rm eff}}$, $\log g$, [Fe/H], continuum location and
photon statistics in the B line, as well as the total (net) error.  We
quantified the dependence of the B abundance on changes in the main
stellar parameters by running multiple syntheses and curves of growth,
although curves of growth are less reliable since the B lines are
blended. From available literature determinations for this
star, we considered $\pm$75~K, $\pm$0.30, and $\pm$0.10~dex to be
typical uncertainties to be associated with ${T_{\rm eff}}$, $\log g$, and [Fe/H]
respectively.
Errors for BD~+3~740 and HD~140283 are also included in Table~2
for comparison.

Uncertainty in stellar metallicity in metal-poor stars causes error due
to lines which blend with the B features.  At our resolution, the
Co\,{\sc i}\ 2496.708 \AA\ feature blends with B\,{\sc i}\ 2496.772 \AA.  Assumption
of a higher metallicity attributes more of the blend to Co, decreasing
the derived B abundance, and vice-versa for a lower assumed
metallicity.  The column ``[Fe/H]'' in Table~2 gives the amount of this
direct effect when the metallicity is uncertain by $\pm$0.1 dex.

McWilliam et~al.\ (1995) studied Co in their spectroscopic
analysis of the most metal-poor stars.  They present evidence that Co
remains scaled-solar to [Fe/H] $\sim -$2.5, but that it
increases by $\sim$0.5 dex as [Fe/H] decreases to $-$3.0.  Since our
spectra do not resolve the Co-B blend, overabundance of Co means that
the B abundances should be decreased.  However, the blend between the
Co and B features is only partial and an increase in Co of a given amount
produces a decrease in B abundance of a smaller amount.
Synthesis testing with [Co/H]
increased by 0.1 and 0.3 dex caused a decrease in B abundance of
about 0.05 and 0.15 dex respectively.  The data are not consistent with
[Co/H] enhanced by 0.5 dex.  Since we have no way
of deblending the lines, this error is not included in Table~2;
it should be considered a systematic
effect which could affect the derived abundance, and at worst case,
our abundance represents an upper limit.
That this effect is unlikely to be large is supported by the good
agreement of the B abundance derived by Duncan et~al.\ (1997) for the
slightly more metal-rich star HD~140283 with that of the higher
resolution study of Edvardsson et~al.\ (1993), which partially resolves
the Co-B blend.  The spectrum for \bd-13 is similar in resolution to
that for HD~140283 in Duncan et~al.\ (1997).

\begin{table}
\caption{Uncertainties affecting the derived B abundance}
\begin{flushleft}
\begin{minipage}{3in}
\begin{tabular}{ccccccc} \hline
star & ${T_{\rm eff}}$ & $\log g$ & [Fe/H] & Continuum & Photon & {\bf Net}  \\
    &$\pm$75~K & $\pm$0.30 & $\pm$0.10 & $\pm$0.5\% & Statistics& {\bf Error}\\
\hline
\bd-13 & $\pm$0.09 & $\pm$0.03 & $\pm$0.10 & $\pm$0.10 & $\pm$0.10 &
	{\bf $\pm$0.20} \\
BD~+3~740\footnote{Final boron abundance 0.22 LTE and 1.04 NLTE.} &
	$\pm$0.08 & $\pm$0.05 & $\pm$0.10 &
	$\pm$0.10\footnote{$\pm$1.5\% in the continuum.} & $\pm$0.23
	& {\bf $\pm$0.29} \\
HD~140283\footnote{Final boron abundance $-$0.10 LTE and 0.35 NLTE.} &
	$\pm$0.08 & $\pm$0.05 & $\pm$0.10 &
	$\pm$0.10\footnote{$\pm$1.5\% in the continuum.} & $\pm$0.06 &
	{\bf $\pm$0.18} \\
\hline
\end{tabular}
\end{minipage}
\end{flushleft}
\end{table}

\subsection{NLTE abundances}

The largest source of systematic errors in our B abundance is likely to
be due to NLTE effects.  As calculated by Kiselman (1994), and Kiselman
\& Carlsson (1996), NLTE effects increase B abundances most
significantly in the most metal-poor stars.  The relatively hot uv
radiation fields in metal-poor stars drive the detailed balance of the
B resonance lines away from statistical equilibrium.  Overionization
and optical pumping tend to weaken the B\,{\sc i}\ 2498 \AA\ doublet, which
leads to an underestimate of the B abundance.  NLTE calculations are
difficult, being limited both by the incomplete understanding of the
atomic physics involved and especially by the uncertainties in the uv
flux values at the precise absorption wavelengths.  The Kiselman \&
Carlsson (1996) study is improved over the Kiselman (1994) one in that
it models the effects of large numbers of lines in the background
radiation field, which tend to decrease NLTE effects.   Results of the
two studies differ by approximately 0.1 dex; present NLTE calculations
are necessarily less accurate than that amount.  The Kiselman \&
Carlsson NLTE correction for \bd-13 is $+0.92$ dex.

\section{Discussion}

Fig.~4 shows the
LTE abundances as a function of [Fe/H] from
stars analyzed in the larger investigation of Duncan
et~al.\ (1997), with the point for \bd-13 from the present investigation
emphasized.  It can be seen that there is an approximately linear
relation between $\log$ $ \epsilon $(B$_{\rm LTE}$) and [Fe/H] over both disk
and halo metallicities, and that \bd-13\ is consistent with this
relationship.  A least-squares fit to all the data of Fig.~4 (allowing
for errors in both coordinates) yields a slope of 0.96$\pm$0.07 and a
reduced chi square, $\chi_{\nu}^2$, of 0.71, indicating an excellent
fit.  If NLTE abundances are used, as shown in Fig.~5, the slope is
0.70$\pm$0.07 and $\chi_{\nu}^2 = 1.63$.  Although it is true that
\bd-13\ is used to determine this line, it is important to note that
this star is quite consistent with the trend defined by the other
stars.  Inspection and $\chi_{\nu}^2$ tests confirm this.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{figs/fig4.eps}}
\caption{LTE B abundances from Duncan et~al.\ (1997) with the program
star highlighted.}
\end{figure}

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{figs/fig5.eps}}
\caption{NLTE B abundances from Duncan et~al.\ (1997) with the program
star highlighted.}
\end{figure}

\subsection{Comparison to standard models and a new model for light
element production}

The slope of close to 1 suggesting a primary process is {\it not}
expected from canonical models of CR spallation in the ISM, which
predict a secondary process and thus a steeper relation.  In a
secondary process the rate of light element production depends on the
product of the abundance of target CNO nuclei and the CR flux, both of
which vary with time.  If SNe are the source of the target nuclei and
the ISM is well-mixed, the ISM metallicity is proportional to the integral
(total) number of SN up to a given time.   If, as is commonly supposed,
SNe also seed the acceleration mechanism which produces CRs, the CR
flux is proportional to the SN rate.  The result is light element
abundances which vary quadratically with the metallicity of the ISM, or
a logarithmic slope of 2 (Prantzos et~al.\ 1993).
Figs.~4 and 5 show that such a slope is certainly not consistent with
our data.
Duncan et~al.\ (1997) discuss these issues in greater detail.

As was discussed by Duncan, Lambert, and Lemke (1992), the data for the
first three metal-poor stars observed for B already seemed to show a
linear (in the log) relationship with [Fe/H], suggesting some primary
production mechanism rather than the secondary mechanism in the ISM
described above. This idea has been modelled in detail by Cass\'e
et~al.\ (1995), Ramaty et~al.\ (1995 and 1997), Lemoine et~al.\ (1997),
and Vangioni-Flam et~al.\ (1996).  In the new scenario, B and Be are
primarily produced by the spallation of C and O onto protons and
$\alpha$ particles.  Such a process could occur near massive star SNe,
where the particle flux would be very non-solar in composition;
depleted in H and He and especially enriched in O and C.  Vangioni-Flam
et~al.\ find that a composition matching either winds from massive
(Wolf-Rayet; WR) stars in star-forming regions or massive star SNe
produce a flux of O and C which, after further acceleration, can
reproduce both the magnitude and slope of B production seen in Figs.~4
and 5 through collisions with protons and $\alpha$ particles.  As
Ramaty et~al.\ point out, production of some additional $^{11}$B by the
neutrino process (Woosley et~al.\ 1990) is not ruled out, and may be
favored on energetic grounds.  Nevertheless the bulk of the B and Be
would be produced from the spallation process.

Although the NLTE correction to the B abundance of \bd-13 is relatively
large and tends to raise the B abundance above the line in Fig.~5, two
other effects not included here would tend to move it closer to the
curve.  One is the effect of the blending Co line discussed above,
which could reduce the B abundance as much as $\sim$0.15 dex.  The
other is the fact that if the spallation producing light elements is
caused by O (and to a lesser extent C), $\log$ $ \epsilon $(B) should be
plotted against [O/H] rather than [Fe/H].  As is well-known, very
metal-poor stars are overabundant in O compared to Fe (moving points
representing the most metal-poor stars to the right in a figure with O
on the x-axis).  Duncan et~al.\ (1997) make such a plot, and demonstrate
that although the measurement errors in O are greater than those for
Fe, when all the metal-poor stars are considered together a straight
line of slope 1.10$\pm$0.14 fits the LTE abundances, and one of slope
0.82$\pm$0.10 the NLTE abundances.
However, oxygen abundance measurements are also
surrounded by greater uncertainty than are iron measurements, and
systematic errors in the oxygen abundance which depend on metallicity will
affect the derived slope.

\section{Summary}

The present investigation has presented an analysis of B in \bd-13.  At
a metallicity of [Fe/H]=$-3.0$, this is the most metal-poor star ever
observed for B.  Careful spectrum synthesis was performed to determine
both the B abundance and its uncertainty. NLTE corrections were also
considered.

We find that \bd-13\ is consistent with the straight line of slope
close to 1 between log $ \epsilon $(B) vs.\ [Fe/H] or log $ \epsilon $(B) vs.\ [O/H]
found by Duncan et~al.\ (1997).  This suggests that the production of
light elements such as B and Be is directly related to the production
of heavier elements.  Canonical models which only include spallation
by the process of protons and $\alpha$ particles impinging
onto the general Galactic ISM predict
a steeper relation which fails to fit all the data now available.

Spallation via the reverse process of energetic C and O nuclei onto
protons and $\alpha$ particles does much better
in explaining the linear relationship observed.  This reverse process
was historically discarded as a minor effect (e.g. 20\% according to
the calculations of Walker et~al.\ 1993) since the canonical theory
assumed a CR and ISM composition similar to that in the solar system.  A
population of low-energy CRs
enriched in C and O and depleted in H and He can produce light
elements by spallation onto the ambient medium.  The
metallicity of the medium would have little effect, and the process
would operate the same throughout the life of the Galaxy, leading to
the observed slope of 1.  A model in which most light element
production comes from CR spallation of C and O nuclei onto protons and
$\alpha$ particles in the vicinity of massive supernovae in
star-forming regions fits the data well.  If particle compositions are
taken from calculations of either massive star SNe or WR stars,
spallation can reproduce both the magnitude and slope of the observed B
evolution (Ramaty et~al.\ 1996) while maintaining consistency with the
other light element abundances.
The observations of B in \bd-13 provide
the strongest evidence of how CR spallation operated in the very early
galaxy.

\begin{acknowledgements}
We thank Michel Cass\'e, Reuven Ramaty, Hubert Reeves, Elizabeth
Vangioni-Flam, and Stan Woosley for their comments.

This research was based on observations obtained with the NASA/ESA {\it
Hubble Space Telescope} through the Space Telescope Science Institute,
which is operated by the Association of Universities for Research in
Astronomy, Inc., under NASA contract NAS5-26555.  This research has
made use of NASA's Astrophysics Data System Abstract Service.

CPD gratefully acknowledges support provided by the University of Hawaii
Foundation, and subsequently by NASA through grant HF-1042.01-93A from the
Space Telescope Science Institute.

\end{acknowledgements}


Powerful relativistic jets in spiral galaxies

Abstract: The discovery of high-energy (E>100 MeV) gamma rays from Narrow-Line Seyfert 1 Galaxies (gamma-NLS1s) has confirmed the presence of powerful relativistic jets in this class of active galactic nuclei (AGN). Although the jet emission is similar to that of blazars and radio galaxies, gamma-NLS1s have some striking differences: relatively small masses (10^6-10^8 M_sun), high accretion rates (0.1-1 times the Eddington limit) and are generally hosted by spiral galaxies. It is now possible to study a rather unexplored range of mass and accretion rates of AGN with relativistic jets. Specifically, in this work I present some results obtained by comparing a sample of blazars and gamma-NLS1s with another sample of Galactic binaries with relativistic jets (stellar mass black holes and neutron stars).

\section{Introduction}
During the latest decades, we have observed powerful relativistic jets from active galactic nuclei (AGN), Galactic binaries (stellar mass black holes, neutron stars and cataclysmic variables), and Gamma Ray Bursts (GRBs). Supersonic jets are also visible in protostellar systems. Despite of the diffusion of jets in almost all types of accreting astrophysical systems, it seemed that galaxies made some preference: early observations of the host galaxies of AGN with relativistic jets displayed only elliptical shapes\cite{para1,para2,para3,para4}. Only a few counterexamples (i.e. jets in spiral galaxies) were found in the past. The turning point occurred a couple of years ago with the discovery -- by means of the \emph{Fermi} satellite -- of GeV $\gamma$ rays from narrow-line Seyfert 1 galaxies ($\gamma$-NLS1)\cite{abdo1,abdo2,fosco1}. This class of AGN is composed of sources generally hosted by spiral (barred) galaxies\cite{crenshaw,deo} and, according to the common paradigm, without jets. Therefore, although some early radio observations had suggested the presence of some jet emission in a handful of cases\cite{grupe,komossa}, the discovery of GeV $\gamma$ rays was somehow a surprise. It confirmed that also this class of AGN can develop powerful relativistic jets, like blazars and radio galaxies (see Ref.~\refcite{fosco2} for a review).

Presently, seven $\gamma$-NLS1s are known and surely one of them is hosted by a spiral galaxy\cite{zhou1,anton,fosco2}. Work is in progress to confirm the spiral hosts for the remaining $\gamma$-NLS1s. The fact that no NLS1 is known to be at high redshift\footnote{The 13th edition of the V\'eron-Cetty \& V\'eron catalogue of quasars and AGN\cite{veron} contains only NLS1s with $z\lesssim 1$. This is also in agreement with the findings that these sources are very young\cite{mathur}.} suggests that it is very unlikely to find this type of AGN hosted by ellipticals.

The issue of the host is linked to how the mass can affect the generation and the power of the jet and was studied in the past particularly with the use of the radio loudness parameter ($RL=f_{\rm 5~GHz}/f_{\rm 440~nm}$), which in turn should indicate the dominance of the synchrotron emission from the jet over the radiation from the accretion disc (e.g. Ref.~\refcite{para3,para4}). Being the blazars and radio galaxies ($RL>>10$) confined in a mass range of $10^{8-10}M_{\odot}$, while Seyferts and LINERs ($RL<10$) occupy the range $10^{6-8}M_{\odot}$, it seemed that the jet needed of a large mass to be generated\cite{para3,para4}.

The discovery of powerful relativistic jets in $\gamma$-NLS1 is now filling the region of small masses and high radio loudness (see Ref.~\refcite{fosco3} for more details). In addition, it opens also new interesting questions on the comparison of jets in AGN and in Galactic binaries. This is the topic of the present work and is complementary to Ref.~\refcite{fosco3}.

\section{Sample selection}
I have used the same sample of AGN in Ref.~\refcite{fosco3}, made of 9 BL Lac Objects, 30 flat-spectrum radio quasars (FSRQs) and 7 $\gamma$-NLS1s. For the present work, I have added a sample of Galactic objects (X-ray Binaries, XRB), made of three stellar mass black holes (GX~339$-$4\cite{corbel1}, H1743$-$322\cite{coriat}, V404~Cyg\cite{corbel2}) and two neutron stars (Aql~X-1\cite{miller}, 4U~1728$-$34\cite{migliari}). For XRB, the X-ray flux measured in different bands ($3-9$~keV in most of the cases), which samples the accretion radiation in Galactic objects, has been converted into the $2-10$~keV flux with {\tt webPIMMS}\footnote{{\tt http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html}} by using a power-law model with photon index equal to 2. Then, the values have been corrected according to Ref.~\refcite{migliari} and scaled to the Eddington luminosity, calculated by assuming a mass of $1.4M_{\odot}$ for neutron stars and $5.8M_{\odot}$, $10M_{\odot}$, and $11.7M_{\odot}$ for GX~339$-$4\cite{gallo}, H1743$-$322\cite{coriat}, and V404~Cyg\cite{gallo}, respectively.

The radio emission has been converted into jet power by assuming the classical relationship\cite{blandford,kording} according to which $P_{\rm jet}\propto L_{\rm core,radio}^{12/17}$. For the sake of simplicity and homogeneity with the AGN sample, I have used the Eq.~(2) of Ref.~\refcite{fosco3}, taking into account that the radio emission from jets has a rather flat spectral index ($\alpha_{\rm radio}\sim 0$). The luminosities have been calculated by assuming the distances and the $N_{\rm H}$ reported in Refs.~\refcite{coriat,gallo,migliari}.

\begin{figure}[!t]
\centerline{\includegraphics[angle=270,scale=0.4]{foschini_1_f1.eps}}
\vspace*{8pt}
\caption{Jet power per unit of $M^{1.4}$ vs accretion. Data from AGN are red, while XRB are blue. Squares indicate BL Lac Objects, circles are for FSRQs, and stars are $\gamma$-NLS1s. Triangles are stellar mass black holes (in different states) and asterisks are for neutron stars. See the text for more details.  \label{f1}}
\end{figure}

\section{Jet power and accretion}
As already noted in Ref.~\refcite{fosco3} (Fig.~5), the addition of $\gamma$-NLS1s to the blazars makes this sample of AGN directly comparable with a sample of jets from Galactic compact objects\footnote{In Ref.~\refcite{fosco3}, I compared Fig.~5 with an analogous figure in Ref.~\refcite{coriat} (incidentally, still Fig.~5!).}. Blazars can be compared with stellar mass black holes, while $\gamma$-NLS1s occupy the region of neutron stars. Therefore, in term of relativistic jets, $\gamma$-NLS1s are the low mass counterpart of blazars so as neutron stars are the small mass version of Galactic black holes. The difference of the jet power is about 12-13 orders of magnitude, more than the mass difference between the two samples. To remove the mass dependence, I have calculated the jet power per unit mass. Then, I have searched for correlations in the two samples by using the {\tt ASURV~v.~1.2} software package\cite{asurv1,asurv2,asurv3}. In the case of AGN, I obtain:

\begin{equation}
{\rm {\bf AGN:}} \quad \log \frac{P_{\rm jet}}{M} = (38.3 \pm 0.1) + (0.55\pm 0.06)\log \frac{L_{\rm disc}}{L_{\rm Edd}}
\label{agn}
\end{equation}

\noindent and the results for the statistical tests are: $Z=4.1$ ($P_{\rm chance}<10^{-4}$) for the Kendall's test, while the Spearman's method gives $\rho =0.6$ ($P_{\rm chance}=2\times 10^{-4}$).

In the case of Galactic binaries:

\begin{equation}
{\rm {\bf XRB:}}\quad \log \frac{P_{\rm jet}}{M} = (34.8 \pm 0.1) + (0.33\pm 0.03)\log \frac{L_{\rm disc}}{L_{\rm Edd}}
\label{xrb}
\end{equation}

\noindent with Kendall's $Z=6.7$ and Spearman's $\rho=0.7$, both with $P_{\rm chance}<10^{-4}$.

\begin{figure}[!t]
\centerline{\includegraphics[angle=270,scale=0.4]{foschini_1_f2.eps}}
\vspace*{8pt}
\caption{Jet power per unit of $M^{1.4}$ and $(L_{\rm disc}/L_{\rm Edd})^{0.5}$ vs accretion. Symbols are the same of Fig.~\ref{f1}.  \label{f2}}
\end{figure}

The difference in normalization of the jet power per unit mass between AGN and XRB is now about three orders of magnitude. This difference, in excess of the difference of mass ranges, has been already noted by several authors (e. g. Refs.~\refcite{heinz,merloni}). Particularly, Ref.~\refcite{heinz} explained it by invoking a non-linear dependence on the mass, founding that $L_{\rm radio}\propto M^{17/12}$, in the hypothesis of the scale invariance of the jet. Indeed, by dividing the jet power by $M^{1.4}$, the difference in the normalization between AGN and XRB is removed (Fig.~\ref{f1}) and the correlation is now:

\begin{equation}
{\rm {\bf AGN+XRB:}}\quad \log \frac{P_{\rm jet}}{M^{1.4}} = (34.70 \pm 0.07) + (0.44\pm 0.03)\log \frac{L_{\rm disc}}{L_{\rm Edd}}
\label{all}
\end{equation}

\noindent with Kendall's $Z=9.1$ and Spearman's $\rho=0.7$, both with $P_{\rm chance}<10^{-4}$.

Interestingly, Fig.~\ref{f1} shows also that the most efficient engines are the $\gamma$-NLS1s, with the greatest jet power per unit $M^{1.4}$ (also per unit mass): small, but nasty. This is agreement with the theory that the central engine of AGN with jets is able to distribute more or less equal power between accretion and ejection\cite{rawlings,ghise3} and, therefore, since $\gamma$-NLS1s have the greatest accretion, they are at the top end of the distribution of specific jet power. In addition, the fact that $\gamma$-NLS1s and FSRQs have very high jet powers and accretion rates, and that the Blandford-Znajek mechanism tends to saturate in a radiation-pressure dominated regime\cite{ghosh} suggest that the jet needs of an additional contribution, perhaps like Blandford-Payne or other hybrid mechanisms\cite{fosco3}. The factor 0.44 of the accretion luminosity displayed in Eq.~(\ref{all}) suggests that it is possible to remove the residual dependence shown in Fig.~\ref{f1} by dividing the specific jet power by the term $(L_{\rm disc}/L_{\rm Edd})^{0.5}$. This is indeed what is shown in Fig.~\ref{f2}. This is a fact. Its interpretation is much more cumbersome.

This term indicates something that is not clear yet. It could point to the efficiency $\eta$ of a spherical accretion flow\cite{frank} and, thus, could emphasize the importance of Bondi accretion\cite{balmaverde,narayan2}. Another option is that it is proportional to the efficiency $\eta$ of the accretion disc, which in turn would point to the presence of some advection. Not necessarily an advection-dominated accretion flow (ADAF), but also a slim disc for the greatest rates. The hypothesis is that AGN could have composite discs, as like as XRB\cite{narayan1}, with a standard disc plus an advective part, which could be ADAF or slim depending on the accretion rate. The possibility that $(L_{\rm disc}/L_{\rm Edd})^{0.5}\propto \eta$ is appealing, because the basic concept is that high disc efficiency means high number of photons that can act both as direct and indirect (e.g. through the ionization of the broad-line region, BLR) source of seed photons to efficiently cool the relativistic electrons of the jet. It is known\cite{elitzur} that the radial scale of the BLR depends on $(L_{\rm disc}/L_{\rm Edd})^{0.5}$ and, then, the slope factor in Eq.~(\ref{all}) could be a direct indicator of the seed richness of the environment close to the central black hole (and, this time, indirectly linked to the radiative efficiency of the disc). Obviously, such an argument could be of some interest for AGN, but not for XRB that do not have a BLR.

On the other hand, it is worth noting that the factor $0.44$ in Eq.~(\ref{all}) is just the average of the factors $0.55$ in Eq.~(\ref{agn}) and $0.33$ in Eq.~(\ref{xrb}). The physical explanation of these factors could be indeed different, indicating a separated dependence on the disc luminosity in Eddington units, although the final aim is the same. In the case of AGN, the factor $0.55$ could be linked to the BLR (and indirectly to the disc), while in the case of XRB, the factor $0.33$ could (directly) be the average for the efficiency of advective discs in different states. It is necessary to take into account that while FSRQs evolve into BL Lacs through gigayears, Galactic black holes change their states on human time scales. $\gamma$-NLS1s are on a different evolutive track as like as neutron stars are not linked to stellar mass black holes. Therefore, if the different slopes indicate different relationships (direct or indirect) with the radiative efficiency, this could be due to the different time scales.

Hypotheses, just hypotheses. How exactly the accretion disc contributes to the jet power is still to be understood. However, despite all these differences and unknowns, Fig.~\ref{f2} confirms the idea that the jets developed by all these systems are basically the same thing.

\section*{Acknowledgments}
I would like to thank J.~M. Paredes, the LOC and the SOC of HEPRO3 for having organized a really pleasant conference and having given me time to express my ideas. Thank you also to F. Tavecchio, G. Ghisellini, and L. Maraschi for amusing and stimulating discussions and to S. Heinz for an interesting comment to my talk.

High Energy Astrophysical Phenomena (astro-ph.HE); Astrophysics of Galaxies (astro-ph.GA)


Parsec-Scale Jet-Environment Interactions in AGN

Abstract: Observations made with the VLBA have led to fundamental advances in our understanding of how radio jets in AGN evolve from parsec-scales out to distances exceeding several hundred kiloparsecs. In this review I discuss current models of young radio source evolution, as well as the observational evidence for a rapid change in jet properties on scales of ~1 kpc. A central topic of current debate is the relative importance of intermittent jet fueling versus jet-environment interactions in causing a drop-off in powerful radio sources at this critical evolutionary stage. Recent 3-D hydrodynamical jet simulations suggest that dense environments and cloud collisions can temporarily stifle, but not completely halt powerful relativistic jets. Several VLBA studies of jet-ISM interactions in both blazars and weak Seyfert jets have indicated that collimated outflows are indeed possible in dense environments. At present, the bulk of the evidence favors intermittent AGN accretion as the dominant factor in determining the evolutionary path of large numbers of AGN jets.

\section{Introduction}

One of the most fundamental questions that can be asked about jets
associated with active galactic nuclei (AGN) is how do they evolve from
their dense, gas-rich parsec-scale environments out to scales of
hundreds of kiloparsecs, well outside their host galaxies.  The
capability of radio wavelength interferometers to penetrate the dense
gas and dust in the centers of AGN host galaxies at high resolution
has brought us tantalizingly close to fully answering this
question. In this review, I briefly describe our current understanding
of young radio jet evolution, and the
relative role played by jet-environment interactions. I begin in
\S2 by discussing what has been learned from statistical
population studies, and devote Sections~3 and 4 to numerical jet
simulations and individual VLBA case studies that have improved our
understanding of interactions between AGN jets and their parsec-scale
environments.

\section{\label{evolution}Evolution of young AGN jets}

Our current knowledge of radio jet evolution owes a great deal to the
gigahertz-peaked spectrum (GPS) class of radio source, which comprise
approximately $\sim 10\%$ of flux-limited samples at
cm-wavelengths. Originally classified in early surveys as 'compact
doubles' by \cite{PM82}, subsequent improvements in VLBI capabilities
revealed weak central components, and in some cases faint bridges of
emission connecting them with bright outer features. It was soon
recognized that these AGN were miniature versions of the classical
kpc-scale lobe-core-lobe radio galaxies, with similar total radio
powers, but over a thousand times smaller in extent.

Based on observed size trends in GPS and compact steep spectrum (CSS)
sources (e.g., \citealt*{JS02}), self-similar expansion models (e.g.,
\citealt*{Beg96,BDO97}) were developed in which the overall linear extent
of the jets grow in proportion with their hotspot diameters. These
hotspots remain in ram pressure equilibrium with the external medium,
which implies that the evolution of the source is strongly dictated by
the density profile of the ISM. Numerical simulations (see \S~3) of
jets expanding into power-law external density profiles confirmed that
a large bow shock forms ahead of the hotspot, allowing the latter to
expand smoothly and propagate outward relatively unimpeded. Unlike the
dentist drill model for kpc-scale lobes, very little side-to-side
motion is expected for the pc-scale hotspot. Spectacular confirmation
of these models came with the first measurements of hotspot proper
motions in GPS radio galaxies (e.g., \citealt*{OCP99}), which
displayed predominantly outward (non-transverse) motion. The derived
kinematic ages, based on constant expansion, were typically $\sim
1000$ y \citep{GTP05}, confirming that these were in fact recently
launched jets.

The first problems with the standard scenario arose with detailed
studies of population statistics. In a steady-state population, one
would expect a rather flat distribution of kinematic ages, but in
fact, the observed one is peaked at young ages
\citep{GTP05}. A similar conclusion had been reached previously
by independent authors who considered the luminosity functions of GPS
sources (e.g., \citealt*{Beg96,RTP96}). Given their high luminosities,
the young radio sources were too numerous compared to their more aged
radio galaxy cousins, implying that must either dim rapidly, or die
out completely before reaching sizes of a few kpc. A lingering issue
of current debate is the relative importance of AGN fueling and
environmental interactions in dictating the evolution of radio jets at
this critical evolutionary stage.

\subsection{AGN fueling and intermittent jet activity}

Although a simple argument for intermittent jet activity in AGN can be
found in the fact that only $\sim10\%$ of all AGN associated with
super-massive black holes are radio loud, yet the lifetimes of
individual AGN are on the order of a few hundred Myr, true 'smoking
gun'-type evidence has become available only relatively recently. The
most compelling has been the discovery of the 'double-double' class of
radio galaxy \citep{SDR00}, of which roughly a dozen are currently
known \citep{MTM06}. These sources contain two sets of nested radio
lobes, which are symmetric with respect to a central component
associated with the active nucleus. The inner double resembles in many
ways a GPS source, with a peaked radio spectrum, bright hotspots, and
fast expansion speed. The outer lobe structures, on the other hand,
have sizes comparable to the those of the largest known radio
galaxies.  The notable gap in radio emission between the two
components is indicative of a long quiescent period, on the order of
$10^6-10^7$ Myr, in which the jet was presumably switched off (e.g.,
\citealt*{OKB01}).

Understanding intermittent jet activity in AGN is undoubtedly an
important factor in building a complete model of jet evolution (e.g.,
\citealt*{RB97}). However, it is still a nascent field in which the
necessary statistical samples (needed because of the long evolutionary
timescales involved) are still being gathered.  As I will describe
below, considerably larger progress has been made in understanding the
role played by jet-ISM interactions in affecting AGN jet evolution.

\subsection{Basic forms of jet-ISM interaction}

Because they are relatively light compared to their external
environments (density contrasts on the order of $10^{-3}$,
e.g., \citealt*{Krause03}), AGN jets are highly susceptible to external
interactions, which can be classified roughly into three main areas:
\begin{itemize}

\item      {\bf Bow shock-hotspot interaction} at the jet terminus, as
in the standard models described above.

\item       {\bf Cloud collisions}, which can cause bending and disruption of the flow.

\item        {\bf Entrainment}, leading to shear layers, deceleration,
instabilities, and possible particle acceleration at the jet boundaries.
\end{itemize}

Although much is known about the physics of entrainment in
kiloparsec-scale jets, progress on parsec-scales has been limited by
several factors. These include the difficulty of observing faint,
diffuse emission at the jet boundaries with limited dynamic-range
VLBI, as well as a paucity of bright, nearby AGN jets which we can
resolve in a transverse direction to the flow. Furthermore, studies of
the crucial 100-1000 milliarcsecond region where jets may undergo
strong internal changes due to entrainment have been hampered by the
lack of a suitable interferometer matching the sensitivity of the VLA
or VLBA. For these reasons I will concentrate hereafter on the issue
of jet interactions with dense clouds in the nuclear region of the
host galaxy.

\section{\label{sims}Numerical jet-cloud simulations}

Numerical simulations continue to play a vital role in understanding
the structure and evolution of AGN jets, by providing the ability to
test various scenarios under controlled conditions. Early numerical
jet-medium interaction studies were able to reproduce classical bow
shock and hotspot structures by propagating supersonic outflows into
external media with uniform density and pressure gradients (e.g.,
\citealt*{HN90}). The extension of MHD codes to the fully
three-dimensional, relativistic regime has made it possible to
robustly examine powerful jet evolution through a more realistic,
non-uniform medium for the first time. I describe here two such
studies (\citealt{CW07}, and \citealt*{SB07}), that are of particular
relevance to young jet evolution.

The simulations of \cite{CW07} employ a fully 3-D, pure hydrodynamic
code to simulate the passage of the relativistic jet through a
two-phase medium. The latter consists of a single dense cloud embedded
in a constant-pressure gas. They examined cases of both high ($\Gamma
= 7$) and low ($\Gamma = 2.29$) Lorentz factor jets striking the cloud
slightly off-axis. During the interaction, an oblique shock forms in
the jet, causing it to bend. Unlike previous non-relativistic studies
(e.g., \citealt*{wang00,higgins99}), the flow itself does not undergo
any significant deceleration or decollimation, and remains stable
after the interaction event. By varying the cloud-to-ambient medium
ratio, the authors find that the highest deflections occur in the case
of low-Mach number jets hitting denser clouds, with cloud density
being the dominant factor. Thicker clouds end up being less
encompassed by the bow shock, allowing earlier interaction with the
Mach disk and stronger oblique shocks in the flow. The clouds
themselves can actually survive the event, provided the cloud/jet
density contrast is high enough to suppress most Kelvin-Helmholtz
instabilities. These regions of shocked gas may be important star
formation sites (see \S~4) and may play a role in creating the
emission-line/jet alignment effect in AGN (e.g.,
\citealt*{MBS87}).

\cite{SB07} investigate the more general case of a
jet propagating through an inhomogeneous medium in the form of a
massive ($10^{10} \;\mathrm{M_{\sun}}$), turbulently supported disk
plus a hot ($10^7$ K) ISM. Like \cite{CW07}, they use a fully 3-D pure
hydrodynamic code, although in this case a non-relativistic one for
which they derive relativistic scaling parameters according to
\cite{KF96}.  In the initial phase of their simulations of a $\sim
10^{43}\; \mathrm {erg\; s^{-1}}$ jet, the morphology looks strikingly
different than those seen in other studies that assume a uniform ISM,
in that the flow attempts to seek out and pass through the
lowest-density locations in the clumpy (fractal) medium. In doing so,
multiple channels are formed and reformed, followed by the formation
of quasi-spherical bubbles around the jet and counter-jets that expand
outward. Making simple assumptions about the gas emissivity, the
authors find that these bubbles should be prominent in hard
X-rays. Once the jet reaches the outer edge of the disk and clears the
last obstruction, a stable, linear outflow develops, containing the
standard re-collimation and bow shock structures. At this point it
pierces the expanding bubble and evolves as in the uniform medium
case.

The authors find a good deal of similarity between the predicted radio
emission from their simulations and the compact symmetric object (CSO)
4C~31.04 \citep{CFG95}. This young radio source is characterized by a
large asymmetry in its jet and counter-jet structure, as well as lobe
spectral index gradients that are difficult to reconcile with standard
models of cocoon backflow \citep{giro03}. Comparison with their
simulations led \cite{SB07} to suggest that the western lobe may be
near the end of the breakout phase, whereas the eastern lobe is at a
slightly earlier stage of evolution. The strong apparent northward
deflection of the western lobe flow at the hotspot is also reminiscent
of structure found in the simulations of \cite{CW07}.

The conclusion that can be drawn from these studies is that powerful
relativistic jets are not likely to be permanently stifled by neither
direct jet-cloud collisions, nor a dense, clumpy external
medium. Instead, it is more likely that they all pass through an
evolutionary stage in which the flow may be bent and not necessarily
well-collimated. The duration of this stage is largely determined by
the power of the jet, and to a lesser extent, the jet/medium density
contrast. Given the good initial agreement with observed jet structure
from these preliminary simulations, it suggests that through careful
study of jet morphologies of young radio sources, it may be possible
to identify the precursors to both high- and low-power radio
galaxies, as well as to characterize their early evolutionary paths.

\section{\label{cases}VLBA studies of jet-environment interactions}

In addition to providing measurements of kinematic expansion speeds,
the VLBA provides a variety of unique tools for studying jet-medium
interactions on parsec scales. These include HI absorption
measurements, Faraday de-polarization and electric vector rotation
measurements at sub-milliarcsecond resolution levels. I discuss here
several recent VLBA studies of ISM interactions in weak Seyfert
jets, as well as in powerful blazars.

\subsection{Seyfert galaxies}

The relative proximity (15-20 Mpc) of Seyfert galaxies makes them
ideal targets for investigating jet-environment effects with the VLBA
at spatial resolutions approaching several thousand A.U. Given that
their jet powers are typically a factor of 100-1000 smaller than
radio-loud quasars (e.g., \citealt*{gold99}), they are much more
subject to entrainment and disruption (e.g.,
\citealt*{deyoung06}). Their sporadic accretion rate also offers the
chance to examine in detail the effects of central engine disruption
on jet structure.

\subsubsection{NGC 4151:}
The nearly-face on Seyfert 1.5 galaxy NGC 4151 has been the subject of
many intensive VLBI studies, due to its well-defined, two-sided,
$\sim100$ parsec-long radio jets, as well as the large quantity of
neutral gas in its nuclear region. HST imaging has revealed numerous
ionized gas clouds in an inner region that is extended about an axis
roughly aligned with the radio jets \citep{hutchings98, kaiser00}. The
spatial geometry of the narrow-line region suggests a thick molecular
torus aligned perpendicular to the jet, which is confirmed by
$\mathrm{H_2}$ measurements \citep{fernan99}. VLBA absorption data
have also provided evidence for an inner HI ring \citep{ulvestad98,
mundell03}. The radio spectral flattening and brightness enhancement
of the jet at this location led \cite{mundell03} to suggest that this
marks a site of jet-ISM interaction. Although the VLBA images lack
sufficient dynamic range to fully examine the extremely weak surface
brightness structure, the jet does undergo an abrupt deviation at this
point, in a manner similar to the jet-cloud simulations of
\cite{CW07}. \cite{mundell03} found the HI absorption line profiles to
vary significantly toward different portions of the jet, indicating a
medium composed of clumpy dense clouds with a variety of
velocities. Although they speculate that some of the other bright
knots in the jet may be the result of jet-cloud encounters, the
authors rule out shock ionization as the main source of the NLR, based
on its imprecise alignment with respect to the radio jet, and the
presence of several low-velocity clouds very near the jet that show no
signs of interaction.

\subsubsection{NGC 3079:}

This is another good example of a Seyfert jet in a dense environment,
albeit in this case the galaxy is viewed nearly edge-on
\citep{sosa01}. Using a series of VLBA measurements over a six year
period, \cite{middel07} have discovered complex kinematics and
variable jet emission in this source. They found one bright jet knot
initially moving at nearly 0.1 c, only to watch it decelerate and
become virtually stationary during the final year of their
observations. During this time its flux density increased and its
spectrum changed to a convex free-free/synchrotron-self absorbed
profile. This behavior is consistent with that expected from the
jet-cloud simulations described in \S~\ref{sims} Furthermore, the
source contains several steeper spectrum features well off the main
jet axis, which could perhaps be remnants of earlier flow channels
as predicted by \citealt*{SB07}. NGC 3079 thus provides an excellent example
of the potential of multi-epoch VLBA studies for exploring the
kinematics of jet-cloud interactions at exceedingly high spatial
resolution.

\subsubsection{PKS 1345+12:}

The ultra-luminous infrared galaxy IRAS 13451+1232 is a recent merger
system with significantly distorted optical morphology and a binary
nucleus, the northwest of which has been classified as a Seyfert 2
(e.g., \citealt*{scoville00}). The latter also contains a spectacular
radio jet (PKS 1345+12), which extends nearly 200 pc in a continuous,
sinusoidal pattern \citep{LKV03}. The counter-jet is also visible, but
only out to $\sim 50$ pc from the nucleus. Although these properties
are consistent with the CSO class, this object is unique in the fact
that \cite{LKV03} measured speeds of 1 c in the innermost
jet region, as well as high fractional polarization at the location of
the southern hotspot. The latter is significant as it implies a
continuous resupply of energy, i.e., the southern jet is not stifled
by this very gas rich galaxy.

By fitting to the apparent ridge line, apparent speeds, and
jet/counter-jet ratio, \cite{LKV03} concluded that the jet follows a
three-dimensional, conical helix aligned 82 degrees from our line of
sight, with an intrinsic flow speed of $\sim 0.8$ c. Similar sinusoidal
ridge lines seen in other CSOs and blazars have led various authors to
conclude that these may be the result of growing Kelvin-Helmholtz
instability modes, driven by small perturbations at the jet nozzle and
excited by interaction with the medium at the jet boundaries. The
northern counter-jet shows a deviation from the predicted best-fit
helical path, and is truncated at the site of dense HI absorption
($>10^{22} \;\mathrm{cm^{-2}}$; \citealt*{morganti05}). This appears
therefore to be a clear case where asymmetries in the external
environment have a strong differential impact on the morphology and
evolutionary rates of the jet and counter-jet of a young radio source.

\subsection{Blazar Jets}

Despite their much larger distances, blazar jets can also serve as
useful probes of parsec-scale jet interactions. First, because they
are viewed directly down the opening in the obscuring torus, there is
much less de-polarization, meaning that the jet polarization and
magnetic field properties can be directly studied. This also means
that any intervening gas can be potentially studied via Faraday
rotation measures (e.g., \citealt*{ZT05}). Second, any slight
deviations in the flow that may be caused by interactions are greatly
magnified by projection effects. Finally, because of Doppler effects,
there are many examples of blazars where over a century of jet evolution is
compressed into a span of only a few years of observing time (e.g.,
\citealt*{KL04}).

\subsubsection{3C 279:}

The powerful jet in the quasar 3C279 was one of the first jets in
which superluminal motion was witnessed, and has been the target of
intensive study in a variety of wave-bands. The jet has been
regularly imaged since 1994 by the 2 cm Survey \citep{KL04} and MOJAVE
\citep{LH05} programs with the VLBA at a wavelength of 2 cm. Shorter
wavelength (7 mm) VLBA monitoring \citep{jorstad04,jorstad05} has
revealed a regular swing in the ejection direction of the jet close to
the nozzle, over a timescale of 3 years. \cite{homan03} describe one
prominent jet feature (C4) that was ejected in late 1984, which moved
steadily along a linear path for over a decade with an apparent speed
of 8 c, before suddenly undergoing an increase in brightness and
change in polarization angle in 1998. These events were followed
shortly thereafter by a rapid apparent acceleration to 13 c, and a
change in trajectory of 26 degrees. Under the most conservative
assumptions, \cite{homan03} found that these changes were consistent
with an intrinsic bend of only 0.5 to 1 degree. Given the fact the
brightening and polarization changed {\it before} the change in
trajectory, the most plausible scenario is one in which C4 is
interacting with the external environment. Furthermore, the direction
of the new trajectory closely matches that of another feature ejected
several decades previously, which rules out a random jet-cloud
collision. The authors suggest instead that the event represents a
collimation of the jet resulting from a jet-boundary interaction at a
de-projected distance $\mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$}}}\hbox{$>$}}} 1$ kpc from the nucleus. Since this is
the first such an event to be witnessed in an AGN jet, it is difficult
to yet draw solid conclusions on the validity of this model. However,
large intensive VLBA monitoring programs such as MOJAVE \citep{LH05}
may soon provide additional examples for further study.

\subsubsection{3C 120:}

Although classified as a Seyfert 1, this nearby (z = 0.033)
broad-lined galaxy shares many properties with blazars, including
superluminal motions of up to 6 c, a one-sided radio jet, and flux
variability. \cite{axon89} found high-velocity emission line
components in the host galaxy that suggested interaction between the
jet and gas clouds in the NLR. The excellent spatial resolution (0.1
pc) achievable by the VLBA at 43 GHz has enabled detailed study of its
jet evolution in both total intensity and linear polarization
\citep{gomez01, jorstad05}. The jet is resolved perpendicular to the
flow direction, and a distinct asymmetry is seen between the northern
and southern edges. In particular, \cite{gomez01} have found a
distinct region in the southern edge, approximately 8 pc
(de-projected) from the base of the jet, where moving jet features
show marked changes as they pass through. These include a brightening
in flux density, and a rotation of their polarization electric vector
position angles (EVPAs). These events are different from that
witnessed in 3C 279, since in this case no accelerations are
seen. \cite{gomez01} conclude the most likely explanation to be
interaction with a cloud, which causes Faraday rotation of the EVPAs,
and shocking of the jet material. There is also an indication of a
slight bend at the interaction site, although the jet remains
well-collimated downstream. Ideally it would be useful to study
additional examples of this type of interaction, but unfortunately
there are still very few known bright jets that are close enough to be
resolved transversely by the VLBA, and yet have viewing angles small
enough not to be heavily de-polarized by foreground nuclear gas.

\section{Summary}

High-resolution radio observations of young radio jets associated with
gigahertz-peaked spectrum AGN have led to considerable insight into
the evolutionary processes of AGN jets. Kinematic and population
studies have shown that these young radio sources undergo a
significant decline in numbers when they reach sizes of $\sim 1$
kpc. VLBA studies of individual jets have provided clear evidence for
interaction with clouds in their external environment, suggesting
stifling by dense gas as a possible cause. However, detailed
numerical simulations of jet-environment interactions indicate that
dense, clumpy environments can only temporarily stifle the flow of
powerful jets, even in the case of direct jet-cloud
collisions. Furthermore, the discovery of 'double-double' galaxies has
provided solid evidence of recurrent jet activity in powerful AGN. It
therefore appears likely that variable accretion rates play a major
role in determining the evolutionary paths of many AGN. The enhanced
resolution and sensitivity of upcoming facilities such as VSOP-II, the
EVLA, and the SKA should provide many new opportunities for studying
the evolution of young radio sources and their interactions with their
external environment.