Abstract: Recent years have seen a tremendous rise of two-dimensional (2D) magnetic materials, several of which were verified experimentally. However, most of the theoretical predictions to date rely on ab initio methods, at zero temperature and fluctuation-free, while one certainly expects detrimental quantum fluctuations at finite temperatures. Here, we present the solution of the quantum Heisenberg model for honeycomb/hexagonal lattices with anisotropic exchange interaction up to third nearest neighbors and in an applied field in arbitrary direction, which answers the question whether long-range magnetization can indeed survive in the ultrathin limit of materials, up to which temperature, and what the characteristic excitation (magnon) frequencies are, all essential to envisaged applications of magnetic 2D materials. We find that long-range magnetic order persists at finite temperature for materials with overall easy-axis anisotropy. We validate the calculations on the examples of monolayers CrI3, CrBr3, and MnSe2. Moreover, we provide an easy-to-use tool to calculate Curie temperatures of new 2D computational materials.

Abstract: 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.

Abstract: Recently, we have used a projection operator to fix the number of particles in a second quantization approach in order to deal with the canonical ensemble. Having been applied earlier to handle various problems in nuclear physics that involve fixed particle numbers, the projector formalism was extended to grant access as well to quantum-statistical averages in condensed matter physics, such as particle densities and correlation functions. In this light, the occupation numbers of the subsequent single-particle energy eigenstates are key quantities to be examined. The goal of this paper is (1) to provide a sound extension of the projector formalism directly addressing the occupation numbers as well as the chemical potential, and (2) to demonstrate how the emerging problems related to numerical instability for fermions can be resolved to obtain the canonical statistical quantities for both fermions and bosons. (C) 2018 Elsevier B.V. All rights reserved.

Abstract: In this contribution, we report a fundamental study of the factors that set the contact resistivity between metals and highly doped n-type 2D and 3D semiconductors. We investigate the case of n-type doped Si contacted with amorphous TiSi combining first principles calculations with Non-Equilibrium Green functions transport simulations. The evolution of the intrinsic contact resistivity with the doping concentration is found to saturate at similar to 2 x 10(-10) Omega.cm(2) for the case of TiSi and imposes an intrinsic limit to the ultimate contact resistance achievable for n-doped Silamorphous-TiSi (aTiSi). The limit arises from the intrinsic properties of the semiconductors and of the metals such as their electron effective masses and Fermi energies. We illustrate that, in this regime, contacting heavy electron effective mass metals with semiconductor helps reducing the interface intrinsic contact resistivity. This observation seems to hold true regardless of the 3D character of the semiconductor, as illustrated for the case of three 2D semiconducting materials, namely MoS2, ZrS2 and HfS2. (C) The Author(s) 2018. Published by ECS.

Abstract: It is generally understood that the resistivity of metal thin films scales with film thickness mainly due to grain boundary and boundary surface scattering. Recently, several experiments and ab initio simulations have demonstrated the impact of crystal orientation on resistivity scaling. The crystal orientation cannot be captured by the commonly used resistivity scaling models and a qualitative understanding of its impact is currently lacking. In this work, we derive a resistivity scaling model that captures grain boundary and boundary surface scattering as well as the anisotropy of the band structure. The model is applied to Cu and Ru thin films, whose conduction bands are (quasi-) isotropic and anisotropic, respectively. After calibrating the anisotropy with ab initio simulations, the resistivity scaling models are compared to experimental resistivity data and a renormalization of the fitted grain boundary reflection coefficient can be identified for textured Ru.

Abstract: In this contribution, we report a fundamental study of the factors that set the contact resistivity between metals and highly doped semiconductors. We investigate the case of n-type doped Si contacted with amorphous TiSi combining first-principles calculations with Non-Equilibrium Green functions transport simulations. The intrinsic contact resistivity is found to saturate at similar to 2x10(-10) Omega.cm(2) with the doping concentration and sets an intrinsic limit to the ultimate contact resistance achievable for n-doped Si vertical bar amorphous-TiSi. This limit arises from the intrinsic properties of the semiconductor and of the metal such as their electron effective masses and Fermi energies. We illustrate that, in this regime, contacting metals with a heavy electron effective mass helps reducing the interface intrinsic contact resistivity.

