Next Article in Journal
Effects of Different Air Particle Abrasion Protocols on the Biaxial Flexural Strength and Fractography of High/Ultra-Translucent Zirconia
Next Article in Special Issue
Influence of the hBN Dielectric Layers on the Quantum Transport Properties of MoS2 Transistors
Previous Article in Journal
Effects of Unidirection/Bidirection Torsional Thermomechanical Processes on Grain Boundary Characteristics and Plasticity of Pure Nickel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bandstructure and Size-Scaling Effects in the Performance of Monolayer Black Phosphorus Nanodevices

Computational Nanoelectronics Group, Faculty of Electrical Engineering and Computing, University of Zagreb, HR 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Materials 2022, 15(1), 243; https://doi.org/10.3390/ma15010243
Submission received: 17 November 2021 / Revised: 22 December 2021 / Accepted: 24 December 2021 / Published: 29 December 2021
(This article belongs to the Special Issue Modeling of Advanced Metal-Oxide-Semiconductor Devices)

Abstract

:
Nanodevices based on monolayer black phosphorus or phosphorene are promising for future electron devices in high density integrated circuits. We investigate bandstructure and size-scaling effects in the electronic and transport properties of phosphorene nanoribbons (PNRs) and the performance of ultra-scaled PNR field-effect transistors (FETs) using advanced theoretical and computational approaches. Material and device properties are obtained by non-equilibrium Green’s function (NEGF) formalism combined with a novel tight-binding (TB) model fitted on ab initio density-functional theory (DFT) calculations. We report significant changes in the dispersion, number, and configuration of electronic subbands, density of states, and transmission of PNRs with nanoribbon width (W) downscaling. In addition, the performance of PNR FETs with 15 nm-long channels are self-consistently assessed by exploring the behavior of charge density, quantum capacitance, and average charge velocity in the channel. The dominant consequence of W downscaling is the decrease of charge velocity, which in turn deteriorates the ON-state current in PNR FETs with narrower nanoribbon channels. Nevertheless, we find optimum nanodevices with W > 1.4 nm that meet the requirements set by the semiconductor industry for the “3 nm” technology generation, which illustrates the importance of properly accounting bandstructure effects that occur in sub-5 nm-wide PNRs.

1. Introduction

Monolayer black phosphorus (BP) or phosphorene, illustrated in Figure 1a,b, is a promising two-dimensional (2D) material for the realization of future electron devices in integrated circuits due to its favorable electronic and transport properties, which include acceptable bandgap and carrier mobility [1,2]. Large-area BP field-effect transistors (FETs) have been experimentally demonstrated [2,3,4], and short-channel and wide-gate phosphorene FETs have been theoretically studied [5,6]. On the other hand, patterning monolayer BP into phosphorene nanoribbons (PNRs) provides a technologically relevant way of adjusting the electronic and transport properties by quantum confinement effects [7,8,9,10]. Ultra-narrow PNRs have been fabricated and characterized recently with widths down to ~0.5 nm [11,12], which further kindles interest in theoretical research of PNR-based nanoelectronic devices.
For a proper assessment of these nanostructured phosphorene FETs, methodology must be based on advanced theoretical formalisms such as quantum transport, e.g., non-equilibrium Green’s function (NEGF) formalism, because the transport physics in nanodevices is governed by quantum effects [13,14,15]. Moreover, these nanostructures must be described by proper atomically-resolved Hamiltonians that consider the complex bandstructure of such materials at the nanoscale [16,17]. The best bandstructure description is naturally provided by ab initio methods such as density-functional theory (DFT) calculations. However, the Hamiltonians obtained by DFT are very large and dense matrices which makes DFT Hamiltonians difficult to implement for quantum transport simulations for realistically sized FETs [5,17,18]. Tight-binding (TB) models are more readily applicable and result in more computationally efficient NEGF simulations. A TB model from the relevant literature (TBL model) for phosphorene, as described in [19], is widely used to explore the properties of BP and PNRs [9,20,21,22]. Nevertheless, in comparison to DFT results [23], the TBL model does not reproduce the intricate multi-valley bandstructure of PNRs with the widths under ~5 nm.
In [24,25], we introduced a new DFT- based TB Hamiltonian model (DFT-TB model) that describes the bandstructure of ultra-narrow PNRs more accurately. The DFT-TB model is more complicated than the TBL model, with larger and denser unit-cell and coupling matrices but is still numerically efficient in comparison to coupled DFT-NEGF simulations. In this paper, we employ NEGF and the DFT-TB model to study size-scaling and bandstructure effects in the electronic and transport properties of PNRs, and the performance of ballistic PNR FETs with nanoribbon widths under 5 nm. We show that PNR width downscaling significantly modifies the dispersion, including the bandgap and effective mass, together with a considerable impact on density of states and transmission through the nanoribbon. For PNR FETs, we demonstrate the deterioration of the ON-state current with decreasing PNR width; however, a surprising maximum of the width-normalized current is reported for ~2.5 nm-wide PNR FET. These findings are further investigated by examining the width- and bias-dependence of channel charge density and average charge velocity. We find that width-scaling-induced bandstructure effects in carrier velocity behavior are the dominant factor that determines the properties of PNR FET current driving capabilities in the ON-state.

2. Methods

2.1. DFT-Based Tight-Binding Hamiltonian

