Nonequilibrium Green Function Research Articles

Enhanced Seebeck effect in gated O-shaped Kekulé graphene junctions

The Seebeck effect in graphene with uniform Kekulé lattice distortion has been studied using the tight-binding model combined with the nonequilibrium Green's function method. The numerical results of this work show that the Seebeck coefficient of the O-shaped Kekulé graphene can be enhanced by modulating the gate voltage. At room temperature (T≈300 K), the Seebeck coefficient for the armchair nanoribbon is Sc≈1370μV/K and for the zigzag nanoribbon is Sc≈1260μV/K. At low temperature (T≈128 K), the Seebeck coefficient for the armchair nanoribbon is Sc≈2500μV/K. In addition, the nano-structuration of the O-shaped Kekulé graphene layer induces an anisotropic behavior in the Seebeck effect. Our work provides an idea for the fabrication of efficient and controllable thermoelectric effect devices.

Impact of Rh, Ru, and Pd Leads and Contact Topologies on Performance of WSe2 FETs: A First Comparative Ab Initio Study.

2D field-effect transistors (FETs) fabricated with transition metal dichalcogenide (TMD) materials are a potential replacement for the silicon-based CMOS. However, the lack of advancement in p-type contact is also a key factor hindering TMD-based CMOS applications. The less investigated path towards improving electrical characteristics based on contact geometries with low contact resistance (RC) has also been established. Moreover, finding contact metals to reduce the RC is indeed one of the significant challenges in achieving the above goal. Our research provides the first comparative analysis of the three contact configurations for a WSe2 monolayer with different noble metals (Rh, Ru, and Pd) by employing ab initio density functional theory (DFT) and non-equilibrium Green's function (NEGF) methods. From the perspective of the contact topologies, the RC and minimum subthreshold slope (SSMIN) of all the conventional edge contacts are outperformed by the novel non-van der Waals (vdW) sandwich contacts. These non-vdW sandwich contacts reveal that their RC values are below 50 Ω∙μm, attributed to the narrow Schottky barrier widths (SBWs) and low Schottky barrier heights (SBHs). Not only are the RC values dramatically reduced by such novel contacts, but the SSMIN values are lower than 68 mV/dec. The new proposal offers the lowest RC and SSMIN, irrespective of the contact metals. Further considering the metal leads, the WSe2/Rh FETs based on the non-vdW sandwich contacts show a meager RC value of 33 Ω∙μm and an exceptional SSMIN of 63 mV/dec. The two calculated results present the smallest-ever values reported in our study, indicating that the non-vdW sandwich contacts with Rh leads can attain the best-case scenario. In contrast, the symmetric convex edge contacts with Pd leads cause the worst-case degradation, yielding an RC value of 213 Ω∙μm and an SSMIN value of 95 mV/dec. While all the WSe2/Ru FETs exhibit medium performances, the minimal shift in the transfer curves is interestingly advantageous to the circuit operation. Conclusively, the low-RC performances and the desirable SSMIN values are a combination of the contact geometries and metal leads. This innovation, achieved through noble metal leads in conjunction with the novel contact configurations, paves the way for a TMD-based CMOS with ultra-low RC and rapid switching speeds.

Performance Projection of Vacuum Gate Dielectric Doping-Free Carbon Nanoribbon/Nanotube Field-Effect Transistors for Radiation-Immune Nanoelectronics.

