Next Article in Journal
Recycled Glass as a Substitute for Quartz Sand in Silicate Products
Next Article in Special Issue
Comparative Research of GaN Growth Mechanisms on Patterned Sapphire Substrates with Sputtered AlON Nucleation Layers
Previous Article in Journal
Dental Implants with Different Neck Design: A Prospective Clinical Comparative Study with 2-Year Follow-Up
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Carbon-Based Band Gap Engineering in the h-BN Analytical Modeling

by
Mohammad Taghi Ahmadi
1,2,*,
Ahmad Razmdideh
3,
Seyed Saeid Rahimian Koloor
4 and
Michal Petrů
4,*
1
Division of Computational Physics, Institute for Computational Science, Ton Duc Thang University, Ho Chi Minh City 758307, Vietnam
2
Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City 758307, Vietnam
3
Nano-electronic Research Group, Physics Department, Faculty of Science, Urmia University, Urmia 57147, Iran
4
Institute for Nanomaterials, Advanced Technologies and Innovation, Technical University of Liberec, Studentska 2, 461 17 Liberec, Czech Republic
*
Authors to whom correspondence should be addressed.
Materials 2020, 13(5), 1026; https://doi.org/10.3390/ma13051026
Submission received: 31 December 2019 / Revised: 14 February 2020 / Accepted: 18 February 2020 / Published: 25 February 2020
(This article belongs to the Special Issue Wide-Bandgap Materials and Applications)

Abstract

:
The absence of a band gap in graphene is a hindrance to its application in electronic devices. Alternately, the complete replacement of carbon atoms with B and N atoms in graphene structures led to the formation of hexagonal boron nitride (h-BN) and caused the opening of its gap. Now, an exciting possibility is a partial substitution of C atoms with B and N atoms in the graphene structure, which caused the formation of a boron nitride composite with specified stoichiometry. BC2N nanotubes are more stable than other triple compounds due to the existence of a maximum number of B–N and C–C bonds. This paper focused on the nearest neighbor’s tight-binding method to explore the dispersion relation of BC2N, which has no chemical bond between its carbon atoms. More specifically, the band dispersion of this specific structure and the effects of energy hopping in boron–carbon and nitrogen–carbon atoms on the band gap are studied. Besides, the band structure is achieved from density functional theory (DFT) using the generalized gradient approximations (GGA) approximation method. This calculation shows that this specific structure is semimetal, and the band gap energy is 0.167 ev.

1. Introduction

