# Electron-hole asymmetry in two-terminal graphene devices

## Abstract

A theoretical model is proposed to describe asymmetric gate-voltage dependence of conductance and noise in two-terminal ballistic graphene devices. The model is analyzed independently within the self-consistent Hartree and Thomas-Fermi approximations. Our results justify the prominent role of metal contacts in recent experiments with suspended graphene flakes. The contact-induced electrostatic potentials in graphene demonstrate a power-law decay with the exponent varying from to . Within our model we explain electron-hole asymmetry and strong Fabri-Perot oscillations of the conductance and noise at positive doping, which were observed in many experiments with submicrometer samples. Limitations of the Thomas-Fermi approximation in a vicinity of the Dirac point are discussed.

###### pacs:

72.80.Vp,73.23.Ad,73.40.Cg## I Introduction

The recent observation of the fractional quantum Hall effect in grapheneDu:2009 () calls for better understanding of two-terminal transport measurements in suspended graphene devices. The two-terminal geometry has also been used to reach the ballistic regime of transport with submicrometer flakes on SiO substrate.Heersche:2007 (); Danneau:2008 () These and other experimentsDu:2008 (); Du:2008:NN () share such common features as the electron-hole asymmetry in the gate-voltage dependence of the conductance and the prominent conductance oscillations at positive doping that are beyond a simplistic theory of ballistic transport.Tworzydlo:2006 (); Titov:2006 () In this paper we attribute these experimental observations to the effect of charge transfer in a vicinity of metal electrodes and to a weak screening in graphene at low doping.

The influence of metal contacts on electron transport in nearly ballistic graphene has been investigated both experimentallyLee:2008 (); Huard:2008 (); Blake:2009 (); Mueller:2009 (); Russo:2010 () and theoretically.Khomyakov:2009 (); Khomyakov:2010 (); Vanin:2010 (); Zhang:2008 (); Datta:2009 (); BarrazaLopez:2010 (); Nouchi:2010 () An agreement between theory and experiment has, however, been limited due to both the complexity of experimental setups and a lack of important ingredients in the considered theoretical models. In this study we suggest a minimal model, which leads to a very good description of experimental data by Xu Du et al.Du:2008:NN () on charge transport in suspended graphene. The gate-voltage dependence of transport properties predicted by the model also agree qualitatively with those observed by Heersche et al.Heersche:2007 () for graphene on SiO substrate.

The electrostatic potential landscape in graphene with metal contacts has been imaged by means of scanning photocurrent microscopyLee:2008 (); Mueller:2009 () to demonstrate the long-range decay of the contact-induced potential (about 1 m). These experiments and the theoryNouchi:2010 () suggest that the charge density in graphene covered by a metal is pinned, i.e. not affected by the gate voltage.

The pinning is much weaker if graphene sheet is oxidized and the metal-graphene contact is resistive. In this case the second minimum in the gate voltage dependence of conductance can emerge which corresponds to the charge density minimum in the metal-covered graphene. In Ref. Nouchi:2010, this mechanism is claimed to be responsible for the electron-hole asymmetry observed in Refs. Nouchi:2008, ; Du:2008, . We argue, however, that the asymmetry seen in Refs. Du:2008, ; Du:2008:NN, ; Du:2009, is entirely due to the contact-induced potential in the free-standing graphene. We also regard the observed saturation of the conductance for large negative gate voltages as the signature of the charge density pinning in the metal-covered graphene. The same mechanism of the electron-hole asymmetry has been put forward recently in Ref. Nouchi:2011, .

The adsorption of graphene on metal substrates has been studied by means of the density functional theory (DFT) in Refs. Khomyakov:2009, ; Vanin:2010, . These studies suggest that the main effect of metal deposition can be modeled by a shift of the Fermi level, i.e. by the chemical doping. This shift can be of any sign. Its precise value, , is comparatively large due to the small density of states of intrinsic graphene, and depends on the workfunctions of graphene and the metal as well as on the direct chemical interaction. The most recent analysisVanin:2010 () predicts, e.g. eV for Al, eV for Ag, eV for Cu, and eV for Au.

In this paper we refer to the metal-covered graphene as the lead, and to the free-standing graphene as the sample. The scalar step-function potential in the Dirac equation provides the simplest model of the sample-lead interface.Tworzydlo:2006 (); Titov:2006 () Remarkably, the exact analytical solution can also be found if the interface is modeled by an exponentially decaying potential.Cayssol:2009 () These models, however, are not self-consistent and cannot describe the large electron-hole asymmetry observed in experiments.

