\section{Introduction} High-quality graphene nanostructures are the most likely building blocks in future electronics, plasmonics and photonics assemblies.
Most of their striking features arise from the strictly two-dimensional, honeycomb lattice structure of the basic material, which causes the nontrivial topology of the electronic wave function and an almost linear low-energy spectrum of the chiral quasiparticles (charge carriers) near the so-called Dirac nodal points~\cite{CGPNG09}.  From an application-technological point of view, the tunability of the transport properties of graphene by external electric and magnetic fields is of particular importance~\cite{Go11}. This allows a controlled modification of selected areas of the sample by gating.  For instance, applying nanoscale top gates on  graphene nanoribbons (GNRs) or bilayer graphene, single or double quantum dots  with particular shape have been produced in recent experiments~\cite{LHV10,AMY12,MKHWPS14}.  

For a large circular dots, the refraction at the boundary  gives rise to two coalescing 
caustics that focus the electron density in the disk-shaped  region~\cite{CPP07,AU13}.  Resonances in the conductance~\cite{BTB09}  and the scattering cross section~\cite{HBF13a}  indicate electron confinement for small quantum dots  in monolayer graphene as well. Apparently the difficulty to localise Dirac electrons by electrostatic potential barriers (which is related to the Klein tunnelling phenomenon) could be overcome: Because of the circular symmetry,  the angle of incidence at the dot's boundary is a constant of motion so that electrons with nonzero angular momentum are confined~\cite{BTB09}. Forward scattering and Klein tunnelling will also be suppressed 
if a Fano resonance arises between the background partition and the resonant contribution to electron scattering~\cite{HBF13a}. These results, obtained within continuum Dirac theory, were numerically confirmed for a lattice model~\cite{PAS11,PHF13}.  The role of resonances have been analysed in various graphene scattering experiments~\cite{HZZ09,MOWBFGGBFM11}.  
Note that the scattering (transport and optics) of massless Dirac fermions on potential barriers substantially differs from that of massive chiral Dirac fermions or Schr\"odinger electrons~\cite{AFkt11}. 



\section{Model and numerical methods} In this work,  we address the scattering and (linear) transport in armchair  GNRs with gate-defined quantum dots in the end-contacted lead-sample geometry most relevant for experiments~\cite{KR12,RG13,LRG13}. The leads can be viewed as semi-infinite quantum wires; the contacts are metallic or gated graphene. 

We describe the sample by the tight-binding Hamiltonian 
\begin{equation}
     {H} =  \sum_{i}V_i^{} {c}_i^{\dag} {c}_i^{} 
           -t \sum_{\langle ij \rangle}({c}_i^{\dag} {c}_j^{} + \text{H.c.})\,,
\label{eq:H}
\end{equation}
where $c_i^{(\dag)}$ annihilates (creates) an electron  at Wannier site $i$ of a honeycomb lattice with carbon-carbon distance $a=1.42$~\AA. Here the hopping amplitude between nearest neighbour sites $\langle ij\rangle$ is given by $t$ ($\simeq$ 3eV) and the site-dependent 
 potentials $V_i=V \sum_{n} \Theta (R-|\vec{r}_i-\vec{r}_{n}|)$ 
model electrostatically defined quantum dots with radius $R$  centred at positions $\vec{r}_n$.

