Spin singlet pairing in the superconducting state of NaxCoO2\cdot1.3H2O: evidence from a ^{59}Co Knight shift in a single crystal

Abstract: We report a ^{59}Co Knight shift measurement in a single crystal of the cobalt oxide superconductor Na_{x}CoO_2\cdot1.3H_2O (T_c=4.25 K). We find that the shift due to the spin susceptibility, K^s, is substantially large and anisotropic, with the spin shift along the a-axis K^s_a being two times that along the c-axis K^s_c. The shift decreases with decreasing temperature (T) down to T\sim100 K, then becomes a constant until superconductivity sets in. Both K^s_a and K^s_c decrease below T_c. Our results indicate unambiguously that the electron pairing in the superconducting state is in the spin singlet form.

\section{Introduction}

The recently discovered superconductivity in cobalt oxide Na$_{x}$CoO$_2$$ \cdot$1.3H$_2$O has attracted considerable attention \cite{Takada}. The new material consists of two-dimensional CoO$_2$ layers that bear similarity to high temperature (high-$T_c$) copper oxide superconductors. A simple but wide-spread view of the electronic state is that increasing the sodium content $x$ adds $x$ electrons to the
spin 1/2 network in the CoO$_2$ layer \cite{Baskaran,WangLeeLee}, which is similar to the case of doping the spin 1/2 network in the CuO$_2$ layer of the high-$T_c$ cuprates. However, there is a large difference between these two classes of transitional-metal oxides. Namely, Co forms a triangular lattice in the cobaltate, instead of square lattice in the cuprates. Because of this crystal structure, exotic superconducting states are expected in the cobaltate. So far, chiral $d$-wave, $p$-wave,  $f$-wave and  $i$-wave  states \cite{Baskaran,WangLeeLee,TanakaHu,Ogata,Lee,Ikeda,Kuroki0,Kuroki}, as well as odd-frequency $s$-wave state \cite{Johanes} have been proposed theoretically. The $d$- and $i$-wave states are spin singlet states, while the remainder are spin triplet states.

Indeed, the nuclear spin-lattice relaxation rate ($1/T_1$) measurements by us indicate that the superconducting state is unconventional; $1/T_1$ is proportional to $T^3$ below the transition temperature $T_c$, suggesting line nodes in the gap function \cite{Fujimoto, Zheng}. However, the symmetry of the spin pairing has been unclear. Previous  measurements of the angle-averaged Knight shift in powder samples have resulted in   controversial conclusions \cite{Kato,Kobayashi}. Kato {\it et al} reported no change of the shift below $T_c$  \cite{Kato}, while Kobayashi {\it et al}  found a decrease of the shift in a powder sample \cite{Kobayashi}. The former suggested  triplet pairing, while the latter concluded to $s$-wave pairing. The difficulty in growing the samples,  which often accompany a  secondary phase that gives rise to a spurious NMR peak overlapping  the main peak in question,  complicates the data analysis in power samples. Therefore, measurements of the spin susceptibility via the Knight shift in single crystals are required.

In particular, the Knight shift in the $c$-axis direction is essential to distinguish various pairing states proposed theoretically  \cite{Baskaran,WangLeeLee,TanakaHu,Ogata,Lee,Ikeda,Kuroki0,Kuroki,Johanes}. This is because even for triplet pairing, the spin susceptibility in the $a$-axis direction can be reduced below $T_c$ but that along the $c$-axis direction should remain unchanged across $T_c$ \cite{Yanase}. Kobayashi {\it et al} have attempted to align small single crystal platelets \cite{Kobayashi2}, but their samples contain multiple phases with almost equal amount of non-superconducting composition that contaminated the NMR signal observed.
 So far,  no work has been reported to determine the absolute value of the $c$-axis shift due to the spin susceptibility which is important in discussing the superconducting state properties.

In this paper, we report the first complete data set of the $^{59}$Co Knight shift  in a high quality {\it single crystal} with $T_c$=4.25 K. We find that the spin contribution to the Knight shift is significantly large and anisotropic. The spin Knight shift along both the $a$-axis and $c$-axis decrease below $T_c$. Our results indicate unambiguously that the spins that form the Cooper pairs  in the superconducting state are anti-paralleled.

Single crystal of Na$_{0.42}$CoO$_2$$\cdot$1.3H$_2$O used in this study was grown by the traveling solvent floating zone (TSFZ) method, as described in a previous publication \cite{Chen}.  AC susceptibility measurement using the in-situ NMR coil indicates that the $T_c$ at zero magnetic field is 4.25 K (data not shown).
NMR experiments were performed using a home-built phase-coherent spectrometer.  The NMR spectra were taken by changing the external magnetic field ($H$) at a fixed RF frequency and recording the echo intensity step by step.

Figure 1 shows a typical example of the NMR spectra for $H\parallel$c-axis and $H\parallel$a-axis. As can been seen in Fig. 1(a), a sharp central transition line accompanied by six satellites due to the nuclear quadrupole interaction, are observed.
There is essentially no signal intensity between the sharp peaks, which indicates good quality of the single crystal.

\begin{figure}
\begin{center}
\includegraphics[scale=0.4]{Fig1.eps}
\caption{(Color online) $^{59}$Co NMR spectra in Na$_{0.42}$CoO$_2$$\cdot$1.3H$_2$O at $T$=4.2 K with $H\parallel$c-axis at 14.1 MHz (upper panel) and with $H\parallel$a-axis at 59.4 MHz (lower panel).}
\label{fig:1}
\end{center}
\end{figure}

Figure 2 shows the temperature dependence of the Knight shift $K_a$ with the field applied along the $a$-axis, and $K_c$ with the field applied along the $c$-axis.
$K_c$ was determined from the central peak as well as from the mid-point between the two first satellites; results obtained by the two procedures agree well \cite{note0}. $K_a$ was determined from the central transition, by taking into account the shift due to quadrupolar interaction effect \cite{Abragam}.
Upon decreasing the temperature from 260 K, both  $K_a$ and $K_c$ decrease, but tend to be saturated below $T$=100 K before superconductivity sets in.
The decrease of the Knight shift with decreasing temperature is consistent with the bulk susceptibility results reported by Foo {\it et al} \cite{Foo} and also for the  sample used in the present study \cite{Alloul}, although a
previous measurement  has failed to detect the temperature dependence \cite{Imai}.
Below, we  show that the spin shift is quite large and can be used to distinguish the pairing symmetry in the superconducting state.

\begin{figure}
\begin{center}
\includegraphics[scale=0.4]{Fig2.eps}
\caption{(Color online) Temperature dependence of the Knight shift $K_a$ with the field applied along the $a$-axis,  and $K_c$ with the field applied along the $c$-axis.}
\label{fig:2}
\end{center}
\end{figure}

The Knight shift is composed of the spin part $K^s$ due to the spin susceptibility ($\chi^s$) and the orbital  part $K^{orb}$ due to the orbital susceptibility ($\chi^{orb}$) which is $T$-independent.
\begin{eqnarray}
K = K^{orb}+K^s
\end{eqnarray}
\begin{eqnarray}
K^s=A^{hf}\chi^s
\end{eqnarray}
where $A^{hf}$ is the hyperfine coupling constant between the nuclear and the electron spins. Figure 3 shows the plot of $K_c$ versus $K_a$. The straight line is a fitting to the data which gives $K_c$=0.5$K_a$+const., which means that $K^s_{a}=2K^{s}_{c}$.
The temperature dependence of the Knight shift  is similar to the case of electron-doped cuprate superconductor (Pr$_{0.91}$La$_{0.09}$)$_2$CeCuO$_4$, in which $K_a$ also decreases with decreasing temperature and becomes a constant below $T$=60 K \cite{Zheng1}.

\begin{figure}
\begin{center}
\includegraphics[scale=0.4]{Fig3.eps}
\caption{(Color online) The plot of $K_c$ vs $K_a$. The straight line is the best fit to the data which gives $K^s_{a}=2K^{s}_{c}$. The arrows indicate the estimated orbital contributions to the shifts (see text).}
\label{fig:3}
\end{center}
\end{figure}

Let us focus  on the temperature dependence of the Knight shift below $T_c$. First, both $K_a$ and $K_c$ decrease below $T_c(H)$, as can be seen in
Figure 4 that shows the enlarged part of the shift below $T$=5 K for $H$(3.82 T)$\parallel a$-axis and $H$(0.93 T)$\parallel c$-axis. $T_c(H\parallel a$)  and $T_c(H\parallel c$) at which the shift decreases is 3.5 K and 3.1 K, respectively, which is in good agreement with the report by Chou {\it et al} \cite{Chou}. Figure 5 compares two of the spectra above and below $T_c$, respectively, for both field configurations.
The reduction of $K_a$ and $K_c$ at $T$=1.6 K amounts to 0.2\% and 0.13\%, respectively. The smaller reduction of $K_c$ is mostly due to  the smaller spin part $K^s_c$.
For $H\parallel c$-axis, we have also carried out measurements at 14.1 MHz ($H\sim$1.38 T), and found detectable shift change between $T_c (H)$ and $T$=1.6 K.  Namely, $\Delta K_c$=$K$ (2.8 K) - $K$ (1.6 K) = 0.05\%. We have further measured the shift at  7.1 MHz ($H\sim$0.69 T) in another sample with $T_c$=3.6 K,  and found $\Delta K_c$=$K$ (3.0 K) - $K$ (1.6 K) = 0.15\%. In all cases,
we have confirmed that the   $^{23}$Na NMR spectrum  has a full width at half maximum (FWHM) of 20 kHz, and does not move appreciablly across $T_c$ within the  experimental accuracy. Assuming that one can detect a change of one tenth  the FWHM in the $^{23}$Na spectrum, we estimate that the contribution to the Knight shift due to the diamagnetism effect ($K^{dia}$) in the vortex state is no larger than 0.02\%  \cite{note}.
Therefore the reduction of  $K_a$ and $K_c$ below $T_c(H)$ is due predominantly to the reduction of the spin susceptibility. This result indicates that the spin pairing in the superconducting state is in the {\it singlet} state.

