Mysterious SiB3: Identifying the Relation between α- and β-SiB3

Possible boron-rich amorphous silicon borides from ab initio simulations

CONTEXT

By means of ab initio molecular dynamics simulations, possible boron-rich amorphous silicon borides (B_nSi_1-n, 0.5 ≤ n ≤ 0.95) are generated and their microstructure, electrical properties and mechanical characters are scrutinized in details. As expected, the mean coordination number of each species increases progressively and more closed packed structures form with increasing B concentration. In all amorphous models, pentagonal pyramid-like configurations are observed and some of which lead to the development of B₁₂ and B₁₁Si icosahedrons. It should be noted that the B₁₁Si icosahedron does not form in any crystalline silicon borides. Due to the affinity of B atoms to form cage-like clusters, phase separations (Si:B) are perceived in the most models. All simulated amorphous configurations are a semiconducting material on the basis of GGA+U calculations. The bulk modulus of the computer-generated amorphous compounds is in the range of 90 GPa to 182 GPa. As predictable, the Vickers hardness increases with increasing B content and reaches values of 25-33 GPa at 95% B concentration. Due to their electrical and mechanical properties, these materials might offer some practical applications in semiconductor technologies.

METHOD

The density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations were used to generate B-rich amorphous configurations.

First-principles calculation of 11B solid-state NMR parameters of boron-rich compounds II: the orthorhombic phases MgB7 and MgB12C2 and the boron modification γ-B28

Based on the work on referencing 11B nuclear magnetic resonance (NMR) spectra for molecular icosahedral boranes and the subsequent transfer to the rhombohedral boron-rich borides of the α-rB12 type, we show that the magic angle spinning (MAS) NMR spectra of boron-rich borides with four or five symmetry-independent boron atoms can also be calculated. The calculations are performed on the level of density functional theory (DFT) using the gauge-including projector-augmented wave (GIPAW) approach. As model compounds o-MgB12C2 and MgB7 are used, for which the experimental spectra could be calculated in excellent agreement with a deviation of 1 to 2 ppm. Based on the calculations, the different B atoms can be assigned to the respective signals, taking into account the quadrupolar coupling constants Cq from computation of the electric field gradient (EFG) with its main axis Vzz. It is shown that due to the specific geometric conditions of icosahedra, the magnitudes of Vzz for the boron atoms involved in exohedral B-B bonds to neighbouring icosahedra depend only on the valence electron density of the bond critical point and the distance. This also applies to the bonds to the interstitial B2 unit in MgB7, but not to bonds to the heteroatom of the C2 dumbbell in o-MgB12C2. Both results are in line with our previous observations for the rhombohedral species (α-rB12; B12X2 with X = P, As, O). Finally, the spectrum of γ-B28 was calculated, whose structure also contains B12 icosahedra and interstitial B2 dumbbells. Here, a very similar bonding situation is found for the icosahedron, but the calculations show that the situation for the B2 unit is clearly different. In general, the only parameter that needs to be varied to fit calculated and measured spectra is the linewidth, as this cannot be calculated. For the cases of o-MgB12C2 and MgB7 signal areas are related to corresponding site multiplicities. A prerequisite for the successful application of the chosen method seems to be the presence of a semiconductor with a sufficiently large band gap, which is the case for the compounds investigated.

First-principles calculation of 11B solid-state NMR parameters of boron-rich compounds I: the rhombohedral boron modifications and B12X2 (X = P, As, O)

In the present study, solid-state nuclear magnetic resonance (NMR) spectra under magic angle spinning conditions of the rhombohedral structures α-B and B12P2 are reported together with the corresponding parameter sets from first principles calculations on α-B B12X2 (X = P, As, O). With the combination of density functional theory (DFT) and the gauge-including projector-augmented wave (GIPAW) approach as the theoretical tools at hand the computed 11B parameters lead to unambiguous explanation of the measurements. Thereby, we overcome common obstacles of processing recorded NMR spectra of solid-state compounds with several crystallographic positions, in particular non-trivial signal assignments and parameter determination due to peak overlap or even unexpected intensity/area ratios. In fact, we find very good agreement between the theoretical results and measured spectra without applying fitting procedures. Using the Perdew-Burke-Ernzerhof (PBE) functional, the results of the common construction types for pseudopotentials and referencing methods for the chemical shift determination are compared. Suggestions and conclusions from experimental 11B NMR studies on parameters according to the icosahedral positions are critically discussed, for instance the early suspected correlation to chemical shifts is not confirmed. Regarding the electric field gradient (EFG) a detailed explanation for obtaining small deviations amongst all investigated structures of the icosahedral polar sites compared to the equatorial sites is given. Our results show an important link between the exohedral bonding situation of compounds with icosahedral structure elements and the main axis of the EFG and therefore, also measurable quadrupole coupling constants if certain geometrical conditions are fulfilled. Finally, this work also contributes to establishing the number of unique sites measured by solid-state NMR methods within the modification of β-B.