Klein-tunneling transistor with ballistic graphene

Today the availability of high mobility graphene up to room temperature makes ballistic transport in nanodevices achievable. In particular, p-n-p transistor in the ballistic regime gives access to the Klein tunneling physics and allows the realization of devices exploiting the optics-like behavior of Dirac Fermions (DF) as in the Vesalego lens or the Fabry P\'erot cavity. Here we propose a Klein tunneling transistor based on geometrical optics of DF. We consider the case of a prismatic active region delimited by a triangular gate, where total internal reflection may occur, which leads to the tunable suppression of the transistor transmission. We calculate the transmission and the current by means of scattering theory and the finite bias properties using Non Equilibrium Green's Function(NEGF) simulation.


I. Introduction
As a high mobility material, graphene is well suited to high frequency electronics [1,2]. Additionally, thanks to the weakness of electron-phonon coupling [3][4][5][6], the high mobility persists at high temperature [7] so that one can envision ballistic graphene electronics in realistic operating conditions. The absence of a band-gap restricts the switching off capability of graphene, the conductivity reaching a shallow minimum at charge neutrality limited by the quantum of conductance e h 4 2 . For logic electronic applications a band-gap can be restored using e.g. geometrical confinement [8], but this is usually at the cost of reduced mobility. In microwave electronics both a high mobility and a significant switching capability are needed to achieve large voltage and power gains at high frequency. To this end we explore here the new possibilities offered by Klein tunneling (KT), a hallmark of Dirac fermions (DFs).
We consider the case of p-n-p transistors where the resistance at the source-gate and gate-drain junctions is dominated by the Klein tunneling effect, while DFs propagate ballistically in the barrier [9]. With graphene being ambipolar, the generalization to n-p-n transistors is straightforward. This is the regime of DF optics where refraction and transmission at the interfaces are determined by Fresnel-like relations [10] and depend on the abruptness of junctions [11,12]. Architectures exploiting the optical properties of DFs have been already proposed, such as the Vaselago lens [11,13], and demonstrated, such as Fabry-Pérot interferometers [14][15][16] or tilted p-n junctions [17,18] (for a review, see reference [19]). The possibility of a 'latticetronics' of Klein tunneling currents in an armchair ribbon using potential barriers has also been recently discussed [20]. Here we consider a different geometry, that exploits total internal reflection in a Klein tunneling prism. The prism is made of an n-doped ballistic triangular domain embedded in a p-doped diffusive area, the latter being controlled electrostatically or chemically. As in light reflectors, an array of KT-prisms can be used to form the active channel of a Klein tunneling transistor; this geometry minimizes the gate length to keep ballistic transport conditions in a way similar to Fresnel lenses which minimize glass weight and light absorption [21]. In this work we calculate the low-energy transmission of such a device using scattering theory, i.e. within an intuitive and physically transparent approach. Our study is supplemented by atomistic simulations using the non-equilibrium Greenʼs function (NEGF) formalism that includes a more complete description of quantum transport. It gives access to finite bias properties while accounting for short gate effects like direct drain-source tunneling, diffraction, dispersion and finite temperature effects. Both approaches predict for the Klein tunneling transistor a strong suppression of conductance that can eventually go below the minimum conductance at charge neutrality. This transistor can be used as a tunable barrier for electrostatic quantum confinement to achieve, e.g. single Dirac fermion pumps working at low temperature. It is also suited for microwave electronics as it cumulates significant resistance in the OFF state with a large conductance in the ON state. A nanoscale variant has been recently proposed that predicts very large ON-OFF ratios for logic applications [22]. Our approach is different: it is based on geometrical DF optics following reference [23]; it is more conservative in accounting for diffusive transport in the leads and targets microwave electronic applications [24,25].
The principle of the Klein tunneling prism is sketched in figure 1. It can be understood by scattering theory and relies on the total internal reflection in a triangular n-doped graphene region (concentration n, angle α), the refractive medium, embedded in a p-doped environment (concentration ≲ p n) playing the role of vacuum. The refraction at the input p-n junction obeys a Snell-Descartes-like relation ϕ ν θ = sin sin 1 1 where ν = − n p / is the (negative) refraction index [10] as sketched in the figure. The transmission ϕ ( )  can be calculated for a sharp junction [11], but also for smooth junctions (see below and reference [12]). The latter case is more realistic and suitable for device modeling. The main feature of Klein tunneling is the enhanced forward scattering for ≫ n p. As shown in the figure, the refracted DF beam is focussed along the junction normal within an angular opening θ θ ⩽ c 1 , where θ ν = arcsin (1/ ) c is the critical angle. At the drain side, which is an n-p junction where DFs are impinging at incidence ϕ α θ = − The reflection can be controlled and transmission restored, on decreasing n-doping, i.e. increasing θ c , as shown in figure 1(b). Note that one obtains an equivalent effect on decreasing the prism angle down to α = 0; it means that the transmission in the open state ≃ n p (or ν = −1) should approach that of a rectangular gate transistor. Taking α = 45 deg, one can estimate the gate doping for full reflection: > n p 6.8 (or ν > 2.6). The transmission of a KT transistor is suppressed here when increasing the gate doping deep in the metallic regime; this is a marked difference with conventional semiconducting transistors where transmission is suppressed on pinching-off the channel. Operating a transistor with a large intrinsic transmission should benefit the dynamical and noise properties of the device. In order to establish the above sketch, a full calculation is needed that accounts for finite transmission and multiple reflection effects.
We rely on the realistic modeling of the p-n junctions. Translational invariance along the junction parallel to the y-axis yields the above-mentioned Snell-Descartes relation between θ and

