 Research
 Open Access
 Published:
Geometric and electronic structures of monolayer hexagonal boron nitride with multivacancy
Nano Convergence volume 4, Article number: 13 (2017)
Abstract
Hexagonal boron nitride (hBN) is an electrical insulator with a large band gap of 5 eV and a good thermal conductor of which melting point reaches about 3000 °C. Due to these properties, much attention was given to the thermal stability rather than the electrical properties of hBN experimentally and theoretically. In this study, we report calculations that the electronic structure of monolayer hBN can be influenced by the presence of a vacancy defect which leads to a geometric deformation in the hexagonal lattice structure. The vacancy was varied from mono to trivacancy in a supercell, and different defective structures under the same vacancy density were considered in the case of an odd number of vacancies. Consequently, all cases of vacancy defects resulted in a geometric distortion in monolayer hBN, and new energy states were created between valence and conduction band with the Fermi level shift. Notably, B atoms around vacancies attracted one another while repulsion happened between N atoms around vacancies, irrespective of vacancy density. The calculation of formation energy revealed that multivacancy including more Bvacancies has much lower formation energy than vacancies with more Nvacancies. This work suggests that multivacancy created in monolayer hBN will have more Bvacancies and that the presence of multivacancy can make monolayer hBN electrically conductive by the new energy states and the Fermi level shift.
Background
Boron nitride (BN) is a material with a superior thermal stability [1,2,3] and exists in various crystal structures such as sphalerite, wurtzite, and hexagonal structures. Among them, hexagonal BN (hBN) is a layered material which has the exactly same structure to graphite. In the case of graphite, C atoms are arranged in the hexagonal lattice structure. Similarly, B and N atoms are orderly occupied in the hexagonal structure of hBN layers. Recently, graphene, or a layer of graphite, [4] has attracted attention due to the electrical properties such as ambipolar field effect [4], high carrier mobility [5,6,7], and quantum Hall effect at room temperature [8,9,10,11,12,13]. And, huge efforts have been made to tune the electronic structure of graphene theoretically [14,15,16] and experimentally [17,18,19]. Different from graphene, hBN is electrically an insulator with a wide band gap of about 5 eV [20,21,22,23,24]. Due to this, less effort has been made to investigate the electrical properties of hBN. Furthermore, converting the electrical characteristics of hBN from an insulator to an electrical conductor is a significant and interesting task.
Several theoretical calculations have been conducted to predict the electrical properties of monolayer hBN by structural modification. According to the calculations, the electronic structure can be affected by substituting B or N atom with different kinds of atoms [15, 21, 24, 25] or by vacancies in monolayer hBN by removing B or N atom [20, 21, 24]. Especially, in the case of theoretical works on a vacancy in monolayer hBN, most of the studies focus on monovacancy such as B or Nvacancy rather than multivacancy defects. However, with the advance in etching hBN [26, 27], creating multivacancy in the layered structure is not a technical barrier any more. Hence, investigating the effect of multivacancy on the electronic structure of hBN layer is a significant work. Specifically, it is because hBN is composed of two kinds of atoms different from graphene and different configurations of multivacancy are formed even under the same vacancy density where the number of vacancy is odd. Therefore, each vacancy configuration will have different electronic structures as well as geometric deformations. However, few efforts have been made to study a change in the electronic structure of hBN by multivacancy defects.
In this study, we theoretically investigated the effect of multivacancy on the geometric and electronic structures of monolayer hBN. The density of multivacancy was varied from mono to trivacancy by removing B or N atoms from monolayer hBN. As mentioned above, a lattice structure of hBN with B and Nvacancy should be considered separately in the case of monovacancy. Trivacancy is also placed in the same situation where 2Bvacancies and Nvacancy should be considered with the vacancy configuration of 2Nvacancies and Bvacancy despite the same vacancy density. Therefore, we conducted calculations on each configuration of defected hexagonal lattice of hBN when the number of vacancy is odd (mono and trivacancy). Also, the formation energy (E_{form}) of hBN with different densities of vacancy was also calculated to clarify which vacancy configuration is preferred to be formed in the structure of hBN monolayer.
Computational method
Our calculations based on dispersioncorrected density functional theory with the local density approximation were implemented in the PWSCF code of the QUANTUM ESPRESSO package. Perdew–Zunger parameterization and norm conserving Troullier–Martins pseudopotentials were used for the electron exchange–correlation interaction. Also, spinpolarization was considered in all cases of the calculations. We employed a 5 × 5 supercell structure of monolayer hBN which contains 50 atoms. And, a defected structure was realized by removing B and N atoms according to the number of a vacancy. Triclinic supercell of 12.5 × 12.5 Å with periodic boundary conditions was used to simulate pristine and defected monolayer hBN. A vacuum region of 15 Å was added to avoid the interactions between periodic images. The equilibrium atomic position was determined by relaxing all B and N atoms in pristine and defected structures. The Brillouin zone of the supercell was sampled by a 12 × 12 × 1 Monkhorst–Pack grid. A plane wave basis set with a kinetic cutoff energy of 30 Ry was used.
The defective hBN was realized by removing B or N atom one by one from pristine hBN structure, depending on the number of vacancy. Figure 1 presents the schematic illustration of modeling hBN with mono, di, and trivacancy. In the case of monovacancy, two cases were considered by eliminating the B or N atom in the position of A or B (hereafter, 1V_{B} and 1V_{N}), respectively. B and N atoms corresponding to A and B sites were removed at the same time to model monolayer hBN with divacancy (hereafter, 2V_{BN}). Two defective configurations of trivacancy were realized by removing atoms in the position of C–A–B and A–B–D (hereafter, 3V_{2NB} and 3V_{2BN}), respectively.
The formation energy of each structure was calculated by using the following equation:
where, E_{V} and E are the total energy of defective and pristine hBN, respectively. And, n and n_{V} are the number of atoms in pristine structure and vacancies in defective hBN.
Results and discussion
Monovacancy
Before investigating the geometric and electronic structures of monolayer hBN with monovacancy, it is necessary to study pristine hBN in order to distinguish geometric and electronic differences between pristine and defective hBN. Thus, the geometric and electronic structures of pristine lattice structure were calculated. Figure 2 shows a relaxed structure, isodensity plot, and electronic band structure of pristine hBN. As shown in Fig. 2a, the bonding length between B and N atoms was calculated to be 1.45 Å from the relaxed structure. As marked with a triangle in the figure, the distances between B atoms and N atoms are all of the same length (2.50 Å). The isodensity plot in Fig. 2b shows that when a vacancy does not exist in monolayer hBN, electron density is evenly distributed throughout the supercell. Figure 2c displays the electronic band structure of pristine hBN which was calculated along the ΓKMΓ path in the Brillouin zone. The band structure reveals that pristine hBN has an indirect band gap of 4.3 eV in the electronic structure, which is consistent with previous calculations [25, 28].
Different from graphene, hBN can have two kinds of defective structures in the case of monovacancy (1V_{B} and 1V_{N}), depending on a vacancy element. Thus, both different configurations should be separately considered to study the effect of monovacancy on the geometric and electronic band structures. Figure 3a and b display relaxed supercell structures with B and N vacancy, respectively. To evaluate the degree of distortion, the distance between atoms around the vacancy site was measured in each case. For this, the atoms around monovacancy were connected by lines which lead to a triangle. When 1V_{B} is created in the structure as seen in Fig. 3a, all distances between the atoms are measured to be 2.61 Å which is greater than 2.50 Å measured from pristine hBN. The relaxed structure of the defected hBN and the measured distance are consistent with the results reported by previous works [21, 24]. The elongated distance indicates that the presence of monovacancy distorted a geometric structure of monolayer hBN. And, the bonding lengths between the atoms at the edge site around B vacancy are changed to 1.41 and 1.44 Å. This also means that a geometric distortion occurred by Bvacancy in comparison with pristine lattice structure. Any planar distortion perpendicular to the monolayer did not happen after relaxation. In the case of 1V_{N}, a similar shape of deformation to Bvacancy appears in the relaxed structure, and the distance between B atoms is measured to be 2.31 Å. This relaxed structure and the distance are also consistent with the previous works [21, 24]. The bonding lengths of atoms existing at the edge of monovacancy ranges from 1.44 to 1.45 Å, indicating that the presence of Nvacancy causes a geometric distortion in the supercell. Also, there was no planar deformation to the vertical axis of hBN lattice structure. Interestingly, the distance between B atoms in Fig. 3b is shorter than the length between N atoms in Fig. 3a. This reason can be found in the isodensity plot displayed in Fig. 3c and d. When B atom is missing in the supercell, the repulsion between N atoms around the vacancy breaks out as shown in Fig. 3c. Meanwhile, B atoms around Nvacancy attract one another, resulting in a shortened distance of 2.31 Å and an elongated bonding length of 1.45 Å. Therefore, the shortened distance can be attributed to the attraction between B atoms around Nvacancy. This different situation gives rise to different electronic structures as shown in Fig. 3e and f. In the case of 1V_{B}, a gap in the band structure has a similar value to that of pristine monolayer hBN. However, the presence of 1V_{B} shifts the Fermi level downward valence band. In the case of 1V_{N}, a new energy state is created along the Fermi level as shown in Fig. 3f, which might help electron jump to conduction band.
Divacancy
Different from an oddnumber vacancy, divacancy in hBN has only one defective configuration which includes all elements. B and N atoms positioned at A and B site in Fig. 1 were removed to model divacancy in the supercell. Figure 4a displays a relaxed structure of monolayer hBN with 2V_{BN}. To examine the degree of geometric distortion, the distance between atoms at the edge site around the vacancies was measured as marked with lines shown in Fig. 4a. As a result, the distance between B atoms around 2V_{BN} was measured to be 1.99 Å while the distance between N atoms at the edge was calculated to be 2.33 Å. These values are less than the distance calculated in pristine lattice structure. In the case of the distance between B and N atoms around 2V_{BN}, the length was calculated to be 3.07 Å in all cases as seen in Fig. 4a. It can be understood that this is due to the repulsion by the localized electrons around N atoms and the attraction between B atoms, resulting in an elongated distance of 3.07 Å greater than 2.50 Å. No planar distortion was observed toward the direction perpendicular to the monolayer. Also, the divacancy influences on the bonding length between atoms at the edge around the vacancies, where the bonding length is distributed from 1.40 to 1.5 Å. Consequently, this leads to a distortion in the geometric structure, and changes the electron density distribution as presented in Fig. 4b. The distribution also shows that a weak covalent bond is formed between B atoms due to the relatively high electron density. This apparently explains why the distance between B atoms was shortened to 1.99 Å in monolayer hBN with 2V_{BN}. The structural change and the presence of divacancy lead to an alteration in the electronic structure. As displayed in Fig. 4c, two new energy states, which do not show in the band structure of pristine hBN, appear between the valence and conduction band. The new energy states can contribute to electron jumps to conduction band.
Trivacancy
Similar to monovacancy, monolayer hBN with trivacancy is present in the form of two defective configurations such as 3V_{2BN} and 3V_{2NB}. Figure 5a shows the relaxed hexagonal lattice structure of monolayer hBN with 3V_{2BN}. To investigate the degree of distortion by the vacancies, the atoms at the edge around the vacancies were connected by lines. As seen in the figure, the lines form a pentagon and have different lengths of 2.43, 2.62, and 3.15 Å. The different distances indicate that the structure was deformed by the presence of 3V_{2BN} in comparison with the pristine structure. This is also supported that the bonding lengths between atoms at the edge sites range from 1.34 to 1.47 Å, which are deviated from the value of 1.45 Å calculated from pristine hBN. Figure 5b displays a deformed structure of monolayer hBN by 3V_{2NB} where the atoms at the edge sites are also connected by lines. The sides of pentagon have different values of 2.06, 2.52, and 2.98 Å, and the bonding lengths of atoms at the edge are in the range between 1.38 and 1.51 Å. As seen in mono and divacancy, the distance between B atoms around the vacancies becomes shorter than 2.50 Å corresponding to the distance between B atoms in pristine hBN. The isodensity plot shown in Fig. 5c and d can explain the reason why the distance between B atoms became short by the vacancies. In the case of 3V_{2BN}, N atoms around the vacancies in Fig. 5c have much higher electron density than B atoms while almost few electrons exist around B atoms at the edge sites in Fig. 5c. However, in the case of 3V_{2NB}, B atoms at the edge site have more electron density and form a weak covalent bond by attraction between two B atoms as seen in Fig. 5d. Consequently, these different electron densities in the lattice structure lead to different electronic structures. As displayed in Fig. 5e and f, new energy states are generated between valence and conduction band in both cases of 3V_{2BN} and 3V_{2NB}. However, in the case of 3V_{2BN}, the new states are evenly distributed in the gap defined by valence and conduction band, and the Fermi level exists between the states. Meanwhile, in the case of 3V_{2NB}, the new states are relatively unevenly positioned, and the Fermi level exists away from the new states, compared with the electronic structure created by 3V_{2BN}. Nevertheless, it can be expected that the new energy states can help electrons jump from valence to conduction band, converting monolayer hBN to an electrically conductive material.
E_{form} of the vacancies
It has been revealed that defective hBN structures have different geometric deformations and electronic structures even under the same vacancy density. This means that each structure has unique E_{form}, depending on the configuration and the number of a vacancy defect. Therefore, calculating E_{form} of each defective structure can give information about which vacancy is an energetically preferable structure. Figure 6 exhibits a plot of E_{form} with respect to the number of a vacancy. In the case of monovacancy, E_{form} of 1V_{B} has much lower value than that of 1V_{N}, meaning that Bvacancy is an energetically preferable defect rather than Nvacancy. The E_{form} of 2V_{BN} has a positive value of 8.9 eV, suggesting that divacancy is not preferable to be energetically formed. A similar trend seen in monovacancy is found in the case of trivacancy. While E_{form} of 3V_{2BN} has a negative value of −75.2 eV, 3V_{2NB} has a positive value of 98.1 eV. This means that trivacancy is energetically preferable to be present in the configuration of 3V_{2BN} rather than of 3V_{2NB}. From this result, it can be said that when hBN has the same or greater number of Bvacancy than that of Nvacancy out of the total number of vacancies, the defective hBN structure is energetically preferable to be formed. Therefore, it is predicted that multivacancy created in the monolayer hBN will have more Bvacancies than Nvacancies.
Conclusions
In this work, we theoretically investigated the effect of multivacancy on the geometric deformation and the electronic structure of monolayer hBN. The presence of vacancy resulted in a geometric deformation which leads to a change in the electronic structure of monolayer hBN irrespective of the number of a vacancy. However, a planar deformation did not break out at any case in the relaxed structure. Specifically, regardless of vacancy density, the repulsion between N atoms around vacancies appeared due to delocalized electrons. And, attraction occurred between B atoms at the edge of vacancies. The presence of vacancies and the geometric deformation in the hexagonal lattice structure created new energy states in the band gap and shifted the Fermi level. Such a change by vacancies can help electrons jump from valence to conduction band, opening a possibility to convert hBN into an electrical conductor.
The calculation of E_{form} revealed that multivacancy with more B vacancies is energetically preferable to be formed even under the same vacancy density. The result suggests that monolayer hBN will have more Bvacancies than Nvacancies when multivacancy is created in the hexagonal lattice structure.
Abbreviations
 BN:

