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

~~~The most unique property of the Galactic microquasar SS 433 is
continuous launching of the highly collimated (with the opening angle
$ \Theta_j\approx1.5$$^\circ~$), mildly relativistic (with the bulk velocity
$ V=\beta c\approx0.26c $), baryonic (i.e. composed of atoms of hydrogen,
helium and heavy elements with the relative proportion close to the
solar one) jets, which reveal themselves via Doppler-shifted optical
and X-ray emission lines detected in the spectrum of the unresolved
core of the system (see \cite{Fabrika2004} for a review). The Doppler
shifts of the lines periodically vary with time in accord with the
kinematic model \citep{Abell1979,Fabian1979,Milgrom1979} postulating
that the jets undergo precessional motion about an axis inclined at $
i=78.05\pm0.05\approx$78$^\circ~$ relative to the line of sight with
half-amplitude $\Theta_p=20.92\pm0.08\approx $21$^\circ~$ and period
$P=162.375\pm0.011\approx162$ days (\citealt{Eikenberry2001}).

~~Apart from the unresolved core, the presence of the jets is clearly
traced by radio emission at distances from $ \sim 10$ milliarcseconds
(e.g. \cite{Marshall2013}) to a few arcseconds, where it forms a
remarkable 'corkscrew' pattern generally consistent with the
prediction of the kinematic model (e.g. \citep{Miller2008}), albeit
with slight deflections that might be caused by the interaction of
the jets with the surrounding medium \citep{Panferov2014}. Further
away, however, the jets stay invisible up to $ \sim 10 $ arcminute
scales, where extended X-ray emission, which probably has relation to
their propagation through the interstellar medium (ISM), is observed
\citep{Brinkmann2005}, and then they appear to be terminated at the
boundary of the surrounding radio nebula W50 \citep{Dubner1998},
giving rise to its severe deformation \citep{Goodall2011} and bright
X-ray emission from the supposed termination region
\citep{Brinkmann1996}.

~~~Thanks to the excellent angular resolution of the \textit{Chandra}
X-ray observatory ($\approx 0.5$ arcsec), \cite{Marshall2002} have
discovered extended X-ray emission at a distance of the order of one
arcsec from the central source, coming from regions coincident with
the jets' propagation trajectory as traced by the radio emission
observed at the same angular scales. Moreover, in further
observations, oppositely Doppler-shifted emission lines were detected
in the spectra of the (on average) approaching (hereafter the 'blue'
jet) and receding (hereafter the 'red' jet) jets, with the lines'
centroids at $\approx 7.3 $ keV and $ \approx 6.4 $ keV, respectively
\citep{Migliari2002}, which indicates that the lines originate in the
material of the jets, still moving at a mildly relativistic speed.

~~~However, the exact identification of these lines remains elusive:
using simulations of the line emission of the extended jets averaged
over a precession period, \cite{Migliari2002} argued that the
rest-frame energy of the lines could be 7.06 keV. In such a case, however,
the line should be identified with the Fe I $  K\beta $ line of
neutral iron rather than with the Fe XXV $K\beta $ line (7.89 keV) of
highly ionized iron, as suggested by \cite{Migliari2002}. While the
presence of the Fe XXV line is a clear indication of hot ($ \gtrsim$
few keV) plasma, the Fe I line arises as a result of fluorescence
following inner shell ionization of neutral iron atoms by X-ray
photons or by interaction with energetic particles. The ambiguity in
the line identification leaves open the question about the mechanism
of production of the extended X-ray emission at arcsec scales,
corresponding to distances of order $\sim 10^{17} $ cm from the
central source.

~~~In this paper, we address the identification issue of the extended
emission lines taking into account time delays between signals coming from
different parts of the precessing jets due to finite light propagation
speed, which prove to be important in predicting the line positions
for radiation emitted at distances of order $ \sim\lambda=\beta
cP=1.1\times 10^{17}$ cm from the central source, in contrast to the
well-known X-ray and optical lines emitted at $\sim10^{11} $ cm and
$\sim10^{15} $ cm, respectively (see \cite{Fabrika2004}). Indeed, at
distances $ \sim\lambda $, the time delay between different regions of
the jets can reach $ \sim2\lambda\tan\Theta_p/c\sim 32$ days, which
corresponds to a significant precession phase shift of $\approx
0.2$ for the observed emission. This effect is also important and has
previously been taken into account for modelling images and
variability of the radio emission observed at the same angular scales
(see e.g. \cite{Bell2011} and references therein).

~~~As a simplest scenario, we may assume that the arcsec X-ray line
emission is produced at some specific distance $ r $ from the central
source or, equivalently, at some specific instant after launch of the
emitting blob from the jet's core. As $r$ gets larger, the observed
line positions will increasingly deviate from the ones predicted by
the kinematic model without accounting for light-travel-time
effects, eventually leading to ambiguity in the predicted line positions
due to simultaneous arrival of emission from different parts of the jet.
This will further lead to broadening and brightening of the observed
emission line if the emission is unresolved spatially and spectrally.

~~~In a more complicated scenario, the X-ray lines are radiated over a
range of distances (or delays after launch) from the central source,
and the observed line positions are determined by properly averaging
over the whole emitting region. In this case, they will still manifest
some modulation with the precession phase even for a brightening
duration comparable to the period of precession, as a result of
Doppler boosting of the emission from the approaching parts of the
jets relative to the receding ones.

~~~We demonstrate that the observed line positions are consistent
either with the Fe XXV K$\alpha$ ($E_0=6.7$  keV) line of helium-like
iron in the former scenario, or with the Fe XXVI Ly$\alpha$
($E_0=6.96$  keV) line of hydrogen-like iron in the latter, hence with
a hot gas origin of the lines in both cases. However, the required
precession phase of the central source in the case of identification
with Fe XXV K$\alpha$ disagrees with the observed one, whereas there
is good agreement in the case of the Fe XXVI Ly$\alpha$ line.

~~~Because quite high ($ \gtrsim 10$ keV) temperature is needed for
the Fe XXVI Ly$\alpha$ line to dominate over the Fe XXV K$\alpha$ line
in the thermal spectrum of hot plasma, the inferred identification
implies that the jets' material experiences efficient heating at
distances $ \sim 1.5\lambda $ from the central source lasting a
significant fraction ($ \sim 0.6$) of the precession period. This may
be caused by partial dissipation of the kinetic energy of the jets via
the propagation of a shock in their matter as a result of their
interaction with the surrounding medium (see,
e.g. \cite{Heavens1990}), or by mutual collisions of discrete dense
blobs within the jets due to more effecient deceleration of the
leading blobs by the interstellar medium with respect to the trailing
ones.

~~~The paper is organised as follows: in the first section, we
describe a model for the variation of the line positions as a function of
the precession phase depending on location and size of their emission
region along the jet. In the second section, we compare predictions of
the model with the observed positions of the emission lines for the
approaching and receding jets and draw a conclusion regarding
their possible identification. The last section is devoted to the
discussion of the results obtained and final conclusions of the
present study.

\section{Predicted line positions}
\label{s:model}
\begin{figure}[t]
\centering
\includegraphics[width=1.0\columnwidth,bb=50 280 550 680]{./figs/sketchx.pdf}
\caption{\small Geometry of the problem for the approaching (on
  average) jet. The precession axis coincides with the $z$ axis, and
  together with the direction to the observer $\vec{n}_{o}$, inclined
  at angle $ i\approx78^\circ~ $ to the former, they form the $xz$
  plane. The direction of the jet $\vec{n}_{j}$ makes an angle $
  \Theta_p\approx21^\circ~ $ with the precession axis, and an angle
  ranging from $ \alpha_{min}=i-\Theta_p\approx58^\circ~$ at precession
  phase $\phi=0$ (shown in blue) and $
  \alpha_{max}=i+\Theta_p\approx99^\circ~$ at $\phi=0.5$ (shown in red)
  with the line of sight. The black dashed line marks the projection
  of the picture plane. The differences in the distance to the
  observer from the approaching and receding parts of the jet relative to
  the distance from the central source are denoted by $c\Delta t_{0}$
  and $c\Delta t_{0.5}$ for the precession phases $\phi=0$ and
  $\phi=0.5$, respectively, and the regions of the 'blue' jet launched
  at $\phi=0.5$ are in fact receding from the observer.}
\label{f:sketch}
\end{figure}

Let us first consider the innermost regions of the precessing
jets, whose observational appearance is not affected by
light-travel-time effects. It is sufficient to consider just one of
the jets, e.g. the blue one, since the results for the other one can
be readily obtained by inversion of the jet's direction vector.