ϕ. A smooth junction, characterized by a length
we rely on the expression by Cayssol et al reference [12], which is the exact solution for a potential step described by a Fermi function, figure 2(b)). We use = d w 4.5 to make contact with the smooth junction formula. The transmission is given by , which is the regime of total internal reflection. The contrast is close to that predicted for a sharp junction but different from that of a smooth junction (not shown) which justifies the use of the Cayssol interpolation formula.
The principle of the Klein tunneling prism can be realized in a Klein tunneling transistor (figure 2(c)) on stacking, head-to-tail, a series of prisms to realize a sawtooth gate with symmetric elementary units made of an isosceles triangle of opening angle π−2α. This geometry allows to implement the refractor principle while keeping the gate length short enough to remain in the ballistic regime. The n-doped barrier is assumed to be ballistic and controlled by the transistor gate. The leads can be either access regions or simply metallic contacts. In the latter case, the carrier concentration is set by electrochemical doping, in the former it can be tuned by a second electrostatic gate. To assess device properties, we calculate the transmission of DF rays across an elementary triangle. We consider a beam with a given incidence angle ϕ 1 and position y along the source side. The trajectory of the ray into the prism is calculated according to above p-n junction refraction and reflection rules iterated up to twenty internal reflections. The beam intensity inside the prism decreases at each reflection, the transmitted intensity outside the prism increases accordingly. We then integrate over y positions and all incident angles (diffusive leads) weighted by the angular density of states to get the overall transmission 〈 〉 shows that of the prism angle while keeping the gate area constant. In standard Klein tunneling with rectangular gates, the transmission saturates at 75% (resp. 45%) for sharp (resp. 10 nm-long) junctions which determines the ON-state current and conductance μ ≃ − g W / 5 mS m ON 1 (10 nm-long junction). As seen in panel (b), the OFF-state current is very sensitive to the prism opening angle; it is zero for α = 45 deg (and ≳ n p 6 ) but increases rapidly whenever α deviates by more than 5 deg from this value (see inset of figure 3(b)). From this analysis, we conclude that the optimal geometry for a KT transistor is α = ± 45 5 deg and

