Multi-frequency radio observations of BAL quasar 1045+352

Abstract: Multi-frequency 1.7, 5 and 8.4-GHz VLBA observations of a radio-loud broad absorption line (BAL) quasar 1045+352 are presented. It is a young compact steep spectrum (CSS) object and its asymmetric, two-sided morphology on a scale of several hundred parsecs, extending in two different directions, may suggest intermittent activity. The young age and unusual morphology of 1045+352 are arguments in favour of an evolution scenario for BAL quasars, in which the BAL features appear at a very early stage of their evolution.

\section{Introduction}

Approximately 10\% of optically selected radio-quiet
quasars display broad absorption lines (BALs) in the blue wings of the
high ionization resonant lines
(C\,IV 1549\AA~-- high-ionization BAL (HiBAL) quasars), caused by the
outflow of gas with velocities up to $\sim0.2\,c$ \cite{hewett03}.
10\% of them also show absorption troughs in low ionization lines
(Mg\,II 2800\AA~-- low-ionization BAL (LoBAL) quasars).
Evidence has accumulated from optically selected BAL quasars, in favour of
an orientation hypothesis to explain their nature.
According to this hypothesis, BAL regions exist in both BAL and non-BAL
quasars, and that BAL quasars are normal quasars seen along a particular
line of sight, e.g. a line of sight skimming the edge of the accretion disk
or torus \cite{weymann91}. This view has been challenged by
the discovery of the existence of a large population of radio-loud BAL quasars
\cite{broth98,becker00, menou01}.
Most of these quasars tend to be compact in the radio
domain with either a flat or steep radio spectrum. It follows
that those belonging to the latter class seem to be related to
Gigahertz-Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS)
sources. This variety of
spectral indices suggests a wide range of orientations, contrary to the
interpretation favoured for optically selected quasars.

The radio morphology of BAL quasars is important because it can
serve as an inclination indicator of BALs, and therefore yields a direct
test of the orientation model. However, information about the radio
structure of BAL quasars is still very limited. Prior to 2006, only
three BAL quasars, FIRST J101614.3+520916 \cite{gregg00}, PKS 1004+13
\cite{wills99}
and LBQS 1138$-$0126 \cite{broth02} were known to have a
double-lobed FR\,II radio morphology on kiloparsec scales, although this
interpretation is doubtful for PKS 1004+13 \cite{gopal00}.
Recently, excluding PKS 1004+13, the population of FR\,II-BAL quasars has increased to ten
objects \cite{gregg06, zhou06}, although
some of them still require confirmation.
Their symmetric structures indicate an ``edge-on'' orientation, which in
turn supports an alternative hypothesis described as ``unification by time'',
with BAL quasars characterised as young or recently refuelled quasars
\cite{becker00, gregg00}. There had been only one attempt (and only at
1.6\,GHz with the EVN) to image radio
structures of the smallest (and probably the youngest) BAL quasars
\cite{jiang03}.
In this paper we present high frequency VLBA images of another very
compact BAL quasar -- 1045+352.

\section{Observations}

1045+352 belongs to the primary sample of 60 candidates for CSS sources
selected from the VLA FIRST catalogue \cite{wbhg97}. Initial observations of all the
candidates were made with MERLIN at 5\,GHz \cite{kun02} and 1.7, 5 and
8.4-GHz VLBA follow-up of 1045+352
was carried out on 13 November 2004 in a snapshot mode with
phase-referencing. The target source
scan was interleaved with a scan on a phase reference source
and the total cycle time (target and phase-reference) was
$\sim$9~minutes including telescope drive times, with $\sim$7\,minutes
actually on the target source per cycle.

The whole data reduction process was carried out using standard AIPS
procedures. For the target source and at each frequency, the corresponding
phase-reference source was mapped and the phase errors so determined were
applied to the target sources, which were then mapped using a few cycles of
phase self-calibration and imaging. IMAGR was used to produce the final
``naturally weighted'' images shown in Fig.~\ref{1045+352_maps}.

\begin{table}[t]
\caption[]{Basic parameters of 1045+352}
\begin{tabular}{@{}l r@{}}
\hline
\hline
~~~~~Parameter & Value~~~~~\\
\hline
~~~~~Source name (B1950)     & 1045+352~~~~~\\
~~~~~Source right ascension (J2000) extracted from FIRST     & $10^{\rm
h}48^{\rm m}34\fs247$~~~~~\\
~~~~~Source declination (J2000) extracted from FIRST     & $+34^{\rm o}
57\farcm24\farcs99$~~~~~\\
~~~~~Total flux density $S_{1.4\,GHz}$~(mJy) extracted from FIRST     &
1051~~~~~\\
~~~~~log$P_{1.4\mathrm{GHz}}$~(W~${\rm Hz^{-1}}$)     & 27.65~~~~~\\
~~~~~Total flux density $S_{4.85\,GHz}$~(mJy) extracted from GB6     &
439~~~~~\\
~~~~~Spectral index $\alpha_{1.4\mathrm{GHz}}^{4.85\mathrm{GHz}}$     &
$-0.70$~~~~~\\
~~~~~Largest Angular Size measured in the 5-GHz MERLIN image     &
$0\farcs5$~~~~~\\
~~~~~Largest Linear Size ($H_0$=100${\rm\,km\,s^{-1}\,Mpc^{-1}}$,
$q_0$=0.5)~($h^{-1}$~kpc) & 2.1~~~~~\\
\hline
\end{tabular}

\label{table1}
\end{table}

\section{Results and discussion}

According to SDSS/DR5, 1045+352 is a galaxy at
RA=\,$10^{\rm h}48^{\rm m}34\fs242$, Dec=\,$+34^{\rm
o}57\farcm24\farcs95$,
which is marked with a cross in the MERLIN map (Fig.~\ref{1045+352_maps}).
However, the spectral observations \cite{willott02} have shown
1045+352 to be a quasar with a redshift
of $z=1.604$. It has also been classified as a HiBAL QSO based
upon the observed very broad
C\,IV absorption, and is a very luminous submillimetre object
with
detections at both 850\,$\mu$m and 450\,$\mu$m \cite{willott02}.

Our MERLIN and VLBA maps of 1045+352 (Fig.~\ref{1045+352_maps}) show this source to be
extended in both NE/SW and NW/SE directions. The central compact feature
with a steep spectrum
visible in all the maps is probably a radio core. The VLBA
image at 1.7\,GHz shows radio structure extended in
a NE/SW direction, with the SW jet being weaker than
that in the NE. The NE/SW emission is visible
in the 5-GHz MERLIN image as a few compact features.
The 5-GHz VLBA image shows also a core and a one-sided jet pointing to the East.
The radio structure in the 8.4-GHz VLBA image is similar to that at
5\,GHz: an extended radio core and a jet pointing in an easterly direction.

The observed radio morphology of 1045+352 could indicate a restart
of activity accompanied by a reorientation of the jet axis with the NE/SW
radio emission being the first phase of activity, now fading away, and
the extension in the NW/SE direction being a signature of the current
phase of activity.

\begin{figure*}[t]
\centering
\includegraphics[width=7cm, height=7cm]{1045+352.merlin.ps}
\includegraphics[width=7cm, height=7cm]{1045+352.18cm.ps}
\includegraphics[width=7cm, height=7cm]{1045+352.6cm.ps}
\includegraphics[width=7cm, height=7cm]{1045+352.4cm.ps}
\caption{The MERLIN 5-GHz (upper left) and VLBA 1.7, 5 and 8.4-GHz maps of
1045+352. Contours increase
by a factor 2 and the first contour level corresponds to $\approx 3\sigma$.
Cross indicates the position of an optical object found using the SDSS/DR5.}
\label{1045+352_maps}
\end{figure*}

A reorientation of the jet axis may result from a merger with another
black hole. It has been shown \cite{merritt02} that a rapid change of jet
orientation can be caused by even a minor merger because of a
spin-flip of the central active black hole arising from the coalescence of
inclined binary black holes.
Also, the Bardeen-Petterson effect \cite{bardeen} can cause a realignment of a rotating
supermassive black hole (SMBH)
and a misaligned accretion disk \cite{liu04}, the timescale of such a realignment
being $t<10^{5}$ years.
The interaction/realignment of a binary and its accretion disk leads to the
development of X-shaped sources \cite{liu04}.
1045+352 is not a typical X-shaped source, however, the
realignment of a rotating SMBH followed by a repositioning of the
accretion disk and jets is a plausible interpretation for misaligned radio
structures, even if they are not conspicuously X-shaped \cite{cohen05}.

The radio luminosity of 1045+352 at 1.4\,GHz (Table~\ref{table1}) is high, making this
source one of the most radio-luminous BAL quasars, with a value
similar to that of the first known radio-loud BAL\,QSO with an FR\,II
structure, FIRST\,J101614.3+520916 \cite{gregg00}.
The radio-loudness parameter,
$R^{\ast}$, defined as the {\it K}-corrected ratio of the 5-GHz radio flux
to 2500\AA~optical flux \cite{stocke92} has been calculated for 1045+352
and amounts to $\log(R^{\ast})=4.9$.
For this, a global radio spectral index, $\alpha_{radio}=-0.8$ and an
optical spectral index, $\alpha_{opt}=-1.0$, have been assumed and the SDSS
$g'$ magnitude has been converted to the
Johnson-Morgan-Cousins {\it B} magnitude using the transformation formula
\cite{smith02}.
Corrections have also been made for intrinsic extinction (local to the
quasar) calculated assuming a Milky-Way extinction curve \cite{willott02}.
Even after correction, $\log(R^{\ast})>1$ and amounts to
$\log(R^{\ast})=4.1$, which means that 1045+352 is still radio-loud object.
The two-sided asymmetric structure of 1045+352 can indicate the source
grows up in an asymmetric environment which is very typical among small scale
CSS sources \cite{saikia01}.
However, it was also found \cite{saikia01,
jeyakumar05} that the radio properties of CSS sources are consistent
with the unified scheme in which the axes of the quasars are observed close
to the line of sight. It is then likely that in the case of the BAL
quasar 1045+352 the orientation effect is substantial.

