Search

Transition metal dichalcogenides (TMDs) are promising candidates for a wide variety of ultrascaled electronic, quantum computation, and optoelectronic applications. The exponential decay of electronic density of states into the bandgap, i.e. the band tail has a strong impact on the performance of TMD applications. In this work, the band tails of various TMD monolayer and multilayer systems when placed on various dielectric substrates is predicted with density functional theory based nonequilibrium Green's functions. Nonlocal scattering of electrons on polar optical phonons, charged impurities and remote scattering on phonons in the dielectric materials is included in the self-consistent Born approximation. The band tails are found to critically depend on the layer thickness, temperature, doping concentration and particularly on the chosen dielectric substrate. The underlying physical mechanisms are studied in high detail and an analytical interpolation formula is given to provide a quick-reference for Urbach parameters in $MoS_2$, $WS_2$ and $WSe_2$.

Short wavelength ultraviolet (UV-C) light deactivates DNA of any germs, including multiresistive bacteria and viruses like COVID-19. Two-dimensional (2D) material-based UV-C light emitting diodes can potentially be integrated into arbitrary surfaces to allow for shadow-free surface disinfection. In this work, we perform a series of first-principles calculations to identify the core components of ultrathin LEDs based on hexagonal boron nitride (hBN). The electrons and holes are predicted to be confined in multiple quantum wells (MQWs) by combining hBN layers with different stacking orders. Various p- and n-doping candidates for hBN are assessed, and the relative p- and n-type metal contacts with low Schottky barrier heights are identified. The findings are summarized in a concrete UV-C LED structure proposal., Comment: 20 pages, 8 figures

The nonequilibrium Green's function (NEGF) method is often used to predict transport in atomistically resolved nanodevices and yields an immense numerical load when inelastic scattering on phonons is included. To ease this load, this work extends the atomistic mode space approach of Ref. [1] to include inelastic scattering on optical and acoustic phonons in silicon nanowires. This work also includes the exact calculation of the real part of retarded scattering self-energies in the reduced basis representation using the Kramers-Kronig relations. The inclusion of the Kramers-Kronig relation for the real part of the retarded scattering self-energy increases the impact of scattering. Virtually perfect agreement with results of the original representation is achieved with matrix rank reductions of more than 97%. Time-to-solution improvements of more than 200$\times$ and peak memory reductions of more than 7$\times$ are shown. This allows for the solution of electron transport scattered on phonons in atomically resolved nanowires with cross-sections larger than 5 nm $\times$ 5 nm., Comment: 11 pages, 10 figures

State-of-the-art industrial semiconductor device modeling is based on highly efficient Drift-Diffusion (DD) models that include some quantum corrections for nanodevices. In contrast, latest academic quantum transport models are based on the non-equilibrium Green's function (NEGF) method that cover all coherent and incoherent quantum effects consistently. Carrier recombination and generation in optoelectronic nanodevices represent an immense numerical challenge when solved within NEGF. In this work, the numerically efficient B\"uttiker-probe model is expanded to include electron-hole recombination and generation in the NEGF framework. Benchmarks of the new multiple-particle B\"{u}ttiker probe method against state-of-the-art quantum-corrected DD models show quantitative agreements except in cases of pronounced tunneling and interference effects.

Simulations are essential to accelerate the discovery of new materials and to gain full understanding of known ones. Although hard to realize experimentally, periodic boundary conditions are omnipresent in material simulations. In this work, we intro-duce ROBIN (recursive open boundary and interfaces), the first method allowing open boundary conditions in material and interface modeling. The computational costs are limited to solving quantum properties in a focus area which allows explicitly discretizing millions of atoms in real space and to consider virtually any type of environment (be it periodic, regular, or ran-dom). The impact of the periodicity assumption is assessed in detail with silicon dopants in graphene. Graphene was con-firmed to produce a band gap with periodic substitution of 3% carbon with silicon in agreement with published periodic boundary condition calculations. Instead, 3% randomly distributed silicon in graphene only shifts the energy spectrum. The predicted shift agrees quantitatively with published experimental data. Key insight of this assessment is, assuming periodici-ty elevates a small perturbation of a periodic cell into a strong impact on the material property prediction. Periodic boundary conditions can be applied on truly periodic systems only. More general systems should apply an open boundary method for reliable predictions.