Abstract: Knowing the exact number of particles N, and taking this knowledge into account, the quantum canonical ensemble imposes a constraint on the occupation number operators. The constraint particularly hampers the systematic calculation of the partition function and any relevant thermodynamic expectation value for arbitrary but fixed N. On the other hand, fixing only the average number of particles, one may remove the above constraint and simply factorize the traces in Fock space into traces over single-particle states. As is well known, that would be the strategy of the grand-canonical ensemble which, however, comes with an additional Lagrange multiplier to impose the average number of particles. The appearance of this multiplier can be avoided by invoking a projection operator that enables a constraint-free computation of the partition function and its derived quantities in the canonical ensemble, at the price of an angular or contour integration. Introduced in the recent past to handle various issues related to particle-number projected statistics, the projection operator approach proves beneficial to a wide variety of problems in condensed matter physics for which the canonical ensemble offers a natural and appropriate environment. In this light, we present a systematic treatment of the canonical ensemble that embeds the projection operator into the formalism of second quantization while explicitly fixing N, the very number of particles rather than the average. Being applicable to both bosonic and fermionic systems in arbitrary dimensions, transparent integral representations are provided for the partition function Z(N) and the Helmholtz free energy F-N as well as for two- and four-point correlation functions. The chemical potential is not a Lagrange multiplier regulating the average particle number but can be extracted from FN+1 – F-N, as illustrated for a two-dimensional fermion gas. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: Heterostructures of III-V materials under a mechanical strain are being actively researched to enhance the performance of the tunnel field-effect transistor (TFET). In scaled III-V device structures, however, the interplay between the effects of strain and quantum confinement on the semiconductor band structure and hence the performance is highly non-trivial. We have therefore developed a computationally efficient quantum mechanical simulator Pharos, which enables self-consistent full-zone k.p-based simulations of III-V TFETs under a general non-uniform strain. We present the self-consistent procedure and demonstrate it on confined staggered bandgap GaAs0.5Sb0.5/In0.53Ga0.47As TFETs. We find a large performance degradation due to size-induced quantum confinement compared to non-confined devices. We show that some performance can be regained either by applying a uniform biaxial tensile strain or through the non-uniform strain profile at a lattice-mismatched heterostructure.

Abstract: A weakly coupled system of two crossed graphene nanoribbons exhibits direct tunneling due to the overlap of the wavefunctions of both ribbons. We apply the Bardeen transfer Hamiltonian formalism, using atomistic band structure calculations to account for the effect of the atomic structure on the tunneling process. The strong quantum-size confinement of the nanoribbons is mirrored by the one-dimensional character of the electronic structure, resulting in properties that differ significantly from the case of inter-layer tunneling, where tunneling occurs between bulk two-dimensional graphene sheets. The current-voltage characteristics of the inter-ribbon tunneling structures exhibit resonance, as well as stepwise increases in current. Both features are caused by the energetic alignment of one-dimensional peaks in the density-of-states of the ribbons. Resonant tunneling occurs if the sign of the curvature of the coupled energy bands is equal, whereas a step-like increase in the current occurs if the signs are opposite. Changing the doping modulates the onset-voltage of the effects as well as their magnitude. Doping through electrostatic gating makes these structures promising for application towards steep slope switching devices. Using the atomistic empirical pseudopotentials based Bardeen transfer Hamiltonian method, inter-ribbon tunneling can be studied for the whole range of two-dimensional materials, such as transition metal dichalcogenides. The effects of resonance and of step-like increases in the current we observe in graphene ribbons are also expected in ribbons made from these alternative two-dimensional materials, because these effects are manifestations of the one-dimensional character of the density-of-states. Published by AIP Publishing.

Abstract: Efficient quantum mechanical simulation of tunnel field-effect transistors (TFETs) is indispensable to allow for an optimal configuration identification. We therefore present a full-zone 15-band quantum mechanical solver based on the envelope function formalism and employing a spectral method to reduce computational complexity and handle spurious solutions. We demonstrate the versatility of the solver by simulating a 40 nm wide In0.53Ga0.47As lineTFET and comparing it to p-n-i-n configurations with various pocket and body thicknesses. We find that the lineTFET performance is not degraded compared to semi-classical simulations. Furthermore, we show that a suitably optimized p-n-i-n TFET can obtain similar performance to the lineTFET. (C) 2015 AIP Publishing LLC.

