The structure and hardness of the highest boride of tungsten, a borophene-based compound

Two-dimensional systems have strengthened their position as a key materials for novel applications. Very recently, boron joined the distinguished group of elements confirmed to possess 2D allotropes, named borophenes. In this work, we explore the stability and hardness of the highest borides of tungsten, which are built of borophenes separated by metal atoms. We show that the WB3+x compounds have Vickers hardnesses approaching 40 GPa only for small values of x. The insertion of extra boron atoms is, in general, detrimental to the hardness of WB3 because it leads to the formation of quasi-planar boron sheets that are less tightly connected with the adjacent tungsten layers. Very high concentrations of boron (x ≈ 1), give rise to a soft (Vickers hardness of ~8 GPa) and unstable hP20-WB4 structure that can be considered to be built of quasi-planar boron α-sheets separated by graphitic tungsten layers. By contrast, we show that the formation of tungsten vacancies leads to structures, e.g. W0.75B3+x, with Vickers hardnesses that are not only similar in value to the experimentally reported load-independent hardnesses greater than 20 GPa, but are also less sensitive to variations in the boron content.

The highest boride of tungsten-often referred to as tungsten tetraboride-is recently best explored for its potential applications as superhard material, however made its first appearance in the literature in 1961, when Chretien and Helgorsky 1 did the first attempt to find its structure. Years later Romans and Krug 2 reported that WB 4 has a hexagonal structure of 20 atoms per unit cell and lattice constants of 5.2 and 6.34 Å for a and c, respectively. The space group of this structure was determined to be P6 3 /mmc. The hP20-WB 4 structure serves now as a reference structure for almost all subsequent experimental studies related to boron-rich materials with a WB 4 -like structure [3][4][5][6][7] . The mechanical properties of WB 4 were first determined by Gu et al. 3 who reported Vickers hardness (H V ) values of 46.2 and 31.8 GPa under applied loads of 0.49 and 4.90 N, respectively, measured by the microindentation technic. Subsequently, Mohammadi et al. 5 also measured the hardness by microindentation method and reported H V values of 43.3 and 28.1 GPa at low (0.49 N) and high (4.90 N) loads, respectively. More recently, the Vickers hardness for W 0.85 B 3 was reported by Tao et al. 8 to be 42.0 and 25.5 GPa under applied loads of 0.098 and 4.90 N, respectively. Finally, Lech et al. 9 determined the maximum nanoindentation hardness of W 0.82 B 3.54 (at a penetration depth of 95.25 nm) to be 41.7 GPa. It is generally accepted that a reliable hardness of a material can be determined from the asymptotic hardness region achieved at high loads 10 . The quite large differences between the H V values reported for WB 4 , especially for high loads, can be attributed to differences in the amount of boron contamination and/or presence of tungsten vacancies, which were experimentally seen in the studied samples 9 . By exploring structures with different compositions, we are able to explain on the theoretical ground the apparent differences between the reported experimental results.
The common description of hP20-WB 4 that can be found in the literature is that this structure consists of graphitic boron layers separated by graphitic layers of W atoms like in the hP16-WB 3 structure but with additional B 2 dimers located between boron sheets and aligned along the c-axis (see Fig. 1a). This description, although very elegant, is completely decoupled from more recent investigations related to 2D boron crystals 11,12 . An 'updated' view to hP20-WB 4 would be that it is a structure consisting of a sequence of quasi-planar boron α-sheets (see Fig. 1b) separated by graphitic W layers. Extensive theoretical investigations have proved, however, that the stoichiometric WB 4 compound in the hP20-WB 4 structure would be thermodynamically and dynamically unstable. Its calculated enthalpy of formation is positive, with a value of 0.4 eV/atom 13 , and from its phonon dispersion the structure was shown to be highly unstable 14 . In fact, we argue that the formation of stable quasi-planar boron layers within the hP20-WB 4 structure is the main reason for the instability of this and related boron-rich structures. Thus, from this viewpoint, the experimentally obtained WB 4 compounds cannot adopt the hP20-WB 4 structure.
In more recent reports 13,15 , the highest borides of tungsten are described as hP16-WB 3 structures contaminated with additional boron atoms. The exact position of the boron atoms in the crystal lattice are difficult to be determined experimentally because of the large mass difference between W and B atoms 15 . Therefore, the combination of theory and experiment is essential for the understanding of the observed findings. Since in the experiment WB 3 is not only contaminated with boron atoms but also, to some extent, possesses tungsten vacancies 8,9,16 , in this work a more precise notation is used when referring to the highest boride of tungsten, namely W 1-y B 3+x , to underline the existence of W vacancies and explore their influence on the stability and properties of WB 3+x .

