Source author record

Kurt Stokbro

Kurt Stokbro appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2016arXiv

General atomistic approach for modeling metal-semiconductor interfaces using density functional theory and non-equilibrium Green's function

Metal-semiconductor contacts are a pillar of modern semiconductor technology. Historically, their microscopic understanding has been hampered by the inability of traditional analytical and numerical methods to fully capture the complex physics governing their operating principles. Here we introduce an atomistic approach based on density functional theory and non-equilibrium Green's function, which includes all the relevant ingredients required to model realistic metal-semiconductor interfaces and allows for a direct comparison between theory and experiments via I-V bias curves simulations. We apply this method to characterize an Ag/Si interface relevant for photovoltaic applications and study the rectifying-to-Ohmic transition as function of the semiconductor doping.We also demonstrate that the standard "Activation Energy" method for the analysis of I-V bias data might be inaccurate for non-ideal interfaces as it neglects electron tunneling, and that finite-size atomistic models have problems in describing these interfaces in the presence of doping, due to a poor representation of space-charge effects. Conversely, the present method deals effectively with both issues, thus representing a valid alternative to conventional procedures for the accurate characterization of metal-semiconductor interfaces.

preprint2016arXiv

Inelastic vibrational signals in electron transport across graphene nanoconstrictions

We present calculations of the inelastic vibrational signals in the electrical current through a graphene nanoconstriction. We find that the inelastic signals are only present when the Fermi-level position is tuned to electron transmission resonances, thus, providing a fingerprint which can link an electron transmission resonance to originate from the nanoconstriction. The calculations are based on a novel first-principles method which includes the phonon broadening due to coupling with phonons in the electrodes. We find that the signals are modified due to the strong coupling to the electrodes, however, still remain as robust fingerprints of the vibrations in the nanoconstriction. We investigate the effect of including the full self-consistent potential drop due to finite bias and gate doping on the calculations and find this to be of minor importance.

preprint2015arXiv

First-principles method for electron-phonon coupling and electron mobility: Applications to 2D materials

We present density functional theory calculations of the phonon-limited mobility in n-type monolayer graphene, silicene and MoS$_2$. The material properties, including the electron-phonon interaction, are calculated from first-principles. We provide a detailed description of the normalized full-band relaxation time approximation for the linearized Boltzmann transport equation (BTE) that includes inelastic scattering processes. The bulk electron-phonon coupling is evaluated by a supercell method. The method employed is fully numerical and does therefore not require a semi-analytic treatment of part of the problem and, importantly, it keeps the anisotropy information stored in the coupling as well as the band structure. In addition, we perform calculations of the low-field mobility and its dependence on carrier density and temperature to obtain a better understanding of transport in graphene, silicene and monolayer MoS$_2$. Unlike graphene, the carriers in silicene show strong interaction with the out-of-plane modes. We find that graphene has more than an order of magnitude higher mobility compared to silicene. For MoS$_2$, we obtain several orders of magnitude lower mobilities in agreement with other recent theoretical results. The simulations illustrate the predictive capabilities of the newly implemented BTE solver applied in simulation tools based on first-principles and localized basis sets.

preprint2015arXiv

Nonequilibrium spin texture within a thin layer below the surface of current-carrying topological insulator Bi$_2$Se$_3$: A first-principles quantum transport study

We predict that unpolarized charge current injected into a ballistic thin film of prototypical topological insulator (TI) Bi$_2$Se$_3$ will generate a {\it noncollinear spin texture} $\mathbf{S}(\mathbf{r})$ on its surface. Furthermore, the nonequilibrium spin texture will extend into $\simeq 2$ nm thick layer below the TI surfaces due to penetration of evanescent wavefunctions from the metallic surfaces into the bulk of TI. Averaging $\mathbf{S}(\mathbf{r})$ over few Å along the longitudinal direction defined by the current flow reveals large component pointing in the transverse direction. In addition, we find an order of magnitude smaller out-of-plane component when the direction of injected current with respect to Bi and Se atoms probes the largest hexagonal warping of the Dirac-cone dispersion on TI surface. Our analysis is based on an extension of the nonequilibrium Green functions combined with density functional theory (NEGF+DFT) to situations involving noncollinear spins and spin-orbit coupling. We also demonstrate how DFT calculations with properly optimized local orbital basis set can precisely match putatively more accurate calculations with plane-wave basis set for the supercell of Bi$_2$Se$_3$.

preprint2014arXiv

Improved initial guess for minimum energy path calculations

