Effect of substitutional defects on resonant tunneling diodes based on armchair graphene and boron nitride nanoribbons lateral heterojunctions

  1. 1,2 ORCID Logo
1Department of Electrical Engineering, Faculty of Engineering, Arak University, Arak, 3815688349, Iran
2Institute of Nanosciences and Nanotechnology, Arak University, Arak, Iran
  1. Corresponding author email
Associate Editor: E. Meyer
Beilstein J. Nanotechnol. 2020, 11, 688–694. https://doi.org/10.3762/bjnano.11.56
Received 11 Nov 2019, Accepted 14 Apr 2020, Published 24 Apr 2020
Full Research Paper
cc by logo
Album

Abstract

A nanometer-scaled resonant tunneling diode based on lateral heterojunctions of armchair graphene and boron nitride nanoribbons, exhibiting negative differential resistance is proposed. Low-bandgap armchair graphene nanoribbons and high-bandgap armchair boron nitride nanoribbons are used to design the well and the barrier region, respectively. The effect of all possible substitutional defects (including BC, NC, CB, and CN) at the interface of graphene and boron nitride nanoribbons on the negative differential resistance behavior of the proposed resonant tunneling diode is investigated. Transport simulations are carried out in the framework of tight-binding Hamiltonians and non-equilibrium Green’s functions. The results show that a single substitutional defect at the interface of armchair graphene and boron nitride nanoribbons can dramatically affect the negative differential resistance behavior depending on its type and location in the structure.

Introduction

2D materials have gained tremendous research interest due to the unique properties that result from their atomic-scale thickness [1-5]. These materials, which include graphene, hexagonal boron nitride, and the large family of transition metal dichalcogenides, have electronic structures exhibiting metallic, semiconducting, and insulating properties. Novel electronic devices have been realized by heterostructures based on vertical stacking or lateral stitching of 2D materials with different electronic properties [6]. Lateral graphene/hexagonal boron nitride (Gr/hBN) heterostructures, due to very low lattice mismatch between graphene and hBN, are most suitable as platforms for fully two-dimensional nanoelectronic devices [7-11].

Resonant tunneling diodes (RTDs) are among various electronic devices realized on the platform of 2D Gr/hBN heterostructures [12-16]. RTDs exhibit negative differential resistance (NDR) and have a wide range of applications including ultra-fast switching devices, oscillators, frequency multipliers, one-transistor static memories and multi-valued memory circuits [12,17-20].

In a RTD, a material with low bandgap energy is sandwiched between two materials with larger bandgaps, i.e., a quantum well between two potential barriers, forming a so-called double-barrier quantum well structure. In the well, the energy of the electrons is quantized due to the quantum confinement of their wave function. Incident electrons with energies equal to the quantized levels of the well pass through the barriers with rather high transmission probabilities. Electrons with other energies have an extremely small chance of passing through. This causes RTDs to exhibit NDR in their current–voltage characteristic.

Conventionally, RTDs are made by vertical stacking of bulk semiconductor materials with different bandgap energies to form a planar 3D structure in which the direction of carrier transport is perpendicular to the interface of stacked materials [21-23]. In recent years, a few RTD structures based on 2D materials have been proposed [24-26]. In such RTDs the bandgap difference needed for normal operation is created by juxtaposing graphene nanoribbons (GNRs) with different widths (utilizing the inverse relation between GNR width and bandgap energy) or by periodically arranging graphene (the well) and boron nitride regions (the barriers). While the performance of conventional RTDs based on bulk semiconductors is degraded by dislocations and lattice mismatch at the interface of different bandgap materials, the RTDs based on heterojunctions between armchair graphene nanoribbons (AGNRs) and armchair boron nitride nanoribbons (ABNNRs) have shown superior performance because of the very low lattice mismatch between graphene and hBN [3]. However, inevitable interfacial defects located at the interface of Gr/hBN heterojunctions, including point defects (single vacancies and substitutional defects) and topological defects can alter the electronic properties of Gr/hBN heterostructures and, consequently, the performance of RTDs based on Gr/hBN heterojunctions [26-33]. Formation energy calculations have revealed that point defects occur preferentially at the interfaces of graphene and hBN domains rather than in the middle of these domains, and that substitutional defects are dominant [30]. Therefore, a reliable and accurate investigation of the electronic behavior in devices based on Gr/hBN heterostructures must consider the effect of such defects.