boron nitride
 hBN:

hexagonal boron nitride
 E_{form} :

formation energy
References
 1.
N. Yang, C. Xu, J. Hou, Y. Yao, Q. Zhang, M.E. Grami, L. He, N. Wang, X. Qu, Preparation and properties of thermally conductive polyimide/boron nitride composites. RSC Adv. 6, 18279–18287 (2016)
 2.
T. Huang, X. Zeng, Y. Yao, R. Sun, F. Meng, J. Xu, C. Wong, Boron nitride@graphene oxide hybrids for epoxy composites with enhanced thermal conductivity. RSC Adv. 6, 35847–35854 (2016)
 3.
J. Yu, X. Huang, C. Wu, X. Wu, G. Wang, P. Jiang, Interfacial modification of boron nitride nanoplatelets for epoxy composites with improved thermal properties. Polymer 53, 471–480 (2012)
 4.
K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, Y. Zhang, S.V. Dubonos, I.V. Grigorieva, A.A. Firsov, Electric field effect in atomically thin carbon films. Science 306, 666–669 (2004)
 5.
K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. Katsnelson, I.V. Grigorieva, S.V. Dubonos, A.A. Firsov, Twodimensional gas of massless dirac fermions in graphene. Nature 438, 197–200 (2005)
 6.
