Structural stability and mechanism of compression of stoichiometric B13C2 up to 68 GPa
Irina Chuvashova, Elena Bykova, Maxim Bykov, Volodymyr Svitlyk, Leonid Dubrovinsky, Natalia Dubrovinskaia
11 Structural stability and mechanism of compression of stoichiometric B C up to 68 GPa Irina Chuvashova a,b* , Elena Bykova b , Maxim Bykov b , Volodymyr Svitlyk c , Leonid Dubrovinsky b , Natalia Dubrovinskaia a a Material Physics and Technology at Extreme Conditions, Laboratory of Crystallography, University of Bayreuth, D-95440 Bayreuth, Germany b Bayerisches Geoinstitut, University of Bayreuth, D-95440 Bayreuth, Germany c European Synchrotron Radiation Facility, BP 220 F-38043 Grenoble Cedex, France *Correspondence to: [email protected], [email protected]
Abstract
Boron carbide is a ceramic material with unique properties widely used in numerous, including armor, applications. Its mechanical properties, mechanism of compression, and limits of stability are of both scientific and practical value. Here, we report the behavior of the stoichiometric boron carbide B C studied on single crystals up to 68 GPa. As revealed by synchrotron X-ray diffraction, B C maintains its crystal structure and does not undergo phase transitions. Accurate measurements of the unit cell and B icosahedra volumes as a function of pressure led to conclusion that they reduce similarly upon compression that is typical for covalently bonded solids. A comparison of the compressional behavior of B C with that of α–B, γ–B, and B C showed that it is determined by the types of bonding involved in the course of compression. Neither ‘molecular-like’ nor ‘inversed molecular-like’ solid behavior upon compression was detected that closes a long-standing scientific dispute.
Introduction
Boron carbide B C was discovered in 1858, and its crystal structure was first established in 1934. It belongs to icosahedral boron compounds, a family of crystalline solids with crystal structures based on various arrangements of B icosahedra, which are considered to be some kind of B “molecules”. In such solids chemical bonding has been rationalized in terms of polycenter bonds on the B closo - cluster (boron icosahedron) and two-electron-two-center ( ) and two-electron-three-center ( ) bonds between the clusters . According to the Wade-Jemmis rule
26 out of 36 valence electrons of a B cluster are accommodated in 13 molecular-orbital-like bonding orbitals to form the cluster. This leaves 10 electrons for external bonding using 12 equivalent external bonding orbitals, thus creating an electron deficiency in the B cluster. For a long time, this ‘electron deficiency’ of the intraicosahedral bonding relative to the intericosahedral one was expected to make the icosahedra more compressible compared to the unit cell, contrary to what is known for typical molecular solids. For this reason the icosahedral boron-rich solids were understood as "inverted-molecular" solids . This gave the origin to the problem of the ratio of the rigidity of the B closo - cluster and the unit cell and the mechanism of compression in icosahedral boron-rich solids that had to be proven experimentally. The first work, which aimed to shed light into the mechanism of compression of icosahedral boron-rich solids in general and boron carbide in particular, was a high-pressure powder neutron diffraction study of boron carbide (B C) to 11 GPa conducted by Nelmes and co-authors . By the linear fit of the experimental ‘pressure versus volume’ (P-V) data obtained for both the icosahedron and the unit cell of B C, the bulk moduli of the icosahedron and the entire structure were determined to be 169(6) GPa and 199(7) GPa, respectively. This led to conclusion that the icosahedron is 23(4) % more compressive than the unit cell of boron carbide, which thus behaved as an "inverted-molecular" solid . An experimental evidence that icosahedral boron materials on compression may not behave as "inverted-molecular" solids came from the single-crystal X-ray diffraction data to 65 GPa obtained for the high-pressure boron allotrope, γ-B . The bulk modulus of the B icosahedron was determined to be 285(6) GPa, whereas the bulk modulus of the entire γ-B structure was found to be 227(3) GPa . A new insight into the situation in boron carbide, compared to the very first measurements , was given due to a compressibility study of ‘nearly stoichiometric boron carbide B C’ . Single crystals of B C were investigated using synchrotron X-ray diffraction in a diamond anvil cell to 74 GPa. The structure of B C was understood as consisting of B C icosahedra interconnected by the C-B-C chains. The parameters of the equation of state (EoS) of boron carbide were K =243(6) GPa (K´=3.6(2)) or 236(8) at fixed K´=4 ( K is the bulk modulus, and K´ is its pressure derivative; the zero pressure unit cell volume was fixed on the ambient pressure experimental value). Dera et al. pointed out that the “icosahedron volume compression did not follow a typical EoS functional behavior” and they did not calculate the bulk modulus of the B C icosahedron. Nevertheless, they reported that the B C icosahedron showed a 13% volume reduction, which was smaller than that of the unit cell volume (18 %). The general conclusion was that B C behaves as a molecular solid that is in accordance with the results obtained for γ-B . Optical properties and evolution of Raman modes of the same B C samples, which were characterized and studied using single-crystal X-ray diffraction in Ref. 10, were also investigated up to 70 GPa by the same research group . Based on their spectroscopy data, the authors reported high-pressure phase transition in boron carbide at 40 GPa, although, according to Dera et al. , no signatures of structural phase transitions were observed in B C by high-pressure XRD studies up to 70 GPa . Verifying these mutually contradictive reports on boron carbide is not only scientifically important on itself, but it is additionally justified in connection to the commercial value of this unique material. Boron carbide is a very hard and, at the same time, lightweight material for applications in personal security (bullet-proof vests) . It possesses the highest Hugoniot elastic limit of ceramic materials (ca. 17–20 GPa), i.e. the maximum uniaxial dynamic stress that the material can withstand elastically, surpassing all its denser competitors such as silicon carbide and alumina by a factor of 2 (Ref. 12). However, it fails just above the Hugoniot elastic limit and the possible source of failure could be clarified through high-pressure experiments. Thus, establishing the mechanism of compression of boron carbide, clarifying its mechanical properties and limits of its stability under loading, are of both scientific and practical interest. Here we report high-pressure investigations of stoichiometric boron carbide B C using high-pressure single-crystal X-ray diffraction up to 68 GPa. Single crystals of B C , which we study here, were characterized in detail in previous work of our group . As it was established that B C is fully ordered and stoichiometric, and carbon atoms occupy a single position (at the ends of the C-B-C chains) , in the present work we could track all changes of the crystal structure, atomic positions, and bond lengths with the high accuracy up to 68 GPa. This enabled us to establish the equations of state of both the B C crystal and the B icosahedron. As the B icosahedra in the stoichiometric boron carbide B C do not contain any carbon atoms, contrary to all previously investigated boron carbides
8, 10, 14-24 , we could compare the compressional behavior of B closo -clusters in B C and boron allotropes α-B and γ-B and conclude not only regarding the ratio of the rigidity of the B closo -cluster and the unit cell in these materials, but also concerning the mechanism of compression of boron-rich icosahedral solids. Results
The equations of state of B C and B icosahedra The single-crystal X-ray diffraction data for B C obtained at variable pressure and some experimental details are presented in Table 1. The quality of the data allowed the refinement of both the lattice parameters and atomic coordinates (Supplementary Table S1). The quality of the structural refinement was good up to the highest pressures achieved that gives evidence that in the course of compression the crystals were maintained in quasihydrostatic environment. In the structure of boron carbide B icosahedra are located in the corners of the rhombohedral cell, and intericosahedral three-atom C-B-C linear chains are oriented along its body diagonal (Fig. 1a, right). As seen in Figure 1, the structure of B C is very similar to those of α-B and γ-B, which can be described in terms of a cubic closest packing ( ccp ) of spheres, where B icosahedra play the role of “spheres” (Figure 1a,b). The unit cell parameters of B C in hexagonal setting are a = 5.5962(3) Å, c = 12.0661(7) Å, (space group R m), as determined in . Boron atoms in the crystal structure of B C (Figure 1) occupy three crystallographically independent positions (B P , B E , and B C ) and the forth position is occupied by carbon atoms in the C-B-C chains . In the present paper we adopt the nomenclature introduced by Ref. 13 : boron atoms forming the boron icosahedra are labeled as B P (polar positions) and B E (equatorial positions), B C designates the boron atom in the center of the C-B-C chain (Figure 1). In this study, all observed reflections perfectly match the B С structure up to the highest pressure reached. Our X-ray diffraction data did not reveal any indication of a phase transition. All the unit cell parameters smoothly decrease on compression (Table 1). Figure 2 presents the dependence of the relative unit cell parameters ( a/a and c/c ) and the relative unit cell ( V/V ) volume of B C on pressure up to 68 GPa. As seen (Figure 2), the structure of B С is slightly more compressible along the c direction. The whole experimental volume-pressure data set was fitted using the third-order Birch-Murnaghan (3BM) equation of state that gave the following EoS parameters: V = 327.25(4) Å , K = 239(7) GPa and K´ = 3.2(3) ( V is the zero pressure unit cell volume, K is the bulk modulus, and K´ is its pressure derivative) (Table 2). The fit with the fixed K´ = 4 resulted in the bulk modulus of K = 222(2) GPa, being lower than that with the free K´. The evolution of the relative volume of the B icosahedron ( V ico /V ico0 ) with pressure is presented in Figure 2. It is very similar to that of the unit cell volume. The whole experimental data set of the icosahedra volumes versus pressure (Table 3) was fitted using the 3BM EoS. The EoS parameters were obtained as follows: V = 12.50(3) Å , K ico = 239(23) GPa, and K ico ´ = 3.8(8) ( V is the zero pressure icosahedron volume, K ico is the bulk modulus of icosahedra, and K ico ´ is its pressure derivative). The fit with the K ico ´ = 4 (fixed) results in a very close value of K ico = 234(7) GPa. A comparison of the bulk moduli found for the B С crystal and the B icosahedron (Table 2) shows that the bulk B С and the icosahedron have the similar rigidity. Evolution of the bonds lengths on compression of B С Recent experimental electron-density study using low-temperature high-resolution single-crystal synchrotron X-ray diffraction data clarified the bonding situation in the stoichiometric boron carbide B С . In the present work, we investigated single crystals from the same batch. For consistency, we adopt here the notations of bonds introduced in Ref. 13. There are seven distinct bonds in the structure of boron carbide (Table 3), which get into three groups (see Figure 1d, bottom, for bonds notations): intra-cluster polycentral bonds (B P –B P , B E –B E , B P –B E and B P –B E ); inter-cluster bonds (B P –B P ), which connect atoms in the polar sites (B P ) of the neighboring icosahedra; and bonds involving C-B-C chains (C–B E –and C– B C ). Single-crystal X-ray diffraction data, collected at eleven pressure points in the interval from 4 GPa to 68 GPa, allowed us to follow changes in the length of each of the seven bonds in the structure of boron carbide B С (Table 3). All bonds gradually shorten under compression (Figure 3); their pressure-dependent evolution does not show any anomalies (Figure 3).
Discussion
As mentioned in the introduction, bulk compressibility of boron carbide compared to compressibility of icosahedra has been a matter of debate. Nelmes et al. reported the crystal structure to be more rigid than the icosahedron cluster, whereas Dera et al. observed an opposite relation. Our results have shown that the rigidity of the crystal structure of stoichiometric boron carbide B C is similar to that of the B icosahedron. Within the standard uncertainty, the bulk moduli of the bulk material and the icosahedra we found to be similar: K = 239(7) GPa (K´ = 3.2(3)) versus K ico = 239(23) (K ico ´ = 3.8(8)) (with the fixed K´ = 4, K = 222(2) GPa versus K ico = 234(7) GPa). The volume reduction of the unit cell of B C and the volume reduction of the icosahedron in the pressure interval between ambient and 60 GPa were calculated to be also similar (18.7% versus for the crystals of B C is in a very good agreement with our result: their 3BM EoS parameters are K = 243(3) GPa, K´ = 3.6(6), and V =330.59(5) Å (V was fixed at experimental value obtained at ambient conditions); with the fixed K´= 4, they got K=236(8) GPa . The volume reduction of the unit cell of B C (18%) matches very well to what we obtained for B C . But interestingly, the volume reduction of the icosahedron in B C in the same pressure interval (between ambient and 60 GPa) in Ref. 10 appeared to be different (13%) (see Table 2). To clarify a reason of the apparent difference in the ratio of the rigidity of the unit cell and icosahedra found for B C and B C , we first compared the icosahedron and the unit cell volume evolution with pressure in the stoichiometric boron carbide B C and boron allotropes α-B and γ-B . In the structures of each of these materials, icosahedra (B ) are built of exclusively boron atoms. Table 2 demonstrates a remarkable observation: in the same pressure range (from ambient up to 60 GPa), the unit cell volume reductions for all these materials are similar (ca. 18% within less than 1% deviation). The B C, containing B C icosahedra, is not an exception (18%) . However, the icosahedra volume reductions are all different and reduce in the raw B C (18.1%), γ-B (16.9%), α-B (14.5%), and B C (13%). This observation is striking enough and desires an explanation through an insight into the compressional behavior of individual bonds in these solids. To visualize the difference in the rates of changes of the bonds length and compare boron carbide B C with α-B and γ-B, experimentally obtained data for the relative changes of the bond lengths ( l P /l P0 ) versus pressure was linearly fitted for all the bonds ( l P is the length of the bond at pressure P; l P0 is the length of this bond at P = 4.0 GPa, the first pressure point available in our experiment in the DAC). We plotted calculated “line slopes” against corresponding interatomic distances l P0 (Figure 4) at the lowest pressure, similarly to how it was done for characterization of the bond lengths’ change under pressure for various boron-rich compounds and α-B . As seen in Figure 4, all the points corresponding to the inter-cluster bonds (C– B C , C–B E , and B P –B P ) in B C (black squares) lie on one line, and the rate of their compression depends on the initial length of the bonds (the C– B C bond between carbon and boron atoms of the C– B C –C chain is the shortest one, see Figure 1a, right panel; Figure 1d, bottom panel; Figure 1b, right panel). Points corresponding to inter-icosahedral bonds in α-B (green circles in Figure 4) drop on the same line: initially shorter B P –B P bonds contract slower than the initially longer B E –B E contacts. All points corresponding to the intra-cluster B-B contacts, i.e. those involved in formation of the polycentral bonds of the B closo -cluster, for both B C and α-B (orange squares and red circles, correspondingly, in Figure 4) appear in a very compact area in Figure 4, indicating similar rates of contraction and their similar lengths at ambient conditions. Why then the B icosahedra in B C and α-B undergo such a dramatically different volume reduction (18.1 vs
13, 27, 28 , the validity of the Wade-Jemmis model (see introduction) was demonstrated for α–B , γ–B , and B C . In all these solids, the molecular-orbital-type bonding on the icosahedral clusters leaves for exo -cluster bonding twelve sp hybrid orbitals perpendicular to the surface of the clusters. Thus, in terms of bonding, we can consider B closo -clusters to be similar in these three substances. Concerning the inter-cluster bonds, they were found to be very different in α–B , γ–B , and B C . In α–B and B C the B P –B P bonds connecting polar B P atoms (those between neighboring close-packed layers of B icosahedra) are strong covalent 2e2c bonds (see Figure 1a-c, left and right panels; Figure 1d, top and bottom panels). The B E –B E contacts in α–B are very flexible, because they are a part of relatively weak 2e3c B E –B E –B E bonds (shown by brown triangles in Figure 1). Instead of these weak bonds, in the same positions (octahedral voids of the ccp ) B C possesses three strong 2e2c C–B E bonds, which are additionally strengthened by the 3e3c bonds of the C–B C –C chains, which impose supplementary negative pressure upon the surroundings, as described in . Now it is clear that the similar contraction of the crystals of α–B and B C (ca. 18%) happens in B C on the expense of the icosahedra, which are set into the very rigid surrounding of strong 2e2c and even stronger 3e3c bonds, but in α–B on the expense of weak 2e3c bonds. This conclusion is justified by the behavior of the unit cell parameters ratio c/a in α–B and B C upon compression: in α–B the c/a ratio increases with pressure, but it decreases in B C (Figure 5). In the both cases, above ca. 40 GPa a clear tendency to the leveling is observed that reflects the more homogeneous compression in the both directions at further pressure increase. In the light of the consideration made for α–B and B C above, it is not a surprise now that the volume reduction of B icosahedra in γ–B (16.9%) is intermediate between those in α–B and B C in the same pressure range. As shown in Ref. 28, in γ–B there is a broad diversity of inter-icosahedral bonds and those involving the boron dumbbell (Figure 1a-d, middle panels): three different kinds of strong 2e2c bonds (B3-B3, B2-B5, and B5-B5); for consistency, we follow here the atoms and bonds notations introduced for γ–B in Refs. 9 and 28. Although initially B5-B5 is slightly longer than B3-B3 and B2-B5, it contracts slightly slower than the latter two. The presence of two types (2e3c and 1e2c) of polar-covalent bonds (2e3c B4-B4-B5, involving two boron atoms of one icosahedron and one atom of a dumbbell, and 1e2c B1-B4, involving atoms of neighboring icosahedra at a distance slightly shorter than that between the B E atoms in α–B) makes the rigidity of the inter-icosahedral space in γ–B to be intermediate between those in α–B and B C (see Figure 1). To compare the pressure evolution of the c/a ratio in α–B and B C with the evolution of the corresponding c´/a´ ratio in the structure of γ–B (c´ is in the direction perpendicular to the close-packed layers and a´ - within the layer), we expressed the c´ and a´ 0 through the a, b and c parameters of the orthorhombic unit cell of γ–B and plotted the c´/a´ ratio vs P in Figure 5. As seen, the c´/a´ ratio in γ–B, similarly to what is observed in α–B, first slightly increases, then saturates. Supplementary Figure S1 shows the pressure evolution of the normalized c/a ratios for all of these three solids: the c´/a´(P) curve of γ–B indeed takes an intermediate position between the c/a(P) curves of α–B and B C . Thus, the compressional behavior of γ–B is likely controlled by the balance between the contraction of icosahedra and the inter-icosahedral bonds. Figure 4 confirms this interpretation: the lengths of intra-icosahedral contacts in γ–B at ambient pressure are quite similar to those of α–B and B C (see positions of purple triangles in Figure 4), whereas rates of their contraction under pressure vary considerably (see their scatter within the black rectangular outlining the intra-icosahedral bonds in Figure 4). Some of them are as compliant as inter-icosahedral contacts (B1-B4 and B5-B5), whereas others are almost as stiff as 2e2c bonds inter-icosahedral and dumbbell bonds of γ–B itself and 2e2c bonds of boron carbide. Coming back to the difference in the compressibility of icosahedra in B C and stoichiometric boron carbide B C , one should take into account different chemical composition, i.e. the boron to carbon ratio, of these two carbides. It is established that boron carbide exists as a single-phase material with a wide homogeneity range of carbon content, from ~7 at. % (B C) to ~20 at. % (B
C), realized by the substitution of boron and carbon atoms for one another within both the icosahedra and intericosahedral chains . Proposed stoichiometries are based on two extreme models, B (CBB) on the boron-rich and B C(CCC) on carbon-rich ends. The location of C atoms in the crystal structure, which is a long-standing debate
10, 14-24 , is difficult to clarify even on single-crystal diffraction data because of the similarity in both electronic and nuclear scattering cross-sections for boron and carbon. For their “nearly stoichiometric boron carbide B C” Dera et al. suggested the presence of 85 atomic % boron and 15 atomic % carbon in the B P positions in B C icosahedra. The presence of carbon likely has to lead to an increase of the rigidity of the icosahedron closo -cluster of B C 1 due to the higher electronegativity of carbon compared to boron. Indeed, the volume reduction of the B C icosahedron (13%) appears to be smallest compared to B icosahedra in boron allotropes α–B (14.5%) and γ–B (16.9%) (see Table 2). The stoichiometric boron carbide B C studied in the present work has been proven to contain carbon only in the C-B-C chains , so that its B icosahedra are much more compliant than B C (Table 2), moreover their compliance is similar to that of the whole structure. This observation gives additional evidence that the stoichiometric boron carbide B C is a compound with true covalent bonding: B icosahedra do not play a role of “molecules”, their conventional separation is surely convenient for geometric presentation of the structure, but the compressional behavior of the stoichiometric boron carbide is governed by a complex interplay of both intra-cluster and inter-cluster bonds, as well as those involving C-B-C chains. To conclude, in the present work we have studied the compressional behavior of the stoichiometric boron carbide B C in the pressure interval up to 68 GPa. Our single-crystal synchrotron X-ray diffraction investigations revealed structural stability of the boron carbide in the studied pressure range. A comparison of the unit cell volume reduction with the reduction of the volume of the B icosahedron upon compression of B C from ambient pressure to 60 GPa revealed their similarity. This confirms that the stoichiometric boron carbide B C is a true covalent compound and does show neither ‘molecular-like’ nor ‘inversed molecular-like’ solid behavior upon compression, as previously disputed. Our analysis has shown that, in agreement with the modern understanding of bonding in α–B, γ–B, and B C based on the experimental electron-density studies
13, 27, 28 , the compressional behavior of these boron allotropes and boron carbide depends on the types of bonding involved in the course of compression, so that the ‘effective compressibility’ of B icosahedra may vary in a broad range, from ca. 14% in α–B to ca. 18% in B C , as compared to ca. 18% of compression of the corresponding crystals. 2 Methods summary
Synthesis of crystals
Single crystals of B С were synthesized at high pressures (8.5–9 GPa) and high temperatures (1873–2073 K) using the large-volume-press technique. The stoichiometric composition B C was determined by single-crystal X-ray diffraction in agreement with energy-dispersive X-ray (EDX) analysis (B C) . The presence of impurities could be excluded. Diamond-anvil cell experiments
The BX90-type diamond anvil cells (DAC) made at Bayerisches Geoinstitut (Bayreuth, Germany) and diamonds with the culet diameters of 250 µm were used in high pressure experiments. Rhenium gaskets were squeezed between the diamonds to make an indentation with the thickness of 30 µm. Then in the center of indentations, round holes of 120 µm in diameter were drilled. The B С crystals were placed into these chambers. Sizes of the crystals were about 10 × 10 × 15 µm and 20×15×10 µm . Neon was used as a pressure transmitting medium. Ruby balls used for pressure calibration were placed into the pressure chamber. Single-crystal X-ray diffraction
Crystals with size of about 10 × 10 × 15 µm were selected for measurements in a DAC at ID27 at the European Synchrotron Radiation Facility (ESRF). Diffraction data were collected at 293 K using the Perkin Elmer XRD1621 flat panel detector. The monochromatic radiation had the wavelength of 0.37380 Å and the crystal-to-detector distance was 383 mm. Pressure in the cell was increased up to 68 GPa with a step of about 6 GPa. 160 frames in the omega scanning range of −40° to +40° were collected (0.5° scanning step size) with an exposure time of 1 s. Integration of the reflection intensities and absorption corrections were performed using CrysAlisPro software
30, 31 . The structure was refined in the anisotropic approximation for all atoms by full matrix least-squares using Jana software . 3 References
1. Thévenot, F. Boron carbide — A comprehensive review.
J. Eur. Ceram. Soc. , 205-225 (1990). 2. Ridgway, RR. Boron carbide: A new crystalline abrasive and wear ‐ resisting product. Trans. Electrochem. Soc. , 117-133 (1934). 3. Albert, B, Hillebrecht, H. Boron: Elementary challenge for experimenters and theoreticians. Angew. Chem. Int. Ed. , 8640-8668 (2009). 4. Jemmis, ED, Balakrishnarajan, MM, Pancharatna, PD. A Unifying electron-counting rule for macropolyhedral boranes, metallaboranes, and metallocenes. J. Am. Chem. Soc. , 4313-4323 (2001). 5. Wade, K. The structural significance of the number of skeletal bonding electron-pairs in carboranes, the higher boranes and borane anions, and various transition-metal carbonyl cluster compounds.
J. Chem. Soc. D: Chem. Commun. , 792-793 (1971). 6. Wade, K. Structural and bonding patterns in cluster chemistry. In:
Adv. Inorg. Chem. Radiochem. (Emeléus, HJ, Sharpe, AG). Academic Press (1976). 7. Emin, D. Icosahedral boron ‐ rich solids. Phys, Today , 55-62 (1987). 8. Nelmes, RJ , et al. Observation of inverted-molecular compression in boron-carbide.
Phys. Rev. Lett. , 2268-2271 (1995). 9. Zarechnaya, E , et al. Pressure-induced isostructural phase transformation in γ-B . Phys. Rev. B , 184111 (2010). 10. Dera, P, Manghnani, MH, Hushur, A, Hu, Y, Tkachev, S. New insights into the enigma of boron carbide inverse molecular behavior. J. Solid State Chem. , 85-93 (2014). 11. Hushur, A, Manghnani, MH, Werheit, H, Dera, P, Williams, Q. High-pressure phase transition makes B
C boron carbide a wide-gap semiconductor.
J. Phys.- Condens. Matter , 12 (2016). 12. Fanchini, G, McCauley, JW, Chhowalla, M. Behavior of disordered boron carbide under stress. Phys. Rev. Lett. , 035502 (2006). 13. Mondal, S , et al. Disorder and defects are not intrinsic to boron carbide.
Sci. Rep. , 19330 (2016). 14. Domnich, V, Reynaud, S, Haber, RA, Chhowalla, M. Boron carbide: Structure, properties, and stability under stress. J. Am. Ceram. Soc. , 3605-3628 (2011). 15. Morosin, B, Kwei, GH, Lawson, AC, Aselage, TL, Emin, D. Neutron powder diffraction refinement of boron carbides - nature of intericosahedral chains. J. Alloys Compd. , 121-125 (1995). 16. Emin, D. Structure and single-phase regime of boron carbides.
Phys. Rev. B , 6041-6055 (1988). 17. Kirfel, A, Gupta, A, Will, G. The nature of the chemical bonding in boron carbide, B C . I. Structure refinement. Acta Crystallogr., Sect. B , 1052-1059 (1979). 18. Bouchacourt, M, Thevenot, F. Analytical investigations in the B-C system. J. Less Common Met , 219-226 (1981). 19. Kwei, GH, Morosin, B. Structures of the boron-rich boron carbides from neutron powder diffraction: Implications for the nature of the inter-icosahedral chains. J. Phys. Chem. , 8031-8039 (1996). 20. Lazzari, R, Vast, N, Besson, JM, Baroni, S, Dal Corso, A. Atomic structure and vibrational properties of icosahedral B С boron carbide.
Phys. Rev. Lett. , 3230-3233 (1999). 21. Tallant, DR, Aselage, TL, Campbell, AN, Emin, D. Boron-carbide structure by Raman spectroscopy. Phys. Rev. B , 5649-5656 (1989). 22. Saal, JE, Shang, S, Liu, Z-K. The structural evolution of boron carbide via ab initio calculations. Appl. Phys. Lett. , 231915 (2007). 23. Stein, H, Aselage, TL, Emin, D. Infrared absorption in boron carbides: Dependence on isotopes on carbon concentration. In: AIP Conf. Proc. (Emin, D, Aselage, TL, Switendick, AC, Morosin, B, Beckel, CL). New York (1990). 24. Morosin, B, Aselage, TL, Emin, D. On the crystal structure of boron carbide. In:
AIP Conf. Proc. (Emin, D, Aselage, TL, Switendick, AC, Morosin, B, Beckel, CL). New York (1990). 25. Chuvashova, I , et al.
High-pressure behavior of α-boron studied on single crystals by X-ray diffraction, Raman and IR spectroscopy.
J. Solid State Chem. , 50-60 (2017).
26. Bykova, E, Gou, H, Bykov, M, Hanfland, M, Dubrovinsky, L, Dubrovinskaia, N. Crystal structures and compressibility of novel iron borides Fe B and Fe x B synthesized at high pressure and high temperature. J. Solid State Chem. , 102-109 (2015). 27. Mondal, S , et al.
Experimental evidence of orbital order in α-B and γ-B polymorphs of elemental boron. Phys. Rev. B , 024118 (2013). 28. Mondal, S , et al. Electron-deficient and polycenter bonds in the high-pressure γ-B phase of boron. Phys. Rev. Lett. , 215502 (2011). 29. Kantor, I , et al.
BX90: A new diamond anvil cell design for X-ray diffraction and optical measurements.
Rev. Sci. Instrum. , 6 (2012). 30. Technologies, A, Technologies, A, UK Ltd., O, UK, CrysAlisPro Software system. CrysAlisPro Software system. V.1.171.37.35. (2013). 31. Farrugia, LJ. WinGX and ORTEP for Windows: An update. J. Appl. Crystallogr. , 849-854 (2012). 32. Petricek, V, Dusek, M, Palatinus, L. Crystallographic Computing System JANA2006: General features. Z. Kristallogr. , 345-352 (2014). 33. Belsky, A, Hellenbrandt, M, Karen, VL, Luksch, P. New developments in the Inorganic Crystal Structure Database (ICSD): Accessibility in support of materials research and design.
Acta Crystallogr., Sect. B , 364-369 (2002). Additional information Accession codes:
The X-ray crystallographic coordinates for structures reported in this article have been deposited at the Inorganic Crystal Structure Database (ICSD) Figure 1. Crystal structures of α-B, γ-B and B C . (a.) Rhombohedral cells of α-B and B C compared to the elemental rhombohedron which can be selected in the structure of γ-B. (b.) Single layers of the cubic closest packing ( ccp ) of icosahedra in the structures of α-B, γ-B and B C . The view is perpendicular to the layers: for α-B and B C the direction of the view coincides with the c -axis of the trigonal unit cell (hexagonal settings), for γ-B it coincides with the c´ direction (see text). (c.) The packing of icosahedra shown in projection along the axes of the rhombohedral cells for α-B and B C and along the a -axis for γ-B. (d.) Fragments of the structures of α-B, γ-B and B C showing different types of bonds in the corresponding structures (see text for details). In all structure drawings, atoms in crystallographically independent positions are marked in different colors: for α-B and B C B P are blue, B E are green; for B C carbon and boron atoms of the C-B-C chains are B C (violet) and C (red). For γ-B B1 are blue, B2 are green, B3 are purple, B4 are orange, B5 are violet. Different types of bonds are shown in different colors: 2e2c bonds are yellow, 1e2c bonds are magenta, 2e3c bonds are brown triangles, 3e3c bonds are black. Figure 2. Pressure-dependent evolution of the relative lattice parameters and the relative unit cell and B icosahedron volumes for B C . Continuous lines show the fit of the corresponding data with the third-order Birch-Murnaghan equation of state . Figure 3. Pressure-dependent evolution of bonds of B C . (a.) Intraicosahedral bonds; (b.)
Intericosahedral bonds. Bonds designations are according to Mondal et al. , corresponding notations from Dera et al. (b1 through b7) are given in brackets to make it easier to compare. Figure 4. Relative change of interatomic distances for α-B, γ-B and B C single crystal plotted against their length at lowest pressure as revealed by in situ single-crystal X-ray diffraction. Circles stay for bonds in α-B, triangles for γ-B, and squares for B C . Intraicosahedral bonds are outlined by the black rectangular; red, purple and orange symbols correspond to α-B, γ-B and B C , respectively. Intericosahedral bonds and those involving B-B dumbbells and C-B-C chains are shown in green, blue, and black colors for α-B, γ-B and B C , respectively. The bonds in γ-B are shown in triangle: intraicosahedral ones (within black rectangular) are designated in purple, between icosahedra are in blue. Figure 5. Evolution of the ratios of the unit cell parameters c/a of α-B and B C compared to that of c´/a´ of γ-B. For α-B the data are taken from Chuvashova et al. for γ-B from Zarechnaya et al. Figure 1 a. Figure 1 b. Figure 1 c. Figure 1 d.
Figure 1. Crystal structures of α-B, γ-B and B C . (a.) Rhombohedral cells of α-B and B C compared to the elemental rhombohedron which can be selected in the structure of γ-B. (b.) Single layers of the cubic closest packing ( ccp ) of icosahedra in the structures of α-B, γ-B and B C . The view is perpendicular to the layers: for α-B and B C the direction of the view coincides with the c -axis of the trigonal unit cell (hexagonal settings), for γ-B it coincides with the c´ direction (see text). (c.) The packing of icosahedra shown in projection along the axes of the rhombohedral cells for α-B and B C and along the a -axis for γ-B. (d.) Fragments of the structures of α-B, γ-B and B C showing different types of bonds in the corresponding structures (see text for details). In all structure drawings, atoms in crystallographically independent positions are marked in different colors: for α-B and B C B P are blue, B E are green; for B C carbon and boron atoms of the C-B-C chains are B C (violet) and C (red). For γ-B B1 are blue, B2 are green, B3 are purple, B4 are orange, B5 are violet. Different types of bonds are shown in different colors: 2e2c bonds are yellow, 1e2c bonds are magenta, 2e3c bonds are brown triangles, 3e3c bonds are black. Figure 2 Pressure-dependent evolution of the relative lattice parameters and the relative unit cell and B icosahedron volumes for B C . Continuous lines show the fit of the corresponding data with the third-order Birch-Murnaghan equation of state. Figure 3 a. Figure 3 b. Figure 3. Pressure-dependent evolution of bonds of B C . (a.) Intraicosahedral bonds; (b.)
Intericosahedral bonds. Bonds designations are according to Mondal et al. ; corresponding notations from Dera et al. (b1 through b7) are given in brackets to make it easier to compare. 3 Figure 4. Relative change of interatomic distances for α-B, γ-B and B C single crystal plotted against their length at lowest pressure as revealed by in situ single-crystal X-ray diffraction. Circles stay for bonds in α-B, triangles for γ-B, and squares for B C . Intraicosahedral bonds are outlined by the black rectangular; red, purple and orange symbols correspond to α-B, γ-B and B C , respectively. Intericosahedral bonds and those involving B-B dumbbells and C-B-C chains are shown in green, blue, and black colors for α-B, γ-B and B C , respectively. The bonds in γ-B are shown in triangle: intraicosahedral ones (within black rectangular) are designated in purple, between icosahedra are in blue. Figure 5. Evolution of the ratios of the unit cell parameters c/a of α-B and B C compared to that of c´/a´ of γ-B. For α-B the data are taken from Chuvashova et al. , for γ-B from Zarechnaya et al. Table 1. Crystallographic details of B C at variable pressure. Formula: B C ; Crystal System: Trigonal; Z=3; Space group: R m , crystal size 10 × 10 × 15µm (crystal 1) and 20 × 15 × 10 µm (crystal 2, designated as “cr2”). P, GPa a, Å c, Å V, Å R int * R [I>3σ(I)]* N m ** N i ** 0.00010(1) *R int and R [I > 3 σ(I)] relate to data collection and structure refinement, respectively. **N m is the number of measured unique reflections; N i - the number of unique reflections with I>3σ(I) Table 2. Parameters of the equations of state of B С , B C, α-B, and γ-B obtained on the basis of single-crystal XRD in comparison with those of icosahedra.
Volume reduction of the unit cells and icosahedra on compression was calculated for the pressure range from ambient to 60 GPa. Material Ref. K , GPa K´ K ico , GPa K´ ico
Volume reduction, % Unit cell Icosahedron boron carbide B C Present study 239(7) 3.2(3) 239(23) 3.8(8) 18.7 18.1 B C Dera et al. (2014) , P<40 GPa 227(3) 2.5(2) 285(6) 1.8(3) 18.9 16.9 7 Table 3.
The B icosahedron volume and bond lengths in boron carbide at variable pressure. “cr2” designates crystal 2. B P , B E and B C are designations from Mondal et al. ; b1 through b7 are designations from Dera et al. P, GPa V ico ,Å Bonds involving CBC Inter-cluster bonds Intra-cluster bonds B E -C (b7) B C -C (b1) B P -B P (b2) B P –B P (b6) B P –B E (b4) B P –B E (b5) B E –B E (b3) 0.00010(1) Supplementary Table S1. Fractional atomic coordinates and equivalent isotropic displacement parameters for B C at variable pressure. (Crystal 2 is designated as “cr2”). *Wyckoff positions are (x, 𝒙� , z) for both B P and B E atoms, (0,0, 0) for B C atom and (0, 0, z) for C atom. P, GPa Atom * x z U eq , Å P E C
0 0 0.0091(19) C 0 0.1189(4) 0.0036(11) 10.0(5) B P E C
0 0 0.0123(17) C 0 0.1193(3) 0.0097(11) 20(1) B P E C
0 0 0.0120(15) C 0 0.1202(3) 0.0074(9) 23(1) B P E C
0 0 0.020(3) C 0 0.1206(6) 0.0162(18) 30(1) cr2 B P E C
0 0 0.0071(15) C 0 0.1208(3) 0.0036(9) 35(1) cr2 B P E C
0 0 0.0144(15) C 0 0.1212(3) 0.0088(9) 43(1) cr2 B P E C
0 0 0.0147(15) C 0 0.1217(3) 0.0079(9) 49(1) cr2 B P E C
0 0 0.0153(15) C 0 0.1226(3) 0.0075(9) 56(1) cr2 B P E C
0 0 0.0151(16) C 0 0.1225(3) 0.0060(9) 65(1) B P E C
0 0 0.016(2) C 0 0.1232(4) 0.0069(12) 68(1) cr2 B P E C
0 0 0.0153(18) C 0 0.1233(3) 0.0036(10) 9
Supplementary Figure S1. Evolution of the normalized ratios of the unit cell parameters c/a of α-B and B C compared to that of c´/a´ of γ-B. For α-B the data are taken from Chuvashova et al. ; for γ-B from Zarechnaya et al.9