The non-equilibrium Green's function (NEGF) method with B\"uttiker probe scattering self-energies is assessed by comparing its predictions for the thermal boundary resistance with molecular dynamics (MD) simulations. For simplicity, the interface of Si/heavy-Si is considered, where heavy-Si differs from Si only in the mass value. With B\"uttiker probe scattering parameters tuned against MD in homogeneous Si, the NEGF-predicted thermal boundary resistance quantitatively agrees with MD for wide mass ratios. Artificial resistances that the unaltered Landauer approach yield at virtual interfaces in homogeneous systems are absent in the present NEGF approach. Spectral information result from NEGF in its natural representation without further transformations. The spectral results show that the scattering between different phonon modes plays a crucial role in thermal transport across interfaces. B\"uttiker probes provide an efficient and reliable way to include anharmonicity in phonon related NEGF. NEGF including the B\"uttiker probes can reliably predict phonon transport across interfaces and at finite temperatures., Comment: 4 pages, 5 figures

High-doping induced Urbach tails and band gap narrowing play a significant role in determining the performance of tunneling devices and optoelectronic devices such as tunnel field-effect transistors (TFETs), Esaki diodes and light-emitting diodes. In this work, Urbach tails and band gap narrowing values are calculated explicitly for GaAs, InAs, GaSb and GaN as well as ultra-thin bodies and nanowires of the same. Electrons are solved in the non-equilibrium Green's function method in multi-band atomistic tight binding. Scattering on polar optical phonons and charged impurities is solved in the self-consistent Born approximation. The corresponding nonlocal scattering self-energies as well as their numerically efficient formulations are introduced for ultra-thin bodies and nanowires. Predicted Urbach band tails and conduction band gap narrowing agree well with experimental literature for a range of temperatures and doping concentrations. Polynomial fits of the Urbach tail and band gap narrowing as a function of doping are tabulated for quick reference.

State of the art quantum transport models for semiconductor nanodevices attribute negative (positive) unit charges to states of the conduction (valence) band. Hybrid states that enable band-to-band tunneling are subject to interpolation that yield model dependent charge contributions. In any nanodevice structure, these models rely on device and physics specific input for the dielectric constants. This paper exemplifies the large variability of different charge interpretation models when applied to ultrathin body transistor performance predictions. To solve this modeling challenge, an electron-only band structure model is extended to atomistic quantum transport. Performance predictions of MOSFETs and tunneling FETs confirm the generality of the new model and its independence of additional screening models., Comment: 8 pages, 8 figures

It is assessed in detail both experimentally and theoretically how the interlayer coupling of transition metal dichalcogenides controls the electronic properties of the respective devices. Gated transition metal dichalcogenide structures show electrons and holes to either localize in individual monolayers, or delocalize beyond multiple layers - depending on the balance between spin-orbit inter- action and interlayer hopping. This balance depends on layer thickness, momentum space symmetry points and applied gate fields. A good quantitative agreement of predictions and measurements of the quantum confined Stark effect in gated MoS2 systems unveils intralayer excitons as major source for the observed photoluminesence., Comment: 11 pages, 15 figures

