Electronic Structure of Boron Flat Holeless Sheet

The electronic band structure, namely energy band surfaces and densities-of-states (DoS), of a hypothetical flat and ideally perfect, i.e., without any type of holes, boron sheet with a triangular network is calculated within a quasi-classical approach. It is shown to have metallic properties as is expected for most of the possible structural modifications of boron sheets. The Fermi curve of the boron flat sheet is found to be consisted of 6 parts of 3 closed curves, which can be approximated by ellipses representing the quadric energy-dispersion of the conduction electrons. The effective mass of electrons at the Fermi level in a boron flat sheet is found to be too small compared with the free electron mass m0 and to be highly anisotropic. Its values distinctly differ in directions Γ–K and Γ–M: mΓ–K/m0 ≈ 0.480 and mΓ–M/m0 ≈ 0.052, respectively. The low effective mass of conduction electrons, mσ/m0 ≈ 0.094, indicates their high mobility and, hence, high conductivity of the boron sheet. The effects of buckling/puckering and the presence of hexagonal or other type of holes expected in real boron sheets can be considered as perturbations of the obtained electronic structure and theoretically taken into account as effects of higher order.


Why should Boron Sheets be Formed?
There are different reasons in favor of the formation of stable 2-D all-boron structures.They can be divided into several groups.Let's consider them separately.

3-D All-Boron Structures
Boron, the fifth element of the Periodic Table, is located at the intersection of semiconductors and metals.Due to a small covalent radius (only 0.84 Å) and number (only 3) of valence electrons, boron does not form simple three-dimensional structures but crystals with icosahedral clusters with many atoms in the unit cell.At least three all-boron allotropes are known-α-and β-rhombohedral and high-pressure γ-orthorhombic phases-for the experimental phase diagram of boron see Reference [18].In addition, αand β-tetragonal and a number of other boron structures, probably stabilized by the presence of impurities/defects, were reported.Theoretical studies of the five boron crystal structures (α, dhcp, sc, fcc, and bcc) were carried out using the LAPW (linearized plane wave) method in Reference [19].The current state of research on the phase diagram of boron from a theoretical point of view is given in Reference [20].It should be noted that, in the last decade, several new structures of boron allotropes were discovered and some have been disproved.Currently, even the number of allotropes of boron is uncertain.The reason for this is that there are many such structures, all of them complex, and some of them are minimally different from others.A pseudo-cubic tetragonal boron recently discovered under high-pressure and high-temperature conditions may also be another form of boron allotropes; however, its structure, studied in Reference [21] using a DFT (density functional theory) calculation, is abnormal compared to other allotropes of boron in many ways.
The almost regular icosahedron B 12 with B-atoms at the vertices (Figure 1) serves as the main structural motif not only of boron allotropes but also of all known boron-rich compounds.In the boron icosahedron, each atom is surrounded by 5 neighboring atoms and, as usual, with one more atom from the rest of the crystal.For this reason, the average coordination number of a boron-rich lattice ranges from 5 to 5+1=6.
Condens.Matter 2019, 4, x 2 of 22 many atoms in the unit cell.At least three all-boron allotropes are known-α-and β-rhombohedral and high-pressure γ-orthorhombic phases-for the experimental phase diagram of boron see Reference [18].In addition, α-and β-tetragonal and a number of other boron structures, probably stabilized by the presence of impurities/defects, were reported.Theoretical studies of the five boron crystal structures (α, dhcp, sc, fcc, and bcc) were carried out using the LAPW (linearized plane wave) method in Reference [19].The current state of research on the phase diagram of boron from a theoretical point of view is given in Reference [20].It should be noted that, in the last decade, several new structures of boron allotropes were discovered and some have been disproved.Currently, even the number of allotropes of boron is uncertain.The reason for this is that there are many such structures, all of them complex, and some of them are minimally different from others.A pseudo-cubic tetragonal boron recently discovered under high-pressure and high-temperature conditions may also be another form of boron allotropes; however, its structure, studied in Reference [21] using a DFT (density functional theory) calculation, is abnormal compared to other allotropes of boron in many ways.
The almost regular icosahedron B12 with B-atoms at the vertices (Figure 1) serves as the main structural motif not only of boron allotropes but also of all known boron-rich compounds.In the boron icosahedron, each atom is surrounded by 5 neighboring atoms and, as usual, with one more atom from the rest of the crystal.For this reason, the average coordination number of a boron-rich lattice ranges from 5 to 5+1=6.However, an isolated regular boron icosahedron is an electron-deficient structure-the total number of valence electrons of 12 boron atoms is not sufficient to fill all the covalent bonding orbitals corresponding to such a cage-molecule.Thus, if it were a stable structure, then intra-icosahedral bonds would be only partially covalent but also to some extent metallic.As for boron icosahedra constituting real crystals, it was clearly demonstrated, for example, for β-rhombohedral boron [22][23][24][25][26][27], that they are stabilized by the presence of point structural defects-vacancies and interstitials, in other words, both partially filled regular or irregular boron sites-at very high concentrations.For example, in the case of β-rhombohedral boron, the total effect of such a stabilization is to increase the average number of boron atoms inside the unit cell from the ideal value 105 (Figure 2) to 106.7 [28], which leads to the saturation of the electron-deficient orbitals and a 5-or 6-coordination number for the majority of constituent boron atoms.However, an isolated regular boron icosahedron is an electron-deficient structure-the total number of valence electrons of 12 boron atoms is not sufficient to fill all the covalent bonding orbitals corresponding to such a cage-molecule.Thus, if it were a stable structure, then intra-icosahedral bonds would be only partially covalent but also to some extent metallic.As for boron icosahedra constituting real crystals, it was clearly demonstrated, for example, for β-rhombohedral boron [22][23][24][25][26][27], that they are stabilized by the presence of point structural defects-vacancies and interstitials, in other words, both partially filled regular or irregular boron sites-at very high concentrations.For example, in the case of β-rhombohedral boron, the total effect of such a stabilization is to increase the average number of boron atoms inside the unit cell from the ideal value 105 (Figure 2) to 106.7 [28], which leads to the saturation of the electron-deficient orbitals and a 5-or 6-coordination number for the majority of constituent boron atoms.many atoms in the unit cell.At least three all-boron allotropes are known-α-and β-rhombohedral and high-pressure γ-orthorhombic phases-for the experimental phase diagram of boron see Reference [18].In addition, α-and β-tetragonal and a number of other boron structures, probably stabilized by the presence of impurities/defects, were reported.Theoretical studies of the five boron crystal structures (α, dhcp, sc, fcc, and bcc) were carried out using the LAPW (linearized plane wave) method in Reference [19].The current state of research on the phase diagram of boron from a theoretical point of view is given in Reference [20].It should be noted that, in the last decade, several new structures of boron allotropes were discovered and some have been disproved.Currently, even the number of allotropes of boron is uncertain.The reason for this is that there are many such structures, all of them complex, and some of them are minimally different from others.A pseudo-cubic tetragonal boron recently discovered under high-pressure and high-temperature conditions may also be another form of boron allotropes; however, its structure, studied in Reference [21] using a DFT (density functional theory) calculation, is abnormal compared to other allotropes of boron in many ways.
The almost regular icosahedron B12 with B-atoms at the vertices (Figure 1) serves as the main structural motif not only of boron allotropes but also of all known boron-rich compounds.In the boron icosahedron, each atom is surrounded by 5 neighboring atoms and, as usual, with one more atom from the rest of the crystal.For this reason, the average coordination number of a boron-rich lattice ranges from 5 to 5+1=6.However, an isolated regular boron icosahedron is an electron-deficient structure-the total number of valence electrons of 12 boron atoms is not sufficient to fill all the covalent bonding orbitals corresponding to such a cage-molecule.Thus, if it were a stable structure, then intra-icosahedral bonds would be only partially covalent but also to some extent metallic.As for boron icosahedra constituting real crystals, it was clearly demonstrated, for example, for β-rhombohedral boron [22][23][24][25][26][27], that they are stabilized by the presence of point structural defects-vacancies and interstitials, in other words, both partially filled regular or irregular boron sites-at very high concentrations.For example, in the case of β-rhombohedral boron, the total effect of such a stabilization is to increase the average number of boron atoms inside the unit cell from the ideal value 105 (Figure 2) to 106.7 [28], which leads to the saturation of the electron-deficient orbitals and a 5-or 6-coordination number for the majority of constituent boron atoms.Thus, all-boron 5-and 6-coordinated regular 3-D lattices cannot exist, but one can naturally imagine 2-D flat or buckled/puckered structures with a triangular arrangement of atoms with and without periodically spaced hexagonal (rarely quadric, pentagonal, or heptagonal) holes.Obviously, most of them are expected to be (semi)metallic.At the moment, a number of different atomic geometries for quasi-planar boron sheets are theoretically proposed [4][5][6][7].
Using the ab initio evolutionary structure prediction approach, a novel reconstruction of the α-boron (111) surface with the lowest energy was discovered [52].In this reconstruction, all single interstitial boron atoms bridge neighboring icosahedra by polar covalent bonds, and this satisfies the electron counting rule, leading to the reconstruction-induced semiconductor-metal transition.The new stable boron sheet, called H-borophene, proposed in Reference [53] and constructed by tiling 7-membered rings side by side, should be especially noted.
As for the irregularly distributed holes, they have to be considered as defects.The research [54] is focused on the formation of local vacancy defects and pinholes in a 2-D boron structure-the so-called γ 3 -type boron monolayer.