\section{Conclusions}

1045+352 is a very radio-luminous BAL quasar,
whose complex structure is suggestive of restarted activity and may have
resulted even from a minor merger.
The radio morphology of 1045+352 shows also that the radio jets and BALs can
coexist, at least during some stage of the quasar lifetime. According to
\cite{lipari06}, in radio-loud systems, the jets then remove the clouds
responsible for the generation of BALs, so the extended radio structures
showing BAL features are very rare.

1045+352 is also a CSS object.
CSSs constitute an intermediate class in an evolutionary sequence of radio
sources between smaller GPS objects and very extended Large Symmetric
Objects (LSO). They are young sources with a typical age of $\sim10^{5}$ years
\cite{mur99}, which is an argument in favour of an evolutionary scenario
for BAL quasars, in which the BAL features appear at a very early stage in
the quasar evolution.

\section*{Acknowledgement}
This work was supported by Polish Ministry of Education and Science
under grant 1 P03D 008 30.


Cycles on the Solar-Type Stars and Cooler Dwarfs

Abstract: Features of the development of activity cycles in the solar-type stars and fast-rotating cool dwarfs have been considered for 65 stars observed in some decades. Cycles with duration of 7-18 years compared to the solar cycle were found for about 50% of the studied stars. In cooler dwarfs with rotation periods of less than 5 days, cyclic changes in brightness occur on longer scales, up to 80 years. Activity of the highest level is produced on K dwarfs; their main cycles are long and have the highest amplitudes. Both old and young solar-type stars show a similar tendency in increasing the cycle length with a slower rotation. No evidence for a relation between the rotation period and duration of cycles was found for cool dwarfs with $P_{rot} < 5$ days.

\section{INTRODUCTION}

The magnetic activity typical for the solar-type stars (spectral classes F8-K2, $T_\textsl{eff} = 6200–4500\;$K) and cooler K--M dwarfs with effective temperatures of up to $2200\;$K determines the diversity of variability observed in different wavelength ranges. The time scale of variability is in the range from seconds (fast flares) to tens of years (activity cycles similar to the solar magnetic cycle); there are also secular cycles of hundreds and thousands of years. On the Sun, manifestations of activity in different layers of the atmosphere are regulated by an 11-year cycle. The development of magnetic activity cycles is determined by the rotation of a star and its internal structure.

According to Nizamov et al.~(2017), cycles in G-type dwarfs are established from the epoch when their rotation period reaches 1.1 days, for K and M dwarfs these periods are longer -- 3.3 and 7.2 days, respectively. The stars of the corresponding spectral classes, rotating faster, demonstrate the saturated activity regime.

A certain number of stars with cycles longer than the solar one are known, but most of the studied objects were investigated from the modern databases on a limited time interval; therefore, cycles of no longer than 10--15 years were found for these stars (Savanov, 2012; Suárez Mascareño et al., 2016; Boro Saikia et al., 2018). The results concerning stellar cycles require further specification and confirmation from the long-term systematic observations.

Oláh et al.~(2016) derived the series of data on time scales up to 36 years for 29 G-K stars from the Mount Wilson HK-Project and determined the duration of their dominant cycles (with a maximum amplitude) and short cycles.

The long-term variations in brightness of several dozens of cool K--M dwarfs were studied from photographic archives and broadband photometry (Bondar’, 1995, 2013; Alekseev and Bondar’, 1998; Alekseev and Kozhevnikova, 2017); for some of them the activity of photospheric starspots similar to the observed on the Sun has been found or suspected.

The task of our work is to form a sample of stars, whose cycles have been found from the data over several decades; to substantially extent the available photometric series for some G--M dwarfs by the results of recent observations; to refine or determine cycles of their photospheric activity.

Besides, we will compare characteristics of dominant cycles of solar-type stars and cooler dwarfs obtained over long-term series, consider how they are related to rotation periods, reveal features of this connection for different type stars, that is important for the dynamo theory.

\section{G--M DWARFS WITH CYCLIC ACTIVITY\\ ON LONG-TERM INTERVALS}

The first results on activity cycles of F--M stars on or near the main sequence were obtained by Wilson (1978) within his monitoring program of 91 stars in the Ca II H and K emission lines. Based on observations in 1966--1977, for some stars there have been suspected cycles (trends) which require certain evidence from the longer data series. The program was continued in the framework of the HK-Project for 111 stars. In the paper by Baliunas et al.~(1995), which became the basis for studying magnetic activity of the solar-type stars, the cycles were searched for and classified from the data series which include the 12-year observations by O.\ Wilson and new data, and in whole cover the 25-year interval.

Oláh et al.~(2016) presented new data about chromospherical activity cycles for 29 HK-Project stars from the 36-year observation series. Most of these stars exhibit cycles with a duration of 7--18 years, i.e. which correspond to the sunspot cycle. Some stars, besides the main cycle, demonstrate one or several shorter cycles (multiple cycles), as well as long-term trends of 25 years or longer.

\begin{figure*}[!th]
\centering
\includegraphics[width=\linewidth]{fig1.jpg}
\caption{
Search for the rotation period of V647 Her. Top panel: V -- light curve on the interval of JD2453128--2453160 (2004) from the SuperWASP data. Bottom panel: periodogram of the presented data, the peak corresponds to a period of 1.087 days (left), the phase fit of V-magnitudes with $P_\textsl{rot} = 1.087$\;days (right).
}
\label{Figure1}
\end{figure*}

For our purpose, we take into account the results of a search for cycles on cool dwarfs from the series longer than 20 years, i.e. an investigating time spans must be comparable with those for the solar-type stars in Oláh et al.~(2016).