In this work, we use a combination of first-principles calculations under the density functional theory framework and heat transport simulations using the atomistic Green's function (AGF) method to quantitatively predict the contribution of the different scattering mechanisms to the thermal interface conductance of epitaxial CoSi$_2$-Si interfaces. An important development in the present work is the direct computation of interfacial bonding from density functional perturbation theory (DFPT) and hence the avoidance of commonly used `mixing rules' to obtain the cross-interface force constants from bulk material force constants. Another important algorithmic development is the integration of the recursive Green's function (RGF) method with B\"{u}ttiker probe scattering that enables computationally efficient simulations of inelastic phonon scattering and its contribution to the thermal interface conductance. First-principles calculations of electron-phonon coupling reveal that cross-interface energy transfer between metal electrons and atomic vibrations in the semiconductor is mediated by delocalized acoustic phonon modes that extend on both sides of the interface, and phonon modes that are localized inside the semiconductor region of the interface exhibit negligible coupling with electrons in the metal. We also provide a direct comparison between simulation predictions and experimental measurements of thermal interface conductance of epitaxial CoSi$_2$-Si interfaces using the time-domain thermoreflectance technique. Importantly, the experimental results, performed across a wide temperature range, only agree well with predictions that include all transport processes: elastic and inelastic phonon scattering, electron-phonon coupling in the metal, and electron-phonon coupling across the interface.

It is critical to capture the effect due to strain and material interface for device level transistor modeling. We introduced a transferable sp3d5s* tight binding model with nearest neighbor interactions for arbitrarily strained group IV and III-V materials. The tight binding model is parameterized with respect to Hybrid functional(HSE06) calculations for varieties of strained systems. The tight binding calculations of ultra small superlattices formed by group IV and group III-V materials show good agreement with the corresponding HSE06 calculations. The application of tight binding model to superlattices demonstrates that transferable tight binding model with nearest neighbor interactions can be obtained for group IV and III-V materials.

Electronic performance predictions of modern nanotransistors require nonequilibrium Green's functions including incoherent scattering on phonons as well as inclusion of random alloy disorder and surface roughness effects. The solution of all these effects is numerically extremely expensive and has to be done on the world's largest supercomputers due to the large memory requirement and the high performance demands on the communication network between the compute nodes. In this work, it is shown that NEMO5 covers all required physical effects and their combination. Furthermore, it is also shown that NEMO5's implementation of the algorithm scales very well up to about 178176CPUs with a sustained performance of about 857 TFLOPS. Therefore, NEMO5 is ready to simulate future nanotransistors.

Empirical Tight Binding (TB) methods are widely used in atomistic device simulations. Existing TB methods to passivate dangling bonds fall into two categories: 1) Method that explicitly includes passivation atoms is limited to passivation with atoms and small molecules only. 2) Method that implicitly incorporates passivation does not distinguish passivation atom types. This work introduces an implicit passivation method that is applicable to any passivation scenario with appropriate parameters. This method is applied to a Si quantum well and a Si ultra-thin body transistor oxidized with SiO2 in several oxidation configurations. Comparison with ab-initio results and experiments verifies the presented method. Oxidation configurations that severely hamper the transistor performance are identified. It is also shown that the commonly used implicit H atom passivation overestimates the transistor performance., Comment: 6 pages, 4 figures

The bandstructures of [110] and [001] Bi2Te3 nanowires are solved with the atomistic 20 band tight binding functionality of NEMO5. The theoretical results reveal: The popular assumption that all topological insulator wire surfaces are equivalent is inappropriate. The Fermi velocity of chemically distinct wire surfaces differs significantly which creates an effective in-surface confinement potential. As a result, topological insulator surface states prefer specific surfaces. Therefore, experiments have to be designed carefully not to probe surfaces unfavorable to the surface states (low density of states) and thereby be insensitive to the TI-effects., Comment: 2015 Applied Physics Letters

In this work, transferable empirical tight binding parameters of strained group IV and III-V semiconductors are generated from ab-initio calculations. The empirical tight binding parameters show good transferability when applied to strained bulk materials as well as ultra-thin superlattices.

