-
Dominant source of disorder in graphene: Charged impurities or ripples?
Authors:
Zheyong Fan,
Andreas Uppstu,
Ari Harju
Abstract:
Experimentally produced graphene sheets exhibit a wide range of mobility values. Both extrinsic charged impurities and intrinsic ripples (corrugations) have been suggested to induce long-range disorder in graphene and could be a candidate for the dominant source of disorder. Here, using large-scale molecular dynamics and quantum transport simulations, we find that the hopping disorder and the gaug…
▽ More
Experimentally produced graphene sheets exhibit a wide range of mobility values. Both extrinsic charged impurities and intrinsic ripples (corrugations) have been suggested to induce long-range disorder in graphene and could be a candidate for the dominant source of disorder. Here, using large-scale molecular dynamics and quantum transport simulations, we find that the hopping disorder and the gauge and scalar potentials induced by the ripples are short-ranged, in strong contrast with predictions by continuous models, and the transport fingerprints of the ripple disorder are very different from those of charged impurities. We conclude that charged impurities are the dominant source of disorder in most graphene samples, whereas scattering by ripples is mainly relevant in the high carrier density limit of ultraclean graphene samples (with a charged impurity concentration < 10 ppm) at room and higher temperatures.
△ Less
Submitted 12 May, 2016;
originally announced May 2016.
-
Silicon and silicon-nitrogen impurities in graphene: structure, energetics and effects on electronic transport
Authors:
Mikko M. Ervasti,
Zheyong Fan,
Andreas Uppstu,
Arkady Krasheninnikov,
Ari Harju
Abstract:
We theoretically study the atomic structure and energetics of silicon and silicon-nitrogen impurities in graphene. Using density-functional theory, we get insight into the atomic structures of the impurities, evaluate their formation energies and assess their abundance in realistic samples. We find that nitrogen, as well as oxygen and hydrogen, are trapped at silicon impurities, considerably alter…
▽ More
We theoretically study the atomic structure and energetics of silicon and silicon-nitrogen impurities in graphene. Using density-functional theory, we get insight into the atomic structures of the impurities, evaluate their formation energies and assess their abundance in realistic samples. We find that nitrogen, as well as oxygen and hydrogen, are trapped at silicon impurities, considerably altering the electronic properties of the system. Furthermore, we show that nitrogen doping can induce local magnetic moments resulting in spin-dependent transport properties, even though neither silicon nor nitrogen impurities are magnetic by themselves. To simulate large systems with many randomly distributed impurities, we derive tight-binding models that describe the effects of the impurities on graphene π electron structure. Then by using the linear-scaling real-space Kubo-Greenwood method, we evaluate the transport properties of large-scale systems with random distribution of impurities, and find the fingerprint-like scattering cross sections for each impurity type. The transport properties vary widely, and our results indicate that some of the impurities can even induce strong localization in realistic graphene samples.
△ Less
Submitted 31 August, 2015;
originally announced September 2015.
-
Electronic and transport properties in geometrically disordered graphene antidot lattices
Authors:
Zheyong Fan,
Andreas Uppstu,
Ari Harju
Abstract:
A graphene antidot lattice, created by a regular perforation of a graphene sheet, can exhibit a considerable band gap required by many electronics devices. However, deviations from perfect periodicity are always present in real experimental setups and can destroy the band gap. Our numerical simulations, using an efficient linear-scaling quantum transport simulation method implemented on graphics p…
▽ More
A graphene antidot lattice, created by a regular perforation of a graphene sheet, can exhibit a considerable band gap required by many electronics devices. However, deviations from perfect periodicity are always present in real experimental setups and can destroy the band gap. Our numerical simulations, using an efficient linear-scaling quantum transport simulation method implemented on graphics processing units, show that disorder that destroys the band gap can give rise to a transport gap caused by Anderson localization. The size of the defect induced transport gap is found to be proportional to the radius of the antidots and inversely proportional to the square of the lattice periodicity. Furthermore, randomness in the positions of the antidots is found to be more detrimental than randomness in the antidot radius. The charge carrier mobilities are found to be very small compared to values found in pristine graphene, in accordance with recent experiments.
△ Less
Submitted 22 December, 2014; v1 submitted 19 December, 2014;
originally announced December 2014.
-
Anderson localization in two-dimensional graphene with short-range disorder: One-parameter scaling and finite-size effects
Authors:
Zheyong Fan,
Andreas Uppstu,
Ari Harju
Abstract:
We study Anderson localization in graphene with short-range disorder using the real-space Kubo-Greenwood method implemented on graphics processing units. Two models of short-range disorder, namely, the Anderson on-site disorder model and the vacancy defect model, are considered. For graphene with Anderson disorder, localization lengths of quasi-one-dimensional systems with various disorder strengt…
▽ More
We study Anderson localization in graphene with short-range disorder using the real-space Kubo-Greenwood method implemented on graphics processing units. Two models of short-range disorder, namely, the Anderson on-site disorder model and the vacancy defect model, are considered. For graphene with Anderson disorder, localization lengths of quasi-one-dimensional systems with various disorder strengths, edge symmetries, and boundary conditions are calculated using the real-space Kubo-Greenwood formalism, showing excellent agreement with independent transfer matrix calculations and superior computational efficiency. Using these data, we demonstrate the applicability of the one-parameter scaling theory of localization length and propose an analytical expression for the scaling function, which provides a reliable method of computing the two-dimensional localization length. This method is found to be consistent with another widely used method which relates the two-dimensional localization length to the elastic mean free path and the semiclassical conductivity. Abnormal behavior at the charge neutrality point is identified and interpreted to be caused by finite-size effects when the system width is comparable to or smaller than the elastic mean free path. We also demonstrate the finite-size effect when calculating the two-dimensional conductivity in the localized regime and show that a renormalization group beta function consistent with the one-parameter scaling theory can be extracted numerically. For graphene with vacancy disorder, we show that the proposed scaling function of localization length also applies. Lastly, we discuss some ambiguities in calculating the semiclassical conductivity around the charge neutrality point due to the presence of resonant states.
△ Less
Submitted 30 June, 2014; v1 submitted 30 October, 2013;
originally announced October 2013.
-
Obtaining localization properties efficiently using the Kubo-Greenwood formalism
Authors:
Andreas Uppstu,
Zheyong Fan,
Ari Harju
Abstract:
We establish, through numerical calculations and comparisons with a recursive Green's function based implementation of the Landauer-Büttiker formalism, an efficient method for studying Anderson localization in quasi-one-dimensional and two-dimensional systems using the Kubo-Greenwood formalism. Although the recursive Green's function method can be used to obtain the localization length of a mesosc…
▽ More
We establish, through numerical calculations and comparisons with a recursive Green's function based implementation of the Landauer-Büttiker formalism, an efficient method for studying Anderson localization in quasi-one-dimensional and two-dimensional systems using the Kubo-Greenwood formalism. Although the recursive Green's function method can be used to obtain the localization length of a mesoscopic conductor, it is numerically very expensive for systems that contain a large number of atoms transverse to the transport direction. On the other hand, linear-scaling has been achieved with the Kubo-Greenwood method, enabling the study of effectively two-dimensional systems. While the propagating length of the charge carriers will eventually saturate to a finite value in the localized regime, the conductances given by the Kubo-Greenwood method and the recursive Green's function method agree before the saturation. The converged value of the propagating length is found to be directly proportional to the localization length obtained from the exponential decay of the conductance.
△ Less
Submitted 24 February, 2014; v1 submitted 28 October, 2013;
originally announced October 2013.
-
Efficient linear-scaling quantum transport calculations on graphics processing units and applications on electron transport in graphene
Authors:
Zheyong Fan,
Andreas Uppstu,
Topi Siro,
Ari Harju
Abstract:
We implement, optimize, and validate the linear-scaling Kubo-Greenwood quantum transport simulation on graphics processing units by examining resonant scattering in graphene. We consider two practical representations of the Kubo-Greenwood formula: a Green-Kubo formula based on the velocity auto-correlation and an Einstein formula based on the mean square displacement. The code is fully implemented…
▽ More
We implement, optimize, and validate the linear-scaling Kubo-Greenwood quantum transport simulation on graphics processing units by examining resonant scattering in graphene. We consider two practical representations of the Kubo-Greenwood formula: a Green-Kubo formula based on the velocity auto-correlation and an Einstein formula based on the mean square displacement. The code is fully implemented on graphics processing units with a speedup factor of up to 16 (using double-precision) relative to our CPU implementation. We compare the kernel polynomial method and the Fourier transform method for the approximation of the Dirac delta function and conclude that the former is more efficient. In the ballistic regime, the Einstein formula can produce the correct quantized conductance of one-dimensional graphene nanoribbons except for an overshoot near the band edges. In the diffusive regime, the Green-Kubo and the Einstein formalisms are demonstrated to be equivalent. A comparison of the length-dependence of the conductance in the localization regime obtained by the Einstein formula with that obtained by the non-equilibrium Green's function method reveals the challenges in defining the length in the Kubo-Greenwood formalism at the strongly localized regime.
△ Less
Submitted 1 July, 2013;
originally announced July 2013.
-
Electronic states in finite graphene nanoribbons: Effect of charging and defects
Authors:
M. Ijäs,
M. Ervasti,
A. Uppstu,
P. Liljeroth,
J. van der Lit,
I. Swart,
A. Harju
Abstract:
We study the electronic structure of finite armchair graphene nanoribbons using density-functional theory and the Hubbard model, concentrating on the states localized at the zigzag termini. We show that the energy gaps between end-localized states are sensitive to doping, and that in doped systems, the gap between the end-localized states decreases exponentially as a function of the ribbon length.…
▽ More
We study the electronic structure of finite armchair graphene nanoribbons using density-functional theory and the Hubbard model, concentrating on the states localized at the zigzag termini. We show that the energy gaps between end-localized states are sensitive to doping, and that in doped systems, the gap between the end-localized states decreases exponentially as a function of the ribbon length. Doping also quenches the antiferromagnetic coupling between the end-localized states leading to a spin-split gap in neutral ribbons. By comparing dI/dV maps calculated using the many-body Hubbard model, its mean-field approximation and density-functional theory, we show that the use of a single-particle description is justified for graphene π states. Furthermore, we study the effect of structural defects in the ribbons on their electronic structure. Defects at one ribbon termini do not significantly modify the electronic states localized at the intact end. This provides further evidence for the interpretation of a multi-peaked structure in a recent scanning tunneling spectroscopy (STS) experiment resulting from inelastic tunneling processes [J. van der Lit et al., Nature Commun., in press (2013)]. Finally, we show that the hydrogen termination at the flake edges leaves identifiable fingerprints on the positive bias side of STS measurements, thus possibly aiding the experimental identification of graphene structures.
△ Less
Submitted 12 June, 2013;
originally announced June 2013.
-
Single- and multi-mode Fabry-Pérot interference in suspended graphene
Authors:
Mika Oksanen,
Andreas Uppstu,
Antti Laitinen,
Daniel J. Cox,
Monica Craciun,
Saverio Russo,
Ari Harju,
Pertti Hakonen
Abstract:
Phase coherence of charge carriers leads to electron-wave interference in ballistic mesoscopic conductors. In graphene, such Fabry-Pérot-like interference has been observed, but a detailed analysis has been complicated by the two-dimensional nature of conduction, which allows for complex interference patterns. In this work, we have achieved high-quality Fabry-Pérot interference in a suspended grap…
▽ More
Phase coherence of charge carriers leads to electron-wave interference in ballistic mesoscopic conductors. In graphene, such Fabry-Pérot-like interference has been observed, but a detailed analysis has been complicated by the two-dimensional nature of conduction, which allows for complex interference patterns. In this work, we have achieved high-quality Fabry-Pérot interference in a suspended graphene device, both in conductance and in shot noise, and analyzed their structure using Fourier transform techniques. The Fourier analysis reveals two sets of overlapping, coexisting interferences, with the ratios of the diamonds being equal to the width to length ratio of the device. We attribute these sets to a unique coexistence of longitudinal and transverse resonances, with the longitudinal resonances originating from the bunching of modes with low transverse momentum. Furthermore, our results give insight into the renormalization of the Fermi velocity in suspended graphene samples, caused by unscreened many-body interactions.
△ Less
Submitted 5 June, 2013;
originally announced June 2013.
-
Ab-initio transport fingerprints for resonant scattering in graphene
Authors:
Karri Saloriutta,
Andreas Uppstu,
Ari Harju,
Martti J. Puska
Abstract:
We have recently shown that by using a scaling approach for randomly distributed topological defects in graphene, reliable estimates for transmission properties of macroscopic samples can be calculated based even on single-defect calculations [A. Uppstu et al., Phys. Rev. B 85, 041401 (2012)]. We now extend this approach of energy-dependent scattering cross sections to the case of adsorbates on gr…
▽ More
We have recently shown that by using a scaling approach for randomly distributed topological defects in graphene, reliable estimates for transmission properties of macroscopic samples can be calculated based even on single-defect calculations [A. Uppstu et al., Phys. Rev. B 85, 041401 (2012)]. We now extend this approach of energy-dependent scattering cross sections to the case of adsorbates on graphene by studying hydrogen and carbon adatoms as well as epoxide and hydroxyl groups. We show that a qualitative understanding of resonant scattering can be gained through density functional theory results for a single-defect system, providing a transmission "fingerprint" characterizing each adsorbate type. This information can be used to reliably predict the elastic mean free path for moderate defect densities directly using ab-initio methods. We present tight-binding parameters for carbon and epoxide adsorbates, obtained to match the density-functional theory based scattering cross sections.
△ Less
Submitted 30 November, 2012;
originally announced November 2012.
-
High-field magnetoresistance revealing scattering mechanisms in graphene
Authors:
Andreas Uppstu,
Ari Harju
Abstract:
We show that the type of charge carrier scattering significantly affects the high-field magnetoresistance of graphene nanoribbons. This effect has potential to be used in identifying the scattering mechanisms in graphene. The results also provide an explanation for the experimentally found, intriguing differences in the behavior of the magnetoresistance of graphene Hall bars placed on different su…
▽ More
We show that the type of charge carrier scattering significantly affects the high-field magnetoresistance of graphene nanoribbons. This effect has potential to be used in identifying the scattering mechanisms in graphene. The results also provide an explanation for the experimentally found, intriguing differences in the behavior of the magnetoresistance of graphene Hall bars placed on different substrates. Additionally, our simulations indicate that the peaks in the longitudinal resistance tend to become pinned to fractionally quantized values, as different transport modes have very different scattering properties.
△ Less
Submitted 13 September, 2012; v1 submitted 2 May, 2012;
originally announced May 2012.
-
Electronic transport in graphene-based structures: an effective cross section approach
Authors:
Andreas Uppstu,
Karri Saloriutta,
Ari Harju,
Martti Puska,
Antti-Pekka Jauho
Abstract:
We show that transport in low-dimensional carbon structures with finite concentrations of scatterers can be modeled by utilising scaling theory and effective cross sections. Our reults are based on large scale numerical simulations of carbon nanotubes and graphene nanoribbons, using a tightbinding model with parameters obtained from first principles electronic structure calculations. As shown by a…
▽ More
We show that transport in low-dimensional carbon structures with finite concentrations of scatterers can be modeled by utilising scaling theory and effective cross sections. Our reults are based on large scale numerical simulations of carbon nanotubes and graphene nanoribbons, using a tightbinding model with parameters obtained from first principles electronic structure calculations. As shown by a comprehensive statistical analysis, the scattering cross sections can be used to estimate the conductance of a quasi-1D system both in the Ohmic and localized regimes. They can be computed with good accuracy from the transmission functions of single defects, greatly reducing the computational cost and paving the way towards using first principles methods to evaluate the conductance of mesoscopic systems, consisting of millions of atoms.
△ Less
Submitted 1 November, 2011;
originally announced November 2011.
-
Quantum confined electronic states in atomically well-defined graphene nanostructures
Authors:
Sampsa Hämäläinen,
Zhixiang Sun,
Mark P. Boneschanscher,
Andreas Uppstu,
Mari Ijäs,
Ari Harju,
Daniël Vanmaekelbergh,
Peter Liljeroth
Abstract:
Despite the enormous interest in the properties of graphene and the potential of graphene nanostructures in electronic applications, the study of quantum confined states in atomically well-defined graphene nanostructures remains an experimental challenge. Here, we study graphene quantum dots (GQDs) with well-defined edges in the zigzag direction, grown by chemical vapor deposition (CVD) on an irid…
▽ More
Despite the enormous interest in the properties of graphene and the potential of graphene nanostructures in electronic applications, the study of quantum confined states in atomically well-defined graphene nanostructures remains an experimental challenge. Here, we study graphene quantum dots (GQDs) with well-defined edges in the zigzag direction, grown by chemical vapor deposition (CVD) on an iridium(111) substrate, by low-temperature scanning tunneling microscopy (STM) and spectroscopy (STS). We measure the atomic structure and local density of states (LDOS) of individual GQDs as a function of their size and shape in the range from a couple of nanometers up to ca. 20 nm. The results can be quantitatively modeled by a relativistic wave equation and atomistic tight-binding calculations. The observed states are analogous to the solutions of the text book "particle-in-a-box" problem applied to relativistic massless fermions.
△ Less
Submitted 20 October, 2011; v1 submitted 19 October, 2011;
originally announced October 2011.