High-pressure and doping studies of the superconducting antiperovskite SrPt3P

Abstract: We report the results of our investigation of SrPt3P, a recently discovered strong-coupling superconductor with Tc = 8.4 K, by application of high physical pressure and by chemical doping. We study hole-doped SrPt3P, which was theoretically predicted to have a higher Tc, resistively, magnetically, and calorimetrically. Here we present the results of these studies and discuss their implications.

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

	In 2012, Takayama et al.\cite{Takayama2012} announced their discovery of $\mathrm{SrPt_3P}$, a centrosymmetric material that crystallizes in an antiperovskite structure similar to that of $\mathrm{CePt_3Si}$, with the notable difference that the distorted $\mathrm{Pt_6P}$ octahedral units are arranged antipolar rather than polar\cite{Takayama2012} which leads to centrosymmetry as opposed to non-centrosymmetry in $\mathrm{CePt_3Si}$\cite{Bauer2004}. This ternary platinum phosphide is a strong-coupling superconductor with a $T_c$ of $\sim$ 8.4 K\cite{Takayama2012}. Following the discovery of $\mathrm{SrPt_3P}$, several conflicting theoretical papers studying the electronic, vibrational, and thermodynamic properties of this material were published\cite{Subedi2013,Szczesniak2014, Kang2013,Chen2012}. One suggested that the structure might be tuneable by an appropriate choice of a substitute for phosphorus\cite{Subedi2013}; another suggested that the $T_c$ might increase with hole doping\cite{Nekrasov2012a}.

    Herein we report results of the chemical doping and high pressure effects on the properties of $\mathrm{SrPt_3P}$. In this study, the samples were prepared as follows. Stoichiometric amounts of  platinum powder, phosphorus powder, and strontium pieces were mixed together in a glove box under an argon atmosphere with total $O_2$ and moisture level less than 1 ppm. The mixture was pressed into a small pellet and loaded into a $Al_2O_3$ crucible. The crucible together with the pellet was then sealed in a clean quartz tube under vacuum. The whole assembly was then put inside a tube furnace for reaction. The reaction sequence is the same as reported by Takayama et al\cite{Takayama2012}. X-ray diffraction characterization of the material was carried out using a Panalytical X’pert diffractometer, magnetic measurements using the Quantum Design MPMS, and specific heat measurements using the Quantum Design PPMS. Resistivity measurements were made using a Linear Research LR-400 AC Resistance Bridge operated at 15.9 Hz and a modified probe in the MPMS and high pressure resistivity measurements were made using a BeCu piston-cylinder cell with a quasi-hydrostatic pressure medium, Fluorinert77.

\begin{figure}[b!]
  \centering
    \includegraphics[width=0.5\textwidth]{XRD_indexed}
  \caption{X-ray diffraction pattern for $SrPt_3P$. Unknown peaks are noted with a red mark.}
\end{figure}

The X-ray powder pattern of $\mathrm{SrPt_3P}$ is shown in Fig. 1. The sample is rather pure, as evidenced by the matching of almost all peaks to the structure proposed by Takayama et al. The relatively good quality of the sample was further verified by the rather sharp resistivity drop of the superconducting transition. The width of the resistive transition at 90\% of the drop is $\approx$ 0.31 K, and the diamagnetic shielding is $\approx$1.3 (without demagnetization corrections) as shown in Fig. 2 and its insets. The large difference between the field-cooled (FC) and zero-field-cooled (ZFC) data suggests type II superconductivity with possible strong field pinning.

We began our investigation of the properties of $\mathrm{SrPt_3P}$ by applying high physical pressure to the material.
\begin{figure}[h!]
  \centering
    \includegraphics[width=0.5\textwidth]{Resistivity_Insets}
  \caption{Temperature dependent resistivity and magnetic susceptibility of $SrPt_3P$. Inset (top left): the narrow resistive transition in an enlarged temperature scale from 0 to 20 K. Inset (bottom right): the magnetic susceptibility in ZFC and FC modes}
\end{figure}
\begin{figure}[h!]
  \centering
    \includegraphics[width=0.5\textwidth]{SrPt3P_HighPressure}
  \caption{(a) Temperature dependent resistivity of $\mathrm{SrPt_3P}$ at different pressures. (b) Enlarged temperature scale of resistivity vs. temperature near $T_c$ for $SrPt_3P$ at different pressures (the sequential order of runs indicated by numbers in parenthesis). Insets: in (a) resistivity vs. temperature above $T_c$ for loading and unloading runs, and (b) $T_c$ vs pressure,}
\end{figure}
The temperature dependent resistivity measurements under pressure are shown in Fig. 3. Taking the midpoint of the transition as the superconducting transition, we can clearly see an increase of the superconducting transition from 8.35 K at ambient pressure to a maximum 8.49 K at a pressure of 9.90 kbar (Fig. 3b). Further application of pressure beyond 9.90 kbar appears to lower the transition temperature, with the highest pressure of 17.53 kbar leading to a $T_c$ of 8.47 K. The change in slope of the superconducting transition at high pressure may be related to its width. Measurements were taken at several pressures upon unloading during the course of the experiment; we find that the trend is largely reversible  with a minor split ($<$0.01 K) between the pressure-increasing and pressure-decreasing branches.
\begin{figure}[t]
  \centering
    \includegraphics[width=0.5\textwidth]{131XRD_doping}
  \caption{XRD results for the $SrPt_3P$ and $SrPt_3P_{0.8}Si_{0.2}$ samples. Inset: shift of the (201) XRD peak with increasing Si content.}
\end{figure} Upon unloading to ambient pressure, the $T_c$ does not recover its original value and instead is slightly higher. Irreversible defects may have been introduced to the sample during the course of applying pressure as evidenced by our observation that, upon unloading pressure, the resistivity above $T_c$ is higher than before loading (Fig. 3a, inset). Previous reports characterize $\mathrm{SrPt_3P}$ as a traditional electron-phonon superconductor in the strong-coupling regime\cite{Takayama2012,Subedi2013}; therefore, the increase in $T_c$ upon application of high pressure can be explained by an increase in the characteristic phonon energy $\Theta_D$\cite{Lorenz2005a}. The observed decrease in $T_c$ after the maximum is reached could be viewed as the result of a competition between an increasing Debye temperature and a decreasing density of states at the Fermi level\cite{McMillan1968, Allen1975}.
\begin{figure}[b]
  \centering
    \includegraphics[width=0.5\textwidth]{131_Tc_v_x}
  \caption{$T_c$ vs x for $\mathrm{SrPt_3P_{1-x}Si_x}$}
\end{figure}

Chemical doping was also carried out for the $\mathrm{SrPt_3P}$ compound. We selected Si as the dopant for two reasons: 1) the size of the Si is very close to that of P, minimizing any possible interference due to size difference with the formation of the material; 2) the doping of Si will likely introduce more hole carriers and will allow the study of hole-doping effects theoretically proposed by Nekrasov et al.\cite{Nekrasov2012a}. Fig. 4 shows the results of X-ray diffraction (XRD) with Si doping. The XRD spectrum of the un-doped $\mathrm{SrPt_3P}$ sample is compared with that of the Si-doped sample (in red). It is clear from this comparison that Si can replace P up to 20\% without inducing an impurity phase, as evidenced by the decrease of the lattice parameter show in the inset of Fig. 4. Higher Si-doping in $\mathrm{SrPt_3P_{1-x}Si_x}$ with x = 0.4 and 0.6 has also been carried out. Our XRD results clearly show a phase separation with the impurity phase increasing with x. This can also be seen in the deviation of $T_c$ from linearity with x (Fig. 5).

Temperature dependent magnetization and resistivity measurements were performed on samples with x = 0 and x = 0.2. Figure 6a shows that the onset transition temperature of the parent compound $\mathrm{SrPt_3P}$ is roughly 8.4 K.

\begin{figure}[t!]
  \centering
   \centerline{\includegraphics[width=0.6\textwidth]{MagneticField_Data}}
  \caption{Top: Magnetic susceptibility vs. temperature at different magnetic fields for (a) $\mathrm{SrPt_3P}$, (b) $SrPt_3P_{0.8}Si_{0.2}$; resistivity vs. temperature with different magnetic fields applied for (c) $SrPt_3P$, (d) $SrPt_3P_{0.8}Si_{0.2}$.}
\end{figure}
From our resistivity measurements, we observe suppression of superconductivity with a magnetic field of 3 T to below the lowest temperature measured. The temperature dependent magnetization of x = 0.2 is shown in Fig. 6b with an onset transition temperature of $\sim$7.8 K. We have also performed magnetization measurements on samples with x = 0.4 and x = 0.6, and observed that the $T_c$ decreases to 7.5 K and 7.0 K with increasing Si doping, respectively.

Bulk superconductivity was established in both the parent compound $\mathrm{SrPt_3P}$ and the sample with 20\% Si, $\mathrm{SrPt_3P_{0.8}Si_{0.2}}$, by measuring the specific heat of both the parent compound $\mathrm{SrPt_3P}$ and the sample with 20\% Si, $\mathrm{SrPt_3P_{0.8}Si_{0.2}}$. We apply the $\alpha$-model\cite{Padamsee1973, Johnston2013} to extract further information from the calorimetric data, which extends BCS theory by introducing a variable parameter $\alpha$ = $\Delta_0$/$k_B T_c$. Subtracting the normal state specific heat $\mathrm{C_N}$ (achieved by applying 7T magnetic field) from the superconducting state, we can arrive at a plot of $\Delta$C/T vs. T (Fig. 7). From this plot (knowing that the actual value of $\gamma = \Delta C / T - C_{el}(T)/T$ must lie between the lowest measured value of $\Delta C/T$ and the 0 K $\alpha$-model value) we can estimate the normal state $\gamma$ values for both samples. For the x = 0 sample, we find $\gamma$ = 6.31 $\pm$ 0.10 $mJ/mol\cdot K^2$, and for the x = 0.2 sample we find $\gamma$ = 7.77 $\pm$ 0.73 $mJ/mol\cdot K^2$. We arrive at $\alpha$ = 2.4 for x=0 and $\alpha$ = 2.2 for x=0.2. From this information we can extract $\Delta_0$ = 1.65 meV for $SrPt_3P$, which is agreeable with the report of $\Delta_0$ = 1.58 meV by Khasanov et al.\cite{Khasanov2014a}; we find $\Delta_0$ = 1.28 meV for $SrPt_3P_{0.8}Si_{0.2}$.
\begin{figure}[htp]
  \centering
    \includegraphics[width=0.38\textwidth]{SpecificHeat}
  \caption{$\Delta$C/T vs. T for (a) $SrPt_3P$ and (b) $\mathrm{SrPt_3P_{0.8}Si_{0.2}}$. (c) $C_{el}/\gamma T_c$ vs. $T/T_c$ for both $SrPt_3P$ (broadened) and $\mathrm{SrPt_3P_{0.8}Si_{0.2}}$. The red lines represents the $\alpha$-model, and the blue lines represents the BCS-theory. Inset: $\Delta$C/T vs. T for $SrPt_3P$ with a broadened transition temperature width.}