Empirical tight binding(ETB) methods are widely used in atomistic device simulations. Traditional ways of generating the ETB parameters rely on direct fitting to bulk experiments or theoretical electronic bands. However, ETB calculations based on existing parameters lead to unphysical results in ultra small structures like the As terminated GaAs ultra thin bodies(UTBs). In this work, it is shown that more reliable parameterizations can be obtained by a process of mapping ab-initio bands and wave functions to tight binding models. This process enables the calibration of not only the ETB energy bands but also the ETB wave functions with corresponding ab-initio calculations. Based on the mapping process, ETB model of Si and GaAs are parameterized with respect to hybrid functional calculations. Highly localized ETB basis functions are obtained. Both the ETB energy bands and wave functions with subatomic resolution of UTBs show good agreement with the corresponding hybrid functional calculations. The ETB methods can then be used to explain realistically extended devices in non-equilibrium that can not be tackled with ab-initio methods., Comment: 11pages, 10 figures

Quantum cascade lasers are unipolar semiconductor lasers covering a wide range of the infrared and terahertz spectrum. Lasing action is achieved by using optical intersubband transitions between quantized states in specifically designed multiple-quantum-well heterostructures. A systematic improvement of quantum cascade lasers with respect to operating temperature, efficiency and spectral range requires detailed modeling of the underlying physical processes in these structures. Moreover, the quantum cascade laser constitutes a versatile model device for the development and improvement of simulation techniques in nano- and optoelectronics. This review provides a comprehensive survey and discussion of the modeling techniques used for the simulation of quantum cascade lasers. The main focus is on the modeling of carrier transport in the nanostructured gain medium, while the simulation of the optical cavity is covered at a more basic level. Specifically, the transfer matrix and finite difference methods for solving the one-dimensional Schr\"odinger equation and Schr\"odinger-Poisson system are discussed, providing the quantized states in the multiple-quantum-well active region. The modeling of the optical cavity is covered with a focus on basic waveguide resonator structures. Furthermore, various carrier transport simulation methods are discussed, ranging from basic empirical approaches to advanced self-consistent techniques. The methods include empirical rate equation and related Maxwell-Bloch equation approaches, self-consistent rate equation and ensemble Monte Carlo methods, as well as quantum transport approaches, in particular the density matrix and non-equilibrium Green's function (NEGF) formalism. The derived scattering rates and self-energies are generally valid for n-type devices based on one-dimensional quantum confinement, such as quantum well structures.

Bi$_{2}$Te$_{3}$ and Bi$_{2}$Se$_{3}$ are well known 3D-topological insulators. Films made of these materials exhibit metal-like surface states with a Dirac dispersion and possess high mobility. The high mobility metal-like surface states can serve as channel material for TI-based field effect transistors. While such a transistor offers superior terminal characteristics, they suffer from an inherent zero band gap problem. The absence of a band gap for the surface states prevents an easy turn-off mechanism. In this work, techniques that can be employed to easily open a band gap for the TI surface states is introduced. Two approaches are described: 1) Coating the surface states with a ferromagnet which has a controllable magnetization axis. The magnetization strength of the ferromagnet is incorporated as an exchange interaction term in the Hamiltonian. 2) An \textit{s}-wave superconductor, because of the proximity effect, when coupled to a 3D-TI opens a band gap on the surface. This TI-superconductor heterostructure is modeled using the Bogoliubov-de Gennes Hamiltonian. A comparison demonstrating the finite size effects on surface states of a 3D-TI and edge states of a CdTe/HgTe/CdTe-based 2D-TI is also presented. 3D-TI nanostructures can be reduced to dimensions as low as 10.0 $ \mathrm{nm} $ in contrast to 2D-TI structures which require a thickness of at least 100.0 $ \mathrm{nm} $. All calculations are performed using the continuum four-band k.p Hamiltonian., Comment: 9 pages,12 figures

Semi-Empirical Tight Binding (TB) is known to be a scalable and accurate atomistic representation for electron transport for realistically extended nano-scaled semiconductor devices that might contain millions of atoms. In this paper an environment-aware and transferable TB model suitable for electronic structure and transport simulations in technologically relevant metals, metallic alloys, metal nanostructures and metallic interface systems is described. Part I of this paper describes the development and validation of the new TB model. The new model incorporates intra-atomic diagonal and off-diagonal elements for implicit self-consistency and greater transferability across bonding environments. The dependence of the on-site energies on strain has been obtained by appealing to the Moments Theorem that links closed electron paths in the system to energy moments of angular momentum resolved local density of states obtained ab-initio. The model matches self-consistent DFT electronic structure results for bulk FCC metals with and without strain, metallic alloys, metallic interfaces and metallic nanostructures with high accuracy and can be used in predictive electronic structure and transport problems in metallic systems at realistically extended length scales, Comment: 34 pages including appendices, parameters and supplemental material in double spaced format