\begin{figure}
\begin{center}
\includegraphics[scale=0.4]{Fig4.eps}
\caption{(Color online) The enlarged part of the $^{59}$Co Knight shift at low temperatures. The vertical arrows indicate $T_c(H)$ at the respective external fields. The horizontal arrows indicated $K_{a,c}^{orb}$, the shift due to orbital susceptibility (see text for the estimation).}
\label{fig:4}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[scale=0.4]{Fig5.eps}
\caption{(Color online) Comparison of the spectra above $T_c$ (broken curves) and below $T_c$ (solid curves) at different field configurations.}
\label{fig:5}
\end{center}
\end{figure}

Among the existing theoretical models, the $p$-wave \cite{TanakaHu}, $f$-wave \cite{Ogata,Lee,Ikeda,Kuroki0} states and the odd-frequency $s$-wave state \cite{Johanes} are spin triplet states, while the chiral $d$-wave state \cite{Baskaran,WangLeeLee} is a spin singlet but fully gapped state. They are not in  agreement with our finding, but  probably the $i$-wave state \cite{Kuroki} is an exception. However, a pure $d$-wave state with line nodes in the gap function appears to be consistent with both the Knight shift data and the $T^3$ dependence of $1/T_1$ \cite{Zheng}. Generally speaking, a state with lower angular momentum ({\it i.e.} $d$-wave) is energetically more favorable compared to its higher angular momentum counterparts ({\it i.e.} $i$-wave).  In addition, there are two conditions pertinent to the present case that favor a pure $d$-wave state. Firstly, although a fully-gapped $d_1+id_2$ state is predicted for an exact triangular lattice \cite{Baskaran,WangLeeLee},
anisotropy (10$\sim$20\%) of the lattice will stabilize a pure $d$-wave state \cite{Watanabe}. Secondly, it has been found that the spin order in Na$_{0.5}$CoO$_2$ is antiferromagnetic, with alternative rows of Co$^{4+}$ (s=1/2) and Co$^{3+}$ (s=0) in the CoO$_2$ plane \cite{Neutron}. In other words, the triangles with regards to the spin configuration are effectively squared. In fact, both  the Knight shift found in this study and the temperature dependence of $1/T_1T$ in the normal state in Na$_{x}$CoO$_2$$\cdot$1.3H$_2$O ($x=0.25\sim0.35$) \cite{Zheng} point to the presence of antiferromagnetic correlations, which would favor $d$-wave pairing.

Next we attempt to separate the $K^{s}$ and $K^{orb}$ from the observed shift data, and argue that the data below $T_c$ can be consistently accounted for by an anisotropic superconducting state with  nodes in the gap function and with spin singlet pairing.
In $s$-wave superconductors with isotropic gap, the spin susceptibility outside the vortex cores vanishes completely at a temperature far below $T_c$. In anisotropic superconductors with nodes in the gap function, however, the spin susceptibility outside the vortex cores can remain a finite value since the quasiparticles state can be extended from the vortex core along the nodal directions \cite{Volovic}. Volovik pointed out that the density of such delocalized quasiparticle state, $N_V$, is proportional to $\sqrt{H/H_{c2}}$ where $H_{c2}$ is the upper critical field  \cite{Volovic}. Indeed, in an overdoped cuprate superconductor, it has been found that  $N_V/N_0\sim0.7\sqrt{H/H_{c2}}$ (ref.\cite{Zheng2}), where $N_0$ is the density of states (DOS) in the normal state. Note that the coefficient 0.7 is less than 1, because NMR only probes the DOS outside the vortex cores. If we use the data of $H_{c2}^a$(0)=8 T and $H_{c2}^c$(0)=4 T by Chou {\it et al} \cite{Chou} and the same $N_V$  as in the cuprate, we obtain $K^{orb}_a$ =2.94\% from the relation $0.7\sqrt{H/H_{c2}}(K_a$(4.2 K)-$K^{orb}_a )=K_a$(1.6 K)-$K^{orb}_a$. Once $K^{orb}_a$ is known, one obtains readily $K^{orb}_c$ =2.02\% from the $K_c$ vs $K_a$ plot in Fig. 3. Such estimate  consistently explains the $K_c$ result; the $H$-induced shift at the lowest temperature due to $N_V(H\parallel c)$ is calculated to be 0.05\%, which is in fair agreement with the experimental result of 0.04$\pm$0.03\%. Note that  $K_a^s\geq$0.72\%, which is about two times the spin shift in cuprate superconductors \cite{Zheng1,Slichter}
Finally, we comment that the estimated $K_{a,c}^{orb}$ agrees fairly well with that independently  estimated from the  Knight shift vs DC susceptibility ($K vs \chi$) plot \cite{MatanoZheng}.

In conclusion,  the $^{59}$Co NMR measurements in a single crystal  of Na$_{x}$CoO$_2$$\cdot$1.3H$_2$O ($T_c$=4.25 K) allowed us, for the first time, to obtain a complete data set of the Knight shift. We found that the shift due to the spin susceptibility is substantially large (0.72\% along the $a$-axis) and anisotropic, with $K^s_a=2K^s_c$. Both $K^s_a$ and $K^s_c$ decrease below $T_c$,  indicating that the spin pairing of the Cooper pairs is in the singlet state. Our work that  settled  the electron pairing symmetry in the superconducting state should help pave the first step toward understanding the mechanism of the superconductivity  in  this new class of materials.

We thank M. Nishiyama for assistance in some of the measurements. This work was  partially supported by  MEXT grants for scientific research.

Superconductivity (cond-mat.supr-con); Strongly Correlated Electrons (cond-mat.str-el)


A stellar endgame -- the born-again Sakurai's object

Abstract: The surface chemical composition of this remarkable star shows that it is hydrogen-deficient, carbon-rich and enriched in the light s-process elements. Spectra taken in May and October 1996 indicate a decrease in the surface hydrogen abundance by 0.7 dex in five months along with an increase in the abundances of Li, Sr, Y and Zr. The abundance changes are in agreement with the hypothesis of the star being a rapidly evolving ``born-again'' AGB star experiencing a final He-shell flash, similar to FG Sge. The C^12/C^13 ratio in October is very low, also suggesting hydrogen ingestion. By chemical composition, Sakurai's object resembles the R Coronae Borealis (R CrB) stars.

\section{Introduction}

The nova-like brightening of Sakurai's object in Sagittarius  has been
attributed to
the  final helium-shell flash of a central star of a planetary
nebula, that returns the star towards the domain of red giants in the
Hertzsprung-Russell diagram
(Duerbeck \& Benetti 1996).
Few  stars have been identified with
this phase of evolution:
examples include
FG Sge, V605 Aql, the planetary nebulae Abell 30, Abell 78 and N66.
It is expected that a
born-again red giant will consume hydrogen and become starkly
hydrogen-deficient, helium- and carbon-rich.
Low-resolution spectra led Duerbeck and Benetti (1996) to
suggest that Sakurai's object is hydrogen-poor.
The presence of strong lines
of neutral carbon and oxygen was also noted.

Changes of the surface chemical composition may be rapid for born-again
AGB stars, as was observed for FG Sge. Sakurai's object offers the
prospect of monitoring such secular changes in another born-again candidate.
Observations, as reported here, are surely crucial for an improved understanding of the
final He-shell flash.

\section{Observations}

Spectra covering 3700-10150\,\AA $ $ at a resolution of about 30,000
were obtained with
the 2.7\,m telescope at McDonald Observatory on May 5 and 6
and on October 7,  1996.
A spectrum was also obtained with the 2.1\,m telescope:
this spectrum from May 9, 1996, covers the region
5720\,\AA\ to 7200\,\AA\ at a resolution of about 60,000.

\section{Chemical composition}

Our analysis is based on
line-blanketed, hydrogen-deficient
model atmospheres, similar to those
described by Asplund et al. (1997) but with a range of hydrogen
abundances.
In estimating the stellar parameters T$_{\rm eff}$, log $g$ and hydrogen
abundance  various ionization
(Fe\,{\sc i}/Fe\,{\sc ii}, Mg\,{\sc i}/Mg\,{\sc ii},
Si\,{\sc i}/Si\,{\sc ii},
Cr\,{\sc i}/Cr\,{\sc ii})
and excitation equilibria
([O\,{\sc i}]/O\,{\sc i}, Fe\,{\sc i}, Fe\,{\sc ii})
together with the H$\beta$ and H$\alpha$
line profiles (with line broadening data following Seaton 1990)
have been used.
The C/He ratio was determined from the C\,{\sc ii} and He\,{\sc i} lines
in the May spectra, which indicate C/He $\simeq 10$ \%.
The same ratio had to be assumed for October when
the
lines were too weak to be utilized.
The microturbulence parameter
was estimated from Ti\,{\sc ii}, Fe\,{\sc i} and Fe\,{\sc ii}
lines of different strengths.
The May spectra are characterized by $T_{\rm eff} = 7500\pm 300$\,K,
log\,$g = 0.0\pm 0.3$, and
$\xi_{\rm t} = 8.0\pm 1.0$\,km\,s$^{-1}$, while it had cooled significantly
in October: $T_{\rm eff} = 6900$\,K, log\,$g = 0.5$, and
$\xi_{\rm t} = 6.5$\,km\,s$^{-1}$.
In fact, the derived parameters are not consistent with a constant stellar
luminosity but rather indicate a decrease by a factor of 4,
which is not supported by the observed photometry.
It could, however,
be that hydrostatic equilibrium is
inapplicable in May due to an expansion of the star or effects of
turbulent pressure:
a dynamical atmosphere can be mimicked by an underestimate of log\,$g$
when assuming hydrostatic equilibrium.
Indeed, with the May parameters the star is located at the
classical Eddington limit
(e.g. Asplund \& Gustafsson 1996).

The analysis of the C\,{\sc i} lines reveals the same inconsistency between
theoretical and observed line strengths as for R\,CrB stars (Gustafsson \& Asplund 1996;
Lambert et al., in preparation): the strengths of weak lines predicted
with the input C abundance are a factor of 4 stronger than observed
(Fig. \ref{f:spectra} and \ref{f:hbeta}).
It should be noted that no agreement
between all $T_{\rm eff}$-log\,$g$ indicators
could be achieved using consistent C abundance for the analysis.
Naturally, this C\,{\sc i} problem makes the absolute abundances uncertain
but relative abundances
are generally expected to be much less affected (Lambert et al., in preparation).