\end{figure}
It can be seen from the specific heat data that the $T_c$ of the x=0.2 doped sample is broader than that of the undoped sample. For x = 0.2, the width is roughly 1.0 K, whereas for x = 0, the width is roughly 0.5 K. To verify that the difference in $\alpha$ value obtained from our fitting is not an artifact due to this difference between the samples, we artificially broadened the data of the x = 0 sample by assuming a 20\% content with 5\% lower $T_c$. We assume the same temperature dependence below $T_c$ for the lower $T_c$ component. This assumption results in a transition with a width that is the same as the x = 0.2 sample (1.0 K). This modified data is still well-fitted by $\alpha$ = 2.4 (Fig 8a. inset). Comparing the scaled specific heat $C_{el}/\gamma T_c$ vs. $T/T_c$ (Fig. 7c) for the x=0 data and the x=0.2 data with equal transition temperature width, we observe a significant decrease in the specific heat jump $\Delta C/ \gamma T_c$ at $T_c$ from $\approx$ 1.87 for the undoped sample to $\approx$ 1.36 for the x = 0.2 doped sample.

	The apparent increase of $N(\epsilon_F)$ with Si doping is in agreement with previous electronic structure calculations showing the Fermi level located on a negative slope in the density of states\cite{Nekrasov2012a}. The decrease in $T_c$ despite this increase in $N(\epsilon_F)$ supports the suggestion that the $T_c$ does not simply scale with $N(\epsilon_F)$\cite{Takayama2012}. The observed decrease in $\alpha$ could explain the observed decrease in $T_c$, and is in agreement with a previous report that decreasing coupling strength suppressed $T_c$ despite an increase in carrier density\cite{Kasahara2009}.

	In conclusion, we investigated the superconductivity in $\mathrm{SrPt_3P}$ resistively, magnetically, and calorimetrically, including under the application of high physical pressure and with the partial replacement of phosphorus with silicon. The high pressure investigation yielded the result that the superconducting transition temperature increases by a maximum of approximately 0.14 K upon application of pressure up to 9.90 kbar, with an apparent decrease in $T_c$ beyond that pressure. Contrary to the previous predictions \cite{Nekrasov2012a}, hole doping with Si results in a systematic decrease of $T_c$, despite an apparent increase in $\gamma_N$. Our specific heat measurements demonstrate the bulk nature of the superconductivity in the parent compound as well as the silicon doped sample with x = 0.2; furthermore, we observe a decrease in $\alpha$ as well as $\Delta C / \gamma T_c$, implying that a decreasing coupling strength may be responsible for the suppression of $T_c$. We suggest that the observed decrease in coupling strength, combined with the suppression of $T_c$, indicates that the vibrational modes and hence the $T_c$ in this system may be highly sensitive to and tunable by modification of the details of the crystal structure.

\begin{acknowledgements}
The work in Houston is supported, in part, by U.S. Air Force Office of Scientific Research Grant FA9550-09-1-0656, the T.L.L. Temple Foundation, the John J. and Rebecca Moores Endowment, and the State of Texas through the Texas Center for Superconductivity at the University of Houston.

\end{acknowledgements}


Very High Redshift Radio Galaxies

Abstract: High redshift radio galaxies (HzRGs) provide unique targets for the study of the formation and evolution of massive galaxies and galaxy clusters at very high redshifts. We discuss how efficient HzRG samples are selected, the evidence for strong morphological evolution at near-infrared wavelengths, and for jet-induced star formation in the z=3.800 HzRG 4C41.17.

\section{Introduction}

Radio sources are surprisingly effective beacons for identifying
galaxies at extremely high redshifts.  Optical/near--IR campaigns
during the past few years by several groups have resulted in the
discovery of more than 120 radio galaxies at $z > 2$, including 17
with $z > 3$, and 3 with $z > 4$.  At low redshifts powerful radio
galaxies are uniquely identified with massive ellipticals. If this is
true also at high redshift, as seems reasonable given the surprisingly
good Hubble $K-z$ relation for radio galaxies at $0 < z < 4.4$, then
we should be able to systematically study the evolution of massive
elliptical galaxies over large lookback times using samples selected
by their radio emission.  While recently developed techniques of
finding very distant star--forming galaxies ({\it e.g.}\ Steidel {\it et al.~}\ 1996)
are yielding substantial galaxy populations at $z \sim 3$, radio
galaxy samples remain the best means of finding the most massive
galaxies at $z \sim 3$, and even higher.

Hierarchical galaxy formation scenarios suggest that these massive
galaxies are assembled from smaller structures at relatively late
cosmic epochs ({\it e.g.}\ Baron and White 1987).  Observations of HzRGs may
thus provide a unique opportunity to study the beginning of this
process.  Furthermore, when they have been found, HzRGs may also be
used for 'color-dropout' searches of galaxy clusters around them, and
begin cluster evolution studies at very high redshift.

\section{How To Find High Redshift Radio Galaxies}

\subsection{Ultra Steep Spectrum Sources}

It has been known for many years that radio sources with very steep
spectra are mostly associated with faint, distant galaxies ({\it e.g.}\ Tielens
{\it et al.~}\ 1979; Blumenthal \& Miley 1979).  An example of a steep radio
spectrum, for the powerful nearby radio galaxy Cygnus A, is shown in
figure \ref{cyga}. One can also see that the radio spectrum steepens at
higher frequency.  If all powerful radio galaxies have such curving
spectra, due to synchrotron and/or inverse Compton losses, then with
increasing redshift we would observe the steeper parts of their
spectra. Analysis of the 3CR radio sample indeed showed such an
expected $\alpha - z$ correlation (van Breugel and McCarthy 1990; Fig
\ref{3crrest} shows an updated version). A similar result has recently
been reported for a sample of ultra-luminous radio galaxies at high
redshifts (Carilli {\it et al.~} 1998).

\begin{figure}[h]
\centerline{
\psfig{figure=wilfig1.ps,width=5.5cm}
}
\caption{\label{cyga}
\small{The radio spectral energy distribution of Cygnus A, using single dish
measurements from the literature ({\it e.g.}\ Kellermann, Pauliny-Toth and
Williams 1969, and others). Note the spectral steepening with
frequency.
}}
\end{figure}

\begin{figure}[h]
\centerline{
\psfig{figure=wilfig2.ps,width=5.7cm}
}
\caption{\label{3crrest}
\small{Spectral steepening with redshift for all 3CR radio galaxies.
}}
\end{figure}

Thus this ultra-steep spectrum (USS) `red radio colors' selection
technique may be used as a very effective tool for identifying very
high redshift radio galaxies.  The first comprehensive USS search for
HzRGs was begun by Chambers and collaborators and resulted in the
discovery of the $z = 3.800$ radio galaxy 4C41.17 (Chambers, Miley and
van Breugel 1990). For many years this was the most distant galaxy
known.  This relatively bright object, with its extended {\rm\,Ly-$\alpha$~} halo,
radio-aligned UV continuum, and evidence for star formation (see
below) in many ways stimulated the continued search for more USS
HzRGs.

With the advent of several new, deep radio surveys it is now possible
to define much better USS samples which reach 10 -- 100 times fainter
flux density levels than previous samples, and which allow much more
accurate spectral index determinations. We have used several such
surveys to define USS samples for study in the northern and southern
hemispheres.  Our primary northern hemisphere sample uses the WENSS
325 MHz survey (Rengelink {\it et al.~} 1997) together with the NVSS 1.4 GHz
survey (Condon {\it et al.~} 1998), to define the spectral indices, and
the FIRST 1.4 GHz survey (Becker {\it et al.~} 1995), to obtain radio maps.
In the southern hemisphere we use the low frequency Texas 365 MHz
(Douglas {\it et al.~} 1996; Dec $> -35$\ifmmode {^{\circ}}\else {$^\circ$}\fi) and MRC 408 MHz (Large {\it et al.~} 1981;
Dec $< -35$\ifmmode {^{\circ}}\else {$^\circ$}\fi) surveys, in combination with the NVSS and PMN 5 GHz
(Griffith \& Wright 1993; Dec $< -35$\ifmmode {^{\circ}}\else {$^\circ$}\fi) surveys.

The sources in our southern hemisphere sample all have spectral
indices $\alpha < -1.2$.  We have obtained radio images of all our
southern $USS$ targets with the Australia Telescope and VLA so that we
can make accurate identifications. A deep near-infrared identification
program of this sample was begun with the CTIO 4m telescope.
Spectroscopic observations of some objects from this sample with the
ESO 3.6 m have shown that USS sample is indeed extremely efficient at
finding HzRGs, and have already resulted in the discovery of the most
distant radio galaxy in the Southern Hemisphere known to date (TN
J1338$-$1942 at $z=4.13$; De Breuck {\it et al.~}, these proceedings).  A
summary of the current status of HzRG identification programs is given
in Table 1.

\begin{table}
\begin{center}
\begin{tabular}{| c || c || c | c | c || c |}\hline
\small
Sample & Definition & \multicolumn{3}{c||}{Known Redshift} & Unknown \\
\hline
 & & $z<2$ & $2<z<3$ & $z>3$ &  \\
\hline
3CR & S$_{178} >$ 10 Jy & 99.5 \% & 0.5 \% & 0 \% & 1\% \\
MRC/1Jy & S$_{408} >$ 0.95 Jy & 93.4 \% & 5.9 \% & 0.7 \% & 41 \% \\
4C USS & $\alpha_{\small 178}^{\small 1414} < -1.0$ & 53 \% & 35 \% & 12 \% & 50 \% \\
New USS & $\alpha_{\small 325}^{\small 1400} < -1.3$ & 35 \% & 35 \% & 30 \% & 45\% (faint) \\
\hline
\end{tabular}
\end{center}
\caption{
\small{The HzRG content for four samples of radio sources as a function of z: 3CR (Spinrad, private communication), MRC/1Jy (McCarthy {\it et al.~}\ 1996), 4C (Chambers {\it et al.~}\ 1996), and our new USS sample.}}
\end{table}
\normalsize

\begin{figure}[h]
\centerline{
\psfig{figure=wilfig3.ps,width=4.8cm}
}
\caption{\label{mphis}
\small{Hystogram of southern USS sources at DEC $<$ -35\ifmmode {^{\circ}}\else {$^\circ$}\fi
from the MRC and PMN surveys.
In total 76 out of 2915 sources have $\alpha^{5 GHz}_{408 MHz} < -1.20$.
}}
\end{figure}

\subsection{The $K-z$ Diagram}

Around the time that 4C41.17 was discovered another systematic search
for HzRGs also proved to be successful.  This method uses the
magnitude -- redshift relationship at infrared wavelengths, the $K-z$
`Hubble' diagram, for powerful radio galaxies galaxies to find
promising targets. This resulted in the discovery of the $z = 3.395$
radio galaxy B2 0902+34 (Lilly 1988). Indeed, the $K-z$ diagram has
since proved to be a very reliable, if little understood,
tool. Despite strong morphological evolution seen at near-infrared
wavelengths, discussed below, the most powerful radio sources continue
to follow the $K-z$ relationship even to z $\sim$ 4.4 (Fig
\ref{kzdiag}).  It is our experience that, by combining the USS and
$K-z$ techniques and using near-infrared identifications of USS
sources which are unidentified on UKST/POSS2 plates, one is virtually
guaranteed to be successful in HzRG hunts.  Without near-IR
identifications and photometry approximately 2/3 of the galaxies are
at $z > 2$. With near-IR identifications and photometry one can select
the redshift range one wishes to study and to date all galaxies
observed at Keck were found to be at the redshift predicted from the
$K-z$ diagram, at least for $z \mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$} 4.4$.

\begin{figure}[h]
\centerline{
\psfig{figure=wilfig4.ps,width=7.5cm}
}
\caption{\label{kzdiag} \small{Hubble K-z diagram for the 3C and 6C
surveys (Eales {\it et al.~} 1997), and NIRC data for HzRGs from van Breugel
{\it et al.~}\ (1998). The new USS $K$-band detections and very high redshift
candidates are shown by horizontal bars at the predicted redshift. The
horizontal arrow indicates the range of predicted redshift for our
faintest $K$-band object. The open star represents TN J1338-1942 (De
Breuck {\it et al.~}, this volume). All magnitudes are corrected to a 64 kpc
metric diameter, assuming H$_{0}$ = 50, q$_0$ = 0.  The inset shows a
blowup of the NIRC data in the $1.5 < z < 5.5$ and $17.5 < K < 21$
region.}}
\end{figure}

