\section{Green's functions under RWA}
\label{correct_formulae}
In Ref.~\onlinecite{Berakdar_12}, the Green's functions of graphene without
connecting the leads $g^{\tau\tau,r}_{s,s'}(\varepsilon)$
given by Eqs.~(A.3)-(A.8) are not the correct solution of Eq.~(A.2) in that
paper. In the following, we present the correct formulae of 
$g^{\tau\tau,r}_{s,s'}(\varepsilon)$ in the cases with $V_g=0$.
\begin{eqnarray}
  g^{\tau\tau,r}_{11}(\varepsilon)&=&-\frac{S}{2\pi v_{\rm F}^2}\bigg\{ \frac{1}{2}
  (\varepsilon+v_0)u_1 - \frac{\Delta^2-\varepsilon(\varepsilon+v_0)}
  {\sqrt{|\Delta^2-{\varepsilon}^2|}}
  \nonumber \\ &&
  {} \times \Big[ u_2\theta(\Delta-|\varepsilon|) 
  -u_3\theta(|\varepsilon|-\Delta) \Big] + D \bigg\} 
  \nonumber \\ && {}-\frac{iS\ {\rm sgn}(\varepsilon)}
  {4v_{\rm F}^2\sqrt{\varepsilon^2-\Delta^2}} \bigg[
  \theta\bigg(\Delta<|\varepsilon|<\sqrt{\frac{\omega_0^2}{4}+\Delta^2}\bigg)
  \nonumber \\ && {}\times
  \Big(\varepsilon-\sqrt{\varepsilon^2-\Delta^2}\Big)v_{-} 
  + \Big(\varepsilon+\sqrt{\varepsilon^2-\Delta^2}\Big)v_{+} 
  \nonumber \\ && {}\times
  \theta\bigg(\Delta<|\varepsilon|<\sqrt{(D-\frac{\omega_0}{2})^2+\Delta^2}
  \bigg) \bigg],
\label{gsmall_11}
\end{eqnarray}
\begin{eqnarray}
  g^{\tau\tau,r}_{22}(\varepsilon)&=&-\frac{S}{2\pi v_{\rm F}^2}\bigg\{ \frac{1}{2}
  (\varepsilon-v_0)u_1 + \frac{\Delta^2-\varepsilon(\varepsilon-v_0)}
  {\sqrt{|\Delta^2-{\varepsilon}^2|}}
  \nonumber \\ &&
  {} \times \Big[ u_2\theta(\Delta-|\varepsilon|) 
  -u_3\theta(|\varepsilon|-\Delta) \Big] - D \bigg\} 
  \nonumber \\ && {}-\frac{iS\ {\rm sgn}(\varepsilon)}
  {4v_{\rm F}^2\sqrt{\varepsilon^2-\Delta^2}} \bigg[
  \theta\bigg(\Delta<|\varepsilon|<\sqrt{\frac{\omega_0^2}{4}+\Delta^2}\bigg)
  \nonumber \\ && {}\times
  \Big(\varepsilon+\sqrt{\varepsilon^2-\Delta^2}\Big)v_{-} 
  + \Big(\varepsilon-\sqrt{\varepsilon^2-\Delta^2}\Big)v_{+} 
  \nonumber \\ && {}\times
  \theta\bigg(\Delta<|\varepsilon|<\sqrt{(D-\frac{\omega_0}{2})^2+\Delta^2}
  \bigg) \bigg],
\label{gsmall_22}
\end{eqnarray}
\begin{eqnarray}
  g^{\tau\tau,r}_{12}(\varepsilon)&=&g^{\tau\tau,r}_{21}(\varepsilon)
  =-\frac{S\Delta}{4\pi v_{\rm F}^2}
  \bigg\{ \frac{\omega_0}{\sqrt{|\Delta^2-{\varepsilon}^2|}} 
  \Big[ u_2\theta(\Delta-|\varepsilon|) 
  \nonumber \\ && {}
  -u_3\theta(|\varepsilon|-\Delta) \Big] + u_1 \bigg\} 
   -\frac{iS\ {\rm sgn}(\varepsilon)\Delta}
  {4v_{\rm F}^2\sqrt{\varepsilon^2-\Delta^2}} 
  \nonumber \\ && {} \times \Bigg[ v_{+}\ \theta\left(\Delta<|\varepsilon|
    <\sqrt{(D-\frac{\omega_0}{2})^2+\Delta^2} \right)
  \nonumber \\ && {} + v_{-}\ \theta\left(\Delta<|\varepsilon|<
      \sqrt{\frac{\omega_0^2}{4}+\Delta^2}\right) \Bigg],
  \label{gsmall_12}
\end{eqnarray}
with $S$ representing the area of sample.
$v_{\nu}=\frac{\omega_0}{2}+\nu\sqrt{|\varepsilon^2-\Delta^2|}$
with $\nu=0,\ \pm$ and
\begin{eqnarray}
  && \hspace{-1.cm}
  u_1={\rm ln}\left|\frac{(D-\frac{\omega_0}{2})^2
      +\Delta^2-{\varepsilon}^2}
    {\frac{\omega_0^2}{4}+\Delta^2-{\varepsilon}^2}\right|, \\
  && \hspace{-1.cm}
  u_2={\rm arctan}\frac{2D-\omega_0}
  {2\sqrt{\Delta^2-{\varepsilon}^2}}
  + {\rm arctan}\frac{\omega_0}
  {2\sqrt{\Delta^2-{\varepsilon}^2}}, \\
  && \hspace{-1.cm}
  u_3=\frac{1}{2}{\rm ln}\left|\frac
     {(D-\frac{\omega_0}{2}+\sqrt{\varepsilon^2-\Delta^2})
       (\frac{\omega_0}{2}+\sqrt{\varepsilon^2-\Delta^2})}
     {(D-\frac{\omega_0}{2}-\sqrt{\varepsilon^2-\Delta^2})
       (\frac{\omega_0}{2}-\sqrt{\varepsilon^2-\Delta^2})} \right|.
\end{eqnarray}

Comparing the above formulae with Eqs.~(A.3)-(A.8) in
Ref.~\onlinecite{Berakdar_12}, one observes
a lot of differences between them.
Specifically, in $g^{\tau\tau,r}_{11}$ and $g^{\tau\tau,r}_{22}$, they missed
$\omega_0/2$ in all $v_\nu$ terms.
The errors in the interband term are even more serious: their formula can be
written as
\begin{eqnarray}
  g^{\tau\tau,r}_{12}(\varepsilon)&=&g^{\tau\tau,r}_{21}(\varepsilon)
  =-\frac{S\Delta}{4\pi v_{\rm F}^2}
  \bigg[ u_1 \theta(\Delta-|\varepsilon|)
  \nonumber \\ && {
 +\frac{u_3}{\sqrt{|\Delta^2-{\varepsilon}^2|}} 
  \theta(|\varepsilon|-\Delta) \Big]  
   -\frac{iS\ {\rm sgn}(\varepsilon)\Delta}
  {4v_{\rm F}^2\sqrt{\varepsilon^2-\Delta^2}} 
  \nonumber \\ && {} \times
  \Bigg[\theta\bigg( \Delta<|\varepsilon|
  <\sqrt{(D-\frac{\omega_0}{2})^2+\Delta^2} \bigg)
  \nonumber \\ && {}
  - \theta\Bigg( \Delta<|\varepsilon|
  <\sqrt{\frac{\omega_0^2}{4}+\Delta^2} \Bigg)\Bigg].
  \label{gsmall_wrong}
\end{eqnarray}
It is seen that all terms in the above equation for
$|\varepsilon|>\Delta$ take the incorrect dimension.
Since there are so many errors in their formulae, 
these errors are not likely to only come from typos.


In order to show their errors more clearly, we calculate the DOS of the isolated
graphene in the field-free case, which has the form
\begin{eqnarray}
  && \hspace{-0.8cm}
  D^{\tau}_0(\varepsilon)= -\frac{1}{\pi}{\rm Im} 
  \sum_{s}\tilde{g}^{\tau,r}_{00,ss}(\varepsilon)\bigg|_{\Delta=0}
  \nonumber \\
  && \hspace{0.2cm} = -\frac{1}{\pi}{\rm Im} 
  \sum_{s}{g}^{\tau,r}_{00,ss}\Big(\varepsilon-s\frac{\omega_0}{2}\Big)\bigg|_{\Delta=0}.
\end{eqnarray}
From our formulae, one obtains
\begin{equation}
  D^{\tau}_0(\varepsilon)=\frac{S}{2\pi v_{\rm F}^2}\;|\varepsilon|\;
  \theta(|\varepsilon|<D),
\end{equation}
which is exactly the well-known formula 
of the DOS in graphene\cite{Neto_rev_09}
and is linear with energy. In contrast, from the formulae by Ding {\em et al.},
 one obtains
\begin{equation}
  D^{\tau}_0(\varepsilon)=\frac{S}{2\pi v_{\rm F}^2}
  \Big(|\varepsilon|-\frac{\omega_0}{2}\Big)
  \theta(|\varepsilon|<D).
\end{equation}
Obviously, the above formula is incorrect, especially it gives the negative DOS when 
$|\varepsilon|<\frac{\omega_0}{2}$.

\end{appendix}