This paper investigates the performance of vacuum gate dielectric doping-free carbon nanotube/nanoribbon field-effect transistors (VGD-DL CNT/GNRFETs) via computational analysis employing a quantum simulation approach. The methodology integrates the self-consistent solution of the Poisson solver with the mode space non-equilibrium Green's function (NEGF) in the ballistic limit. Adopting the vacuum gate dielectric (VGD) paradigm ensures radiation-hardened functionality while avoiding radiation-induced trapped charge mechanisms, while the doping-free paradigm facilitates fabrication flexibility by avoiding the realization of a sharp doping gradient in the nanoscale regime. Electrostatic doping of the nanodevices is achieved via source and drain doping gates. The simulations encompass MOSFET and tunnel FET (TFET) modes. The numerical investigation comprehensively examines potential distribution, transfer characteristics, subthreshold swing, leakage current, on-state current, current ratio, and scaling capability. Results demonstrate the robustness of vacuum nanodevices for high-performance, radiation-hardened switching applications. Furthermore, a proposal for extrinsic enhancement via doping gate voltage adjustment to optimize band diagrams and improve switching performance at ultra-scaled regimes is successfully presented. These findings underscore the potential of vacuum gate dielectric carbon-based nanotransistors for ultrascaled, high-performance, energy-efficient, and radiation-immune nanoelectronics.

Theoretical study on the structure, spectroscopic, and current–voltage behavior of 11-Cis and Trans retinal isomers in rhodopsin

In this study, the electronic transport properties of 11-Cis and Trans retinal, components of rhodopsin, were investigated as optical molecular switches using the nonequilibrium Green’s function (NEGF) formalism combined with first-principles density functional theory (DFT). These isomers, which can be reversibly converted into each other, were examined in detail. The structural and spectroscopic properties, including infrared (IR), Raman, nuclear magnetic resonance (NMR), and ultraviolet (UV) spectra, were analyzed using the hybrid B3LYP/6–311 + + G** level of theory. Complete vibrational assignments were performed for both forms utilizing the scaled quantum mechanical force field (SQMFF) methodology. To evaluate the conductivity of these molecules, we utilized current–voltage (I-V) characteristics, transmission spectra, molecular projected self-consistent Hamiltonian (MPSH), HOMO–LUMO gap, and second-order interaction energies (E2). The trendline extrapolation of the current–voltage plots confirmed our findings. We investigated the effect of different electrodes (Ag, Au, Pt) and various connection sites (hollow, top, bridge) on conductivity. The Ag electrode with the hollow site exhibited the highest efficiency. Our results indicate that the Cis form has higher conductivity than the Trans form.

Spin-Orbit Torque in Single-Molecule Junctions from ab Initio.

The use of electric fields applied across magnetic heterojunctions that lack spatial inversion symmetry has been previously proposed as a nonmagnetic means of controlling localized magnetic moments through spin-orbit torques (SOT). The implementation of this concept at the single-molecule level has remained a challenge, however. Here, we present first-principles calculations of SOT in a single-molecule junction under bias and beyond linear response. Employing a self-consistency scheme invoking density functional theory and nonequilibrium Green's function theory including spin-orbit interaction, we compute the change of the magnetization with the bias voltage and the associated current-induced SOT. Within the linear regime our quantitative estimates for the SOT in single-molecule junctions yield values similar to those known for magnetic interfaces. Our findings contribute to an improved microscopic understanding of SOT in single molecules.

Time‐Resolved Optical Response of the Dicke's Model via the Nonequilibrium Green's Function Approach

Due to their conceptual appeal and computational convenience, two‐level systems (TLSs) and generalizations are used to investigate nonlinear behavior in quantum optics and assess the applicability of theoretical methods. Herein, the focus is on second‐harmonic generation (SHG) and, as system of interest, on the Dicke model, which consists of several TLSs inside an optical cavity. The main aspect addressed is the scope of nonequilibrium Green's function (NEGF) to describe the effect of inhomogeneities and electron–electron (e–e) interactions on the SHG signal. For benchmarking purposes, exact diagonalization (ED) results are also presented and discussed. SHG spectra obtained with NEGF and ED are found to be in very good mutual agreement in most situations. Furthermore, inhomogeneity in the TLS and e–e interactions reduce the SHG signal, and the reduction is stronger with inhomogeneity than with interactions. This trend is consistently noted across different (small to large) system sizes. Finally, a modified NEGF approach is proposed to account for cavity leakage, where the quantum photon fields are coupled to a bath of classical oscillators. As expected, within this mixed quantum‐classical scheme, a decrease in the intensity of the fluorescent spectra takes place depending on the entity of leakage.