\section{Morphological evolution of HzRGs}

Near-infrared images obtained with the W.\ M.\ Keck I Telescope of
HzRGs with $1.9 < z < 4.4$ show strong morphological evolution at {\it
rest--frame optical} ($\lambda_{\rm rest} > 4000$ \AA) wavelengths (van
Breugel {\it et al.~} 1998; Fig \ref{morphev}).  At the highest redshifts, $z
> 3$, the rest--frame visual morphologies exhibit structure on at
least two different scales: relatively bright, compact components with
typical sizes of $\sim$10~{\rm\,kpc}\ surrounded by large--scale ($\sim$ 50
$-$ 100 {\rm\,kpc}) diffuse emission. The brightest components are often
aligned with the radio sources, and their {\it individual}
luminosities are $M_B \sim -20$ to $-22.$ For comparison,
present--epoch L$_\star$ galaxies and, perhaps more appropriately,
ultra-luminous infrared starburst galaxies, have, on average, $M_B \sim
-21.0$. The {\it total, integrated} rest--frame B--band luminosities
are $3 - 5$ magnitudes more luminous than present epoch $L_\star$
galaxies.

At lower redshifts, $z < 3$, the rest--frame optical morphologies
become smaller, more centrally concentrated, and less aligned with the
radio structure.  Galaxy surface brightness profiles for the $z < 3$
HzRGs are much steeper than those of at $z > 3$.  We attempted to fit
the $z < 3$ surface brightness profiles with a de Vaucouleurs
r$^{1/4}$ law and with an exponential law, the forms commonly used to
fit elliptical and spiral galaxy profiles, respectively.  We
demonstrate the fitting for our best resolved object at $z < 3$, 3C
257 at $z = 2.474$ (Fig \ref{morphev}).  Within the limited dynamical
range of the data, both functional forms fit the observed
profiles---neither is preferred.  Interestingly, despite this strong
morphological evolution the $K - z$ `Hubble' diagram for the most
luminous radio galaxies remains valid even at the highest redshifts,
where a large fraction of the K-band continuum is due to a
radio--aligned component.

\begin{figure}[h]
\centerline{
\psfig{figure=wilfig5.ps,width=7.8cm}
}
\caption{\label{morphev}
\small{Selected near--IR images of HzRGs, presented in order of decreasing
redshift, and the surface brightness profile of 3C257.
}}
\end{figure}

\section{Jet-induced star formation in 4C41.17: HST imaging}

Deep HST images of 4C41.17 at `R'-band (rest-frame UV) and in {\rm\,Ly-$\alpha$~}
show that the 4C41.17 system consists of two components: 4C41.17-North
with a bright string of UV knots and {\rm\,Ly-$\alpha$~} emission along the radio
axis, and 4C41.17-South with several much fainter UV knots,
distributed in random fashion throughout a low surface brightness
halo. The brightest radio knot in 4C41.17-North is associated with the
brightest UV knot and arc-like {\rm\,Ly-$\alpha$~} emission. One of the field objects
in the HST images was also seen at near--IR (Graham {\it et al.~} 1994; object
\# 16) and radio wavelengths (Carilli {\it et al.~} 1994). This was used to
align the HST and radio frames with an estimated relative accuracy of
$\sim$ 0.1\ifmmode {^{\prime\prime}~}\else $^{\prime\prime}~$\fi.  The central, radio--aligned UV and {\rm\,Ly-$\alpha$~} emission
is shown in Figure \ref{4c4117hst} with the 0.21\ifmmode {^{\prime\prime}~}\else $^{\prime\prime}~$\fi resolution
radio X-band image from Carilli {\it et al.~} overlaid. Figure
\ref{4c4117smooth} shows the HST continuum image of the entire 4C41.17
system smoothed to 0.3\ifmmode {^{\prime\prime}~}\else $^{\prime\prime}~$\fi resolution to enhance low surface
brightness features.

\begin{figure}[ht]
\centerline{
\psfig{figure=wilfig6.ps,width=7.7cm}
}
\caption{\label{4c4117hst}
HST WFPC2 images of 4C41.17-North with the X-band map of Carilli {\it et al.~}
1994 superimposed.
}
\end{figure}

Downstream from the bright radio/UV/{\rm\,Ly-$\alpha$~} knot the radio source curves
towards a faint, very steep spectrum NE lobe (see Carilli {\it et al.~} 1994),
while upstream from this knot, towards the central radio core, the UV
continuum appears edge--brightened. These morphological features
suggest a strong interaction between the radio jet and dense ambient
gas and, in fact, are as expected in jet--induced star formation
models where sideways shocks induce star formation in the dense medium
of forming galaxies ({\it e.g.}\ De Young 1989; Bicknell {\it et al.~} 1998).

The star formation rates ($SFR$) of the various components as
deduced from the rest--frame UV HST photometry range from 5 -- 40
\ifmmode {\rm\,M_\odot}\else ${\rm\,M_\odot}$\fi/yr for the knots in 4C41.17-North, to 5 -- 10 \ifmmode {\rm\,M_\odot}\else ${\rm\,M_\odot}$\fi/yr in
4C41.17-South.  Here we have assumed L$_{1500\AA} \sim 10^{40.1}$
\ifmmode {\rm\,erg\,s^{-1}}\else ${\rm\,erg\,s^{-1}}$\fi\AA$^{-1}$ for a $SFR$ = 1 \ifmmode {\rm\,M_\odot}\else ${\rm\,M_\odot}$\fi/yr (Conti {\it et al.~} 1996) and no
dust reddening. The derived values are surely lower limits, given the
detection of dust at sub---mm wavelengths in 4C41.17 by Dunlop {\it et al.~}
(1994).

The entire 4C41.17 system is embedded in a common halo of diffuse, low
surface brightness emission which extends over a very large area of
54$h_{50}^{-1}$~kpc $\times$ 76$h_{50}^{-1}$~kpc (5\ifmmode {^{\prime\prime}~}\else $^{\prime\prime}~$\fi $\times$
7\ifmmode {^{\prime\prime}~}\else $^{\prime\prime}~$\fi). This includes a faint region, 4C41.17-South, with half a
dozen compact knots distributed in random fashion. Spectroscopic
observations have shown that 4C41.17-South is indeed at the same
redshift as 4C41.17-North (Dey {\it et al.~} 1999).  The range of UV
luminosities and $SFR$ rates for the individual knots in
4C41.17-South is lower than in 4C41.17-North, and very similar to the
`normal' (radio-quiet) Lyman-break galaxies discovered by Steidel {\it et al.~}
(1996). The random distribution and on average lower $SFR$ in the
4C41.17-South knots suggests that star formation here is unaided by
bowshocks from the radio jet.  The total star formation rate, integrated
over the entire 4C41.17 system and including the low surface brightness
emission, is $\sim$ 660 \ifmmode {\rm\,M_\odot}\else ${\rm\,M_\odot}$\fi/yr.  Of this perhaps as much as 2/3 of
the star formation may be occurring in the inter-knot regions. If the
total star formation would continue at this rate for $2 \times 10^{8} -
2 \times 10^{9}$ yrs an entire massive elliptical galaxy of $10^{11}
\ifmmode {\rm\,M_\odot}\else ${\rm\,M_\odot}$\fi - 10^{12} \ifmmode {\rm\,M_\odot}\else ${\rm\,M_\odot}$\fi$ might be assembled between $z \sim 4$ and $z
\sim 2.5$, consistent with the morphological evolution for HzRGs seen
in the near--IR Keck observations.

\begin{figure}[ht]
\centerline{
\psfig{figure=wilfig7.ps,width=10.5cm}
}
\caption{\label{4c4117smooth}
Smoothed version of the HST WFPC2/F702W image
showing the entire 4C41.17 system, including the clumpy companion system
4C41.17-South.
}
\end{figure}

\section{Jet-Induced Star Formation In 4C41.17: Keck Spectroscopy}

Deep spectropolarimetric observations with the Keck II telescope by
Dey {\it et al.~} (1997) have provided strong evidence in support of the
jet--induced star formation model for 4C41.17-North suggested above on
the basis of the HST and radio morphologies. These observations showed
that the bright, radio--aligned rest--frame UV continuum is unpolarized
($P_{UV}(2\sigma) < 4\%$).  This implies that scattered AGN light,
which is generally the dominant contributor to the rest-frame UV
emission in $z\sim 1$ radio galaxies, is unlikely to be a major
component of the UV flux from 4C 41.17. Instead, the total light
spectrum shows absorption lines and P--Cygni--like features that are
similar to those detected in the spectra of the recently discovered
population of star forming galaxies at slightly lower ($z\sim2-3$)
redshifts (Fig \ref{4c4117spec}).  The detection of the S V$\lambda 1502$ stellar
photospheric absorption line, the shape of the blue wing of the Si IV
profile, the unpolarized continuum emission, the inability of other
AGN--related processes to account for the UV continuum flux, and the
overall similarity of the UV continuum spectra of 4C 41.17 and the
nearby star forming region NGC 1741B strongly suggest that the light
from 4C 41.17 is dominated by young, hot stars.  The presence of
radio--aligned features in many of the $z > 3$ HzRGs suggests, by
analogy to 4C41.17, that jet--induced star formation may be a common
phenomenon at these very high redshifts.

\begin{figure}[ht]
\centerline{
\psfig{figure=wilfig8.ps,width=13cm}
}
\caption{\label{4c4117spec}
Keck spectrum (Dey {\it et al.~} 1997) of the radio-aligned component
4C41.17-North, compared with a UV spectrum of the Wolf-Rayet starburst
galaxy NGC 1741 (Conti {\it et al.~} 1996).
}
\end{figure}

\noindent {\bf Acknowledgments}

{\small Part of the work described here was performed while WvB was on
sabbatical leave during January 1997 -- April 1997 at the
Anglo-Australian Observatory, the Australian National Telescope
Facility, and the Mount Stromlo and Sidings Springs Observatories. He
appreciates the support provided by these institutes. He thanks the
Kookaburra's for their wake--up calls, and his Australian colleagues
for their warm hospitality and invigorating discussions, with special
thanks to Drs.\ J.\ Bland--Hawthorn, G.\ Bicknell, M.\ Dopita, and R.\
Sutherland.  We gratefully acknowledge S.\ Rawlings for advance
information regarding 6C~0140+326 and 8C~1435+635 and his work on the
$K - z$ diagram, and C.\ Carilli for providing high quality radio
images of 4C41.17 which allowed to improve on the relative
radio/optical astrometry for this source. The research by WvB, CDB and
AS at IGPP/LLNL is performed under the auspices of the US Department
of Energy under contract W--7405--ENG--48.}

Very High Redshift Radio Galaxies


Large magneto-thermal effect and the spin-phonon coupling in a parent insulating cuprate Pr_{1.3}La_{0.7}CuO_4

Abstract: The magnetic-field (H) dependence of the thermal conductivity \kappa of Pr_{1.3}La_{0.7}CuO_4 is found to show a pronounced minimum for in-plane fields at low temperature, which is best attributed to the scattering of phonons by free spins that are seen by a Schottky-type specific heat and a Curie-Weiss susceptibility. Besides pointing to a strong spin-phonon coupling in cuprates, the present result demonstrates that the H-dependence of the phonon heat transport should not be naively neglected when discussing the \kappa(H) behavior of cuprates, since the Schottky anomaly is ubiquitously found in cuprates at any doping.

\section{Introduction}

To elucidate the nature of the superconducting ground state, the
quasiparticle (QP) transport properties are widely studied in
high-$T_c$ cuprates at low temperature ($T$), and the discussions
are often based on the magnetic-field ($H$) dependence of the
thermal conductivity $\kappa$.\cite{Krishana, Aubin, Ando1, Chiao,
Sun1, Ando2, Hawthorn} The most striking findings include a
``plateau" in the $\kappa(H)$ isotherm at low
temperature,\cite{Krishana, Aubin, Ando1} an increase in the QP
population due to the supercurrents around vortices,\cite{Chiao}
and the magnetic-field-induced QP localization in underdoped
cuprates.\cite{Sun1, Ando2, Hawthorn} A common assumption in these
previous works is that the phononic thermal conductivity
$\kappa_p$ is {\it independent} of $H$ and all the $H$-dependence
of $\kappa$ is of electronic origin. This assumption may seem
reasonable, because the vortex scattering of phonons (which is
important in conventional superconductors) somehow becomes
insignificant at low temperature in cuprates.\cite{Franz} However,
phonons could interact also with magnetic excitations and be
scattered, which would introduce some $H$-dependence in
$\kappa_p$. Therefore, it is important to elucidate how the
low-$T$ phonon heat transport is affected by magnetic excitations
in cuprates, particularly in the superconducting doping regime.
Remember, even though antiferromagnetic (AF) magnon excitations
are likely to be irrelevant in the superconducting regime,
Schottky anomalies in the specific heat, which are usually
attributed to excitations of localized free spins, are commonly
found in superconducting cuprates.\cite{Moler, Revaz, Wright,
Chen, Brugger, Hien, Schottky}

To gain insight into the possible $H$-dependence of $\kappa_p$, it
is useful to study parent insulating cuprates, where it is
expected that the electronic contribution is absent and $\kappa_p$
can be unambiguously studied; however, it has been discussed that
in AF insulating materials the magnon excitations can also act as
heat carriers or scatterers,\cite{Sun2, Hess, Hofmann, Takeya}
which brings complications to the analysis of the $\kappa(H)$ data
of insulating samples. In particular, Jin {\it et al.}\cite{Jin}
recently reported a large magnetic-field-induced enhancement of
$\kappa$ in Nd$_2$CuO$_4$ (NCO), a parent insulating cuprate in
the tetragonal $T'$-phase, and they proposed that the increase in
$\kappa$ gives evidence for the magnon heat transport becoming
active due to the closing of a magnon gap in the magnetic field.
Note that in NCO the magnons for the Cu spins should not
contribute to carrying heat because of a large ($\sim$5 meV) gap,
and thus it must be the Nd magnons that is relevant;\cite{Li} the
existence of the Nd-magnon contribution to the heat transport has
been backed up by a more recent study by Li {\it et al.}\cite{Li}
Considering these findings, to unambiguously investigate the $H$
dependence of $\kappa_p$, we turn to the
Pr$_{1.3}$La$_{0.7}$CuO$_4$ (PLCO) system where there is
essentially no magnetic contribution from the rare-earth ions.
Although it was found that a small magnetic moment is induced in
Pr$^{3+}$ ions (whose ground state is singlet and would be
effectively nonmagnetic),\cite{Sumarlin, Lavrov1} in
Pr$_{1.3}$La$_{0.7}$CuO$_4$ as much as 35\% of the Pr$^{3+}$ sites
are diluted by nonmagnetic La$^{3+}$ ions and thus the Pr magnons,
even if exist, should be strongly damped and give negligible
contribution to the heat transport. Also, the single crystals of
PLCO were found to show very clean phonon heat transport ({\it
e.g.}, $\kappa_p$ of PLCO is up to 5 times larger at 20 K compared
to La$_2$CuO$_4$,\cite{Sun3}) so that exotic scattering mechanisms
of phonons are expected to be rather easily seen in PLCO without
being masked by scatterings by disorder.

In this paper, we report detailed thermal conductivity, specific
heat, and magnetic susceptibility measurements of high-quality
Pr$_{1.3}$La$_{0.7}$CuO$_4$ (PLCO) single crystals. At low
temperature, the $\kappa(H)$ curve shows a striking dip feature,
which is different from previous data on similar $T'$-phase
compounds, Pr$_{2}$CuO$_4$ (PCO) and NCO.\cite{Jin,Sales} Our new
result on $\kappa(H)$ indicates that excitations of free spins,
evidenced by the specific heat and the magnetic susceptibility
data, are responsible for the peculiar $H$ dependence of the
phonon heat transport. Since the Schottky-type specific heat that
is indicative of the existence of free spins is found almost
ubiquitously in high-$T_c$ cuprates,\cite{Moler, Revaz, Wright,
Chen, Brugger, Hien, Schottky} the present result suggests that
the phonon heat transport can be $H$ dependent at low temperature
in a wide range of cuprate samples and one should be careful upon
discussing the QP properties based on $\kappa(H)$ behavior.
Physically, the present result points to a strong spin-phonon
coupling in cuprates, which might bear intriguing implications on
the mechanism of high-$T_c$ superconductivity.

\section{Experiments}

High-quality Pr$_{1.3}$La$_{0.7}$CuO$_4$ single crystals are grown
by the traveling-solvent floating-zone technique in flowing
oxygen.\cite{Sun3} The partial substitution of La for Pr not only
stabilizes the crystal growth\cite{Sun3, Fujita} but also disturbs
the possible Pr magnons and makes them irrelevant as heat
carriers, which is useful for the purpose of this work. For the
heat transport measurements, the crystals are cut into rectangular
platelets with a typical size of $2.5 \times 0.5 \times 0.1$
mm$^3$ (where the edges are made to be parallel to the
crystallographic axes within 1$^{\circ}$) and annealed in flowing
Ar to remove excess oxygen.\cite{Sun3}

The $a$- and $c$-axis thermal conductivities ($\kappa_a$ and
$\kappa_c$) are measured by a steady-state technique in a $^3$He
refrigerator.\cite{Sun1} Two RuO$_2$ chip sensors are used for the
``top" and ``bottom" thermometers on the sample. For the
magnetic-field dependent $\kappa(H)$ measurements, the base
temperature of heat sink is precisely controlled by an ac
resistance bridge (Linear Research LR-700) within 0.01\% accuracy
with another RuO$_2$ thermometer or a Cernox thermometer (for
above 10 K), and the magnetic-field-dependence of the thermometers
were carefully calibrated beforehand using a SrTiO$_3$ capacitance
sensor and a high-resolution capacitance bridge (Andeen Hagerling
2500A), which allow to keep the temperature unchanged with a very
high precision ($\pm$ 0.1 mK at 0.3 K and $\pm$ 1 mK above 1 K)
during the magnetic field sweeps. The relative uncertainty of the
$\kappa(H)$ measurement was confirmed to be typically less than
1\%. Specific heat and magnetic susceptibility measurements are
carried out using a Physical Properties Measurement System and a
SQUID magnetometer (Quantum Design), respectively.

\section{Results and discussions}

\subsection{Thermal conductivity}

\begin{figure}
\includegraphics[clip,width=7.5cm]{fig1.eps}
\caption{Temperature dependences of $\kappa_a$ and $\kappa_c$ of
PLCO single crystals in zero field. Inset: $\kappa$ vs $T^3$ plot
of the low $T$ data below 0.43 K. Thin solid lines show the $T^3$
dependence expected in the boundary-scattering regime.}
\end{figure}

Figure 1 shows the temperature dependences of $\kappa_a$ and
$\kappa_c$ of PLCO in zero field down to 0.30 K, where the data
above 4 K are a reproduction from Ref. \onlinecite{Sun3}. As
already discussed, one can assume that both $\kappa_a$ and
$\kappa_c$ show purely phononic thermal conductivity in this
insulating compound, and the large magnitude of the phonon peak
evidences a high quality of our PLCO crystals.\cite{Sun3} The
inset of Fig. 1 shows that the boundary-scattering regime of
phonons (where $\kappa_p \sim T^3$) is achieved below $\sim$ 0.35
K.

\begin{figure}
\includegraphics[clip,width=8.5cm]{fig2.eps}
\caption{(a)--(e) Magnetic field dependences of $\kappa$ in PLCO
single crystals for various configurations of the heat current
$J_H$ and the magnetic field. (f) Comparison of $H_{\rm dip}$
(where the $\kappa(H)$ curve shows a minimum) for $J_H \parallel
a$, $H \parallel$ [110] to the spin-reorientation field $H_{\rm
sr}$ for the same field direction (Plakhty {\it et al.}, Ref.
\onlinecite{Plakhty}).}
\end{figure}

Figure 2 shows our main result, the peculiar magnetic-field
dependences of $\kappa$ at temperatures from 0.36 to 18 K, for
several different configurations of heat current and magnetic
field. Notably, a pronounced dip in the $\kappa(H)/\kappa(0)$
curves is observed for in-plane fields at low tempratures. At
first sight, this may seem a typical $\kappa(H)$ behavior in a
spin-flop system where magnons act as phonon
scatterers:\cite{Buys} The zero-field $\kappa$ is purely phononic,
because there are no low-energy magnons due to the anisotropy gap
in the magnon dispersion; when the magnetic field is applied, the
Zeeman energy causes the magnon excitations to become gapless at
the spin-flop transition field $H_{\rm sf}$, but the gap opens
again at higher field; consequently, the magnon scattering of
phonons is the strongest at $H_{\rm sf}$ where the magnons are the
most populated, and causes a dip in the $\kappa(H)$
curve.\cite{Buys}

In the $T'$-phase compounds of NCO, PCO, and PLCO, the magnetic
structures are essentially identical and spin-flop (or
spin-reorientation) transition of the long-range-ordered Cu spins
has been studied by the neutron scattering:\cite{Sumarlin,
Lavrov1, Skanthakumar, Plakhty, Plakhty2, Petitgrand1,
Petitgrand2, Cherny} The Cu spins form a noncollinear structure
with the spin easy axes of [100] and [010] in zero-field;
depending on the magnetic-field direction, a first-order
transition (for $H \parallel$ [100]) or a continuous
spin-reorientation transition (for $H \parallel$ [110]) is
induced, and in both cases the high-field state has a collinear
spin structure. If one tries to interpret our $\kappa(H)$ data in
terms of this well-studied spin reorientation, several problems
become apparent: Most notably, the characteristic field for the
spin reorientation, $H_{\rm sr}$, is known to decrease upon
heating, while the ``dip" field, $H_{\rm dip}$, where the
$\kappa(H)$ curve shows a minimum, increases upon heating. We show
quantitative comparison between $H_{\rm sr}$ and $H_{\rm dip}$ for
$H \parallel$ [110] in Fig. 2(f), where it becomes clear that
$H_{\rm dip}$ is not likely to be related to the spin
reorientation. Also, $H_{\rm dip}$ is observed to be very close
for $H \parallel$ [100] and $H \parallel$ [110], while $H_{\rm
sr}$ is much larger for $H \parallel$ [100] than for $H
\parallel$ [110].\cite{Lavrov1, Plakhty2, Cherny} Furthermore,
the $H$-dependence of $\kappa_c$ is very similar to that of
$\kappa_a$, which is difficult to understand in the spin-flop
scenario because the $c$-axis phonons can hardly couple with the
magnons of the two-dimensional spin system. Therefore, the
peculiar $\kappa(H)$ behavior in PLCO is clearly {\it not} due to
the magnon scattering of phonons associated with the spin
reorientation.

It is known from previous neutron experiments\cite{Sumarlin,
Lavrov1} that the exchange field of the ordered Cu spins induces a
small ($\sim$ 0.1 $\mu_B$) ordered moment on Pr$^{3+}$ ions. One
may expect that the Pr magnons should contribute to the
$\kappa(H)$ behavior. Since the weak Pr moments are induced by the
interaction between Cu and Pr ions, the Pr magnons, if exist and
play a role (as either heat carriers or scatterers) in determining
the $\kappa(H)$ behavior, should affect $\kappa(H)$ in a way that
corresponds to the spin reorientation transition. Therefore, the
lack of correspondence between $H_{dip}$ and $H_{sr}$ means that
not only Cu magnons but also Pr magnons are irrelevant to the
observed $\kappa(H)$ behavior, which is likely due to the
considerable dilution of the Pr sites by nonmagnetic La$^{3+}$
ions.

\subsection{Specific heat}

\begin{figure}
\includegraphics[clip,width=8.5cm]{fig3.eps}
\caption{(a) Temperature dependences of the specific heat of a
PLCO single crystal in various magnetic fields along [110]. (b)
Schottky anomaly for $H \parallel$ [110], obtained by subtracting
the phonon contribution from the raw data in panel (a). (c)
Schottky anomaly at 4 T for different orientations of $H$,
together with the theoretical curves expected for the Schottky
anomaly from 1\% of $s$ = 1/2 moments with the anisotropic
Land$\acute{e}$ factor $g$ = 2 (for [100] and [110]) and $g$ =
1.25 (for [001]). (d) Magnetic field dependences of the peak
temperature in the Schottky anomaly for the three directions.}
\end{figure}

To investigate the actual scatterers of phonons and the origin of
the peculiar $\kappa(H)$ behavior, the specific heat of PLCO is
measured at low temperature using a large single crystal ($\sim$
34 mg, Ar-annealed) from the same batch. As shown in Fig. 3, there
is a pronounced $H$-dependent enhancement, whose temperature and
field dependences are fully consistent with the two-level Schottky
anomaly due to free $s$ = 1/2 moments\cite{Schottky, Moler}
\begin{equation} C_{Schottky}(T,H) = n(\frac{g\mu_BH}{k_BT})^2
\frac{e^{g\mu_BH/k_BT}}{(1+e^{g\mu_BH/k_BT})^2}, \label{Schottky}
\end{equation}
where $g$ is the Land$\acute{e}$ factor. The concentration of free
spins is $n/R$, where $R$ is the universal gas constant. In our
data, the peak position of the Schottky anomaly, $T_{peak}$,
increases linearly with $H$ [see Fig. 3(d)], which precisely
follows Eq. (\ref{Schottky}) and allows us to determine the $g$
factor. Moreover, fitting of the data to Eq. (\ref{Schottky}) [see
Fig. 3(c)] gives the concentration of the free spins to be $\sim$
1\%, which is too large to be associated with impurities and
suggests the existence of some intrinsic mechanism to produce free
moments.

\subsection{Magnetic susceptibility}

\begin{figure}
\includegraphics[clip,width=8.5cm]{fig4.eps}
\caption{(a) Temperature dependences of susceptibility $\chi_a$
and $\chi_c$ of a PLCO crystal in oxygenated (as-grown) and
reduced (Ar-annealed) states. Panels (b) and (c) magnify the
low-$T$ part; solid lines show fittings with the formula $\chi(T)
= \chi_0 + C/(T-\theta)$, where the first and second terms are
$T$-independent background and additional Curie-Weiss
contribution, respectively; the $C$ parameter for the reduced
sample is consistent with $\sim$ 1\% of $s$ = 1/2 free spins with
$g$ = 2 ($g$ = 1.25) for $H \parallel a$ ($\parallel c$).}
\end{figure}

Another evidence for the existence of free moments in PLCO is the
magnetic susceptibility ($\chi$) data. Figure 4 shows the
temperature dependences of $\chi$ for $H \parallel a$ and $H
\parallel c$ ($\chi_a$ and $\chi_c$) in 0.5 T field. Our data are
essentially consistent with those for PCO single
crystals,\cite{Hundley} including the low-$T$ upturn that is
observed in both $\chi_a$ and $\chi_c$. As theoretically discussed
by Aharony's group,\cite{Sachidanandam} the magnetic
susceptibility of R$_2$CuO$_4$ (R = Nd, Pr, or Sm) is governed by
the magnetic states of rare-earth ions. The ground multiplet of
Pr$^{3+}$ ion is $^3H_4$ in the tetragonal crystalline electric
field (CEF), so that Pr$^{3+}$ ions have a singlet ground state
with a well separated ($\sim$ 18 meV) first-excited
state.\cite{Sachidanandam} Therefore, the low-$T$ susceptibility
is predicted to be $T$-independent,\cite{Sachidanandam} and the
observed low-$T$ upturn, which is well described by the
Curie-Weiss law, must come from some free moments not considered
in the mean-field theory. Note that the sensitivity of the low-$T$
upturn in $\chi$ to the oxygen treatment, shown in Fig. 4,
suggests that the free moments do not originate from impurity
ions, but are produced by some mechanism inherent to PLCO.
Candidates for such a mechanism include an appearance of Pr$^{4+}$
ions in some particular crystal environment or some ``free"
Cu$^{2+}$ spins possibly produced in the AF domain
boundaries.\cite{Vaknin, Lavrov3} Incidentally, the analyses of
the Schottky anomaly and the Curie-Weiss contribution are
consistent with the same anisotropic $g$ factor, which gives
confidence that the two phenomena are of the same origin.

\subsection{Phonon scattering by paramagnetic moments}

Whatever the origin of the free moments, they are clearly
correlated with the $H$-dependence of $\kappa$ for in-plane
fields. Although there has been no report in the cuprate context,
it has been known that phonons can be scattered by free
paramagnetic moments that form a two-level system,\cite{Berman}
where the energy levels of the two states split in magnetic fields
by $\Delta E$ due to the Zeeman energy. In this situation, phonons
are typically scattered in the following way: A phonon with energy
$\Delta E$ excites the lower-level state of a magnetic moment to
the higher-level state and is absorbed by a direct process; then,
another phonon with the same energy (but with a different wave
vector) is subsequently emitted by the higher-level state as it
relaxes. Since the phonon spectrum has a (broad) maximum at $\sim
4k_{B}T$, this ``paramagnetic scattering" of phonons is the
strongest when the Zeeman energy $\Delta E$ is equal to $4k_{B}T$.
Naturally, this causes a dip at $H_{\rm dip}$ in the $\kappa(H)$
isotherm with $H_{\rm dip}$ roughly proportional to
$T$,\cite{Berman} which is actually the case in our data below 2 K
[see Fig. 2(f)]. Note that the free spins themselves do not carry
heat because there is no dispersive collective excitations due to
the negligibly small interactions between these diluted spins.

In passing, we note that the peculiar $\kappa(H)$ behavior in NCO
(Ref. \onlinecite{Jin}) may also be understood in terms of the
paramagnetic scattering of phonons without employing the magnon
heat transport, at least for $T >$ 1.5 K: It has been discussed
that the doublet ground state of the Nd$^{3+}$ ions in NCO is
split slightly ($\sim$ 3 K) in zero field due to the interaction
with Cu spins,\cite{Sachidanandam} which causes the Nd$^{3+}$ ions
to form a two-level system; hence, the phonons may well be
scattered by paramagnetic Nd$^{3+}$ ions at $T >$ 1.5 K, since the
Nd ions show Ne\'el order only below $\sim$ 1.5
K,\cite{Sachidanandam} and a sufficiently large $H$ can remove the
paramagnetic scatterings and enhance $\kappa$
significantly.\cite{Berman} Thus, although the magnons are likely
to be carrying heat below $\sim$ 1.5 K,\cite{Li} their role at
higher temperature had better be scrutinized.

We should point out that we do not currently have a good
understanding for the enhancement of $\kappa$ with $H
\parallel c$ shown in Fig. 2(e), but it may be that for this field
direction some additional heat carriers are induced and add a
major contribution on top of $\kappa_p(H)$. One possibility is a
remaining contribution of Pr magnons, although as much as 35\% of
the Pr$^{3+}$ sites are diluted by nonmagnetic La$^{3+}$ ions in
PLCO. Remember that La-free PCO single crystals had been reported
to exhibit a large enhancement of $\kappa$ with $H$,\cite{Sales}
which might share a common origin to the present behavior shown in
Fig. 2(e).

Although the interaction between phonons and the free spins has
been known for a long time,\cite{Berman} such phenomenon and its
impact on the heat transport have never been revealed for cuprate
materials. The present results demonstrate that the phonon heat
transport can be $H$ dependent at low temperature in high-$T_c$
cuprates whenever excitations of free spins (that cause a Schottky
anomaly in the specific heat) are present. The effect of the
paramagnetic scattering is particularly pronounced in PLCO,
probably because in this parent compound the phononic heat
transport is very clean and also the Schottky anomaly is large. In
the superconducting doping regime of cuprates, while the Schottky
anomaly is almost ubiquitously observed,\cite{Moler, Revaz,
Wright, Chen, Brugger, Hien, Schottky} its contribution is much
smaller than that in PLCO and also the phonons are additionally
scattered by doped carriers, both of which cause the contribution
of the paramagnetic scattering to be reduced. Nevertheless, it is
likely that some fraction of the $H$ dependence of $\kappa$ is due
to the paramagnetic scattering of phonons even in superconducting
samples, and therefore an analysis of $\kappa(H)$ should be done
very carefully. Note that, while {\it quantitative analyses} of
$\kappa(H)$ would be difficult in light of the paramagnetic
scattering of phonons, {\it qualitative comparisons} of the
$\kappa(H)$ behavior for different dopings \cite{Sun1,Ando2} would
still be legitimate, because the Schottky anomaly is observed
across the whole doping range without any clear doping dependence.

Lastly, the present result suggests, together with the curious
magnetic shape-memory effect in
La$_{2-x}$Sr$_{x}$CuO$_4$,\cite{Lavrov2} that the spin-phonon
coupling is rather strong in cuprates. The significance of the
spin-phonon coupling in {\it doped} cuprates, although still
remains to be verified, is rather naturally expected based on the
behavior of the parent compound. Because of this coupling, AF
fluctuations and phonons may couple in an unusual way in these
materials and this could be the reason why the mechanism of
high-$T_c$ superconductivity is so difficult to understand. In
this regard, it is useful to note that possible roles of the
spin-phonon coupling in the peculiar electronic state of doped
cuprates have been theoretically discussed.\cite{Normand,
Jarlborg}

\section{Summary}

We measure the magnetic-field dependence of the thermal
conductivity of a parent insulating cuprate
Pr$_{1.3}$La$_{0.7}$CuO$_4$ to elucidate the possible field
dependence of phonon heat transport in cuprates. A striking dip
feature is found in $\kappa(H)$ isotherms for in-plane fields at
low temperature. The ``dip" field is found to increase linearly
with $T$ at low temperature, which is indicative of the phonon
scattering by free spins. Both the low-$T$ specific heat and the
magnetic susceptibility are measured to give supplemental evidence
for the existence of free spins. Since the Schottky-type specific
heat that signifies the existence of free spins is found
ubiquitously in high-$T_c$ cuprates in a wide doping range, the
present result suggests that the $H$-dependence of the phonon heat
transport should not be naively neglected, particularly when one
tries to extract the QP information from the quantitative analyses
of $\kappa(H)$. In addition, the present result points to a strong
spin-phonon coupling in high-$T_c$ cuprates.

\begin{acknowledgments}
We thank Y. Aoki, A. N. Lavrov, and J. Takeya for helpful
discussions. This work was supported by the Grant-in-Aid for
Science provided by the Japan Society for the Promotion of
Science.
\end{acknowledgments}


An Adaptive Nature-inspired Fog Architecture

Abstract: During the last decade, Cloud computing has efficiently exploited the economy of scale by providing low cost computational and storage resources over the Internet, eventually leading to consolidation of computing resources into large data centers. However, the nascent of the highly decentralized Internet of Things (IoT) technologies that cannot effectively utilize the centralized Cloud infrastructures pushes computing towards resource dispersion. Fog computing extends the Cloud paradigm by enabling dispersion of the computational and storage resources at the edge of the network in a close proximity to where the data is generated. In its essence, Fog computing facilitates the operation of the limited compute, storage and networking resources physically located close to the edge devices. However, the shared complexity of the Fog and the influence of the recent IoT trends moving towards deploying and interconnecting extremely large sets of pervasive devices and sensors, requires exploration of adaptive Fog architectural approaches capable of adapting and scaling in response to the unpredictable load patterns of the distributed IoT applications. In this paper we introduce a promising new nature-inspired Fog architecture, named SmartFog, capable of providing low decision making latency and adaptive resource management. By utilizing novel algorithms and techniques from the fields of multi-criteria decision making, graph theory and machine learning we model the Fog as a distributed intelligent processing system, therefore emulating the function of the human brain.

\section{Introduction}\label{sec:introduction}}
\else
\section{Introduction}
\label{sec:introduction}
\fi

\IEEEPARstart{T}{his} demo file is intended to serve as a ``starter file''
for IEEE Computer Society journal papers produced under \LaTeX\ using
IEEEtran.cls version 1.8b and later.
I wish you the best of success.

\hfill mds

\hfill August 26, 2015

\subsection{Subsection Heading Here}
Subsection text here.

\subsubsection{Subsubsection Heading Here}
Subsubsection text here.

\section{Conclusion}
The conclusion goes here.

\appendices
\section{Proof of the First Zonklar Equation}
Appendix one text goes here.

\section{}
Appendix two text goes here.

\ifCLASSOPTIONcompsoc

  \section*{Acknowledgments}
\else

  \section*{Acknowledgment}
\fi

The authors would like to thank...

\ifCLASSOPTIONcaptionsoff
  \newpage
\fi

\section{Introduction}
This demo file is intended to serve as a ``starter file''
for IEEE Computer Society conference papers produced under \LaTeX\ using
IEEEtran.cls version 1.8b and later.
I wish you the best of success.

\hfill mds

\hfill August 26, 2015

\subsection{Subsection Heading Here}
Subsection text here.

\subsubsection{Subsubsection Heading Here}
Subsubsection text here.

\section{Conclusion}
The conclusion goes here.

\ifCLASSOPTIONcompsoc

  \section*{Acknowledgments}
\else

  \section*{Acknowledgment}
\fi

The authors would like to thank...

\section{Introduction}
\IEEEPARstart{T}{his} demo file is intended to serve as a ``starter file''
for IEEE Communications Society journal papers produced under \LaTeX\ using
IEEEtran.cls version 1.8b and later.
I wish you the best of success.

\hfill mds

\hfill August 26, 2015

\subsection{Subsection Heading Here}
Subsection text here.

\subsubsection{Subsubsection Heading Here}
Subsubsection text here.

\section{Conclusion}
The conclusion goes here.

\appendices
\section{Proof of the First Zonklar Equation}
Appendix one text goes here.

\section{}
Appendix two text goes here.

\section*{Acknowledgment}

The authors would like to thank...

\ifCLASSOPTIONcaptionsoff
  \newpage
\fi

\section{Introduction}
\IEEEPARstart{T}{his} demo file is intended to serve as a ``starter file''
for IEEE journal papers produced under \LaTeX\ using
IEEEtran.cls version 1.8b and later.
I wish you the best of success.

\hfill mds

\hfill August 26, 2015

\subsection{Subsection Heading Here}
Subsection text here.

\subsubsection{Subsubsection Heading Here}
Subsubsection text here.

\section{Conclusion}
The conclusion goes here.

\appendices
\section{Proof of the First Zonklar Equation}
Appendix one text goes here.

\section{}
Appendix two text goes here.

\section*{Acknowledgment}

The authors would like to thank...

\ifCLASSOPTIONcaptionsoff
  \newpage
\fi

\section{Introduction}
\IEEEPARstart{T}{his} demo file is intended to serve as a ``starter file''
for IEEE \textsc{Transactions on Magnetics} journal papers produced under \LaTeX\ using
IEEEtran.cls version 1.8b and later.
I wish you the best of success.

\hfill mds

\hfill August 26, 2015

\subsection{Subsection Heading Here}
Subsection text here.

\subsubsection{Subsubsection Heading Here}
Subsubsection text here.

\section{Conclusion}
The conclusion goes here.

\appendices
\section{Proof of the First Zonklar Equation}
Appendix one text goes here.

\section{}
Appendix two text goes here.

\section*{Acknowledgment}

The authors would like to thank...

\ifCLASSOPTIONcaptionsoff
  \newpage
\fi

\section{Introduction}
This document is a model and instructions for \LaTeX.
Please observe the conference page limits.

\section{Ease of Use}

\subsection{Maintaining the Integrity of the Specifications}

The IEEEtran class file is used to format your paper and style the text. All margins,
column widths, line spaces, and text fonts are prescribed; please do not
alter them. You may note peculiarities. For example, the head margin
measures proportionately more than is customary. This measurement
and others are deliberate, using specifications that anticipate your paper
as one part of the entire proceedings, and not as an independent document.
Please do not revise any of the current designations.

\section{Prepare Your Paper Before Styling}
Before you begin to format your paper, first write and save the content as a
separate text file. Complete all content and organizational editing before
formatting. Please note sections \ref{AA}--\ref{SCM} below for more information on
proofreading, spelling and grammar.

Keep your text and graphic files separate until after the text has been
formatted and styled. Do not number text heads---{\LaTeX} will do that
for you.

\subsection{Abbreviations and Acronyms}\label{AA}
Define abbreviations and acronyms the first time they are used in the text,
even after they have been defined in the abstract. Abbreviations such as
IEEE, SI, MKS, CGS, ac, dc, and rms do not have to be defined. Do not use
abbreviations in the title or heads unless they are unavoidable.

\subsection{Units}
\begin{itemize}
\item Use either SI (MKS) or CGS as primary units. (SI units are encouraged.) English units may be used as secondary units (in parentheses). An exception would be the use of English units as identifiers in trade, such as ``3.5-inch disk drive''.
\item Avoid combining SI and CGS units, such as current in amperes and magnetic field in oersteds. This often leads to confusion because equations do not balance dimensionally. If you must use mixed units, clearly state the units for each quantity that you use in an equation.
\item Do not mix complete spellings and abbreviations of units: ``Wb/m\textsuperscript{2}'' or ``webers per square meter'', not ``webers/m\textsuperscript{2}''. Spell out units when they appear in text: ``. . . a few henries'', not ``. . . a few H''.
\item Use a zero before decimal points: ``0.25'', not ``.25''. Use ``cm\textsuperscript{3}'', not ``cc''.)
\end{itemize}

\subsection{Equations}
Number equations consecutively. To make your
equations more compact, you may use the solidus (~/~), the exp function, or
appropriate exponents. Italicize Roman symbols for quantities and variables,
but not Greek symbols. Use a long dash rather than a hyphen for a minus
sign. Punctuate equations with commas or periods when they are part of a
sentence, as in:
\begin{equation}
a+b=\gamma\label{eq}
\end{equation}

Be sure that the
symbols in your equation have been defined before or immediately following
the equation. Use ``\eqref{eq}'', not ``Eq.~\eqref{eq}'' or ``equation \eqref{eq}'', except at
the beginning of a sentence: ``Equation \eqref{eq} is . . .''

\subsection{\LaTeX-Specific Advice}

Please use ``soft'' (e.g., \verb|\eqref{Eq}|) cross references instead
of ``hard'' references (e.g., \verb|(1)|). That will make it possible
to combine sections, add equations, or change the order of figures or
citations without having to go through the file line by line.

Please don't use the \verb|{eqnarray}| equation environment. Use
\verb|{align}| or \verb|{IEEEeqnarray}| instead. The \verb|{eqnarray}|
environment leaves unsightly spaces around relation symbols.

Please note that the \verb|{subequations}| environment in {\LaTeX}
will increment the main equation counter even when there are no
equation numbers displayed. If you forget that, you might write an
article in which the equation numbers skip from (17) to (20), causing
the copy editors to wonder if you've discovered a new method of
counting.

{\BibTeX} does not work by magic. It doesn't get the bibliographic
data from thin air but from .bib files. If you use {\BibTeX} to produce a
bibliography you must send the .bib files.

{\LaTeX} can't read your mind. If you assign the same label to a
subsubsection and a table, you might find that Table I has been cross
referenced as Table IV-B3.

{\LaTeX} does not have precognitive abilities. If you put a
\verb|\label| command before the command that updates the counter it's
supposed to be using, the label will pick up the last counter to be
cross referenced instead. In particular, a \verb|\label| command
should not go before the caption of a figure or a table.

Do not use \verb|\nonumber| inside the \verb|{array}| environment. It
will not stop equation numbers inside \verb|{array}| (there won't be
any anyway) and it might stop a wanted equation number in the
surrounding equation.

\subsection{Some Common Mistakes}\label{SCM}
\begin{itemize}
\item The word ``data'' is plural, not singular.
\item The subscript for the permeability of vacuum $\mu_{0}$, and other common scientific constants, is zero with subscript formatting, not a lowercase letter ``o''.
\item In American English, commas, semicolons, periods, question and exclamation marks are located within quotation marks only when a complete thought or name is cited, such as a title or full quotation. When quotation marks are used, instead of a bold or italic typeface, to highlight a word or phrase, punctuation should appear outside of the quotation marks. A parenthetical phrase or statement at the end of a sentence is punctuated outside of the closing parenthesis (like this). (A parenthetical sentence is punctuated within the parentheses.)
\item A graph within a graph is an ``inset'', not an ``insert''. The word alternatively is preferred to the word ``alternately'' (unless you really mean something that alternates).
\item Do not use the word ``essentially'' to mean ``approximately'' or ``effectively''.
\item In your paper title, if the words ``that uses'' can accurately replace the word ``using'', capitalize the ``u''; if not, keep using lower-cased.
\item Be aware of the different meanings of the homophones ``affect'' and ``effect'', ``complement'' and ``compliment'', ``discreet'' and ``discrete'', ``principal'' and ``principle''.
\item Do not confuse ``imply'' and ``infer''.
\item The prefix ``non'' is not a word; it should be joined to the word it modifies, usually without a hyphen.
\item There is no period after the ``et'' in the Latin abbreviation ``et al.''.
\item The abbreviation ``i.e.'' means ``that is'', and the abbreviation ``e.g.'' means ``for example''.
\end{itemize}
An excellent style manual for science writers is \cite{b7}.

\subsection{Authors and Affiliations}
\textbf{The class file is designed for, but not limited to, six authors.} A
minimum of one author is required for all conference articles. Author names
should be listed starting from left to right and then moving down to the
next line. This is the author sequence that will be used in future citations
and by indexing services. Names should not be listed in columns nor group by
affiliation. Please keep your affiliations as succinct as possible (for
example, do not differentiate among departments of the same organization).

\subsection{Identify the Headings}
Headings, or heads, are organizational devices that guide the reader through
your paper. There are two types: component heads and text heads.

Component heads identify the different components of your paper and are not
topically subordinate to each other. Examples include Acknowledgments and
References and, for these, the correct style to use is ``Heading 5''. Use
``figure caption'' for your Figure captions, and ``table head'' for your
table title. Run-in heads, such as ``Abstract'', will require you to apply a
style (in this case, italic) in addition to the style provided by the drop
down menu to differentiate the head from the text.

Text heads organize the topics on a relational, hierarchical basis. For
example, the paper title is the primary text head because all subsequent
material relates and elaborates on this one topic. If there are two or more
sub-topics, the next level head (uppercase Roman numerals) should be used
and, conversely, if there are not at least two sub-topics, then no subheads
should be introduced.

\subsection{Figures and Tables}
\paragraph{Positioning Figures and Tables} Place figures and tables at the top and
bottom of columns. Avoid placing them in the middle of columns. Large
figures and tables may span across both columns. Figure captions should be
below the figures; table heads should appear above the tables. Insert
figures and tables after they are cited in the text. Use the abbreviation
``Fig.~\ref{fig}'', even at the beginning of a sentence.

\begin{table}[htbp]
\caption{Table Type Styles}
\begin{center}
\begin{tabular}{|c|c|c|c|}
\hline
\textbf{Table}&\multicolumn{3}{|c|}{\textbf{Table Column Head}} \\
\cline{2-4}
\textbf{Head} & \textbf{\textit{Table column subhead}}& \textbf{\textit{Subhead}}& \textbf{\textit{Subhead}} \\
\hline
copy& More table copy$^{\mathrm{a}}$& &  \\
\hline
\multicolumn{4}{l}{$^{\mathrm{a}}$Sample of a Table footnote.}
\end{tabular}
\label{tab1}
\end{center}
\end{table}

\begin{figure}[htbp]
\centerline{\includegraphics{fig1.png}}
\caption{Example of a figure caption.}
\label{fig}
\end{figure}

Figure Labels: Use 8 point Times New Roman for Figure labels. Use words
rather than symbols or abbreviations when writing Figure axis labels to
avoid confusing the reader. As an example, write the quantity
``Magnetization'', or ``Magnetization, M'', not just ``M''. If including
units in the label, present them within parentheses. Do not label axes only
with units. In the example, write ``Magnetization (A/m)'' or ``Magnetization
\{A[m(1)]\}'', not just ``A/m''. Do not label axes with a ratio of
quantities and units. For example, write ``Temperature (K)'', not
``Temperature/K''.

\section*{Acknowledgment}

The preferred spelling of the word ``acknowledgment'' in America is without
an ``e'' after the ``g''. Avoid the stilted expression ``one of us (R. B.
G.) thanks $\ldots$''. Instead, try ``R. B. G. thanks$\ldots$''. Put sponsor
acknowledgments in the unnumbered footnote on the first page.

\section*{References}

Please number citations consecutively within brackets \cite{b1}. The
sentence punctuation follows the bracket \cite{b2}. Refer simply to the reference
number, as in \cite{b3}---do not use ``Ref. \cite{b3}'' or ``reference \cite{b3}'' except at
the beginning of a sentence: ``Reference \cite{b3} was the first $\ldots$''

Number footnotes separately in superscripts. Place the actual footnote at
the bottom of the column in which it was cited. Do not put footnotes in the
abstract or reference list. Use letters for table footnotes.

Unless there are six authors or more give all authors' names; do not use
``et al.''. Papers that have not been published, even if they have been
submitted for publication, should be cited as ``unpublished'' \cite{b4}. Papers
that have been accepted for publication should be cited as ``in press'' \cite{b5}.
Capitalize only the first word in a paper title, except for proper nouns and
element symbols.

For papers published in translation journals, please give the English
citation first, followed by the original foreign-language citation \cite{b6}.

\section{Introduction}
\input{sections/Introduction.tex}
\section{Terminology}
\label{sec:terminology}
\input{sections/Terminology.tex}
\section{Related work}
\label{sec:related}
\input{sections/Related.tex}
\section{Background}
\label{sec:background}
\input{sections/Background.tex}
\section{Architectural model}
\label{sec:architecture}
\input{sections/Architecture.tex}
\section{Implementation}
\label{sec:implementation}
\input{sections/Implementation.tex}
\section{Experimental evaluation}
\label{sec:evaluation}
\input{sections/Evaluation.tex}
\section{Conclusion}
\label{sec:conclusion}
\input{sections/Conclusion.tex}

\section*{Acknowledgment}
This work is being accomplished as a part of the Tiroler Cloud Project: ”Energiebewusste F{\"o}derierte Cloud f{\"u}r Anwendungen aus Industrie und Forschung”, funded by the Bridge programme under grant agreement No 848448.

\bibliographystyle{unsrt}

\subsection{Communication gateways selection}
The process of Communication gateways selection is of utmost importance for reducing the decision making latency and increasing the resources utilization. Essentially, the Communication gateways act as a high-bandwidth aggregation points between the edge and the Cloud. Furthermore, they perform essential operations for clustering of the Fog devices in functional areas and mapping of the IoT applications to the physical resources.
The selection of the Communication gateways is performed by the Cloud layer, which sorts the available Fog devices based on dominance in multi-dimensional search space. The approach that we take to perform the non-domination sorting is based on the Fast Elitist Non-dominated Sorting Genetic Algorithm (NSGA-II) \cite{deb2000fast}. In this approach, every Fog device from the full set of available devices, is checked with a partially filled set $F'$. To begin with, the first Fog device from the list of devices is initially kept in the set $F'$. Afterwards, each Fog device $f$ is compared with all members of the set $F'$ one by one. If the Fog device $f$ dominates any other fog device $q$ in the set $F'$, then the device $q$ is removed from the set. On the other hand, if device $f$ is dominated by all member of $F'$, then $f$ is ignored. Moreover, if device $f$ is not dominated by any device in $F'$, then it is included in $F'$. When all Fog devices known to the cloud are checked, the remaining members of $F'$ constitute the non-dominated set.

The non-domination sorting algorithm is performed by considering the following conflicting criteria: (i) betweenness centrality, (ii) computational performance, and (iii) communication latency to the Cloud. The main purpose of the non-domination soring is to identify the Fog devices with highest number of available resources and the most optimal position in the Fog system. The betweenness centrality, as one of the objectives in the non-domination sorting algorithm, provides information on how well the each Fog device is connected in relation to the other Fog devices, the Cloud layer above and the IoT layer below. Respectively, the computational performance objective provides information on the computational capacity of the Fog devices. For the purposes of the SmartFog system, the computational performance is represented through the number of instructions that the processor of the Fog device can execute in one second - MIPS. The communication latency objective is represented through the virtual link latency between a given Fog device and any other device. In our implementation, historical data on the previous data transfers is being preserved to calculate the average latency between any two devices in the Fog environment.

The process of non-domination sorting algorithm results with a set of multiple optimal solutions, visualized in the form of Pareto front. As all of the Fog devices represented in the Pareto front are optimal, we utilize automated decision making strategy, which considers a-priory input information to select only the required number of Fog devices to act as Communication gateways. The decision making algorithm considers the number $N_f$ and the type $F$ of the functional areas to give priority which Fog devices to be selected. For example, if a given functional area should be optimized for processing, then the decision making algorithm will select the Fog device with the highest computational performance from the set $F'$ of Pareto optimal Fog devices. The selection process of the Communication gateways is presented in Algorithm \ref{alg:sort}.
\begin{algorithm}
\label{alg1}
\caption{Communication gateways selection}\label{alg:sort}
\begin{algorithmic}[1]
\INPUT
\Statex $N$ \Comment{Number of Fog devices}
\Statex $L = (L_1, L_2, ... , L_N)$ \Comment{List of Fog devices}
\Statex $C = (C_1, C_2, ... , C_N)$ \Comment{CPUs per Fog device}
\Statex $M = (M_1, M_2, ... , M_N)$ \Comment{Memory per Fog device}
\Statex $T$ \Comment{Topology of the Fog layer}
\Statex $N_f$ \Comment{Number of functional areas}
\Statex $F = (F_1, F_2, ..., F_{N_f} )$ \Comment{Type of the functional area}
\While{$i < N$}
\State $B_i \gets evaluate\_betweenness\_centrality(L_i, T)$
\State $E_i \gets evaluate\_Fog\_device(C_i, M_i,B_i)$
\State $i \gets i+1$
\EndWhile\label{mainloop}
\State $S \gets non\_domination\_sorting(E)$
\State $D \gets automated\_decision\_making (S, N_f,F)$
\end{algorithmic}
\end{algorithm}

\subsection{Functional areas clustering}
The concept of functional areas, introduced by the SmartFog architecture, evolves around the notion of Fog devices' grouping into distinctive clusters optimized for special types of IoT applications, such as compute or memory intensive. Essentially, the clustering process enables grouping of logically similar Fog devices by considering various resource related criteria.

The clustering process and the formation of the functional areas is performed in distributed manner among the communication gateways. Every communication gateway performs Spectral clustering by considering the available computing and memory resources of the available Fog devices. A-priory input from the Cloud layer is used to steer the clustering process in the preferred direction, therefore creating clusters of Fog devices with similar amount of available processing and memory resources. For example, if a communication gateway is intended to create compute optimized functional area, the Spectral clustering will identify all similar Fog devices, which have sufficient amount of processing resources. The process of functional areas formation is presented in Algorithm \ref{alg:cluster}. It is essential to be noted, that every communication gateway creates each own cluster, therefore a given Fog device can be member of multiple functional areas, provided it has sufficient resources.

\begin{algorithm}
\caption{Functional areas clustering}\label{alg:cluster}
\begin{algorithmic}[1]
\INPUT
\Statex $N$ \Comment{Number of Fog devices}
\Statex $k$ \Comment{Number of clusters}
\Statex $C = (C_1, C_2, ... , C_N)$ \Comment{CPUs per Fog device}
\Statex $M = (M_1, M_2, ... , M_N)$ \Comment{Memory per Fog device}
\Statex $G$ \Comment{Value of the Gaussian filter}
\State $L \gets create\_dense\_matrix(C, M)$
\State $Y \gets find\_eigen\_vectors(L, k)$
\State $K_m \gets apply\_k\_means(Y, k)$
\State $F_k \gets create\_functional_clusters (K_m)$

\end{algorithmic}
\end{algorithm}

\subsection{Multi-objective optimization}
In this work we extend and utilize essential concepts from the area of multi-criteria optimization and decision making with a main goal to enable efficient communication between the IoT devices and the Fog/Cloud systems above. In the most general sense, optimization is a process of identifying one or multiple solutions, which correspond to the extreme values of two or more objective functions within given constraints set. In the cases in which the optimization task utilizes only a single objective function it results in a single optimal solution. Moreover, the optimization can also consider multiple conflicting objectives simultaneously. In those circumstances, the process will usually result in a set of optimal trade-off solutions, so-called Pareto front. The task of finding the optimal set of Pareto solutions in the form of Pareto front is known in the literature as a multi-objective optimization~\cite{branke2008multiobjective}~\cite{kimovski2016multi}. The Pareto front is an essential tool for decision support and preference discovery, whose shape provides new insights and allows scientists to explore the space of non-dominated solutions, possibly revealing regions of interest that are impossible to see otherwise.

In general, the multi-objective optimization problem involves two or more objective functions which have to be either minimized or maximized. The problem of optimization can be formulated as:
\begin{equation}
min/max (f_1(Y), f_2(Y),\dots,f_n(Y))
\end{equation}
where $n \geq 2$ is the number of objectives functions $f$ that we want to minimize or maximize, while $Y=(y_1,y_2,\dots,y_k)$ is a region enclosing the set of feasible decision vectors.

\subsection{Betweenness centrality}
In the research field of graph theory, the  \textit{betweenness centrality} is used as a measure of node centrality within a graph, in relation to the number of shortest paths passing through the node (vertex) \cite{barthelemy2004betweenness}. For each pair of vertices in a connected unweighted graph, there exists at least one shortest path, such that the number of edges that the route passes through can be minimized. Respectively, for weighted graphs the betweenness centrality utilizes the the sum of the edges' weights to find the shortest paths. The betweenness centrality for each vertex is the number of the shortest paths that pass through the vertex and can be represented through the following equation:
\begin{equation}
g(n)= \sum_{s \neq n \neq d} \frac{\sigma_{sd}(n)}{\sigma_{sd}}
\end{equation}
where $\sigma_{sd}$ is the total number of shortest path routes from vertex $s$ to vertex $d$, while $\sigma_{sd}(n)$ represents the number of shortest paths from vertex $s$ to vertex $d$ passing through vertex $n$.

High value of the betweenness centrality usually implies that a given vertex can reach others with lowest possible latency. Moreover, high betweenness centrality can also indicate that a given node lies on the path of many shortest routes. Consequently, if one removes a node with large centrality it can lengthen the paths between many other pairs of nodes. This concept can be utilized in Fog computing in order to find the most central Fog devices, therefore enabling more efficient communication between the IoT devices and the Fog.

\subsection{Spectral clustering}
In the Big Data era the vital tool for dealing with large data-sets is the concept of classification or grouping of data objects into a set of categories or clusters. The classification of the objects is conducted based on the similarity or dissimilarity of multiple features that describe them. Essentially, the classification methods can be divided into two categories, namely supervised and unsupervised \cite{cluster}. In supervised classification, the features' mapping from a set of input data vectors is classified to a finite set of discrete labeled classes and it is modeled in terms of some mathematical function. On the other hand, in unsupervised classification, called clustering, no labeled data-sets are available. The aim of the clustering is to separate a finite unlabeled data-sets into a finite and discrete set of clusters.
For the purpose of efficient clustering of the Fog devices into functional groups we utilize unsupervised clustering technique called Spectral clustering \cite{spectral2007tutorial}. Spectral clustering is a technique that make use of the eigenvalues \cite{edelman1988eigenvalues} of the similarity matrix $S_{ij}=s(x_i,x_j)$ of the data points $x_i,x_j$ to perform dimensionality reduction before clustering the data in relation to a reduced number of dimensions. The commonly used similarity measures in Spectral clustering are based on the Euclidean distance and the Gaussian kernel. After performing dimensionality reduction, Spectral clustering relies on more simple clustering techniques, such as k-means \cite{kmeans}, to group the similar data points in distinctive similarity based sub-graphs. Spectral clustering is suitable for finding resemblance between different data points by utilizing the concept of graph similarity, making it easily applicable for sub-diving of network graphs, such as the Fog environment.

\subsection{Communication gateway non-domination sorting and decision making}
\label{cg}
The essential feature of the SmartFog architecture is the introduction of the communication gateways, which act both as a aggregation points between the Edge and the Cloud, and are actively involved in the creation of the functional areas. For those reasons, it is important to evaluate the behavior, efficiency and scalability of the implemented non-domination sorting and decision making algorithms. For the evaluation purposes we utilize network overlays with sizes varying from 20 to 40 vertices, i.e. Fog devices. In the simulation test-bed we assumed that the Fog devices have limited computation resources, ranging from 800 to 1200 MIPS. In relation to the memory and storage resources, we assume that each device can have between 1 and 4 GB of operating memory and small amount of solid-state based storage. All experiments were repeated 100 times and the median value, together with the standard deviation, are presented in the Figures below.

Initially, the evaluation activities were focused towards validation of the concepts, introduced from the areas of graph theory and multi-criteria decision making. Therefore, we first explore the connectivity characteristics, expressed as betweenness centrality, of every Fog device, which are utilized in the multi-criteria non-domination sorting. The exploration of the scalability potential of the betweenness centrality algorithm is essential, as it should be executed each time there are changes in the highly volatile IoT/Fog network overlay. Moreover, this process requires calculation of all shortest path routes in a given topology, therefore increasing the complexity of the algorithm. The evaluation of the the betweenness centrality algorithm is presented on Figure \ref{fig:Between}.

\begin{figure}[h]
    \centering
    \includegraphics[width=0.4\textwidth]{images/Fig_3.png}
    \caption{Computation scalability of the betweenness centrality algorithm}
    \label{fig:Between}
\end{figure}

As it can be observed from the simulation results, the execution time of the betweenness centrality calculation rises linearly with the size of the network overlay. Nevertheless, even for 40 Fog devices, which can potentially serve thousands of Edge and IoT devices, the execution time is below 40 $ms$.

Furthermore, on Figure \ref{fig:Nondom}, we provide evaluation data of the non-domination sorting algorithm. The experimental scenario is identical with the one described above. The execution time presented in the Figure includes the time required for the decision making to be performed, while omitting the betweenness centrality calculation. We can observe that for our testbed, the non-domination sorting was very efficient and induced latencies below 5 $ms$.

\begin{figure}[h]
    \centering
    \includegraphics[width=0.4\textwidth]{images/Fig_2.png}
    \caption{Computation scalability of the Non-domination sorting algorithm}
    \label{fig:Nondom}
\end{figure}

\subsection{Functional areas clustering}
The other concept, exploited by the SmartFog architecture, is the Spectral clustering algorithm, which is being utilized by the communication gateways to create specialized functional areas. The current implementation of the SmartFog architecture supports the creation of computational and memory optimized functional areas. For the purpose of Spectral clustering we again consider varying sizes of the network, which ranges from 20 to 40 Fog devices. The evaluation results, presented on Figure \ref{fig:Spectral}, clearly show that the algorithm scales very well, with latencies in the range of 300 $ms$ and very low standard deviation values. Moreover, on Figure \ref{fig:clus} sample clustering result for three functional areas is presented.

\begin{figure}[h]
    \centering
    \includegraphics[width=0.4\textwidth]{images/Fig_1.png}
    \caption{Computation scalability of the Spectral clustering algorithm}
    \label{fig:Spectral}
\end{figure}

\begin{figure}[h]
    \centering
    \includegraphics[width=0.4\textwidth]{images/Fig_4.png}
    \caption{Sample results from the functional areas clustering}
    \label{fig:clus}
\end{figure}

\subsection{Simulation experiment}
In order to explore the influence of the introduced concepts on the overall efficiency of the Fog architectures we have conducted extensive simulation-based analysis. We have simulated a Sense-Process-Actuate Model (SPAM) with Edge-ward placement strategy. In SPAM, the sensors continuously emit data in the form of tuples. Tuple is basic unit of communication between different entities in iFogSim, which is specified by the following parameters: data source, data destination and amount of required computation and network resources. These tuples are transmitted to the Fog devices as data streams. The IoT applications modules, executed on the Fog devices, process this data and transmit it back to the actuators for performing a specific action. For the Edge-ward placement strategy, specific application's modules are deployed at edge of the network, while the remaining modules are placed to act as communication gateways to the Cloud. Therefore part of the data is being processed by Fog devices placed at edge of network and the remaining data is transmitted to the Cloud for further processing.

In our simulation model, the sensors and actuators are attached to the edge devices (mobiles). These mobiles are connected with the Fog devices through edge-gateways. Whereas on the other end, the Fog devices are connected with the Cloud through the communication gateways. We have used edge-ward placement strategy in our simulation. Therefore we have defined two process-control loops to calculate the end-to-end latency or IoT application loop delay. The first is a sense-process-actuate (SPA) loop in which the Fog devices receive the data from the sensors, and then process it and transmit it back to the actuators for further action. The second one is the process-control (PC) loop in which, through communication gateways, the Cloud receives data from the Fog devices, and then process it and transmits it back to the source Fog devices.  Therefore, we conducted a series of experiments to evaluate the performance of the SmartFog against unoptimized-Fog in relation to the SPA and PC loop delays. These series consisted of three experimental scenarios with three overlay topologies with sizes of 20, 30, and 40 Fog devices respectively. The Fog devices in the overlays are connected randomly with one another creating a mesh kind of topology. Furthermore, the Fog devices are connected with the sensors and actuators, which are continuously emitting tuples. The Fog devices are also attached to Cloud through the communication gateways. For calculating the SPA loop delay, we have used tuples with a requirement of 1000-8000 MIPS and 100 bytes network width. Moreover, for the PC loop delay, the tuples have minimal requirement of 40000 MIPS and 100 bytes. In our experiments the tuples are sent randomly from the sensors to the Fog devices, while SmartFog manages the data transfers from the Fog to the Cloud. The simulation results from the loop delay evaluation are presented on Figures \ref{SPA} and \ref{PC}.

\begin{figure}[h]
    \centering
    \includegraphics[width=0.4\textwidth]{images/Fig_7.png}
    \caption{SPA loop latency in Fog environment}
    \label{SPA}
\end{figure}

The evaluation results, both for the SPA and the PC loop latency, clearly show that SmartFog can reduce the decision making latency by up to 8\%, which is significant improvement in volatile environments such as the Fog. In the cases when the data needs to be processed in the Cloud, more concretely within the PC loop, we can conclude that for smaller overlays SmartFog can reduce the communication latency, while for larges overlays the latency remains the same. Furthermore, in the case of the SPA loop, SmartFog can significantly reduce the latency, both for small and large overlays.

\begin{figure}[h]
    \centering
    \includegraphics[width=0.4\textwidth]{images/Fig_6.png}
    \caption{PC loop latency in Fog environment}
    \label{PC}
\end{figure}

Moreover, we have further focused our experimental evaluation on the network load induced by the IoT applications. For this reason, we have evaluated the total network traffic over the overlays for the SmartFog architecture and unoptimized Fog. The results, presented on Figure \ref{Nload}, clearly show that for all network sizes the communication load can be reduced by up to 13\%.

\begin{figure}[h]
    \centering
    \includegraphics[width=0.4\textwidth]{images/Fig_5.png}
    \caption{Network load in Fog environment}
    \label{Nload}
\end{figure}

\subsection{Simulation environment} \label{sec:over}
We simulated an unstructured Fog overlay network system on top of the iFogSim \cite{gupta2017ifogsim}, which is an efficient toolkit for modeling and simulating resource management techniques in IoT and Fog computing environments. To accommodate the novel concepts introduced by SmartFog we have extended iFogSim to support the automated gateways selection and devices' clustering.

Within the simulation environment we assume that every Fog node in the overlay network have a specific processing capacity and system architecture, such as ARM or x86. Furthermore, every node has limited operating memory and storage disk capacity with a fixed number of stored data items of varying size. In addition a specific scheme has been defined to allow every Fog node to leave and join the network after a random time interval, which is an important scenario supported by the SmartFog architecture.

\subsection{Multi-criteria decision making and spectral clustering}
\label{sec:optimization}
We perform the gateway selection based on multiple conflicting objectives, which have been modeled as a part of the information collected from the initial Fog overlay simulation. We utilize modified  NSGA-II multi-objective optimization algorithm to perform the non-domination sorting of the Fog devices. In addition to this, we also utilize low-latency automated decision making algorithm. To instantiate NSGA-II as a main component of the communication gateway selection module, we have extended the jMetal \cite{durillo2011jmetal} object-oriented Java framework for multi-objective optimization problems. We have implemented particular modifications in the jMetal framework to deal with the specific characteristics of the proposed architecture. More concretely, we have modified the non-domination sorting algorithm, as the standard operators of the jMetal framework do not support independent multi-criteria sorting. Furthermore, jMetal was extended to support automated decision-making by utilizing the algorithm proposed in \cite{dm}. This algorithm performs low-latency decision making by dividing the Pareto front in multiple regions based on the a-priory knowledge that gives priority to specific criteria. Therefore, this algorithm is suitable for the SmartFog architecture. Moreover, we modified the jMetal framework and extended its API to support integration with the iFogSim simulation environment.

Lastly, the Spectral clustering was implemented by extending the WEKA data mining framework \cite{hall2009weka}. The implementation of the spectral clustering was based on the algorithm proposed in \cite{ng2002spectral}. Additionally, the Spectral clustering implementation was integrated within the jMetal framework, thus fully supporting the novel concepts introduced in the SmartFog.

Distributed, Parallel, and Cluster Computing (cs.DC)