In the limit of vanishing bias voltage, the conductance between the left (L) and right (R) leads can be obtained within the limits of the Landauer-B\"uttiker approach: 
\begin{equation}
G_{\rm LR}= \frac{\mathrm{e}^2}{\mathrm{h}} \sum_{n \in L, m \in R} | S_{n,m} |^2\,,
\label{eq:cond}
\end{equation}
where  $S_{n,m}$ is the scattering matrix between all open (i.e., active) lead channels~\cite{Dat95}. We furthermore exploit the local density 
of states (LDOS) at a certain site $i$ of the sample, 
\begin{equation}
\text{LDOS} (E) = \sum_{l} |\langle i | l\rangle |^2 \delta (E-E_l)\,,
\label{eq:ldos}
\end{equation}
as a probe of a possible electrostatic confinement
of the charge carrier~\cite{PHWF14,SB14}.
In Eq.~(\ref{eq:ldos}), $E$ is the particle's energy
and the sum extends over all single-electron eigenstates
$|l\rangle=c_l^\dagger |0\rangle$ of $H$ with energy $E_l$. We use the kernel polynomial method~\cite{WWAF06} and the `Kwant' software package~\cite{Kwant14}
for evaluating the conductance and LDOS numerically. 
\begin{figure}[t]
\includegraphics[width=.91\linewidth]{Q.pdf}\\
\includegraphics[width=.95\linewidth]{G_single_dot.pdf}\\
\includegraphics[width=\linewidth]{singledot_ldos.pdf}
\caption{Total and partial scattering efficiencies  (top), conductance (middle) and LDOS (bottom) of a single quantum dot 
(the black circle  specifies the dimension of the dot).   For further explanation see text.}
\label{fig1}
 \end{figure} 
\section{Results and discussion}
\subsection{Single quantum dot}
We first consider a single dot, see Fig.~\ref{fig1}.
To begin with we determine  the scattering efficiency $Q$, which is 
the scattering cross section divided by the geometric cross section,
of a plane Dirac electron wave hitting on a dot embedded 
in an infinite graphene sheet, within the Dirac approximation
(for details of such a  calculation cf. 
Ref.~\cite{HBF13a}). Since the energy is conserved in our setup 
(i.e., there is no coupling to bath degrees of freedom etc.), the scattering is always elastic (coherent). 
The data displayed in the upper panel for $E=0.05 \,t$ and $R=5$~nm
shows that a series of resonances appears when the gate voltage is varied (besides the rather broad $a_0$ and $a_1$ modes, the 
$a_{m>1}$ resonances are very sharp),
which can be attributed to the excitations of higher dot normal modes~\cite{HBF13a}. 
So each maximum marks a combination of $E$, $R$, and $V$,
where a certain normal mode fits particularly well into the dot,
while all other  modes are suppressed.
The middle panel indicates that these signatures also show up in the conductance,
which this time, however, was computed for a 20~nm wide armchair GNR placed between leads.
That is why we can assign the signal at $V\simeq - 0.09 t$ to the $a_1$ mode.
The value for $G[V=0]$  tells us that the system exhibits three open transport channels.
The lowermost panel presents the LDOS close to the $a_1$ resonance marked
in the $G$~vs.~$V$ plot by the blue diamond.
Obviously the electron is almost completely localised at the quantum dot,
where the superposition of the states with $j=\pm(m+1/2)$ leads to
the characteristic vortex pattern of the $a_1$ mode. There are six vortices close to the boundary of the dot which are reflected by three preferred 
scattering directions in the far field~\cite{HBF13a}.  
\begin{figure}[t]
\centering
\includegraphics[width=0.8\linewidth]{G_chain_dot.pdf}
\includegraphics[width=\linewidth]{chain_ldos.pdf}
\caption{Conductance (top) and LDOS (lower panels) of a GNR with a linear array of quantum dots. The colour code is the same as in Fig.~\ref{fig1}.}
\label{fig2}
 \end{figure}  
 
 
\subsection{Qunatum dot chain}
We now study quantum dots arranged like a string of pearls.
Figure~\ref{fig2}  presents the conductance and the LDOS 
for a contacted GNR sample with five dots at intervals of 20~nm. 
Again we set $R=5$~nm, $E=0.05 t$, and the width of the armchair GNR is 20~nm. 
Note that the conductance undergoes a dramatic change 
if $V$ is slightly varied near the resonance points (see red and blue diamonds).
In this way the system may act as a switch.
The results for the LDOS provide further insight into
how the system behaves locally for  $V=-0.088\,t$ ($V=-0.080\,t$),
i.e., at the local minimum (maximum) of the conductance.
At $V=-0.088\,t$, $E=0.05 t$ (middle panel), the quantum dots are off-resonance from an eigenmode
of the particle: As a result the scattering is enhanced and the conductance is 
low (note that the figure shows a stationary situation). 
For an  applied bias $V=-0.080\,t$ (lower panel) all dots largely support the $a_1$ mode.
Since the wave functions of quasi-bound $a_1$ modes at neighbouring dots
overlap resonantly, electrons are no longer permanently confined 
to a particular dot but will be transferred through the GNR 
by an effective inter-dot hopping process. Remarkably, in this case the conductance is almost as large as 
for the GNR without quantum dots ($G[V=0]$) where the LDOS is the same (i.e., uniform) everywhere between the leads.
\begin{figure}[t]
\centering
\includegraphics[width=0.8\linewidth]{G_array_dot.pdf}
\includegraphics[width=\linewidth]{array_ldos.pdf}
\caption{Conductance (top) and LDOS (middle/bottom) of an armchair GNR with a square dot superlattice.}
\label{fig3}
 \end{figure} 
\subsection{Quantum dot matrix}
Figure~\ref{fig3} illustrates that the same behaviour 
could be achieved for a square quantum dot superlattice 
made on an armchair GNR (now the width of the GNR is 60~nm,
the other parameters are the same as in Fig.~\ref{fig2}). 
For $V=0$ the system has nine open transport channels. 
Again the relatively large conductance at $V=-0.080\,t$ 
(compared to $G$ for $V=-0.065\,t$) originates 
from an effective (coherent) inter-dot  transfer of the electron.
This hopping process is expected to take place on a reduced time scale~\cite{PHF13}.
As can be seen from the lowermost panel, at $V=-0.065\,t$,
the LDOS notably shrinks if we move to the centre of the GNR;
this depletion of the local particle density 
reduces the inter-dot hopping and consequently the conductance.     
\begin{figure}[t]
\centering
\includegraphics[width=\linewidth]{minimal_ldos.pdf}
\caption{LDOS for a linear (top) and a square (bottom)  quantum dot array on a GNR with almost vanishing conductance.}
\label{fig4}
 \end{figure} 
 
While graphene has a gapless band structure and also the armchair GNRs
considered in Figs.~\ref{fig1}-\ref{fig3} have  $G >0$  for all $E$
at $V=0$  (i.e., there are always open transport channels),
the conductance of the GNRs with imprinted quantum dots may
tend to zero in certain ranges of $V$; see the black squares in Figs.~\ref{fig2} and~\ref{fig3}.
Apparently the linear and square lattice periodic arrays of gate-defined
quantum dots give rise to new superlattice band structures,
which exhibit energy (pseudo) gaps  just as for ordinary solids~\cite{LWSYL09}. 
This is reflected in the LDOS depicted in Fig.~\ref{fig4}.

\section{Conclusions} To summarise, using exact numerics,
we have studied the electronic and transport properties 
of contacted graphene nanoribbons with gate-defined quantum dots of realistic size (i.e., 
with a few thousand carbon atoms that feel the electrostatic potential step). 
The conductivity of such systems is very sensitive to the model parameters and can even vanish. 
This provides an opportunity to use them as switches by varying the top gate potential.
A detailed analysis of the local density of states shows 
that an electrostatic electron confinement  is possible for single quantum dots
in the quantum regime due to resonant scattering. 
Thereby the incident electron is fed into vortices, which trap the particle. That means for small dots the 
particle confinement is not by total internal reflection.  In addition, for periodic quantum dot arrays a  superlattice
band structure evolves that allows for an effective inter-dot  
transport on a reduced energy scale. The phenomena detected in this work 
might be exploited in actual graphene-based devices. From a theoretical point of view, it seems worthwhile 
to study a modulation of the quantum dot barriers---via the top gates---by oscillating fields, 
which induces both inelastic tunnelling and resonant transport through excited states\cite{SHF15b}.
Another issue is the electrical and magneto transport trough quantums dots on folded graphene samples
that exhibit superlattice structures in itself~\cite{SRSH14}. We finally like to emphasize that such graphene quantum 
dots might used to construct logic gates and ``quantum computing'' devices~\cite{WL11}, which would be also 
an interesting topic for future studies. 

\begin{acknowledgement}
This work was funded by the Deutsche Forschungsgemeinschaft through 
Priority Program SPP1459 `Graphene' and the Competence Network 
for Scientific High-Performance Computing 
in Bavaria (KONWIHR III, project  PVSC-TM). 
\end{acknowledgement}

\bibliographystyle{pss}