Steric Effects on Single-Molecule Conductance in Flat-Lying Phenanthrene.

A previous combined experimental and theoretical study found that the position of anchoring groups on a phenanthrene (PHE) backbone played a large role in determining the single-molecule conductance of the PHE derivative. However, a consistent 0.1 G0 feature was found across all PHE derivatives. To understand this, the previously investigated PHE derivatives were placed flat on a simulated Au substrate with a scanning tunneling microscope (STM) tip over PHE and conductance was calculated using the non-equilibrium Green's function technique in conjunction with density functional theory (NEGF-DFT). The location of the tip was varied to find the most conductive and most energetically favorable arrangements, which did not coincide. Furthermore, the variation in conductance found in erect junctions was not present when PHE derivatives were lying flat, with all derivatives calculated to have conductance values around 0.1 G0.

Jiezi: An open-source Python software for simulating quantum transport based on non-equilibrium Green's function formalism

We present a Python-based open-source library named Jiezi, which provides the means of simulating the electronic transport properties of nanoscaled devices on the atomistic level. The key feature of Jiezi lies in its core algorithm, i.e., self-consistent orchestration between the non-equilibrium Green's function (NEGF) method and a Poisson's equation solver. Beyond the construction of the tight-binding (TB) Hamiltonian with empirical parameters for conventional materials, the package offers a comprehensive framework for constructing the Wannier-based Hamiltonian matrix, enabling the investigation of novel materials and their heterostructures. To expedite the solution of NEGF systems, a methodology based on renormalization theory is proposed for reducing the dimension of the Hamiltonian matrix. Additionally, we adopt a non-linear Poisson equation solver with no analytical approximation in this software. The software facilitates seamless integration with external tools for geometry and mesh generation and post-processing. In this paper, we present the main capabilities and workflow by demonstrating with a simulation for the carbon nanotube field-effect transistor (CNTFET). Program summaryProgram Title: JieziCPC Library link to program files:https://doi.org/10.17632/nk79kbtww4.1Developer's repository link:https://github.com/Jiezi-negf/JieziLicensing provisions: GPLv3Programming language: PythonNature of problem: Simulates the quantum transport property of nano-scaled transistors based on the predefined device structure and the material composition.Solution method: Solves the coupled Schrödinger equation and Poisson equation by NEGF and finite element method.

Molecular heat transport across a time-periodic temperature gradient.

The time-periodic modulation of a temperature gradient can alter the heat transport properties of a physical system. Oscillating thermal gradients give rise to behaviors such as modified thermal conductivity and controllable time-delayed energy storage that are not present in a system with static temperatures. Here, we examine how the heat transport properties of a molecular lattice model are affected by an oscillating temperature gradient. We use analytical analysis and molecular dynamics simulations to investigate the vibrational heat flow in a molecular lattice system consisting of a chain of particles connected to two heat baths at different temperatures, where the temperature difference between baths is oscillating in time. We derive expressions for heat currents in this system using a stochastic energetics framework and a nonequilibrium Green's function approach that is modified to treat the nonstationary average energy fluxes. We find that emergent energy storage, energy release, and thermal conductance mechanisms induced by the temperature oscillations can be controlled by varying the frequency, waveform, and amplitude of the oscillating gradient. The developed theoretical approach provides a general framework to describe how vibrational heat transmission through a molecular lattice is affected by temperature gradient oscillations.

Carrier Transport Properties of the Orthorhombic Phase Boron Nitride Nanoribbons and Rectifying Device Design