S.V. Morozov, K.S. Novoselov, M.I. Katsnelson, F. Schedin, D.C. Elias, J.A. Jaszczak, A.K. Geim, Giant intrinsic carrier mobilities in graphene and its bilayer. Phys. Rev. Lett. 100, 016602 (2008)
 7.
M. Han, B. Ozyilmaz, Y. Zhang, P. JarilloHerero, P. Kim, Electronic transport measurements in graphene nanoribbons. Phys. Status Solidi (b) 244, 4134–4137 (2007)
 8.
K.S. Novoselov, E. McCann, S.V. Morozov, V.I. Falko, M.I. Katsnelson, U. Zeitler, D. Jiang, F. Schedin, A.K. Geim, Unconventional quantum Hall effect and Berry’s phase of 2π in bilayer graphene. Nat. Phys. 2, 177–180 (2006)
 9.
Z. Jiang, Y. Zhang, Y.W. Tan, H.L. Stormer, P. Kim, Quantum Hall effect in graphene. Solid State Commun. 143, 14–19 (2007)
 10.
Z. Jiang, Y. Zhang, H.L. Stormer, P. Kim, Quantum Hall states near the chargeneutral dirac point in graphene. Phys. Rev. Lett. 99, 106802 (2007)
 11.
Y. Zhang, Y.W. Tan, H.L. Stormer, P. Kim, Experimental observation of the quantum Hall effect and Berry’s phase in graphene. Nature 438, 201–204 (2005)
 12.