Nowadays, frontier technologies, especially in the semiconductor and battery industries, are improving material scaling to enhance the efficiency of their products. In particular, foiling unwanted transport between anode and cathode in the batteries is engaged by material performance. One the other hand, from the semiconductor-based point of view, one of the available techniques to progress the artifact is hidden in the BN material performance, which can be controlled by supplementary carbon atoms. For example, in the form of BC2N, the band gap can be controlled effortlessly. On the other hand, the BC2N allotropes as negative electrodes and anchoring materials for Li-S has been reported [1,2]. Since an increase in the numbers and efficiency of transistors diminishes the cost and leads to computers with faster processing ability, the number of such transistors on a microchip is commonly on the rise. Given this, it seems very cost effective to reduce the size of transistors [3,4]. However, this minimizing process will be finally stopped due to the need for developments in electronic commerce; thus, alternative technologies are needed [3]. Hence, the search for a matter with good nanoelectronics and optoelectronics functions is an important issue [5]. Alternative technology should solve past problems and have economic justification. Hereupon, to deal with this problem, the nanotechnology was integrated with electronic industry, and therefore, nanoelectronics was born [6].
Over the past 50 years, the use of Moore’s law has been effective in relieving silicon-based electronic components for nanoelectronics [7] and optoelectronics [8]. The two-dimensional graphene, due to its extraordinary electronic properties including high carrier velocity and mechanical strengths, is a good candidate for electronic devices [9,10]. However, the lack of a gaping band prevents it from being used in field effect transistors (MOSFET) [11], and this has led to the use of two-dimensional minerals such as graphite or single-layer boron nitride. Moreover, boron nitride has a structural similarity to carbon (C) and appears in different crystalline forms [12]. Besides, unlike graphene, due to the specific accumulation of monolayers by boron nitride, it is rarely available. In addition, because of the polarity of the bonds between B–N atoms in forming a multilayer of boron nitride, the whole system will be stable [13]. Boron nitride is known as a semiconductor with a large gap of 5 to 6 electron volts due to its wideband semiconducting behavior [14,15]. Furthermore, an interesting possibility is the substitution of the C atoms with B and N atoms in the sites of the honeycomb structure of graphene [16], which leads to the formation of a triplet-stoichiometric compound.
In 2006, the structure and electronic properties of BXCYNZ with stoichiometry ratios and various formations of the Ab Initio calculations were studied [17]. The results show that the amount of stability and energy gap depends on the array of atoms and their stoichiometry properties and is not related to the dimensions of the cell [17]. The structural similarity between graphite and hexagonal boron nitride became a motive for synthesizing the alloys of this compound [18]. On the one hand, the electronic properties of BXCYNZ are similar to graphite in a semimetal manner [19], and the insulator attitude is similar to that of boron nitride based on the doping with boron and nitrogen atoms [20]; therefore, the location and the amount of impurities affect the structure’s stability. Meanwhile, the structure of BC2N is zigzag-bonded, and it makes separate units of graphite and island-like boron nitride [21]. Therefore, determining the preferable structure of the BC2N monolayer, configurational strain, and the atomic bond has a considerable role [22].
It was theoretically predicted that through the use of first-principles calculation and CVD (chemical vapor deposition) synthesis in the presence of BC13, CC14, N2, and H2 based on XRD analysis, the existence of hexagonal layers by a 2.44 angstrom lattice constant with a 3.40 angstrom interlayer space has been confirmed [23]. Although the X-ray data fail to provide the atomic arrangement in planar form, they do designate that the boron, carbon, and nitrogen atom arrangement in the BC2N is in the form of sp2-type bonding [24]. As a result, it has been anticipated that the atoms will be arranged in planes of hexagonal rings similar to graphite, with weak interplanar interactions [25]. At this point, it remains vague whether BC2N functions as a semiconductor or a semimetal [26]. Furthermore, presumably the struggle between strain and bonding energy leads to the single-layered BC2N structure. Hence, whether BC2N is a semiconductor or a semimetal is an issue of controversy [26]. Therefore, band-gap engineering is a traditional and an influential technique in adjusting a new energy gap to a desirable value and designing new semiconductor materials [27]. To this end, there are diverse methods, including atom arrangement change, doping, and applying electric fields, among others.
The present study is concerned with the analytical investigation of energy dispersion through the use of the Nearest Neighbor Tight Binding (NNTB) approximation method [28] for BC2N material in which boron and nitrogen are both bonded with three carbons in a way that none of the carbons are bonded with each other, as shown in Figure 1. Another feature of this structure is that it has no symmetry in the space group, but the angle between the straight lines of combined boron and nitrogen atoms is 60 degrees.

2. Modeling