Due to the limitations of the TBL Hamiltonian model, especially in describing ultra-narrow PNRs, this work uses the recently developed DFT-TB model, which reproduces all important size-scaling effects on the bandstructure of PNRs [24,25]. In the development of the new DFT-TB model, the DFT simulations were performed using the OpenMX package [26,27], employing generalized gradient approximation (GGA) with Perdew–Burke–Ernzerhof (PBE) exchange-correlation (XC) functional. The DFT results were used as inputs into TBStudio [28], a new software package that implements the Slater–Koster (SK) method, and allows the choice of orbital number and type for fitting the DFT data [29]. In the new DFT-TB model, 4 orbitals are included for each phosphorus atom (s, px, py and pz), and all relevant SK overlap integrals are accounted for (ssσ, spσ, ppσ, ppπ). We showed in [24,25] that our new DFT-TB model achieves excellent agreement with DFT bandstructure in the energy range of interest, i.e., within ~1 eV from the conduction band minimum (CBM) and valence band maximum (VBM). The model data needed for the construction of PNR Hamiltonians can be found in Supplementary Materials. Electronic properties, transport properties, and PNR FET performance obtained using the new DFT-TB model are compared to those obtained with a simpler widely-used TBL model [19], in order to demonstrate the strong impact of bandstructure effects. These two TB Hamiltonians are used to study width-dependent dispersion of ultra-scaled PNRs and as inputs in the NEGF equations, which enables the investigation of the transport properties of PNRs and the performance of PNR FETs.

2.2. Quantum Transport with NEGF

The NEGF formalism is utilized in this work to solve the Schrödinger’s equation with open boundary conditions (OBCs). As will be explained later in Section 2.3, we need the transmission function and density of states of the PNR, so that only the equilibrium part of our in-house NEGF code for 2D material nanostructures [30,31] is needed. The central quantity of the NEGF formalism is the retarded Green’s function of the device and is obtained by
G R ( E ) = E + i 0 + I H Σ S R ( E ) Σ D R ( E ) 1
where H is the device Hamiltonian constructed using either the new DFT-TB or the existing TBL model. The size and sparsity of the Hamiltonian matrix depends on the model used, as well as on nanoribbon width (W) and length (L). The ΣR matrices are the retarded contact self-energies that consider the OBCs to the two contacts (source, S, and drain, D). These self-energy matrices are found by the iterative and numerically efficient Sancho–Rubio method [32]. The NEGF calculations in this work assume ideal contacts, meaning that the S/D extensions or reservoirs are semi-infinite semiconducting PNRs. This choice is common in the literature as it eliminates noncoherent effects at the channel-contact interfaces, and introduces no additional contact resistance into the nanostructure [33]. The retarded (GR) and advanced Green’s function (GA) of the device, where GA = GR, are then used to find the transmission function, T(E), and density of states, DOS(E), according to expressions given in e.g., [9,34,35].
Figure 1. Illustration of a phosphorene nanoribbon with (a) side and (b) top view, as well as of (c) PNR FET cross-sectional view. ToB model description with (d) capacitive model, (e) bandstructure along the channel length, and (f) example of a PNR bandstructure calculated by the DFT-TB model.
Figure 1. Illustration of a phosphorene nanoribbon with (a) side and (b) top view, as well as of (c) PNR FET cross-sectional view. ToB model description with (d) capacitive model, (e) bandstructure along the channel length, and (f) example of a PNR bandstructure calculated by the DFT-TB model.
Materials 15 00243 g001

2.3. Top-of-the-Barrier Device Model

The characteristics of ultra-scaled PNR FETs, illustrated in Figure 1c, and relevant device performance metrics are obtained using the top-of-the-barrier (ToB) device model, which self-consistently solves the Poisson equation that exploits the NEGF results and provides ballistic device characteristics [22,36,37]. The ToB model does not include tunneling, so its relevance is limited to the above-threshold operation region and can be utilized to adequately assess the ON-state performance in devices with channel lengths larger than ~10 nm [36,37].
In the ToB model, the central parameter is the top-of-the-barrier potential (Usc), located at the maximum of the source-drain barrier. Assuming a grounded source, the capacitive model defined in Figure 1d results in the following potential:
U s c = U s c 0 q α G V G S + α D V D S + q C o x Q s Q s ( V G S = 0 )
where the α parameters describe the capacitive coupling between the electrodes and the channel, e.g., gate coupling is defined as αG = CG/(CS + CD + CG). We set αG = 1 and αD = 0, assuming that the gate electrode exhibits ideal control over the atomically-thin nanoribbon channel. In (2), Cox is the gate oxide capacitance, and Qs is the inversion charge density in the channel at top-of-the-barrier. The Qs depends on the positions of Usc and Fermi levels in source and drain regions (EF,S and EF,D, respectively), which are shown in Figure 1e, and is found as follows:
Q s = Q s , S + Q s , D
Q s , S ( D ) = 1 2 + D O S ( E U s c ) f ( E E F , S ( D ) ) d E
with EF,D = EF,SqVDS. After convergence is achieved for Usc and Qs, the drain current is calculated using the Landauer formula [34]. In all calculations, we assume that devices operate at room temperature, i.e., T = 300 K; studying temperature-related effects, including self-heating, is beyond the scope of this work.
The charge density defined in (4) and the Landauer’s current formula take as inputs the DOS(E) and T(E), respectively, which are found by the NEGF simulations. Therefore, despite its relative simplicity, the ToB device model inherently includes all size-scaling and bandstructure-related effects with atomistic and orbital resolution. A simplified single-band effective-mass bandstructure in the conduction band is shown in Figure 1e, whereas our results are based on TB Hamiltonians, with an example of PNR dispersion obtained with the DFT-TB model shown in Figure 1f. Therefore, the described approach will allow us to explore accurately the impact of bandstructure and size scaling on the performance of ballistic PNR FETs.

3. Results and Discussions

3.1. Electronic and Transport Properties of Ultra-Narrow PNRs