Let the $z$ axis be aligned with the precession axis and the $x$ axis
lie in the plane formed by the $z$ axis and the direction from the
source to the observer. Then, the $y$ axis and the line in the $xz$
plane normal to the line of sight form the picture plane (see
Fig. \ref{f:sketch}). Hence, the line-of-sight and precession axis
directions are given by vectors
\begin{equation}
\vec{n}_{o}=(\sin i,~0,~\cos i),~ \vec{n}_{p}=(0,~0,~1),
\end{equation}
where $i$ is the inclination angle of the precession axis relative to the line of sight.

The direction of the blue jet varies periodically in time and depends
on the precession phase as
\begin{equation}
 \vec{n}_{j}(\phi)=(\sin\Theta_p\cos 2\pi\phi,~\sin\Theta_p\sin 2\pi\phi,~\cos\Theta_p),
\end{equation}
where $ \Theta_p $ is the amplitude of precession, i.e. the
half-opening angle of the precession cone, $ P $ is the precession
period, and the zero precession phase is defined here as the moment
corresponding to the minimum angle $ \alpha_{min}=i-\Theta_p $ between
the jet's direction and the line of sight\footnote{Hereafter, the
  angle between the jet's direction and the line of sight is measured
  in the rest frame of the observer, so it does not need to be
  corrected for relativistic aberration, as would be the case for the
  angle measured in the reference frame comoving with the jet's
  matter. Time is also measured in the observer's reference frame
  and is thus not affected by relativistic dilation.}. The maximum
angle $ \alpha_{max}=i+\Theta_p $ then corresponds to $\phi=0.5$, and
since $\Theta_p>\frac{\pi}{2}-i$, the blue jet is actually receding at this
moment (see Fig. \ref{f:sketch}).
\begin{equation}
\mu(\phi)=\cos i~\cos\Theta_p+\sin i~\sin\Theta_p\cos 2\pi\phi,
\end{equation}
reaching the maximum $\mu_{max}=\cos(i-\Theta_p)$ at $ \phi=0$ and minimum $ \mu_{min}=\cos(i+\Theta_p)$ at $ \phi=0.5$.

The Doppler factor for the jet's matter moving ballistically with the bulk velocity $ V=\beta c$ is given by
\begin{equation}
\delta(\mu)=\frac{1}{\gamma(1-\mu(\phi)\beta)},
\end{equation}
with $ \gamma=1/\sqrt{1-\beta^2}$. Photons with energy $ E_0 $ emitted by the jet will be detected by the observer with the energy
\begin{equation}
E(\phi)=\frac{E_0}{\gamma(1-\mu(\phi)\beta)},
\label{eq:e}
\end{equation}
and this equation indeed provides an excellent fit to the observed line positions for the compact (unresolved) jets in SS 433 (with the line positions for the red jet obtained by substituting $ \mu $ by $ -\mu $).

Suppose now that the jet's matter becomes bright (e.g. after having been reheated) at some instant $\tau P$ (expressed in units of the precession period $P$ through the dimensionless time-delay parameter $\tau$) after its launch from the central object. The corresponding distance covered by the jet by this time is $ \beta c~\tau P=\lambda\tau $, where $ \lambda=\beta c P$ is the distance covered by the jets over one precession period. Photons from the reheating region of the jets will reach the observer with a delay
\begin{equation}
\Delta t(\tau,\mu)=(1-\mu\beta)\tau P
\end{equation}
relative to photons emitted by the central source, where the $-\mu\beta\tau P$ component is caused by the finite light travel time.

Thus, if the observer cannot spatially separate the emission of the
central source from the radiation produced at distance $ \lambda\tau $
from it, the detected signal will contain both photons coming from the
central source at its current precession phase $ \phi_0 $ and from
blobs launched from the central source at an earlier phase $
\phi_e=\phi_0-\Delta\phi$, where
\begin{equation}
\Delta\phi(\tau,\mu_e)=(1-\mu_e\beta)\tau,
\end{equation}
and $\mu_e(\tau,\phi_0)$ (cosine of the angle between the line of
sight and the velocity direction of the matter emitting at
distance $ \lambda\tau $) is determined implicitly by the non-linear
equation
\begin{equation}
\mu_e=\mu(\phi_0-\Delta\phi(\tau,\mu_e)).
\label{eq:mu}
\end{equation}

\begin{figure}
\centering
\includegraphics[width=0.4\textwidth,bb=80 200 530 680]{./figs/eTau05fnl.pdf}
\includegraphics[width=0.4\textwidth,bb=80 200 530 650]{./figs/eTau1fnl.pdf}
\includegraphics[width=0.4\textwidth,bb=80 180 530 650]{./figs/eTau15fnl.pdf}
\caption{\small Cosine of the angle between the observed jet's direction and the line of sight as a function of the current precession phase of the central source $ \phi_0 $ shifted by $ \tau $ for $ \tau=0.5,1.0,1.5 $ (from top to bottom). Model predictions for the blue and red jets are shown with blue and red solid lines, respectively. Variation of the line position for the compact jets is shown by the dotted line, while the dashed line shows the 'linear' approximation of the exact solution for the extended jets. The vertical dashed lines mark the stationary points corresponding to the jet's propagation in the picture plane.}
\label{f:mu}
\end{figure}

We see that the total phase shift $ \Delta\phi(\tau,\mu_e) $ is the sum of two components. It is the $\mu_e\tau\beta$ term that makes Eq. \ref{eq:mu} implicit, while the $\tau$ component just shifts the observed precession phase $ \phi_0 $ backward. When $ \beta\tau $ is small, one can linearise Eq. \ref{eq:mu} assuming that the non-linear phase shift is small, i.e. by substituting $ \mu_e $ in the argument of $ \Delta \phi(\tau,\mu_e) $ by $ \mu_0=\mu(\phi_0-\Delta\phi(\tau,\mu_0))$, so that
\begin{equation}
\mu_{e,lin}=\mu(\phi_0-\tau(1-\mu(\phi_0-\tau)\beta))).
\label{eq:mulin}
\end{equation}
This approximation works quite well up to $\beta\tau\lesssim0.1$ (see Fig. \ref{f:mu}).

The precession phases corresponding to blobs propagating in the picture plane (i.e. $ \mu=0 $), given by
\begin{equation}
\phi_{\perp,\pm}=\pm\arccos\left(\frac{-\cos\Theta_p\cos i}{\sin\Theta_p\sin i}\right)/2\pi\approx\pm0.34,
\end{equation}
are obviously not affected by light-travel-time effects at all and can
be considered stationary points for the considered
transformation. Interestingly, if the blobs get re-brightened at a
sufficiently long distance from the central source, it is possible
that photons emitted by an approaching ($ \mu>0 $) blob will reach the
observer before photons produced by a blob moving in the picture plane
($\mu=0$) despite the former having been launched later than the
latter. This can cause ambiguity in the expected position of the
observed X-ray line for time delays larger than a critical value that
can be found by minimising the expression
\begin{equation}
\tau=\frac{1}{\beta\mu}\left[\arccos\left(\frac{\mu-\cos\Theta_p\cos i}{\sin\Theta_p\sin i}\right)/2\pi-\phi_{\perp,+}\right]
\end{equation}
over $ \mu$ in the range from $ 0 $ to $\mu_{max}=\cos(i-\Theta_p)$. For the parameters of the SS 433 jets, we find $ \tau_1\approx1.8$ at $ \mu\approx0.28 $.

The distance for which photons emitted by a blob launched at the phase
$ \phi=0 $ (and having the maximum $ \mu $) begin to catch up with
photons from a blob launched earlier at the phase $\phi_{\perp,-} $ is given by
\begin{equation}
\tau_2=\frac{1}{\beta\mu_{max}}\left[\arccos\left(\frac{\mu_{max}-\cos\Theta_p\cos i}{\sin\Theta_p\sin i}\right)/2\pi-\phi_{\perp,+}\right],
\end{equation}
which yields $ \tau_2\approx2.35$ for SS 433 parameters, and this
marks the instant when half of all approaching blobs are involved in
the ambiguity of the predicted line position.

Further, $ \tau $ for which the arrival time of photons from the
$ \mu=\mu_{max} $ jet's regions coincides with the arrival time of
photons emitted by a blob launched at the phase $\phi_{\perp,+}$ of
the preceding precession cycle equals
\begin{equation}
\tau_3=\frac{1}{\beta\mu_{max}}\left[\mu_{max}-\arccos\left(\frac{\cos\Theta_p\cos i}{\sin\Theta_p\sin i}\right)/2\pi-1+\phi_{\perp,+}\right],
\end{equation}
{after which more than half of all points of the initial precession
  curve are involved in the ambiguity of the predicted line position), and $
  \tau_3\approx4.4$ for SS 433. All these situations are illustrated
  in Fig. \ref{f:munl}. Finally, full 'phase mixing' takes place at $
  \tau\approx2\tau_3\sim9$.