A method is presented for generating a good initial guess of a transition path between given initial and final states of a system without evaluation of the energy. An objective function surface is constructed using an interpolation of pairwise distances at each discretization point along the path and the nudged elastic band method then used to find an optimal path on this image dependent pair potential (IDPP) surface. This provides an initial path for the more computationally intensive calculations of the true minimum energy path using some method of choice for evaluating the energy and atomic forces, for example by ab initio or density functional theory. The optimal path on the IDPP surface is significantly closer to the true minimum energy path than a linear interpolation of the Cartesian coordinates and, therefore, reduces the number of iterations needed to reach convergence and averts divergence in the electronic structure calculations when atoms are brought too close to each other in the initial path. The method is illustrated with three examples: (1) rotation of a methyl group in an ethane molecule, (2) an exchange of atoms in an island on a crystal surface, and (3) an exchange of two Si-atoms in amorphous silicon. In all three cases, the computational effort in finding the minimum energy path with DFT was reduced by a factor ranging from 50 % to an order of magnitude by using an IDPP path as the initial path. The time required for parallel computations was reduced even more because of load imbalance when linear interpolation of Cartesian coordinates was used.

preprint2013arXiv

Electron transport across a metal-organic interface

We simulate the electron transport across the Au(111)-pentacene interface using non-equilibrium Green's functions and density-functional theory (NEGF-DFT), and calculate the bias-dependent electron transmission. We find that the electrical contact resistance is dominated by the formation of a Schottky barrier at the interface, and show that the conventional semiconductor transport models across Schottky barriers need to be modified in order to describe the simulation data. We present an extension of the conventional Schottky barrier transport model, which can describe our simulation results and rationalize recent experimental data.

preprint2013arXiv

First principles study of pentacene on Au(111)

We investigate the atomic and electronic structure of a single layer of pentacene on the Au(111) surface using density functional theory. To find the candidate structures we strain match the pentacene crystal geometry with the Au(111) surface, in this way we find pentacene overlayer structures with a low strain. We show that the geometries obtained with this approach has lower energy than previous proposed surface geometries of pentacene on Au(111). We also show that the geometry and workfunction of the obtained structures are in excellent agreement with experimental data.

preprint2012arXiv

Atomic-scale model for the contact resistance of the nickel-graphene interface

We perform first-principles calculations of electron transport across a nickel-graphene interface. Four different geometries are considered, where the contact area, graphene and nickel surface orientations and the passivation of the terminating graphene edge are varied. We find covalent bond formation between the graphene layer and the nickel surface, in agreement with other theoretical studies. We calculate the energy-dependent electron transmission for the four systems and find that the systems have very similar edge contact resistance, independent of the contact area between nickel and graphene, and in excellent agreement with recent experimental data. A simple model where graphene is bonded with a metal surface shows that the results are generic for covalently bonded graphene, and the minimum attainable edge contact resistance is twice the ideal edge quantum contact resistance of graphene.

preprint2010arXiv

First-principles modelling of molecular single-electron transistors

We present a first-principles method for calculating the charging energy of a molecular single-electron transistor operating in the Coulomb blockade regime. The properties of the molecule are modeled using density-functional theory, the environment is described by a continuum model, and the interaction between the molecule and the environment are included through the Poisson equation. The model is used to calculate the charge stability diagrams of a benzene and C$_{60}$ molecular single-electron transistor.

preprint2010arXiv

Semi-Empirical Model for Nano-Scale Device Simulations

We present a new semi-empirical model for calculating electron transport in atomic-scale devices. The model is an extension of the Extended Hückel method with a self-consistent Hartree potential. This potential models the effect of an external bias and corresponding charge re-arrangements in the device. It is also possible to include the effect of external gate potentials and continuum dielectric regions in the device. The model is used to study the electron transport through an organic molecule between gold surfaces, and it is demonstrated that the results are in closer agreement with experiments than ab initio approaches provide. In another example, we study the transition from tunneling to thermionic emission in a transistor structure based on graphene nanoribbons.

preprint2009arXiv

Efficient wave function matching approach for quantum transport calculations

The Wave Function Matching (WFM) technique has recently been developed for the calculation of electronic transport in quantum two-probe systems. In terms of efficiency it is comparable with the widely used Green's function approach. The WFM formalism presented so far requires the evaluation of all the propagating and evanescent bulk modes of the left and right electrodes in order to obtain the correct coupling between device and electrode regions. In this paper we will describe a modified WFM approach that allows for the exclusion of the vast majority of the evanescent modes in all parts of the calculation. This approach makes it feasible to apply iterative techniques to efficiently determine the few required bulk modes, which allows for a significant reduction of the computational expense of the WFM method. We illustrate the efficiency of the method on a carbon nanotube field-effect-transistor (FET) device displaying band-to-band tunneling and modeled within the semi-empirical Extended Hückel theory (EHT) framework.