Results
The structure of W 1-y B 3+x . The highest borides of tungsten are obtained starting from WB 3 in the hP16-WB 3 structure by adding additional boron atoms at the positions shown in red in Fig. 1c and/or by selective removal of W atoms. The fully 'packed' structure is the hP20-WB 4 structure, shown in Fig. 1a, that has buckled boron α-sheets, shown in Fig. 1b, separated by W layers. The buckling height is 1.49 Å and is larger than that of the freestanding triangular boron sheet (0.82 Å) 17 . The complete removal of all the W atoms from hP20-WB 4 leads to the all-boron structure that consists of 16 atoms per unit cell and is shown in Fig. 1d. This structure is nothing more than a sequence of quasi-planar boron α-sheets arranged in such a way that the boron atoms that stick out of the graphitic frames face each other forming dimers. We can describe the hP20-WB 4 structure in the same manner except that the boron α-sheets are intercalated by tungsten layers. The hP16-B structure is less stable than the α-rhombohedral boron (hR12-B) by 0.47 eV/atom. Interestingly enough, by removing one of the boron dimers in hP16-B a slightly more stable (by 12 meV/atom) hP14-B structure is obtained. The hP16-B and hP14-B structures have P6 3 /mmc and P6/mmc space groups, respectively, and a = 5.034 Å, c = 6.166 Å and a = 5.081 Å, c = 5.195 Å lattice constants, respectively. Stability of the compounds. To explore the relative stability of the generated structures, we calculate for each structure its enthalpy of formation per atom, ΔE. The ΔE values are calculated relative to the enthalpies of tungsten and boron solids in the body-centered cubic (cI2-W) and α-rhombohedral (hR12-B) structures, respectively. The results for ΔE versus boron atomic content, (3 + x)/(4 + x − y), are summarized in Fig. 2a and b. All the structures that have enthalpies of formation above the horizontal dashed lines in Fig. 2a and b are, in principle, thermodynamically unstable. It is instructive, however, to draw also a line that connects cI2-W with hP16-B, what is shown in Fig. 2a and b by a dashed line that is above the horizontal dashed line. The enthalpies of formation of almost all the considered W 1-y B 3+x structures are located bellow the cI2-W ↔ hP16-B line. This means that the incorporation of W atoms in between boron sheets is energetically favorable. The enthalpies of formation versus boron atomic content are presented in two ways. In Fig. 2a, we organize the results according to the number of extra boron atoms, n B , in W 1-y B 3+x relative to hP16-WB 3 . In Fig. 2b, the same results have been organized emphasizing the number of W vacancies in W 1-y B 3+x also relative to hP16-WB 3 . It is clear from Fig. 2a that negative enthalpies of formation or positive ΔE but close to 0, have structures with none or no more than 2 extra B atoms. The cases with 3 and 4 extra B atoms have ΔE  0.2 eV/atom or larger. This also includes the highly debated hP20-WB 4 structure (ΔE = 0.36 eV/atom). From Fig. 2b, we can learn that the only relevant cases are those for which the number of W vacancies, n V(W) , is 0 or 1, since among those cases we can find structures with negative or close to 0 enthalpies of formation. Combining all the information coming from Fig. 2, we choose 8 structures that in principle can be important to understand experimental results. Six of those structures are shown in Fig. 3. For each relevant boron atomic content, we choose the structure with the lowest enthalpy of formation. The relevant structures (the highest boride of tungsten) are those with boron contents ranging from 0.75 (WB 3 ) to 0.83 (W 0.75 B 3.75 ). The highest boron content is chosen following ref. 9.
The lattice constants and symmetry of the structures shown in Fig. 3 are summarized in Table 1. In this table, we also include, for each structure, the occupations of the B and W atoms relative to the P6 3 /mmc space group. It is important to notice that if we take the average of all the lattice constants listed in Table 1, we get 5.192 and 6.303 Å for a and c, respectively, that is, values that match quite well those reported in the experiment (5.2 and 6.34 Å for a and c, respectively) 9 . This may suggest that in the experiment is observed a non-stoichiometric structure with a random distribution of both the extra boron atoms and W vacancies and may farther explain the difficulties in the interpretations of X-ray and neutron diffraction data 9 .
Mechanical properties of the compounds. The elastic properties of the studied structures are summarized in Table 2, whereas the plot of the Vickers hardness versus boron atomic content is shown in Fig. 4. From this figure, we see that only stoichiometric WB 3 (hP16-WB 3+x ) can be considered as superhard material. The hardness is however affected by contamination by extra B atoms. This is clearly seen in Fig. 4 for WB 3+x , for which the Vickers hardness changes from ~40 to ~8 GPa for an increase of the B content by 5%. For the tungsten-deficient W 0.75 B 3+x structures the picture is different, namely, we obtain smaller variations of the Vickers hardness with the increase of B content (see Fig. 4). Most of the considered structures have Vickers hardnesses larger or equal to 20 GPa, what means that the highest boride of tungsten is a hard material but not superhard (at least in the range of considered boron contents). A particularly soft structure is hP20-WB 4 , which has a comparable bulk modulus to that of hP16-WB 3 but much smaller shear modulus (see Table 2). The softening of WB 4 may be attributed to the formation of stable 2D boron layers (α-sheets) which are less tightly bound to the tungsten layers. The average nearest neighbor W-B distance is 2.324 and 2.383 Å in WB 3 and WB 4 , respectively, what reflects the weakening of the W-B bond in WB 4 with respect to WB 3 .
Another way of explaining the lack of stability of the hP20-WB 4 structure could be by using the concept of self-doping in boron sheets introduced in ref. 18. By adding extra boron atoms to the hP16-WB 3 structure, we essentially transfer extra negative charge into the boron graphitic layers present in hP16-WB 3 . This extra   Table 1. Phases, occupations relative to the P6 3 /mmc space group, lattice constants, and resulting space groups of the most stable structures at a given boron atomic content.
charge influences the B-B distances in these layers. In Fig. 1a,c, and d, we show the B-B nearest neighbor distances between boron atoms for the structures present there. The B-B distances (1.739 Å and 1.794 Å) for the in-plane boron atoms in hP20-WB 4 are clearly longer than the respective B-B distances (1.708 Å and 1.732 Å) in hP16-WB 3 . This means that the graphitic layers in hP20-WB 4 are overcharged with electrons. Also the inclusion of tungsten atoms into the hP16-B structure seems to do the same effect as boron doping, since the B-B distances (1.696 Å and 1.669 Å) for the in-plane atoms in hP16-B are shorter than the respective distances in hP20-WB 4 . The out-of-plane boron atoms in hP20-WB 4 stay as distant as possible from the graphitic boron layers in order to minimize the transfer of negative charge to the in-plane boron atoms. As a consequence, we have a formation of B 2 dimers that are loosely bound to the graphitic frames.
In summary, we show that the insertion of extra boron atoms into the WB 3 structure is, in general, energetically unfavorable and lowers its shear modulus while maintaining a high value for the bulk modulus, effectively leading to a softer material. A high degree of boron contamination leads to the formation of quasi-planar boron α-sheets separated by graphitic tungsten layers in WB 4 . Structures of the W 1−y B 3+x type, in which boron contamination is accompanied by the presence of tungsten vacancies, are more stable and harder than WB 4 . Finally, the formation of tungsten vacancies gives rise to structures (e.g. W 0.75 B 3+x ) with Vickers hardnesses that are less sensitive to variations in the boron content and are similar in value to the experimentally reported load-independent values, which are greater than 20 GPa. Our results should provide guidance for the development of new WB 4 synthesis strategies.

Methods
Our first principles calculations were based on density functional theory (DFT) and the projector augmented wave (PAW) method as implemented in the Quantum ESPRESSO simulation package 19 . For the exchange and correlation functional, we used the revised Perdew-Burke-Ernzerhof spin-polarized generalized gradient approximation (PBEsol-GGA) functional. The plane-wave basis set was converged using a 40 Ry energy cutoff. A 8 × 8 × 8 k-point mesh and a Gaussian smearing of 0.005 Ry was used in the Brillouin zone integration. The calculations were carried out using supercells containing up to 20 atoms. For each considered structure, a full atomic position and lattice parameter relaxation was preformed.
A total of 60 low-energy W 1−y B 3+x structures with high boron atomic content were selected using the cluster-expansion method implemented in the Alloy-Theoretic Automated Toolkit (ATAT) 20 . The enthalpies of formation were calculated from the formula:  E tot are the calculated total energies per atom and c = (3 + x)/(4 + x − y) is the boron atomic content. The elastic properties of the most stable structures were calculated using the ElaStic code 21 . To compute the Vickers hardness, we used the semi-empirical hardness model proposed by Chen et al. 22  , where G and k are the shear modulus and the Pugh modulus ratio (k = G/B, where B is the bulk modulus), respectively.