\begin{figure}
\centering
\includegraphics[width=0.4\textwidth,bb=80 200 530 680]{./figs/eTau18fnl.pdf}
\includegraphics[width=0.4\textwidth,bb=80 200 530 650]{./figs/eTau24fnl.pdf}
\includegraphics[width=0.4\textwidth,bb=80 180 530 650]{./figs/eTau44fnl.pdf}
\caption{\small The same as Fig. \ref{f:mu} but for the regime when photons from several blobs launched at different precession phases arrive at the observer at the same time, giving rise to ambiguity in the predicted X-ray line positions. The values of $ \tau$ are 1.8, 2.4 and 4.4 (from top to bottom).}
\label{f:munl}
\end{figure}

The steepening of the $\mu(\phi_0)$ curve at $ \phi_{\perp,-}\approx0.66 $ for  $\tau$ close $ \tau_1\approx1.8 $ (and at $ \phi_{\perp,+}\approx0.34$ for the red jet) means that significant brightening of the jet emission must occur at this phase (if the emission is spatially unresolved), since blobs that were launched over $\approx20\%$ of the precession period will become bright at this moment simultaneously (see Fig. \ref{f:munl}). The arising X-ray line will be significantly broadened. The same effect should also take place for other $ \tau $ but the associated 'amplification factor',
\begin{equation}
A(\tau,\phi_0)=|\left(\frac{d\phi_e(\tau,\phi_0)}{d\phi_0}\right)^{-1}|=|\left(1-\frac{d\Delta\phi(\tau,\phi_0)}{d\phi_0}\right)^{-1}|,
\end{equation}
is smaller (see Fig. \ref{f:ampl}).

In the regime when a set or range of $\phi_e$ values, $ [\phi_e] $, correspond to the same $\phi_0$, one needs to perform summation (integration) over this set (range) with the weights of different $ \phi_e $ determined by the square of the Doppler boosting factor $ \delta(\mu_e)$ (for the photon intensity from a blob with given $ \mu_e $).

 Thus, some property $S$ of the observed X-ray line can be calculated as
\begin{equation}
\tilde{S}(\tau,\phi_0)=\int_{[\phi_e]}d\phi_e~A_{\phi_e}(\tau,\phi_0)~\delta^2(\mu_e(\phi_e))~S(\phi_e)/w_0,
\end{equation}
where
\begin{equation}
w_0(\tau,\phi_0)=\int_{[\phi_e]}d\phi_e~A_{\phi_e}(\tau,\phi_0)~\delta^2(\mu_e(\phi_e))
\end{equation}
is the normalising factor.

\begin{figure}
\centering
\includegraphics[width=0.4\textwidth,bb=70 200 530 680]{./figs/ampl.pdf}
\caption{\small Amplification factor for the intensity of the
  unresolved emission due to simultaneous arrival of photons emitted
  by parts of the jet launched at different precession phases, as a
  function of the precession phase of the central source, for
  different values of the time-delay parameter $ \tau $ as indicated
  by the numbers above each curve's peak (coloured alternately blue
  and red for better visibility). A characteristic double-peaked
  structure emerges for $ \tau>1.8$ reflecting two 'turnover' points
  in the line position curve (see Fig. \ref{f:munl}). }
\label{f:ampl}
\end{figure}

Let us now consider a more general situation when the jet's matter
becomes X-ray bright not at one particular $ \tau $, but for $ \tau $
ranging from $ \tau_s$ to $ \tau_f$, and one cannot spatially resolve
the signal corresponding to different $\tau$ and $ \phi_e $ and
measures some integral characteristics of the line $<S> $ (primarily,
position and the total intensity) in the spectrum of such emission,
which can be calculated as
\begin{equation}
<S>_{ \tau_s,\tau_f}(\phi_0)=\int^{\tau_f}_{\tau_s}d\tau w_0(\tau,\phi_0)\tilde{S}(\tau,\phi_0)/W_0,
\label{eq:long}
\end{equation}
with the normalising coefficient
\begin{equation}
W_0=\int^{\tau_f}_{\tau_s}d\tau w_0(\tau,\phi_0).
\end{equation}
Figure \ref{f:mut} shows the resulting effect for various $ [\tau_s,\tau_f]$ intervals defined by their width $ \Delta\tau=\tau_f-\tau_s $ and the central point $<\tau>=(\tau_s+\tau_f)/2$.

As could be expected, the amplitude of the observed modulation with
precession phase decreases with increasing width of the averaging
window $ \Delta\tau $. However, the modulation remains noticeable even
for $\Delta\tau=0.9$, due to the Doppler boosting on the segment of the
modulation curve with the steep variation in the predicted line
centroid energy. For the receding jet, the situation is inverted,
which results in a significantly smaller amplitude of the modulation
of the line centroid with precession phase.

Since $\tilde{S}(\tau,\phi_0)$ is given by integration over $ \phi_e
$, one can change the integration order in Eq. (\ref{eq:long}), namely
first average the amplification factor  $ A_{\phi_e}(\tau,\phi_0) $
over $\tau$ and thus significantly weaken its influence on the net
result.  For the duration of the brightening period $ \Delta\tau >
(\mu_{max}-\mu_{min})\beta\tau\approx 0.2\tau $, the effect of the
relative amplification is already quite small, so the exact solution
is well approximated by the prediction that only includes Doppler
boosting when integrating over the whole emitting region (see the
dash-dotted line in Fig. \ref{f:mut}).

Before comparing our predictions with the observed line positions, it
is worth mentioning that since $ \mu\beta  <<1 $, $ \delta(\mu) $ is
well approximated by the linear function $(1+\beta\mu)/\gamma $, which
implies that the predicted variation of the line position $
E=\delta(\mu_e)E_{0}$ must closely follow the dependence of $ \mu_e $
on the precession phase.

\begin{figure}
\centering
\includegraphics[width=0.4\textwidth,bb=80 200 530 680]{./figs/et15dt03wx.pdf}
\includegraphics[width=0.4\textwidth,bb=80 200 530 650]{./figs/et15dt06wx.pdf}
\includegraphics[width=0.4\textwidth,bb=80 180 530 650]{./figs/et15dt09wx.pdf}
\caption{\small Averaged cosine of the angle between the jet's
  direction and the line of sight for the scenario of long-duration
  brightening as a function of the central source's precession phase
  shifted by $<\tau>$. The blue and red solid lines show the
  predictions for the blue and red jets respectively, taking into
  account both Doppler boosting and arrival time overlaps, while the
  dash-dotted lines show the prediction accounting for Doppler
  boosting only. The predicted line positions for $\tau=<\tau>$ in the
  short flare scenario with and without allowance for the
  light-travel-time effects are shown with dashed and dotted lines,
  respectively.
}
\label{f:mut}
\end{figure}
\section{Comparison with the data}
\label{s:comparison}

\cite{Migliari2002} reported the detection of X-ray emission lines at
energies $ E_{b,o}=7.28^{+0.02}_{-0.23} $ keV and $
E_{r,o}=6.39^{+0.12}_{-0.15} $ keV, apparently associated with the
blue and red jets, respectively, at distances $ \approx
1.5\times10^{17}$ cm $\approx1.4\lambda$ (in projection) from the
central source (see also \cite{Migliari2005}), which translates to $
\tau =1.4-1.6$ depending on $ \mu_e $. Since the duration of the
relevant (as well as subsequent ones) \textit{Chandra} observation
of the arcsec-scale X-ray emission is a few tens of kiloseconds
(see Table 1 in \cite{Migliari2005}), i.e. less than $ \sim0.1 $\% of
the precession period, these observations can be considered
instantaneous snaphots of the extended jets. The SS 433 precession
phase for the \cite{Migliari2002} observation was $\phi_0=0.63 $ (see
Table 1 in \cite{Migliari2005}, but note the 0.5 shift between the
definition of the zero phase in their work and the one adopted here).

Below, we first check whether the observed X-ray line positions can be reconciled with the predictions of the kinematic model taking into account light-travel-time effects. We consider four possible line identifications: Fe XXV K$\alpha$ (6.7 keV) and Fe XXVI Ly$\alpha$ (6.96 keV) -- the two brightest iron lines from hot optically thin plasma; and Fe I K$\alpha$ (6.4 keV) and Fe I K$\beta$ (7.06 keV) -- the brightest fluorescent lines of neutral iron.

\begin{figure}
\centering
\includegraphics[width=0.35\textwidth,bb=50 200 550 670]{./figs/eeTau15e67nl.pdf}
\includegraphics[width=0.35\textwidth,bb=50 200 550 650]{./figs/eeTau15e696nl.pdf}
\includegraphics[width=0.35\textwidth,bb=50 200 550 650]{./figs/eeTau15e64nl.pdf}
\includegraphics[width=0.35\textwidth,bb=50 180 550 650]{./figs/eeTau15e706nl.pdf}
\caption{\small Predicted line positions (solid lines) for the blue
  and red jets in the short brightening scenario with $ \tau=1.5$. The
  line's rest-frame energy $ E_0 $ equals to 6.7, 6.96, 6.4 and 7.06
  keV (from top to bottom). The dotted line shows the prediction of
  the kinematic model without allowance for light-travel-time
  effects. The hatched regions depict the 1$ \sigma$ uncertainty
  regions for the line centroids measured by \cite{Migliari2002}
  (solid horizontal lines) at the precession phase marked by the
  vertical dash-dotted line.}