K.S. Novoselov, Z. Jiang, Y. Zhang, S.V. Morozov, H.L. Stormer, U. Zeitler, J.C. Maan, G.S. Boebinger, P. Kim, A.K. Geim, Roomtemperature quantum Hall effect in graphene. Science 315, 1379 (2007)
 13.
B. Özyilmaz, P. JarilloHerrero, D. Efetov, D.A. Abanin, L.S. Levitov, P. Kim, Electronic transport and quantum Hall effect in bipolar graphene pnp junctions. Phys. Rev. Lett. 99, 166804 (2007)
 14.
L. Liu, Z. Shen, Bandgap engineering of graphene: a density functional theory study. Appl. Phys. Lett. 95, 252104 (2009)
 15.
J. da Rocha Martins, H. Chacham, Disorder and segregation in B−C−N graphenetype layers and nanotubes. ACS Nano 5(1), 385–393 (2010)
 16.
S.M. Kozlov, F. Viñes, A. Görling, Bandgap engineering of graphene by physisorbed adsorbates. Adv. Mater. 23, 2638–2643 (2011)
 17.
R. Balog, B. Jørgensen, L. Nilsson, M. Andersen, E. Rienks, M. Bianchi, M. Fanetti, E. Lægsgaard, A. Baraldi, S. Lizzit, Z. Sljivancanin, F. Besenbacher, B. Hammer, T.G. Pedersen, P. Hofmann, L. Hornekær, Bandgap opening in graphene induced by patterned hydrogen adsorption. Nat. Mater. 9, 315–319 (2010)
 18.
