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: We propose a simple quantum mechanical model describing the time dependent diffusion current between two fermion reservoirs that were initially disconnected and characterized by different densities or chemical potentials. The exact, analytical solution of the model yields the transient behavior of the coupled fermion systems evolving to a final steady state, whereas the long-time behavior is determined by a power law rather than by exponential decay. Similar results are obtained for the entropy production which is proportional to the diffusion current. (C) 2014 Elsevier B.V. All rights reserved.

Abstract: In the realm of Ehrenfests theorem, classical trajectories obeying Newtons laws have been proven useful to construct explicit solutions to the time-dependent WignerLiouville equation. Whereas previous works have particularly focused on the initial distribution function as a vehicle found to carry the signatures of quantum statistics into the time-dependent solution, the present paper shows that the LagrangeCharpit method based on classical trajectories can be successfully invoked as well to tackle quantum mechanical features with no classical counterpart, such as tunneling.

Abstract: 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. (C) 2017 Elsevier Inc. All rights reserved.

Abstract: For many decades the Boltzmann distribution function has been used to calculate the non-equilibrium properties of mobile particles undergoing the combined action of various scattering mechanisms and externally applied force fields. When the latter give rise to the occurrence of inhomogeneous potential profiles across the region through which the particles are moving, the numerical solution of the Boltzmann equation becomes a highly complicated task. In this work we highlight a particular algorithm that can be used to solve the time dependent Boltzmann equation as well as its quantum mechanical extension, the WignerBoltzmann equation. As an illustration, we show the calculated distribution function describing electrons propagating under the action of both a uniform and a pronouncedly non-uniform electric field.

Abstract: Since Ehrenfests theorem, the role and importance of classical paths in quantum dynamics have been examined by several means. Along this line, we show that the classical equations of motion provide a solution to quantum dynamics, if appropriately incorporated into the Wigner distribution function, exactly reformulated in a type of Boltzmann equation. Also the quantum-mechanical features of the canonical ensemble can be studied in this framework of Newtonian dynamics, if the initial distribution function is appropriately constructed from the statistical operator.

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: 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: Ando's model provides a rigorous quantum-mechanical framework for electron-surface roughness scattering, based on the detailed roughness structure. We apply this method to metallic nanowires and improve the model introducing surface roughness distribution functions on a finite domain with analytical expressions for the average surface roughness matrix elements. This approach is valid for any roughness size and extends beyond the commonly used Prange-Nee approximation. The resistivity scaling is obtained from the self-consistent relaxation time solution of the Boltzmann transport equation and is compared to Prange-Nee's approach and other known methods. The results show that a substantial drop in resistivity can be obtained for certain diameters by achieving a large momentum gap between Fermi level states with positive and negative momentum in the transport direction. (C) 2015 AIP Publishing LLC.

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.

Abstract: A general framework to calculate the Zener current in an indirect semiconductor with an externally applied potential is provided. Assuming a parabolic valence and conduction band dispersion, the semiconductor is in equilibrium in the presence of the external field as long as the electron-phonon interaction is absent. The linear response to the electron-phonon interaction results in a non-equilibrium system. The Zener tunneling current is calculated from the number of electrons making the transition from valence to conduction band per unit time. A convenient expression based on the single particle spectral functions is provided, enabling the evaluation of the Zener tunneling current under any three-dimensional potential profile. For a one-dimensional potential profile an analytical expression is obtained for the current in a bulk semiconductor, a semiconductor under uniform field, and a semiconductor under a non-uniform field using the WKB (Wentzel-Kramers-Brillouin) approximation. The obtained results agree with the Kane result in the low field limit. A numerical example for abrupt p-n diodes with different doping concentrations is given, from which it can be seen that the uniform field model is a better approximation than the WKB model, but a direct numerical treatment is required for low bias conditions.

Abstract: Recently, a renewed interest in Zener tunneling has arisen because of its increasing impact on semiconductor device performance at nanometer dimensions. In this paper we evaluate the tunnel probability under the action of a nonuniform electric field using a two-band model and arrive at significant deviations from the commonly used Kanes model, valid for weak uniform fields only. A threshold on the junction bias where Kanes model for Zener tunneling breaks down is determined. Comparison with Kanes model particularly shows that our calculation yields a higher tunnel probability for intermediate electric fields and a lower tunnel probability for high electric fields. When performing a current calculation comparing to the WKB approximation for the case of an abrupt p-n junction significant differences concerning the shape of the I-V curve are demonstrated.

Abstract: We present comprehensive calculations of the low-field hole mobility in Ge p-channel inversion layers with SiO2 insulator using a six-band k·p band-structure model. The cases of relaxed, biaxially, and uniaxially (both tensily and compressively) strained Ge are studied employing an efficient self-consistent methodmaking use of a nonuniform spatial mesh and of the Broyden second methodto solve the coupled envelope-wave function k·p and Poisson equations. The hole mobility is computed using the KuboGreenwood formalism accounting for nonpolar hole-phonon scattering and scattering with interfacial roughness. Different approximations to handle dielectric screening are also investigated. As our main result, we find a large enhancement (up to a factor of 10 with respect to Si) of the mobility in the case of uniaxial compressive stress similarly to the well-known case of Si. Comparison with experimental data shows overall qualitative agreement but with significant deviations due mainly to the unknown morphology of the rough Ge-insulator interface, to additional scattering with surface optical phonon from the high- insulator, to Coulomb scattering interface traps or oxide chargesignored in our calculationsand to different channel structures employed.