Boron Quasi-Planar Clusters
Indirectly, the reality of boron sheets can be proved by the presence of various quasi-planar boron clusters, i.e., finite fragments of sheets, in gaseous state and also boron nanotubes, which are the fragments of boron sheets wrapped into cylinders (see for example, the review from Reference [55] and the references therein).Experimental and theoretical evidences that small boron clusters prefer planar structures were reported in Reference [8].
In addition, recently, a highly stable quasi-planar boron cluster B 36 of hexagonal shape with a central hexagonal hole [9], which is viewed as a potential basis for an extended 2-D boron sheet, and boron fullerene B 40 [10], which can be imagined as the fragment of a boron sheet wrapped into the sphere, were discovered experimentally.Photoelectron spectroscopy in combination with ab initio calculations have been carried out to probe the structure and chemical bonding of the B 27 − cluster [56].
A comparison between the experimental spectrum and the theoretical results reveals a 2-D global minimum with a triangular lattice containing a tetragonal defect and two low-lying 2-D isomers, each with a hexagonal vacancy.

Liquid Boron Structure
There are also evidences [57][58][59] that liquid boron does not consist of icosahedra but mainly of quasi-planar clusters.Ab initio MD (molecular-dynamics) simulations of the liquid boron structure yields that at short length scales, B 12 icosahedra, a main structural motif of boron crystals and boron-rich solid compounds, are destroyed upon melting.Although atoms form an open packing, they maintain the 6-coordination.
According to measurements of the structure factor and the pair distribution function, the melting process is associated with relatively small changes in both the volume and the short-range order of the system.Results of a comprehensive study of liquid boron with X-ray measurements of the atomic structure and dynamics coupled with ab initio MD simulations also show that there is no evidence of survival of the icosahedral arrangements into the liquid, but many atoms appear to adopt a geometry corresponding to the quasi-planar pyramids.