Figure 2 reports the dispersion curves for various PNRs, with W ranging from 0.49 nm to 4.41 nm, calculated using the TBL and DFT-TB Hamiltonians. We observe that the DFT-TB Hamiltonian results in multi-valley bandstructure in both the conduction and valence band, which agrees with DFT studies of ultra-narrow PNRs reported in [23,38]. In contrast, the TBL Hamiltonian produces a symmetric single valley for both electrons and holes in the conduction and valence bands, respectively, irrespective of PNR width. Therefore, we expect a more accurate analysis of electronic, transport and device properties of PNRs and PNR FETs by using the DFT-TB model.
With downscaling the PNR width from 4.41 nm to 0.49 nm, bandgap (EG) increases considerably for both TB models, i.e., from 1.57 eV to 2.61 eV (TBL) and from 0.71 eV to 1.61 eV (DFT-TB). The TBL model provides wider bandgaps due to different XC functionals used for the development of the two Hamiltonian models. Namely, PBE was used for our DFT-TB model [24,25], which results in a lower EG in comparison to the HSE functional used for the development of the TBL model [19]. The PBE functionals are known to underestimate the bandgap, so the realistic EG value is expected to be between those obtained by PBE and HSE DFT simulations. Nevertheless, our DFT-TB model provides more accurate dispersion properties for both electrons and holes. Since the ToB device model relies on the bandstructure properties not far away from CBM and VBM, we avoid artificial EG adjustment in the DFT-TB model.
Due to the larger number of orbitals considered in the DFT-TB model, the related bandstructure plots in Figure 2 contain a larger number of subbands and, hence, conducting modes than those obtained by the TBL model. Consequently, NEGF simulations using the DFT-TB model should result in higher transmission probabilities at the same energy away from the CBM or VBM in comparison to the TBL results. When PNR width is scaled down, the number of subbands decreases for both models. However, improved bandstructure description by the DFT-TB model shows that even the narrowest PNRs, with W of 0.49 nm (Figure 2a) and 1.47 nm (Figure 2b), exhibit a much richer dispersion so that higher transmission probabilities and more enhanced DOS are expected in these devices than predicted by the TBL model.
Regarding carrier effective masses (m*), in wider PNRs, both TB models provide similar m* values near the CBM and VBM, while the difference increases considerably when the width is scaled down. For W = 4.41 nm, m* in the first subband equals ~0.28m0 (TBL) and ~0.21m0 (DFT-TB). For W = 0.49 nm, the DFT-TB model results in much heavier electrons in comparison to the TBL model. Namely, the DFT-TB model gives m*~1.9m0 in the first subband and m*~0.5m0 in the second subband, while the TBL model results in m*~0.5m0 in the first subband. The lower dispersion curvature, i.e., heavier carriers, in either the conduction or valence band generally leads to increased DOS that benefits inversion charge density in the channel. At the same time, heavier carriers exhibit diminished carrier velocities, which negatively impacts the current drivability of PNR FETs. Due to interplay of different phenomena, from subband number to effective mass change, it is difficult to qualitatively estimate how the PNR width downscaling will impact device performance purely from analyzing the dispersions.
Size-scaling effects observed in dispersion characteristics in Figure 2 consequently have a strong impact on the DOS and transmission of ultra-narrow PNRs. Figure 3 shows the DOS obtained with DFT-TB and TBL models for PNR widths of 0.49 nm, 1.47 nm, 2.45 nm, and 4.41 nm. When W is scaled down, DFT-TB DOS generally decreases away from the CBM and VBM due to reduced number of subbands in narrower nanoribbons. The two TB models provide approximately equal DOS in the vicinity of the CBM and VBM, within ~0.1 eV to ~0.2 eV. Further away from this energy range, the DFT-TB model gives much higher DOS values, which means that PNRs described by a more detailed Hamiltonian model will provide a higher amount of inversion charge that can be induced by appropriate adjustment of the quasi-Fermi level in the channel by the gate electrode, for both n- and p-type devices. For the narrowest PNR, the two DOS characteristics diverge considerably even near the CBM and VBM, which is caused by the large difference in m* between the two models. Therefore, we can expect the largest discrepancy in PNR FET performance between the two TB models for 0.49-nm-wide PNR FETs.
Connection between the bandstructure and transmission through the nanoribbon is illustrated in Figure 4a–c for the 1.47 nm-wide PNR. Namely, transmission function counts the number of conducting modes with a positive velocity at a certain energy and, hence, the DFT-TB model gives a generally higher transmission due to richer bandstructure than the TBL model. For W = 1.47 nm in Figure 4c, and other PNR widths from 0.49 nm to 4.41 nm reported in Figure 4d–f, we again observe that the TBL model reproduces the transmission well only near the CBM and VBM, while at higher energies E > CBM + 0.2 eV, the DFT-TB model provides a more complex characteristic and higher transmission values. Due to the higher number of conducting modes, the improved bandstructure description by the DFT-TB model results in up to ~4× higher transmission in comparison to the TBL model.
As shown in Figure 4, downscaling of W leads to an increased transmission gap and a reduced transmission maximum. In the examined energy range, from −2 eV to 2 eV, the transmission maxima decrease from 27, over 16 and 10, to 4 when PNR width decreases from 4.41 nm, over 2.45 nm and 1.47 nm, down to 0.49 nm. At the same time, the transport gap extracted as the energy range where T(E) < 0.01 increases from 0.71 eV for W = 4.41 nm up to 1.61 eV for the 0.49 nm-wide PNR. This increase makes narrower PNRs more immune to tunneling effects, which improves the performance of narrow PNR FETs in the OFF-state. It is expected that the reduced transmission in narrower PNRs will lead to poorer ON-state performance when W decreases, but device simulations are needed for a proper assessment due to several competing factors that play a role in device operation. For example, in the case of W = 0.49 nm, electrons are considerably heavier than in wider PNRs (Figure 2a), which results in high DOS (Figure 3a), and a high inversion charge density can be expected in 0.49 nm-wide PNR FETs. In addition, the transmission is boosted near the CBM in comparison to the TBL results (Figure 4d), which suggest a high drain current. However, the high effective mass (m*~1.9m0) points toward a lower carrier velocity and overall poorer current-driving capabilities of 0.49 nm-wide PNR FETs.