Microscopic DFT-modeling of the sample-lead interface has been undertaken in Ref. BarrazaLopez:2010, to investigate the contact-induced potential in short (nm) graphene samples, and its effect on the two-terminal conductance. In this study a gate electrode is absent and the conductance is studied as a function of Fermi energy using electrostatic potential for undoped graphene. The reported electron-hole asymmetry is, therefore, of a different type than the one occurring with pinned charge density in the leads. It is also hard to scale the results of Ref. BarrazaLopez:2010, to realistic system sizes ( nm).

Below we propose an effective model, which takes into account both the charge density pinning in the leads and the self-consistent contact potential in the sample. In formulating the model we focus on the two-terminal geometry and neglect the resistance of the sample-lead interfaces. We employ the effective Dirac Hamiltonian and take into account the large density of states in a metal. The gate electrode is incorporated into the model, while the system size poses no serious restriction. The presence or absence of a substrate is taken into account by an appropriate choice of the permittivity constant. We calculate the self-consistent potential in graphene using two different approximations, namely the self-consistent Hartree (SCH) and the Thomas-Fermi (TF) approximation. From the Kubo formula we, then, calculate the two-terminal conductance and noise as a function of gate voltage.

We find that the different approximations result in a qualitatively similar form of the effective potential. The potential penetrates deeply into the sample and reveals a power-law decay with the exponent varying from to depending on the electron concentration. This behavior agrees with the TF analysis of Ref. Khomyakov:2010, . We demonstrate that the slow decay of the potential is responsible for the asymmetric gate-voltage dependence of conductance and noise, increasing in intensity with the lead doping . For small positive gate voltages the Dirac point crosses the chemical potential twice inside the sample so that two n-p interfaces are formed. The slow potential decay ensures that the transmission through the interfaces is strongly selective with respect to the momentum direction.Cheianov:2006 () The electron scattering at the n-p interfaces is the reason for the enhanced Fabry-Pérot oscillations of the conductance as a function of electron concentration at moderate positive doping. We note that the oscillations transform to the pronounced resonances in the ribbon geometry .Darancet:2009 () When describing gate-voltage asymmetries we refer to negatively doped leads, . For positively doped leads the polarity of the effect is reversed. One can further argue that Fabry-Pérot oscillations are easily distorted by disorder. The presence and shape of the Fabry-Pérot oscillations can, therefore, be used as a direct indicator of the sample quality.

## Ii Charge and potential profiles

The charge carriers in the graphene sheet (defining the -plane) are described by the Dirac Hamiltonian

(1) |

where is a vector of Pauli matrices and is the spatially dependent effective potential. In Eq. (1) we neglect any effects arising from the finite band width and assume translational invariance in .

Doping of the graphene sheet is caused by the contact with the metal electrodes and by applying a voltage to the global back gate. We model the lead by graphene with the chemical potential measured from the Dirac point as shown in Fig. 1(a).

We use the externally compensated charge density in graphene as an input parameter for a self-consistent calculation. We focus on a two-terminal setup for which the metal-covered and free-standing parts of the graphene sheet correspond to (leads) and (sample), respectively. Accordingly, the compensated charge density has the spatial dependence

(2) |

where is proportional to the gate voltage, since the distance to the gate electrode is typically such that the quantum corrections to capacitance can be neglected. We assume that the charge density in the lead, , is pinned,Nouchi:2010 () i.e. not influenced by the gate-voltage. The screening inside the gate electrode can be neglected for , with the relative permittivity of the gate dielectric. Finally, the inter-valley coupling at the metal-graphene interface is also disregarded.

The relation between the charge densities and and their respective chemical potentials and (cf. Fig. 1(a)) is given by

(3) |

where is the Fermi distribution function, is the Heaviside step function, and

(4) |

is the density of states (DOS) in infinitely extended ballistic graphene described by Hamiltonian (1) including the spin and valley degeneracy. At zero-temperature Eq. (5) becomes

(5) |

Screening in the leads is strongly enhanced by the metal electrodes deposited on top of the graphene sheet. This effect, which ensures the charge-density pinning, is taken into account by means of an additional DOS , which is assumed to be energy independent. The transport properties do not depend on the value of under these conditions.