Growing of Boron Sheets
Currently, some of the 2-D materials beyond graphene also are used [60].But for non-layer structured 3-D materials such as boron, it is a real challenge to fabricate the corresponding 2-D nanosheets due to the absence of the driving force of anisotropic growth.There are rare examples of some 2-D metal nanosheets; see for example, the recent report [61] on single-crystalline Rh nanosheets with a 3-5 atomic layers thickness.Boron sheets are expected to be metallic as well.Thus, this should increase the chances of their actual formation.
In this regard, we have to mention the recent report [62] in which large-scale single-crystalline ultrathin boron nanosheets have been fabricated via the thermal decomposition of diborane.
It is obvious that an infinite boron sheet does not exist in nature and that its finite pieces are not stable compared to bulk and/or nanotubular structures of boron.To grow boron sheets, one needs a substrate which binds boron atoms strongly to avoid bulk phases while, at the same time, provides sufficient mobility of boron atoms on the substrate.Possible candidates for substrates are surfaces of (close-packed transition) metals.The feasibility of different synthetic methods for 2-D boron sheets was assessed [47,[63][64][65][66] using ab initio calculations, i.e., "synthesis in theory" approach.A large-scale boron monolayer has been predicted with mixed hexagonal-triangular geometry obtained via either depositing boron atoms directly on the surface or soft landing of small planar B-clusters.
Recently, a series of planar boron allotropes with honeycomb topology has been proposed [67].Although the free-standing honeycomb B allotropes are higher in energy than α-sheets, these calculations show that a metal substrate can greatly stabilize these new allotropes.
The atomically thin, crystalline 2-D boron sheets, i.e., borophene, were actually synthesized [16] on silver surfaces under ultrahigh-vacuum conditions (Figure 3).An atomic-scale characterization, supported by theoretical calculations, revealed structures reminiscent of fused boron clusters with multiple scales of anisotropic, out-of-plane buckling.Unlike bulk boron allotropes, borophene shows metallic characteristics that are consistent with predictions of a highly anisotropic 2-D metal.In this regard, we have to mention the recent report [62] in which large-scale single-crystalline ultrathin boron nanosheets have been fabricated via the thermal decomposition of diborane.
It is obvious that an infinite boron sheet does not exist in nature and that its finite pieces are not stable compared to bulk and/or nanotubular structures of boron.To grow boron sheets, one needs a substrate which binds boron atoms strongly to avoid bulk phases while, at the same time, provides sufficient mobility of boron atoms on the substrate.Possible candidates for substrates are surfaces of (close-packed transition) metals.The feasibility of different synthetic methods for 2-D boron sheets was assessed [47,[63][64][65][66] using ab initio calculations, i.e., "synthesis in theory" approach.A large-scale boron monolayer has been predicted with mixed hexagonal-triangular geometry obtained via either depositing boron atoms directly on the surface or soft landing of small planar B-clusters.
Recently, a series of planar boron allotropes with honeycomb topology has been proposed [67].Although the free-standing honeycomb B allotropes are higher in energy than α-sheets, these calculations show that a metal substrate can greatly stabilize these new allotropes.
The atomically thin, crystalline 2-D boron sheets, i.e., borophene, were actually synthesized [16] on silver surfaces under ultrahigh-vacuum conditions (Figure 3).An atomic-scale characterization, supported by theoretical calculations, revealed structures reminiscent of fused boron clusters with multiple scales of anisotropic, out-of-plane buckling.Unlike bulk boron allotropes, borophene shows metallic characteristics that are consistent with predictions of a highly anisotropic 2-D metal.The experimental work in Reference [17] shows that 2-D boron sheets can be grown epitaxially on a Ag(111) substrate.Two types of boron sheets, β12 and χ3, both exhibiting a triangular lattice but with different arrangements of periodic holes, were observed by scanning tunneling microscopy.DFT simulations indicate that both sheets are planar without obvious vertical undulations.
According to the ab initio calculations [68], the periodic nanoscale 1-D undulations can be preferred in borophenes on concertedly reconstructed Ag(111).This "wavy" configuration is more stable than its planar form on flat Ag(111) due to an anisotropic high bending flexibility of borophene.An atomic-scale ultrahigh vacuum scanning tunneling microscopy characterization of a borophene grown on Ag(111) reveals such undulations, which agree with the theory.Although the lattice is coherent within a borophene island, the undulations nucleated from different sides of the island form a distinctive domain boundary when they are laterally misaligned.
Recently, borophene synthesis monitored in situ by low-energy electron microscopy, diffraction, and scanning tunneling microscopy and modeled using ab initio theories has been reported in Reference [69].By resolving the crystal structure and phase diagram of borophene on Ag(111), the domains are found to remain nanoscale for all growth conditions.However, by growing borophene on Cu(111) surfaces, large single-crystal domains (up to 100 μm) are obtained.The crystal structure is a novel triangular network with a concentration of hexagonal vacancies of η= 1/5.These experimental data together with ab initio calculations indicate a charge-transfer coupling to the substrate without significant covalent bonding.
Boron on a Pb(110) surface was simulated [70] by using ab initio evolutionary methodology and found that 2-D Pmmn structures can be formed because of a good lattice matching.By increasing the The experimental work in Reference [17] shows that 2-D boron sheets can be grown epitaxially on a Ag(111) substrate.Two types of boron sheets, β 12 and χ 3 , both exhibiting a triangular lattice but with different arrangements of periodic holes, were observed by scanning tunneling microscopy.DFT simulations indicate that both sheets are planar without obvious vertical undulations.
According to the ab initio calculations [68], the periodic nanoscale 1-D undulations can be preferred in borophenes on concertedly reconstructed Ag(111).This "wavy" configuration is more stable than its planar form on flat Ag(111) due to an anisotropic high bending flexibility of borophene.An atomic-scale ultrahigh vacuum scanning tunneling microscopy characterization of a borophene grown on Ag(111) reveals such undulations, which agree with the theory.Although the lattice is coherent within a borophene island, the undulations nucleated from different sides of the island form a distinctive domain boundary when they are laterally misaligned.
Recently, borophene synthesis monitored in situ by low-energy electron microscopy, diffraction, and scanning tunneling microscopy and modeled using ab initio theories has been reported in Reference [69].By resolving the crystal structure and phase diagram of borophene on Ag(111), the domains are found to remain nanoscale for all growth conditions.However, by growing borophene on Cu(111) surfaces, large single-crystal domains (up to 100 µm) are obtained.The crystal structure is a novel triangular network with a concentration of hexagonal vacancies of η = 1/5.These experimental data together with ab initio calculations indicate a charge-transfer coupling to the substrate without significant covalent bonding.
Boron on a Pb(110) surface was simulated [70] by using ab initio evolutionary methodology and found that 2-D Pmmn structures can be formed because of a good lattice matching.By increasing the thickness of 2-D boron, the three-bonded graphene-like P2 1 /a boron was revealed to possess lower molar energy, indicating the more stable 2-D boron.
The influence of the excess negative charge on the stability of borophenes-2-D boron crystals-was examined in Reference [71] using an analysis of the decomposition of the binding energy of a given boron layer into contributions coming from boron atoms that have different coordination numbers to understand how the local neighborhood of an atom influences the overall stability of the monolayer structure.The decomposition is done for the α-sheet related family of structures.It was found a preference for 2-D boron crystals with very small or very high charges per atom.Structures with intermediate charges are not energetically favorable.It has been also found a clear preference in terms of binding energy for the experimentally observed γ-sheet and δ-sheet structures that is almost independent on the considered excess of negative charge of the structures.
Two-dimensional boron monolayers have been extensively investigated using ab initio calculations [72].A series of boron bilayer sheets with pillars and hexagonal holes have been constructed.Many of them have a lower formation energy than an α-sheet boron monolayer.However, the distribution and arrangement of hexagonal holes can cause a negligible effect on the stability of these structures.
Recently, an ab initio study [73] of the effect of electron doping on the bonding character and stability of borophene for the neutral system has revealed previously unknown stable 2-D structures: ε-B and ω-B.The chemical bonding characteristic in this and other boron structures is found to be strongly affected by an extra charge.Beyond a critical degree of electron doping, the most stable allotrope changes from ε-B to a buckled honeycomb structure.Additional electron doping, mimicking a transformation of boron to carbon, causes a gradual decrease in the degree of buckling of the honeycomb lattice.

Applications
In general, the formation of a boron sheet would have wide applications in techniques because the boronizing of metal surfaces is known as an effective method of formation of protective coatings [74].In particular, quasi-planar bare boron surfaces can serve as lightweight protective armor.
Boron sheets are expected to be very good conductors with potential applications in nanoelectronics, e.g., in high-temperature nanodevices.Boron sheets could have potential as metallic interconnects and wiring in electronic devices and IC (integrated circuits) [41].
A theoretical investigation [32] of both the molecular physisorption and dissociative atomic chemisorption of hydrogen by boron sheets predicts physisorption as the leading mechanism at moderate temperatures and pressures.Further calculations performed on hydrogen-storage properties showed that the decoration of pristine sheets with the right metal elements provide additional absorption sites for hydrogen [47].Thus, boron sheets can serve for good nanoreservoirs of fuel hydrogen used in green-energy production.
Due to the high neutron-capture cross section of 10 B nuclei, solid-state boron allotropies, as well as boron-rich compounds and composites, are good candidates to be used as neutron-protectors.Boron sheets will be especially useful as an absorbing component in composite neutron shields [75].Materials with the high bulk concentration of B-atoms usually are nonmetals and, therefore, not suitable for electromagnetic shielding purposes.However, frequently, the simultaneous protection against both the neutron irradiation and electromagnetic waves is needed, in particular, because neutron absorption by 10 B nuclei is accompanied by a gamma-radiation.For this reason, in the boron-containing nanocomposites designed for neutron-protection, it is necessary to introduce some foreign components with metallic conductivity.Utilizing of the metallic boron sheet as a component may resolve this problem [76].
Recently, the mechanical properties of 2-D boron-borophene-have been studied by ab initio calculations [77].The borophene with a 1/6 concentration of hollow hexagons is shown to have the Föppl-von Kármán number per unit area over twofold higher than graphene's value.The record high flexibility combined with excellent elasticity in boron sheets can be utilized for designing advanced composites and flexible devices.The transfer of undulated borophene onto an elastomeric substrate would allow for high levels of stretchability and compressibility with potential applications for emerging stretchable and foldable devices [68].
The boron sheets are quite inert to oxidization and interact only weakly with their substrate.For this reason, they may find applications in electronic devices in the future [17].
In large-scale single-crystalline ultrathin boron nanosheets fabricated [62] via the thermal decomposition of diborane, the strong combination performances of low turn-on field-of-field emissions, favorable electron transport characteristics, high sensitivity, and fast response time to illumination reveal that the nanosheets have high potential as applications in field-emitters, interconnects, IC, and optoelectronic devices.
Some other applications of borophene are described in recent reviews [78,79].

Available Electron Structure Calculations
Because boron sheets are of great academic and practical interests, their electronic structure is studied intensively.Most of them are found to be metallic.
Let us note that there are some indirect evidences for metallic conduction in boron sheets.The absence of icosahedra in liquid boron affects its properties including electrical conductivity [57,59], and it behaves like a metal.
The very stable quasi-planar clusters of boron B n for n up to 46, considered to be fragments of bare boron quasi-planar surfaces, have to possess a singly occupied bonding orbital [29].Assuming that conduction band of the infinite surface is generated from the HOMO (highest-occupied-molecular-orbital) of a finite fragment, it means the partial filling of the conduction band, i.e., the metallic mechanism of conductance.
Diamond-like, metallic boron crystal structures were predicted in Reference [80] employing so-called decoration schemes of calculations, in which the normal and hexagonal diamond-like frameworks are decorated with extra atoms across the basal plane.They should have an overly high DoS near the Fermi level.This result may provide a plausible explanation for not only the anomalous superconductivity of boron under high pressure but also the nonmetal-metal transition in boron structures.
In Reference [45], it was presented the results of a theoretical study of the phase diagram of elemental boron showing that, at high pressures, boron crystallizes in quasi-layered bulk phases characterized by in-plane multicenter bonds and out-of-plane bonds.All these structures are metallic.
Band structures of a series of planar boron allotropes with honeycomb topologies recently proposed in Reference [67], exhibit Dirac cones at the K-point, the same as in graphene.In particular, the Dirac point of honeycomb boron sheet locates precisely on the Fermi level, rendering it as a topologically equivalent material to graphene.Its Fermi velocity is of 6•10 5 m/s, close to that of graphene.However, in H-borophene [53] constructed by tiling 7-membered rings side by side, a Dirac point appeared at about 0.33 eV below the Fermi level.
According to some theoretical results [36,37,47,48], boron sheets can be not only metal but in some cases also an almost zero band-gap semiconductor depending on its atomistic configuration.Probably, the semiconducting character is related to the nonzero thickness of buckled/puckered 2-D boron sheets or double-layered structures.
Some of borophenes can be magnetic.Based on a tight-binding model of 8-Pmmn borophene developed in Reference [81], it is confirmed that the crystal hosts massless Dirac fermions and the Dirac points are protected by symmetry.Strain is introduced into the model, and it is shown to induce a pseudomagnetic field vector potential and a scalar potential.The 2-D antiferromagnetic boron, designated as M-boron, has been predicted [82] using an ab initio evolutionary methodology.M-boron is entirely composed of B 20 clusters in a hexagonal arrangement.Most strikingly, the highest valence band of M-boron is isolated, strongly localized, and quite flat, which induces spin polarization on either cap of the B 20 cluster.This flat band originates from the unpaired electrons of the capping atoms and is responsible for magnetism.M-boron is thermodynamically metastable.
Boron sheets grown on metal surfaces are predicted [63] to be strongly doped with electrons from the substrate, showing that a boron sheet is an electron-deficient material.As mentioned, by simulating [70] boron on Pb(110) surface using ab initio evolutionary methodology, it was found that the 2-D Dirac Pmmn boron can be formed.Unexpectedly, by increasing the thickness of 2-D boron, the three-bonded graphene-like structure P2 1 /a was revealed to possess double anisotropic Dirac cones.It is the most stable 2-D boron with particular Dirac cones.The puckered structure of P2 1 /a boron results in the peculiar Dirac cones.
This present work aims to provide more detailed calculations on the electronic structure of boron sheet including not only DoS but also band structure, electron effective mass, Fermi curve, etc.

Theoretical Approach
We use an original theoretical method of the quasi-classical type [83] based on the proof that the electronic system of any substance is a quasi-classical system; that is, its exact and quasi-classical energy spectra are close to each other.
As for the determination of the materials' electron structure, the quasi-classical method is reduced to the LCAO (linear combination of atomic orbitals) method with a basis set of quasi-classical atomic orbitals.Within the initial quasi-classical approximation, the solution of the corresponding mathematical problem consists of two main stages: 1.
The construction of matrix elements for secular equation, which, within the initial quasi-classical approximation, reduces to a geometric task of determining the volume of the intersection of three spheres [89], and 2.
The solving of the secular equation, which determines the crystalline electronic energy spectrum [90].
This method has been successfully applied to electronic structure calculations performed for various modifications of boron nitride, BN, one of the most important boron compounds [91][92][93][94], as well as metal-doped β-rhombohedral boron [95].
The maximal relative error of a quasi-classical calculation itself, i.e., without errors arisen from input data, is estimated as approx.4%.
As for the input data, the quasi-classical method of band structure calculations requires them to be in the form of quasi-classical parameters of constituent atoms: the inner and outer radii of the classical turning points for electron states in atoms, the radii of layers of the quasi-classical averaging of potential in atoms, and the averaged values of the potential within corresponding radial layers of atoms.These quantities for an isolated boron atom (as well as for other atoms) in the ground state were pre-calculated in Reference [96] on the basis of ab initio theoretical, namely Hartree-Fock (HF), values of electron levels [97].Thus, in our case, the accuracy of the quasi-classical parameters is determined by that of the HF approach.
As is known, the electronic structure of any atomic system is influenced by its geometric structure and vice versa.Often, one starts with the question of how to find the most stable idealized atomic configuration.Despite this, here, we will directly begin with the electron band structure of a flat triangular boron sheet, neglecting the buckling/puckering effects and hexagonal holes (see references above), assuming that in real sheets (e.g., grown on metal surfaces), their buckled/puckered or vacant parts would not be arranged in a periodic manner, and thus, they should be regarded as perturbations which can be taken into account within a higher-order perturbation theory.
Multi-layered (buckled) boron sheets can be imagined by substituting metal Me atoms again with boron B atoms in the layered structure of a metal boride MeB 2 .Analyses of the isolated layer instead of multilayered structure also seems to be quite sufficient for the initial approximation because, in such structures, only intra-layer conductivity is metallic, while interlayer conductivity is nonmetallic due to larger interlayer bond lengths if compared with those in layers.
The 2-D unit cell of the perfectly flat boron sheet without hexagonal holes (Figure 4) is a rhomb with an acute angle of β = π/3, i.e., with a single lattice constant a (Figure 5).Let There is a number of different values for the lattice constant of a boron sheet suggested theoretically.For the self-consistency of calculations, we use a = 3.37 a.u. of length, i.e., 1.78 Å, which corresponds to the B-B pair interatomic potential in the same quasi-classical approximation [91].
Condens.Matter 2019, 4, x 8 of 22 instead of multilayered structure also seems to be quite sufficient for the initial approximation because, in such structures, only intra-layer conductivity is metallic, while interlayer conductivity is nonmetallic due to larger interlayer bond lengths if compared with those in layers.The 2-D unit cell of the perfectly flat boron sheet without hexagonal holes (Figure 4) is a rhomb with an acute angle of 3 / π β = , i.e., with a single lattice constant a (Figure 5).Into the basis set of a simple LCAO formalism, it has been included core (1s), fully (2s), and partially filled (2p) valence and empty excited (2p) atomic orbitals.
Experimentally, there are 10 detected different existing states in the boron atom.To minimize the calculation errors related to the approximation of the crystalline potential by the superposition of atomic potentials, we choose orbitals with the same symmetry as the partially filled valence orbital, i.e., 2p, with the closest energy level and, consequently, with the closest classical turning point radii of electrons.
Taking into account the degeneracy of atomic energy levels by magnetic and spin quantum numbers of 2, 2, 6, and 6, respectively, we can state that this set of 4 orientation-averaged orbitals replaces 16 angularly dependent atomic orbitals.
The secular equation takes the form instead of multilayered structure also seems to be quite sufficient for the initial approximation because, in such structures, only intra-layer conductivity is metallic, while interlayer conductivity is nonmetallic due to larger interlayer bond lengths if compared with those in layers.The 2-D unit cell of the perfectly flat boron sheet without hexagonal holes (Figure 4) is a rhomb with an acute angle of 3 / π β = , i.e., with a single lattice constant a (Figure 5).  .Then, the radius-vectors of the lattice sites are    Into the basis set of a simple LCAO formalism, it has been included core (1s), fully (2s), and partially filled (2p) valence and empty excited (2p) atomic orbitals.
Experimentally, there are 10 detected different existing states in the boron atom.To minimize the calculation errors related to the approximation of the crystalline potential by the superposition of atomic potentials, we choose orbitals with the same symmetry as the partially filled valence orbital, i.e., 2p, with the closest energy level and, consequently, with the closest classical turning point radii of electrons.
Taking into account the degeneracy of atomic energy levels by magnetic and spin quantum numbers of 2, 2, 6, and 6, respectively, we can state that this set of 4 orientation-averaged orbitals replaces 16 angularly dependent atomic orbitals.
The secular equation takes the form instead of multilayered structure also seems to be quite sufficient for the initial approximation because, in such structures, only intra-layer conductivity is metallic, while interlayer conductivity is nonmetallic due to larger interlayer bond lengths if compared with those in layers.The 2-D unit cell of the perfectly flat boron sheet without hexagonal holes (Figure 4) is a rhomb with an acute angle of 3 / π β = , i.e., with a single lattice constant a (Figure 5).  .Then, the radius-vectors of the lattice sites are    Into the basis set of a simple LCAO formalism, it has been included core (1s), fully (2s), and partially filled (2p) valence and empty excited (2p) atomic orbitals.
Experimentally, there are 10 detected different existing states in the boron atom.To minimize the calculation errors related to the approximation of the crystalline potential by the superposition of atomic potentials, we choose orbitals with the same symmetry as the partially filled valence orbital, i.e., 2p, with the closest energy level and, consequently, with the closest classical turning point radii of electrons.
Taking into account the degeneracy of atomic energy levels by magnetic and spin quantum numbers of 2, 2, 6, and 6, respectively, we can state that this set of 4 orientation-averaged orbitals replaces 16 angularly dependent atomic orbitals.
The secular equation takes the form Into the basis set of a simple LCAO formalism, it has been included core (1s), fully (2s), and partially filled (2p) valence and empty excited (2p) atomic orbitals.
Experimentally, there are 10 detected different existing states in the boron atom.To minimize the calculation errors related to the approximation of the crystalline potential by the superposition of atomic potentials, we choose orbitals with the same symmetry as the partially filled valence orbital, i.e., 2p, with the closest energy level and, consequently, with the closest classical turning point radii of electrons.
Taking into account the degeneracy of atomic energy levels by magnetic and spin quantum numbers of 2, 2, 6, and 6, respectively, we can state that this set of 4 orientation-averaged orbitals replaces 16 angularly dependent atomic orbitals.
The secular equation takes the form where S(α Formally, these expressions contain infinite series.However, within the initial quasi-classical approximation, due to the finiteness of quasi-classical atomic radii, only a finite number of summands differs from zero.Thus, the series are terminated unambiguously.
The input data in a.u. in the form of quasi-classical parameters of boron atoms are shown in Tables 1 and 2. As it was mentioned above, the parameters of electron states fully or partially filled with electrons in the ground state were calculated on the basis of the theoretical, namely HF, values of electron levels, while for the excited state, we use the experimental value [98], which, however, is modulated by the multiplier of order of 1, 0.984151, leading to the coincidence between experimental and HF-theoretical first ionization potentials for an isolated boron atom: 0.304945 and 0.309856 a.u., respectively.Note that for the ground state, the 1s 2 2s 2 2p configuration is considered, not the 1s 2 2s2p 2 configuration, from which the ground state and first excited states of some boron-like ions arise [99].All the matrix elements and electron energies are calculated in points α 1 k 2 of the reciprocal space with parameters −1/2 ≤ α 1 , α 2 ≤ +1/2, i.e., within a rhombic unit cell of the reciprocal lattice (Figure 7).The first Brillouin zone for a boron flat sheet has a hexagonal shape.Of course, the areas of hexagonal and rhombic unit cells are equal (Figure 8).The unit cell is covered evenly by 1,002,001 points, at which the energy is found as a solution to the generalized eigenvalue problem.
The calculation has been performed in atomic units, a.u.Then, the results have been converted according to the relations: 1 a.u. of energy = 27.212eV and 1 a.u. of length = 0.52918 Å.
Based on the resulting data set, we have constructed the electron band surfaces, the distribution of DoS in the bands, and the Fermi curve, emphasizing that, instead of the Fermi surface, the characteristics of 3-D crystals, 2-D crystals are characterized by Fermi curves.

Results and Discussion
In the quasi-classically calculated electronic structure of the flat boron sheet, we resolve four bands of energy.We have to emphasize that for simplicity, the band surfaces below are shown over a rhombic (not hexagonal) domain.
The lowest energy band 1 E surface is found to be almost a plane placed at the level of    The unit cell is covered evenly by 1,002,001 points, at which the energy is found as a solution to the generalized eigenvalue problem.
The calculation has been performed in atomic units, a.u.Then, the results have been converted according to the relations: 1 a.u. of energy = 27.212eV and 1 a.u. of length = 0.52918 Å.
Based on the resulting data set, we have constructed the electron band surfaces, the distribution of DoS in the bands, and the Fermi curve, emphasizing that, instead of the Fermi surface, the characteristics of 3-D crystals, 2-D crystals are characterized by Fermi curves.

Results and Discussion
In the quasi-classically calculated electronic structure of the flat boron sheet, we resolve four bands of energy.We have to emphasize that for simplicity, the band surfaces below are shown over a rhombic (not hexagonal) domain.
The lowest energy band 1 E surface is found to be almost a plane placed at the level of   The unit cell is covered evenly by 1,002,001 points, at which the energy is found as a solution to the generalized eigenvalue problem.
The calculation has been performed in atomic units, a.u.Then, the results have been converted according to the relations: 1 a.u. of energy = 27.212eV and 1 a.u. of length = 0.52918 Å.
Based on the resulting data set, we have constructed the electron band surfaces, the distribution of DoS in the bands, and the Fermi curve, emphasizing that, instead of the Fermi surface, the characteristics of 3-D crystals, 2-D crystals are characterized by Fermi curves.

Results and Discussion
In the quasi-classically calculated electronic structure of the flat boron sheet, we resolve four bands of energy.We have to emphasize that for simplicity, the band surfaces below are shown over a rhombic (not hexagonal) domain.
The lowest energy band E 1 surface is found to be almost a plane placed at the level of E 1min = E 1max = −276.21eV.Thus, the chemical shift against the core 1s atomic level E 1s = −209.41eV equals to δE 1 = E 1s − E 1 = 66.80 eV.Dispute the shift of the B 1s atomic level, it retains an order of magnitude after transforming in an electronic band of the boron flat sheet.The lowest-lying band E 1 is fully filled with electrons.
The band E 2 is the highest fully filled band (Figure 9) with bottom at E 2min = −37.21eV and top at E 2max = −19.85eV, i.e., with a width of ∆E 2 = E 2max − E 2min = 17.36 eV.Note that, this range of energies is comparable in order of magnitude with a valence 2s atomic level of E 2s = −13.46eV.
The band E 3 (Figure 10) is partially filled, i.e., partially empty, with a bottom at E 3min = −23.08eV and top at E 3max = −17.16eV, i.e., with a width of ∆E 3 = E 3max − E 3min = 5.92 eV.Note that this range of energies is comparable in order of magnitude with a valence 2p atomic level E 2p = −8.43eV.
The band E 4 (Figure 11) is empty, with a bottom at E 4min = −17.65 eV and top at E 4max = −8.08 eV, i.e., with a width of ∆E 4 = E 4max − E 4min = 9.57 eV.Note that this range of energies is comparable in order of magnitude with the modulated value of the excited 2p level E 2p = −5.84eV.
Note that this range of energies is comparable in order of magnitude with a valence 2p atomic level     The band 4 E (Figure 11) is empty, with a bottom at eV 65 .17 The band 4 E (Figure 11) is empty, with a bottom at eV 65 .17  To compare easily our results with the literature data, in addition to the presentation of the band structure using the contour plots of the whole Brillouin zone in Figures 12 and 13, we plot the band energies (as well as their second derivatives and corresponding parabolic approximations) along the main lines of symmetry.
Our quasi-classical calculation of the crystalline band structure, like any other approach also utilizing HF parameters of constituting atoms, cannot determine the absolute values of energy parameters with a high accuracy.By this reason, the above mentioned value E Fermi cannot be used directly to determine the electron work function of the boron sheet.This goal can be achieved only after corrections are made to include the electron-correlation and to exclude the electron-self-interacting effects, which have to allow an accurate determination of the position for the vacuum level of energy E = 0.However, the shifting of the reference point at the energy axis does not affect the energy differences, which are credible as are determined with quite an acceptable accuracy.They are collected in Table 3.
conduction electrons at the Fermi level, are bounded inside the 2-D crystal.This result once more evidences the correctness of the calculations performed in this work.The total width of valence and conduction bands equal to eV 79 .17 .As expected, it is negligible if compared with that of a conduction band.
To compare easily our results with the literature data, in addition to the presentation of the band structure using the contour plots of the whole Brillouin zone in Figures 12 and 13, we plot the band energies (as well as their second derivatives and corresponding parabolic approximations) along the main lines of symmetry.    .As expected, it is negligible if compared with that of a conduction band.
To compare easily our results with the literature data, in addition to the presentation of the band structure using the contour plots of the whole Brillouin zone in Figures 12 and 13, we plot the band energies (as well as their second derivatives and corresponding parabolic approximations) along the main lines of symmetry.The Fermi curve of a boron flat sheet is found to be consisting of parts of a number of closed curves including concentric ones, the centre of which can be approximated by ellipse with long and short axes along the Γ-K and Γ-M directions, respectively (Figure 14).The Fermi curve of a boron flat sheet is found to be consisting of parts of a number of closed curves including concentric ones, the centre of which can be approximated by ellipse with long and short axes along the Γ-Κ and Γ-Μ directions, respectively (Figure 14).As it is known, for semiconductors, the effective mass conception referring to band zone curvatures is used to approximate the wave-vector dependence of electron energies near the band gap.As for metals, the Fermi surface curvature can be used to estimate the effective mass of conduction electrons and hence their mobility.
The ellipse representing a branch of intersection between the 3 E -band surface with the Fermi E -plane can be described by the equation  As it is known, for semiconductors, the effective mass conception referring to band zone curvatures is used to approximate the wave-vector dependence of electron energies near the band gap.As for metals, the Fermi surface curvature can be used to estimate the effective mass of conduction electrons and hence their mobility.
The ellipse representing a branch of intersection between the E 3 -band surface with the E Fermi -plane can be described by the equation where k Γ-K and k Γ-M are wave-number components along perpendicular axes Γ-K and Γ-M and F = E Fermi − E 3min = 3.66 eV is the Fermi energy.The effective masses m Γ-K and m Γ-M can be estimated from this equation if it is rewritten in the form of a normalized ellipse equation where k Γ-K 0 and k Γ-M 0 are half-axes in the directions Γ-K and Γ-M, respectively: Then, one can calculate the effective mass of the conduction electrons m σ , i.e., electrons placed at the Fermi level, from the relation The effective electron mass at the Fermi level reveals a significant anisotropy.For the central ellipse, the effective masses are m Γ-K /m 0 ≈ 0.480 and m Γ-M /m 0 ≈ 0.052, with m σ /m 0 ≈ 0.094, where m 0 is the free electron mass.
The Fermi curve of a boron flat sheet is found to consist of 6 parts of 3 ellipses representing the quadric energy-dispersion of the conduction electrons; see Figure 15., where 0 m is the free electron mass.
The Fermi curve of a boron flat sheet is found to consist of 6 parts of 3 ellipses representing the quadric energy-dispersion of the conduction electrons; see Figure 15., where 0 m is the free electron mass.
The Fermi curve of a boron flat sheet is found to consist of 6 parts of 3 ellipses representing the quadric energy-dispersion of the conduction electrons; see Figure 15.The overall shapes of DoSs obtained by us and previously by others, especially, in References [31,33,47] are rather similar but with some differences.It is understandable as these structures are buckled/puckered or flat variants of the same triangular lattice with or without hexagonal holes (Figure 17).The discrepancies may be attributed with the perturbations related to the mentioned structural changes and differences between the computing methods utilized, as well as the The overall shapes of DoSs obtained by us and previously by others, especially, in References [31,33,47] are rather similar but with some differences.It is understandable as these structures are buckled/puckered or flat variants of the same triangular lattice with or without hexagonal holes (Figure 17).The discrepancies may be attributed with the perturbations related to the mentioned structural changes and differences between the computing methods utilized, as well as the difference between projected onto in-plane or out-of-plane orbitals' densities-of-states (PDoSs) from the total DoS of the sheet.e overall shapes of DoSs obtained by us and previously by others, especially, in References 47] are rather similar but with some differences.It is understandable as these structures are d/puckered or flat variants of the same triangular lattice with or without hexagonal holes 17).The discrepancies may be attributed with the perturbations related to the mentioned ral changes and differences between the computing methods utilized, as well as the nce between projected onto in-plane or out-of-plane orbitals' densities-of-states (PDoSs) from al DoS of the sheet.The Fermi curve of the monolayer flat boron sheet approximated by parts of concentric c ellipse-like curves could be considered as a certain kind of topological analog of the Fermi su (Figure 18) in the form of a half-torus and distorted cylinder of magnesium diboride MgB2 [ which is believed to be a structural analog of the hypothetical multilayered boron sheet, where m Me atoms in a metal diboride MeB2 structure are replaced by B-atoms themselves.The low effective mass of conduction electrons at the Fermi level indicates a high mobil electrons and, hence, a high conductivity of the flat boron sheet.

Conclusions
In summary, we can conclude that the electronic band structure of a boron flat triangular has been calculated within a quasi-classical approach for the quasi-classical structural param (B-B bond length) of   [31], (b) from Reference [33], and (c) modified from Reference [37].
The Fermi curve of the monolayer flat boron sheet approximated by parts of concentric closed ellipse-like curves could be considered as a certain kind of topological analog of the Fermi surface (Figure 18) in the form of a half-torus and distorted cylinder of magnesium diboride MgB 2 [100], which is believed to be a structural analog of the hypothetical multilayered boron sheet, where metal Me atoms in a metal diboride MeB 2 structure are replaced by B-atoms themselves.The Fermi curve of the monolayer flat boron sheet approximated by parts of concentric closed ellipse-like curves could be considered as a certain kind of topological analog of the Fermi surface (Figure 18) in the form of a half-torus and distorted cylinder of magnesium diboride MgB2 [100], which is believed to be a structural analog of the hypothetical multilayered boron sheet, where metal Me atoms in a metal diboride MeB2 structure are replaced by B-atoms themselves.The low effective mass of conduction electrons at the Fermi level indicates a high mobility of electrons and, hence, a high conductivity of the flat boron sheet.

Conclusions
In summary, we can conclude that the electronic band structure of a boron flat triangular sheet has been calculated within a quasi-classical approach for the quasi-classical structural parameter (B-B bond length) of  The low effective mass of conduction electrons at the Fermi level indicates a high mobility of electrons and, hence, a high conductivity of the flat boron sheet.

Conclusions
In summary, we can conclude that the electronic band structure of a boron flat triangular sheet has been calculated within a quasi-classical approach for the quasi-classical structural parameter (B-B bond length) of a = 1.78Å.It is shown to have metallic properties like most of other modifications of boron sheets.
The Fermi curve of a boron flat sheet consists of parts of 3 ellipses with semimajor and semiminor axes along the Γ-K and Γ-M directions, respectively.The effective electron mass at the Fermi level reveals a distinct anisotropy: m Γ−K /m 0 ≈ 0.480 and m Γ−M /m 0 ≈ 0.052, with m σ /m 0 ≈ 0.094 for conduction mass.The low effective mass of conduction electrons indicates a high mobility of electrons and, hence, a high conductivity of flat boron sheets.
The shapes of density-of-states obtained here for flat holeless boron sheets and previously calculated ones are rather similar, which is understandable as these structures are buckled/puckered or flat but with hexagonal holes, variants of the same triangular lattice.The remaining discrepancies may be attributed to the perturbations associated with the mentioned structural changes and differences in the used models.

Figure 1 .
Figure 1.A regular icosahedron B12 with B-atoms at the vertices.

Figure 2 .
Figure 2.An idealized unit cell of a β-rhombohedral boron crystal.

Figure 1 .
Figure 1.A regular icosahedron B 12 with B-atoms at the vertices.

Figure 1 .
Figure 1.A regular icosahedron B12 with B-atoms at the vertices.

Figure 2 .
Figure 2.An idealized unit cell of a β-rhombohedral boron crystal.

Figure 2 .
Figure 2.An idealized unit cell of a β-rhombohedral boron crystal.

Figure 3 .
Figure 3.The borophene structure on a silver substrate: the top and side views of the monolayer structure (unit cell indicated by the box) [16].

Figure 3 .
Figure 3.The borophene structure on a silver substrate: the top and side views of the monolayer structure (unit cell indicated by the box) [16].

1 .
Then, the radius-vectors of the lattice sites are a number of different values for the lattice constant of a boron sheet suggested theoretically.For the self-consistency of calculations, we use 37 .3 = a a.u. of length, i.e., 1.78 Å, which corresponds to the B-B pair interatomic potential in the same quasi-classical approximation [91].

Figure 4 .
Figure 4.A boron perfect flat sheet without hexagonal holes.

Figure 5 .
Figure 5.A 2-D rhombic unit cell of a boron flat sheet.

Figure 6 .
Figure 6.The vectors of a reciprocal lattice of a boron flat sheet.

Figure 4 .
Figure 4.A boron perfect flat sheet without hexagonal holes. Let a number of different values for the lattice constant of a boron sheet suggested theoretically.For the self-consistency of calculations, we use 37 .3 = a a.u. of length, i.e., 1.78 Å, which corresponds to the B-B pair interatomic potential in the same quasi-classical approximation [91].

Figure 4 .
Figure 4.A boron perfect flat sheet without hexagonal holes.

Figure 5 .
Figure 5.A 2-D rhombic unit cell of a boron flat sheet.

Figure 6 .
Figure 6.The vectors of a reciprocal lattice of a boron flat sheet.

Figure 5 .
Figure 5.A 2-D rhombic unit cell of a boron flat sheet. Let a number of different values for the lattice constant of a boron sheet suggested theoretically.For the self-consistency of calculations, we use 37 .3 = a a.u. of length, i.e., 1.78 Å, which corresponds to the B-B pair interatomic potential in the same quasi-classical approximation [91].

Figure 4 .
Figure 4.A boron perfect flat sheet without hexagonal holes.

Figure 5 .
Figure 5.A 2-D rhombic unit cell of a boron flat sheet.

Figure 6 .
Figure 6.The vectors of a reciprocal lattice of a boron flat sheet.

Figure 6 .
Figure 6.The vectors of a reciprocal lattice of a boron flat sheet.

Figure 7 .
Figure 7.The transform from the ) , ( 2 1 α α domain to the ) , ( y x k k domain of reciprocal space.

Figure 8 .
Figure 8.The hexagonal (first Brillouin zone) and rhombic unit cells of a reciprocal lattice of a flat boron sheet.

.
Dispute the shift of the B 1s atomic level, it retains an order of magnitude after transforming in an electronic band of the boron flat sheet.The lowest-lying band 1 E is fully filled with electrons.The band 2 E is the highest fully filled band (Figure9) with bottom at , this range of energies is comparable in order of magnitude with a valence 2s atomic level of eV

Figure 8 .
Figure 8.The hexagonal (first Brillouin zone) and rhombic unit cells of a reciprocal lattice of a flat boron sheet.

.
Dispute the shift of the B 1s atomic level, it retains an order of magnitude after transforming in an electronic band of the boron flat sheet.The lowest-lying band 1 E is fully filled with electrons.The band 2 E is the highest fully filled band (Figure9) with bottom at that, this range of energies is comparable in order of magnitude with a valence 2s atomic level of eV

Figure 8 .
Figure 8.The hexagonal (first Brillouin zone) and rhombic unit cells of a reciprocal lattice of a flat boron sheet.

Figure 9 . 3 E
Figure 9.The band 2 E energy surface (a) and contour plots (b) over a rhombic unit cell.

Figure 9 .
Figure 9.The band E 2 energy surface (a) and contour plots (b) over a rhombic unit cell.

Figure 9 .
Figure 9.The band 2 E energy surface (a) and contour plots (b) over a rhombic unit cell.

he band 3 E
(Figure10) is partially filled, i.e., partially empty, with a bottom at

Figure 10 .
Figure 10.The band 3 E energy surface (a) and contour plots (b) over a rhombic unit cell.

Figure 10 .Figure 10 .
Figure 10.The band E 3 energy surface (a) and contour plots (b) over a rhombic unit cell.

Figure 11 .
Figure 11.The band 4 E energy surface (a) and contour plots (b) over a rhombic unit cell.Between bands 1 E and 2 E , there is a very wide energy gap of eV 00 .239 max 1 min 2 12 = − = Δ E E E , while pairs of bands 2 E and 3 E , and 3 E and 4 E overlap each with other, i.e. there are obtained pseudo-gaps of eV 23 .3 max 2 min 3 23 − = − = Δ E E E and eV 49 .0 max 3 min 4 34 − = − = Δ E E E .

Figure 11 .
Figure 11.The band E 4 energy surface (a) and contour plots (b) over a rhombic unit cell.Between bands E 1 and E 2 , there is a very wide energy gap of ∆E 12 = E 2min − E 1max = 239.00eV, while pairs of bands E 2 and E 3 , and E 3 and E 4 overlap each with other, i.e. there are obtained pseudo-gaps of ∆E 23 = E 3min − E 2max = −3.23 eV and ∆E 34 = E 4min − E 3max = −0.49eV.The Fermi level is found at E Fermi = −19.42eV, within the part of the band E 3 without overlapping with other bands.This result confirms the metallicity of the boron sheet.Thus, all the electron energies are found to be negative.It means that all electrons, including conduction electrons at the Fermi level, are bounded inside the 2-D crystal.This result once more evidences the correctness of the calculations performed in this work.The total width of valence and conduction bands equal to ∆E V = E Fermi − E 2min = 17.79 eV and ∆E C = E 4max − E Fermi = 11.34 eV, respectively.The upper valence band width is ∆E VU = E Fermi − E 2max = 0.43 eV.As expected, it is negligible if compared with that of a conduction band.To compare easily our results with the literature data, in addition to the presentation of the band structure using the contour plots of the whole Brillouin zone in Figures12 and 13, we plot the band energies (as well as their second derivatives and corresponding parabolic approximations) along the main lines of symmetry.Our quasi-classical calculation of the crystalline band structure, like any other approach also utilizing HF parameters of constituting atoms, cannot determine the absolute values of energy parameters with a high accuracy.By this reason, the above mentioned value E Fermi cannot be used directly to determine the electron work function of the boron sheet.This goal can be achieved only after corrections are made to include the electron-correlation and to exclude the electron-self-interacting effects, which have to allow an accurate determination of the position for the vacuum level of energy

Figure 12 .
Figure 12.The section of the conduction band surface along the main diagonal of a rhombic unit cell (direction Γ-Κ) of reciprocal space (in atomic units).

Figure 13 .
Figure 13.The section of the conduction band surface along a small diagonal of a rhombic unit cell (direction Γ-Μ) of reciprocal space (in atomic units).

Figure 12 .
Figure 12.The section of the conduction band surface along the main diagonal of a rhombic unit cell (direction Γ-K) of reciprocal space (in atomic units).

Figure 12 .
Figure 12.The section of the conduction band surface along the main diagonal of a rhombic unit cell (direction Γ-Κ) of reciprocal space (in atomic units).

Figure 13 .
Figure 13.The section of the conduction band surface along a small diagonal of a rhombic unit cell (direction Γ-Μ) of reciprocal space (in atomic units).

Figure 13 .
Figure 13.The section of the conduction band surface along a small diagonal of a rhombic unit cell (direction Γ-M) of reciprocal space (in atomic units).

Figure 14 .
Figure 14.Curves of the intersection of the band surface with the Fermi plane in neighboring rhombic unit cells of reciprocal space.
number components along perpendicular axes Γ-Κ and Γ-Μ and eV equation if it is rewritten in the form of a normalized ellipse equation axes in the directions Γ-Κ and Γ-Μ, respectively:

Figure 14 .
Figure 14.Curves of the intersection of the band surface with the Fermi plane in neighboring rhombic unit cells of reciprocal space.
The effective electron mass at the Fermi level reveals a significant anisotropy.For the central ellipse, the effective masses are 480

Figure 15 .Figure 15 .
Figure 15.The Fermi curve of a boron flat sheet.DoS within the bands 2 E , 3 E , and 4 E against electron energy renormalized to the Fermi level Fermi E E E − → are presented in Figure 16 in two different scales for the convenient consideration.As for the band 1 E , DoS within this band is proportional to the Diracfunction, [eV]) 42 .256 ( ~+ E δ , with the accordingly renormalized argument eV 79 .256 eV) 42 .19 ( eV 21 .276 eV 21 .276 = − − − → − .

Figure 15 .Figure 16 .
Figure 15.The Fermi curve of a boron flat sheet.DoS within the bands 2 E , 3 E , and 4 E against electron energy renormalized to the Fermi level Fermi E E E − → are presented in Figure 16 in two different scales for the convenient consideration.As for the band 1 E , DoS within this band is proportional to the Diracfunction, [eV]) 42 .256 ( ~+ E δ , with the accordingly renormalized argument eV 79 .256 eV) 42 .19 ( eV 21 .276 eV 21 .276 = − − − → − .

Figure 16 .
Figure 16.The density-of-electron-states renormalized to the Fermi level in a valence band and the lower and upper conduction bands of a boron flat sheet in two different scales: general view (a) and in Fermi level vicinity (b).

gure 16 .
The density-of-electron-states renormalized to the Fermi level in a valence band and the wer and upper conduction bands of a boron flat sheet in two different scales: general view (a) and Fermi level vicinity (b).
is shown to have metallic properties like most of modifications of boron sheets.
is shown to have metallic properties like most of other modifications of boron sheets.
1 , α 2 ) and H(α 1 , α 2 ) are 16 × 16 matrices of overlapping integrals and single-electron Hamiltonian, respectively, reducible to 4 × 4 matrices.E(α 1 , α 2 ) is a required electron energy band.About the parameters α 1 and α 2 , see below.This equation has 4 different real and negative roots E m (α 1 , α 2 ), m = 1, 2, 3, 4. It can be demonstrated that they exhibit all the different solutions of the corresponding secular equation with 16 × 16 matrices.Within the initial quasi-classical approximation, these matrix elements can be found from the relations shown in the Appendix A.

Table 1 .
The inner and outer classical turning point radii r i and r i of electrons in boron atom.

Table 2 .
The radii r λ of radial layers of quasi-classical averaging of potential in boron atoms and the averaged values of potential ϕ λ .

Table 3 .
The band widths and (pseudo)gaps between bands.