The Semi-Empirical TB model developed in part I is applied to metal transport problems of current relevance in part II. A systematic study of the effect of quantum confinement, transport orientation and homogeneous strain on electronic transport properties of Cu is carried out. It is found that quantum confinement from bulk to nanowire boundary conditions leads to significant anisotropy in conductance of Cu along different transport orientations. Compressive homogeneous strain is found to reduce resistivity by increasing the density of conducting modes in Cu. The [110] transport orientation in Cu nanowires is found to be the most favorable for mitigating conductivity degradation since it shows least reduction in conductance with confinement and responds most favorably to compressive strain., Comment: 13 pages with 8 figures in double spaced manuscript format

In this work, the low rank approximation concept is extended to the non-equilibrium Green's function (NEGF) method to achieve a very efficient approximated algorithm for coherent and incoherent electron transport. This new method is applied to inelastic transport in various semiconductor nanodevices. Detailed benchmarks with exact NEGF solutions show 1) a very good agreement between approximated and exact NEGF results, 2) a significant reduction of the required memory, and 3) a large reduction of the computational time (a factor of speed up as high as 150 times is observed). A non-recursive solution of the inelastic NEGF transport equations of a 1000 nm long resistor on standard hardware illustrates nicely the capability of this new method., Comment: 9 pages, 7 figures

As metal-oxide-semiconductor field-effect transistors (MOSFET) channel lengths (Lg) are scaled to lengths shorter than Lg<8 nm source-drain tunneling starts to become a major performance limiting factor. In this scenario a heavier transport mass can be used to limit source-drain (S-D) tunneling. Taking InAs and Si as examples, it is shown that different heavier transport masses can be engineered using strain and crystal orientation engineering. Full-band extended device atomistic quantum transport simulations are performed for nanowire MOSFETs at Lg<8 nm in both ballistic and incoherent scattering regimes. In conclusion, a heavier transport mass can indeed be advantageous in improving ON state currents in ultra scaled nanowire MOSFETs., Comment: 6 pages, 7 figures, journal

The topological insulator properties of CdTe/HgTe/CdTe quantum wells are theoretically studied. The CdTe/HgTe/CdTe quantum well behaves as a topological insulator beyond a critical well width dimension. It is shown that if the barrier(CdTe) and well-region(HgTe) are altered by replacing them with the alloy Cd$_{x}% $Hg$_{1-x}$Te of various stoichiometries, the critical width can be changed.The critical quantum well width is shown to depend on temperature, applied stress, growth directions and external electric fields. Based on these results, a novel device concept is proposed that allows to switch between a normal semiconducting and topological insulator state through application of moderate external electric fields., Comment: 8 pages, 13 figures

The Empirical Tight Binding(ETB) method is widely used in atomistic device simulations. The reliability of such simulations depends very strongly on the choice of basis sets and the ETB parameters. The traditional way of obtaining the ETB parameters is by fitting to experiment data, or critical theoretical bandedges and symmetries rather than a foundational mapping. A further shortcoming of traditional ETB is the lack of an explicit basis. In this work, a DFT mapping process which constructs TB parameters and explicit basis from DFT calculations is developed. The method is applied to two materials: GaAs and MgO. Compared with the existing TB parameters, the GaAs parameters by DFT mapping show better agreement with the DFT results in bulk band structure calculations and lead to different indirect valleys when applied to nanowire calculations. The MgO TB parameters and TB basis functions are also obtained through the DFT mapping process., Comment: single column, 9 pages, 5 figures