In this work, a nanometer-scaled RTD based on lateral AGNR/ABNNR heterojunctions is proposed and the effect of all possible types of substitutional defects at the interface of AGNR/ABNNR heterojunctions on the electronic behavior of the proposed RTD is investigated.

The effect of substitutional defects on the electronic transport properties of zigzag graphene nanoribbons symmetrically decorated with BN is described in [34]. However, it considers the interfaces of the Gr/hBN regions are parallel to the transport direction so that there are no bandgap variations in the transport direction. Since the electron transport is in zigzag direction, spin-polarized transport calculations are utilized. Instead, in the study presented here, the electron transport is in armchair direction, which is not spin-polarized. Therefore, nonequilibrium Green’s functions with tight-binding Hamiltonians (without considering spin degree of freedom) are utilized for electronic transport calculations. Furthermore, in the proposed RTD there are four GNR/BNNR heterojunctions perpendicular to the transport direction (Figure 1), which are required to construct the double-barrier quantum well structure.

Calculation Method

Figure 1a schematically shows the structure of the proposed RTD. An AGNR of 19 carbon atoms in width (19-AGNR) and four hexagonal carbon rings in length (1.7 × 1.6 nm2) is sandwiched between two ABNNRs of the same width but of only one hexagonal hBN ring in length (1.7 × 0.42 nm2). Two spacer regions are considered at both sides of the structure to exclude outward interfacial defects from the contact regions. The whole structure is assumed to be connected to semi-infinite AGNR contacts at both sides. The energy band edge diagram of the proposed RTD along with the quantized energy level of the well are shown in Figure 1b. Due to small lattice mismatch between graphene and hBN the edge bond relaxation correction for carbon atoms at the interface of C and hBN domains could be ignored [35]. Therefore, the electronic structure of the proposed RTD could be modeled via a tight-binding Hamiltonian, HD, including first nearest-neighbor interactions [36]. The hopping energies between C, B and N atoms and the on-site energies of B and N atoms are listed in Table 1 [37].

[2190-4286-11-56-1]

Figure 1: (a) Schematic representation of the proposed RTD. Carbon, boron and nitrogen atoms, are depicted in grey, red and yellow, respectively. The contacts are assumed to be heavily doped semi-infinite AGNRs. (b) Energy diagram of the proposed RTD (solid lines) and resonant energy of the well (dashed line).

Table 1: Tight-binding parameters for the CBN composite structurea.

εB εC εN tCC tCB tCN tBN
2.76 0.00 −1.64 2.65 2.25 1.70 2.40

aε and t denote the on-site and hopping energies, respectively (both in eV).

In nonequilibrium Green’s functions the retarded Green’s function of the device (scattering region) is written as [38]:

[2190-4286-11-56-i1]
(1)

where η is an infinitesimal positive number and the non-Hermitian self-energy matrix, Σ(E) = Σ1(E) + Σ2(E), represents the escape rate of electrons from the device into the semi-infinite contacts. The self-energy matrices are calculated through a highly convergent recursive method [39]. Then the transmission as a function of the energy is obtained via [39]:

[2190-4286-11-56-i2]
(2)

in which Ga(E) = (Gr(E)) is the advanced Green's function and Γj (j = 1, 2), represent the level broadening due to the coupling between device and contacts:

[2190-4286-11-56-i3]
(3)

Finally, the current through the device is calculated via [38]:

[2190-4286-11-56-i4]
(4)

in which q is the electron charge, f(E) is the Fermi–Dirac distribution function at the contacts and V is the bias voltage. The Fermi level at both contacts is assumed to be 0.03 eV above the conduction band edge [40]. Room temperature (300 K) is considered in all simulations.

Results and Discussion

Substitutional carbon atoms in boron or nitrogen sublattices (CB and CN) as well as boron or nitrogen atoms on carbon sites (BC and NC) are considered (Figure 2). Each defect could occur inside the well or in the contact regions. Since the typical defect concentration for real samples is one defect per ca. 10 Å, one defect of each type is considered at the Gr/hBN heterojunctions [41].