The density profile corresponds to a charge-neutral setup with disconnected leads. Bringing the sample into electric contact with the leads causes the charge redistribution and the band bending shown in Fig. 1(b). The excess charge density in the entire device is defined as

(6) |

where is the total charge density in the graphene sheet relative to intrinsic graphene, and

(7) |

is the additional charge density induced in the metal.

The potential profile and the total charge density have to be determined self-consistently. In this study we ignore the exchange interaction and restrict ourselves to the scalar Hartree potential

(8) |

where

(9) |

The coupling constant is given by , where is the fine structure constant in ballistic graphene, and is the effective background dielectric constant in the graphene plane ( for suspended graphene and for graphene on SiO with the other side exposed to air/vacuum). The distance depends on the “geometry” of the graphene sheet. In a planar geometry, , one has to regularize the interaction term by choosing a finite screening length , which requires that the extension of the leads in is much smaller than in order to fulfill the charge neutrality condition. Alternatively one can let by considering a graphene sheet in the form of a cylinder with circumference . This choice yields the same results for with a computational advantage.

The self-consistent scheme is closed by relating the total charge density to the potential profile . The quantum-mechanical expression for in the self-consistent Hartree (SCH) approximation reads

(10) |

The local density of states (LDOS), , is calculated by solving the equation for the retarded Green’s function, ,

(11) |

where is a small positive parameter. We note that the dependence of the LDOS on the effective potential is non-local. Taking advantage of the translational invariance in we introduce the Fourier transform

(12) |

where the summation runs over the discrete values of conserved transverse momentum, , with integer (periodic boundary conditions in direction).

The boundary conditions in for the Green’s function in the channel representation are obtained from the exact analytical solution for and , where is such that . Selecting the decaying solution ( for ) we take the limit . Still, since the spectrum is partially discrete (some modes are confined to the sample), it is necessary to keep finite in order to maintain computational stability of the scheme. We consider the limit such that coherence is preserved on the scale .

The numerical computation of the Green’s function for is demanding since a large energy range has to be considered. For the lower energy bound at which the integrand in Eq. (10) becomes negligible is roughly . To ensure the charge conservation one has to keep the parameters and in Eq. (10) identical for both and , so that

(13) |

instead of Eq. (4). The expression (13) converges to Eq. (4) in the limit , .

The relation between charge density and effective potential is local in the TF approximation, which dramatically improves the computational efficiency. Replacing in Eq. (10) by , setting the temperature to zero, and using the ideal DOS (4) leads to the semiclassical expression for ,

(14) |

which we regard as the charge density in the TF approximation.

The set of Eqs. (8),(6), along with Eq. (10) for SCH and Eq. (14) for TF approximations, is solved by means of an iterative algorithm. It is run until reaches a self-consistency with accuracy of 10 relative to . For strong interactions () and large potential steps the procedure requires the use of strong damping in each iteration to ensure convergence (the weight of a new iteration is roughly of the order of ).

Potential profiles for different values of but fixed , calculated at zero temperature, are shown in Fig. 2. The dimensionless potential in the TF approximation depends on the two parameters and . Such scaling is only approximate for the SCH potential. We also note that the TF potential transforms exactly as under the global transformation , which is not the case for the SCH potential.

Our results apply to large system sizes which are not accessible by the DFT models. For instance, for Al contacts ( eV) and nm one finds . The sample width is of minor influence as far as . Although the parameter determines the potential decay in the leads (but not the charge profile), it has no significant effect on transport properties as long as .

The validity of the TF approximation is governed by the semi-classical criterion , where is the de Broglie wavelength. This is equivalent to the condition ,Khomyakov:2010 () which is violated for the upper curve in Fig. 2 (). The notable difference between the TF and SCH potential in this case indicates the importance of non-local quantum effects: evanescent modes and Friedel oscillations. Evanescent modes strongly increase the charge density in a vicinity of the sample-lead interface for energies close to the Dirac point.Titov:2010 () This effect, however, is partially compensated in the considered geometry by the Friedel oscillations. The latter suppress the LDOS near the interface for energies far from the Dirac point. To illustrate the compensation we plot in Fig. 3 the quantum correction for energies (dominated by the evanescent modes) and (dominated by the Friedel oscillations). The partial compensation of the non-local quantum corrections makes the TF approximation reliable even outside its applicability range. Still the mentioned discrepancy between the TF and SCH potentials shows the limitations of local density approximations such as TF in positively doped samples.comment ()