The NNTB approximation method is used to model the band structure of BC2N, based on Figure 1a [28]. In this regard, one-way is used to investigate how the electron–electron interaction in the structure of crystalline solves a separable Schrödinger equation, and this model is achieved by means of a time-independent Schrödinger equation. The outcome of the time-independent equation (Schrödinger) is written in the form of
E 0 = [ h ( k ) ] 0
where 0 is the wave function, and h(k) is the matrix equation as
h ( k ) = m = 1 4 H m n e i k ( d m d n )
In Equation (2),   H m n stands for the Hamiltonian matrix equation, k is the wave vector, and   d m and d n are the displacements of the mth and nth unit cells from the origin, respectively. On the other hand, for this special model, no chemical bonds between carbon–carbon bonds are considered. Now, in this case, B assigns E o C , E o B , and E o N to the onsite energies of carbon, boron, and nitrogen atoms, respectively, and the overlapping energy between the two atoms of B–C and N–C are designated by t b c = t   and t c n = t p , respectively. Thus, Equation (2) can be presented via:
h ( k )   =   [ h 1 0 0 h 1 t p h 2 h 2 h 4 t p h 3 h 2 h 4 h 1 0 0 h 1 ]
where h1 = (E0 + 2 × E0 × (cos(k.r1) + cos(k.r2)))E, h2 = t + t × (cos(k.r1)), h3 = tp + tp × (cos(k.r1)), h4 = t × (cos(k.r2)), k = kx i + ky j, r1 = 3 × a0 i, and r2 = 3a0i +   3 a0j. For interpretation purposes, it is assumed that Figure 1b is an ideal hexagonal shape and has an equal bond length of a0 = 1.446 Å between boron and nitrogen atoms. Therefore, to define the energy eigenvalues in the matrix form, the energy matrix is diagonalized by
det [ h 1 0 0 h 1 t p h 2 h 2 h 4 t p h 3 h 2 h 4 h 1 0 0 h 1 ]   =   0 .
The energy eigenvalues based on the Schrödinger equation solution lead to the following four energy dispersion relations for this particular BC2N structure due to the introduced unit cell, as shown in Figure 1a.
E 1 , 2 = ( 2 cos ( a ) + 2 cos ( b ) + 1 ) E 0 ± 2 2 ( ( 2 t 3 t cos ( a ) 2 cos ( b ) 2 + 4 t 3 t cos ( a ) 2 cos ( b ) + 4 t 3 t cos ( a ) cos ( b ) 2 + 8 t 3 t cos ( a ) cos ( b ) + 4 t 2 t 2 cos ( a ) 2 cos ( b ) + 8 t 2 t 2 cos ( a ) cos ( b ) + 2 t t 3 + t 4 + 3 t 2 t 2 + t 4 + t 4 ( cos ( a ) 4 × 4 t 4 cos ( a ) 3 + 6 t 4 cos ( a ) 2 + 4 t 4 cos ( a ) + t 4 cos ( b ) 4 + 2 t 4 cos ( b ) 2 2 t 3 t + 4 t t 3 cos ( a ) + 2 t t 3 cos ( a ) 2 t 3 + 4 t 2 t 2 cos ( b ) + 8 t 2 t 2 cos ( a ) 2 t 2 t 2 cos ( b ) 2 + 4 t 2 t 2 ( cos ( a ) ) 2 + 8 t 2 t 2 ( cos ( a ) ) 2 + 2 t t 2 ( cos ( b ) ) 2 8 t 3 t cos ( a ) + 4 t 3 t cos ( b ) + t 2 t 2 ( cos ( a ) ) 4 2 t 3 t ( cos ( a ) ) 4 8 t 3 t ( cos ( a ) ) 3 12 t 3 t ( cos ( a ) ) 2 + 2 t 4 ( cos ( a ) ) 2 ( cos ( b ) ) 2 + 2 t 4 ( cos ( a ) ) 2 ( cos ( b ) ) 2 ) 1 2 + t t + t 2 ( cos ( a ) ) 2 ) + t 2 ( cos ( b ) ) 2 + t 2 + t 2 + 2 t 2 cos ( a ) + 2 t t cos ( a ) + t t ( cos ( a ) ) 2 ) 1 2
E 3 , 4 = ( 2 × cos ( a ) + 2 cos ( b ) + 1 ) × E 0 ± 2 2 × ( t t ± 2 t 3 t cos ( a ) 2 cos ( b ) 2 + 4 t 3 t cos ( a ) 2 cos ( b ) + 4 t 3 t cos ( a ) cos ( b ) 2 + 8 t 3 t cos ( a ) cos ( b ) + 4 t 2 t 2 cos ( a ) 2 cos ( b ) + 8 t 2 t 2 cos ( a ) cos ( b ) + 2 t t 3 + t 4 + 3 t 2 t 2 + t 4 + t 4 ( cos ( a ) ) 4 × 4 t 4 ( cos ( a ) ) 3 + 6 t 4 cos ( a ) 2 + 4 t 4 cos ( a ) + t 4 cos ( b ) 4 + 2 t 4 cos ( b ) 2 2 t 3 t + 4 t t 3 cos ( a ) + 2 t t 3 cos ( a ) 2 + 4 t 2 t 2 cos ( b ) + 8 t 2 t 2 cos ( a ) 2 t 2 t 2 cos ( b ) 2 + 4 t 2 t 2 ( cos ( a ) ) 2 + 8 t 2 t 2 ( cos ( a ) ) 2 + 2 t t 2 ( cos ( b ) ) 2 8 t 3 t cos ( a ) + 4 t 3 t cos ( b ) + t 2 t 2 ( cos ( a ) ) 4 2 t 3 t ( cos ( a ) ) 4 8 t 3 t ( cos ( a ) ) 3 12 t 3 t ( cos ( a ) ) 2 + 2 t 4 ( cos ( a ) ) 2 ( cos ( b ) ) 2 + 4 t 4 cos ( a ) ( cos ( b ) ) 2 ) 1 2 + t t + t 2 ( cos ( a ) ) 2 + t 2 ( cos ( b ) ) 2 + t 2 + t 2 + 2 t 2 cos ( a ) + 2 t t cos ( a ) + t t ( cos ( a ) ) 2 ) 1 2
To elucidate the variation of the electronic band structure in BC2N, we generated an E–K relation according to Figure 1. The variation in the band gap together with solving Equations (5) and (6), which are named dispersion relation equations, is shown in Figure 2.