Using density functional theory combined with nonequilibrium Green's function, the electronic structures and carrier transport properties of the orthorhombic phase boron nitride nanoribbons (BNNRs) with different edges and different widths are investigated. The calculated results show that both armchair‐ and zigzag‐edged BNNRs are direct bandgap semiconductors. The quantum confinement gives rise to a distinct bandgap for the 1D BNNRs, and the edged atoms play an important role on their electronic structures. The bandgaps of the armchair‐edged BNNRs approach 0.5–0.8 eV when their widths are sufficiently large. However, the electronic structures of the zigzag‐edged BNNRs exhibit more variations dependent on their edged atoms. The J–V characteristics of the BNNRs reveal that the zigzag‐edged BNNRs have superior J–V performance to that of the armchair‐edged counterparts. Eventually, a computational nanodevice design is proposed to achieve rectifying behavior by constructing an in‐plane heterojunction based on the BNNRs without any doping, and the left and right electrodes are formed by differently edged atoms and a small strain is performed on the left electrode. The rectification ratio of about 77.36 is realized. The findings provide comprehensive understandings of the orthorhombic BN monolayer, which can be helpful to its potential application for nanoelectronics.

Numerical investigation on the convergence of self-consistent Schrödinger-Poisson equations in semiconductor device transport simulation

Semiconductor devices at the nanoscale with low-dimensional materials as channels exhibit quantum transport characteristics, thereby their electrical simulation relies on the self-consistent solution of the Schrödinger-Poisson equations. While the non-equilibrium Green’s function (NEGF) method is widely used for solving this quantum many-body problem, its high computational cost and convergence challenges with the Poisson equation significantly limit its applicability. In this study, we investigate the stability of the NEGF method coupled with various forms of the Poisson equation, encompassing linear, analytical nonlinear, and numerical nonlinear forms Our focus lies on simulating carbon nanotube field-effect transistors (CNTFETs) under two distinct doping scenarios: electrostatic doping and ion implantation doping. The numerical experiments reveal that nonlinear formulas outperform linear counterpart. The numerical one demonstrates superior stability, particularly evident under high bias and ion implantation doping conditions. Additionally, we investigate different approaches for presolving potential, leveraging solutions from the Laplace equation and a piecewise guessing method tailored to each doping mode. These methods effectively reduce the number of iterations required for convergence.

The enhancement of infrared characterization of passivated InAs nanowires

Indium arsenide (InAs) nanowires (NWs) show significant advantages in optoelectronic devices due to their excellent electronic and optical properties. While the surface states affect such outstanding properties, and the performance of optoelectronic devices can be further optimized by modulating the surface state. Herein, the adsorption and substitution properties and the optoelectronics response of passivated InAs NWs are studied using first-principles calculations combined with nonequilibrium Green's function (NEGF). The results show that O atom adsorption can reduce the band gap and improve infrared absorption of InAs NWs. Meanwhile, for the O atom substitution model, the impurity energy level appears near the Fermi level of InAs NWs, and for a higher concentration O atoms substitution, the semiconduction properties of the InAs NWs are destroyed. We also choose the Au, Cs, Bi, Sn, and Y atoms as decoration atoms to passivate InAs NWs, and compared with the intrinsic InAs NW surface model, the energy bands gap of the metal-atoms adsorption models all show a decreasing trend in addition to Au-atom adsorption, and the absorption coefficients of InAs NWs adsorbed by several atoms such as Cs, Sn, Bi, and Y in the infrared region have redshifts in different degrees. The absorption coefficients of InAs NWs adsorbed by several metal atoms, Sn, Bi, and Y, are higher in the infrared region compared to that of InAs NW adsorbed by O atom. After the screening of materials, we select InAs NWs adsorbed with Cs, Bi, Sn, and Y atoms with good infrared absorption characteristics for the performance study of optoelectronic devices and it is shown that the InAs NW devices with adsorbed Bi and Sn atoms have the optimal performance with low dark current and high photocurrent.

Exploring Intertwined quantum and cryogenic behaviour in ultra-scaled 10 nm MOSFET: a NEGF quantum ballistic simulation