3.2. Performance of Ultra-Scaled PNR FETs

In this work, we investigate the performance of 15 nm-long PNR FETs, for which we assume SiO2 as gate dielectric with a thickness of 1 nm and S/D doping of m = 0.001, where m is the molar fraction of the areal density of P atoms in PNRs, resulting in a doping density of ~4 × 1012 cm−2. A common threshold voltage (VTH) of 0.24 V, as projected in the International Roadmap for Devices and Systems (IRDS) at the “3 nm” CMOS node [39], is set for all devices by automatically adjusting the gate work function. Setting the same VTH, and consequently the same OFF-state current (IOFF), allows a meaningful and fair comparison of PNR FETs with different nanoribbon widths. The resulting IOFF is ~1 nA/µm, defined for VDS = 0.7 V and VGS = 0 V. The supply voltage is 0.7 V, so the ON-state current (ION) is extracted from I-V characteristics for VDS = VGS = 0.7 V. Average charge velocity (vavg) at ToB is calculated from the drain current and the obtained ToB charge density. The ON-state ToB charge density (Qs,ON) and ON-state velocity (vON) are also extracted at the same bias point as ION, i.e., with gate and drain biased at the supply voltage.
Figure 5 reports the dependence of ION on PNR width for both the absolute magnitude of the current (Figure 5a) and for the width-normalized current (Figure 5b). The absolute ION monotonically decreases with the downscaling of PNR width, from 8.7 µA in the 4.41 nm-wide PNR FET down to 0.25 µA for W = 0.49 nm. In comparison to the TBL model, DFT-TB Hamiltonians with a more accurate PNR bandstructure provide higher driving currents, except for the narrowest device. These ION values in single PNRs are too low to be practically relevant, so several PNRs must be connected in parallel to provide sufficiently high ION. The plausibility of utilizing PNRs in extremely scaled FETs is quantified by the width-normalized ION reported in Figure 5b, which allows an assessment of PNR FET performance against IRDS requirements. As shown previously [22], ION obtained by the TBL model exhibits a generally decreasing trend with a weak modulation by W downscaling. Moreover, none of the examined PNR FETs using the TBL model fulfills the IRDS requirement for ION at the “3 nm” node, i.e., ION > 1.9 mA/µm [39]. In contrast, device simulations using DFT-TB Hamiltonians provide significant qualitative and quantitative changes. First, the ION vs. W characteristic is non-monotonic and the width-normalized ION exhibits a local maximum of 2.17 mA/µm for W = 2.45 nm, which means that 2.45 nm-wide PNR FETs are the most area-efficient devices in terms of current drivability. Regarding IRDS requirements, DFT-TB model reveals that PNR FETs with W > 1.4 nm can surpass the ION target set by the IRDS. The DFT-TB model reveals the severity of bandstructure effects in the narrowest devices because ION for W = 0.49 nm is only 0.51 mA/µm, which is ~2.7× lower than obtained by the TBL model. Going towards the widest examined nanoribbons, both Hamiltonian models converge to the same ION, which is expected as both approaches describe large-area phosphorene equally well.
The size-scaling bandstructure effects on the device performance can be more easily understood by exploring the size- and bias-dependent properties of charge density (Figure 6) and average charge velocity (Figure 7) in the channel of PNR FETs. As shown in Figure 6a, Qs increases with increasing VGS up to ~7 × 1012 cm−2, and generally increases with the downscaling of PNR width, although Qs-VGS curves are closely spaced for W > 1.5 nm. This Qs behavior is attributed to the increasing DOS near the CBM when W decreases (see Figure 3). Figure 6b compares the bias-dependence of Qs for the two TB models for 0.49 nm and 2.45 nm-wide phosphorene nanodevices. Clearly, the TBL model overestimates the channel charge density for W = 2.45 nm and underestimates it for W = 0.49 nm. The impact of width downscaling on Qs in the ON-state (Qs,ON) is reported in Figure 6c, and we observe that both models result in monotonic increase of Qs,ON when W decreases. The DFT-TB model provides a somewhat stronger modification of charge density with W, and the Qs,ON from DFT-TB simulations surpasses the TBL model results only for W < 1.26 nm. This result for W = 0.49 nm is expected because, from Figure 2a, it is clear that DFT-TB bandstructure exhibits heavier electrons in the lowest subbands and a generally higher number of subbands near the CBM when compared to the simpler TBL model. On the other hand, wider PNRs defined by DFT-TB Hamiltonians have slightly lighter carriers, which results in Qs,ON being somewhat lower in the 1.5 nm to 3.5 nm width range. Since the Qs behavior is in disagreement with ION trends, the reasons for improvement introduced by the DFT-TB model reported in Figure 5 must come from the gate capacitance (CG) or the average charge velocity.
Regarding CG, it can be assessed as a series of Cox and quantum capacitance of the channel (Cq) defined as Cq = d(qQs)/dVGS. Since gate oxide (SiO2) thickness is 1 nm for all devices, Cox = 34.5 fF/µm2 and is independent of PNR width. On the other hand, Cq increases from ~102 fF/µm2 for W = 4.41 nm up to ~248 fF/µm2 (TBL) and ~298 fF/µm2 (DFT-TB) for the 0.49 nm-wide device. Therefore, the dependence of Cq is qualitatively the same as Qs in Figure 6c. In addition, Cq is much larger than Cox, so the total CG only slightly increases with the downscaling of PNR width, and the absolute values of CG are 26–31 fF/µm2. If the oxide thickness were to decrease, Cox would rise, which in turn would enhance the relative impact of Cq in the total CG of PNR FETs. Consequently, bandstructure effects reported using the DFT-TB model would be more pronounced, but a more detailed investigation is beyond the scope of the current paper.
Figure 7a plots VGS-dependence of vavg for VDS = 0.7 V obtained using the DFT-TB model for various PNR widths. Generally, vavg increases with increasing VGS, while for W = 0.49 nm the velocity is independent of gate bias. The minimum vavg of 0.44 × 107 cm/s is reported for the 0.49 nm-wide PNR FET, while the highest velocity is reached for W = 2.45 nm irrespective of VGS. For W = 2.45 nm, vavg equals 1.79 × 107 cm/s at threshold and grows to 2.43 × 107 cm/s in the ON-state. The changes introduced by the improved Hamiltonian model are illustrated in Figure 7b that compares vavg in 0.49 nm and 2.45 nm-wide nanoribbons calculated using the two TB models. While both models give identical qualitative vavg behavior with respect to VGS, the TBL model underestimates the velocity for W = 2.45 nm and overestimates it for the 0.49 nm-wide PNR FET. In turn, these characteristics lead to weak W-dependence of ION for the TBL model reported in Figure 5b. Finally, the influence of PNR width downscaling on vON is illustrated in Figure 7c, and the curve exhibits a monotonic vON decrease in case of the TBL model. In contrast, a non-monotonic vON behavior is observed in case of the DFT-TB model with a maximum vON of 2.43 × 107 cm/s recorded in the 2.45 nm-wide PNR FET. This local vON maximum for W = 2.45 nm is a consequence of two competing mechanisms. As shown in Figure 2 and discussed in the related text, the first-subband effective mass monotonically increases in narrower PNRs, so a monotonic velocity decrease should occur with W downscaling. On the other hand, the average charge velocity considers all populated subbands, including higher subbands that exhibit larger effective masses. Therefore, the characteristics of the first subband are important, but not the only one responsible for the overall vON behavior. When the PNR width decreases, the separation between subbands increases, which in turn decreases their population and their relative impact on vON. Hence, 2.45 nm presents an optimum width where these two mechanisms jointly result in maximum average charge velocity in the ON-state. Finally, by comparing the results in Figure 5b and Figure 7c, we conclude that the ION characteristics reported in Figure 5b are dominantly caused by the bias- and width-dependence of carrier velocity in the channel.