\begin{figure}[t]
\centerline{
\psfig{figure=figure1.ps,height=6.cm}}
\caption{A selected piece of spectrum in May (solid) and October (dashed)
showing the increase of some elements, e.g. Sc, Ti, and Y. The dotted
curve is the synthetic spectrum with the stellar parameters of
October but with the May abundances.
Note also that all predicted C\,{\sc i} lines are too strong
}
         \label{f:spectra}
\end{figure}

\begin{figure}[t]
\centerline{
\psfig{figure=figure2.ps,height=9.8cm}}
\caption{{\bf a} H$\beta$ in October (thick solid)
compared with predicted
line profiles for solar H (dotted) and H-deficient by
3.0\,dex (dashed).
{\bf b} C$_2$ (1-0) Swan band for $^{12}$C/$^{13}$C = 2 (dotted),
5 (dashed) and 10 (dash-dotted), together with the
observed October spectrum
(thick solid).
Also shown in both figures are the May spectra (solid) but
displaced upwards by 0.2 for clarity
}
         \label{f:hbeta}
\end{figure}

The derived LTE
abundances for May and October are summarized in Table \ref{t:abund}.
More details on the analysis and atomic data (lines, gf, hfs, etc) as
well as a comparison with V854\,Cen
will be given elsewhere.
The weak Balmer lines certainly
rule out a solar hydrogen abundance (Fig. \ref{f:hbeta}).
Note that the absolute abundances of most elements
are effectively unchanged
from May to October within the uncertainties (typically $\leq 0.3\,$dex).
Some elements, however, exhibit a marked change, for example, hydrogen declined
as lithium and the light $s$-process elements increased
in abundance
by a factor of about 4 (Fig. \ref{f:spectra}).
Also Sc, Ti, Cr and Zn seem to have increased during the time\-span
(Fig. \ref{f:spectra}).
The general agreement between the May and October abundances for most
elements suggests that the stellar parameters are not seriously in
error, which could
otherwise
have resulted in spurious abundance
effects.
Besides Li, the abundances of
elements showing variations are not very sensitive
to the stellar parameters:
the required $\Delta T_{\rm eff} \approx 1000$\,K for
either May or October to annul the abundance variations
would be inconsistent with the $T_{\rm eff}$--log\,$g$ indicators
and introduce other as severe changes (e.g. for Ca)
less easily explainable;
a different log\,$g$ can not simultaneously explain
all changes.
It would also only aggravate the luminosity discrepancy.
Hence, the few
changes seem to be real. Furthermore, they are
limited to elements expected to show alterations due
to a final flash.

