Abstract

In this computational study, a novel borophene polymorph, Ogee-Borophene, characterized by irregular decagon-shaped hollows was reported. The structure involves a deviation from hexagonal configurations found in all other known borophene polymorphs. The decagon-shaped hollow is highly related to the anisotropy of the structure, which results in three types of boron atoms with different electronic properties in the structure. In the study, the electronic structure and density of states of this novel structure were investigated with the help of density functional theory calculations. The electronic structure of Ogee-Borophene shows Dirac cone formations near the Fermi level. The discovery of a novel borophene polymorph, Ogee-Borophene, with irregular-shaped hollows represents a different point of view in the 2D materials field. The unique electronic properties of this material suggest that Ogee-Borophene has the potential to be used in a variety of applications, including transistors, and selective sensor applications without the need for additional doping.

1. Introduction

With the discovery of graphene, especially in the last decade, interest in other 2-dimensional (2D) materials has increased [1, 2]. In the studies that emerged as a result of this interest, it was seen that 2D materials have different and rich properties from each other [3]. When these 2D materials are classified, it is seen that a few species such as silicene, germanene, arsenene, phosphorene, antimonene, stanene, bismuthene, tellurene, and borophene are composed of a single element [412].

For the last 20 years, 2D boron nanosheets, 0D nanoclusters, and 1D nanotubes have attracted great interest as research topics [1325]. Borophene, which is a 2D allotrope of boron, was synthesized in 2015 by Mannix et al. [26] on an Ag(111) substrate under vacuum conditions. The experimental realization of borophene attracted tremendous interest due to its extraordinary electronic and mechanical properties [27]. In terms of electronic structure, while most 2D materials are as investigated as insulators, semimetals, or semiconductors, borophene, and other boron allotropes are found to be metals [28]. Therefore, borophene and other boron allotropes provide a strong base for future nanoelectronic applications as a 2D metal [27]. In terms of mechanical properties, borophene shows an anisotropic structure and structural polymorphism, which also gives tunable electronic properties to borophene and its polymorphs [29]. Normally, it can be predicted that flat boron sheets may form a triangular bonding configuration. However, this configuration is unstable due to excess electrons occupying the antibonding orbitals [30]. The structure of borophene can only be stabilized by the introduction of vacancies. A hollow-hexagon-shaped vacancy-mediated superstructural motif structural polymorphism is suggested and the structural stability and diversity of elemental boron layers are evaluated by treating them as a pseudo alloy B1–xVax, where Va is a vacancy in the close-packed triangular B lattice [29]. Therefore, many different studies report structures with different hexagonal hollow concentrations and are investigated to find stable free-standing and grown-on structures [3135].

With the increasing experimental studies, it is found that the substrate is affecting the ground state of the 2D borophenes [32]. The preferred 2D-grown polymorph is found to be a flat v1/6 type sheet on Ag, Cu, and Ni substrates. But, on an Au substrate, a v1/12 type buckled sheet is found to be preferred [36]. On experimental Ag(111) growths and its calculations [36], hollow-hexagons chains are found to be aligned as separated by hexagonal boron rows (named ) or narrower zigzag boron rows separated by hole arrays (named ), which shows the alignments of these vacancies can be changed by different substrates and it is important to study other possible alignments and vacancy types. However, to our knowledge, no study has been reported about the investigation of any other hollow vacancies in borophenes.

In this study, we utilize first-principle density functional theory (DFT) to suggest a new polymorph with irregular decagon-shaped hollows. Because this vacancy is not regular-shaped, which is more like an ellipse than a circle, the irregularity of the vacancy’s shape will induce different electronic properties to different atoms in the vicinity of the vacancy. Our study demonstrates the different electronic properties of different boron atoms provide tuning or selective interaction functionality to this newly suggested 2D metallic polymorph with the shape of an Ogee pattern.

2. Computational Details

The structure considered in this work was built and relaxed using the Atomic Simulation Environment (ASE) software package [37]. All the calculations were carried out using DFT where the wave functions are expanded in plane waves as implemented in the grid-based projector-augmented wave (GPAW) code [38, 39]. In calculations, we have used the gpaw-tools scripts to interact with ASE and GPAW, create, read, manipulate, and analyze output data and other user interface needs [40]. The crystal structures and electron density maps were visualized using the VESTA software [41]. The exchange–correlation energy is described by the generalized gradient approximation (GGA) using the Perdew–Burke–Ernzerhof (PBE) functional. A vacuum space of 20 Å in the normal orientation of the Ogee-Borophene sheet was used to avoid interactions between the neighboring sheets. The calculations were performed on a C2/m (12) space group cell. The plane-wave cutoff energy was 700 eV. The k-point-mesh for Brillouin zone sampling of primitive cells, based on the Monkhorst–Pack scheme was 9 × 9 × 1. All atoms were allowed to relax until forces were below 0.05 eV Å−1. Phonon dispersion calculation with a 3 × 3 × 1 supercell was performed using Phonopy code [42].

3. Results and Discussion

Figure 1 shows the atomic structure of Ogee-Borophene used in the calculations for a unit cell. The unit cell of Ogee-Borophene consists of eight boron atoms. Figure 1(b) shows the side view of the structure shows an S-shaped buckled behavior. The calculated bond lengths and angles are listed in Tables 1 and 2, respectively. Because of the irregular bond lengths distributed in the structure, the hollow has 10 sides instead of six. Therefore, it is clear that the structure can be suggested to be a polymorph of a Borophene with irregular decagon-shaped hollows.

Figure 2 shows a 3 × 3 supercell of the investigated structure. Different colors of boron atoms represent similar atom types in the structure. The mentioned irregular hollow decagon is shown in Figure 2. For the rest of the article, we will use the α, β, and γ naming for the red, blue, and green represented boron atom types, respectively.