Abstract: The envelope function method traditionally employs a single basis set which, in practice, relates to a single material because the k.p matrix elements are generally only known in a particular basis. In this work, we defined a basis function transformation to alleviate this restriction. The transformation is completely described by the known inter-band momentum matrix elements. The resulting envelope function equation can solve the electronic structure in lattice matched heterostructures without resorting to boundary conditions at the interface between materials, while all unit-cell averaged observables can be calculated as with the standard envelope function formalism. In the case of two coupled bands, this heterostructure formalism is equivalent to the standard formalism while taking position dependent matrix elements. (C) 2015 Elsevier B.V. All rights reserved

Abstract: The tunneling current between two crossed graphene ribbons is described invoking the empirical pseudopotential approximation and the Bardeen transfer Hamiltonian method. Results indicate that the density of states is the most important factor determining the tunneling current between small (similar to nm) ribbons. The quasi-one dimensional nature of graphene nanoribbons is shown to result in resonant tunneling.

Abstract: Ando's surface roughness model is applied to metallic nanowires and extended beyond small roughness size and infinite barrier limit approximations for the wavefunction overlaps, such as the Prange-Nee approximation. Accurate and fast simulations can still be performed without invoking these overlap approximations by averaging over roughness profiles using finite domain distribution functions to obtain an analytic solution for the scattering rates. The simulations indicate that overlap approximations, while predicting a resistivity that agrees more or less with our novel approach, poorly estimate the underlying scattering rates. All methods show that a momentum gap between left- and right-moving electrons at the Fermi level, surpassing a critical momentum gap, gives rise to a substantial decrease in resistivity.

Abstract: We study the resistivity scaling in nanometer-sized metallic wires due to surface roughness and grain-boundaries, currently the main cause of electron scattering in nanoscaled interconnects. The resistivity has been obtained with the Boltzmann transport equation, adopting the relaxation time approximation of the distribution function and the effective mass approximation for the conducting electrons. The relaxation times are calculated exactly, using Fermi's golden rule, resulting in a correct relaxation time for every sub-band state contributing to the transport. In general, the relaxation time strongly depends on the sub-band state, something that remained unclear with the methods of previous work. The resistivity scaling is obtained for different roughness and grain-boundary properties, showing large differences in scaling behavior and relaxation times. Our model clearly indicates that the resistivity is dominated by grain-boundary scattering, easily surpassing the surface roughness contribution by a factor of 10. (C) 2014 AIP Publishing LLC.

Abstract: Group IV based tunnel field-effect transistors generally show lower on-current than III-V based devices because of the weaker phonon-assisted tunneling transitions in the group IV indirect bandgap materials. Direct tunneling in Ge, however, can be enhanced by strain engineering. In this work, we use a 30-band k.p method to calculate the band structure of biaxial tensile strained Ge and then extract the bandgaps and effective masses at Gamma and L symmetry points in k-space, from which the parameters for the direct and indirect band-to-band tunneling (BTBT) models are determined. While transitions from the heavy and light hole valence bands to the conduction band edge at the L point are always bridged by phonon scattering, we highlight a new finding that only the light-holelike valence band is strongly coupling to the conduction band at the Gamma point even in the presence of strain based on the 30-band k.p analysis. By utilizing a Technology Computer Aided Design simulator equipped with the calculated band-to-band tunneling BTBT models, the electrical characteristics of tensile strained Ge point and line tunneling devices are self-consistently computed considering multiple dynamic nonlocal tunnel paths. The influence of field-induced quantum confinement on the tunneling onset is included. Our simulation predicts that an on-current up to 160 (260) mu A/mu m can be achieved along with on/off ratio > 10(6) for V-DD = 0.5V by the n-type (p-type) line tunneling device made of 2.5% biaxial tensile strained Ge. (C) 2014 AIP Publishing LLC.