4. Conclusions

Our results are obtained for ballistic transport across the channel, but in reality, carriers would likely experience scattering by intrinsic or extrinsic scattering centers, such as acoustic and optical phonons, Coulomb centers, defects, etc. [40,41]. Moreover, we neglect contact resistance in this study although it is a severe performance limiter in all 2D material-based FETs [42,43], including monolayer BP [4,44,45] and nanoribbon-based FET structures [22,33,46,47]. Nevertheless, our work illustrated the importance of using a proper Hamiltonian in the simulation of PNR nanodevices, quantified the magnitude of size-scaling and bandstructure effects, and demonstrated a considerable improvement of PNR FET figures-of-merit in comparison to the simpler TBL model. We found that for W = 2.45 nm (best device), the PNR FET is able to tolerate 13% current loss due to scattering and still meet the IRDS ION target at the “3 nm” technology node. The ballisticity level of 87% seems attainable because in a large-area phosphorene FET with a 15 nm-long channel, the ballisticity of ~90% was reported after including phonon scattering [48]. We note that the impact of electron-phonon scattering on the transport properties and performance of phosphorene nanodevices is under investigation, and the ballisticity level reported in [48] might be lower [18].
Monolayer BP or phosphorene is an attractive candidate for future electronic devices, while nanostructured BP in the form of PNRs offers an additional avenue for adjusting the electronic, transport, and device properties by quantum confinement effects. The nature of nanoscale BP devices demands advanced theoretical approaches based on quantum transport and appropriate atomically-resolved device Hamiltonians. In this work, we explored ultra-narrow PNRs and PNR FETs with the widths under ~5 nm using NEGF simulations based on a recently developed DFT-TB Hamiltonian model that accurately reproduces multi-valley dispersion and valence-conduction band asymmetry observed in PNRs by ab initio calculations. We explored the dispersion, DOS, and transmission through PNRs of various widths and reported a strong impact of width-scaling on the number and shape of conducting subbands, bandgap, and carrier effective mass. Focusing on 15 nm-long ballistic PNR FETs, we found that the ON-state inversion charge density increases up to ~7 × 1012 cm−2, and that the average charge velocity decreases considerably (from ~2.2 × 107 cm/s to ~0.4 × 107 cm/s) with the downscaling of PNR width. Velocity decrease was found to be the dominant factor in current-driving properties, so the ON-state current in PNR FETs also declines in narrower nanoribbons, from ~2 mA/µm (W = 4.41 nm) to ~0.5 mA/µm (W = 0.49 nm). Nevertheless, using an improved bandstructure description with the DFT-TB model revealed that ballistic PNR FETs with W > 1.4 nm can meet the IRDS requirement for ION at the “3 nm” CMOS technology node. Moreover, an optimum PNR FET with W = 2.45 nm was found, which exhibits ION ~ 2.2 mA/µm and which can operate at 87% of the ballistic limit and still meet the IRDS target for the ON-state current.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ma15010243/s1, File S1: compressed folder containing the tight-binding model data (matrices) for the construction of nanoribbon Hamiltonians.

Author Contributions