Figure 3 shows the band structure and the density of states (DOS) results of the unstrained Ogee-Borophene structure. Like the other allotropes and polymorphs of borophene, Ogee-Borophene also shows metallic behavior. There are several Dirac cone formations observed in the band structure however, the ones near the Fermi level in between region and region at 0.3 or −0.02 eV can be attributed as interesting, respectively. At this point, state density becomes zero which can be further changed by applying strain or using a substrate. However, these further possible manipulations are out of the scope of this study and can be the subject of another study.

To better understand the electronic structure of the investigated structure, projected DOS (PDOS) behavior is investigated in Figure 4. Figure 4(a) shows that most of the states show p-orbital behavior. In Figure 4(b) the percentage of p-orbital behavior of the states is shown. Near the Fermi level, p-orbital behavior shows a maximum, where we can conclude a possible conductivity of this metal is near a complete p-orbital related. Near the Dirac cone, the graph is not continuous because of zero states at that energy level. To understand better the source of the states near the Fermi level, we investigated the PDOS behaviors of each atom type, which are named α, β, and γ in Figures 4(c) and 4(d). Near the Fermi level, s-orbital states of β-boron and p-orbital states of α-boron are dominated. Because the p-orbital is dominant at the Fermi level, it can be suggested that the p-orbital states of α-boron are deciding the nature of the conduction in the investigated structure.

In Figure 5, PDOS of the px, py, and pz orbitals are shown for α, β, and γ types of boron atoms which are mentioned before. Because of the structural anisotropy of the Ogee-Borophene, it is not possible to any similarities between px and py orbitals. For α and β types, valance electrons are highly py orbital dominated. And for β and γ types, conduction electrons are highly pz orbital dominated. Although the conduction electrons in the α-type boron atom are generally in the px orbital characteristic, it is possible to observe that the above-mentioned electrons, which are near the Fermi level and have a dominant role in conduction, have pz orbital behavior. Also from the figure, it is possible to understand the origin of Dirac cone formations observed in the band structure near the Fermi level in between region and region. Dirac cone formation between region is mostly dominated by pz orbital of α type of boron, however Dirac cone formation between region is dominated by pz orbital of all boron types and px orbital of α type of boron. Therefore, it can be argued that the α type of boron has an important role in explaining the Dirac formations, which are important because of having unique electronic properties, near the Fermi level. It may be possible to argue that α type of boron atoms is also important because the structural anisotropy of the Ogee-Borophene is caused by them. By removing the α type of boron atoms from the structure, the hollows of the system will become isotropic.

In Figure 6, isosurfaces of electron density distributions with blue and yellow isosurfaces are shown for densities of 0.1 and 0.15 , respectively. Anisotropy in a bond between two boron atoms can be observed in yellow isosurfaces. Furthermore, it can be suggested that the variations in the bond lengths within the structure, gave a potential for different interaction strengths along the decagon-shaped hollows compared to the bridges between the gaps.

Figure 7 shows the calculated phonon dispersion curves and the phonon DOS for the Ogee-Borophene structure. The calculated phonon dispersion curves contain many phonon branches with imaginary frequencies across the entire k-space. Imaginary frequencies are mostly regarded as the inherent structural instability of the investigated structure. However, the presence of imaginary frequencies in the phonon spectrum can provide valuable insight into the nature and properties of a material [43]. The investigation of Ogee-Borophene in this study suggests that the “free-standing” form of this polymorph is not thermodynamically stable. The term “free-standing” refers to the isolated Ogee-Borophene structure without any substrate or support. Decagonal-shaped hollows and anisotropy of the structure may contribute to its potential instability in free-standing form. However, it is worth considering the possibility of achieving stability by investigating the Ogee-Borophene on known metallic or semiconductor surfaces, which we may call a substrate. By introducing a substrate, that can provide structural support and influence the interaction between Ogee-Borophene and its surroundings, it may be possible to stabilize the structure. Further DFT calculations and molecular dynamics simulations can help elucidate the energetic, adhesion, and electronic coupling between Ogee-Borophene and the possible substrates. These studies can guide the selection of optimal substrate materials and conditions for stabilizing Ogee-Borophene.

4. Conclusion

This study presents first-principle DFT calculations of a new hollow irregular decagon-shaped vacancy-mediated polymorph of borophene, named Ogee-Borophene. The irregular shape of the hollow decagon gives rise to different electronic properties for different atoms near the hollow decagons. Based on the DFT calculations, three types of boron atoms with distinct electronic properties, denoted as α, β, and γ, are found in this hypothetical boron polymorph. The electronic band structure shows metallic behavior with a Dirac cone far above the Fermi level. DOS and PDOS of the investigated structure show that conductivity is mostly dominated by the pz orbital of the α-type boron atom. Our study demonstrates the different electronic properties of different boron atoms provide tuning or selective interaction functionality to this hypothetical 2D metallic polymorph. The unique electronic properties of this material, arising from the anisotropy of the structure and the presence of boron atoms with different electronic properties, suggest that Ogee-Borophene has the potential to be used in a variety of applications, including transistors, selective sensor applications without the need for doping.

The study also found that Ogee-Borophene cannot be found in a free-standing stable form. However, the possibility of obtaining stability by investigating the interaction with substrates or building heterostructured 2D materials provides a promising direction for further research.

Overall, the present study makes a significant contribution to the scientific community by providing new insights and approaches into the fundamental properties of borophene and demonstrating the potential of this material for a variety of applications.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The numerical calculations reported in this paper were partially performed at TUBITAK ULAKBIM, High Performance and Grid Computing Center (TRUBA resources), and the National Center for High-Performance Computing of Turkey (UHeM) under grant number 1012802022.