To make a direct comparison of our self-consistent calculation to the TF analysis of Ref. Khomyakov:2010, , we calculate the TF potential for a single sample-lead interface at ( for and for ). We find that the TF potential decays for as . The exponent is given by for and otherwise (which means that can vary with position). A similar behavior is found in Ref. Khomyakov:2010, when the doping inside the sample, , is due to the charged impurities. Our results for gated graphene are different since no charge density pinning in the leads is assumed in Ref. Khomyakov:2010, .

To compare our results to the full DFT treatment of Ref. BarrazaLopez:2010, we plot in Fig. 4 the self-consistent TF potentials calculated from our model for the same set of parameters eV, , , and 3.4, 6.8, and 13.6 nm. Despite the simplicity of our model (notably, the absence of exchange interactions) the potentials agree well with those calculated in Ref. BarrazaLopez:2010, .

## Iii Transport properties

For a given potential profile one may calculate the zero-temperature conductance from the Landauer formula , where the sum extends over the transverse momenta . The shot noise is quantified by the Fano factor .

The transmission probability for a given channel at the Fermi energy is related to the Green’s function by the Kubo formulaOstrovsky:2007 ()

(15) |

where is an optional cross-section. We choose such that for .

The conductance in the TF approximation is shown in Fig. 5 as a function of the charge density in the sample, , for different values of (using the SCH approximation gives very similar results). We note that the conductance scales approximately with the parameter for . A similar plot of the Fano factor is given in Fig. 6.

The most evident consequence of the charge transfer between the sample and the leads is an electron-hole asymmetry in the dependence of the conductance and noise on the charge density (or the gate voltage). To illustrate the key role of the slow potential decay we plot in Figs. 5,6, for comparison, the results obtained from the step-function model ( for and for ) with thin solid lines. The step-function model also leads to an asymmetry, which is, however, negligible in the experimentally relevant regime , .

We find that the conductance is enhanced for negative doping since in the entire sample and suppressed for sufficiently large positive doping since in this case. The minimal conductance is gradually increasing for large , while the position of the conductivity minimum is shifted towards the positive doping.

Another consequence of the slow potential decay is the specific form of the Fabry-Pérot oscillations, which are only visible on one side of the conductance minimum (in the Fano factor one can also observe weak oscillations for negative doping). As demonstrated by thin solid lines in Fig. 6, a small asymmetry already exists in the step-function model with finite doping of the leads. However, the oscillation amplitude is strongly enhanced (reduced) for positive (negative) doping as a consequence of the slow potential decay. The enhancement is strongest when the position of the Dirac point, , coincides with the chemical potential in a spatially extended region, which (as illustrated in Fig. 7) leads to the selective transmission of the carries with small transverse momenta at the n-p interfaces.Cheianov:2006 ()

In a particular set of experiments it is the sample length, , rather than the lead doping, that is varied. In Fig. 8 we rearrange the data to show the dependence of conductance on the sample doping for different sample lengths. As expected the contact-induced electron-hole asymmetry is more pronounced for shorter samples. The shift of the conductance minimum qualitatively agrees with recent measurements reported in Ref. Nouchi:2011, .

## Iv Conclusions

In conclusion we have modeled two-terminal electron transport through ballistic graphene samples by taking into account the effects of charge transfer at the metal-contact/graphene interfaces. Our analysis explains the electron-hole asymmetry in the conductance and the Fabry-Pérot conductance oscillations at positive doping, which have been observed in many experiments.Heersche:2007 (); Du:2008 (); Du:2008:NN (); Du:2009 () These phenomena are most clearly resolved in Ref. Du:2008:NN, (Fig. 2 b,c) for two different sample lengths. The period of the oscillations, found by transforming the positions of the peaks into , corresponds to , which unambiguously confirms that the oscillations are of the Fabry-Pérot type. The amplitude of the oscillations in this experiment is somewhat stronger than in Fig. 5, which is likely due to a modification of graphene by the metal leads that has resulted in an increased reflectivity of the sample-lead interfaces.

We are grateful to I. V. Gornyi, P. M. Ostrovsky, H. Schomerus, and J. Cayssol for discussions. WRH acknowledges funding by the Scottish Universities Physics Alliance. MJ acknowledges partial support from the Swedish VR and the Korean WCU program funded by MEST/NFR (Grant No. R31-2008-000-10057-0). MT acknowledges the hospitality of Karlsruhe Institute of Technology where part of the work was done.