Conceptualization, M.P.; methodology, M.P. and M.M.; data curation, M.M. and M.P.; writing—original draft preparation, M.P. and M.M.; writing—review and editing, M.P.; project administration, M.P.; funding acquisition, M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Croatian Science Foundation under the project CONAN2D (Grant No. UIP-2019-04-3493).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are contained within the article and are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, H.; Neal, A.T.; Zhu, Z.; Luo, Z.; Xu, X.; Tománek, D.; Ye, P.D. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033–4041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Das, S.; Demarteau, M.; Roelofs, A. Ambipolar Phosphorene Field Effect Transistor. ACS Nano 2014, 8, 11730–11738. [Google Scholar] [CrossRef] [PubMed]
  3. Li, L.; Yu, Y.; Ye, G.J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X.H.; Zhang, Y. Black Phosphorus Field-Effect Transistors. Nat. Nanotechnol. 2014, 9, 372–377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Haratipour, N.; Robbins, M.C.; Koester, S.J. Black Phosphorus P-MOSFETs with 7-Nm HfO2 Gate Dielectric and Low Contact Resistance. IEEE Electron Device Lett. 2015, 36, 411–413. [Google Scholar] [CrossRef]
  5. Klinkert, C.; Szabó, Á.; Stieger, C.; Campi, D.; Marzari, N.; Luisier, M. 2-D Materials for Ultrascaled Field-Effect Transistors: One Hundred Candidates under the Ab Initio Microscope. ACS Nano 2020, 14, 8605–8615. [Google Scholar] [CrossRef]
  6. Cao, X.; Guo, J. Simulation of Phosphorene Field-Effect Transistor at the Scaling Limit. IEEE Trans. Electron Devices 2015, 62, 659–665. [Google Scholar] [CrossRef]
  7. Guo, H.; Lu, N.; Dai, J.; Wu, X.; Zeng, X.C. Phosphorene Nanoribbons, Phosphorus Nanotubes, and van Der Waals Multilayers. J. Phys. Chem. C 2014, 118, 14051–14059. [Google Scholar] [CrossRef] [Green Version]
  8. Taghizadeh Sisakht, E.; Zare, M.H.; Fazileh, F. Scaling Laws of Band Gaps of Phosphorene Nanoribbons: A Tight-Binding Calculation. Phys. Rev. B 2015, 91, 085409. [Google Scholar] [CrossRef] [Green Version]
  9. Poljak, M.; Suligoj, T. Immunity of Electronic and Transport Properties of Phosphorene Nanoribbons to Edge Defects. Nano Res. 2016, 9, 1723–1734. [Google Scholar] [CrossRef]
  10. Poljak, M.; Suligoj, T. The Potential of Phosphorene Nanoribbons as Channel Material for Ultrascaled Transistors. IEEE Trans. Electron Devices 2018, 65, 290–294. [Google Scholar] [CrossRef]
  11. Watts, M.C.; Picco, L.; Russell-Pavier, F.S.; Cullen, P.L.; Miller, T.S.; Bartuś, S.P.; Payton, O.D.; Skipper, N.T.; Tileli, V.; Howard, C.A. Production of Phosphorene Nanoribbons. Nature 2019, 568, 216–220. [Google Scholar] [CrossRef]
  12. Zhang, W.; Enriquez, H.; Tong, Y.; Mayne, A.J.; Bendounan, A.; Smogunov, A.; Dappe, Y.J.; Kara, A.; Dujardin, G.; Oughaddou, H. Flat Epitaxial Quasi-1D Phosphorene Chains. Nat. Commun. 2021, 12, 5160. [Google Scholar] [CrossRef]
  13. Luisier, M.; Klimeck, G. Atomistic Full-Band Simulations of Silicon Nanowire Transistors: Effects of Electron-Phonon Scattering. Phys. Rev. B 2009, 80, 155430. [Google Scholar] [CrossRef] [Green Version]
  14. Pizzi, G.; Gibertini, M.; Dib, E.; Marzari, N.; Iannaccone, G.; Fiori, G. Performance of Arsenene and Antimonene Double-Gate MOSFETs from First Principles. Nat. Commun. 2016, 7, 12585. [Google Scholar] [CrossRef] [Green Version]
  15. Afzalian, A.; Pourtois, G. ATOMOS: An ATomistic MOdelling Solver for Dissipative DFT Transport in Ultra-Scaled HfS2 and Black Phosphorus MOSFETs. In Proceedings of the 2019 International Conference on Simulation of Semiconductor Processes and Devices (SISPAD), Udine, Italy, 4–6 September 2019; pp. 1–4. [Google Scholar]
  16. Marin, E.G.; Perucchini, M.; Marian, D.; Iannaccone, G.; Fiori, G. Modeling of Electron Devices Based on 2-D Materials. IEEE Trans. Electron Devices 2018, 65, 4167–4179. [Google Scholar] [CrossRef]
  17. Afzalian, A.; Akhoundi, E.; Gaddemane, G.; Duflou, R.; Houssa, M. Advanced DFT–NEGF Transport Techniques for Novel 2-D Material and Device Exploration Including HfS2/WSe2 van Der Waals Heterojunction TFET and WTe2/WS2 Metal/Semiconductor Contact. IEEE Trans. Electron Devices 2021, 68, 5372–5379. [Google Scholar] [CrossRef]
  18. Afzalian, A. Ab Initio Perspective of Ultra-Scaled CMOS from 2D-Material Fundamentals to Dynamically Doped Transistors. npj 2D Mater. Appl. 2021, 5, 5. [Google Scholar] [CrossRef]
  19. Rudenko, A.N.; Katsnelson, M.I. Quasiparticle Band Structure and Tight-Binding Model for Single- and Bilayer Black Phosphorus. Phys. Rev. B 2014, 89, 201408. [Google Scholar] [CrossRef] [Green Version]
  20. Poljak, M. Electron Mobility in Defective Nanoribbons of Monoelemental 2D Materials. IEEE Electron Dev. Lett. 2020, 41, 151–154. [Google Scholar] [CrossRef]
  21. Yin, D.; Yoon, Y. Design Strategy of Two-Dimensional Material Field-Effect Transistors: Engineering the Number of Layers in Phosphorene FETs. J. Appl. Phys. 2016, 119, 214312. [Google Scholar] [CrossRef]
  22. Poljak, M.; Matić, M.; Zeljko, A. Minimum Contact Resistance in Monoelemental 2D Material Nanodevices with Edge-Contacts. IEEE Electron Device Lett. 2021, 42, 1240–1243. [Google Scholar] [CrossRef]
  23. Liang, L.; Wang, J.; Lin, W.; Sumpter, B.G.; Meunier, V.; Pan, M. Electronic Bandgap and Edge Reconstruction in Phosphorene Materials. Nano Lett. 2014, 14, 6400–6406. [Google Scholar] [CrossRef]
  24. Poljak, M.; Matić, M. Bandstructure Effects in Phosphorene Nanoribbon MOSFETs from NEGF Simulations Using a New DFT-Based Tight-Binding Hamiltonian Model. In Proceedings of the 2021 International Conference on Simulation of Semiconductor Processes and Devices (SISPAD), Dallas, TX, USA, 27–29 September 2021; pp. 180–183. [Google Scholar]
  25. Poljak, M.; Matić, M. DFT-Based Tight-Binding Model for Atomistic Simulations of Phosphorene Nanoribbons. In Proceedings of the 2021 44th International Convention on Information, Communication and Electronic Technology (MIPRO), Opatija, Croatia, 27 September–1 October 2021; pp. 80–84. [Google Scholar]
  26. Ozaki, T. Variationally Optimized Atomic Orbitals for Large-Scale Electronic Structures. Phys. Rev. B 2003, 67, 155108. [Google Scholar] [CrossRef]
  27. OpenMX Website. Available online: http://www.openmx-square.org/ (accessed on 9 February 2021).
  28. Nakhaee, M.; Ketabi, S.A.; Peeters, F.M. Tight-Binding Studio: A Technical Software Package to Find the Parameters of Tight-Binding Hamiltonian. Comput. Phys. Commun. 2020, 254, 107379. [Google Scholar] [CrossRef]
  29. Slater, J.C.; Koster, G.F. Simplified LCAO Method for the Periodic Potential Problem. Phys. Rev. 1954, 94, 1498–1524. [Google Scholar] [CrossRef]
  30. Poljak, M.; Song, E.B.; Wang, M.; Suligoj, T.; Wang, K.L. Influence of Edge Defects, Vacancies, and Potential Fluctuations on Transport Properties of Extremely Scaled Graphene Nanoribbons. IEEE Trans. Electron Devices 2012, 59, 3231–3238. [Google Scholar] [CrossRef]
  31. Poljak, M.; Glavan, M.; Kuzmić, S. Accelerating Simulation of Nanodevices Based on 2D Materials by Hybrid CPU-GPU Parallel Computing. In Proceedings of the 2019 42nd International Convention on Information and Communication Technology, Electronics and Microelectronics (MIPRO), Opatija, Croatia, 20–24 May 2019; pp. 51–56. [Google Scholar]
  32. Sancho, M.P.L.; Sancho, J.M.L.; Sancho, J.M.L.; Rubio, J. Highly Convergent Schemes for the Calculation of Bulk and Surface Green Functions. J. Phys. F Met. Phys. 1985, 15, 851. [Google Scholar] [CrossRef]
  33. Poljak, M.; Matić, M. Metallization-Induced Quantum Limits of Contact Resistance in Graphene Nanoribbons with One-Dimensional Contacts. Materials 2021, 14, 3670. [Google Scholar] [CrossRef] [PubMed]
  34. Datta, S. Quantum Transport: Atom to Transistor, 2nd ed.; Cambridge University Press: New York, NY, USA, 2005; ISBN 0-521-63145-9. [Google Scholar]
  35. Poljak, M. Intrinsic Capabilities for Digital Switching of Silicene Nanoribbons with Edge Defects. IEEE Trans. Electron Devices 2020, 67, 354–359. [Google Scholar] [CrossRef]
  36. Rahman, A.; Guo, J.; Datta, S.; Lundstrom, M.S. Theory of Ballistic Nanotransistors. IEEE Trans. Electron Devices 2003, 50, 1853–1864. [Google Scholar] [CrossRef] [Green Version]
  37. Kaneko, S.; Tsuchiya, H.; Kamakura, Y.; Mori, N.; Ogawa, M. Theoretical Performance Estimation of Silicene, Germanene, and Graphene Nanoribbon Field-Effect Transistors under Ballistic Transport. Appl. Phys. Express 2014, 7, 035102. [Google Scholar] [CrossRef]
  38. Carvalho, A.; Rodin, A.S.; Neto, A.H.C. Phosphorene Nanoribbons. Europhys. Lett. 2014, 108, 47005. [Google Scholar] [CrossRef]
  39. IEEE International Roadmap for Devices and Systems (IRDS), 2021 Edition. Available online: https://irds.ieee.org/ (accessed on 1 October 2021).
  40. Gaddemane, G.; Vandenberghe, W.G.; Van de Put, M.L.; Chen, S.; Tiwari, S.; Chen, E.; Fischetti, M.V. Theoretical Studies of Electronic Transport in Monolayer and Bilayer Phosphorene: A Critical Overview. Phys. Rev. B 2018, 98, 115416. [Google Scholar] [CrossRef] [Green Version]
  41. Gaddemane, G.; Van de Put, M.L.; Vandenberghe, W.G.; Chen, E.; Fischetti, M.V. Monte Carlo Analysis of Phosphorene Nanotransistors. J. Comput. Electron. 2021, 20, 60–69. [Google Scholar] [CrossRef]
  42. Driussi, F.; Venica, S.; Gahoi, A.; Gambi, A.; Giannozzi, P.; Kataria, S.; Lemme, M.C.; Palestri, P.; Esseni, D. Improved Understanding of Metal–Graphene Contacts. Microelectron. Eng. 2019, 216, 111035. [Google Scholar] [CrossRef]
  43. Gahoi, A.; Kataria, S.; Driussi, F.; Venica, S.; Pandey, H.; Esseni, D.; Selmi, L.; Lemme, M.C. Dependable Contact Related Parameter Extraction in Graphene–Metal Junctions. Adv. Electron. Mater. 2020, 6, 2000386. [Google Scholar] [CrossRef]
  44. Du, Y.; Liu, H.; Deng, Y.; Ye, P.D. Device Perspective for Black Phosphorus Field-Effect Transistors: Contact Resistance, Ambipolar Behavior, and Scaling. ACS Nano 2014, 8, 10035–10042. [Google Scholar] [CrossRef] [Green Version]
  45. Telesio, F.; le Gal, G.; Serrano-Ruiz, M.; Prescimone, F.; Toffanin, S.; Peruzzini, M.; Heun, S. Ohmic Contact Engineering in Few–Layer Black Phosphorus: Approaching the Quantum Limit. Nanotechnology 2020, 31, 334002. [Google Scholar] [CrossRef]
  46. Poljak, M.; Matić, M. Quantum Transport Simulations of Phosphorene Nanoribbon MOSFETs: Effects of Metal Contacts, Ballisticity and Series Resistance. In Proceedings of the 2020 International Conference on Simulation of Semiconductor Processes and Devices (SISPAD), Kobe, Japan, 23 September–6 October 2020; pp. 371–374. [Google Scholar]
  47. Poljak, M. Impact of Width Scaling and Parasitic Series Resistance on the Performance of Silicene Nanoribbon MOSFETs. IEEE Trans. Electron Devices 2020, 67, 4705–4708. [Google Scholar] [CrossRef]
  48. Szabo, A.; Rhyner, R.; Carrillo-Nunez, H.; Luisier, M. Phonon-Limited Performance of Single-Layer, Single-Gate Black Phosphorus n- and p-Type Field-Effect Transistors. In Proceedings of the 2015 IEEE International Electron Devices Meeting (IEDM), Washington, DC, USA, 7–9 December 2015; pp. 297–300. [Google Scholar]