3. Result and Discussion

It is found that there are four answers due to the variables of wave vector and energy. The first answer is E1 (Figure 2a), which indicates a gap of 2.416 ev; the second answer is E2 with a 1.212 ev band gap (Figure 2b), the band gap due to the third answer, E3, is equal to 0.7182 ev (Figure 2c), and finally, the fourth answer E4 with a 1.3254 ev band gap (Figure 2d). The band structure that is plotted using E1, E2, E3, and E4, resulted from Equations (5) and (6), which demonstrates the suggested model, as shown in Figure 3. In the other words, the dispersion relation (Equation (6)) leads to the formation of the band structure, as shown in Figure 2 and Figure 3.
As shown in Figure 3, the possible band gaps based on the presented mode are analyzed and plotted in three-dimensional (right-hand side figures) and two-dimensional (left-hand side figures) forms in comparison with each other starting with E1 as shown in Figure 2a, followed by the accumulation of band structures E2 to E4.
In the technology application, the ultimate lattice structure is under unwanted strain and stresses that cause band energy variation. To analyze this annoying effect, one of the imaginable methods is the overlap energy variation investigation, which can be realized in the form of lattice parameter variation. Therefore, the effect of overlap energy between carbon–nitrogen and carbon–boron atoms as a disparity between the highest and lowest energy levels is considered as shown in Figure 4. In the other words, the reduction in the lattice parameter could cause an increase in the spatial overlap of the orbitals. On the other hand, the increased anti-bonding appears to be larger than the increased bonding. For simplification purposes, the hopping energy of boron–carbon and nitrogen–carbon is shown in terms of t and tp, respectively. As indicated in Figure 3, the energy overlaps between B–C, namely, tp, is fixed at 2.59 ev, and the energy overlaps between C–N, namely t, is changed. Therefore, band gap energies of about a) E1 = 2.594 ev, b) E1 = 4.266 ev, c) E1 = 5.11 ev, and d) E1 = 5.6 ev are reported. These results show that in accordance with the increase in tp, the band gap increases.
On the other hand, via changing the overlapping energy between nitrogen and carbon, namely t, an overlap between energy bands is reported. Besides, the effect of overlap energy between carbon and boron is plotted, as shown in Figure 5.
It is concluded that any variation in the overlap energy indicates a direct effect on the BC2N band structure, which can be explained by an applied quantum confinement effect in the energy matrix. In the presented model, it can be assumed that the quantum confinement effect is buried in the overlap energy variation. On the other hand, density functional theory (DFT) and tight binding (TB) methods have been employed intensively in the material property investigation [29,30,31,32]. The DFT illustrates good agreement with experimental results; however, it is computationally very expensive. Therefore, its application has been limited to a small amount of atom calculations [32]. In contrary, the TB method in the band structure calculation without self-consistent progressions needs smaller amounts of computational possessions. Therefore, TB models have been applied in the large structure (up to millions of atoms) investigations [29]. In addition, the TB model often leads to analytical terminologies that improve the logical investigation of material properties; consequently, in this research, a TB model is being implemented. On the other hand, to simulate the DFT, we use an OPENMX3.8.5-open source computer code with a lattice constant of about a = b = 5.0503 Å and c = 15.2997 Å. In addition, for the linear combination of the pseudo-atomic localized basis set, 150 Ry cutoff energy and 10 × 10 × 1 k-point are used. The basis set function is B7.0-s2p2d1, N6.0-s2p2d1, and C6.0-s2p2d1 for boron, nitrogen, and carbon, respectively. The first symbol shows the chemical name together with the cutoff radius and the initial orbitals, which is shown via the last set. All data are achieved using generalized gradient approximations (GGA) and the PBE exchange-correlation functional.
The result of the simulated structure is closer to the results of Figure 5 in response to the analytical model (comparison between Figure 5 and Figure 6), which indicates that the band structure is about 0.167 ev. The comparison study between the presented model and DFT simulation result is carried out as shown in Figure 7.
An acceptable agreement—especially on the k point—is detected, which confirms the accuracy of the proposed model. Therefore, the proposed analytical calculation could represent the prediction of electrical performance of BC2N, which also provides very fast modeling and simulation tools for band-gap investigation. On the other hand, the lattice parameter effect on the BC2N gap performance can be explored by the overlap energy gradient due to the stress and strain associated with the temperature or device fabrication limits.