\label{f:ee}
\end{figure}

Figure~\ref{f:ee} shows the predicted variations of the positions of
these four lines in the short brightening scenario. Evidently, for $
\tau\approx1.5$, the observed X-ray line positions are consistent only
with the Fe XXV K$\alpha$ line at a phase close to $\phi_0-\tau
\approx\phi_{\perp,-}$, i.e. $ \phi_0\approx0.2 $. For Fe I K$\alpha$,
agreement is possible only for $ \tau>\tau_2=2.4$, which, however,
contradicts the observed location of the X-ray emission at $ \tau<2
$. Similarly, for the Fe XXVI Ly$\alpha$ and Fe I K$\beta$ lines, a
match is possible for $ \tau\approx3.5 $ or $ \tau\approx0.2$, thus also
inconsistent with the observations.

In order to quantify the quality of the model fit to the observed line
positions, let us introduce the statistics
\begin{equation}
\chi^2(\tau,E_0)=\left(\frac{E_{b,m}-E_{b,o}}{\sigma(E_{b,o})}\right)^2+\left(\frac{E_{r,m}-E_{r,o}}{\sigma(E_{r,o})}\right)^2,
\label{eq:chi}
\end{equation}
where $E_{b,m}$ and $E_{r,m}$ are the model predictions for the blue
and red jets, and  $ \sigma(E_{b,o}) $ and $ \sigma(E_{r,o}) $ are the
upper and lower uncertainties for the measured values.

Projection of the parameter space region encompassed by the $
\Delta\chi^2(\tau,E_0)=1$ contour on the $E_0$ axis with $ \tau $ in
the range from 1.1 to 1.8 gives $ E_0\approx 6.7\pm0.2$ keV, where
the uncertainty is $ 1\sigma$. Thus, in the short brightening
scenario, the observed line positions are consistent with the
prediction for the Fe XXV K$\alpha$ line, whereas the other considered
lines cannot provide fits of comparable quality. However, for $
E_0=6.7 $ keV, the best fitting phase $ \phi_0\approx 0.2 $ is in
disagreement with the precession phase of the central source for the
\textit{Chandra} observation. This motivates us to consider the more
complicated long brightening scenario.

The predicted line positions in the long brightening scenario are
shown in Fig. \ref{f:eet} for $ <\tau>=1.5$ (i.e. $r\approx 1.6\times
10^{17}$ cm from the central source) and $ \Delta\tau=0.3,0.6,$ and $
0.9$ (corresponding to the emission region's size of $\approx
3.3\times 10^{16}$, $6.55\times 10^{16}$ and $9.8\times 10^{16}$ cm,
respectively). We see that for $ <\tau>=1.5$, almost no dependence on
$\Delta\tau$ is expected for the predicted line position for the blue
jet at the precession phase corresponding to the observation
discussed. This allows us to exclude both the $ E_0=6.4$ keV and $
E_0=6.7$ keV lines as possible identification candidates in this
scenario. Our predictions for the $ E_0=6.96$ keV and $ E_0=7.06$ keV
lines are broadly consistent with the observations for $\Delta\tau$
from $ \approx $0.45 to $\approx$1 in the former case and $\Delta\tau$
from $ \approx $0.4 to $\approx$0.75 in the latter. The best fitting
$\Delta\tau$ is $\approx0.6$ in both cases.

The predicted line widths calculated as the dispersion of the
predicted line positions for emission coming from different blobs
observed simultaneously equals $ \sigma(E)\approx 0.3$ keV for both
the blue and red jets if $\Delta\tau=0.6$. For $\Delta\tau=0.9$, $
\sigma(E)\approx 0.5$ keV and $ \sigma(E)\approx 0.4$ keV for the blue
and red jets, respectively. Our results are thus broadly consistent
with the results of the simulation performed by \cite{Migliari2002},
who found $ E_0=7.06$ keV for the $\Delta\tau=1$ scenario, with the
predicted line widths of $ \sim $ few $ 0.1 $ keV.

\begin{figure}
\centering
\includegraphics[width=0.35\textwidth,bb=50 200 550 670]{./figs/eet67t15dt03wx.pdf}
\includegraphics[width=0.35\textwidth,bb=50 200 550 650]{./figs/eet696t15dt03wx.pdf}
\includegraphics[width=0.35\textwidth,bb=50 200 550 650]{./figs/eet64t15dt03wx.pdf}
\includegraphics[width=0.35\textwidth,bb=50 180 550 650]{./figs/eet706t15dt03wx.pdf}
\caption{\small Same as Fig.\ref{f:ee} but for the scenario with a long brightening period centered at  <$\tau$>=1.5  and having duration $d\tau=\tau_f-\tau_f$ equal to 0.3 (short-dashed line), $0.6$ (long-dashed line), and $0.9$ (solid line). For comparison, the prediction for the short flare scenario with $ \tau=<\tau>$ is shown with the dotted line.}
\label{f:eet}
\end{figure}

\section{Discussion and conclusions}
\label{s:discussionconslusions}

\cite{Migliari2002} proposed that the observed X-ray lines are
associated with the 7.06 keV, i.e. Fe I K$\beta$, line. However, such
a line can arise as a result of fluorescence of cold gas irradiated by
X-rays, rather than from a hot ionized plasma as was suggested by
\cite{Migliari2002}. Moreover, if the observed lines were indeed
associated with Fe I K$\beta$, there must also be $ \sim 5-10$ times
stronger Fe I K$\alpha$ lines present in the spectrum, which are not
observed. We have demonstrated that the observed line positions are in
fact consistent with the Fe XXVI Ly$\alpha$ (6.96 keV) line, with
different Doppler shifts for the blue and red jets. Physically,
this implies that the emitting plasma should be re-heated up to
temperatures $ \gtrsim 10$ keV at $\sim 10^{17}$~cm from the central
source in a region of size $\sim 6\times 10^{16}$ cm. If the
temperature of the emitting plasma were somewhat lower, the centroids
of the observed lines would be shifted to lower energies as a result
of the contribution from the Fe XXV K$\alpha$ (6.7 keV) line to the
spectrally unresolved broad emission line.

\cite{Migliari2002} proposed that internal shocks arising due to variation in the speeds of blobs launched in approximately the same direction, i.e. with a time delay of $ \approx P$, can provide a mechanism for the jets' re-heating. Based on the detected short-term (on the time-scale of hours and days) variability, \cite{Migliari2005} further hypothesised the existence of an underlying outflow of highly energetic particles driving a fast shock wave. Another mechanism can be related to heating of the jet's matter as a result of interaction with the surrounding cocoon, which causes an oblique shock propagation inside the jets \citep{Heavens1990}.

At distances of $ \sim 10^{17} $ cm from the central source, the
matter of the jets could in fact be composed of a large number of
dense blobs, similarly to what is believed to occur in the region of
SS 433 optical line emission \citep{Fabrika2004}. From the energetic
considerations, the number density of the gas responsible for the
extended X-ray emission needs to be quite high:
\begin{equation}
n\sim 1\times 10^6~L_{X,33.5}/N_{50}/\Lambda_{-22.5} ~ {\rm cm^{-3}},
\end{equation}
where $ L_{X,33.5}\approx1$ is the observed luminosity of the extended
X-ray emission in units of $ 10^{33.5} $ erg/s \citep{Migliari2002},
$\Lambda_{-22.5}\approx1$ is the plasma emissivity in the same spectral
band for $ T\sim10  $ keV in units of $ 10^{-22.5} $ erg/s cm$^3$, and
$ N_{50}=k\dot{M}_j\Delta\tau P/m_p/10^{50}\sim k\Delta\tau
$ is the total number of hydrogen nuclei (divided by $10^{50}$) in the emitting region of
size $ \Delta\tau $ (along the jet) for the mass flux through the jet
$ \dot{M}_j\sim10^{-6.5}~M_{\odot}/$yr \citep{Fabrika2004}, and $k$ is
the mass fraction of the emitting gas.