Silicon-based spin qubits have emerged as promising candidates for scalable quantum information processing. This study first time investigates the behaviour of ultra-scaled 10 nm gate length and 3 nm channel thickness nanoscale double gate metal-oxide semiconductor field-effect transistors (MOSFETs) over a broad temperature range, from deep cryogenic (4 K) to room temperature (300 K). Employing the Non-Equilibrium Green’s Function (NEGF) method, the research explores the intertwined quantum and cryogenic behaviours of the various quantum phenomena, including eigen energies, eigen-functions, electron concentration, current characteristics, and more. This comprehensive analysis sheds light on the intricate interplay of quantum effects in nanoscale transistors under deep cryogenic conditions, offering valuable insights into the development of cryo-CMOS circuits for quantum computing.

Numerical simulation of Rashba spin–orbit coupling effect on quantum Hall resistivity within tight-binding approximation

A numerical simulation of the effect of Rashba spin–orbit coupling (RSO) on quantum transport properties of electrons confined to a two-dimensional electron gas (2DEG) in the integer quantum Hall regime is presented. The mechanism of quantization and the dependence of the Hall resistivity (HR) on the system built-in and external parameters are studied using the non-equilibrium Green’s function (NEGF) technique. Relevant Green’s functions are computed numerically using spectral decomposition method in combination with recursive techniques. The tight-binding-like representation of the system Hamiltonian including magnetic field is obtained within the finite difference method, using Dirichlet boundary conditions. Corresponding formula linking the Hall resistivity with the transmission function of the system has been derived analytically within the NEGF formalism in a general form, in the sense that the formula can also be used beyond the linear response limit and for nonzero temperatures. In particular, it was shown explicitly that in the presence of the RSO coupling the even- and odd plateaus emerging in the magnetic-field dependence of the quantum Hall resistivity are related to different physical mechanisms.

Study of carbon monoxide adsorption on B-N doped hydrogenated Sn nanoribbons towards nanosensor applications

This study utilized Van der Waals corrected density functional theory (DFT) and nonequilibrium Green's function (NEGF) to consider the potential of hydrogenated Armchair stanene nanoribbons of width 8 (HASnNR(8)) for CO gas detection. The study analyzed the most stable adsorption configuration, charge transfer, electronic transport, adsorption energy, sensitivity and electronic properties of CO gas molecules on HASnNR(8), both before and after B-N doping. The results showed that B-N doping led to a significant increase in adsorption energy and a gap opening at the Gamma point, resulting in semiconductor behavior. Analysis of charge transfer indicated that CO molecules acted as donors to HASnNR(8). Moreover, there was a notable change in the transmission functions of CO adsorbed on B-N doped HASnNR(8) compared to pure HASnNR(8). Based on the I-V characteristics, the study concluded that CO molecules could be easily detected after B-N doping. Finally, the study performed a recovery time analysis and sensitivity CO adsorbed on bare and B-N doped HASnNR(8) surfaces, with the results indicating that B-N doped HASnNR is a promising material for CO molecule sensing applications.

Super-Suppression of Long-Wavelength Phonons in Constricted Nanoporous Geometries.

In a typical semiconductor material, the majority of the heat is carried by long-wavelength, long-mean-free-path phonons. Nanostructuring strategies to reduce thermal conductivity, a promising direction in the field of thermoelectrics, place scattering centers of size and spatial separation comparable to the mean free paths of the dominant phonons to selectively scatter them. The resultant thermal conductivity is in most cases well predicted using Matthiessen's rule. In general, however, long-wavelength phonons are not as effectively scattered as the rest of the phonon spectrum. In this work, using large-scale molecular-dynamics simulations, non-equilibrium Green's function simulations, and Monte Carlo simulations, we show that specific nanoporous geometries that create narrow constrictions in the passage of phonons lead to anticorrelated heat currents in the phonon spectrum. This effect results in super-suppression of long-wavelength phonons due to heat trapping and reductions in the thermal conductivity to values well below those predicted by Matthiessen's rule.

Transport properties of Mn-doped/adsorbed zigzag boron nitride nanoribbon based nanodevices