[2190-4286-11-56-2]

Figure 2: All possible substitutional defects at the interfaces of RTD based on AGNR/ABNNR heterojunctions. Carbon, boron and nitrogen atoms are depicted in grey, red and yellow, respectively.

Figure 3 shows the effect of substitutional defects on the energy bandstructure of 19-AGNR and 19-ABNNR. Both BC and NC defects increase the bandgap of AGNR (Figure 3b,c) due to breaking the symmetry between two graphene sublattices. A BC defect, due to the p-type character, shifts the conduction band edge toward higher energies, while an NC defect, due to the n-type character, shifts it toward lower energies. In contrast, both CB and CN defects decrease the bandgap of ABNNR. However, a CB defect in ABNNR, due to the n-type character, shifts the conduction band edge toward lower energies while the CN defect, due to the p-type character, shifts it toward higher energies (Figure 3e,f).

[2190-4286-11-56-3]

Figure 3: Effect of substitutional defects on the energy bandstructure of 19-AGNR and 19-ABNNR. Carbon, boron and nitrogen atoms are depicted in grey, red and yellow, respectively.

Figure 4a depicts the effect of CB and NC substitutional defects at the interface between the left ABNNR barrier and the well on the current–voltage characteristic of the proposed RTD. Both defects lower the on-site energy at the place of the substituted atom due to the n-type character, which shifts all energy levels including the conduction band edge toward lower energies. The NC defect inside the well also lowers the resonant energies of the well. Therefore, the transmission peaks move downward (Figure 4b), which causes the peak current (Ip) to increase and the peak voltage (Vp) to decrease (Figure 4a). A CB defect in the left barrier lowers the conduction band edge (Figure 3e). This also, by reducing the height of the left energy barrier, lowers the resonant energies of the well, which in turn moves the transmission peaks downward (Figure 4b). Lower transmission peaks translate to decreased Vp and increased Ip (Figure 4a).

[2190-4286-11-56-4]

Figure 4: Effect of CB and NC defects at the interface between the left ABNNR barrier and the well on (a) the current–voltage characteristic and (b) the transmission trough the proposed RTD as a function of the electron energy.

Both BC and CN substitutional defects have p-type character (adding an extra hole to the system) because the added atom has one electron less than the removed atom. Therefore, a BC defect in the well shifts all allowed energy levels (including discrete resonant energies) upward (Figure 3b), which causes Vp to increase and Ip to decrease (Figure 5a). A CN defect at the interface between the right barrier and the well (inside the barrier) also moves the conduction band edge upward, which in turn shifts the resonant energies of the well upward causing Vp to increase and Ip to decrease (Figure 5a).

[2190-4286-11-56-5]

Figure 5: Effect of BC and CN defects at the interface between the right barrier and the well on (a) current–voltage characteristic and (b) the transmission trough the proposed RTD as a function of the electron energy.

CB and NC defects may occur at the interface between the right barrier and the right spacer (see Figure 2). Because these substitutional defects are outside the well region they have no effect on the resonant energies of the well (Figure 6b). Therefore, in both cases Vp remains unchanged as illustrated in Figure 6a. Nevertheless, a CB defect in the right barrier lowers the conduction band edge and consequently the barrier height. Therefore, the transmission probabilities over the right barrier at all energies are increased compared to a defect-free RTD (red curve in Figure 6b). The NC defect in the right spacer also lowers the conduction band edge. This creates a local quantum well in the right spacer region, which, by localizing electron wave functions, reduces the transmission coefficients (yellow curve in Figure 6b). This causes Ip to decrease compared to a defect-free structure.

[2190-4286-11-56-6]

Figure 6: Effect of CB and NC defects at the interface between the right barrier and the right spacer on (a) current-voltage characteristic and (b) the transmission trough the proposed RTD as a function of the electron energy.