Figure 2. Comparison of dispersions obtained by the DFT-TB (panels on the right) and TBL (panels on the left) models for PNR widths of (a) 0.49 nm, (b) 1.47 nm, (c) 2.45 nm, and (d) 4.41 nm.
Figure 2. Comparison of dispersions obtained by the DFT-TB (panels on the right) and TBL (panels on the left) models for PNR widths of (a) 0.49 nm, (b) 1.47 nm, (c) 2.45 nm, and (d) 4.41 nm.
Materials 15 00243 g002
Figure 3. Density of states in PNRs obtained using the DFT-TB and TBL models for PNR widths of (a) 0.49 nm, (b) 1.47 nm, (c) 2.45 nm, and (d) 4.41 nm.
Figure 3. Density of states in PNRs obtained using the DFT-TB and TBL models for PNR widths of (a) 0.49 nm, (b) 1.47 nm, (c) 2.45 nm, and (d) 4.41 nm.
Materials 15 00243 g003
Figure 4. Dispersion obtained by (a) TBL and (b) DFT-TB model, and (c) transmission for the 1.47 nm-wide PNR. Comparison of transmission functions obtained with the two TB Hamiltonian models for W of (d) 0.49 nm, (e) 2.45 nm, and (f) 4.41 nm.
Figure 4. Dispersion obtained by (a) TBL and (b) DFT-TB model, and (c) transmission for the 1.47 nm-wide PNR. Comparison of transmission functions obtained with the two TB Hamiltonian models for W of (d) 0.49 nm, (e) 2.45 nm, and (f) 4.41 nm.
Materials 15 00243 g004
Figure 5. Impact of PNR width downscaling on (a) absolute ION and (b) width-normalized ION in PNR FETs. The plots compare the results obtained by DFT-TB and TBL models.
Figure 5. Impact of PNR width downscaling on (a) absolute ION and (b) width-normalized ION in PNR FETs. The plots compare the results obtained by DFT-TB and TBL models.
Materials 15 00243 g005
Figure 6. (a) Charge density vs. gate bias characteristics in PNR FETs. Comparison of (b) Qs-VGS curves for 0.49 nm and 2.45 nm-wide devices, and (c) Qs-W characteristics in the ON-state obtained by the two TB models.
Figure 6. (a) Charge density vs. gate bias characteristics in PNR FETs. Comparison of (b) Qs-VGS curves for 0.49 nm and 2.45 nm-wide devices, and (c) Qs-W characteristics in the ON-state obtained by the two TB models.
Materials 15 00243 g006
Figure 7. (a) Average charge velocity vs. VGS in PNR FETs for the DFT-TB model. Comparison of (b) vavg-VGS curves for 0.49 nm and 2.45 nm-wide PNR FETs, and (c) ON-state vavg-W characteristics obtained by the DFT-TB and TBL models.
Figure 7. (a) Average charge velocity vs. VGS in PNR FETs for the DFT-TB model. Comparison of (b) vavg-VGS curves for 0.49 nm and 2.45 nm-wide PNR FETs, and (c) ON-state vavg-W characteristics obtained by the DFT-TB and TBL models.
Materials 15 00243 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Poljak, M.; Matić, M. Bandstructure and Size-Scaling Effects in the Performance of Monolayer Black Phosphorus Nanodevices. Materials 2022, 15, 243. https://doi.org/10.3390/ma15010243

AMA Style

Poljak M, Matić M. Bandstructure and Size-Scaling Effects in the Performance of Monolayer Black Phosphorus Nanodevices. Materials. 2022; 15(1):243. https://doi.org/10.3390/ma15010243

Chicago/Turabian Style

Poljak, Mirko, and Mislav Matić. 2022. "Bandstructure and Size-Scaling Effects in the Performance of Monolayer Black Phosphorus Nanodevices" Materials 15, no. 1: 243. https://doi.org/10.3390/ma15010243

APA Style

Poljak, M., & Matić, M. (2022). Bandstructure and Size-Scaling Effects in the Performance of Monolayer Black Phosphorus Nanodevices. Materials, 15(1), 243. https://doi.org/10.3390/ma15010243

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop