Abstract: To design fast memory devices, we need material combinations that can facilitate fast read and write operations. We present a heterostructure comprising a two-dimensional (2D) magnet and a 2D topological insulator (TI) as a viable option for designing fast memory devices. We theoretically model the spin-charge dynamics between 2D magnets and 2D TIs. Using the adiabatic approximation, we combine the nonequi-librium Green's function method for spin-dependent electron transport and a time-quantified Monte Carlo method for simulating magnetization dynamics. We show that it is possible to switch a magnetic domain of a ferromagnet using the spin torque from spin-polarized edge states of a 2D TI. We show further that the switching of 2D magnets by TIs is strongly dependent on the interface exchange (Jint), and an opti-mal interface exchange, is required for efficient switching. Finally, we compare experimentally grown Cr compounds and show that Cr compounds with higher anisotropy (such as CrI3) result in a lower switching speed but a more stable magnetic order.

Abstract: We compare three first-principles methods of calculating the Curie temperature in two-dimensional (2D) ferromagnetic materials (FM), modeled using the Heisenberg model, and propose a simple formula for estimating the Curie temperature with high accuracy that works for all common 2D lattice types. First, we study the effect of exchange anisotropy on the Curie temperature calculated using the Monte Carlo (MC), the Green's function, and the renormalized spin-wave (RNSW) methods. We find that the Green's function method overestimates the Curie temperature in high-anisotropy regimes compared to the MC method, whereas the RNSW method underestimates the Curie temperature compared to the MC and the Green's function methods. Next, we propose a closed-form formula for calculating the Curie temperature of 2D FMs, which provides an estimate of the Curie temperature that is greatly improved over the mean-field expression for magnetic material screening. We apply the closed-form formula to predict the Curie temperature 2D magnets screened from the C2DB database and discover several high Curie temperature FMs, with Fe2F2 and MoI2 emerging as the most promising 2D ferromagnets. Finally, by comparing to experimental results for CrI3, CrCl3, and CrBr3, we conclude that for small effective anisotropies, the Green's-function-based equations are preferable, while for larger anisotropies, MC-based results are more predictive.

Abstract: Using first-principles calculations, we investigate the magnetic order in two-dimensional (2D) transition-metal-dichalcogenide (TMD) monolayers: MoS2, MoSe2, MoTe2, WSe2, and WS2 substitutionally doped with period four transition-metals (Ti, V, Cr, Mn, Fe, Co, Ni). We uncover five distinct magnetically ordered states among the 35 distinct TMD-dopant pairs: the non-magnetic (NM), the ferromagnetic with out-of-plane spin polarization (Z FM), the out-of-plane polarized clustered FMs (clustered Z FM), the in-plane polarized FMs (X-Y FM), and the anti-ferromagnetic (AFM) state. Ni and Ti dopants result in an NM state for all considered TMDs, while Cr dopants result in an anti-ferromagnetically ordered state for all the TMDs. Most remarkably, we find that Fe, Mn, Co, and V result in an FM ordered state for all the TMDs, except for MoTe2. Finally, we show that V-doped MoSe2 and WSe2, and Mn-doped MoS2, are the most suitable candidates for realizing a room-temperature FM at a 16-18% atomic substitution.

Abstract: Transition metal dichalcogenides, intercalated with transition metals, are studied for their potential applications as dilute magnetic semiconductors. We investigate the magnetic properties of WSe2 doped with third-row transition metals (Co, Cr, Fe, Mn, Ti and V). Using density functional theory in combination with Monte Carlo simulations, we obtain an estimate of the Curie or Neel temperature. We find that the magnetic ordering is highly dependent on the dopant type. While Ti and Cr-doped WSe2 have a ferromagnetic ground state, V, Mn, Fe and Co-doped WSe2 are antiferromagnetic in their ground state. For Fe doped WSe2, we find a high Curie-temperature of 327 K. In the case of V-doped WSe2, we find that there are two distinct magnetic phase transitions, originating from a frustrated in-plane antiferromagnetic exchange interaction and a ferromagnetic out-of-plane interaction. We calculate the formation energy and reveal that, in contrast to earlier reports, the formation energy is positive for the intercalated systems studied here. We also show that in the presence of W-vacancies, it becomes favorable for Ti, Fe, and Co to intercalate in WSe2.

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: In tunnel field-effect transistors, trap-assisted tunneling (TAT) is one of the probable causes for degraded subthreshold swing. The accurate quantum-mechanical (QM) assessment of TAT currents also requires a QM treatment of phonon-assisted tunneling (PAT) currents. Therefore, we present a multi-band PAT current formalism within the framework of the quantum transmitting boundary method. An envelope function approximation is used to construct the electron-phonon coupling terms corresponding to local Frohlich-based phonon-assisted inter-band tunneling in direct-bandgap III-V semiconductors. The PAT current density is studied in up to 100 nm long and 20 nm wide p-n diodes with the 2- and 15-band material description of our formalism. We observe an inefficient electron-phonon coupling across the tunneling junction. We further demonstrate the dependence of PAT currents on the device length, for our non-self-consistent formalism which neglects changes in the electron distribution function caused by the electron-phonon coupling. Finally, we discuss the differences in doping dependence between direct band-to-band tunneling and PAT current. Published under license by AIP Publishing.

Abstract: In the quantum mechanical simulation of exploratory semiconductor devices, continuum methods based on a k.p/envelope function model have the potential to significantly reduce the computational burden compared to prevalent atomistic methods. However, full-zone k.p/envelope function simulation approaches are scarce and existing implementations are not self-consistent with the calculation of the electrostatic potential due to the lack of a stable procedure and a proper normalization of the multi-band envelope functions. Here, we therefore present a self-consistent procedure based on a full-zone spectral k.p/envelope function band structure model. First, we develop a proper normalization for the multi-band envelope functions in the presence of transmitting boundary conditions. This enables the calculation of the free carrier densities. Next, we construct a procedure to obtain self-consistency of the carrier densities with the electrostatic potential. This procedure is stabilized with an adaptive scheme that relies on the solution of Poisson's equation in the Gummel form, combined with successive underrelaxation. Finally, we apply our procedure to homostructure In0.53Ga0.47As tunnel field-effect transistors (TFETs) and staggered heterostructure GaAs0.5Sb0.5/In0.53Ga0.47As TFETs and show the importance of self-consistency on the device predictions for scaled dimensions. Published by AIP Publishing.

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 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: Promising predictions are made for III-V tunnel-field-effect transistor (FET), but there is still uncertainty on the parameters used in the band-to-band tunneling models. Therefore, two simulators are calibrated in this paper; the first one uses a semi-classical tunneling model based on Kane's formalism, and the second one is a quantum mechanical simulator implemented with an envelope function formalism. The calibration is done for In0.53Ga0.47As using several p+/intrinsic/n+ diodes with different intrinsic region thicknesses. The dopant profile is determined by SIMS and capacitance-voltage measurements. Error bars are used based on statistical and systematic uncertainties in the measurement techniques. The obtained parameters are in close agreement with theoretically predicted values and validate the semi-classical and quantum mechanical models. Finally, the models are applied to predict the input characteristics of In0.53Ga0.47As n- and p-lineTFET, with the n-lineTFET showing competitive performance compared to MOSFET.

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.