Figure 7 compares the current–voltage characteristic of RTDs including BC and CN defects at the interface between the left barrier and the left spacer with that of a defect-free RTD. If there is no defect in the left spacer, by increasing the bias voltage a quantum well is formed in this region, which localizes incident electron waves. However, a BC defect in the left spacer region shifts the conduction band edge upward. Thus, a potential barrier with a very small height (with respect to the ABNNR potential barrier) is formed. Although by increasing the bias voltage the height of this potential barrier is lowered, a quantum well is never formed in the spacer region in this case (at least up to a bias voltage of 250 mV, which is considered in this work). The absence of a local quantum well in the left spacer region (for bias voltages ranging from 0 to 250 mV) causes the transmission peak to vanish at larger bias voltages with respect to the defect-free structure. Therefore, the peak values of current and voltage are increased with respect to defect-free structure (Figure 7a). A CN defect in the left barrier region, besides reducing the bandgap, slightly lowers the conduction band edge (Figure 3f). This has a twofold effect: a small downward shift of the resonant energies of the well and a slight increase in transmission probabilities over the left barrier (Figure 7b), which in turn cause Ip to slightly increase and Vp to slightly decrease, as illustrated by the yellow curve in Figure 7a. Thus, a single substitutional defect can severely alter the NDR behavior of RTDs based on AGNR/ABNNR heterojunctions. This implies that by intentionally introducing such defects in devices based on AGNR/ABNNR heterojunctions, nanoelectronic devices with desired performance characteristics can be designed.

[2190-4286-11-56-7]

Figure 7: Effect of BC and CN defects at the interface between the left barrier and the left spacer on (a) current–voltage characteristic and (b) the transmission trough the proposed RTD as a function of the electron energy.

Conclusion

The effect of substitutional defects on the NDR behavior of a nanometer-scaled RTD based on 2D heterojunctions of AGNR/ABNNR was investigated. It was shown that a single substitutional defect, depending on its type and position, could severely alter the NDR behavior of the proposed RTD. While a BC defect inside the well region decreases Ip and increases Vp, it increases both Ip and Vp if located at the interface between left barrier and spacer. An NC defect inside the well increases Ip and decreases Vp, while it decreases Ip and does not affect Vp if located at the interface between right barrier and spacer. Although a CN defect inside the well has the same effect as a BC defect, it increases Ip and decreases Vp if located at the interface between left barrier and spacer. A CB defect inside the well has the same effect as an NC defect, while it decreases Ip and does not alter Vp if located at the interface between right barrier and spacer. Substitutional defects can be intentionally incorporated in 2D heterojunctions in a controllable manner. Therefore, since peak current and peak voltage of the proposed RTD depend on the position of the resonant energy of the well, and the resonant energy itself is highly sensitive to the type and position of substitutional defects, the intentional introduction of such defects can be utilized to design nanoscale RTDs with desired NDR characteristic and RTD-based strain or pressure sensors with improved sensitivity.