The mean density of a continuous jet at this distance would be $ \sim
10^{2} $ cm$^{-3}$ (assuming the X-ray jets in the compact core are
characterized by the size $\sim 10^{11}$ cm and number density $\sim
10^{14}$ cm$^{-3}$, e.g. \citealt{Marshall2002, KMS2016}), i.e. less than the
required density by at least 4 orders of magnitude. If the length of
the X-ray bright period $ \Delta\tau P$ is determined by the cooling
time of the emitting matter $\sim T/(n\Lambda)$, then it is easy to
find that the required $ k$ is $\approx 0.03$, i.e. 3 \% of the total
mass of the jet's matter in this region. This mass could be provided
by a number of dense blobs experiencing collisions as a result of more
efficient deceleration of leading blobs by the interstellar medium
compared to trailing ones in a regime where the boundary of the jet
becomes virtually perpendicular to the local gas velocity direction,
i.e. at distances comparable with $\lambda/\pi\tan\Theta_p\approx
0.8\lambda$ \citep{Heavens1990}.

The characteristic collision velocities are $\sim\delta V\sim \Theta
_{j}\beta~c\sim 1000$ km/s, which should result in heating of the blobs
up to temperatures $ \sim 10 $ keV, similar to the observed
temperatures at the base of the X-ray jets near the central source
\citep{Marshall2002}. Blob collisions are possible only during the
period $ \delta t_{rc}\sim r\Theta_{j}/\delta V\approx P $
until a complete lateral re-collimation of the precessing jet takes
place and its shape becomes that of a band. After that, the X-ray
emission should gradually cease at the cooling time-scale.

Based on the constraints on the origin of the observed lines obtained
here, namely the requirements of high temperature of the emitting
plasma and relatively long duration of the brightening episode $
\Delta \tau \sim 0.6$, and taking into account the observed rapid
(down to scales of several hours, \citealt{Migliari2005}) stochastic
variability of the X-ray emission, we suggest that dissipation of the
jets' kinetic energy via collisions of their discrete dense blobs due
to their deceleration by the medium of the jets' cocoons, provides a
natural explanation of the observed characteristics of the
arcsec-scale X-ray emission.

\section{acknowledgements}
The research was supported by the Russian Science Foundation (grant 14-12-01315).


\section{{\normalsize{}Introduction}}

After the synthesis of graphene \cite{graphene_synth}, and the discovery
of its unique electronic, optical and thermal properties, research
on low-dimensional materials has increased tremendously. Graphene
is attractive because of its unique properties such as its high electronic
conductivity, and high mobility at room temperature \cite{gn1,gn2,gn3}.
However its applications are limited in semiconducting devices because
of its zero band gap. This gaplessness of graphene has motivated the
discovery of alternate nanomaterials that possess finite band gap,
and have the potential to replace the silicon in semiconductor technology.
Recently two-dimensional (2D) materials which have finite band gap
such as monolayers of hexagonal boron-nitride (h-BN) \cite{bn_synth},
transition metal chalcogenides (MoS$_{2}$, WS$_{2}$ etc.)\cite{mos2,opt_tmdc},
and phosphorene have been successfully synthesized\cite{phospho}.

Quasi one-dimensional materials such as nanoribbons, nanotubes, nanorods,
and nanowires have also attracted attention because of their interesting
photochemical, photophysical, and transport properties\cite{nr1,nr2}.
Nanoribbons (NRs), in particular, have received great attention from
researchers because of their unique electronic properties which can
be modified on the basis of their edge configuration, and width. Quantum
confinement because of the presence of edges, and their different
possible passivations, lead to interesting optical, electronic and
magnetic properties. Graphene nanoribbons (GNRs) and boron-nitride
nanoribbons (BNNRs) have been experimentally synthesized by unwrapping
of carbon nanotubes, and BN nanotubes, respectively\cite{gnr_synthl,bnnr_synth}.
Armchair GNRs (AGNRs) are non magnetic semiconductors for all widths,
and have oscillating band gaps over families, approaching the zero
band gap of a 2D sheet, for large widths\cite{AGNR_Eg_osc}. Zigzag
GNRs (ZGNRs), on the other hand, have tunable band gaps ranging from
metal to semiconductor, depending on the width, and passivation edges\cite{O_zgnr,pass_zgnr}.
Furthermore, ZGNRs also exhibit half-metallic behavior with spin-polarized
band gaps when an electric field is applied along the direction of
the width of the ribbons, and their gaps can be controlled by the
strength of the electric field\cite{ZGNR_Efield}.

\textcolor{black}{C and Si have same number of electrons, therefore
SiC structures are stable in bulk and have electronic properties which
are intermediate in between bulk Si and carbon structures. Silicon
carbide (SiC) crystallizes in several forms such as hexagonal, rhombohedral,
and cubic Bravais lattices\cite{edge_H_zig}. Bulk SiC is a wide band
gap semiconductor, and has been used in high-temperature, high-pressure
and high-frequency device applications\cite{SiC_app,SiC_Eg}. The
graphene-like 2D SiC monolayer has not been experimentally fabricated
yet, but has been extensively studied using theoretical methods}\cite{SiC_mono,ASiCNR_Eg_osc,sic_bare,sic2d_eg,SiX,sic_bulk_gw,SiC_BN_BeO,group4,opt_SiCNT,tight_bind}\textcolor{black}{.
Unlike graphene, it is a direct band gap semiconductor with a band
gap of 2.5 eV. An ultra-thin SiC nanosheet which has a thickness of
$0.5-1.5\text{ nm}$ has been fabricated, and is used in light emitting
applications }\cite{SiC2D}\textcolor{black}{. Various SiC structures
such as SiC nanotubes, SiC nanowires, microribbons, and crystalline
and bicrystalline nanobelts have also been successfully synthesized\cite{SiCNT_synth,sicnr_micro,3C_nanobelts,bicrystaline_nanobelt}.
Furthermore multilayer SiC nanoribbons in lengths of micrometers,
and thicknesses of nanometers have also been synthesized}\cite{SiCNR_synth}\textcolor{black}{. }

As far as theory is concerned\textcolor{black}{,{} SiC nanoribbons
have been extensively studied using }first-principles calculations,
based on the density functional theory approach (DFT)\textcolor{black}{{}
}\cite{ASiCNR_Eg_osc,ZSiCNR_DFT,ZSiCNR_Efield,assymetry_H,edge_H_zig,sic_bare,O_S_zig,mod_eg,B_N_SiCNR,half_bare_zig,vac_SiCNR}\textcolor{black}{.
}Sun et al. \cite{ASiCNR_Eg_osc} \textcolor{black}{reported that
hydrogen- passivated armchair SiCNRs (ASiCNRs) are direct band gap
semiconductors, with band gaps in the range of 2.3 - 2.4 eV. Furthermore,}
\textcolor{black}{they reported that band gaps of zigzag SiCNRs (ZSiCNRs)
are spin polarized, and exhibit intrinsic half-metallic behavior,
without the application of an external electric field.}\cite{ASiCNR_Eg_osc}
Lou \emph{et al}.\cite{ZSiCNR_DFT} studied narrow ZSiCNRs with widths
in the range 0.6\textendash 1.6 nm, and showed that the ferrimagnetic
state is more stable,\textcolor{cyan}{{} }with an anti parallel spin
orientation between the two edges. ZSiCNRs that have a half-bare-edge
were studied by Tang \emph{et al.}\cite{half_bare_zig} using hybrid
DFT approach. They reported that half-bare-edge ZSiCNRs with a bare
carbon edge are magnetic semiconductors, while the ones with bare
Si edge atoms are magnetic metals. The influence of doping with B
and N impurities has been studied by Costa \emph{et al.},\cite{B_N_SiCNR}
who showed that B doped ZSiCNRs retain half metallic behaviour, while
N doped ZSiCNRs become metallic. Lopez-Benzanilla\emph{ et} \emph{al.}\cite{O_S_zig}
employed a local spin-density approximation (LSDA) approach to study
sulfur and oxygen passivated ZSiCNRs to discover that, ZSiCNRs turned
from half-metallic to semiconductors, or metals, as a result of passivation.
The effects of substitution of edge atoms, with B and N atoms, has
been studied by Zheng \emph{et al.}\cite{mod_eg} who reported that
modified ZSiCNRs are semiconductors. Bekaroglu \emph{et al.,}\cite{sic_bare}
and Morbec \emph{et al.},\cite{vac_SiCNR} studied the effects of
vacancies on the electronic structure of ZSiCNRs by using the DFT-GGA
approach, and found that double (Si and C) vacancies induce magnetism.\textcolor{black}{{}
Ultra narrow ZSiCNRs were studied by Ping et al.\cite{lou_short}
including short-range exact exchange effects, and they obtained lower
edge energies, and higher band gaps, as compared to the ones computed
using GGA. }

