Efficient solution of the Wigner-Liouville equation using a spectral decomposition of the force field

The Wigner-Liouville equation is reformulated using a spectral decomposition of the classical force field instead of the potential energy. The latter is shown to simplify the Wigner-Liouville kernel both conceptually and numerically as the spectral force Wigner-Liouville equation avoids the numerical evaluation of the highly oscillatory Wigner kernel which is nonlocal in both position and momentum. The quantum mechanical evolution is instead governed by a term local in space and non-local in momentum, where the non-locality in momentum has only a limited range. An interpretation of the time evolution in terms of two processes is presented; a classical evolution under the influence of the averaged driving field, and a probability-preserving quantum-mechanical generation and annihilation term. Using the inherent stability and reduced complexity, a direct deterministic numerical implementation using Chebyshev and Fourier pseudo-spectral methods is detailed. For the purpose of illustration, we present results for the time-evolution of a one-dimensional resonant tunneling diode driven out of equilibrium

Theoretical study of scattering in graphene ribbons in the presence of structural and atomistic edge roughness

We investigate the diffusive electron-transport properties of charge-doped graphene ribbons and nanoribbons with imperfect edges. We consider different regimes of edge scattering, ranging from wide graphene ribbons with (partially) diffusive edge scattering to ribbons with large width variations and nanoribbons with atomistic edge roughness. For the latter, we introduce an approach based on pseudopotentials, allowing for an atomistic treatment of the band structure and the scattering potential, on the self-consistent solution of the Boltzmann transport equation within the relaxation-time approximation and taking into account the edge-roughness properties and statistics. The resulting resistivity depends strongly on the ribbon orientation, with zigzag (armchair) ribbons showing the smallest (largest) resistivity and intermediate ribbon orientations exhibiting intermediate resistivity values. The results also show clear resistivity peaks, corresponding to peaks in the density of states due to the confinement-induced subband quantization, except for armchair-edge ribbons that show a very strong width dependence because of their claromatic behavior. Furthermore, we identify a strong interplay between the relative position of the two valleys of graphene along the transport direction, the correlation profile of the atomistic edge roughness, and the chiral valley modes, leading to a peculiar strongly suppressed resistivity regime, most pronounced for the zigzag orientation.Comment: 13 pages, 7 figure

Monte Carlo Study of Electronic Transport in Monolayer InSe

The absence of a band gap in graphene makes it of minor interest for field-effect transistors. Layered metal chalcogenides have shown great potential in device applications thanks to their wide bandgap and high carrier mobility. Interestingly, in the ever-growing library of two-dimensional (2D) materials, monolayer InSe appears as one of the new promising candidates, although still in the initial stage of theoretical studies. Here, we present a theoretical study of this material using density functional theory (DFT) to determine the electronic band structure as well as the phonon spectrum and electron-phonon matrix elements. The electron-phonon scattering rates are obtained using Fermi’s Golden Rule and are used in a full-band Monte Carlo computer program to solve the Boltzmann transport equation (BTE) to evaluate the intrinsic low-field mobility and velocity-field characteristic. The electron-phonon matrix elements, accounting for both long- and short-range interactions, are considered to study the contributions of different scattering mechanisms. Since monolayer InSe is a polar piezoelectric material, scattering with optical phonons is dominated by the long-range interaction with longitudinal optical (LO) phonons while scattering with acoustic phonons is dominated by piezoelectric scattering with the longitudinal (LA) branch at room temperature (T = 300 K) due to a lack of a center of inversion symmetry in monolayer InSe. The low-field electron mobility, calculated considering all electron-phonon interactions, is found to be 110 cm2V−1s−1, whereas values of 188 cm2V−1s−1 and 365 cm2V−1s−1 are obtained considering the long-range and short-range interactions separately. Therefore, the calculated electron mobility of monolayer InSe seems to be competitive with other previously studied 2D materials and the piezoelectric properties of monolayer InSe make it a suitable material for a wide range of applications in next generation nanoelectronics