3.1 Optimized structure
Two newly discovered B2CO structures have been proposed and their crystal structures schematic at atmospheric pressure are shown in Fig. 1. For the first phase with SG Pmn21 is primitive orthorhombic structure containing 8 atoms per unit cell. In order to distinguish from the previously proposed oP8-B2CO, here denoted it as oP8’-B2CO (oP8’ for short). As shown in Fig. 1a, there are two classes of B atoms (B1 and B2) that are not equivalent in oP8’, the Y coordinates of B1 atom are 0.049 or 0.951, which almost located at the cell boundary, while the Y coordinates of B2 atom are 0.380 or 0.620. For B1 atoms, they just bind with three C atoms by sp2 hybridization and form BC3 triangular configuration, however for B2 atoms, they connect with three O atoms and one C atom by sp3 hybridization and form tetrahedron. The C atoms all have 4 coordination with four neighboring atoms as three B1 atoms and one B2 atom, while for O atoms all bonding with three B2 atoms by sp2 hybridization. In other word, there are sp2-sp3 hybridization forms coexist in oP8’, which is quite different from oP8, which is just sp3 hybridization. There are holes in the oP8’ structure that run along the X-axis, result in its low density of 2.813 g/cm3.
For the second B2CO phase, it belongs to the primitive monoclinic with the SG of P21/c, and the unit cell contains 16 atoms, hence denoted as mP16-B2CO (mP16 for short). There are also two kinds of B atoms (B1 and B2) that not equivalent in mP16, the X coordinates of B1 atom is 0.180 or 0.820, and the X coordinates of B2 atom are 0.338 or 0.662. Different from the oP8' structure, all the atoms in the mP16 structure are quadruple coordination relations, thereinto all O atoms are bonded to two B1 atoms and two B2 atoms, all C atoms are also connected with two B1 atoms and two B2 atoms, all of B atoms are bonded to two O atoms and two C atoms. Hence there is no any case about C atom directly connected with O atoms. At the same time, from the Z-axis, the structure is similar to a six-cell honeycomb, while the composition of B-C-O six membered ring is diverse and the six atoms are not in the same plane. Other detail structural information about the two new structures is exhibited in Table 1.
Table 1 The cell parameters a, b, c (Å), ρ (g/cm3) and Atomic Wyckoff positions for two newly B2CO phases.
oP8ʹ
|
mP16 (β=88.6°)
|
ρ (g/cm3)
|
a (Å)
|
b (Å)
|
c (Å)
|
ρ (g/cm3)
|
a (Å)
|
b (Å)
|
c (Å)
|
2.813
|
2.617
|
5.804
|
3.859
|
3.052
|
4.509
|
5.246
|
4.567
|
AWP
|
x
|
y
|
z
|
AWP
|
x
|
y
|
z
|
B1
|
0.5
|
0.049
|
0.118
|
B1
|
0.820
|
0.354
|
0.410
|
B2
|
0.5
|
0.620
|
0.303
|
B2
|
0.662
|
0.625
|
0.915
|
C
|
0.5
|
0.876
|
0.443
|
C
|
0.304
|
0.869
|
0.401
|
O
|
0.5
|
0.437
|
0.598
|
O
|
0.186
|
0.137
|
0.910
|
3.2 Stability analysis
The elastic stabilities of oP8’- and mP16-B2CO at AP are evaluated by the calculation of the independent elastic constants Cijs, as presented in Table 2.
For oP8’, it belongs to orthorhombic system, Born criteria are listed in Eq. (1) [34]:
![](https://myfiles.space/user_files/83062_751fab6dfaef2446/83062_custom_files/img1625047193.png)
For mP16-B2CO with Laue class as 2/m, belong to monoclinic system, it has 13 independent elastic constants. The generic necessary and sufficient criterion that all eigenvalues of elastic matrix be positive is easy to check with simple linear algebra routines. [34]
Table 2 The elastic matrices of two newly B2CO phases at AP, the numerical unit is unified with GPa.
711.6
|
39.2
|
68.3
|
|
|
|
483.5
|
211.8
|
132.1
|
|
-3.0
|
|
……
|
338.2
|
120.4
|
|
|
|
……
|
409.5
|
69.9
|
|
-10.3
|
|
……
|
……
|
195.4
|
|
|
|
……
|
……
|
621.1
|
|
52.1
|
|
|
|
|
177.7
|
|
|
|
|
|
212.7
|
|
1.2
|
|
|
|
|
184.4
|
|
……
|
……
|
……
|
|
222.1
|
|
|
|
|
|
|
147.1
|
|
|
|
……
|
|
188.4
|
oP8’-B2CO
|
mP16-B2CO
|
The elastic matrices are symmetric about the main diagonal, and these symmetric nonzero data are represented by apostrophe symbol “……”, the blank cell represents data that is 0. There is no doubt that from the above elastic matrix in Table 2, the Cijs satisfy the criteria above, indicating oP8’-B2CO and mP16-B2CO are elastic stable at AP.
The imaginary frequency will lead to crystal distort, indicates the dynamical instability. Phonon dispersion spectrum and the related phonon density of states (DOS) of two newly predicted B2CO phases at AP are researched and plotted in Fig. 2. There are no negative phonon modes in entire Brillouin zone of their unit cells, suggesting their dynamically stable.
The stability of oP8’ and mP16 B2CO were also intensively researched at high pressure. Here, taking 10 GPa, 50 GPa and 100 GPa as examples, the elastic constant and phonon scattering of oP8’ and mP16 B2CO under high pressure are studied. As displayed in Table S1, the elastic constant Cijs of oP8’ and mP16 B2CO under high pressure all satisfied the Born criterion of elastic stability. And the phonon scattering spectra and their density of states in Fig. S1 all indicate that there are no virtual frequency. That is to say, in the range of pressures that we studied, both oP8’ and mP16 B2CO have high pressure stability.
3.3 Potential synthesis reaction
For providing guidance of experimental synthesis, it’s of great significance and extremely urgent to deploying the potential synthesis of novel IED B-C-O phases. Based on the common reactants such as boron trioxide B2O3, boron-rich oxides B6O, graphite C, and α-boron B, construct two different reaction paths, then balanced the synthesis equations of these novel IED B2CxO phases, and therewith the relationships between formation enthalpy ΔHf and pressure re revealed according the following Eq. (2)~(3). The ΔHf1 and ΔHf2 represent the formation enthalpy via path 1 and path 2, respectively.
![](https://myfiles.space/user_files/83062_751fab6dfaef2446/83062_custom_files/img1625047158.png)
Fig. 3 summarizes the relationships between formation enthalpy (normalized on a per formula) and pressure for all IED B-C-O compounds. The enthalpy of formation is negative, indicating the likelihood that the target IED B-C-O compounds will be synthesized by the selected reaction, and more negative the formation enthalpy, more reaction driven energy, more possible to be synthesized. Hence the critical pressure for formation enthalpy reaching zero is the key to determine the HP synthesis.
As shown in Fig. 3a, for B2CO constituted by sp3 hybridization, no matter what kind of synthesis path, the critical pressure of formation enthalpy from positive to negative is the minimum for tI16 and the maximum for mP16. And the critical pressures of sp3 hybridization B2CO fall into two categories: mP16 and oC16 with high critical pressures, others with low critical pressures. For B2CO contained sp3 and sp2 hybridization displaying in Fig. 3b, oI16 is not a pressure driven structure due to the formation enthalpy increases gradually with the increase of pressure and changes from negative to positive, hence it is unlikely to be synthesized by high pressure of common reactants such as B2O3, B6O, C and B. Compared with tP16, oP8’ has significantly smaller critical pressure, which indicates that oP8’ has the advantage of synthetic conditions. As exhibited in Fig. 3c, for B2CXO with similar structure type, with the increase of carbon content, the critical pressure gradually decreases in either synthesis path. In other words, the increase of carbon is conducive to the synthesis of B2CXO (X=2, 3, 5).
Except that oI16 is not a pressure-driven structure, the rest are all pressure-driven structures, and the critical pressure is within 100 GPa, which can be synthesized by high pressure experiment with diamond anvil cell (DAC). Some of these IED B-C-O phases have critical pressures of less than 30 GPa and may be synthesized by high pressure experiments in LVP. Hence, the formation enthalpies as a function of pressure at 0 to 30 GPa with sampling intervals of 2 GPa has been studied detailed. As shown in Fig. S2a, for the path 1 of B2CO, their critical pressures are between 12~20 GPa, and following the sequence as tI16 < tP4 < oP16 < oP8 < oP8’. The path 2 critical pressures of tI16 and tP4 are in the range of 27~30 GPa, which has obvious disadvantages in terms of synthetic pressure than path 1. As exhibited in Fig. S2b, for path 1 of B2CXO (X=2, 3, 5), their critical pressures are 29.2 GPa, 20.0 GPa and 14.9 GPa, respectively. In other words, with the increase of carbon content of B2CXO (X=2, 3, 5), the critical pressure of formation enthalpy gradually decreases from positive to negative. For the path 2 of B2C3O and B2C5O, their critical pressures are 27.4 GPa and 19.3 GPa, both are higher than that of path1. As the working pressure and temperature of LVP can reach tens of GPa and thousands of degrees centigrade, and the high temperature will benefit to reduce the reaction barrier, hence synergistic high temperature and HP experiments in LVP maybe a potentially efficient synthesis technology for the target B-C-O product at a lower pressure than the theoretical critical pressure value at zero temperature. Of course, special technologies need to be combined, such as rapid pressure relief, rapid cooling and so on, which will help to preserve these IED B-C-O phases.
3.4 High pressure deformation
It is well known that pressure is an important parameter to regulate the state of matter. In general, solid materials undergo a continuous densification during HP compression. As exhibited in Fig. 4, for these B2CO phases we studied, during the pressure from AP to HP as 100 GPa, the volumes per molecular formula of these phases decrease monotonically. For the six pure sp3 hybridization connected phases, their volume compressions are basically similar, the largest volume compression appears on mP16, while the minimum volume compression appears on tP4, the difference between the two is very small. For three sp2 and sp3 hybridization coexist phases, all have large pressure induced volume change. Detailed research has revealed that the B2CO phases with sp3 hybrid connection have lower volume compression ratio than that of sp2 and sp3 hybrid coexist. This is due to the existence of holes in the sp2 and sp3 hybridization coexisting phases, which results in their weaker resistance to deformation than the dense phases formed by pure sp3 hybridization.
In addition, the volume per molecular formula of B2CXO (X=2, 3, 5) as a function of pressure was studied in detail, and shown in Fig. S3a. For B2CXO (X=2, 3, 5), their volumes vary with pressure between AP to 100 GPa high pressure as follows: B2C2O 6.388 Å3/20.034 %; B2C3O 7.154 Å3/19.144 %; B2C5O 8.855 Å3/18.213 %. That is to say, with the increase of carbon content, their volumes shrinkage gradually increase, while under the effect of volume size at AP, their volume compression rate decrease.
In fact, the phonon vibration frequency is affected by the strength of the chemical bond. The max phonon vibration frequency is an indicator of the strength of the chemical bond. Based on the phonon scattering of oP8’ and mP16 B2CO under different pressures, it is found that their max phonon vibration frequencies are 40.162 THz/AP, 41.547 THz/10 GPa, 45.483 THz/50 GPa, 47.512 THz/100 GPa for oP8’ and 39.411 THz/AP, 40.754 THz/10 GPa, 44.691 THz/50 GPa, 48.765 THz/100 GPa for mP16. With the increase of pressure, the maximum vibration frequency increases, indicating that pressure will induce the increase of bond energy, which may be due to the decrease of atomic spacing under the action of pressure. Furthermore, zero point vibrational energy (ZPVE) is the energy that a quantum will vibrate under the zero point of absolute temperature. The amplitude of zero point vibration increases with increasing temperature, and the lighter the atom (that is, the smaller the element number), the more pronounced the zero point vibration. Based on the phonon scattering spectrum and its DOS, ZPVE can be calculated based on
is the phonon density of states,
is Planck constant. ZPVE per B2CO chemical formula of oP8’ are 0.541 eV/AP, 0.560 eV/10 GPa, 0.614 eV/50 GPa, 0.666 eV/100 GPa. And ZPVE (eV) per formula of mP16 are 0.522 /AP, 0.547 /10 GPa, 0.615 /50 GPa, 0.676 /100 GPa. The researched relationship of ZPVE and pressure has revealed that there is also a positive correlation between pressure and ZPVE of oP8’ and mP16.
3.5 Electrical property
Based on GGA-PBE functional, the electronic band structures along the symmetry points of the Brillouin zone of all known B2CO phases at AP are investigated and exhibited in Fig. 5, thereinto a~f are the B2CO phases formed by sp3 hybridized, and g~i are the B2CO phases constituted by sp2-sp3 hybridized.
For tP4 and tI16, there is no band crossing Fermi level, and the forbidden bands are determined by valence band max VBM (trefoil) and conduction band min CBM (heart). As the VBM and CBM locates at different symmetry points, both tP4 and tI16 are indirect bandgap semiconductors with gaps as 1.666 eV and 2.971 eV, respectively. For three orthorhombic phases oP8, oP16 and oC16, they all have semiconductor properties with the indirect band gap as 3.522 eV, 3.244 eV and 3.461 eV, respectively. For sp2-sp3 hybridized phases oI16 and tP16, they both have semiconductor properties with an indirect band gap as 3.658 eV and 1.880 eV, respectively. Among all the B2CO, mP16 has the widest band gap as 4.681 eV, oP8’ has the narrowest band gap as 1.607 eV. The proposed two phases greatly expand the band-gap properties of B2CO as semiconductor materials and enrich the potential applications of B2CO in the electronic industry field.
Here the electrical properties of the carbon-rich type B2CXO (X=2, 3, 5) at atmospheric pressure were revealed based on PBE. As displayed in Fig. 6, as the VBM (trefoil) and CBM (heart) of B2C2O and B2C3O locates at the same symmetry points G, both B2C2O and B2C3O are direct semiconductors, which is consistent with the conclusions of Zhang et al [18]. As for B2C5O, the VBM and CBM locates at symmetry points G and M, indicating that B2C5O is an indirect semiconductor, which is also in well agreement with previous studies [18]. The band gap values of B2CXO (X=2, 3, 5) are 1.888 eV, 2.485 eV and 2.280 eV, respectively.
The relationship of band gap values and pressure of these IED B-C-O phases are researched. As shown in Fig. 7, there are two variation trends of band gaps with pressure. For these B2CO phases with high symmetry and composed of pure sp3 hybridization, such as tP4, tI16, oP8, oP16 and oC16, their band gaps all increases monotonically as the pressure increases. In the pressure range of 100 GPa, their band gaps increased 0.166 eV/9.96%, 0.187 eV/6.29%, 0.639 eV/18.14%, 0.252 eV/7.77%, and 0.515 eV/15.04%, respectively. As for these B2CO phases composed by sp3-sp2 hybridization coexist (oI16, oP8’ and tP16) or with low symmetry like mP16, their band gaps increased first and then decreased with the increase of pressure. In different phases, the inflection point of rising first and falling second is different. For mP16 and oP8’, the inflection point is at 10 GPa, and their max band gaps are 4.706 eV and 1.643 eV, respectively. As for oI16 and tP16, their inflection points are where the pressure is higher, and the max gap values are 4.807 eV and 2.662 eV, respectively.
Here the relationship between the band gap of the carbon-rich system B2CXO (X=2, 3, 5) and the pressure has been investigated. As displayed in Fig. S3, in the pressure range of ambient pressure up to high pressure as 100 GPa, the band gap value of B2C3O undergo near invariance and then slow decline, among the band gap between ambient pressure and 30 GPa varies by only several meV, then the gap drops about 30 meV within 30 GPa to 100 GPa. For B2C2O, its gap increases at first and then decrease, the max gap value 2.416eV appears in 70~80 GPa. The gap of B2C2O goes through sustained decline with pressure increasing, and the 100 GPa pressure result in the gap decreasing amplitude as 0.323 eV/17.11%.
3.6 Mechanical property
Here the mechanical properties as bulk modulus B and shear modulus G of two newly predicted B2CO phases and their relationship with pressure are detailed studied. The independent elastic constant Cijs of oP8’ at different pressure are listed in Table S2, and the elastic matrix of mP16 are displayed in Table S3. Because mP16 is composed of pure sp3 hybridization, which is denser than oP8' formed by the coexistence of sp2 and sp3 hybridization. Generally, B indicates the ability to resist volume deformation by loading pressure,[35] and can be calculated by independent elastic parameter. [36] Hence mP16 has larger elasticity moduli than oP8’ across the whole studied pressure range, as shown in Fig. 8a. The pressure has a positive correlation to the bulk modulus B of oP8’ and mP16, therein the linear fitting shows that the bulk modulus B of mP16 has a perfect linear relationship with pressure P (AP to 100 GPa) as B=268.3+3.323P with R2=0.998.
Generally, shear modulus G can be acquired by independent elastic parameter, [36] represents the ability to resist deformation upon shear stress. [35] Due to the structural differences, the shear modulus G of oP8’ is smaller than that of mP16 at AP. As exhibited in Fig. 8b, the shear modulus G of oP8’ and mP16 is also positively affected by pressure. When the pressure is above 80 GPa, the shear modulus G of oP8’ will be larger than that of mP16. The data fitting results show that the shear modulus G of oP8’ has a perfect linear relationship with pressure P (AP to 100 GPa) as G=144.9+2.073P with R2=0.993, and G of mP16 has the relationship as G=193.3+3.281P-0.037P2+1.715×10-4P3, R2=0.996.
The Young modulus E represents the rigidity of the condensed material, the greater E is, the less likely to deform. Based on the relationship for B, G and E as E=9BG/(3B+G), [36] the Young modulus E of oP8’ and mP16 are also detailed studied. As displayed in Fig. 8c, mP16 has higher E than oP8' at AP. However, with the increase of pressure, the Young's modulus E of oP8’ increases rapidly, which is nearly the same as that of mP16 at 80 GPa, and then exceeds that of mP16 with the increase of pressure. The positive correlation between pressure P (AP to 100 GPa) and Young's modulus E can be expressed as mP16: E=464.2+7.839P-0.078P2+3.481×10-4P3 with R2=0.998, and oP8’: E=337.4+5.271P with R2=0.992.
As reflect the brittleness or toughness of materials, Poisson's ratio μ is also an important mechanical property of materials and can be calculated based on μ=0.5(3B-2G)/(3B+G).[36] We also studied in detail the effect of pressure on the Poisson's ratio μ of oP8’ and mP16. As demonstrated in Fig. 8d, they both have value of μ smaller than 0.333, which is the critical value for a criterion of ductility or brittleness. The results of Poisson's ratio show that both of them are brittle materials. Ignored the first point at AP, the pressure P (10 GPa~100 GPa) to Poisson's ratio of mP16 satisfies the following relationship as: μ=0.190+8.549×10-4P with R2=0.997.
The hardness is extensively adopted to evaluate the mechanical properties of condensed matter. Here, the Vickers hardness (Hv) of two newly discovered B2CO phases were calculated based on Eq. (4). [1] The directional dependence of material’s elastic property can be quantified as anisotropy. The mechanical behavior of materials such as plastic deformation, elastic instability and crack behavior are dominantly affected by the elastic anisotropy. Hence, it's essential for us to systematically investigate their elastic anisotropy for potential industrial applications. The universal elastic anisotropy index Au [37] is employed to evaluate the integrated anisotropy as Eq. (5).
![](https://myfiles.space/user_files/83062_751fab6dfaef2446/83062_custom_files/img1625046884.png)
We also studied in detail the effect of pressure on their hardness and elastic anisotropy index Au. As displayed in Fig. 8e, both oP8’ and mP16 have high hardness than 20 GPa in all the pressure range we studied, suggesting that oP8’ and mP16 are hard materials. The hardness of mP16 increases during the pressure from AP reaches to 20 GPa, and then the hardness continues to decrease in the process of increasing the pressure to 100 GPa. The relationship between hardness Hv and pressure P (30 GPa ~ 100 GPa) in the descending stage is satisfied as Hv=34.2-0.0756P with R2=0.998. For oP8’, its hardness starts with a small increase and then a rapid increase and ends with a near constant change during the process from chamber pressure to 100 GPa. The value of Au for isotropic crystal is 0. Any departure from 0 indicates the anisotropic degree which considers the contributions of both the bulk and shear moduli. As exhibited in Fig. 8f, with the increase of pressure, the anisotropy of oP8’ experienced a process of first decreasing, then nearly unchanged, and then decreasing again. For mP16, it has the degree of anisotropy lower than oP8’. The elastic anisotropy index Au of mP16 first decrease and then increase with pressure increase from AP to high pressure as 100 GPa. In the pressure range from 20 to 100 GPa, the Au of mP16 has a relationship with pressure P as Au=0.071+2.79×10-3P with R2=0.998.
3.7 Simulated tensile deformation
The material’s ideal strength, which is defined as the critical stress when a perfect crystal lost its elastic stability, equal to the upper limit for material strength. The atomistic mechanism for structural deformation and failure models can be clearly interpretated with the investigation of strain-stress relations and bond-breaking processes. Here the tensile deformation of oP8’ with high symmetry has been studied in detail.
Fig. 9a presents the simulated tensile stress–strain relationships along the specific directions including [100], [010] and [001] for oP8’. It can be seen that tI16-B8C6O2 has the max tensile stress of 64.7 GPa in the [100] direction with the strain at 0.21, with the strain exceed 0.24 in [100] direction, oP8’ will failure. For [010] and [001] direction, overall tensile stress increases synchronously with the strain, the max tensile stress and max strain appear simultaneously. Once the stress or strain exceeds the maximum ([010]: about 37.2 GPa /0.16; [001]: about 12.0 GPa /0.13), oP8’ will be destroyed. As shown in Fig. 9b, with the appearance of structural strains in different directions, the band gap of oP8’ will mutate, then as the strain of [010] direction increases, the band gap continues to decrease. As the strain of [100] direction increases, the band gap continues to increase to 2.455 eV until the strain reach 0.17, and then decreases with strain increase. For [001] direction, the gap decreases during the strain up to 0.08 and then increase with strain increase.