\textcolor{black}{It is well known that the DFT-based approaches underestimate
the band gap in semiconductors, because they do not include electron
correlation effects.\cite{gw_acc} It has been demonstrated that in
nanoribbons, not including correlation effects causes severe errors
as compared to their bulk counterpart, because in low-dimensional
systems these effects are enhanced \cite{GNR_GW,GNR_GW2,GNR_GW3,GNR_GW4,GNR_GW5,AMoS2,Nalaal}.
It will be also be interesting to study many-body electron-electron
interactions in evaluating accurate band gaps and electron-hole effects
in order to calculate optical absorption spectra in ZSiCNRs; this
knowledge will enable the exploration of their potential in optical
device applications. In the present work, we calculate the quasiparticle
band structures and the optical absorption spectra of ZSiCNRs by using
GW approximation, and by solving the Bethe-Salpeter equation (BSE),
respectively. We find that ultra-narrow ZSiCNRs, with widths narrower
than 1 nm, are non-magnetic semiconductors. GW corrections open up
significant band gaps in ZSiCNRs that are wider than 1 nm, and transform
them from half-metals to semiconductors. We find that narrowest ZSiCNR,
with the width of 0.6 nm, has strongly bound excitons, with a binding
energy of 2.1 eV, while other nanoribbons have weakly bound excitons
because their band gaps are too small. Therefore, the narrowest ribbon
can be utilized in optoelectronic applications in the infrared frequency
range. Broader nanoribbons exhibit distinct band gaps for the two
spin orientations, and, therefore, they can be utilized in spintronic
devices such as spin valves.}

\textcolor{black}{The remainder of this paper is organized as follows.
Section \ref{sec:theory} presents details of computations employed
in this work, while section \ref{sec:results} presents results of
our calculations on ZSiCNRs of varying widths. Finally, in section
\ref{sec:conclusions}, we present our conclusions.}

\section{{\normalsize{}Computational Details}}

\label{sec:theory}

\textcolor{black}{We carried out first-principles many-body electronic
structure calculations for ZSiCNRs using a three-step procedure. Firstly,
the DFT calculations are performed using a plane wave approach, in
the generalized gradient approximation (GGA),\cite{gga} with Perdew-Burke-Enzerhof
(PBE) pseudopotentials.\cite{psp} Geometry optimizations were carried
out using the computer program Vienna ab initio simulation package
(VASP) \cite{VASP}. All the nanoribbons considered here were fully
relaxed until the force on each atom was lower than 0.01 eV/\AA.
The threshold energy convergence was set to $10^{-4}$ eV. A Monokhorst
k-point grid of $1\times1\times11$ was used for Brillouin zone integration
during the relaxation. After the relaxation, we used optimized atomic
positions and lattice parameters to calculate Kohn-Sham band gaps
of ZSiCNRs using the DFT-GGA approach as implemented in the software
package ABINIT \cite{abinit1,abinit2}. GGA wave functions are a good
starting point for the many-body calculations, aimed at including
electron correlation effetcs.}

\textcolor{black}{In the second step we calculate quasiparticle corrections
for the DFT-GGA band gaps using the many-body approach G$_{0}$W$_{0}$
approximation\cite{g0w0}, which is a single-shot GW approach for
computing self-energy corrections in terms of the Green's function
(G), and the screened Coulomb potential (W). The screening effects
were included by using the plasmon pole model\cite{ppm}. A kinetic
energy cutoff of 26 Ha has been used for all the systems considered
here. Cutoff energies of 5 Ha, and 12 Ha have been used for the correlation,
and the exchange, parts of self-energy calculations, respectively.
A total number of 200 bands were used, out of which more than 120
were unoccupied. In the last step, the BSE is solved within the Tamm-Dancoff
approximation\cite{tda} in order to obtain the optical response inclusive
of the contribution of the electron-hole interactions, retaining only
the resonant part of the Bethe-Salpeter Hamiltonian: This is because
the inclusion of the coupling has been shown to have a negligible
effect. This approximation has been also used to study similar systems
such as graphene nanoribbons,\cite{GNR_GW,GNR_GW6} and boron-nitride
nanoribbons\cite{GNR_GW4}. In the GW, and the BSE calculations, a
Coulomb truncation scheme was also used to avoid the inclusion of
long range interactions between the periodic images. A $1\times1\times30$
k-point grid has been used in the the GW and BSE calculations for
the systems considered here. Five valence band and five conduction
bands were included to obtain the optical absorption spectra. Nanoribbons
were considered to be periodic along the $z$ direction, and more
than 13 \AA\ wide vacuum layers were taken in non-periodic directions
($x$ and $y$ directions), in order to represent an isolated system.
Width of a given ZSiCNR was specified by prefixing it with the number
of zigzag lines ($N_{z}$) across its width. Thus ZSiCNR-3 denotes
a ribbon with $N_{z}$=3 zigzag lines across the width. In the present
work, we have studied ribbons with $2\leq N_{z}\leq8$, which corresponds
to widths between 0.6 $\text{nm}$ and 2.2 $\text{nm}$.}

\section{{\normalsize{}Results And Discussion}}

\label{sec:results}

\subsection{Formation energy and stability}

\begin{figure}
\begin{centering}
\includegraphics{fig1}
\par\end{centering}
\caption{(a) Geometric structure of ZSiCNR-4 unit cell shown in the red dashed
box is considered for the calculation\textcolor{red}{. }(b) Calculated
spin-density of ZSiCNR-4. Blue and yellow colored spin densities denote
opposite spin orientations. (c) Computed edge formation energy of
ZSiCNRs as a function of width. }

\label{str_zsicnr}
\end{figure}

Fig. \ref{str_zsicnr} (a) shows the geometric structure of ZSiCNR-4.
Optimized bond lengths corresponding to Si-C, Si-H, and C-H were obtained
to be 1.77 \AA, 1.49 \AA, and 1.09 \AA, respectively, while the
lattice constant was found to be 3.11 \AA, for all the structures
considered here. These parameters are in good agreement with previous
studies\cite{ASiCNR_Eg_osc,lou_short}. Here, we calculate the edge
formation energies of these nanoribbons in order to understand their
thermodynamic stability, during edge formation. We evaluated the formation
energy by using the following formula:

\textcolor{black}{
\[
E_{f}=\left(E_{T}-N_{SiC}E_{SiC}-0.5N_{H}E_{H_{2}}\right)/2l,
\]
where $E_{T}$ is the total energy of the hydrogen passivated zigzag
SiC nanoribbon , $E_{SiC}$ is the total energy of the 2D SiC sheet,
$E_{H_{2}}$ is the total energy of an isolated hydrogen molecule,
$N_{SiC}$ is the number of SiC dimers in the unit cell, $N_{H}$
is the number of H atoms in the unit cell, and $l$ is the periodic
length of an edge. In general, the lower the value of $E_{f}$ of
a material, the higher its stability}. The formation energy of ZSiCNRs
as a function of width has been plotted in Fig. \ref{str_zsicnr}
(c). Unlike ASiCNRs,\cite{Nalaal} formation energy of these nanoribbons
is strongly dependent on their width, and it increases with the width
of the nanoribbons. From the figure it is obvious that ribbons which
have narrow widths have high stability, suggesting that they will
be relatively easier to synthesize in the laboratory.

\subsection{Nonmagnetic ZSiCNRs}

\subsubsection{Quasiparticle energies}

First we discuss the electronic structure of ZSiCNRs with widths $N_{z}=\text{2}$
and $\text{3}$, which do not exhibit any magnetic behavior. Ping
\emph{et al.} \cite{lou_short} reported that these two nanoribbons
are non-magnetic semiconductors, which do not exhibit half-metallic
behaviour within the GGA or HSE approaches, unlike ribbons which are
wider than 1 nm. Their DFT-GGA study predicted that ZSiCNR-2 is a
non-magnetic semiconductor with a band gap of 0.97 eV, and ZSiCNR-3
is essentially metallic with a band gap of 0.07 eV. Our calculated
DFT-GGA and GW band structures of ZSiCNR-2 and ZSiCNR-3 are shown
in Figs. \ref{band_z2z3} (a) and (b ) respectively. The blue lines
correspond to the DFT-GGA, while the red lines denote the GW band
structure. Our DFT results agree with those of Ping et al \cite{lou_short}.
Fig. \ref{band_z2z3} (a) shows that ZSiCNR-2 is a direct band gap
material at the DFT-GGA level, and the gap occurs at the $Z$ point
located at the edge of Brillioun zone. After including the many-body
effects, computed band structure is depicted by the red lines. We
found that quasiparticle corrections increase the band gap from 0.97
to 2.4 eV, which is an indirect one (quasi-direct) at the GW level.
In the GW band structure of ZSiCNR-2, the valence band maximum (VBM)
occurs small distance away from the $Z$ point, and there is no change
in the conduction band minimum (CBM). Fig. \ref{band_z2z3} (b) shows
the band structures computed by GGA and GW approaches, for ZSiCNR-3.
Self-energy corrections widened the band gap from 0.07 eV to 0.45
eV, which is a six fold increase due to electron correlation effects.
This increase shows that many-body corrections transform ZSiCNR-3
from a nearly metallic system, to a semiconductor.