, of voltage gain
and of transit frequency GHz V 1 (for a gate length = L 100 nm). We discuss now the physical limits of our geometrical optics description which relies on the assumption that ≲ ≪ The refraction principle is scale independent but in practice the characteristic length of the gate L should be smaller than the ballistic length, while remaining larger than the electronic wave length to avoid deleterious effects of diffraction. The principle of the reflector being very sensitive to refraction angles (see the inset of figure 2), an in angular orientation, due to finite wavelength, may drastically affect the total internal reflection. Another limitation for the OFF state resistance is the direct drain source tunneling which is also strongly scale dependent. Finally we have neglected for simplicity the dispersion effects arising from broad sources such as finite bias DF emission. In order to quantify these effects we have performed a numerical simulation of the KT transistor by means of NEGF method. Our model is based on a tight-binding Hamiltonian to describe the electron state in the graphene honeycomb lattice [26].
NEGF simulation assumes a transistor of ballistic graphene whose size is scaled down by a factor 2.5 with respect to a realistic device to make the calculation achievable with available computational resources. Additionally, the self-consistent solution of Poissonʼs equation has been deactivated here to save CPU time. The gate-induced potential in graphene was modelled as a 'square' barrier of height reduced by V 2 ds at the source-end and enhanced by V 2 ds at the drain-end. This simple description of potential has been shown realistic from self-consistent simulation of usual graphene field effect transistors (GFETs) in the p-n-p or n-p-n regimes [27]. The transistor gate is taken either as rectangular or triangular but with the same gate area. The triangular gate is an elementary unit of figure 2(c). The channel width is 80 nm; the gate length is 40 nm for the rectangular device. It is modulated between 20 nm and 60 nm for the triangular device. We set the source side doping to = × − p 2.3 10 cm 12 2 which corresponds to a Fermi energy T is displayed for the rectangular geometry (red lines) and the triangular geometry (blue lines) for the ON state of the transistor (dotted lines) and the OFF state (solid lines). In the ON state (small gate doping) two minima of transmission appear corresponding to the energy position of the CNP in the lead and in the barrier. A transmission maximum is observed that lies mid-way between both CNPs. As ∼ n p 2 , the maximum coincides with the Fermi energy, irrespective of the gate geometry. The transmission is large as expected for an ON state. In the OFF state (large gate doping) the transmission maximum is shifted toward higher energy and scaled up. The transmission at the Fermi level becomes sensitive to the gate geometry: it is large for the rectangular gate and small for the triangular one. Oscillation of ε ( ) T observed for rectangular gates reflect Fabry-Pérot interferences; they are blurred for triangular gates. This effect has been neglected in the geometrical optics approach.
From the transmission, one calculates the current flowing through the device on integrating transmission in the Fermi energy interval between source and drain. Finally we plot in  figure 5(b)). To compare with the scattering theory results we also plot in 5(a) the sharp junction rectangular (resp. triangular) device transmission (red dotted line, resp. blue dotted line) from figure 3(a).
First of all we notice that NEGF simulation for the rectangular device is in close agreement with scattering theory. Regarding the triangular device, the agreement is more qualitative. The maximum current in the ON state is smaller and the OFF state current larger than geometrical optics predict. We assign the difference to the diffraction effect that is more pronounced in As an illustration, we have plotted in 5(a), the current for a triangular device with an opening angle of α 1.15 that compares favorably with the numerical calculation. The high bias response in figure 5(b) shows a comparable behavior. In both cases, numerical simulations show a current suppression below the CNP minimum which is the signature of the Klein tunneling transistor. Our models correspond to two extreme cases of a macroscopic Klein tunneling transistor on one side and a nanometric device on the other one. An actual Klein tunneling transistor, with a gate length ≲ L 100 nm limited by the elastic mean free path, should have transport properties bracketed by these two limits. From a technical point of view, abrupt junctions can be realized using local back gates, split by a nanometer thin gap and insulated from graphene by atomically thin hexagonal boron nitride layers [7].
In conclusion, we have introduced a Klein tunneling transistor architecture that takes advantage of anomalous refraction properties of Dirac fermions in graphene to realize a tunable electrostatic barrier for Dirac fermions. We have used a geometrical optics model to explain the principle of the transistor. We have performed extensive numerical simulations of nanoscale variants of the device that confirm the transistorʼs effect while taking a full account of finite size effects, in particular Dirac fermion diffraction. Our modeling will prove useful for the design of actual devices and to evaluate their potential in terms of single electron pumps and microwave electronics.