4. Conclusions

BC2N is one of the most stable structures of CX (BN) y configuration. This structure is an in-between combination of hexagonal graphene and hexagonal boron nitride. The calculations showed that the band-gap energy of CX (BN) y is between pristine graphene with a band gap of zero and boron nitride with a band-gap energy of 5.5 eV. The band gap firmly depends on either the formation energy or binding energy. Conceding and doping in the graphene sheet with boron and nitrogen atoms firmly depends on its location and amount, which affects its stability. In the present study, we focused on the amazing and notable structures of the BC2N family. In this structure, no bonding is found between C–C atoms, which caused the inspiration of unique electronic and optical properties as well as its application in designing electronic devices. This study focused on the dispersion relation investigation using an analytical calculation based on the nearest neighbor’s tight-binding method. This method supposes that the band gap of the main structure of BC2N alters from 2.41 in E1 to 1.12 in E2 and from 0.7182 eV to 1.3254 eV in E4. On the other hand, using figured alternating wave vector versus energy, it is found that the changes in overlap energy between boron–carbon and nitrogen–carbon can engender the band gap.

Author Contributions

The main idea are been suggested by M.T.A. which is been simulated by A.R. In the final form, the results are been discussed and supported by S.S.R.K. and M.P. All authors were in the administration team and contributed in the manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by Ton Duc Thang University, Ho Chi Minh City, Vietnam. It is also supported by the Ministry of Education, Youth and Sports of the Czech Republic and the European Union (European Structural and Investment Funds–Operational Programme Research, Development and Education) in the frames of the project “Modular platform for autonomous chassis of specialized electric vehicles for freight and equipment transportation”, Reg. No. CZ.02.1.01/0.0/0.0/16_025/0007293.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shao, Y.; Wang, Q.; Hu, L.; Pan, H.; Shi, X. BC2N monolayers as promising anchoring materials for lithium-sulfur batteries: First-principles insights. Carbon 2019, 149, 530–537. [Google Scholar] [CrossRef]
  2. Matsuda, Y.; Morita, M.; Hanada, T.; Kawaguchi, M. A new negative electrode matrix, BC2N, for rechargeable lithium batteries. J. Power Sources 1993, 43, 75–80. [Google Scholar] [CrossRef]
  3. Keyes, R.W. Physical limits of silicon transistors and circuits. Rep. Prog. Phys. 2005, 68, 2701–2746. [Google Scholar] [CrossRef] [Green Version]
  4. Goodman, J.R. Using cache memory to reduce processor-memory traffic. ACM SIGARCH Comput. Archit. News 1983, 11, 124–131. [Google Scholar] [CrossRef]
  5. Chircu, A.M.; Kauffman, R.J. Limits to value in electronic commerce-related IT investments. J. Manag. Inf. Syst. 2000, 17, 59–80. [Google Scholar] [CrossRef]
  6. Goser, K.F.; Pacha, C.; Kanstein, A.; Rossmann, M.L. Aspects of systems and circuits for nanoelectronics. Proc. IEEE 1997, 85, 558–573. [Google Scholar] [CrossRef] [Green Version]
  7. Tarafdar, J.C.; Raliya, R. Nanotechnology; Scientific Publishers: Jodhpur, India, 2012. [Google Scholar]
  8. Brazhkin, V.; Lyapin, A. Hard and superhard carbon phases synthesised from fullerites under pressure. Сверхтвердые материалы 2012, 6, 75–105. [Google Scholar]
  9. Welser, J.; Wolf, S.A.; Avouris, P.; Theis, T. Applications: Nanoelectronics and nanomagnetics. In Nanotechnology Research Directions for Societal Needs in 2020; Springer: Dordrecht, The Netherlands, 2011; pp. 375–415. [Google Scholar]
  10. Koloor, S.S.R.; Rahimian-Koloor, S.M.; Karimzadeh, A.; Hamdi, M.; Petrů, M.; Tamin, M.N. Nano-level damage characterization of graphene/polymer cohesive interface under tensile separation. Polymers 2019, 11, 1435. [Google Scholar] [CrossRef] [Green Version]
  11. Rahimian-Koloor, S.M.; Moshrefzadeh-Sani, H.; Hashemianzadeh, S.M.; Shokrieh, M.M. The effective stiffness of an embedded graphene in a polymeric matrix. Curr. Appl. Phys. 2018, 18, 559–566. [Google Scholar] [CrossRef]
  12. Geim, A.K.; Novoselov, K.S. The rise of graphene. In Nanoscience and Technology: A Collection of Reviews from Nature Journals; World Scientific: London, UK, 2010; pp. 11–19. [Google Scholar]
  13. Bertolazzi, S.; Bondavalli, P.; Roche, S.; San, T.; Choi, S.Y.; Colombo, L.; Bonaccorso, F.; Samorì, P. Nonvolatile Memories Based on Graphene and Related 2D Materials. Adv. Mater. 2019, 31, 1806663. [Google Scholar] [CrossRef] [Green Version]
  14. Tang, C.; Bando, Y.; Sato, T.; Kurashima, K. A novel precursor for synthesis of pure boron nitride nanotubes. Chem. Commun. 2002, 12, 1290–1291. [Google Scholar] [CrossRef] [PubMed]
  15. Chimene, D.; Alge, D.L.; Gaharwar, A.K. Two-dimensional nanomaterials for biomedical applications: Emerging trends and future prospects. Adv. Mater. 2015, 27, 7261–7284. [Google Scholar] [CrossRef] [PubMed]
  16. Xu, M.; Liang, T.; Shi, M.; Chen, H. Graphene-like two-dimensional materials. Chem. Rev. 2013, 113, 3766–3798. [Google Scholar] [CrossRef] [PubMed]
  17. Pierret, R.F. Semiconductor Device Fundamentals; Pearson Education India: Noida, India, 1996. [Google Scholar]
  18. Yoder, M.N. Wide bandgap semiconductor materials and devices. IEEE Trans. Electron Devices 1996, 43, 1633–1636. [Google Scholar] [CrossRef]
  19. Rani, P.; Jindal, V. Designing band gap of graphene by B and N dopant atoms. RSC Adv. 2013, 3, 802–812. [Google Scholar] [CrossRef] [Green Version]
  20. Mazzoni, M.S.; Nunes, R.; Azevedo, S.; Chacham, H. Electronic structure and energetics of BxCyNz layered structures. Phys. Rev. B 2006, 73, 073108. [Google Scholar] [CrossRef]
  21. Golberg, D.; Bando, Y.; Tang, C.; Zhi, C. Boron nitride nanotubes. Adv. Mater. 2007, 19, 2413–2432. [Google Scholar] [CrossRef]
  22. Kawaguchi, M. B/C/N materials based on the graphite network. Adv. Mater. 1997, 9, 615–625. [Google Scholar] [CrossRef]
  23. Ci, L.; Song, L.; Jin, C.; Jariwala, D.; Wu, D.; Li, Y.; Srivastava, A.; Wang, Z.; Storr, K.; Balicas, L. Atomic layers of hybridized boron nitride and graphene domains. Nat. Mater. 2010, 9, 430. [Google Scholar] [CrossRef]
  24. Park, C.; Chadi, D. Stability of deep donor and acceptor centers in GaN, AlN, and BN. Phys. Rev. B 1997, 55, 12995. [Google Scholar] [CrossRef]
  25. Liu, A.Y.; Wentzcovitch, R.M.; Cohen, M.L. Atomic arrangement and electronic structure of BC2N. Phys. Rev. B 1989, 39, 1760. [Google Scholar] [CrossRef] [PubMed]
  26. Terrones, M.; Botello-Méndez, A.R.; Campos-Delgado, J.; López-Urías, F.; Vega-Cantú, Y.I.; Rodríguez-Macías, F.J.; Elías, A.L.; Munoz-Sandoval, E.; Cano-Márquez, A.G.; Charlier, J.-C. Graphene and graphite nanoribbons: Morphology, properties, synthesis, defects and applications. Nano Today 2010, 5, 351–372. [Google Scholar] [CrossRef]
  27. Yao, B.; Chen, W.; Liu, L.; Ding, B.; Su, W. Amorphous B–C–N semiconductor. J. Appl. Phys. 1998, 84, 1412–1415. [Google Scholar] [CrossRef]
  28. Tang, Q.; Zhou, Z. Graphene-analogous low-dimensional materials. Prog. Mater. Sci. 2013, 58, 1244–1315. [Google Scholar] [CrossRef]
  29. Tran, V.T.; Saint-Martin, J.; Dollfus, P.; Volz, S. Third nearest neighbor parameterized tight binding model for graphene nano-ribbons. AIP Adv. 2017, 7, 075212. [Google Scholar] [CrossRef] [Green Version]
  30. Simão, C.D.; Reparaz, J.S.; Wagner, M.R.; Graczykowski, B.; Kreuzer, M.; Ruiz-Blanco, Y.B.; García, Y.; Malho, J.-M.; Goñi, A.R.; Ahopelto, J.; et al. Optical and mechanical properties of nanofibrillated cellulose: Toward a robust platform for next-generation green technologies. Carbohydr. Polym. 2015, 126, 40. [Google Scholar] [CrossRef] [Green Version]
  31. Bena, C.; Montambaux, G. Remarks on the tight-binding model of graphene. New J. Phys. 2009, 11, 095003. [Google Scholar] [CrossRef]
  32. Papaconstantopoulos, D.A.; Mehl, M.J. The Slater–Koster tight-binding method: A computationally efficient and accurate approach. J. Phys. Condens. Matter 2003, 15, R413. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Unit cell of BC2N; (b) ultimate Structure of BC2N.