\begin{figure}[H]
\begin{centering}
\includegraphics[scale=1.2]{fig2}
\par\end{centering}
\caption{(a) and (b) GGA and GW band structures of ZSiCNR-2 and ZSiCNR-3. Blue
and red lines represent GGA and GW band structures, respectively,
and black arrows denote interband transitions. (c) and (d) Quasiparticle
self-energy corrections (differences between quasiparticle energies
and GGA energy values) for various bands to the GGA Kohn-Sham energies
for ZSiCNR-2, and ZSiCNR-3, respectively. Corrections to the states
corresponding to the edge atoms in valence, and conduction, bands
are enclosed in the red dashed ellipses.}

\label{band_z2z3}
\end{figure}

Self-energy corrections for various bands for both the ribbons are
shown in Figs. \ref{band_z2z3} (c) and (d). In these nanoribbons,
$\pi$ band states corresponding to the $p_{x}$ orbital extend into
the vacuum perpendicular to the ribbon plane, while the $\sigma$
states are composed of $p_{z}$ and $p_{y}$ orbitals lying in the
plane of the sheet. The self-energy corrections to the $\pi$ states
are larger than those to the $\sigma$ states because of the screening
effects due to other $\pi$ electrons . The corrections enclosed in
the red ellipses correspond to the edge states which experience enhanced
Coulombic interactions. When compared to GGA band structure, GW bands
not only shift but also stretch by an average of 20\%-30\% because
of non-uniform quasiparticle corrections. \textcolor{black}{We find
that corrections to the states which derive main contribution from
the interior atoms are less linear when compared to states which are
predominantly localized on edge atoms}\textcolor{black}{\uline{.}}

\begin{figure}[H]
\begin{centering}
\includegraphics[scale=0.6]{fig3}
\par\end{centering}
\caption{\textcolor{black}{Projected density of states for ZSiCNR-2 and ZSiCNR-3.
Red and blue curves show states which derive contributions from the
edge C and Si atoms, while black and green curves represent states
which derive main contributions from internal C and Si atoms, respectively.
Fermi level is represented by black dashed line, which has been set
to 0 eV.}}

\label{fig:pdos_z2_z3}
\end{figure}
 \textcolor{black}{We present projected density of states (PDOS) for
ZSiCNR-2 and ZSiCNR-3 in Fig. \ref{fig:pdos_z2_z3}, in order to illustrate
the contribution of various atoms from the edges, and the interior,
to different orbitals. PDOS plots show that the VBM in both the ribbons
arises from edge carbon atoms, while the CBM is derived from edge
Si atoms. The states located on internal C and Si atoms are away from
the Fermi-level in both the ribbons.}

\subsubsection{Optical absorption spectra}

\begin{figure}[H]
\begin{centering}
\includegraphics[scale=0.5]{fig4}
\par\end{centering}
\caption{\textcolor{black}{Optical absorption spectra of (a) ZSiCNR-2, (b)
ZSiCNR-3. The blue curve shows the spectrum computed without the contribution
of electron-hole interactions, while the red curve represents spectrum
calculated by including electron-hole interactions.}}
\centering{}\label{fig:bse_zsicnr}
\end{figure}

\textcolor{black}{Using the quasiparticle band structure, we calculated
the optical absorption spectrum of the ZSiCNR-2 and ZSiCNR-3, with,
and without electron-hole interactions as is shown in Figs. \ref{fig:bse_zsicnr}
(a) and (b). Electron-hole interactions were incorporated by using
the BSE, and in all the calculations the polarization direction of
incident light was considered to be along the length of the nanoribbons.
We use $E_{g,d}$ notation to label the direct transitions which occur
between first valence band and first conduction band as shown in Fig.
\ref{band_z2z3}. The blue curve represents absorption due to interband
transitions at the GW-RPA level, while the red curve is a result of
solving the BSE, including electron-hole effects. In both the ribbons,
the first prominent peak comes from the transition between the first
valence band, and the first conduction band. In GW-RPA absorption
spectra we can observe that there are several peaks due to different
interband transitions in ZSiCNR-2. However, there is only one prominent
peak in ZSiCNR-3 which is a result of the transition between the VBM
and the CBM. The inclusion of electron-hole interactions changes the
spectrum significantly because, not only the positions of the peaks
are shifted, but also their shape looks different, when compared with
the GW-RPA absorption spectrum. In the absorption spectra of both
the ribbons, which include electron-hole interactions, we observe
only one intense excitonic peak due to the transition between first
valence and first conduction band, $E_{g,d}$. We did not find any
other transitions among other bands, even though we included five
valence and five conduction bands in our calculations. In Fig. \ref{fig:bse_zsicnr}
it is obvious that ZSiCNR-2 has a strongly bound exciton, of binding
energy of 2.1 eV ($E_{g,d}$), but with relatively weaker intensity,
corresponding to the first absorption peak. We also observed another
strongly bound excitonic peak ($E_{g,d}^{1}$), with a very high intensity,
and binding energy of 1.5 eV, originating at distinct $k$-points
as compared to the first one. These binding energies are significantly
larger than the calculated exciton binding energy, 1.17 eV, for the
2D SiC sheet,\cite{sic_bulk_gw} suggesting that the reduced dimensionality,
and quantum confinement, lead to an enhancement in the exciton binding
energies of nanoribbons. From among all the ribbons considered here,
it is clear that ZSiCNR-2 is a small gap semiconductor which can be
used in optical applications in the low energy region of spectrum.
Unlike ZSiCNR-2, strongly bound excitons are not observed in ZSiCNR-3
because the band gap is much smaller for this ribbon. Band to band
absorption spectrum peak at GW level gets blue shifted upon inclusion
of electron-hole effects, as shown in Fig. 4. The final location of
the blueshifted absorption peak clearly implies negative exciton binding
energy. Thus, the nature of optical transitions contributing to various
peaks is also modified in ZSiCNR-3, once the many-body effects are
included.}

\subsection{Spin-polarized ZSiCNRs}

\begin{figure}[H]
\begin{centering}
\includegraphics[scale=0.5]{fig5}
\par\end{centering}
\caption{\textcolor{black}{(a) GGA and (b) quasiparticle band structures of
ZSiCNR-4. Red and blue lines represent spin-up and spin-down states,
respectively. The atoms contributing to the states of the valence
and conduction band on either side of the band gap are marked for
ZSiCNR-4. Towards the zone boundary, the states are localized at the
ribbon edges.}}

\label{band_znr4}
\end{figure}

\textcolor{black}{Next we present the results of quasiparticle band
structure calculations of spin-polarized ZSiCNR-4, because $N_{z}=4$
to $N_{z}=8$ have similar band structures, and at mean-field DFT
level exhibit intrinsic half-metallic behavior.\cite{ASiCNR_Eg_osc,ZSiCNR_DFT}
In Figs. \ref{band_znr4} (a) and (b) we compare the GGA and GW quasiparticle
band structures of ZSiCNR-4, and show the widened quasiparticle band
gap, when electron-correlation effects are included. While DFT-GGA
predicts a very small gap of 30 meV, the GW quasiparticle gap is seen
to have widened to 0.33 eV, clearly showing that the ZSiCNRs are narrow-gap
semiconductors (Table \ref{tab-gaps1}). Although, after including
many-body effects through GW calculations, half-metallicity disappears,
but still the band gaps for two spin orientations are different. The
atoms contributing to various states are marked in Fig. \ref{band_znr4}
(a). Towards the zone boundary, the states of the first valence and
conduction bands of both spins are localized at the ribbon edges.\cite{ZSiCNR_DFT}
From the spin-density plot presented in Fig. \ref{str_zsicnr} (b),
it is obvious that the spins localized on the two edges have opposite
spin orientations. The states to the left of the VBM/CBM extend over
the entire ribbon width, except for the spin-down conduction band
which mainly comprises of the Si edge states. The screening at the
edges shows a complicated dependence on the states, which is enhanced
in ZSiCNRs because the charge and spin distributions differ at the
Si and C edges. Fig. \ref{fig:structure} shows that the self energy
corrections are highly state-dependent. The states which are enclosed
in rectangular boxes, away from the main region, represent the spin-up
edge states (valence states at the Si edge, and conduction states
at the C edge). As the wave vector approaches the zone boundary (Z
point), the states get more localized at the edges,\cite{ZNR_edge_states}
and the deviations between up- and down-spin states, as also the energy
gap for up-spin states, become larger (}\textcolor{black}{\emph{cf.}}\textcolor{black}{{}
Fig. \ref{band_znr4}(a)). Thus, the nature of spin-up bands is very
similar to the bands in spin degenerate zigzag graphene} nanoribbons.\cite{GNR_GW3}

\begin{table}
\begin{tabular}{ccccccc}
\hline
 &  & spin-up  &  &  & spin-down  & \tabularnewline