H. Chang, J. Cheng, X. Liu, J. Gao, M. Li, J. Li, X. Tao, F. Ding, Z. Zheng, Facile synthesis of widebandgap fluorinated graphene semiconductors. Chem. Eur. J. 17, 8896–8903 (2011)
 19.
Y.C. Chen, D.G. de Oteyza, Z. Pedramrazi, C. Chen, F.R. Fischer, M.F. Crommie, Tuning the band gap of graphene nanoribbons synthesized from molecular precursors. ACS Nano 7, 6123–6128 (2013)
 20.
M.S. Si, J.Y. Li, H.G. Shi, X.N. Niu, D.S. Xue, Divacancies in graphitic boron nitride sheets. EPL 86, 46002 (2009)
 21.
S. Azevedo, J.R. Kaschny, C.M.C. de Castilho, F. de Brito Mota, Electronic structure of defects in a boron nitride monolayer. Eur. Phys. J. B 67, 507–512 (2009)
 22.
L.C. Yin, H.M. Cheng, R. Saito, Triangle defect states of hexagonal boron nitride atomic layer: density functional theory calculations. Phys. Rev. B 81, 153407 (2010)
 23.
M. Topsakal, E. Aktürk, S. Ciraci, Firstprinciples study of two and onedimensional honeycomb structures of boron nitride. Phys. Rev. B 79, 115442 (2009)
 24.
B. Huang, H. Lee, Defect and impurity properties of hexagonal boron nitride: a firstprinciples calculation. Phys. Rev. B 86, 245406 (2012)
 25.
N. Berseneva, A. Gulans, A.V. Krasheninnikov, R.M. Nieminen, Electronic structure of boron nitride sheets doped with carbon from firstprinciples calculations. Phys. Rev. B 87, 035404 (2013)
 26.
P. Sutter, J. Lahiri, P. Albrecht, E. Sutter, Chemical vapor deposition and etching of highquality monolayer hexagonal boron nitride films. ACS Nano 5, 7303–7309 (2011)
 27.
N. Alem, R. Erni, C. Kisielowski, M.D. Rossell, W. Gannett, A. Zettl, Atomically thin hexagonal boron nitride probed by ultrahighresolution transmission electron microscopy. Phys. Rev. B 80, 155425 (2009)
 28.
K.T. Park, K. Terakura, N. Hamada, Bandstructure calculations for boron nitrides with three different crystal structures. J. Phys. C: Solid State Phys. 20, 1241–1251 (1987)
Authors’ contributions
DK conducted theoretical calculations and manuscript preparation. HK, MS, SL, and SYL contributed to discussion and manuscript preparation. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The authors have no data to share since all data are shown in the submitted manuscript.
Funding
S. Lee and M. W. Song gratefully acknowledge support from Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2016R1D1A1B03932999).
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Kim, DH., Kim, HS., Song, M.W. et al. Geometric and electronic structures of monolayer hexagonal boron nitride with multivacancy. Nano Convergence 4, 13 (2017). https://doi.org/10.1186/s4058001701070
Received:
Accepted:
Published:
Keywords
 Boron nitride
 Vacancy
 Defect
 Deformation
 Band structure