Figure 1. (a) Unit cell of BC2N; (b) ultimate Structure of BC2N.
Materials 13 01026 g001
Figure 2. (a) One of the acceptable respondents from the energy determinant with a three-dimensional band structure scheme mainly consisting of BC2N; (b) The band gap resulting from the Equation (3) is due to E1 is 2.416 ev; (c) The band gap corresponding to Equation (4), namely E2, is about 1.212 ev; (d) On the other hand, the band gap resulting from Equation (5) E3 is 0.7182 ev. Finally, as shown, through solving Equation (6), it is found that the fourth possible band gap E4 is 1.3254 ev.
Figure 2. (a) One of the acceptable respondents from the energy determinant with a three-dimensional band structure scheme mainly consisting of BC2N; (b) The band gap resulting from the Equation (3) is due to E1 is 2.416 ev; (c) The band gap corresponding to Equation (4), namely E2, is about 1.212 ev; (d) On the other hand, the band gap resulting from Equation (5) E3 is 0.7182 ev. Finally, as shown, through solving Equation (6), it is found that the fourth possible band gap E4 is 1.3254 ev.
Materials 13 01026 g002
Figure 3. (a) Dispersion relation plotted for E1, which indicates a band gap = 2.312 ev; (b) dispersion relation plotted for both E1, and E2, in which the minimum band gap is found to be 1.1988 ev; (c) dispersion relations for both E1 and E3 are compared, and a minimum of about 0.7182 ev is reported; (d) dispersion relations for both E1 and E4 are compared, and a minimum band gap of about 1.304 ev is shown; (e) dispersion relation results for E1, E2, E3, and E4 are compared, and a minimum band gap of about 0.7182 ev is addressed. Note: In the right-hand side of each figure, a two-dimensional view is presented with the same label.
Figure 3. (a) Dispersion relation plotted for E1, which indicates a band gap = 2.312 ev; (b) dispersion relation plotted for both E1, and E2, in which the minimum band gap is found to be 1.1988 ev; (c) dispersion relations for both E1 and E3 are compared, and a minimum of about 0.7182 ev is reported; (d) dispersion relations for both E1 and E4 are compared, and a minimum band gap of about 1.304 ev is shown; (e) dispersion relation results for E1, E2, E3, and E4 are compared, and a minimum band gap of about 0.7182 ev is addressed. Note: In the right-hand side of each figure, a two-dimensional view is presented with the same label.
Materials 13 01026 g003aMaterials 13 01026 g003b
Figure 4. (a) t = 2.59 ev, tp = 1 ev; (b) for t = 2.59 ev, tp = 2 ev; (c) t = 2.59 ev, tp = 2.5 ev; (d) t = 2.5 ev, tp = 2.83 ev.
Figure 4. (a) t = 2.59 ev, tp = 1 ev; (b) for t = 2.59 ev, tp = 2 ev; (c) t = 2.59 ev, tp = 2.5 ev; (d) t = 2.5 ev, tp = 2.83 ev.
Materials 13 01026 g004
Figure 5. (a) t = 1 ev, tp = 1 ev; (b) t = 1 ev, tp = 2.83 ev; (c) t = 2 ev, tp = 2.83 ev; (d) t = 2.5 ev, tp = 2.83 ev.
Figure 5. (a) t = 1 ev, tp = 1 ev; (b) t = 1 ev, tp = 2.83 ev; (c) t = 2 ev, tp = 2.83 ev; (d) t = 2.5 ev, tp = 2.83 ev.
Materials 13 01026 g005
Figure 6. The band structure of BC2N.
Figure 6. The band structure of BC2N.
Materials 13 01026 g006
Figure 7. The comparison between the presented model (red dots) and density functional theory (DFT) simulation (black stars).
Figure 7. The comparison between the presented model (red dots) and density functional theory (DFT) simulation (black stars).
Materials 13 01026 g007

Share and Cite

MDPI and ACS Style

Ahmadi, M.T.; Razmdideh, A.; Rahimian Koloor, S.S.; Petrů, M. Carbon-Based Band Gap Engineering in the h-BN Analytical Modeling. Materials 2020, 13, 1026. https://doi.org/10.3390/ma13051026

AMA Style

Ahmadi MT, Razmdideh A, Rahimian Koloor SS, Petrů M. Carbon-Based Band Gap Engineering in the h-BN Analytical Modeling. Materials. 2020; 13(5):1026. https://doi.org/10.3390/ma13051026

Chicago/Turabian Style

Ahmadi, Mohammad Taghi, Ahmad Razmdideh, Seyed Saeid Rahimian Koloor, and Michal Petrů. 2020. "Carbon-Based Band Gap Engineering in the h-BN Analytical Modeling" Materials 13, no. 5: 1026. https://doi.org/10.3390/ma13051026

APA Style

Ahmadi, M. T., Razmdideh, A., Rahimian Koloor, S. S., & Petrů, M. (2020). Carbon-Based Band Gap Engineering in the h-BN Analytical Modeling. Materials, 13(5), 1026. https://doi.org/10.3390/ma13051026

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