References

  1. Liu, B.; Abbas, A.; Zhou, C. Adv. Electron. Mater. 2017, 3, 1700045. doi:10.1002/aelm.201700045
    Return to citation in text: [1]
  2. Shim, J.; Park, H.-Y.; Kang, D.-H.; Kim, J.-O.; Jo, S.-H.; Park, Y.; Park, J.-H. Adv. Electron. Mater. 2017, 3, 1600364. doi:10.1002/aelm.201600364
    Return to citation in text: [1]
  3. Sanaeepur, M.; Jafari, M. R.; Esmaeili, M. IEEE Trans. Electron Devices 2020, 67, 725–729. doi:10.1109/ted.2019.2958421
    Return to citation in text: [1] [2]
  4. Sanaeepur, M. IEEE Trans. Nanotechnol. 2019, 18, 1181–1187. doi:10.1109/tnano.2019.2945354
    Return to citation in text: [1]
  5. Sanaeepur, M. IEEE Trans. Electromagn. Compat. 2019, 61, 532–537. doi:10.1109/temc.2018.2830182
    Return to citation in text: [1]
  6. Li, M.-Y.; Chen, C.-H.; Shi, Y.; Li, L.-J. Mater. Today 2016, 19, 322–335. doi:10.1016/j.mattod.2015.11.003
    Return to citation in text: [1]
  7. Fiori, G.; Betti, A.; Bruzzone, S.; Iannaccone, G. ACS Nano 2012, 6, 2642–2648. doi:10.1021/nn300019b
    Return to citation in text: [1]
  8. Sanaeepur, M.; Yazdanpanah Goharrizi, A.; Sharifi, M. J. Beilstein J. Nanotechnol. 2014, 5, 1569–1574. doi:10.3762/bjnano.5.168
    Return to citation in text: [1]
  9. Goharrizi, A. Y.; Sanaeepur, M.; Sharifi, M. J. Superlattices Microstruct. 2015, 85, 522–529. doi:10.1016/j.spmi.2015.04.035
    Return to citation in text: [1]
  10. Ebrahimi, M.; Horri, A.; Sanaeepur, M.; Tavakoli, M. B. J. Appl. Phys. 2020, 127, 084504. doi:10.1063/1.5130777
    Return to citation in text: [1]
  11. Ebrahimi, M.; Horri, A.; Sanaeepur, M.; Tavakoli, M. B. J. Comput. Electron. 2020, 19, 62–69. doi:10.1007/s10825-019-01442-z
    Return to citation in text: [1]
  12. Pereira, J. M., Jr.; Vasilopoulos, P.; Peeters, F. M. Appl. Phys. Lett. 2007, 90, 132122. doi:10.1063/1.2717092
    Return to citation in text: [1] [2]
  13. Nguyen, V. H.; Mazzamuto, F.; Bournel, A.; Dollfus, P. J. Phys. D: Appl. Phys. 2012, 45, 325104. doi:10.1088/0022-3727/45/32/325104
    Return to citation in text: [1]
  14. Song, Y.; Wu, H.-C.; Guo, Y. Appl. Phys. Lett. 2013, 102, 093118. doi:10.1063/1.4794952
    Return to citation in text: [1]
  15. Fallahazad, B.; Lee, K.; Kang, S.; Xue, J.; Larentis, S.; Corbet, C.; Kim, K.; Movva, H. C. P.; Taniguchi, T.; Watanabe, K.; Register, L. F.; Banerjee, S. K.; Tutuc, E. Nano Lett. 2015, 15, 428–433. doi:10.1021/nl503756y
    Return to citation in text: [1]
  16. Xiao, Z.; Ma, C.; Huang, J.; Liang, L.; Lu, W.; Hong, K.; Sumpter, B. G.; Li, A.-P.; Bernholc, J. Adv. Theory Simul. 2019, 2, 1800172. doi:10.1002/adts.201800172
    Return to citation in text: [1]
  17. Mazumder, P.; Kulkarni, S.; Bhattacharya, M.; Sun, J. P.; Haddad, G. I. Proc. IEEE 1998, 86, 664–686. doi:10.1109/5.663544
    Return to citation in text: [1]
  18. Pan, W.; Lyo, S. K.; Reno, J. L.; Simmons, J. A.; Li, D.; Brueck, S. R. J. Appl. Phys. Lett. 2008, 92, 052104. doi:10.1063/1.2840996
    Return to citation in text: [1]
  19. Hikita, Y.; Kourkoutis, L. F.; Susaki, T.; Muller, D. A.; Takagi, H.; Hwang, H. Y. Phys. Rev. B 2008, 77, 205330. doi:10.1103/physrevb.77.205330
    Return to citation in text: [1]
  20. Dalgleish, H.; Kirczenow, G. Nano Lett. 2006, 6, 1274–1278. doi:10.1021/nl060040m
    Return to citation in text: [1]
  21. Kovacic, S. J.; Ojha, J. J.; Simmons, J. G.; Jessop, P. E.; Mand, R. S.; SpringThorpe, A. J. IEEE Trans. Electron Devices 1993, 40, 1154–1160. doi:10.1109/16.214743
    Return to citation in text: [1]
  22. Seabaugh, A. C.; Luscombe, J. H.; Randall, J. N.; Colter, P. C.; Dip, A.; Eldallal, G. M.; Bedair, S. M. Thin Solid Films 1993, 225, 99–104. doi:10.1016/0040-6090(93)90135-c
    Return to citation in text: [1]
  23. Reed, M. A. Superlattices Microstruct. 1988, 4, 741–747. doi:10.1016/0749-6036(88)90206-6
    Return to citation in text: [1]
  24. Dubois, S. M.-M.; Declerck, X.; Charlier, J.-C.; Payne, M. C. ACS Nano 2013, 7, 4578–4585. doi:10.1021/nn401322t
    Return to citation in text: [1]
  25. Xu, Y.; Guo, Z.; Chen, H.; Yuan, Y.; Lou, J.; Lin, X.; Gao, H.; Chen, H.; Yu, B. Appl. Phys. Lett. 2011, 99, 133109. doi:10.1063/1.3643899
    Return to citation in text: [1]
  26. Li, M.; Zheng, B.; Duan, K.; Zhang, Y.; Huang, Z.; Zhou, H. J. Phys. Chem. C 2018, 122, 14945–14953. doi:10.1021/acs.jpcc.8b02750
    Return to citation in text: [1] [2]
  27. Liu, X.; Zhang, G.; Zhang, Y.-W. Nano Lett. 2016, 16, 4954–4959. doi:10.1021/acs.nanolett.6b01565
    Return to citation in text: [1]
  28. Ding, N.; Lei, Y.; Chen, X.; Deng, Z.; Ng, S.-P.; Wu, C.-M. L. Comput. Mater. Sci. 2016, 117, 172–179. doi:10.1016/j.commatsci.2015.12.052
    Return to citation in text: [1]
  29. Zhang, J.; Xie, W.; Xu, X.; Zhang, S.; Zhao, J. Chem. Mater. 2016, 28, 5022–5028. doi:10.1021/acs.chemmater.6b01764
    Return to citation in text: [1]
  30. Pruneda, J. M. Phys. Rev. B 2012, 85, 045422. doi:10.1103/physrevb.85.045422
    Return to citation in text: [1] [2]
  31. Li, Q.; Liu, M.; Zhang, Y.; Liu, Z. Small 2016, 12, 32–50. doi:10.1002/smll.201501766
    Return to citation in text: [1]
  32. Ding, N.; Chen, X.; Wu, C.-M. L. Sci. Rep. 2016, 6, 31499. doi:10.1038/srep31499
    Return to citation in text: [1]
  33. Wang, J.; Ma, F.; Liang, W.; Sun, M. Mater. Today Phys. 2017, 2, 6–34. doi:10.1016/j.mtphys.2017.07.001
    Return to citation in text: [1]
  34. Leon, C.; Costa, M.; Chico, L.; Latgé, A. Sci. Rep. 2019, 9, 3508. doi:10.1038/s41598-019-39763-5
    Return to citation in text: [1]
  35. Mukherjee, S.; Kaloni, T. P. J. Nanopart. Res. 2012, 14, 1059. doi:10.1007/s11051-012-1059-2
    Return to citation in text: [1]
  36. Sanaeepour, M.; Abedi, A.; Sharifi, M. J. IEEE Trans. Electron Devices 2017, 64, 1300–1304. doi:10.1109/ted.2017.2656944
    Return to citation in text: [1]
  37. Zhao, R.; Wang, J.; Yang, M.; Liu, Z.; Liu, Z. J. Phys. Chem. C 2012, 116, 21098–21103. doi:10.1021/jp306660x
    Return to citation in text: [1]
  38. Sanaeepur, M.; Goharrizi, A. Y.; Sharifi, M. J. IEEE Trans. Electron Devices 2014, 61, 1193–1198. doi:10.1109/ted.2013.2290049
    Return to citation in text: [1] [2]
  39. Sancho, M. P. L.; Sancho, J. M. L.; Sancho, J. M. L.; Rubio, J. J. Phys. F: Met. Phys. 1985, 15, 851–858. doi:10.1088/0305-4608/15/4/009
    Return to citation in text: [1] [2]
  40. Teong, H.; Lam, K.-T.; Liang, G. Jpn. J. Appl. Phys. 2009, 48, 04C156. doi:10.1143/jjap.48.04c156
    Return to citation in text: [1]
  41. Huang, B.; Liu, F.; Wu, J.; Gu, B.-L.; Duan, W. Phys. Rev. B 2008, 77, 153411. doi:10.1103/physrevb.77.153411
    Return to citation in text: [1]
Other Beilstein-Institut Open Science Activities