\hline
\textcolor{black}{N$_{z}$ } & $E_{GGA}^{gap}$  & $E_{GW}^{gap}$  & $\Delta E^{gap}$  & $E_{2\ GGA}^{gap}$  & $E_{2\ GW}^{gap}$  & $\Delta E_{2}^{gap}$\tabularnewline
\hline
4  & 0.03  & 0.33  & 0.30  & 0.55  & 0.59 & 0.04\tabularnewline
5  & 0.06  & 0.82  & 0.76  & 0.72  & 0.77 & 0.05\tabularnewline
6  & 0.006  & 0.38  & 0.37  & 0.64 & 0.74  & 0.1\tabularnewline
7 & 0.002 & 0.384 & 0.38 & 0.64 & 0.92 & 0.28\tabularnewline
8 & 0.08 & 0.77 & 0.69 & 0.76 & 0.93 & 0.14\tabularnewline
\hline
\end{tabular}

\caption{\textcolor{black}{The values of GGA and GW band gaps, and their difference,
for spin-polarized ZSiCNRs for widths N$_{z}$= 4 to 8. Columns 2-4
(Columns 5-7) correspond to spin-up (spin-down) channel. All energies
are in eV. \label{tab-gaps1}}}
\end{table}

\vspace{0.2in}

\begin{figure}[H]
\begin{centering}
\includegraphics[scale=0.7]{fig6}
\par\end{centering}
\caption{\textcolor{black}{Quasiparticle self-energy corrections to the GGA
Kohn-Sham energies of ZSiCNR-4. Blue and red empty circles represent
spin-up, and spin-down, states respectively. Corrections to the spin-up
edge states are shown in the dashed rectangular boxes. The dashed
ellipses enclose the $\sigma$ states.}}

\label{fig:structure}
\end{figure}

\vspace{0.2in}

The larger spin-down gap at the zone boundary exhibits a smaller quasi-particle
correction. Fig. \ref{fig:structure} shows the valence spin-down
states bending away from the spin-up states, and, similarly, spin-down
states of the conduction band too have smaller corrections than the
corresponding spin-up states. These states correspond to the valence
C, and the conduction Si, states localized at the edges. The spin-down
edge states show smaller self energy corrections when compared to
those localized in the interior, while the spin-up edge states exhibit
exactly the opposite behavior. The charge distributions are different
at the two edges, leading to varying amounts of screening. For example,
in the ground state, there is a partial valence electron transfer
from Si to C.\cite{SiCnr_polar} As a consequence, the edge C states
are screened more, and have smaller self energy corrections, while
the lesser screened edge valence Si states experience enhanced Coulomb
interactions. Therefore, the spin-down band gap does not change as
much as the spin-up band gap.

A few points enclosed by dashed circles lying outside the main region
are $\sigma$ states. Apart from the first valence and conduction
bands, the next few bands comprise both $\sigma$ and $\pi$ states
(shown in Fig. \ref{fig:structure}). Similar to ASiCNRs,\cite{Nalaal}
and narrow ZSiCNRs, the $\sigma$ states have smaller-self energy
corrections than the $\pi$ states.

\textcolor{black}{For the sake of completeness, we also performed
calculations on the non-magnetic (i.e., non-spin polarized) states
of ZSiCNRs, and found that spin polarized state does not exist for
narrow ribbons ($N_{z}=2,\:3$), and spin-polarized calculations for
those ribbons yield the same results as non-magnetic calculations.
However, for broader ribbons, both spin-polarized and non-magnetic
states exist, and quite expectedly, total energies for the spin-polarized
cases are lower as compared to the non-magnetic ones (see Table \ref{tab:mvsnm}).
We also performed GGA and GW calculations for the non-magnetic states
of broader ribbons and results are presented in Table \ref{tab:nonmag}.
From the table it is obvious that: (a) at DFT-GGA level band gaps
of the ribbons decrease rapidly with their width, and (b) quasiparticle
corrections are large for the ultra narrow ZSiCNRs, but they also
reduce in magnitude rapidly with the increasing width. Combining these
two results, we conclude that non-magnetic calculations predict that
broader ZSiCNRs will exhibit metallic behavior. Therefore, an experimental
measurement of the band gaps of broader ZSiCNRs can settle the issue
whether these systems exhibit a magnetic ground state, or a nonmagnetic
one.}

\begin{table}
\begin{tabular}{c|cc|c}
\hline
\multirow{2}{*}{\textcolor{black}{N$_{z}$}} & \multicolumn{2}{c|}{\textcolor{black}{Total energy (eV) }} & \multirow{2}{*}{\textcolor{black}{$\Delta E$ (eV)}}\tabularnewline
\cline{2-3}
 & \textcolor{black}{NM} & \textcolor{black}{SP} & \tabularnewline
\hline
\textcolor{black}{4 } & \textcolor{black}{-62.106675} & \textcolor{black}{-62.115140} & \textcolor{black}{0.008465}\tabularnewline
\textcolor{black}{5 } & \textcolor{black}{-76.101117} & \textcolor{black}{-76.120842} & \textcolor{black}{0.019725}\tabularnewline
\textcolor{black}{6 } & \textcolor{black}{-90.119364} & \textcolor{black}{-90.144896} & \textcolor{black}{0.025532}\tabularnewline
\hline
\end{tabular}

\caption{\textcolor{black}{Total energies per unit cell for non-magnetic (NM)
and spin-polarized (SP) ZSiCNRs for widths N$_{z}$= 4 to 6. Last
column gives energy difference between both the configurations}\textcolor{green}{.
\label{tab:mvsnm}}}
\end{table}

\begin{table}
\begin{tabular}{cccc}
\hline
\textcolor{black}{N$_{z}$ } & $E_{GGA}^{gap}$  & $E_{GW}^{gap}$  & $\Delta E^{gap}$ \tabularnewline
\hline
2 & 0.97 & 2.4 & 1.43\tabularnewline
3 & 0.07 & 0.45 & 0.36\tabularnewline
4  & 0.011 & 0.22 & 0.209\tabularnewline
5  & 0.015 & 0.15  & 0.135\tabularnewline
6  & 0.002 & 0.07  & 0.068\tabularnewline
\hline
\end{tabular}

\caption{\textcolor{black}{The values of GGA and GW band gaps, and their difference,
for non-magnetic states of ZSiCNRs for widths N$_{z}$= 2 to 6 . All
energy values are in eV.}\textcolor{black}{\uline{\label{tab:nonmag}}}}
\end{table}

\textcolor{black}{We have not presented optical absorption spectra
of spin-polarized ZSiCNRs as they undergo different quasiparticle
self-energy corrections because of two spin band gaps. We could expect
weakly bound excitons in ribbons that are wider than 1 $\text{nm}$.
Small band gap of ZSiCNRs which are wider than 0.6 $\text{nm}$ makes
them potentially useful in IR devices.}

\section{Conclusions}

\label{sec:conclusions}

\textcolor{black}{We performed first-principles many-body calculations
in order to investigate quasi particle band structures and optical
absorption of hydrogen-passivated zigzag SiC nanoribbons whose widths
ranged from 0.6 nm to 2.2 nm. For the study, computationally intensive
GW approximation was employed to compute the self-energy corrected
quasiparticle band gaps. We also used the BSE to calculate optical
absorption spectra including eletron-hole effects. It was found that
the many-body effects are significant due to the reduced dimensionality
of nanoribbons. Self-energy corrections transformed nearly half-metallic
zigzag SiC nanoribbons with width larger than 1 $\text{nm}$, to semiconductors.
The inclusion of electron-hole effects changed optical absorption
spectra in a significant manner, both qualitatively, and quantitatively.
The excitonic binding energies and luminescence properties are dependent
on width of the ribbon. We found that nanoribbon with a width of 0.6
$\text{nm}$ has strongly bound excitons with binding energy of 2.1
eV, suggesting its possible utility in optoelectronic applications.
We also computed the edge formation energy and showed that narrower
ribbons are more stable when compared to wider ones. The unique Coulombic
interactions at the spin-polarized edges make zigzag SiC nanoribbons
interesting systems to study. As far as behavior of their band gaps
with respect to the width is concerned, because the band gaps in half-metallic
ZSiCNRs arise due to the presence of magnetic edge states, therefore,
broader ribbons are expected to have spin-polarized non-zero band
gaps, hinting at their possible utility in spintronic applications,
such as spin valves. }
\begin{acknowledgments}
NA and NM gratefully acknowledge the support from Monash HPC, National
Computing Infrastructure of Australia, and the Pawsey Supercomputing
facility. This research was partially supported by the Australian
Research Council Centre of Excellence in Future Low-Energy Electronics
Technologies (project number CE170100039) and funded by the Australian
Government. AS acknowledges the financial support from Department
of Science and Technology, Government of India, under project no.
SB/S2/CMP-066/2013.
\end{acknowledgments}