We present novel methods implemented within the non-equilibrium Green function code (NEGF) transiesta based on density functional theory (DFT). Our flexible, next-generation DFT-NEGF code handles devices with one or multiple electrodes (Ne ≥ 1) with individual chemical potentials and electronic temperatures. We describe its novel methods for electrostatic gating, contour optimizations, and assertion of charge conservation, as well as the newly implemented algorithms for optimized and scalable matrix inversion, performance-critical pivoting, and hybrid parallellization. Additionally, a generic NEGF "post-processing" code (tbtrans/phtrans) for electron and phonon transport is presented with several novelties such as Hamiltonian interpolations, Ne ≥ 1 electrode capability, bond-currents, generalized interface for user-defined tight-binding transport, transmission projection using eigenstates of a projected Hamiltonian, and fast inversion algorithms for large-scale simulations easily exceeding 10 6 atoms on workstation computers. The new features of both codes are demonstrated and bench-marked for relevant test systems. * The transport of charge, magnetic moments and, in general, any sort of excitation is a fascinating fundamental physical problem that has demanded attention for a long time [1]. Today, the interest is enhanced by the technological needs of an industry increasingly based on devices whose detailed atomistic structure matters [2], but the treatment of transport is still a formidable open task. Spurred by the fast developments of the microelectronic industry, the first attempts to understand electronic transport at the atomic scale where based on scattering theory [3]. The electron transmission between two semi-infinite reservoirs was treated in a time-independent fashion solving the scattering matrix connecting the reservoirs. At this stage, transport was described as one-electron scattering by a static contact region and this granted access to many concepts and to devising new experiments [4][5][6]. However, the problem is fundamentally a non-equilibrium one that requires evolving many-body states [7][8][9][10].Density functional theory (DFT) has been one method to address some aspects of this problem. Conceptually, DFT is a mean-field many-body theory of the ground state. As such, it can in principle give exact results for the linear conductance because the linear response is a property of the ground state [11]. Beyond linear conductance, not even ideal DFT works because of the need to describe excited states and dynamics of the system. Such limitations may be mitigated by using time-dependent DFT [12,13], but going beyond the linear regime is highly nontrivial. A main issue of a DFT description stems from the approximations made to compute the ground state. Indeed, it has been recently shown that cases where strong correlations rein, such as the Coulomb blockade regime, the commonly used exchange-and-correlation functionals fail and new ones have to be used [14].Probably the most significant conceptu...

A review of the present status, recent enhancements, and applicability of the SIESTA program is presented. Since its debut in the mid-nineties, SIESTA's flexibility, efficiency and free distribution has given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of SIESTA combines finite-support pseudoatomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here we describe the more recent implementations on top of that core scheme, which include: full spin-orbit interaction, non-repeated and multiple-contact ballistic electron transport, DFT+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, densityfunctional perturbation theory, efficient Van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as WANNIER90 and the second-principles modelling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering SIESTA runs, and various postprocessing utilities. SIESTA has also been a) Electronic mail:

The electron microscope has been a powerful, highly versatile workhorse in the fields of material and surface science, micro and nanotechnology, biology and geology, for nearly 80 years. The advent of two-dimensional materials opens new possibilities for realizing an analogy to electron microscopy in the solid state. Here we provide a perspective view on how a two-dimensional (2D) Dirac fermion-based microscope can be realistically implemented and operated, using graphene as a vacuum chamber for ballistic electrons. We use semiclassical simulations to propose concrete architectures and design rules of 2D electron guns, deflectors, tunable lenses and various detectors. The simulations show how simple objects can be imaged with well-controlled and collimated in-plane beams consisting of relativistic charge carriers. Finally, we discuss the potential of such microscopes for investigating edges, terminations and defects, as well as interfaces, including external nanoscale structures such as adsorbed molecules, nanoparticles or quantum dots.

Designing platforms to control phase-coherence and interference of electron waves is a cornerstone for future quantum electronics, computing or sensing. Nanoporous graphene (NPG) consisting of linked graphene nanoribbons has recently been fabricated using molecular precursors and bottom-up assembly [Moreno et al, Science 360, 199 (2018)] opening an avenue for controlling the electronic current in a two-dimensional material. By simulating electron transport in real-sized NPG samples we predict that electron waves injected from the tip of a scanning tunneling microscope (STM) behave similarly to photons in coupled waveguides, displaying a Talbot interference pattern. We link the origins of this effect to the band structure of the NPG and further demonstrate how this pattern may be mapped out by a second STM probe. We enable atomistic parameter-free calculations beyond the 100 nm scale by developing a new multi-scale method where first-principles density functional theory regions are seamlessly embedded into a large-scale tight-binding. arXiv:1811.07576v1 [cond-mat.mes-hall]

Bottom-up prepared carbon nanostructures appear as promising platforms for future carbon-based nanoelectronics, due to their atomically precise and versatile structure. An important breakthrough is the recent preparation of nanoporous graphene (NPG) as an ordered covalent array of graphene nanoribbons (GNRs). Within NPG, the GNRs may be thought of as 1D electronic nanochannels through which electrons preferentially move, highlighting NPG's potential for carbon nanocircuitry. However, the π-conjugated bonds bridging the GNRs give rise to electronic cross-talk between the individual 1D channels, leading to spatially dispersing electronic currents. Here, we propose a chemical design of the bridges resulting in destructive quantum interference, which blocks the cross-talk between GNRs in NPG, electronically isolating them. Our multiscale calculations reveal that injected currents can remain confined within a single, 0.7 nm wide, GNR channel for distances as long as 100 nm. The concepts developed in this work thus provide an important ingredient for the quantum design of future carbon nanocircuitry. *