### References

- X. Du, I. Skachko, F. Duerr, A. Luican, and E. Y. Andrei, Nature 462, 192 (2009); K. I. Bolotin, F. Ghahari, M. D. Shulman, H. L. Stormer, and P. Kim, Nature 462, 196 (2009).
- H. B. Heersche, P. Jarillo-Herrero, J. B. Oostinga, L. M. K. Vandersypen, and A. F. Morpurgo, Nature 446, 56 (2007); Solid State Commun. 143, 72 (2007).
- R. Danneau, F. Wu, M. F. Craciun, S. Russo, M. Y. Tomi, J. Salmilehto, A. F. Morpurgo, and P. J. Hakonen, Phys. Rev. Lett. 100, 196802 (2008).
- X. Du, I. Skachko, and E. Y. Andrei, Int. J. Mod. Phys. B 22, 4579 (2008).
- X. Du, I. Skachko, A. Barker, and E. Y. Andrei, Nature Nanotech. 3, 491 (2008).
- J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and C. W. J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006).
- M. Titov and C. W. J. Beenakker, Phys. Rev. B 74, 041401(R) (2006).
- E. J. H. Lee, K. Balasubramanian, R. T. Weitz, M. Burghard, and K. Kern, Nature Nanotech. 8, 486 (2008).
- B. Huard, N. Stander, J. A. Sulpizio, and D. Goldhaber-Gordon, Phys. Rev. B 78, 121402(R) (2008).
- P. Blake, R. Yang, S. V. Morozov, F. Schedin, L. A. Ponomarenko, A. A. Zhukov, R. R. Nair, I. V. Grigorieva, K. S. Novoselov, and A. K. Geim, Solid State Commun. 149, 1068 (2009).
- T. Mueller, F. Xia, M. Freitag, J. Tsang, and Ph. Avouris, Phys. Rev. B 79, 245430 (2009).
- S. Russo, M.F. Craciun, M. Yamamoto, A.F. Morpurgo, and S. Tarucha. Physica E 42, 677 (2010).
- P. A. Khomyakov, G. Giovannetti, P. C. Rusu, G. Brocks, J. van den Brink, and P. J. Kelly, Phys. Rev. B 79, 195425 (2009); G. Giovannetti, P. A. Khomyakov, G. Brocks, V. M. Karpan, J. van den Brink, and P. J. Kelly, Phys. Rev. Lett. 101, 026803 (2008).
- P. A. Khomyakov, A. A. Starikov, G. Brocks, and P. J. Kelly, Phys. Rev. B 82, 115437 (2010).
- M. Vanin, J. J. Mortensen, A. K. Kelkkanen, J. M. Garcia-Lastra, K. S. Thygesen, and K. W. Jacobsen, Phys. Rev. B 81, 081408(R) (2010).
- L. M. Zhang and M. M. Fogler, Phys. Rev. Lett. 100, 116804 (2008).
- R. Golizadeh-Mojarad and S. Datta, Phys. Rev. B 79, 085410 (2009).
- S. Barraza-Lopez, M. Vanević, M. Kindermann, and M. Y. Chou, Phys. Rev. Lett. 104, 076807 (2010).
- R. Nouchi and K. Tanigaki, Appl. Phys. Lett. 96, 253503 (2010).
- R. Nouchi, M. Shiraishi, and Y. Suzuki, Appl. Phys. Lett. 93, 152104 (2008)
- R. Nouchi, T. Saito, and K. Tanigaki, arXiv:1102.2954 (2011).
- J. Cayssol, B. Huard, and D. Goldhaber-Gordon, Phys. Rev. B 79, 075428 (2009).
- V. V. Cheianov and V. I. Fal’ko, Phys. Rev. B 74, 041403 (2006).
- P. Darancet, V. Olevano, and D. Mayou, Phys. Rev. Lett 102, 136803 (2009).
- M. Titov, P. M. Ostrovsky, and I. V. Gornyi, Semicond. Sci. Technol. 25, 034007 (2010).
- Strictly speaking, there exists a purely technical reason for the difference between TF and SCH potentials since the latter is calculated for a finite broadening . This effect is, however, negligible in Fig. 2.
- P. M. Ostrovsky, I. V. Gornyi, and A. D. Mirlin. Eur. Phys. J. Special Topics 148, 63 (2007).