We have used the results derived for 40 dK--dM stars from photographic data series (Bondar’, 1995, 2013), from photometric researches by Alekseev and Kozhevnikova (2017) and Messina et al.~(2002). For 9 stars (Bondar’ 2019) the available data series were added by data from photometric catalogs. Data from the Hipparcos catalog (\url{http://cds.u-strasbg.fr/cgi-bin/Dic-Simbad?HIP}) cover\\ the span of 1990–1993, data from the ASAS catalog (Pojmanski, 1997) were obtained in 2002–2009, data from SuperWASP (\url{https://wasp.cerit-sc.cz/form}) cover the interval between 2004–2008, data obtained in 2008\-–\-2018 were taken from the KWS database (Kamogata / Kiso / Kyoto Wide-field Survey, \url{http://kws.cetus-net.org/~maehara/VSdata.py?object}).

\begin{table*}[!th]
\centering
\vspace*{0.5cm}
\caption{
G- and M-type dwarfs with active cycles determined from new photometric data series.
}
\vspace*{0.5cm}

\begin{tabular}{l|c|l|c|r|c|c|l|l}
Star & B-V & Sp type & $P_\textsl{rot}$, d & $P_\textsl{cyc}$, yr & $A_\textsl{cyc}$, mag & Time span & Variable type & Sources\\
\hline
BE Cet & 0.66 & G2.5 & 7.76 & 6.8 $\pm$ 0.1 & $<$ 0.04 & 1986–2018 & BY Dra & 2, 4, 5\\
DX Leo & 0.77 & G9 & 5.42 & 3.3 $\pm$ 0.1 & 0.03 & 1989–2018 & BY Dra & 1, 2, 4, 5\\
V1005 Ori & 1.37 & M0 & 4.40 & 20.4 $\pm$ 0.2 & 0.12 & 1976–2018 & BY Dra & 1, 2, 4, 7\\
DT Vir & 1.48 & M0 & 2.88 & 31.6 $\pm$ 0.4 & 0.15 & 1911–2018 & fl & 1, 2, 4, 6\\
V647 Her & 1.64 & M3.5 & 1.09 & 76.1 $\pm$ 0.6 & 0.30 & 1939–2018 & fl & 1, 2, 3, 4, 6\\
V577 Mon & 1.69 & M4.5 & 8.00 & 40 $\pm$ 0.4 & 0.50 & 1909–2018 & fl & 1, 2, 4, 6\\
DX Cnc & 2.06 & M6.5 & 0.46 & 2.7 $\pm$ 0.1 & 0.10 & 1984–2018 & fl & 2, 4, 6\\
V374 Peg & 1.26 & M3.5 & 0.44 & var & 0.10 & 1989–2017 & fl & 1, 2, 3, 4\\
HU Del & 1.65 & M4.5 & 0.35 & var & 0.32 & 1940–2017 & fl & 2, 4, 6\\
\end{tabular}
\vspace*{0.5cm}

{
Remarks. Spectral type, B-V and type of variability are taken from the SIMBAD astronomical database; numbers in column “Sources” mark databases Hipparcos (1), ASAS (2), SuperWASP (3), KWS (4); Messina and Guinan, 2002 (5), Bondar’, 2013 (6), Alekseev and Kozhevnikova, 2017 (7).
The type of stellar variability is marked according to the SIMBAD astronomical database: BY Dra-type stars (BY Dra), flare stars (fl).
}
\label{Table1}
\end{table*}

The list of stars and used data sources are given in Table~1. Two stars from the list -- BE Cet and DX Leo -- belong to G-dwarfs. On the interval between 1986--2000 Messina and Guinan (2002) found only short-term cycles for these stars. We obtained data series covering a time span twice as long, until 2018; this allows us to consider the stability of short cycles and to search for possible long-term cycles. The cooler stars in our sample have spectral classes M0V--M6.5V; for some of them the long-term photographic measurements of brightness were used from Bondar’ (1995). Two stars, V374 Peg and HU Del, are of interest as fast-rotating fully convective dwarfs. Light curves of V374 Peg from the photometry in 1998--2013 were studied by Vida et al.~(2016), but the authors found no cycles on this star. The parameters of cycles (the cycle length and amplitude) given in Table 1 were obtained as a result of the analysis of composite light curves formed by us, including recent data of 2017 and 2018. We made a search for the cycles from periodograms derived by the Yurkevich, Hartley and Scargle methods using the AVE program (\url{http://www.gea.cesca.es}) as described in (Bondar’ 2019).

Note that the errors in the $P_\textsl{cyc}$ values show the accuracy of their determination, but not in the duration of cycles, which may vary up to several years at different epochs.

Variations in light curves of V374 Peg and HU Del are not cyclic, and the parameter $A_\textsl{cyc}$ shows the magnitude of the maximum changes in their yearly mean brightness over the studied time interval. No long-term cycles were found on the BE Cet and DX Leo stars. We confirm that only the short cycles detected by Messina and Guinan (2002) have been found from the 30-year photometric series. Considering parameters of the cycle, we took information about the duration of cycles for V833 Tau (78 yr) on the interval of 1899–2009 from Bondar’ (2015) and for YZ CMi (27.5 yr) from data on the time span of 1926–2009 (Bondar’, 2019).

The rotation period of the star V647 Her was determined according to the SuperWASP data-base obtained with high temporal resolution. A search for the period was performed using the AVE program by the Hartley and Scargle methods on the JD2453128 – 2453160 (2004) interval containing 1292 points for periods of 0.1–5 days and 0.1–10 days. The periodogram in Fig.\ 1 (the right panel) shows the presence of a peak corresponding to a period of 1.09 days determined by the Scargle method (Scargle, 1982).

We also confirmed the rotation periods of the stars DT Vir and V374 Peg. According to our analysis of data from the ASAS catalog since 2003 to 2008, the rotation period of the star DT Vir $P_\textsl{rot} = 2.88\;$d is consistent with
$P_\textsl{rot} = 2.85\;$d derived by Donati et al.~(2008); for V374 Peg, according to data in 2004 from the SuperWASP catalog,
we found $P_\textsl{rot} = 0.448\;$d to agree with $P_\textsl{rot} = 0.446\;$d acquired by Vida et al.~(2016).

\section{ROTATION PERIODS AND\\ PARAMETERS OF ACTIVITY CYCLES}

Here we consider 65 G–M dwarfs, for which a search for cycles has been performed on long-term intervals. A sample includes 36 stars studied from the photometric series and 29 stars, including the Sun, from the Mount Wilson HK-Project (Oláh et al., 2016). The number distribution of objects by spectral types is as follows: 14 G, 28 K, and 23 M dwarfs.

\begin{figure*}[!th]
\centering
\vspace*{-0.9cm}
\includegraphics[width=\linewidth]{fig2.jpg}
\vspace*{0.2cm}
\caption{
Distribution of rotation periods and duration of cycles for stars from our sample.
}
\label{Figure2}
\end{figure*}

The rotational periods are known for 59 stars of the sample; cycles were detected or suspected for 62 stars. Fig.\ 2 (left panel) shows the distribution of periods in bins, each of them is equal to 5\ d, the spectral classes are marked with different symbols, the dashed line indicates the total number of stars in the bin. Most of G-type dwarfs (10 objects) have periods from 5 to 20 days and only the star HD 81809 (G2V) rotates slower than the Sun, with a period of 39.3 days. Fast-rotating dwarfs with periods of less than 5 days include 25 stars from the sample; mainly, these are cool dwarfs. Three M-stars have rotation periods of less than 0.5 day. 13 K--M dwarfs fall into the range of rotation periods of 30–60 days, the slowest ones are two objects HD 95735 (M2V) and HD 16160 (K3V) with rotation periods of 54 and 57 days, respectively.

The distribution on the length of cycles is shown in Fig.\ 2 (right panel), bins correspond to 7 years. The cycle durations of 50\%\ of stars are from 7 to 21 years that is comparable to the solar spot cycle. The long cycles, up to 80 years, are detected in fast-rotating K--M dwarfs. The length of dominant cycles for 15 stars is short, less than 7 years. The long-term brightness variations of fully convective M-dwarfs HU Del and V374~Peg with $P_\textsl{rot} <0.5$~d showed no certain periodicity. The light curves for some stars and the results of a search for the cycles are given in (Bondar', 2019).

\begin{figure*}[!th]
\centering
\vspace*{-0.9cm}
\includegraphics[width=\linewidth]{fig3.jpg}
\vspace*{0.2cm}
\caption{
The amplitude and duration of cycles of G–M dwarfs with different rotation periods.
}
\label{Figure3}
\end{figure*}

The amplitude of the dominant cycle in the solar-type stars does not exceed 0.2$^m$, and for the fast-rotating cool dwarfs it can increase up to 0.5$^m$ or more; that indicates on the large areas of active regions. As follows from Fig.~3, K-dwarfs with any rotation periods are distinguished from G and M dwarfs by such parameters as cycle duration and amplitude; this leads to the conclusion that their level of activity is the highest among dwarfs on or near the main sequence.

\section{FEATURES OF THE DEVELOPMENT\\ OF CYCLES IN G–M DWARFS}

The development of activity cycles in the solar-type stars mainly depends on their rotational velocity. On the diagram $P_\textsl{cyc}$ and $P_\textsl{rot}$ (Bondar’, 2019) stars are divided into two sequences near a value of $P_\textsl{rot} = 5\;$d. For stars with $P_\textsl{rot} > 5$~d, the cycle length increases if a rotation period becomes slower. This sequence is not homogeneous, since it includes objects with periods from several days to tens of days, i.e. having different ages. For the older stars, rotating slower than the Sun, the above tendency is less pronounced than for the young G and K dwarfs with $5 < P_\textsl{rot} < 20$~days. For these stars cycles of more than 22 years have not been detected yet, but among young K dwarfs there are stars with cycles of about 40 years. The Sun is located between groups of old and young stars.

The youngest stars with $P_\textsl{rot} <5$~d show no any definite connection between the rotation and cycle length. Baliunas et al.~(1996) derived a relation between the observed parameters $P_\textsl{cyc}$ and $P_\textsl{rot}$ and the theoretical dynamo number $D$ expressed as $P_\textsl{cyc}/P_\textsl{rot} \sim D^i$. Later, such a relation was considered in several papers for different samples of stars (Oláh et al., 2009; Savanov, 2012; Vida et al., 2014; Katsova et al., 2015; Suárez Mascareño et al., 2016). Vida et al.~(2014) have found for all the long and short cycles the value of the slope $i \sim 0.77$, which agrees with Baliunas et al.~(1996), but M dwarfs have been excluded from the analysis.

Suárez Mascareño et al.~(2016) and Savanov (2012) have considered only M stars; the stellar cycles were determined from observations on intervals of less than 10 years. The obtained large values of $i$ indicate that for the cooler dwarfs such a relation is weaker or absent.

\begin{figure*}[!th]
\centering
\vspace*{-0.9cm}
\includegraphics[width=\linewidth]{fig4.jpg}
\vspace*{0.2cm}
\caption{
Relation between the number of revolutions in the cycle and the rotation period for G--M dwarfs. The spectral classes of stars are indicated as in Fig.\ 2. The position of the Sun is shown by the corresponding symbol. The solid regression line fits all the stars of our sample, the dashed line shows the regression line from Oláh et al.~(2016).
}
\label{Figure4}
\end{figure*}

Figure 4 shows the essential difference between the slopes of regression lines for the entire sample of stars (solid line, $i = 1.25 \pm 0.068$) and for the dominant cycles of stars from Oláh et al.~(2016) (dashed line, $i = 0.76 \pm 0.15$). The differences in slopes are caused not only by possible uncertainties in the determination of cycle parameters, but are related to a large fraction of rapidly rotating dwarfs in our sample. However, there is still no a sufficient number of stars with $P_\textsl{rot} <5$~d to analyze such a relation for each spectral group.

\section{CONCLUSION}

The stars on or near the main sequence are related by structure (they have convective shells) and the observed manifestations of activity provided by the presence of magnetic fields. Different methods of observations on long-term time intervals reveal slow quasi-periodic changes in brightness of active dwarfs on intervals compared to the duration of the solar cycle. Such a pattern of changes is caused by the evolution of surface spots. For a number of stars, the analogy between the 11-year Schwabe cycle and the dominant cycles of stars is confirmed by the long-term observations. Multiple cycles, the butterfly effects were found in some stars studied in detail. The longer cycles, by several decades, are assumed to be caused by a change in spottedness of the stellar surface. Confident evidence for the evolution of spots on such time spans may be obtained from the long-term observations in different photometric bands. Modern methods for measuring magnetic fields make it possible to detect magnetic spots of different polarities, and from the many-year data one can identify cycles similar to the Hale cycle. The possibility of detecting such cycles on the stars YZ CMi, DT Vir, DX Cnc and others is discussed by Shulyak et al.~(2015). Differences in the duration of dominant cycles are associated with differences in the inner structure of stars, their physical parameters, and rotational periods. Certain results on the relation between characteristics of stars and their activity level were obtained for the solar-type stars (Oláh et al., 2016; Fabbian et al., 2017) from many-year observations of indicators of activity in the chromosphere and photosphere. Later stars are mainly investigated on time intervals of about 10 years.

On the basis of long-term observation series, features in the development of activity cycles were considered for 65 G--M dwarfs with rotation periods from fractions of a day to several tens of days. No cycles of more than 40 years were found in the solar-type stars; for most of them, activity cycles are in the range from 7 to 21 years, i.e. compared to the 11-year solar cycle, also having du-ration from 7 to 18 years at different epochs. Cycles of the fast-rotating dwarfs can reach 80 years. So far, such long cycles have been suspected in two stars -- V833Tau and V647 Her. About 30\%\ of the considered stars show only short dominant cycles of less than 7 years.

For rapidly rotating dwarfs, there is no definite relationship between the cycle length and the rotational period, which is typical for stars with rotation periods longer than 5 days. Their cycle length increases with a deceleration of rotation. But in this group of stars, the longer cycles are observed in K-dwarfs rotating faster than the Sun.

\section*{ACKNOWLEDGEMENTS}

The authors are grateful to the referee for useful remarks and suggestions.

We used data from the following databases: ASAS and SuperWASP catalogs, the International Variable Star Index (VSX) database, the Kamogata – Kiso – Kyoto Wide-Field Survey data. The authors are thankful to all the staff for supporting these databases and providing access to data.

\section*{FUNDING}
This work was carried out with the partial support of the Russian Foundation for Basic Research (grants 18-52-06002 Az$_a$ and 19-02-00191a).

\section*{CONFLICT OF INTEREST}
The authors declare that they have no conflicts of interest.

\phantomsection
\bibliographystyle{unsrt}

Solar and Stellar Astrophysics (astro-ph.SR)


AGN Feedback in groups and clusters of galaxies

Abstract: The lack of very cool gas at the cores of groups and clusters of galaxies, even where the cooling time is significantly shorter than the Hubble time, has been interpreted as evidence of sources that re-heat the intergalactic medium. Most studies of rich clusters adopt AGN feedback to be this source of heating.From ongoing GMRT projects involving clusters and groups, we demonstrate how low-frequency GMRT radio observations, together with Chandra/XMM-Newton X-ray data, present a unique insight into the nature of feedback, and of the energy transfer between the AGN and the IGM.

\section{Introduction}

Simulations of structure formation, using gravitational interactions
between dark matter particles, reproduce, to reasonable accuracy, the
large-scale properties of the observed web of galaxies and clusters
(e.g. Springel et al. 2005). However, when gas physics are added, the
simulated galaxies fail to match some of the basic observed properties
(e.g., Croton et al. 2005).  To produce the observed luminosity function
of galaxies, the rate of gas cooling has to be set too high, leading
to a severe overestimation of the fraction of baryons in stars, and an
underestimate of the fraction in the X-ray emitting gas. Thus,
additional sources of heating are required:
gravitational heating, supernovae, mergers and conduction are among plausible
candidates, but most studies of rich clusters adopt AGN heating to be
predominant (see McNamara \& Nulsen 2007 and references therein).

The nature of this feedback, vital to our understanding of galaxy and
structure evolution, is one of the most important unresolved questions
in extragalactic astronomy. In the cores of groups and clusters of
galaxies, radiative cooling times are shorter than the Hubble time
(see Fig.~1a). Thus, {\it cooling flows} should result, yet the
spectral signature of gas cooling to low temperatures is not seen
(Peterson \& Fabian 2006). This should be linked to the same
source of feedback that is required by the simulations.  Furthermore, the
raised entropy of the gas within galaxy groups, compared to that
expected from scaling cluster properties (``symmetry breaking''), shows
evidence of the heating of baryons at modest redshifts before they are
assembled into clusters, providing more support for feedback on
the scale of galaxies (Ponman, Cannon \& Navarro 1999).

\section{Observational study of AGN feedback}

Most studies of feedback have hitherto focused on AGN jet/cavity
systems in the most massive, X-ray luminous galaxy clusters, where the
disturbed structures are most clearly seen. However, most galaxies and
most of the baryonic matter in the Universe reside in substantially
smaller groups, where feedback in principle has the greatest impact on
galaxy formation and evolution. It is vital to study feedback in a
sample with a wide variety of groups {\it and} clusters, if we are to
understand how feedback has influenced the thermal history of galaxies
and the intergalactic medium (IGM), and thus of most of the baryons in
the Universe.

 \begin{figure}
\centering
\includegraphics[width=5.0cm]{somak-fig1a.eps}
\includegraphics[width=5.9cm]{somak-fig1b.eps}
 \caption{{\itshape (a, left)\/} The cooling time of the
X-ray emitting intergalactic medium at the cores of galaxy clusters
is less than the Hubble time. The dashed lines represent
clusters in which the expected cool core is not found (from
Sanderson et al. 2006). This is true of X-ray emitting
groups of galaxies as well.
 {\itshape (b, right)\/} Cavities, created in the IGM
as a result of the  plasma injected by the central radio source,
are confined by the pressure of the IGM. They buoyantly
rise through the IGM,
transferring a significant part of their enthalpy.
Here, in the Chandra X-ray image of the cluster Abell 2597
(McNamara et al. 2001), one such detached (``ghost'') cavity
is marked.
}
 \end{figure}

\subsection{The synergy of X-ray and radio observations}

The IGM in groups and clusters emits X-rays, while the relativistic
particles from the AGN produce radio emission through synchrotron
radiation.  However, no simple relation has been established between
radio power and the energy required to destroy cool cores.  A unified
study of observations in multiple wavebands is required to examine the
detailed physics of this process of transfer of energy between the AGN
and the IGM.  The high resolution of the Chandra observatory reveals
highly disturbed structures in the cores of many clusters, including
shocks, cavities and sharp density discontinuities.  Comparing with
radio maps, it is clear that AGN jets are associated with many of
these disturbances.

The most typical configuration is for jets from
the central dominant elliptical of a group or cluster to extend
outwards and inflate lobes of radio-emitting plasma, pushing aside the
X-ray emitting gas of the cluster halo to create cavities in the hot
IGM (Fig.~1b). The enthalpy of a cavity filled with relativistic
plasma is $4pV$, where $p$ and $V$ are the pressure and volume
respectively (McNamara \& Nulsen 2007), and a large fraction of this
can be transferred to the IGM as the buoyant cavities travel outwards
(see, e.g., McNamara \& Nulsen 2007, Jetha et
al. 2008).

 \begin{figure}
\centering
\includegraphics[width=6.0cm]{somak-fig2a.eps}
\includegraphics[width=6.0cm]{somak-fig2b.eps}
 \caption{
In the NGC 741 fossil group, the jets are emitted from
a non-central dwarf galaxy, NGC 742.
{\itshape (a, left)\/} Our 610~MHz GMRT contours, superposed
on an optical image, and
{\itshape (b, right)\/} Our 235~MHz GMRT contours, superposed
on the 0.5-2~keV Chandra X-ray image, show far more
extended radio emission than previous VLA maps. Lowest contours are at
0.15 and 0.9 mJy/beam respectively.}
 \end{figure}

\subsection{GMRT studies of galaxy groups}
\label{Obsgrp}

Groups have crucial advantages over richer clusters for studies of
feedback -- the gas in their relatively shallow potential wells is
expected to respond more readily to energy input, and the lower
temperature of the IGM would lead to stronger metal emission lines.
Our in-depth study of 18 elliptical-dominated galaxy groups, which
show structures indicative of AGN interaction in deep Chandra/XMM
X-ray images, at 3 GMRT frequencies (235-610 MHz), demonstrates
(Vrtilek et al. 2009, O'Sullivan et al. 2009) that by combining deep
X-ray and GMRT observations, the extent of the AGN outbursts can be
better ascertained, and the chance of detecting earlier outbursts
enhanced, than in $>$1~GHz observations (Giacintucci et al. 2008).
These observations (Figs.~2 \& 3a) allow one to measure the spectral
index as a function of position along the jets, and thus to carefully
model the extended emission, and, using equipartition values for the
ambient magnetic field, to derive the age of the outburst.  This leads
to an estimate of the duty cycle of AGN activity, and helps to
decouple multiple outbursts, if present.  From all of this, in turn,
an accurate estimate of the energy of the AGN outburst can be made.

\subsection{GMRT studies of feedback in galaxy clusters}
\label{Obsclus}
Previous studies of feedback in clusters (e.g. Birzan et al. 2004,
Dunn et al. 2005, Fabian et al. 2006, McNamara \& Nulsen 2007) have
concentrated on the core regions of clusters where X-ray images reveal
the presence of cavities and related features.  To assess how
different kinds of radio outbursts contribute to cluster heating, and
how common it is to find AGN activity related to their brightest
galaxies, a statistical sample needs to be studied. We have obtained
GMRT 610~MHz images of all 26 GMRT-accessible clusters ($z\!<\! 0.2$)
of the REXCESS sample (B\"ohringer et al. 2007). Data reduction is in
progress: we show a very preliminary map of one of these clusters in
Fig.~3b.  These clusters have extensive multiwavelength observations
(including deep X-ray observations), which we will be able to use in
conjunction with the radio observations to model the behaviour of the
AGN and its effect on the cluster IGM.

 \begin{figure}
\centering
\includegraphics[width=5.5cm]{somak-fig3a.eps}
\includegraphics[width=6.8cm]{somak-fig3b.eps}
\caption{
{\itshape (a, left)\/}
Smoothed Chandra X-ray image of the HCG~62 compact group
of galaxies, in which two
  $\sim$8~kpc cavities are detected.  Our GMRT observations (here
  235~MHz contours on the Chandra 0.3--2 keV image) show far more
  extended emission than previous VLA observations (Vrtilek et al. 2002),
   and evidence of
  progressive spectral ageing along the jets.
{\itshape (b, right)\/} (Preliminary) GMRT 610~MHz contours superposed on
the XMM-Newton image of the 2.5~keV
REXCESS cluster RXJ2023.0-2056 at redshift $z=0.056$,
showing clear detection of jets
associated with the central galaxy.
}
 \end{figure}

\section{Conclusions}
The most likely sources of feedback in galaxy groups and clusters are
active galactic nuclei.
However, the radio power of AGN in central galaxies measured at
high radio frequencies is not sufficient to explain the evidence of
heating in the IGM of most groups and clusters.
Low frequency GMRT observations, along with high-resolution
Chandra X-ray observations, are playing a crucial role in
the detailed study of the history of AGN outbursts
and the mechanism of energy transfer, in an ongoing series of
investigations.

\acknowledgements
We thank the staff of the GMRT (run by NCRA-TIFR)
for making these observations possible.

Cosmology and Nongalactic Astrophysics (astro-ph.CO); Astrophysics of Galaxies (astro-ph.GA)


Enhanced Superconductivity in Close Proximity to the Structural Phase Transition of Sr1-xBaxNi2P2

Abstract: The structural evolution and superconductivity of a 122-type solid solution Sr1-xBaxNi2P2 were studied. We found that an orthorhombic-tetragonal structural phase transition takes place at x = 0.5, and is characterized by the P-P dimers breaking. The superconducting transition temperature exhibited its highest value of 2.85 K at x = 0.4.

\section{Introduction}
The transition-metal compounds AT$_2$X$_2$ exhibit a variety of crystal structures, and the most common is the tetragonal ThCr$_2$Si$_2$-type (space group $I4/mmm$, $D_{4h}^{17}$, No. 139). \cite{Hoffmann}
This structure can further be classified into two types:
the ``collapsed" tetragonal that is characterized by the formation of molecule-like X-X dimers along the crystallographic $c$-axis, and the ``uncollapsed" tetragonal in which the X-X dimers are broken.
Interestingly, some compounds exhibit structural phase transitions from the uncollapsed to collapsed tetragonal phases when hydrostatic pressure or chemical doping is applied.
Typical examples include CaFe$_2$As$_2$ under pressure, \cite{Kreyssig,Goldman} CaFe$_2$(As$_{1-x}$P$_{x}$)$_2$, \cite{Kasahara} and Ca(Fe$_{1-x}$Rh$_{x}$)$_2$As$_2$, \cite{Danura} in which the formation of X-X bonds results in the loss of the Fe magnetic moment and the disappearance of superconductivity.
Other notable examples are Sr$_{1-x}$Ca$_{x}$Co$_2$P$_2$,\cite{Jia} LaCo$_2$(Ge$_{1-x}$P$_{x}$)$_2$,\cite{Jia1} and SrCo$_2$(Ge$_{1-x}$P$_{x}$)$_2$,\cite{Jia2} in which the ferromagnetic quantum critical point is induced by the X-X dimer breaking. \cite{Jia2}

A rare but intriguing structure among AT$_2$X$_2$ appears in SrNi$_2$P$_2$.
This compound crystallizes in an orthorhombic structure with the space group $Immm$ ($D_{2h}^{25}$, No. 71) \cite{Keimes}.
In this structure, the unit cell is made up of 12 phosphorus atoms, of which four phosphorus atoms form dimers while the remaining eight do not. This results in a superstructure along the crystallographic $b$-axis, and hence causes orthorhombic distortion, as shown in Fig. 1(a).
On the other hand, BaNi$_2$P$_2$ crystallizes in an uncollapsed tetragonal structure ($I4/mmm$), \cite{Keimes} in which no P-P dimers are formed, as shown in Fig. 1(b).
Both compounds exhibit superconductivity, and the superconducting transition temperatures $T_{\rm c}$ are 1.4 K in SrNi$_2$P$_2$\cite{Ronning} and 2.5 K in BaNi$_2$P$_2$\cite{Mine,Tomioka}.
This raises questions regarding which composition causes the P-P dimer breaking to occur in the solid solution Sr$_{1-x}$Ba$_{x}$Ni$_2$P$_2$, and how $T_{\rm c}$ varies as a function of the composition $x$.

In this paper, we present the results of our synthesis of the SrNi$_2$P$_2$--BaNi$_2$P$_2$ solid solution along with the crystallographic and superconducting properties of the samples.
Our observations suggest that there was interplay between the structural transitions and superconductivity in Sr$_{1-x}$Ba$_{x}$Ni$_2$P$_2$.

\begin{figure}[t]
\begin{center}
\includegraphics[width=8cm]{fig1_Kudo.eps}
\caption{
(Color online) The crystal structure of (a) SrNi$_2$P$_2$ (space group $Immm$, $D_{2h}^{25}$, No. 71) and (b) BaNi$_2$P$_2$ ($I4/mmm$, $D_{4h}^{17}$, No. 139).
}
\end{center}
\end{figure}

\begin{figure}[t]
\begin{center}
\includegraphics[width=7cm]{fig2_Kudo.eps}
\caption{
(Color online) The $x$ dependencies of the lattice parameters and the interlayer P-P distances in Sr$_{1-x}$Ba$_x$Ni$_2$P$_2$ at room temperature.
The subscripts T and O denote the tetragonal and orthorhombic phases, respectively. See Fig. 1 for P(1) and P(2) of the orthorhombic phase.
}
\end{center}
\end{figure}

Polycrystalline samples of Sr$_{1-x}$Ba$_x$Ni$_2$P$_2$ with nominal Ba content of 0.00 $\le$ $x$ $\le$ 1.00 were synthesized by means of a solid-state reaction. The stoichiometric amounts of the starting materials Ba, Sr, Ni, and P were mixed, ground, placed in an alumina crucible, and then sealed in an evacuated quartz tube. The ampoule was heated at 300 $^\circ$C for 5 h, at 700 $^\circ$C for 3 h, and at 1000 $^\circ$C for 24 h, followed by natural cooling in the furnace. Powder X-ray diffraction (XRD) studies confirmed that the resulting samples were a single phase of Sr$_{1-x}$Ba$_x$Ni$_2$P$_2$. The lattice parameters and interlayer P-P distances were estimated by the Rietveld refinement using the RIETAN-FP program. \cite{Izumi} The magnetization $M$ was measured using the Magnetic Property Measurement System (MPMS) by Quantum Design at temperatures above 1.8 K.

\begin{figure}[t]
\begin{center}
\includegraphics[width=8.5cm]{fig3_Kudo.eps}
\caption{
(Color online) The temperature dependence of the magnetization $M$ divided by the magnetic field $H$ for Sr$_{1-x}$Ba$_x$Ni$_2$P$_2$ with (a) $x$ $<$ 0.50 and (b) $x$ $\ge$ 0.50 at $H=$ 30 Oe in the zero-field and field cooling conditions. No correction for the diamagnetizing field was made.
}
\end{center}
\end{figure}

\begin{figure}[t]
\begin{center}
\includegraphics[width=7cm]{fig4_Kudo.eps}
\caption{
(Color online) The $x$ dependencies of (a) the superconducting transition temperature $T_{c}$ and (b) the superconducting volume fraction determined from the magnetization data measured in the field cooling conditions in Sr$_{1-x}$Ba$_x$Ni$_2$P$_2$ (blue circles). The $T_{\rm c}$ and volume fraction for the end member SrNi$_2$P$_2$ (red triangles) were determined from specific heat data in Ref. \citen{Ronning}.
}
\end{center}
\end{figure}

Based on the XRD data measured at room temperature, the structure at $x <$ 0.5 was confirmed to be the orthorhombic structure ($Immm$), whereas at $x >$ 0.5 the structure was confirmed to be the uncollapsed tetragonal ($I4/mmm$).
The orthorhombic--tetragonal structural transition can be seen in the $x$ dependence of the $c$ parameter, which exhibited a jump at $x$ = 0.5, as shown in Fig.~2(a).
On the other hand, the variation in the $a$ parameter was very small, as shown in Fig.~2(b).
The orthorhombic phase at $x <$ 0.5 can be characterized by the two distinct P-P distances, P(1)-P(1) and P(2)-P(2), as shown in Fig.~2(c).
The P(2)-P(2) distance was comparable to the P-P distances of 2.251--2.710 {\AA} in the collapsed tetragonal phase of CaTM$_2$P$_2$ (TM = Fe, Co, Ni, Cu),  \cite{Hoffmann} and thus indicates a covalent bond.
A large P(1)-P(1) distance suggests the absence of a covalent bond between them.
In contrast, in the tetragonal phase at $x >$ 0.5, the P-P distance was uniform and much longer than that of the covalent P-P bond, as shown in Fig.~2(c), and the covalent P-P bonds were completely broken.

Figure 3 shows the temperature dependence of the magnetization $M$ for Sr$_{1-x}$Ba$_x$Ni$_2$P$_2$ at a magnetic field of 30 Oe in the zero-field and field cooling conditions.
Using this data, we estimated the superconducting transition temperature $T_{\rm c}$ and superconducting volume fraction as a function of the doping $x$, as summarized in Fig.~4.
Bulk superconductivity with a substantial superconducting volume fraction was observed at $x$ $\geq$ 0.4 for Sr$_{1-x}$Ba$_x$Ni$_2$P$_2$.
The superconducting transition temperature exhibited the highest value of $T_{\rm c}$ = 2.85 K at $x$ = 0.4.
This $T_{\rm c}$ value was twice as high as the $T_{\rm c}$ = 1.4 K of the end member SrNi$_2$P$_2$. \cite{Ronning}
$T_{\rm c}$ decreased monotonically with increasing $x$, and
the end member BaNi$_2$P$_2$ exhibited $T_{\rm c}$ = 2.53 K, which is in good agreement with previous reports. \cite{Mine,Tomioka}
In contrast, the superconductivity at 0.15 $\leq$ $x$ $<$ 0.4 is not bulk in nature.
The superconducting volume fraction decreased with decreasing $x$ below $x$ = 0.4, and the $T_{\rm c}$ was almost independent of $x$.
This behavior suggests the coexistence of superconducting and non-superconducting phases, and their phase boundary is located at $x$ $\simeq$ 0.4.
However, the XRD data suggested the formation of a continuous solid solution in this doping range. The orthorhombic--tetragonal structural phase boundary was located at $x$ $\simeq$ 0.5, as shown in Fig. 2.
To account for this discrepancy, another phase transition, presumably to the collapsed tetragonal structure, should be invoked below room temperature.
Incidentally, SrNi$_{2}$(P$_{1-x}$As$_{x}$)$_{2}$ is known to exhibit three structural phases, namely the orthorhombic ($Immm$), collapsed, and uncollapsed tetragonal ($I4/mmm$) structures, depending on the temperature and composition. \cite{Keimes2}

In conclusion, we have demonstrated that the SrNi$_2$P$_2$-BaNi$_2$P$_2$ solid solution exhibits a structural phase transition that is characterized by the breaking of the P(2)-P(2) dimers of the orthorhombic structure.
The superconducting transition temperature was enhanced up to $T_{\rm c}$ = 2.85 K when the composition was in the vicinity of the structural transition.
This $T_{\rm c}$ was higher than those of the end members SrNi$_2$P$_2$ (1.4 K) and BaNi$_2$P$_2$ (2.5 K).
Thus, the tuning of the X-X bonding state in AT$_2$X$_2$ with the ThCr$_2$Si$_2$-type structure is an effective means to develop superconductors with enhanced $T_{\rm c}$.

\acknowledgment
A part of this work was performed at the Advanced Science Research Center, Okayama University.
This work was partially supported by Grants-in-Aid for Scientific Research (No. 26287082, 15H01047, 15H05886, and 16K05451) provided by the Japan Society for the Promotion of Science (JSPS)
and the Program for Advancing Strategic International Networks to Accelerate the Circulation of Talented Researchers from JSPS.


Light propagation and fluorescence quantum yields in liquid scintillators

Abstract: For the simulation of the scintillation and Cherenkov light propagation in large liquid scintillator detectors a detailed knowledge about the absorption and emission spectra of the scintillator molecules is mandatory. Furthermore reemission probabilities and quantum yields of the scintillator components influence the light propagation inside the liquid. Absorption and emission properties are presented for liquid scintillators using 2,5-Diphenyloxazole (PPO) and 4-bis-(2-Methylstyryl)benzene (bis-MSB) as primary and secondary wavelength shifter. New measurements of the quantum yields for various aromatic molecules are shown.

\section{Introduction}
In many particle physics experiments using organic liquid scintillators, precise knowledge on the parameters affecting the propagation of scintillation light inside the detector is required. This knowledge is needed to estimate the detector response and to simulate detector signals under various conditions. In particular, for large liquid scintillator (LS) detectors in the ton-scale basic properties of the LS molecules as the absorption and emission spectra or the fluorescence quantum yields need to be known. Such detectors are for example used in neutrino physics, since the cross sections for neutrino detection reactions are tiny and therefore large target masses are needed. The understanding of light propagation in combination with calibration data is paving the way for the determination of the energy deposited inside the detector by ionizing particles.

In chapter 2 we describe the light absorption and emission properties of some substances typically used in LS neutrino experiments. In databases these spectra are normally given assuming low concentrations and an inert medium around the fluorescent molecules. Both conditions are not fulfilled for standard liquid scintillators which can be affected by solvent and concentration effects. Therefore the spectra as measured in the corresponding solvent and at a given concentration might be more relevant. In addition, the spectrum can change as the light propagates through the medium.

In chapter 3 we present results of our new measurements on the fluorescence quantum yields for several molecules used in various LS neutrino detectors. It is hard to control the systematic uncertainties in such measurements, which might be the reason why the literature values of the absolute numbers are spread over a wide range with significant deviations from each other. The agreement between quantum yield ratios of molecules measured in our setup compared to the corresponding ratios in other publications is found to be more robust.

\section{Scintillator absorption and emission}

Ionizing particles traversing liquid scintillators mainly excite the solvent molecules of the liquid. Deexcitation can occur essentially via the emission of fluorescence light or non-radiative energy transfer to another molecule. Since the aromatic solvent molecules are not transparent for their own flourescence light, in particular in large scale LS detectors, primary and sometimes secondary wavelength shifters are added to the liquid. Most modern experiments use PPO (2,5-Diphenyloxazole) as primary and bis-MSB (4-bis-(2-Methylstyryl)benzene) as secondary wavelength shifter.

For the modelling of light propagation inside the LS it is important to know on which molecule the photon is absorbed. UV light in wavelength regions below 280 nm is mainly absorbed by the solvent molecules. Benzene derivatives typically have absorption bands around 260~nm and reemit light around 300~nm~\cite{Ber}. Ideally the region of solvent emission matches the absorption peak of the primary wavelength shifter which dominates the absorption in the mixture from around $280-350$~nm. The scintillation light seen by the photomultipliers in a large scale LS detector is above 350~nm. In this region the secondary wavelength shifter is typically dominating the absorption. Therefore, it is very important to have a high quantum yield for these molecules, since this parameter determines the reemission probability after absorption of scintillation light in the liquid. The absorption of secondary wavelength shifters like bis-MSB becomes negligible in the region above 430~nm. There the impurity levels in the chemicals are most relevant for the attenuation length. Light absorbed by the impurities is typically not reemitted. The regions of main absorption for the individual LS components are illustrated in Figure~\ref{Fig1}.

\begin{figure}[tb]
\begin{center}
\includegraphics[scale=0.50]{Figure_01.png}
\end{center}
\caption[]{In an organic liquid scintillator system with a primary and secondary wavelength shifter the absorption bands of the individual components should be reasonably separated. The bars in the plot show for each component the wavelength region for which its own absorbance is dominating in the scintillator mixture.}
\label{Fig1}
\end{figure}

The emission spectra of PPO and bis-MSB diluted in cyclohexane at concentrations of 0.5~mg/l are shown in Figure~\ref{Fig2}. In addition, the emission spectrum for a LS mixture as measured in a 2~mm cell in the so-called front face geometry is plotted. For this geometry, emission spectrum of the sample is recorded at the same side where the excitation takes place. The sample contains 3~g/l PPO and 20~mg/l bis-MSB dissolved in pure ortho-phenylxylylethane (o-PXE), a high flashpoint scintillator solvent used e.g.~in the Borexino CTF~\cite{PXE} and in the Double Chooz experiment~\cite{DC1}. This sample was excited at 260~nm, where the light mainly excites the o-PXE.

\begin{figure}[tb]
\begin{center}
\includegraphics[scale=0.50]{Figure_02.jpg}
\end{center}
\caption[]{Emission spectra of PPO and bis-MSB samples diluted in cyclohexane are compared to a o-PXE based scintillator (3~g/l PPO, 20~mg/l bis-MSB) spectrum measured in a geometry with negligible self absorption.}
\label{Fig2}
\end{figure}

The LS emission spectrum in Figure~\ref{Fig2} shows the basic features of a pure PPO spectrum shifted by 8~nm to longer wavelengths. This wavelength shift can be explained by the combination of a solvent and concentration effect. In the spectrum, there is also some contribution of the bis-MSB (less than 30\%) and a small PXE emission peak around 290~nm. We mainly see the PPO emission when exciting the PXE due to the fact that there is non-radiative energy transfer from the excited PXE to the PPO whereas the energy transfer from PPO to bis-MSB is mainly radiative. Since propagation through the liquid is minimized for the detected light in the front face geometry, radiative transfer is suppressed in this case. The rather poor non-radiative transfer between the two wavelength shifters due to a slight mismatch in the main PPO emssion and bis-MSB absorption could be improved by using a different primary wavelength shifter like e.g.~p-terphenyl~\cite{PXE} which on the other hand has the disadvantage of limited solubility.

In detector simulations it is normally better to use the emission spectrum after non-radiative transfer instead of the ones of the single molecules~\cite{SWD}. Otherwise the simulation might assume radiative energy transfer, which differs from non-radiative transfer in several aspects as e.g.~the efficiency. Non-radiative transfer is much harder to simulate due to contribution of multiple processes and complicated microphysics~\cite{Ent}.

As the light propagates through the liquid and gets absorbed and reemitted by the secondary wavelength shifter, the emission spectrum changes and shifts to longer wavelengths, since the absorbed photons can only be reemitted at equal or lower energies. Figure~\ref{Fig3} shows again the LS spectrum of Figure~\ref{Fig2}. In addition, the spectrum of the same PXE based scintillator mixture is plotted for a case when the light had to pass through the liquid for few millimetres. In this spectrum the PPO component has essentially disappeared and the light is mainly emitted in the bis-MSB region above 400~nm. At longer pathlengths of several centimetres or even metres the spectrum will be even more shifted to longer wavelengths. To predict the response of light sensors such as photomultiplier tubes (PMTs) the wavelength dependent quantum efficiency of the device as well as the light spectrum at the PMT position needs to be known. As we see from Figure~\ref{Fig3}, the distance between the point of light creation and PMT position is relevant for the LS emission spectrum.

\begin{figure}[tb]
\begin{center}
\includegraphics[scale=0.50]{Figure_03.jpg}
\end{center}
\caption[]{Scintillator emission in front face geometry (small self-absorption) and in a 1~cm triangular cell (some self-absorption).}
\label{Fig3}
\end{figure}

Scintillation photons from $350-430$~nm are mainly absorbed by the secondary wavelength shifter as described above. At shorter wavelength the amount of scintillation light is negligible. Nevertheless, the reemission probability below 350~nm is of interest as well, in case of Cherenkov light production. Good modelling of the Cherenkov contribution can be of crucial importance to understand the non-linearity in the energy scale of a detector~\cite{SWD}. Photons emitted after light absorption below 350~nm mainly originate from excited PPO molecules. Here the PPO molecules are excited either directly ($280 - 350~$nm) or indirectly via non-radiative energy transfer from the solvent ($< 280$~nm). Nevertheless, the reemission probability does not only depend on the fluorescence quantum yield of PPO, in particular for wavelengths lower than 300~nm. It is also influenced by the interaction mechanisms between molecules. The reemission probability in the LS has a rising trend with increasing wavelength and converges to the fluorescence quantum yield of PPO before bis-MSB absorption and emission take over~\cite{SWD}. Therefore the quantum yield not only of the secondary, but also of the primary wavelength shifter is an important parameter for the modelling of Cherenkov light. The secondary wavelength shifter is added typically at low concentrations to the scintillator and the emission spectrum is shifted relative to the absorption bands of the other molecules. Therefore, the reemission probability above 350~nm is less affected by self-interactions or solvent effects and can be approximated by the bis-MSB fluorescence quantum yield.

\section{Fluorescence quantum yields}

The literature values of fluorescence quantum yields, in particular for PPO and bis-MSB, vary over a wide range. In many databases the values of Berlman~\cite{Ber} are quoted. However, these values might overestimate the true numbers. For the quantum yield of one of the samples ``Berlman has rather arbitrarily assigned a value of 1.00''~\cite{Dem}. This sample was used as a standard and all other numbers in the Berlman~\cite{Ber} handbook refer to it. The literature values we found for PPO range from 0.71 to 1.00, those for bis-MSB from 0.75 to 0.96 \cite{Ber, XIA, Mas, Joh}.

The reasons why it is so hard to reproduce those results are manifold. The determination of the quantum yield strongly depends on many parameters such as sample concentration, solvent choice, oxygen content, excitation/emission wavelength dependencies, sample geometry, self-absorption or instrumental effects. Since the values are typically determined in relative measurements, a rather large uncertainty on the knowledge of the quantum yield standard also enters the measurement. Therefore it is useful to have several independent measurements using different setups.

For our measurements we diluted the samples in cyclohexane. This rather inert solvent has the advantage of high transparency and optical purity. The PPO and bis-MSB samples were diluted to a concentration of 0.5~mg/l. The absorption spectrum was measured using a Varian Cary400 UV/Vis photospectrometer. To get the emission spectrum a Varian CaryEclipse fluorimeter was used. The measurements were done in a 1~cm quartz cell. Oxygen was removed before the measurements by nitrogen bubbling. After bubbling the cells were closed with a sealed cap. To avoid any influence due to variation in the light intensity at different excitation wavelengths we extracted the quantum yield using the same excitation wavelength in the sample and the standard. As a standard we were using quinine sulfate ($1.5\cdot 10^{-5}$~M) in 0.1~M sulfuric acid.

This quinine sulfate reference fulfills several criteria of a adequate standard as good separation of broad absorption and emission bands, which are in similar wavelength regions as those of the samples. Furthermore, it is rather insensitive to oxygen or concentration quenching and stable in solution. Its quantum yield $\Phi_r$ of 0.55~\cite{Mel} is rather accurately known and constant over the investiagted excitation wavelengths (290 -- 380~nm)~\cite{Gil}. However, there is a temperature dependence of the yield~\cite{Mel}. The measurements were done at a room temperature of about $21 - 25^\circ$C. In this range corrections to the reported quantum yield of the standard should be less than 1~\%.

Since the solvent in the sample and the solvent in the standard are different, a refractive index correction had to be applied. The quantum yields of the samples $\Phi_x$were calculated using
\begin{equation}
\Phi_x = \Phi_r \frac{B_r}{B_x} \frac{I_x}{I_r} \frac{n_x^2}{n_r^2}.
\end{equation}
In this equation $B_x$ and $B_r$ correspond to the fraction of light absorbed in the sample ($x$) and the standard ($r$). The values for $B$ were kept below $20$~\% in sample and reference. The  wavelength integrated intensities of the emission spectra $I_x$ and $I_r$ were corrected for the instrumental wavelength dependent detection efficiency. For the refractive index of cyclohexane we used $n_x = 1.43$ and for the one of the standard $n_r = 1.34$.

The PPO quantum yield was determined at an excitation wavelength of 290~nm and from 300 to 330 nm in 5~nm steps. An average value of 0.842 was found. For bis-MSB the result we obtain is 0.863 which is also an average in a wavelength range from 300 to 380 nm (5~nm steps from 300 -- 350~nm, 10~nm steps from 350 -- 380~nm). Both values are found to be stable within the wavelength range covered by the measurements. The variations from the average value was within 3~\% for all wavelengths. As an additional data point the PPO yield was checked at an excitation wavelength of 350~nm. To obtain a reasonable absorption at this wavelength the PPO-concentration in cyclohexane was increased by more than 2 orders of magnitude to $4.5\cdot 10^{-4}$~M. Even at these high concentrations no indication for any decrease in the quantum yield could be found in the tail of the absorption spectrum.

Since the emission of the bis-MSB light is mainly above 400~nm it would be interesting to test the yield also in that region to get the reemission probability when the photons are self-absorbed. However, the bis-MSB absorption drops rapidly in that area and much higher concentrations are needed to determine those values. This might induce systematic effects influencing the precision of the results. In addition, the quantum yield of the standard is not reliable anymore above 400~nm and one would need to compare it to data taken at lower excitation wavelengths, which requires additional correction factors increasing the uncertainty further. For those reasons we restricted ourselves in the determination of the quantum yields to wavelengths up to 380~nm.

To check systematic effects the measurements were repeated varying some of the relevant parameters. To investigate effects of concentration quenching or self-absorption the yield was determined at lower concentrations (factor 2). In addition, the cell geometry was modified using also a thinner cell of 4~mm length and 1~cm width. The results obtained all agree within a $1 \sigma$ error of 5~\%, which was estimated from the systematic contributions as listed above. In particular, this 5~\% relative error covers the uncertainties related to the precision of the reference yield, knowledge of the refractive indices and quenching effects due to inefficient oxygen removal. After oxygen removal in our samples the quantum yield increased in the case of bis-MSB by 7~\%. A weaker effect was found for PPO.

The ratio of the bis-MSB to PPO quantum yield is 1.02 for our results. This is in reasonable agreement to the ratio of Berlman which is 0.94~\cite{Ber} assuming the Berlman ratio to also have an uncertainty in the 5~\% range. However, our absolute values are significantly lower, which might be explained by the fact that Berlman overestimated the yield of his reference, as already suggested by Demas and Crosby~\cite{Dem}. On the other hand our results are higher than the ones measured by members of the Borexino Collaboration~\cite{Joh, Mas}, but also here we have reasonable relative agreement for the $\Phi_\textrm{\scriptsize bis-MSB}/\Phi_\textrm{\scriptsize PPO}$ ratio. If we compare the results presented here with the values of Xiao et al.~\cite{XIA} we find good agreement for the absolute value of PPO, however, the yield for bis-MSB is about 10~\% lower in our case. Part of the reason could be a strong solvent effect, since the solvent for bis-MSB in reference~\cite{XIA} is linear alkyl benzene (LAB) instead of cyclohexane.

\vspace{0.2 cm}

To study such solvent effects we also did the measurements replacing cyclohexane by \linebreak \mbox{n-dodecane}, LAB and o-PXE. The values in n-dodecane (the refractive index $n_\textrm{\scriptsize dodecane}$ at 405~nm, 18$^\circ$C is 1.42), the main component in the KamLAND~\cite{Kam} and Double Chooz detectors~\cite{LSP}, were within the precision of the measurement consistent with the ones determined in cyclohexane. The aromatic solvent LAB ($n_\textrm{\scriptsize LAB} = 1.49$) which is e.g.~used in the targets of the Daya Bay~\cite{DB} and RENO~\cite{RENO} reactor neutrino experiments has a rather strong self emission when excited below 340~nm making the measurement complicated. Since the PPO absorption is already very weak above 340~nm we only determined the yield for bis-MSB in LAB from 340 to 380~nm. At a bis-MSB concentration of about $1\cdot 10^{-6}$~M the yields were 5 to 10~\% higher in LAB compared to cyclohexane which is consistent with the bis-MSB result reported by Xiao et al.~\cite{XIA}. For o-PXE ($n_\textrm{\scriptsize PXE}$ at 405~nm, 18$^\circ$C is 1.604) which is used in Double Chooz there are similar difficulties at lower wavelengths as for the LAB, so also here the measurement was only done for bis-MSB at 370 and 380~nm. Slightly lower yields were recorded in o-PXE compared to cyclohexane, but the effect was only 5~\% or less and therefore not significant. In reference~\cite{Joh} quite large solvent effects were reported in pseudocumene (1,2,4-trimethylbenzene, PC) which is the basis of the Borexino~\cite{Bor} scintillator and also contained in the KamLAND detector. Compared to the cyclohexane solvent the value increased by 13~\% for the case of PPO and it decreased by 13~\% for bis-MSB.

\vspace{0.2 cm}

In Table~\ref{QY_sum} we summarize all our results obtained in cyclohexane for several solvent or wavelength shifting molecules and other fluorescence materials of interest. Anthracene and 9,10-di-phenylathracene (DPA) are sometimes used as standards in this type of measurements. Both of them have the disadvantage that they are very sensitive to oxygen quenching and that there absorption bands are very narrow and spiky making them sensitive to instrumental parameters as the excitation bandwidth. Besides, the separation of absorption and emission bands is rather poor increasing self-absorption processes inside the sample cell. So we observed in the DPA sample a rather high spread in the quantum yield numbers for the different excitation wavelengths and found a relative difference in the integrated light emission with and without nitrogen bubbling of almost 20~\%.

In many studies anthracene is diluted in ethanol, so we determined the yield also in this solvent ($n = 1.36$). Values were measured between 300 and 370~nm in 5~nm steps and the yields were found to be stable within errors in this region. For anthracene in ethanol ($1.1\cdot 10^{-5}$~M) $\Phi = 0.297\pm0.022$ was found which is in agreement to the literature value of 0.27~\cite{Dem}. Dissolved in cyclohexane ($1.1\cdot 10^{-5}$~M) we get a higher value of $\Phi = 0.323\pm0.021$. This number can be compared to the one in the handbook of Berlman which is 0.36. Our lower result is consistent with the suspicion that all of Berlman's values should be reduced.

\begin{table}[h]
\caption[Composition]{Fluorescence quantum yield of aromatic molecules used in liquid scintillators. All samples were diluted in cyclohexane and measured at room temperature inside a 1~cm fluorescence cell (for DPA, BPO, PBD and butly-PBD the yields determined in the 4~mm cell were used). Besides the range of excitation wavelengths used to determine the quantum yields the wavelengths of the absorption and emission peaks are quoted.} \label{QY_sum}
\begin{center}
\begin{tabular}{lccccc}
Molecule & concentration & exc.~range~[nm] & Quantum yield & abs.~max. & em.~max.\\
\hline
bis-MSB &  $1.6\cdot 10^{-6}$~M & 300 - 380 & $0.863\pm0.043$ & 345~nm & 418~nm \\
PPO	&  $2.3\cdot 10^{-6}$~M &  290 - 330 & $0.842\pm0.042$ & 303~nm & 358~nm \\
Anthracene &  $1.1\cdot 10^{-5}$~M &  300 - 370 & $0.323\pm0.021$ & 357~nm & 400~nm\\
DPA & $1.5\cdot 10^{-6}$~M & 350 - 380 & $0.91\pm0.05$ & 373~nm & 407~nm \\
BPO & $1.7\cdot 10^{-6}$~M & 320 - 350 & $0.91\pm0.05$ & 320~nm & 384~nm \\
PBD & $1.7\cdot 10^{-6}$~M & 300 - 320 &  $0.84\pm0.05$  & 302~nm & 358~nm \\
butyl-PBD & $1.4\cdot 10^{-6}$~M & 300 - 320 &  $0.89\pm0.05$ & 302~nm & 361~nm \\
POPOP &  $1.4\cdot 10^{-6}$~M  & 350 - 380 & $0.90\pm0.05$ & 360~nm & 411~nm \\
o-PXE & $1.2\cdot 10^{-4}$~M & 260 - 280 & $0.33\pm0.03$ & 269~nm & 290~nm \\
PC & $8.3\cdot 10^{-5}$~M & 250 - 270 & $0.41\pm0.04$ & 267~nm & 290~nm \\
LAB & $4.1\cdot 10^{-4}$~M & 250 - 270 & $0.20\pm0.02$ & 260~nm & 284~nm \\
DIN & $4.7\cdot 10^{-5}$~M & 250 - 270 & $0.32\pm0.03$ & 279~nm & 338~nm \\
\hline
\end{tabular}
\end{center}
\end{table}

Among all the molecules studied by Berlman as well as in our studies the highest fluorescence quantum yield is found for DPA. This yield was set to 1.00 in the Berlman book and used as reference for all other molecules. At excitation around the main absorption peak (350~nm to 380~nm), we find a quantum yield of 0.91 measured in a 4~mm fluorescence cell. This means if compared to the results presented here all of Berlman's values should be reduced by 9~\%.

Alternative options for PPO or bis-MSB in organic LS were measured for comparison. The primary wavelength shifters BPO (2-(4-biphenyl)-5-phenyloxazole), PBD (2-(4-biphenyl)-5-phenyl-1,3,4-oxadiazole) and butyl-PBD (2-(4-biphenyl)-5-(4-tert-butyl-phenyl)-1,3,4-oxadiazole) were already investigated in the context of Indium loaded LS \cite{JL}. For PBD and butyl-PBD the reference sample was the diluted PPO-solution which absorbs and emits in very similar wavelength regions. For the primary wavelength shifter candidates the highest values are found for BPO and butyl-PBD. Liquid scintillators containing these molecules typically also provide high light yields~\cite{JL}. Although there is reasonable relative agreement between quantum yield ratios for most molecules investigated in this work with the available literature values of Berlman~\cite{Ber} and most discrepancies in the absolute values could be explained by a common correction factor of about 10~\%, there is a notable difference between the PBD and the PPO quantum yield. We find identical numbers whereas Berlman finds a PBD quantum yield which is 17~\% below the one of PPO.

A possible replacement of the secondary wavelength shifter bis-MSB is POPOP (5-Phenyl-2-[4-(5-phenyl-1,3-oxazol-2-yl)phenyl]-1,3-oxazole) which seems to have a slightly higher yield.

For the solvent molecules in Table~\ref{QY_sum} the estimated relative uncertainty is higher than the one of the wavelength shifters, since absorption and emission bands are shifted to lower wavelengths compared to the standard. In addition, those solvent molecules are very sensitive to oxygen quenching. Whereas several neutrino experiments in the past used PC-based scintillators providing high light yields, more recent experiments prefer high flash-point solvents which are advantageous in terms of safety aspects as o-PXE or LAB. The yield of o-PXE (Dixie Chemicals) was first determined relative to the quinine sulfate reference. The other solvents were then measured relative to PXE. For PC (Aldrich, 98~\%) the highest yield of all tested solvents was found. The LAB (Petresa 550-Q) had significantly lower yields than the other molecules. Another solvent with properties which are promising for large scale neutrino detectors, Diisopropyl naphtalene (DIN, Ruetasolv DI-S, mixture of isomers), has a similar yield as o-PXE.

\section{Conclusion}

For a good modelling of the light propagation in large scale liquid scintillator detectors a detailed knowledge of the light absorption properties, the photon emission behaviour and the fluorescence quantum yields of the wavelength shifters is essential. For the commonly used molecules PPO and bis-MSB the literature values for the quantum yield vary by 30~\%. Therefore new measurements were required with an emphasis on the control over systematic effects as presented in this article. The quantum yields obtained in our studies are $\Phi_\textrm{\scriptsize PPO} = 0.842\pm0.042$ and $\Phi_\textrm{\scriptsize bis-MSB} = 0.863\pm0.043$. Similar values ranging from $84-91$~\% were found in alternative primary and secondary wavelength shifters when dissolved in cyclohexane. Variations of the scintillator solvent can change the yield by about 10~\% in some cases.

Instrumentation and Detectors (physics.ins-det)