\begin{table}[t]
\caption{ Chemical compositions of Sakurai's object, the
R\,CrB stars and the Sun
(normalized to log\,($\Sigma \mu_i \epsilon_i$) = 12.15)
\label{t:abund}
}
\begin{tabular}{lccccc}
 \hline
  Element & Sun$^{\rm a}$ & \multicolumn {2} {c} {Sakurai's object} &
\multicolumn {2} {c} {R\,CrB$^{\rm b}$}   \\
 && May & October & majority & minority  \\
\hline \\
H  & 12.0 &  9.7 & 9.0 &   &   $<4.1 - 10.8$   \\
He & 11.0 & 11.4$^{\rm c}$ & 11.4$^{\rm c}$ & 11.5$^{\rm c}$ & 11.5$^{\rm c}$   \\
Li &  3.3$^{\rm a}$ &  3.6 & 4.2 &       &   \\
C  &  8.6 &  9.7$^{\rm d}$ & 9.8$^{\rm d}$ & 8.9$^{\rm d}$ & 8.6 -- 9.5$^{\rm d}$ \\

N  &  8.0 &  8.9 &  8.9 &  8.6 &  7.6 -- 8.6 \\
O  &  8.9 &  9.5 &  9.4 &  8.2 &  7.5 -- 8.8 \\
Ne &  8.1 &  9.3 &      &      &  7.9 -- 9.6 \\
Na &  6.3 &  6.7 &  6.8 &  6.1 &  5.8 -- 5.9 \\
Mg &  7.6 &  6.6 &  6.5 &      &  6.1 -- 7.3 \\
Al &  6.5 &  6.6 &  6.3 &  6.0 &  5.3 -- 5.6 \\
Si &  7.5 &  7.1 &  7.5 &  7.1 &  7.3 -- 8.1 \\
S  &  7.3 &  6.6 &  6.9 &  6.9 &  6.7 -- 7.6 \\
K  &  5.1 &  4.8 &  5.0 &      &            \\
Ca &  6.4 &  5.6 &  5.5 &  5.4 &  5.0 -- 5.3 \\
Sc &  3.2 &  3.1 &  3.9 &      &            \\
Ti &  5.0 &  4.1 &  4.6 &      &            \\
Cr &  5.7 &  4.5 &  5.1 &      &            \\
Fe &  7.5 &  6.3 &  6.6 &  6.5 &  5.0 -- 5.8 \\
Ni &  6.2 &  6.1 &  6.2 &  5.9 &  5.2 -- 5.8 \\
Cu &  4.2 &  4.9 &  5.0 &      &            \\
Zn &  4.6 &  4.7 &  5.4 &  4.3 &  3.8 -- 4.1 \\
Rb &  2.6 & $<3.7$ &  4.6 &      &            \\
Sr &  3.0 &  4.9 &  5.4: &      &            \\
Y  &  2.2 &  3.3 &  4.2 &  2.1 &  0.6 -- 2.8 \\
Zr &  2.6 &  3.0 &  3.5 &      &            \\
Ba &  2.1 &  1.5 &  1.9 &  1.6 &  0.7 -- 1.3 \\
La &  1.2 & $<1.6$ &  1.5 &      &  \\
\hline
\end{tabular}

\begin{list}{}{}
\item[$^{\rm a}$] From Grevesse et al. (1996). For Li the meteoritic
value is adopted.
\item[$^{\rm b}$] From Rao \& Lambert (1996) and Jeffery \& Heber (1993).
The majority is an average of 14 stars while
the minority consists of V\,CrA, VZ\,Sgr, V3795\,Sgr and DY\,Cen.
\item[$^{\rm c}$] Input C/He ratio for model atmospheres: C/He=1\% assumed for
R\,CrB stars and 10\% estimated for Sakurai's object from the 1996 May spectra.
\item[$^{\rm d}$] Spectroscopically determined C\,{\sc i}
abundance, see text.
\end{list}

\end{table}

The metallicity of Sakurai's object is, judging from Fe, slightly below
solar by 0.2\,dex in mass fraction (0.9\,dex if the input rather than
the spectroscopic C abundance is adopted).
The quantities [Si/Fe], [S/Fe], [Ca/Fe],
and [Ti/Fe], which are 0.8, 0.6, 0.3 and 0.4
respectively, are, if unchanged
from the star's birth, also
indicative of a metal-poor star (Edvardsson et al. 1993).
An isotopic ratio $1.5\leq^{12}$C/$^{13}$C$ \leq 5$ is determined from
the strong C$_2$ (1-0) and (0-1) Swan bands (Fig. \ref{f:hbeta}).
The strengthening of the C$_2$ bands due to the change
in stellar parameters is clearly illustrated in Fig. \ref{f:hbeta}.

It is of considerable interest to compare the compositions of Sakurai's object
and FG\,Sge, another born-again candidate which has
recently experienced R\,CrB-like visual declines. FG\,Sge resembles Sakurai's
object in that it is
strongly $s$-element enriched (Langer et al. 1974), as well as
carbon-rich and poor in iron-group elements, except for Sc
(Kipper \& Kipper 1993).  In FG Sge, however,
the heavy $s$-elements are as
overabundant as the light, and it has not yet been shown to be
hydrogen-deficient.
FG\,Sge may therefore have experienced a late shell flash as a luminous post-AGB
star rather than a final flash as a white dwarf
(Bl\"ocker \& Sch\"onberner 1996).

Two of the outstanding aspects of the chemical composition of Sakurai's
object are hallmarks of the R\,CrBs: H-deficiency and a high C content,
but also other similarities in relative abundances exist
(Lambert et al., in preparation; Rao \& Lambert 1996;
Lambert \& Rao 1994).
Except for the high Y/Fe other observed X/Fe ratios are similar
to those found in R\,CrB stars.
In particular it resembles the (relatively) H-rich V854\,Cen
(Asplund et al., in preparation).
If, however, C/He is correctly estimated, it may sooner be related
to objects such as V605\,Aql, Abell 30 and 78 and the hot
R\,CrB star V348\,Sgr, which are also surrounded by planetary
nebulae and have been proposed to be final flash candidates
(Renzini 1990).

Similar abundance patterns as presented here for Sakurai's object
have also been obtained in less detailed analyses by
Shetrone \& Keane (1997) and Kipper \& Klochkova (1997).
Shetrone \& Keane's finding of a
near normal H abundance is however very puzzling.

\section{Abundance variations and nucleosynthesis}

In broad terms, the composition of Sakurai's object shows evidence of
severe contamination by material exposed to hydrogen and helium burning
and associated nuclear reactions. Close examination provides
some interesting
constraints on the nucleosynthesis experienced by the star.

The present atmosphere is not
a simple mix of initial unprocessed gas, gas
run through the H-burning CNO-cycles, and  H-exhausted gas exposed to He-burning,
but must have been accompanied by further processing.
This is
demonstrated by the low observed $^{12}$C/$^{13}$C ratio, which
encompasses the equilibrium value of 3.5 for CNO-cycling.
As the equilibrium abundance of $^{13}$C is very low following He-burning,
the observed ratio suggests that $^{12}$C from He-burning has been
exposed to hot protons. It would seem that
C-rich material from He-burning has been mixed with
ingested hydrogen such that
the proton supply is effectively exhausted in converting
inhibited (see Renzini 1990).

Not all protons are consumed in He-rich regions.
Production of lithium is ascribable to the
Cameron-Fowler (1971) mechanism.
Here  $^3$He synthesised in a low mass main sequence star is
converted to $^7$Li
in an  envelope that convects $^7$Li to low temperatures
where it survives
until re-exposed to high temperatures.
Production of lithium implies
H-burning in regions not previously exposed to
H-burning temperatures; $^3$He
which is destroyed in regions that have undergone H-shell or H-core burning
can hardly be resynthesised.
The observed Li is not a fossil from an
earlier stage as a  Li-rich AGB star:
the predicted Li/H ratio for
AGB stars
which have undergone hot-bottom burning
is 10$^{-8}$ while the
observed ratio
is 10$^{-5}$ to 10$^{-6}$ and hydrogen consumption
necessarily destroys fossil lithium.
The overabundant Na and Al have likely been synthesised through
$^{22}$Ne(p,$\gamma$)$^{23}$Na and $^{25}$Mg(p,$\gamma$)$^{26}$Al.
As in other H-deficient stars Ne is very high,
which can not be explained by $\alpha$-captures on N from initial
CNO, but must be due to products of He-burning and possibly
additional CNO-cycling.

A remarkable feature of Sakurai's object is the large overabundance of
light $s$-process elements and the high ratio of light to heavy $s$-process
elements. Probably, $^{13}$C($\alpha,n)^{16}$O is the neutron source.
The $s$-processing may be characterized by the neutron exposure $\tau$.
We find a good fit to the
abundances from Ni to La for October with $\tau =  0.2 \pm 0.1$\,mb$^{-1}$
using Malaney's (1987) predictions for a single exposure.
The Rb abundance indicates a low neutron density of
$N_{\rm n} \approx 10^8$\,cm$^{-3}$ (Malaney 1987), while
no useful limit could be set on the Tc abundance.
An exponential distribution of exposures provides less
good agreement with the observed abundances.
For the R\,CrB star
U Aqr, which also shows light $s$-element enhancements,
Bond et al. (1979) obtained $\tau \simeq 0.6$\,mb$^{-1}$. Such exposures
imply that about 10 neutrons were captured by each Fe seed nucleus.
Given that the observed ratio $^{13}$C/Fe $\simeq 10^3$, the exposure, even
in the presence of neutron poisons such as $^{14}$N, seems an achievable
goal. The fact that Ni, Cu and Zn are well fit by the predictions indicates
that the envelope consists mostly of material exposed to neutrons. This fact also
likely explains the anomalous high ratios of K/Fe and Sc/Fe.

The final He-shell flash may occur in a luminous post-AGB star or the
white dwarf that evolves from the post-AGB star.
In the latter case, hydrogen may be mixed  with deep layers of He
and C and consumed. In contrast, the H-burning layer in the post-AGB star
prevents deep mixing. About 10\% of all AGB stars may experience their final
He-shell flash as a white dwarf and, if H consumption is severe, may
convert the born-again AGB star to an R\,CrB star
(Renzini 1990; Iben et al. 1996).
Iben \& MacDonald (1995) have presented a model in which mixing and
nucleosynthesis were followed: their chosen model of a 0.6\,$M_{\odot}$ star
ended with an outer layer having the abundance ratios (by number of atoms)
H/He $\simeq 10^{-0.8}$,
C/He $\simeq 10^{-1.2}$,
N/C $\simeq 10^{-0.5}$, and O/C $\simeq 10^{-1.3}$.
This resembles the
composition of Sakurai's object,
apart from the predicted H deficiency
not being as severe as observed.
The model O/C ratio is lower than observed but might
be raised by adjustment of the uncertain rate for the reaction
$^{12}$C($\alpha,\gamma$)$^{16}$O.
In summary, final flash models offer a
tantalising
prospect of accounting for
the observed composition of Sakurai's object.

Life as a born-again AGB star is brief: the model by Iben \& MacDonald
(1995) brightens by a factor of 10 and cools from T$_{\rm eff}$ of 40,000K
to 6300K in just 17 yr. Evolution over the narrow temperature range
covered by Sakurai's object between May and October is, of course, much
faster.
The evolutionary timescale seems similar to that of V605\,Aql
(Lundmark 1921).
The  timescale for compositional changes for Sakurai's object is
likely even shorter; processed material rising from below will
mix very quickly with the atmosphere.
The other
final flash candidate FG\,Sge has also
showed rapid abundance alterations, e.g. some $s$-process elements
increased by 0.8 dex in 7 years (Langer et al. 1974).

It is now important  to extend the few available calculations
of the final flash to a wider range of
initial conditions, and to include Li-production and $s$-processing. The hints
that the surface composition is evolving rapidly must be pursued by
continued spectroscopic observations,
which may shed further light on its evolutionary
status and relation to the R\,CrB stars;
we may have witnessed the birth of an R\,CrB star.
Furthermore, a determination of the nebular composition would reveal the original
composition of Sakurai's object prior to the final flash.
Monitoring the visual variability of the star
searching for R\,CrB-like declines is naturally of
importance.

\begin{acknowledgements}
We are grateful to Craig Wheeler for bringing the star to our attention, to
Guillermo Gonzalez for a spectrum of Sakurai's object,
and to Jim MacDonald and Icko Iben for helpful comments.
Nikolai Piskunov and Sveneric Johansson are thanked for help with the
hydrogen profiles and atomic data, and
V.P. Arkhipova and V.P. Goranskij for photometric information.
Financial support from the Swedish
Natural Research Council, NSF (grant
AST93-15124) and the Robert A. Welch Foundation of Houston, Texas
is acknowledged.
The analysis has made use of the VALD database.
The referee Simon Jeffery is thanked for helpful comments.
\end{acknowledgements}


Active Galactic Nuclei in polarized light

Abstract: Due to the compactness active galactic nuclei (AGNs) are still unresolved with optical observations. However, structure and physical conditions of the matter in their central parts are especially important to study the processes of the matter accretion to supermassive black holes and eventually these investigations are essential to understand the galaxy evolution. Polarization contains information about the interaction of electromagnetic waves with the environment and provides information about the physical processes in the central regions of the AGNs that could not be found with the help of other optical observations. In this paper, the importance of applying polarimetry methods to the study of geometry, kinematics, and physical processes in active galactic nuclei (AGN) in polarized light is discussed. An overview of the mechanisms of polarization formation, their connection with different structures and scales are provided. Also, we overview the polarimetric investigations based on different assumptions that are done using the observations conducted in Special Astrophysical Observatory of the Russian Academy of Sciences.

\section{Introduction}
Active galactic nuclei (AGN) are bright, compact regions emitting up to 90\% of the entire galaxy energy. Now the high rate of energy release in AGN is associated with the accretion processes onto the supermassive black hole (SMBH). The activity of nuclei is manifested by a number of extraordinary observed properties. Even in the first observations of AGN, the variability in all spectral ranges was observed both on long-term scales of few years (the first works in optics were Sharov\& Efremov 1963, Matthews \& Sandage 1963, Smith \& Hoffleit 1963a,b; in the radio band Dent 1965) and few hours (Wagner \& Witzel 1995). Atypical for normal galaxies emission lines in AGN spectra are also variable (see e.g. Ili{\'c} et al. 2015, Shapovalova et al. 2016, Shapovalova et al. 2017, Shapovalova et al. 2019). Moreover, the amount of observational data on changing-look quasars where broad lines have appeared or disappeared in the spectrum has recently increased (e.g., Khachikian \& Weedman 1971, Cohen et al. 1986, Goodrich 1989 and more recent e.g. MacLeod et al. 2019 and even in polarized light e.g. Hutsemékers et al. 2019). In general, active galaxies are brighter than normal ones, which is especially evident in significant energy excesses in radio, IR, UV bands and at higher frequencies. The spectral energy distributions (SED) of AGN are now considered to be composed of many components (see e.g. Collins et al. 2014), and the major part of radiation is contained in the optical-gamma bands produced in the central region of AGN with an average size of only 1 pc.

However, an important property, which will be in a central place in this article, will be the polarization of active nuclei. As soon as non-thermal synchrotron spectra were detected in AGN, it was assumed by analogy with the radiation of the Crab Nebula the radiation of AGN should be polarized (Dibai \& Shakhovsky 1966). Soon linear polarization was indeed discovered (e.g. Efimov et al. 1979, Hough et al. 1987). As it was predicted by Antonucci (2002) and is absolutely clear today polarization provides additional parameters of the radiation of the studied source and is essential in AGN investigation. The central regions of AGN are now spatially unresolved in optical band, and it is impossible to study them using direct images, so the models of the AGN inner structure have only indirect evidences. Yet, using interferometric methods in radio and IR bands as GRAVITY or EHT it was recently obtained that the principle assumptions of AGN theory are correct such as the accreting SMBH or rotating matter in broad line region (BLR) (Castelvecchi 2019 and Gravity Collaboration et al. 2018, respectively). At the same time, polarization preserves information about both the spatial structure of the nucleus and the physical properties within.

This paper is a short review of polarimetric studies of AGN in the optical range, carried out by observations on the 6-m BTA telescope at the Special Astrophysical Observatory of the Russian Academy of Sciences (SAO RAS), and aims to emphasize the importance of AGN studies in polarized light.

The paper is organized as following. In Sec. 2 we describe the mechanism producing polarization in AGNs in general. Then the individual cases (continuum polarization, broad line polarization and its rapid variability) are analyzed in Sec. 3, 4 and 5, respectively. Finally, in Sec. 6
we outline our conclusions.

\section{Polarization in AGNs}

\begin{figure*}
    \centering
    \includegraphics[scale=0.6]{figure.pdf}
    \caption{"Unified Model in polarized light". It is clear from the figure that the AGN orientation relative to the observer has a strong influence on the dominant polarization mechanism and its orientation. All major AGN components and their sizes are denoted. }
    \label{PUM}
\end{figure*}

As it is known, radiation becomes polarized if it passes through an anisotropic medium. Anisotropy can be caused by various mechanisms, which are usually divided into reflection and transmission processes. The degree of polarization of the reflected radiation is high (up to 100\%), while it usually causes a smaller polarization when passing through the medium. Concerning active nuclei, it is more convenient to divide the polarization mechanisms into internal (inside the central parsec) and external.

The main internal mechanisms of polarization are:

\begin{itemize}
    \item[(i)] the polarization of radiation due to the influence of the magnetic field near the SMBH event horizon within the frames of general relativity (the Lense-Thirring effect);
    \item[(ii)] the polarization of accretion disk radiation due to radiation transfer (Beloborodov 1998) and Faraday rotation in magnetized medium (e.g. Gnedin et al. 2006, Silant'ev et al. 2009);
    \item[(iii)] the optical jet synchrotron radiation;
    \item[(iv)] the scattering by electrons in a hot corona.
\end{itemize}

The internal mechanisms of polarization are summed up, and it is impossible to separate them unambiguously though they have the different direction and they vary on different time scales.

Among the external mechanisms are only two of them:

\begin{itemize}
    \item[(i)] the equatorial scattering at the dusty torus;
    \item[(ii)] the polar scattering at the ionization cone.
\end{itemize}

These mechanisms are clearly observed in Seyfert galaxies (Sy), which are classified into two types due to the features in the spectrum: in Sy 1 there are both broad and narrow emission components, whereas in Sy 2 there are only narrow ones. Antonucci (1982) found that two types of Sy galaxies are also characterized by different orientation polarization: in Sy 1 polarization tends to be parallel to the radio axis, and in Sy 2 to be perpendicular. Within the Unified model (UM, Antonucci 1993; Urry \& Padovani 1995) these observational properties are explained by the orientation of the active nuclei relative to the observer, so the dusty torus can obscure the central regions or not, and therefore this determines the dominant scattering mechanism. Therefore, the different polarization mechanisms dominate in different AGN types. The "unified model in polarized light" is presented in Fig. \ref{PUM}. It describes the types of polarization mechanism depending on the AGN orientation relative to the observer. The AGN orientation relative to the observer has a strong influence on the dominant polarization mechanism and its orientation. At different angles an observer could see different AGN components. At the angles close to $45^\circ$ - i.e. in case of Type 1 AGNs - the central parsec are observed (accretion disk, BLR, region of the jet formation) together with the dusty torus. So, in their spectra both broad and narrow lines are seen. Note here that within the Type 1 AGNs one can also find significant difference as it is concluded in the paper by Popović (2018). At lower angles closer to the equatorial plane - i.e. Type 2 AGNs - the dusty torus obscure the inner regions from the observer and only narrow lines are detectable as NLR is located much farther. The detection of broad emission lines in polarized light in the galaxy NGC 1068, previously considered a typical representative of galaxies of Sy 2 type (Miller et al. 1991), was a good proof of the described assumption of obscuration and the UM in general.  Then, observing the AGN at very high angles one watch approximately into the jet cone, where the radiation is so energetic that its contribution "swamp"\ the details.

Moreover, the polarimetric studies of NGC 1068 were conducted by {\it HST} with the Faint Object Camera in the ultraviolet and with the WF/PC-I in the visual (Capetti et al. 1997 and references therein). Observations have shown the polarization changing across the field. A few years after Kishimoto (1999) reanalyzed the {\it HST} polarimetric data and constructed a model under the assumption of light coming from a point source and scattering at clouds (in this case, NLR), and the two-dimensional polarization distribution was de-projected into a three-dimensional matter distribution in the AGN.

According to the above considerations it can be seen that being associated with the physical state of matter and its spatial distribution polarization becomes an essential tool to study the spatially unresolved areas of AGN. To illustrate this, the results of AGN studies in polarized light in the continuum and broad emission lines, as well as studies of the rapid polarization variability, will be presented below.

Note here that the Stokes parameters of polarization are used in paper. This formalism was introduced by Stokes (1852). Within the frames of this formalism $I$ parameter is the total object intensity, $Q$ and $U$ describes the linear polarization and $V$ is for the circular one. to find a mathematical description see (Walker 1954).

\section{Polarization in continuum}

It is known that in AGN polarization depends on the wavelength, and it s due to different mechanisms of polarization formation. Webb et al. (1993) did the first attempt to separate the different mechanisms of continuum polarization, one of which could be an accretion disk. However, even with the anisotropic distribution of the hot gas in the accretion disk and hence the anisotropic Thomson scattering, the observed polarization should not depend on the wavelength. The situation changes if there is a magnetic field in the accretion disk. Gnedin \& Silant'ev (1984, 1997) found that for weak ($B < 10^3$ G)\footnote{The expected theoretical range for magnetic field strength in SMBH accretion disks is $10^{-2}-10^4$ G (see e.g. Silant'ev et al. 2009 and references therein).} magnetic fields in the accretion disk, optical anisotropy occurs due to the Faraday rotation of the polarization plane at the photon free path, which leads to a strong dependence of the degree of linear polarization on the wavelength $\lambda$. Spectropolarimetric observations of a sample of AGN at the 6-m BTA telescope (Afanasiev et al. 2011) shown that the dependence $P(\lambda) \propto \lambda^n$ exists, where $P$ is the polarization degree and $\lambda$ is wavelegth, and the polarization is formed at small spatial scales – as it will be shown later, less than the BLR size, and can be associated with the processes in the accretion disk. Silant'ev et al. (2009) shown that if the temperature changes along the radius according to the $T(R) \propto R^{-p}$ law (as for the standard Shakura-Sunyaev disk (Shakura \& Sunyaev 1973)) and the magnetic field strength also changes radially $B(R) \propto B_H (R_H/R)^s$, where $R_H$ and $B_H$ are the event horizon radius and the magnetic field strength at it, respectively, then the degree of polarization will change with the wavelength as $P \propto \lambda^{(s/p-2)}$. Comparing this ratio with the observed one, the magnetic field strength ranging from 0.37 to 300 G was obtained for 14 AGN. As it was discussed in the paper (Afanasiev et al. 2011) the magnetic field strength decreases with
increasing increasing SMBH mass, and it is in a agreement with the magnetic coupling model.

Besides, it is possible to determine the radiative efficiency of the accretion disk, depending on the spin of the SMBH. Du et al. (2014) shown that the radiative efficiency could be calculated knowing mass, luminosty, Eddington ratio and the inlination angle of AGN. The determination of inclination angle is the hardest problem here as the systems are not optically resolved. Afanasiev et al. (2018) shown that obtaining the polarization of the source and comparing it with the predicted by Sobolev–Chandrasekhar theory the inclination could be found. Such calculations were made for 47 AGN, and 90\% of the sample show spins >0.9, which corresponds to Kerr black holes.

Also, the polarization in the continuum varies with the amplitude of about 1-2\%. The spectropolarimetric monitoring was conducted at BTA telescope for Sy 1.5 Mrk 6 (Afanasiev et al. 2014) and Sy 1 3C 390.3 (Afanasiev et al. 2015). According to the obtained observational data it was found for Mrk 6 that the delay in broad H$_{\alpha}$ line relative to the continuum at 5100\AA\ is 22 days, which corresponds to the size of the BLR region - 22 lt days (it is consistent with the reverberation mapping results by Bentz \& Katz (2015)). As well, the polarized continuum at 5100\AA\ shows only 2 days delay\footnote{Yet the observations sampling was greater than 2 days (there were 13 epochs at $\sim$1000 days, the method of interpolated cross-correlation function was applied. The linear interpolation was used. Comparison of the linear interpolation with interpolation by exponential continuum model shown that the results were corrected.}, that corresponds to the size by an order smaller than BLR. A similar result was obtained for 3C 390.3: the size of the BLR in H$_{\alpha}$ was 120 lt days (for H$_{\beta}$ was 2 times less, 60 lt days), and the size of the region where polarized continuum forms – 10 lt days. Thus, the region of the polarized continuum, approximately 10 times smaller than BLR, and it is consistent with the idea of the polarized radiation coming from the accretion disk.

\section{Polarization in broad lines}

\begin{figure}
    \centering
    \includegraphics[scale=0.8]{mrk6.pdf}
    \caption{The spectrum and its polarization parameters observed on 2012 Feb 14. From top to bottom: the observed spectrum of Mrk 6, parameters Q and U and corresponding degree $P$ and angle $\varphi$ of polarization. The thin lines represent observed and the thick lines represent corrected spectra for the ISM polarization. The figure is taken from (Afanasiev et al. 2014).}
    \label{mrk6}
\end{figure}

Discussing in this section, Sy 1 galaxies are characterized by the broad emission lines corresponding to the allowed and semi-allowed transitions (as the Balmer series, MgII, CIV, CIII], etc.) in their spectra. It is considered that the lines are formed in the BLR region located at an average distance of 0.01-0.1 pc from the central source. The broadening of the lines is due to the Keplerian rotation of the BLR clouds as a thick disk and shows the velocities up to $10^4$ km/s. Due to the small spatial sizes BLR could not be optically resolved even in close AGN, so only indirect methods are used to study it. The most popular is the reverberation mapping (Blandford \& McKee 1982), which allows us to determine the distance between the BLR and the central source (more precisely – the accretion disk, where exciting radiation is formed) by the delay between the continuum radiation and re-emitted one in the emission line. The BLR stratification was revealed (Peterson 1994; Baldwin et al. 2003) and the characteristic sizes in more than a hundred galaxies were determined. The leading role now plays SDSS-RM project (Shen et al. 2015) aimed to measure time delays for 849 galaxies at the redshifts up to 4.5. Here we will note that the most important achievement of RM programs is recently the relation of AGN luminosity versus the BLR size (Bentz et al. 2009) that indirectly proves the model of the matter accretion onto SMBH. To find more details on RM see (Bentz 2015) However, the reverberation mapping does not allow to conclude about the character of the clouds motion, and also depends on the angle of the system inclination, not determined directly.

Spectropolarimetry is a powerful tool for studying the kinematic properties of BLR. In Sy 1 AGNs, the dominant polarization mechanism in the lines is external equatorial scattering at the dusty torus (the internal mechanisms make a small contribution, and the emission line can be assumed unpolarized before the scattering). Smith et al. (2005) found that due to the Keplerian rotation of the clouds in BLR within equatorial scattering the polarization angle $\varphi$ has an S-shaped profile depending on the wavelength $\lambda$ along the emission line. The example of polarized observations of H$\alpha$ line in Mrk 6 is demonstrated in Fig. \ref{mrk6}. A more detailed simulation using the {\it STOKES} package (Marin 2018) showed that the equatorial scattering does not affect the line profile in total flux, but in polarized light the two-hump structure and the dramatic change of the polarization angle along the profile become visible. Afanasiev et al. (2014) shown that if the motions are Keplerian ones, then the velocity can be expressed as:
$$
\log(v_{i}/c) = a - b\cdot \log(\tan(\varphi_i)),
$$
$$
a=0.5 \log \Big(\frac{GM_{BH} \cos^2(\theta)}{c^2R_{sc}}\Big),
$$
where $v_i$ and $\varphi_i$ is the velocity and the polarization angle for the definite $\lambda_i$,$c$ is the speed of light, $G$ is the gravitational constant, $M_{BH}$ is a central black hole mass, $\theta$  the angle between the disc and polarization plane (see Fig. 1 in (Afanasiev \& Popovi\'c 2015)) and $R_{sc}$ is the distance between the central source and the scattering region. Thus, observing clouds rotating at Keplerian orbits, the obtained dependence $\log(v_{i}/c)$ vs. $\log(\tan(\varphi_i))$ must be linear, which was discovered firstly by modeling, then by observations in the broad H$_{\alpha}$ line in Mrk 6 (Afanasiev et al. 2014) and in 3C 390.3 (Afanasiev et al. 2015), and for for 35 Sy 1 galaxies (Afanasiev et al. 2019).

It is important to note that the parameter $a$ obtained by linear approximation of the dependence $\log(v_{i}/c)$ vs. $\log(\tan(\varphi_i))$ also depends on two important parameters of AGNs – the radius of the scattering region, or dusty torus, $R_{sc}$ and the mass of the central SMBH $M_{BH}$. This fact is the basis of the method of the SMBH mass determination in Sy 1 galaxies with spectropolarimetric observations, described by Afanasiev \& Popovi{\'c} (2015).

Let’s compare the proposed method of the mass determination with reverberation mapping. Assuming the gas is virialized and the size of the BLR is identically equal to the observed signal delay $R_{BLR}=c\tau$, the relation $M_{SMBH} = f {v^2 R_{BLR}}/{G}$ could be used. This estimation is affected by a number of assumptions. The observed gas velocity depends on the inclination angle of the system: $v = v_{obs} / \sin(i)$, and the angle $i$ is not determined from direct observations, since BLR is spatially unresolved. The factor $f$ is also determined empirically. The coefficient $f$ is considered to describe the inclination angle of the system as well as the structure of BLR (e.g. Yu et al. 2019), however, there is no clarity on this issue yet. Thus, on the one hand, reverberation mapping has a significant advantage due to the fact that it uses the velocities at scales much closer than can be obtained from star velocity dispersions. However, a large number of unknown parameters are used for the estimation, which is why the accuracy of the mass measurement lies within the order of magnitude. In addition, to measure the key parameter – the signal delay between the continuum and the line – it requires a long series of observational data, which is implemented for a small number of galaxies and only at small redshifts.

At the same time, mass estimation by spectropolarimetry requires only 1 observation epoch, from which the dependence of the logarithm of the velocity on the logarithm of the polarization angle is obtained. As it can be seen from the expression for $a$, the only value expected additional measurements is the dusty torus size $R_{sc}$. Now there are two ways to determine it: IR reverberation mapping (Minezaki et al. 2004; Suganuma et al. 2006; Koshida et al. 2009, 2014; Mandal et al. 2018) and near-IR interferometry (e.g. Kishimoto et al. 2009, 2011). Comparison of $R_{sc}$ and $R_{BLR}$ (Afanasiev et al. 2019) showed that there is a linear dependence between these quantities, and $R_{sc}$ is on average 2 times more than $R_{BLR}$. However, the measurements of the dust sublimation radius is still a hard question as the theoretical predictions are not in agreement with IR observation. Also it should be mentioned here that except the difficulty with $R_{sc}$ estimation spectropolarimetric method of mass estimation could be applied only for the systems where the rotation motion dominates the radial velocities and out/inflows are weak. For more details of the method application see Savi\'c et al. 2018.

It is important that the spectropolarimetric method of mass estimation does not depend on the orientation of the system. Then having two methods for determining the mass of BH, dependent and independent on the inclination, it is possible to determine the angle $i$. A sample of galaxies shows a linear relation between the inclination angle of the galaxy and the inclination angle of the BLR (Afanasiev et al. 2019). Accordingly, if we know the true angle of the BLR, the distribution of brightness on the disk could be estimated. From observations the distribution is close to the classical Sakura-Sunyaev accretion disk (Afanasiev et al. 2019).

\section{Short-term variability of polarization}

BL Lac type objects (or blazars)\footnote{Though the terms "BL Lacs" and "blazars" are not equal to each other, within this paper we would assume it interchangeably.} are a special type of AGNs with the jet directed at a small angle to the observer’s line of sight. Because of this orientation, the blazars look like star-like sources, their brightness can exceed the host galaxy luminosity by up to 4 mag. It is also known that their luminosity and polarization has a synchrotron origin and are highly variable even within the night. Impey et al. (2000) mentioned for the first time that the polarization of the blazar S5 0716+714 on scales of several hours walks on the $QU$-plane in arcs, that is equivalent to the rotation of the polarization vector. Shablovinskaya \& Afanasiev (2019) shown that on a longer time series the polarization vector in the blazar S5 0716+714 forms arcs and loops on the $QU$-plane, while the change of direction occurs within appoximately 1.5 hours. A similar, though less clear, picture is seen from the analysis of the 2014 microflash data described in Bhatta et al. (2015a,b) and for blazar BL Lac (Covino et al. 2015). Since observing BL Lac type objects one looks inside the jet and polarization appears due to the plasma rotation in the magnetic field, the rotations observed on the $QU$-plane correspond to the physical plasma motion in the jet.
According to the model of the jet structure (Marscher et al. 2008), optical radiation is formed in a helical magnetic field at a distance of $<10^{-2}$ pc from the central source. Based on this fact, the observed plasma motion should be directed along the arcs, and the time of switching the direction of the polarization vector is equivalent to the size of the region where the radiation is formed. Since the region of the optical jet is not spatially unresolved, it is possible to use the polarization vector rotation on the $QU$-plane as a plasma motion test. Thus, similar arcs described by the polarization vector were found in the radio band polarization for the blazar CTA 102 (Li et al. 2018). However, as shown in (Shablovinskaya \& Afanasiev 2019), model of stable helical magnetic field cannot describe the observed rotation of the polarization, therefore the magnetic field precession was suggested, that for S5 0716+714 is approximately 15 days.

\section{Conclusion}

Polarimetry is a powerful tool to investigate the central regions of active nuclei that are unresolved by the direct observations. Being associated with the physical state of matter at different distances from the central SMBH polarization stores information about the inner structure of the nucleus. Observations of polarization carried out in recent years in SAO RAS allowed to obtain the following results.

\begin{itemize}
    \item     The polarization in {\it continuum} is produced in magnetized AD (0.001-0.01 pc) and depends on the magnetic field strength in accreation disk $B(R)$, SMBH mass $M_{SMBH}$ and spin. Spectropolarimetric investigations revealed that $B(R)=0.3-300$ G and spin $\sim$1 that corresponds to Kerr black holes.
    \item The polarization in {\it broad lines} resolves the gas kinematics in BLR (0.01-0.1 pc), and it leads to the more accurate SMBH mass estimation in Sy 1 galaxies, independent from the inclination angle. Also, comparison between these estimations and ones obtained by reverberation mapping the inclination angles were calculated for the number of galaxies and the luminosity distribution of the BLR disks with known inclination were shown.
    \item     {\it Short-term variability} of the polarization vector in BL Lac type objects marks the plasma kinematics inside the jet. Comparing the observational results with the models of the helical magnetic field in jet it was offered that the field is precessing and the region where the optical polarization is formed is of the 1.5 lt hour size.
\end{itemize}

\footnotesize
\section*{Acknowledgments}
We sincerely thank the organizers of 12 Serbian Conference on Spectral Line Shapes in Astrophysics for the special session "Broad lines in AGNs: The physics of emission gas in the vicinity of super-massive black hole" organized in memory of Dr Alla Ivanovna Shapovalova who was our dear colleague and friend.


A teaching guide of nuclear physics: the concept of bond

Abstract: We propose discussions and hands-on activities for GCSE and A-level students, covering a fundamental aspect of nuclear physics: the concept of bond and the energy released (absorbed) when a bond is created (broken). This is the first of the series of papers named "A teaching guide of nuclear physics", whose main goal is to provide teaching tools and ideas to GCSE and A-level teachers, within a consistent and complete curriculum.

\section{Introduction}
\label{sect:intro}

This work is part of a series of papers where we present activities and discussions that teachers can use to
introduce nuclear physics to GCSE and A-level students, from the basic ideas to the most advanced concepts.
Each aims at guiding teachers and students step by step through
the concepts of nuclear physics.
Using analogies with known situations, hands-on activities
and logics, the students will
be able to understand fundamental ideas and mechanisms, as
the principle of minimum energy, endothermic and exothermic reactions,
the creation or breaking of bonds.
Stimulating questions that could spark the students' interest
and foster useful discussions are also provided.

As the students progress through the series of papers, these basic concepts
will allow them to understand a wide variety of phenomena,
as radioactivity, nuclear reactions, nuclear fusion and fission
\cite{0031-9120-52-2-024001, 0031-9120-52-5-054001}.

The present paper focuses on the concept of bond and the energy that is released or absorbed when
a bond is created or broken, respectively.

\section{Learning objectives}
\label{sect:objectives}

After participating in the activities and discussions described in the next sections,
the students will be able to:
\begin{enumerate}[label=(\alph*)]
\item Describe the concept of bond and explain why
energy is released or absorbed when a bond is formed or broken, respectively.
\item Apply the concept of bond and its associated energy to real-life situations
and explain or predict the outcome of an experiment connecting the macroscopic
observables to the microscopic behaviour of atoms and molecules.
\end{enumerate}

Along the paper, the symbol \textbf{(Q)} will be used to denote questions for the students.
These questions aims at stimulating further discussions and stimulating the curiosity of students.
The discussions between them and with the teacher are
essential to achieve the overarching goal of these activities: having the students discuss
concepts of nuclear physics and understand them through logics and experiment.
The teacher acts as a facilitator, promoting the good ideas, giving a direction
to the discussions and encouraging the students to dare making
hypothesis.

We want to stress a very important point here: during the discussions,
the teacher should promote the usage of a
proper scientific language, with correct terms and verbs throughout the activities.
A correct scientific language is at the basis of the communication among scientists.

As these activity could be conducted with both GCSE and A-level students,
additional material containing more advanced discussions and more
detailed studies of some concepts are also presented
along the paper. These additional material is denoted as \textbf{(EXTRA)} and can be skipped
when working with GCSE students.

The paper is organized as follows: in section \ref{sect:bonds}, the definition of bond is discussed.
Section \ref{visual_kit} introduces the visual kit that will be used in
sections \ref{visualize_bonds} to visualize the concept of bond.
Two hands-on experiments will be conducted in sections \ref{hands-on-exp-1} and  \ref{hands-on-exp-2},
where the students will be able to see first hand the effects of creating and breaking bonds at the molecular level.
In sections \ref{evaporation} and \ref{sport_bags}, the concepts that have been learned in the previous sections will be applied to other real-life situations.

\section{Bonds}
\label{sect:bonds}

Understanding the concept of (chemical) bond may result complicated for students
\cite{barker2000students, boo1998students}.
For such a reason, we have developed a simple hands-on activity that may
take place in the classroom and help fixing basic concepts \cite{brown1992using}.

\textbf{(Q)} What is a bond?

\textbf{(Q)} What does it mean to bond something?

\textbf{(Q)} Can you give some examples of bonds in everyday's life?

The word \emph{bond} has of course many meanings. The ones we are interested
in are well described by the following definitions taken from the
Cambridge dictionary:

\begin{enumerate}
\item to bond: to stick materials together, especially with glue, or to be stuck together like this.
(example: \emph{This new adhesive can bond metal to glass})
\item bond: a place where single parts of something are joined together, especially with glue,
or the type of join made (examples: \emph{a strong/weak/permanent bond; When the glue has
set, the bond formed is watertight})
\end{enumerate}

\subsection{The visual kit to visualize bonds}\label{visual_kit}

The following kit (see Fig.\ref{visual_kit_content}) will be used in section \ref{sect:bonds}
to help the students visualize bonds and to understand the concept of the energy associated
with bonds:

\begin{enumerate}[label=-]
\item Two cubes labeled with the letters A and B, well visible. Each cube has Velcro on the sides.

\item A sphere: a spherical wax candle would do.

\item A metal or a wooden spurtle.

\end{enumerate}

\begin{figure}[!h]
\begin{center}
\includegraphics[width=0.70\textwidth,angle=0]{figure/Visual_kit_1b.jpg}
\end{center}
\caption{(Colors online). The visual kit used to understand the concept of the energy associated to a bond.}
\label{visual_kit_content}
\end{figure}

\subsection{Visualizing the energy of a bond}\label{visualize_bonds}
\textbf{(Q)} What happens when you break a bond?

\textbf{(Q)} What happens when a new bond is formed?

\textbf{(Q)} Think of an example of bond. Now imagine breaking this bond, separating the
two objects that are hold together. What does it take to separate them?

Help the discussion toward the concept that breaking a bond
requires energy. Once again, the visual kit can be used to visualize this concept:

\begin{enumerate}[label=-]
\item Place the two cubes on the table, joined to each other with Velcro (see Fig. \ref{visual_kit__separating}).
\item Take the spurtle
and gently start rubbing the spurtle on the Velcro, slowly
separating the two cubes. Explain that the friction between the spurtle, the Velcro and the cubes
is raising the temperature of the system AB.
You could also use an intermediate activity to convince the students of this point:
ask them to rub their hand quickly against the surface of a table. The temperature
of their hand will soon raise, because of friction.
\item The temperature of the system AB
raises, even if just by a tiny amount, due to the friction with the spurtle.
\item While separating A from B, i.e. breaking the bond, you are transferring
energy to the system AB.

\begin{figure}[!h]
\begin{center}
\includegraphics[width=0.70\textwidth,angle=0]{figure/Visual_kit_separating2b.jpg}
\end{center}
\caption{(Colors online). Breaking the bond between A and B using a spurtle. The process gives energy to the AB system. This energy is used to break the bond between A and B.}
\label{visual_kit__separating}
\end{figure}

\item Once the two cubes are separated, set the spurtle aside and summarize what happened
using a third object: a sphere representing the energy
that you transferred to the system AB in order to break the bond.
You may label the sphere with the letter E.
Reestablish the bond between A and B, leaving only the cubes A and B on the table.
Then bring in the energy E and put it next to the bond.
The result is that A and B are separated (see Fig. \ref{visual_kit_equation}).
Hence, when the two objects A and B are separated, they have more energy than before breaking the
bond.
\end{enumerate}

\begin{figure}[!h]
\begin{center}
\includegraphics[width=0.70\textwidth,angle=0]{figure/Visual_kit_equaton.jpg}
\end{center}
\caption{(Colors online). The process of breaking a bond, visualized with the kit. The energy (red sphere) is added
to the AB bound system and, as a result, the bond is broken and A and B are separated. Energy has been absorbed
to break a bond.}
\label{visual_kit_equation}
\end{figure}

\textbf{(Q)} What happens then in the inverse process, when two separated objects are brought
together and a bond is formed?

Let the students play with the visual kit to find an answer.

The conclusion that they will soon draw is that:

\begin{enumerate}
\item Energy is absorbed by the system when a bond is broken (process represented in Fig. \ref{visual_kit_equation})
\item Energy is released by a system when a bond is created (process represented in Fig. \ref{visual_kit_equation_inverse})
\end{enumerate}

\begin{figure}[!h]
\begin{center}
\includegraphics[width=0.70\textwidth,angle=0]{figure/Visual_kit_equaton_inverse.jpg}
\end{center}
\caption{(Colors online). The process of creating a bond, visualized with the kit. The A and B components
of the system get bound and energy (red sphere) is released.}
\label{visual_kit_equation_inverse}
\end{figure}

The reaction that has been represented could be written as follows:

\begin{center}
  AB + energy $\rightleftharpoons$ A + B
\end{center}

where the double arrow indicates that the process can happen in both senses.

Hence, breaking and creating bonds is a way that a system can use
to gain or loose energy. This is a fundamental concept that lies
behind more complex phenomena such as the binding energy of an
atomic nucleus, the stability of a nucleus, radioactivity, nuclear fusion and fission and many others.
A more detailed discussion on this point can be found in Section \ref{conclusions}.

\subsection{Hands-on activity: the energy associated with a bond (part I)}\label{hands-on-exp-1}
The goal of this activity is to experimentally observe the energy associated
with bonds. The students will be able to apply the concepts that they learned
in the previous discussions to a concrete case and explain the observed phenomena.

The activity is made of two experiments. The first experiment proceeds as follows:

\begin{enumerate}
\item Split the students into small groups
\item Each group of students is given an experimental kit
\item The kit for this first experiment contains (see Fig. \ref{hands-on_material_1}): a large bowl, a beaker containing 500 mL of room-temperature water (as an alternative any container with the same volume and a scale to measure 500 mL of water), about 100 g of kitchen salt (NaCl), a thermometer (minimum range required: -50 $^{\circ}$C to +50 $^{\circ}$C), a rod or a spoon to stir the water
\footnote{In case that a graduated beaker or container is not available, a scale could be used to weigh the 500 mL of water.}.
\item Once the students have poured 500 mL of water into the beaker, they have to
measure the temperature of the water (allowing enough time to the
temperature of the water to stabilize after a thermal equilibrium with the container is reached).
\item After that, they have to dissolve salt in the water, stirring the water in order to facilitate the dissolution of the salt
\footnote{WARNING: stirring too quickly could raise the temperature of the system and
alter the results of the experiment by introducing an external source of energy. It is fundamental that one stirs
the water gently.}.
\item It is important that each group measures the temperature at all times and takes notes of what is happening: is the temperature changing? Does all the salt dissolve in the water?
There is no need to explain what they have to focus on: let them experiment with the quantity of salt
that they like and give them enough time to discuss and draw some conclusion.
\end{enumerate}

\begin{figure}[!h]
\begin{center}
\includegraphics[width=0.70\textwidth,angle=0]{figure/hands_on_salt_material_1_crop.jpg}
\end{center}
\caption{(Colors online). The material for the first part of the hands-on activity. From the top left, in clockwise order: a large bowl, a beaker with 500 mL of room-temperature water, about 100 g of kitchen salt (NaCl), a thermometer, a rod or a spoon to stir the water, a scale.}
\label{hands-on_material_1}
\end{figure}

\begin{figure}[!h]
\begin{center}
\includegraphics[width=0.48\textwidth,angle=0]{figure/hands_on_water_18_degrees.jpg}
\includegraphics[width=0.48\textwidth,angle=0]{figure/hands_on_water_17_degrees.jpg}
\end{center}
\caption{(Colors online). The execution of the first part of the hands-on activity. (Photo on the left) Measuring the initial temperature of the water and then starting pouring salt in it. (Photo on the right) Stirring the water with a spoon in order
to dissolve the salt and, at the same time, measuring the temperature at all times. The temperature decreases by
one or two Celsius degrees.}
\label{hands-on_execution_1}
\end{figure}

The temperature of the water typically drops by a couple of Celsius degrees (see Fig. \ref{hands-on_execution_1}).

This first experiment ends with a discussion among all students and the teacher. Guide them through
an analysis of the outcome of the experiment with questions:

\textbf{(Q)} What happened when you dissolved the salt in the water?

\textbf{(Q)} Why did the temperature drop? How much did it drop?

\textbf{(Q)} Did it depend on the quantity of salt?

\textbf{(Q)} Why didn't the temperature raise instead of dropping?

\textbf{(Q)} And why did the temperature change in the first place?

Remind them of the conclusion on the energy associated to a bond and that a system can gain or lose
energy by breaking or creating bonds.

\textbf{(Q)} What bonds could have been created or broken while dissolving the salt in the water?

\textbf{(Q)} If the temperature dropped, does it mean that bonds have been created or broken?

Guide them through the complexity of the process: when a grain of salt is thrown to water,
water molecules surround the grain of salt; bonds are created between water and salt molecules,
while bonds are broken between water molecules to allow the water molecules to join the salt ions;
more and more bonds are created between water molecules and the salt ions
until the water molecules pull salt atoms away from the grain, one by one; bonds are created between
water molecules and salt atoms (Na$^{+}$ cations and Cl$^{-}$ anions) while the bonds between Na and Cl atoms and the bonds between water molecules are broken.
When the salt is dissolved, there is a huge number of Na$^{+}$ and Cl$^{-}$ ions in the solution, each of them
surrounded by many water molecules. Each bond between a water molecule and a salt ion is weak, but the
number of bonds is large.

\textbf{(EXTRA)} If the students are interested and if their level of chemistry allows it,
the discussion can temporarily divert to the kind of bonds that keep the salt ions together (ionic bonds)
and the kind of bonds that water molecules can form with each other and with the salt ions
(hydrogen bonds) \cite{klein2017organic}.

\textbf{(Q)} Some bonds are created, then, while some other bonds are broken. Which one is larger: the energy released
by the bonds that are created or the energy that is absorbed by the bonds that are broken? How can you tell?

\textbf{(Q)} How could we make this effect larger? How could we make the temperature drop much more dramatically?

The stronger the broken bonds are, the more the energy they absorb and the more the temperature drops.
The bonds that were broken in the first experiment were: the bonds between water molecules and the bonds
between the salt ions. We can not do much about the bonds between the salt ions, but we can surely make
the bonds between water molecules much stronger: we could use ice instead of water.

\textbf{(EXTRA)} Energy is absorbed to break the bonds and this makes the temperature drop. The link between the two is
given by the atomic theory and the kinetic theory of gases \cite{boltzmann1896gases, kauzmann1914kinetic}.
The energy to break the bonds is taken from the kinetic energy
of the surrounding molecules, which then slow down. Slower molecules correspond, on a macroscopic point of view, to a lower temperature.

\subsection{Hands-on activity: the energy associated with a bond (part II)}\label{hands-on-exp-2}

\begin{figure}[!h]
\begin{center}
\includegraphics[width=0.70\textwidth,angle=0]{figure/hands_on_salt_material_2_crop.jpg}
\end{center}
\caption{(Colors online). The material for the second part of the hands-on activity. From the top left, in clockwise order: a large bowl, a bag of ground ice (at least 400 g), a beaker with 500 mL of room-temperature water, about 100 g of kitchen salt (NaCl), a thermometer, a rod or a spoon to stir the water, a scale.}
\label{hands-on_material_2}
\end{figure}

The second experiment proceeds as follows:
\begin{enumerate}
\item Keep the same groups of students as in the first experiment
\item Each group of students is given an experiment kit
\item The kit for this second experiment contains (see Fig. \ref{hands-on_material_2}): a large bowl, a bag of ground ice (at least 400 g), a beaker containing 500 mL of room-temperature water (or, as an alternative, any container with the same volume and a scale to weigh 500 mL of water), about 100 g of kitchen salt (NaCl), a thermometer (minimum range required: -50 $^{\circ}$C to +50 $^{\circ}$C), a rod or a spoon to stir the water
\item The students have to pour 500 mL of water into the beaker and add about 10 cubes of ice. It is important
that all groups add approximately the same amount of ice, so that possible differences cannot be associated
to a different setup of the experiment.
\item The students have to measure the temperature of the water+ice (wait for the
temperature to stabilize after a thermal equilibrium with the container has been reached)
\item After a thermal equilibrium is reached, they can dissolve the salt in the water, stirring the water to
facilitate the dissolution of the salt.
\item It is important that each group measures the temperature at all times and takes notes of what is happening: is the temperature changing? Does all the salt dissolve in the water? Let them experiment with the quantity of salt
that they like and give them enough time to discuss and draw some conclusion.
\item Once the temperature seems not to drop anymore, let the students add more ice. Let them experiment freely. You could even transform this part of the experiment into a competition: which group will reach the lowest temperature?
\end{enumerate}

\begin{figure}[!h]
\begin{center}
\includegraphics[width=0.48\textwidth,angle=0]{figure/hands_on_ice_0_degrees.jpg}
\includegraphics[width=0.48\textwidth,angle=0]{figure/hands_on_ice_minus_18_degrees.jpg}
\end{center}
\caption{(Colors online). The execution of the second part of the hands-on activity. (Photo on the left) Measuring the initial temperature of the water+ice system (around 0 degrees centigrades) and then starting pouring salt in it. (Photo on the right) Stirring the solution with a spoon in order
to dissolve the salt and, at the same time, measuring the temperature at all times. The temperature decreases very quickly and reaches about -18 $^{\circ}$C.}
\label{hands-on_execution_2}
\end{figure}

Before adding salt, the temperature of the water+ice system is about zero degrees centigrade, which
is the melting temperature of the ice. After adding the salt, the temperature
gradually drops to -18 $^{\circ}$C (see Fig. \ref{hands-on_execution_2}). The more salt is dissolved, the lower the temperature gets.

\textbf{(Q)} How can you explain what happened in terms of bond creation or destruction?

Compared to the first experiment, there is an additional bond that is broken: the bond between water molecules in the ice. As this bond is much stronger than the bond between water molecules in water, the energy absorbed when this bond is broken
is much higher. Hence the temperature drops much more with respect with the first experiment.
This also explains why the temperature of the water+ice+salt system does not drop so strongly if the ice
is floating on the surface of the water, without getting in contact with the salt on the bottom of the container.
The most dramatic drop of temperature is attained if the ice cubes touch the salt on the bottom, as this causes
a large number of bonds to be broken between water molecules in the ice. Hence, a good suggestion to the students, in
case they were not able to reach the -18 $^{\circ}$C, would be to get rid of some of the water, so that the ice cubes
could get in contact with the salt.

\subsection{Hands-on activity: evaporation of alcohol}\label{evaporation}
In this activity, the students will be able to apply what they learned in the previous sections
to the case of the evaporation of a liquid.

\begin{enumerate}
\item Keep the same groups of students as in the previous experiments
\item Give a bottle of room-temperature denatured alcohol to each group and ask the student to pour
some drops of alcohol on the back of their hands and to wait for a couple of minutes.
\item Ask them to be patient and to register any change that they may feel.
\end{enumerate}

\textbf{(Q)} What happened after you poured some alcohol on your skin?

\textbf{(Q)} Did it feel hotter or colder?

\textbf{(Q)} Write down a possible explanation for the fact that the interested area
felt colder after some time.

\textbf{(Q)} What does it mean, at the molecular level, that the alcohol evaporated?

The molecules of denatured alcohol are short carbon chains made of only two carbon atoms
(CH3-CH2-OH) \cite{klein2017organic_b}. This facilitates their evaporation.
Evaporation means that the bonds between the alcohol molecules at the surface of the liquid are broken
and these molecules can then abandon the liquid and enter the gas phase of the surrounding air.
To break the bonds between the alcohol molecules in the liquid phase, energy is needed, as we
learned in the previous sections. This energy is taken from the kinetic energy of the molecules
that the skin is made of and this in turn corresponds, at a macroscopic level, to a lower temperature.

\subsection{Hands-on activity: sport bags}\label{sport_bags}
In this activity, the students will be able to apply what they learned in the previous sections
to the case of the well-known bags that are used in sports to cool down or heat up very quickly an athlete's body.

\begin{enumerate}
\item Keep the same groups of students as in the previous experiments
\item Give one bag to each group. Some of the bags, when shacked, will cool down quickly. Some others
will heat up.
\item Tell the students that they cannot open the bags, but they can shake them or fold them. Let them experiment
freely until they discover that the temperature of the bag raises or drops after shaking it.
\end{enumerate}

\textbf{(Q)} What happened with your bag? Did the same happen for all groups?

\textbf{(Q)} Write down a possible explanation for what happened with your bag.

\textbf{(Q)} Write down a possible explanation for the fact that other groups' bag experienced
the opposite change of temperature.

After the discussion, have the different groups compare their explanation and agree on a common conclusion.
The difference between the two kind of bags is that some salts allow a large number of water molecules
to surround each salt ion. This results in a larger number of bonds that are formed and hence
a larger energy that is released. If the energy released is larger than the energy that is required
to break the water-water and ion-ion bonds, then the system's temperature raises.
With some other salts, the number of water molecules surrounding each ion is not large enough
and the energy released is less than the energy required to break the water-water and ion-ion bonds.
Hence the temperature of the system drops.

\section{Conclusions}\label{conclusions}

It is fundamental to help the students to draw conclusions at the end of the activities and to summarize
what has been achieved. This helps the students to fix the ideas and to put them into context.

The most effective way to summarize the activity is doing it interactively. Once again, the teacher acts
as a facilitator, promoting the good ideas and giving a direction to the discussions.

\textbf{(Q)} What fundamental concept have we touched upon in these activities?

\textbf{(Q)} Why is the concept of bond so important?

Guide the students toward the idea that creating and breaking bonds allows a system to exchange
energy with the exterior, and that energy exchanges are essential for a system to evolve. Without entering
into the details, you could mention the principle of minimum energy that states that for a closed system
the internal energy will tend to decrease and to approach a minimum at equilibrium. A closed system is a system that can exchange energy with the exterior, but not matter. The principle of minimum energy, which will be covered in
the following papers of this series, implies that a system naturally evolves toward a minimum energy configuration
through energy exchanges with the exterior.

We saw examples where the effect of a huge number of bonds
leads to a macroscopic increase or decrease of the temperature of the system. The students can be guided through
this line of reasoning with questions as:

\textbf{(Q)} Why did the temperature drop when we added salt to the water with ice?

\textbf{(Q)} What energy transformations are taking place in the process?

The total energy absorbed by the broken bonds is larger than that released by the creation of new bonds. This energy
is taken from the kinetic energy of the surrounding molecules and this is perceived as a drop in temperature.
In conclusion, have the students bring a clear message home: the bonds  and their associated energy are two fundamental
concepts in physics, and as the activities and the discussions in the following papers of this series will show,
these concepts are at the basis of many important phenomena in nuclear physics.
\section*{Acknowledgments}
A.P. acknowledges a support from the STFC Grant No. ST/P003885/1 and ST/P006213/1.