We are conducting research into the spin transport characteristics of models that are built from zigzag Boron Nitride Nanoribbons (ZBNNRs) that have been doped and adsorbed with Mn. These investigations are made possible by using the non-equilibrium Green’s function (NEGF) formalism in conjunction with Density Functional Theory (DFT). We have scrutinized the electronic structures of six different models, and have chosen two for a deeper exploration of its’ transport characteristics. Our computational findings highlight an unusual occurrence where a Mn atom substitutes a B atom in ZBNNRs, resulting in ZB(Mn)NNRs that exhibit significant Giant Magnetoresistance (GMR) and Rectification ratio (RR) - peaking at 106 and 105 respectively, along with the presence of negative differential conductance (NDC) effects. The detected GMR, NDC, and rectification characteristics in ZB(Mn)NNRs offer an optimistic outlook for the conceptualization of future nanodevices.

Exploring the Odd–Even Effect, Current Stabilization, and Negative Differential Resistance in Carbon-Chain-Based Molecular Devices

The transport properties of molecular devices based on carbon chains are systematically investigated using a combination of non-equilibrium Green’s function (NEGF) and density functional theory (DFT) first-principle methods. In single-carbon-chain molecular devices, a distinct even–odd behavior of the current emerges, primarily influenced by the density of states (DOS) within the chain channel. Additionally, linear, monotonic currents exhibit Ohmic contact characteristics. In ladder-shaped carbon-chain molecular devices, a notable current stabilization behavior is observed, suggesting their potential utility as current stabilizers within circuits. We provide a comprehensive analysis of the transport properties of molecular devices featuring ladder-shaped carbon chains connecting benzene-ring molecules. The occurrence of negative differential resistance (NDR) in the low-bias voltage region is noted, with the possibility of manipulation by adjusting the position of the benzene-ring molecule. These findings offer a novel perspective on the potential applications of atom chains.

Surface-dominated conductance scaling in Weyl semimetal NbAs

Protected surface states arising from non-trivial bandstructure topology in semimetals can potentially enable advanced device functionalities in compute, memory, interconnect, sensing, and communication. This necessitates a fundamental understanding of surface-state transport in nanoscale topological semimetals. Here, we investigate quantum transport in a prototypical topological semimetal NbAs to evaluate the potential of this class of materials for beyond-Cu interconnects in highly-scaled integrated circuits. Using density functional theory (DFT) coupled with non-equilibrium Green’s function (NEGF) calculations, we show that the resistance-area RA product in NbAs films decreases with decreasing thickness at the nanometer scale, in contrast to a nearly constant RA product in ideal Cu films. This anomalous scaling originates from the disproportionately large number of surface conduction states which dominate the ballistic conductance by up to 70% in NbAs thin films. We also show that this favorable RA scaling persists even in the presence of surface defects, in contrast to RA sharply increasing with reducing thickness for films of conventional metals, such as Cu, in the presence of surface defects. These results underscore the potential of topological semimetals as future back-end-of-line (BEOL) interconnect metals.

First-Principles Simulation and Materials Screening for Spin-Orbit Torque in 2D van der Waals Heterostructures.

Recent advancements in spin-orbit torque (SOT) technology in two-dimensional van der Waals (2D vdW) materials have not only pushed spintronic devices to their atomic limits but have also unveiled unconventional torques and novel spin-switching mechanisms. The vast diversity of SOT observed in numerous 2D vdW materials necessitates a screening strategy to identify optimal materials for torque device performance. However, such a strategy has yet to be established. To address this critical issue, a combination of density functional theory and non-equilibrium Green's function is employed to calculate the SOT in various 2D vdW bilayer heterostructures. This leads to the discovery of three high SOT systems: WTe2/CrSe2, MoTe2/VS2, and NbSe2/CrSe2. Furthermore, a figure of merit that allows for rapid and efficient estimation of SOT is proposed, enabling high-throughput screening of optimal materials and devices for SOT applications in the future.