Patent Publication Number: US-2023147042-A1

Title: Carbon-based clathrate compounds

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
     This application is a continuation of U.S. Application No. 16/810,388 filed Mar. 5, 2020, which further claims priority to U.S. Provisional Application No. 62/814,024 filed Mar. 5, 2019, both filed in the United States Patent &amp; Trademark Office, the disclosures of which are incorporated herein by reference in their entireties. 
    
    
     STATEMENT OF INTEREST 
     This invention was made with Government support under Grant W31 P4Q-13-1-0005 awarded by DARPA and Grant DE-SC0001057 from the Department of Energy. The U.S. Government has certain rights in the invention. 
    
    
     FIELD OF THE INVENTION 
     The present invention relates to carbon-based clathrate compounds such as SrB 3 C 3  and LaB 3 C 3  as well as methods of making the same. 
     BACKGROUND OF THE INVENTION 
     As a fundamental building block of nature, carbon is unrivalled in its diversity to form stable structures with other elements and itself. One-dimensional (1D) carbon-based materials (e.g., polymers) have thoroughly reshaped society over the past century, in addition to providing the building blocks for life. In recent years, two dimensional (2D) materials, such as graphene, have attracted much attention due to remarkable properties that promise to advance technology. Three dimensional (3D) sp 3  carbon-based structures, such as diamond, exhibit many superlative properties including hardness, strength, thermal conductivity and electron mobility. But aside from diamond, only a handful of materials in this class are actually known. Lonsdaleite (hexagonal diamond), B-doped diamond, SiC and BC 2 N, are valuable high-tech ceramics, but they all maintain the basic structure of diamond. Boron carbide also contains sp 3  hybridized carbon atoms, but these atoms serve as links between B 12  icosahedra, rather than establishing the overall structural framework. From the organic perspective, 3D covalent organic frameworks (COFs), which are formed by linking sp 2 -hybridized molecular building blocks, have attracted much interest as materials for gas storage and separations. Compared with the exquisite synthetic control over porous COF materials, the experimental progress in denser sp 3  carbon-based structures lags far behind. 
     Numerous 3D carbon allotropes and compounds have been predicted to have feasible energies and interesting properties. However, it remains unclear whether any of these materials can be produced in the laboratory. Aside from the diamond structure, almost no other sp 3  carbon-based frameworks are known or can be stabilized at atmospheric pressure. For example, longstanding predictions of 3D sp 3 -bonded C 3 N 4  networks have not been realized thus far, and high-pressure polymeric phases of CO 2 , which consist of a network of CO 4  tetrahedra, decompose into molecular CO 2  phases when decompressed. 
     There is thus a need to predict and synthesize new 3D sp 3  carbon-based structures that have useful properties. Accordingly, the principal object of the invention is to provide carbon-based structures, principally in the form of a carbon-based clathrate compound. Other objects will also be apparent from the detailed description of the invention. 
     SUMMARY OF THE INVENTION 
     Broadly stated, the objects of the invention are realized, accordingly to one aspect of the invention, through the prediction and synthesis of carbon-based clathrate compounds that include a clathrate lattice, guest atoms, and, optionally, substitution atoms for the atoms that form the clathrate lattice. More specifically, the carbon-based clathrate compounds may include (i) a clathrate lattice with atoms of at least one element selected from the group consisting of carbon and boron as a host cage structure; (ii) guest atoms encapsulated within the clathrate lattice; and (iii) substitution atoms that may be substituted for at least one portion of the carbon and boron atoms that constitute the lattice. 
     According to one aspect of the invention, the guest atoms are lanthanum (La). In another aspect of the invention, the guest atoms may also include Ba, Ca, Sr, or other atoms with a similar ionic radius. 
     According to one aspect of the invention, the substitution atoms comprise nitrogen. In another aspect of the invention, the nitrogen atoms are substituted for at least some boron atoms. 
     According to one aspect of the invention, the structure of the carbon-based clathrate compound is cubic bipartite sodalite. 
     According to one aspect of the invention, the clathrate lattice includes cages, each cage including 24 atoms with six four-sided faces and eight six-sided faces. 
     According to another aspect, the invention is a compound of the formula LaB 3 C 3 . 
     According to another aspect of the invention, the LaB 3 C 3  compound is a carbon-based clathrate. In another embodiment, the structure of the carbon-based clathrate is cubic bipartite sodalite. In another embodiment of the invention, the compound comprises a clathrate lattice with atoms selected from the group consisting of carbon and boron as a host cage structure. In another embodiment of the invention, the clathrate lattice is formed of sp 3  hybridized carbon and boron. 
     According to one aspect of the invention, the Lanthanum is a guest atom encapsulated within the clathrate lattice. According to another aspect of the invention, the clathrate lattice comprises cases, each cage comprising 24 atoms with six four-sided faces and eight six-sided faces. 
     Advantageously, the compound of the formula LaB 3 C 3  may have semiconductor properties. Additionally, other carbon-based clathrate compounds such as SrB 3 C 3  may be superconductors. 
     According to one aspect of the invention, the compound of the formula LaB 3 C 3  has a bulk modulus of 255 GPa. In another embodiment of the invention, the compound of the formula LaB 3 C 3  has a has a calculated hardness of about 30 GPa. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee. 
       The invention is more fully described by reference to the following detailed description and the accompanying drawings wherein: 
         FIG.  1 A  is a ternary phase diagram of a Sr-B-C system at 50 GPa. Circles (green) represent thermodynamically stable compounds while squares (orange) represent additional metastable compositions used in the search. 
         FIG.  1 B  is a ternary convex hull for the Sr-B-C system at 50 GPa based on formation enthalpies. Compounds with enthalpy data represented by points (red) are on the convex hull and thermodynamically stable against decomposition. Black points (the points not on the hull) show the formation enthalpies of metastable structures found in the structure searches. 
         FIG.  2    provides the structure of an SrB 3 C 3  clathrate. The structure also may represent the LaB 3 C 3  Clathrate, with La as the guest atom. The cubic structure is composed of face-sharing boron-carbon cages that encapsulate Sr or La atoms. Each cage contains 24 atoms with six four-sided faces and eight six-sided faces (4 6 6 8 ). Different color cages are used to emphasize the stacking of cages that tile 3D space. 
         FIG.  3 A  provides experimental x-ray diffraction data collected at 57(3) GPa with Rietveld refinement of the SrB 3 C 3  phase. More specifically, the experimental X-ray diffraction data (black points) was collected at 57(3) GPa with Rietveld refinement (blue line) of the SrB 3 C 3  phase. Green ticks indicate contributions from Ne with Le Bail refinement. The 2D diffraction cake image aligned with the integrated pattern shows nearly complete powder averaging with sharp peaks for the SrB 3 C 3  phase. The inset shows a magnified view at high angle with sharp SrB 3 C 3  peaks to a limiting resolution of 0.75 Å. 
         FIG.  3 B  provides experimental EoS and calculated EoS data for the SrB 3 C 3  phase. Experimental EoS (solid blue line) with B0 = 249(3) GPa, B0′=4.0 (fixed) and calculated EoS (dashed lines) with B0 (DFT-LDA) = 257 GPa, B0′=4.0 (fixed); B0 (DFT-GGA) = 225 GPa, B0′=4.0 (fixed). Different colored symbols represent data points from six independent experimental runs. 
         FIG.  4 A  provides electronic properties of SrB 3 C 3  and in particular a two-dimensional electron localized function (ELF) for the compound. 
         FIG.  4 B  provides electronic band structures for SrB 3 C 3  at 0 GPa projected onto atomic orbitals, where the width of each band is proportional to the weight of the corresponding orbital character. The projected density of states and Crystal Orbital Hamilton Population (COHP) between adjacent B and C atoms in SrB 3 C 3  at 0 GPa are shown to the right. The Fermi energy is set to 0 eV (dashed line). 
         FIG.  5    provides an image showing in situ synthesis of SrB 3 C 3  with laser heating in a diamond anvil cell. 
         FIG.  6 A  shows phonon dispersion, phonon density of states and electron-phonon integral of SrB 3 C 3  at 0 GPa. 
         FIG.  6 B  shows Tc as a function of pressure for SrB 3 C 3 . 
         FIGS.  7 A- 7 D  show calculated stress-strain relations of SrB 3 C 3 .  FIG.  7 A  and  FIG.  7 B  show calculated tensile stress-strain relations for SrB 3 C 3  with GGA ( FIG.  7 A ) and LDA ( FIG.  7 B ).  FIG.  7 C  and  FIG.  7 D  show calculated shear stress-strain relations for SrB 3 C 3  in the (110) easy cleavage plane with GGA ( FIG.  7 C ) and LDA ( FIG.  7 D ). 
         FIG.  8    provides an XRD pattern showing LaB 3 C 3  recovered to low pressure. The Powder XRD pattern of LaB 3 C 3  is shown during decompression at 5 GPa (top), and after opening the cell in air (bottom). Bragg positions from cubic LaB 3 C 3  are marked by blue lines. The diffraction intensity was limited by the opening angle above 20=14°. When exposed to air, LaB 3 C 3  degrades into unknown, presumably, hydrolysis products. Diffraction from unreacted LaB 3 C 3  can still be observed in air several minutes after the cell was opened. 
         FIG.  9 A  shows computed stability of LaB 3 C 3  via different reaction pathways indicating exothermic formation from pure elements and from stable binaries above ca. 20 GPa. 
         FIG.  9 B  shows Rietveld refinement of LaB 3 C 3  at 57(2) GPa with a = 4.66307(8) Å, wRp-Bknd = 0.042, χ 2  = 3.8. The sample is nearly phase pure with a trace impurity of unreacted LaB 6 . Inset shows well-resolved peaks to a limiting resolution of 0.80 Å. The 2D cake image, presented above the integrated pattern, indicates uniform powder averaging. The pattern obtained before heating (grey trace) shows only cubic LaB 6  and Ne; ballmilled LaC 2  and glassy carbon appear amorphous. 
         FIG.  9 C  provides unit cell volume of LaB 3 C 3  as a function of pressure. Experimental data were fitted using a third-order Birch-Murnaghan EoS with B 0  = 254(9) GPa and B 0 ’ = 3.6(4). DFT results yield B 0  = 248 GPa with B 0 ’ = 3.8 (PBEGW) and B 0  = 268 GPa with B 0 ’ = 3.9 (LDA). 
         FIG.  10    shows phonon dispersion for LaB 3 C 3  at 1 atm showing dynamic stability. 
         FIG.  11    provides stress-strain calculations for LaB 3 C 3 with GGA and LDA., indicating a high ideal strength with hardness near 30 GPa. 
         FIG.  12    provides raw XRD patterns of SrB 3 C 3  before and after heating. 
         FIG.  13    provides raw XRD patterns of SrB 3 C 3  with increasing pressure and heating. 
         FIG.  14    provides Sr-B-C phase identification at 155(6) GPa with Ne medium. 
         FIG.  15    provides Sr-B-C phase identification at 142(10) GPa with AI 2 O 3  medium. 
         FIG.  16    presents data showing the stability of SrB 3 C 3  and analysis of possible stoichiometries. 
         FIG.  17    shows X-ray diffraction of SrB 3 C 3  at atmospheric pressure. 
         FIG.  18    provides electronic band structure for SrB 3 C 3  at 200 GPa projected onto atomic orbitals. 
         FIG.  19 A  and  FIG.  19 B  provide energetic stabilities as a function of pressure, with  FIG.  19 A  showing calculated enthalpies per SrC 2  unit as a function of pressure for various structures with respect to 14 /mmm structure, and  FIG.  19 B  showing calculated enthalpies per SrB 6  unit as a function of pressure for various structures with respect to Pm 3 m structure. 
         FIG.  20    provides the structure of LaB 3 C 3  viewed along different crystallographic directions. C and B atoms make up the framework, while La is trapped in the cages. The single 4 6 6 8  cage has one unique B-C distance of 1.74 Å at atmospheric pressure. 
         FIG.  21 A  shows the difference in charge density for LaB 3 C 3  (crystal density minus isolated atomic density) viewed in the (100) plane at 0 GPa. 
         FIG.  21 B  shows the electronic band structure (HSE) of LaB 3 C 3  at 0 GPa with density of states. Contributions from different orbitals are indicated. 
         FIG.  22    provides an EDX map of recovered sample showing showing homogeneous elemental distribution in heated region (25×25 µm 2 ) with average composition La 1.00±0.07 B 2.95±0.06 C 4.20±0.09 O 0.46±0.07 . The carbon content is biased by adventitious sources, also observed on pure standards. The color composite map is layered over the SEM image and individual elements from the same region are shown below. 
     
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     Carbon clathrate is an impressive 3D sp 3  material. Carbon-based clathrates are open-framework structures composed of host cages that trap guest atoms in which all host atoms are linked by four-coordinate bonds. As sp 3 -bonded frameworks, carbon-based clathrates represent strong and lightweight materials that also offer tunable properties through manipulation of the occupancy and type of guest atoms within the cages. Despite their prominence in other systems with tetrahedral coordination, carbon-based clathrates have not yet been reported due to tremendous challenges associated with their synthesis. 
     Attempts to synthesize carbon clathrates go back at least 50 years since they were postulated following the formation of inorganic silicon clathrates, and their possible structures and properties are of longstanding interest. However, carbon clathrates have not been successively synthesized yet. Some proposed but unrealized carbon clathrates are expected to exhibit exceptional mechanical properties with tensile and shear strengths exceeding diamond, while large electron-phonon coupling is predicted to give rise to conventional superconductivity with high transition temperatures. If produced, these materials may represent a class of diamond-like compounds wherein the electronic structure is tunable by adjusting the occupancy of electron-donating (or withdrawing) atoms within the cages. 
     The inventors have performed a substantial amount of research to answer the persisting question of whether carbon clathrate structures are accessible by experiment. First-principles DFT calculations indicate that both filled and guest-free carbon clathrates are energetically unfavorable but by energies as low as 0.07 eV/atom relative to diamond (for reference, commercially produced C 60  is metastable by nearly six times that energy). Synthesis of carbon clathrates might therefore proceed through a non-equilibrium pathway (e.g., formation from a high-energy precursor or deposition method) or through a chemical substitution/doping strategy to modify the intrinsic thermodynamic stability. No successful metastable pathways to carbon clathrates have been established yet, although three-dimensional polymers of C 60  have been suggested to resemble carbon clathrate-like structures. 
     While non-equilibrium synthesis pathways remain feasible in concept, another strategy is to substitute boron for carbon atoms within the cage frameworks of carbon clathrates. The electron deficient nature of boron creates the ability to form complex chemical bonding with itself or carbon to stabilize polyhedra, such as the icosahedral units in molecular carborane clusters. Zeng et al. calculated that boron substitution can improve the intrinsic thermodynamic stability of carbon clathrate frameworks. Nevertheless, no thermodynamically stable carbon-clathrate was predicted after examination of a small subset of possible B substitution schemes in Li-filled carbon clathrates. A broad search of potential B substitution schemes was therefore needed to validate this chemical stabilization principle. Here, the inventors established the first thermodynamically stable carbon-based clathrate using automatic structure searching methods and validated the prediction via high-pressure and high-temperature experiments, thus expanding known sp 3  carbon materials to a new class with tunable properties. 
     The inventors conducted an extensive computational search in the Sr-B-C system (including pure elements, binaries and ternaries from SrB x C y  with 0 ≤ x, y ≤ 6) at pressures from 0-200 GPa after broader searching in ternary B-C systems with a variety of metals including Li, Na, Mg and Ca. For the Sr-B-C system, several high-pressure compounds were determined to be thermodynamically stable with respect to elemental mixtures (see  FIG.  1 A  and  FIG.  1 B ). At 50 GPa, the hexagonal P6 3 /mmc and y-B structures were calculated to be the most stable forms of Sr and B, respectively, while diamond is the most stable structure for C. The compounds Sr 5 C 2 , SrC, SrC 2 , Sr 5 C 2 , SrB, SrB 2 , SrB 4 , and SrB 6  were found to be the stable binaries on the convex hull at 50 GPa (discussed below in more detail). The inventors found no energetically stable B-C binary compounds above 50 GPa, in agreement with a previous computational study. 
     Two stable ternary compounds were predicted at 50 GPa. The first, hexagonal SrBC (space group P6 3/ mmc), exhibits two-dimensional layers of six-membered B-C rings stacked between layers of Sr atoms, similar to intercalated graphite. This new SrBC phase is isostructural with LiBC found at ambient pressure. The second ternary compound is cubic (space group Pm 3 n)) with composition 2Sr@B 6 C 6  (SrB 3 C 3 ) and has the type-VII clathrate structure, known for the clathrate hydrate HPF 6· 6H 2 O and related to inorganic compounds such as BaPd 2 P 4  through a small structural distortion. The topology of SrB 3 C 3  is that of bipartite sodalite (sod-b), which is distinguished from the sodalite structure (sod) in that carbon atoms are only bonded to boron atoms and vice versa. This is the first prediction of a thermodynamically stable carbon-based clathrate. 
     The SrB 3 C 3  clathrate framework ( FIG.  2   ) is composed of a single truncated octahedral cage with six four-sided faces and eight six-sided faces (4 6 6 8 ). The cages are comprised of 24 vertices with alternating C and B atoms and each cage contains a single Sr atom at the center. This type-VII clathrate phase is predicted to be thermodynamically stable from 50 to at least 200 GPa. The material does not exhibit imaginary phonon frequencies at any pressure indicating dynamic stability and favorable conditions for metastable recovery to ambient conditions. At zero pressure, the optimized lattice parameter is 4.88 Å, and the structure contains one unique B-C bond length of 1.73 Å. The boron-doped clathrate is much more stable than its pure carbon counterpart: at 50 GPa, 2Sr@C 12  is metastable by 0.667 eV/atom, while 2Sr@B 6 C 6  lies on the convex hull. 
     DFT calculations show that at 50 GPa SrB 3 C 3  is a stable product of exothermic reactions of the pure elements and of readily accessible binary compounds. Therefore, the inventors conducted diamond-anvil cell (DAC) experiments using homogeneous fine-grained mixtures of SrB 6 , SrC 2  and glassy C targeting the stoichiometric reaction SrB 6  + SrC 2  + 4C → 2SrB 3 C 3 . Additional experiments were conducted using only mixtures of binary compounds where the most energetically favorable reaction was calculated as SrB 6  + 3SrC 2  → SrB 3 C 3  + 3SrBC. Mixtures of the powders were compressed in Ne or AI 2 O 3  media and heated near ~2500 K using an infrared fiber laser while synchrotron X-ray diffraction patterns were collected to monitor structural changes in situ. 
     The starting SrB 6  possesses the LaB 6  (Pm 3 m) structure, whereas SrC 2  takes on the acetylide structure of CaC 2  (14 /mmm). When compressed at room temperature, SrB 6  remains in the starting cubic phase and SrC 2  transforms to the R 3 m structure (BaC 2  type) above 14 GPa, eventually appearing amorphous to X-rays above 50 GPa. Upon heating near 2500 K, the intensities of diffraction peaks from the starting compounds vanish and a series of new reflections appear. At 57(3) GPa, these sharp lines were initially indexed to a phase-pure BCC cubic lattice with a = 4.5972(2) Å, in excellent correspondence with the type-VII SrB 3 C 3  clathrate (Pm 3 n) with a calculated lattice parameter of a = 4.593 Å at the same pressure. 
     The calculated X-ray diffraction pattern of the SrB 3 C 3  clathrate is compared with experimental scattering data in  FIG.  3 A  and  FIG.  3 B . Given the nearly complete experimental powder averaging statistics, the quantitative diffraction intensities are representative of atomic positions and are in excellent agreement with the calculated pattern of SrB 3 C 3  clathrate to the experimental resolution limit of 0.75 Å. All allowed reflections with calculated intensity are observed to this limit, consistent with the formation of SrB 3 C 3  clathrate. 
     Given the large contrast in electron density between Sr and the framework atoms, the heavier element dominates the intensity of scattered X-rays. While the formation of SrB 3 C 3  is strongly supported by the stoichiometric conversion of the starting materials, the intensities of the allowed reflections that differentiate the primitive bipartite structure from the BCC sodalite version are in fact negligible, and the possibility for another type of cubic Sr lattice must be considered. The inventors exclude a new elemental Sr lattice that could potentially describe these data based on the requirement of large Sr-Sr distances near 3.8 Å, which are &gt;1.1 Å larger than the distances between Sr atoms in the known metallic phases at this pressure. Other Sr compounds with cubic sublattices are excluded by the experimental P-V EoS, which uniquely distinguishes clathrate cage structures based on agreement with DFT calculations over a broad pressure range from 0-150 GPa. SrB 3 C 3  and other higher-energy cage variants are much less compressible than all other stable and metastable elemental structures and binary/ternary compounds calculated except pure boron allotropes and diamond. 
     The SrB 3 C 3  framework exhibits strong covalent bonding between sp 3 -hybridized B and C atoms, and weak interactions with the Sr guest. This strong sp 3 -hybridized covalent framework guarantees a high value for the bulk modulus (B 0  = 249(3) GPa) and incompressibility of SrB 3 C 3  clathrate. Based on electron count, SrB 3 C 3  should be a hole conductor, and calculations show it is. The reasoning is as follows: All all-carbon, four-coordinate zeolites are insulators at low pressures, closed-shell systems analogous to diamond. A sodalite all-carbon clathrate would be that, and so would isoelectronic [C 3 B 3   ]   3- . SrB 3 C 3  is one electron per formula unit short of this magic (insulator) electron count. Indeed, the band structure ( FIG.  4 B ) shows this-a good gap for one electron more than SrB 3 C 3 . See also  FIG.  4 A . 
     SrB 3 C 3  clathrate is likely the first member of a new class of strong and lightweight sp 3 -bonded carbon-based frameworks with tunable properties. Because boron anions are isoelectronic to carbon atoms in the B-C framework, the bipartite structure should exhibit similar properties to hypothetical pure carbon cages. Furthermore, the ability to trap various kinds of guest atoms in the sp 3 -bonded cages allows the carbon-based clathrates to possess diamond-like mechanical properties with a tunable electronic structure. While SrB 3 C 3  is metallic, the electronic structure may be modified by substituting different guest atoms. For example, LaB 3 C 3 , which the inventors also predict to be thermodynamically stable at high pressure, possesses a band gap due to the balanced electron count. The removal of guest atoms from the cages offers an additional space to explore with the possibility for guest-free structures. Although SrB 3 C 3  is only thermodynamically stable at high pressure, the inventors find that it is recoverable to atmospheric pressure, similar to diamond, which is produced on an industrial scale. At 1 atm, SrB 3 C 3  persists when kept under an inert atmosphere, but it begins to degrade when exposed to the moisture in air over prolonged time (hours). The extent to which alternative clathrate structure types are achievable through different boron substitution/guest schemes are being investigated, and the formation pressure and stability may be optimized. The results may also shed light on the role of element substitution in stabilizing entirely different framework structures. It is likely that the framework stabilization principle observed here with boron is applicable to other elements, including more electronegative ones like nitrogen. There are prospects to obtain a much richer landscape for new carbon-based materials with advanced properties by applying this stabilization principle. 
       FIG.  5    provides an image showing in situ synthesis of SrB 3 C 3  with laser heating in a diamond anvil cell. The central hole contains the sample with a laser hot spot that is rastered over the samples. The reflective box is sample converted to SrB 3 C 3 ; the bottom portion is unconverted. 
       FIG.  6 A  and  FIG.  6 B  demonstrate the superconducting transition temperature of SrB 3 C 3 .  FIG.  6 A  shows phonon dispersion, and phonon density of states (PHDOS) and electron-phonon integral λ(ω) of SrB 3 C 3  at 0 GPa.  FIG.  6 B  shows T c  as a function of pressure. The insets show the evolution of λ and ω log  with pressure. The superconducting transition temperature, T c , was estimated from the Allen~Dynes modified McMillan equation, and a typical value of the Coulomb pseudopotential µ* = 0.1 was used. The calculated T c  is 42 K at ambient pressure. 
       FIGS.  7 A- 7 D  show calculated stress-strain relations of SrB 3 C 3 . More specifically,  FIG.  7 A  and  FIG.  7 B  show calculated tensile stress-strain relations for SrB 3 C 3  with GGA ( FIG.  7 A ) and LDA ( FIG.  7 B ).  FIG.  7 C  and  FIG.  7 D  show calculated shear stress-strain relations for SrB 3 C 3  in the (110) easy cleavage plane with GGA (C) and LDA (D). The stress response along different deformation paths under tensile and shear strains, combined with the lowest peak stress defines the corresponding ideal strength, which is the maximum stress that a perfect crystal can sustain before yielding to a plastic deformation. The established method was applied to determine the stress-strain relations for SrB 3 C 3  under tensile strains in three principal crystallographic directions. The lowest value of calculated peak stress is 24 GPa (31 GPa with LDA) along the &lt;110&gt; direction, which indicate that the &lt;110&gt; direction is the weakest tensile direction, and thus the (110) planes are the easy cleavage planes. The inventors next evaluate the shear stress response in the (110) “easy cleavage plane” of SrB 3 C 3 . The lowest peak shear stress of 25 GPa (34 GPa with LDA) in the (110)[ 1 10] shear direction. These strength values place SrB 3 C 3  as a very hard material with a hardness between 24-34 GPa. 
     According to one aspect of the invention, the inventors have predicted and synthesized the carbon-based clathrate structure with La as the guest atom. The clathrate compound LaB 3 C 3  is different in many respects from other carbon clathrate compounds discussed above. 
     The inventors report a carbon-based clathrate with composition 2La@B 6 C 6  (LaB 3 C 3 ). Like SrB 3 C 3 , the La version crystallizes in the cubic bipartite sodalite structure (type VII clathrate) with La atoms encapsulated within truncated octahedral cages comprised of alternating carbon and boron atoms. The covalent nature of B-C bonding results in a rigid, incompressible framework, and due to the balanced electron count, La 3+ [B 3 C 3 ] 3-  is a semiconductor with an indirect band gap estimated near 1.29 eV. Given that different guest atoms can be substituted within the clathrate cages, it is possible that a broad range of boron-stabilized carbon clathrates with varying properties may be synthesized with ranging physical properties. 
     To probe the possibility of forming the bipartite sodalite-type carbon-boron clathrate with a trivalent guest, the inventors first performed first-principles structure searching computations using the CALYPSO method. Like the case of SrB 3 C 3 , the inventors predicted that LaB 3 C 3  becomes stable under high-pressure conditions. After conducting a computational survey of binary and ternary compounds in the La-B-C system, the inventors estimate that LaB 3 C 3  is the exothermic product of stable binary compounds and/or single components above about 20 GPa ( FIG.  9 A ). Compared with calculations for SrB 3 C 3 , the pressure required for LaB 3 C 3  is reduced by about 20 GPa—a much more favorable synthetic condition that we tentatively attribute to the overall balanced electron count. 
     According to one aspect of the invention, LaB 3 C 3  was synthesized according to the reaction LaB 6  + LaC 2  + 4C → 2LaB 3 C 3  at high pressure using diamond anvil cells. After heating near 2500 K with an infrared fiber laser, the starting precursor was converted into an essentially phase-pure cubic phase ( FIG.  9 A ), as expected for the clathrate (space group Pm 3 n). At 57(2) GPa, the experimental lattice parameter a = 4.66 Å is in good agreement with a DFT optimized LaB 3 C 3  structure with a = 4.63 Å at 60 GPa (PBEGW), and Rietveld refinement of the powder data show excellent agreement with the predicted clathrate structure ( FIG.  9 B ). In addition, the stoichiometric conversion from the homogeneous precursor to a single phase strongly supports the 1 La:3B:3C composition, and a uniform single-phase region near this composition was also confirmed using energy-dispersive X-ray spectroscopy (EDX) mapping on the recovered product. Experimentally, the inventors confirmed the formation of LaB 3 C 3  at 38 GPa but did not observe conversion at 35 GPa, placing a tentative bound on the minimum required synthesis pressure. Compared with SrB 3 C 3  synthesized near 50 GPa, this represents a significant pressure reduction and suggests promise for future stabilization efforts at even milder conditions. 
     Synchrotron X-ray diffraction patterns collected during decompression confirm the clathrate framework of LaB 3 C 3  based on compressibility. The decompression data were fitted using a third-order Birch-Murnaghan equation of state (EoS) with the zero-pressure bulk modulus B 0  = 254(9) GPa and it’s derivative B 0 ′= 3.6(4). The bulk modulus compares favorably with calculations based on DFT ( FIG.  11   ), noting that LDA significantly underestimates the volume for a given pressure. LaB 3 C 3  is very incompressible with B 0  in the range of technical ceramics used for armor, such as B 4 C (B 0  = 243(3) GPa), and the observed bulk modulus cannot be explained by any structures other than the clathrate. Similar to SrB 3 C 3 , the La counterpart is recoverable to ambient conditions, but must be preserved carefully under an inert atmosphere. When samples are exposed to air/moisture, they degrade into what is presumed to be hydrolysis products. Several attempts to prepare samples for transmission electron microscopy by focused ion beam milling were unsuccessful and resulted in oxidized amorphous material. 
     Like SrB 3 C 3 , LaB 3 C 3  is predicted to crystallize in the cubic bipartite sodalite structure ( FIG.  20   ). This structure type is known for the clathrate hydrate HPF 6 .6H 2 O and is related to inorganic compounds such as BaPd 2 P 4  through a small structural distortion, as well as superhydrides such as Ca/YH 6  and CaYH 12 . The clathrate structure is comprised of a single truncated octahedral B 12 C 12  cage with six four-sided faces and eight six-sided faces (4 6 6 8 ). All faces are shared between cages so that the single cage type can fully tile three-dimensional space without any voids. All atoms are on special crystallographic positions so that there is only one unique C-B distance and one unique C(B)-Sr distance. The square and hexagonal faces result in a mixture of 90° and 120° C-B-C bond angles. The incorporation of B within the framework reduces the energetic penalty for the 90° bond angles. 
     Given the strong contrast in electron density between framework and guest atoms, powder diffraction can neither distinguish the positions and occupancies of boron and carbon (i.e., the “coloring” of the clathrate lattice) nor the weak intensity reflections that differentiate the primitive (bipartitie) and body-centered cubic (sodalite) lattices. The inventors therefore conducted single-crystal diffraction on suitable grains produced by laser heating at high temperatures above 3000 K. Similar to the case of SrB 3 C 3 , the observation of {120} reflections with significant intensity confirms the primitive space groupPm 3 n and the ordered bipartite structural model produces the lowest reliability factor between experimental observations compared with other possible coloring schemes. 
     At 0 GPa, the LaB 3 C 3  lattice is slightly expanded compared with that of SrB 3 C 3  (a = 4.934 Å vs. 4.868 Å). This is somewhat surprising as the ionic radius of La 3+  is expected to be smaller than that of Sr 2+  (1.36 Å vs. 1.44 Å for 12-fold coordination), see http://abulafia.mt.ic.ac.uk/ shannon/radius.php. The increased cell size and thus B-C bond length in LaB 3 C 3  (1.74 Å vs 1.72 Å SrB 3 C 3 ) can be understood as additional charge transferred from guest to host. Indeed, charge density analysis ( FIG.  21 A  and  FIG.  21 B ) show that 1.82 e -  are transferred from La to the B/C host lattice, whereas 1.35 e -  are transferred for the case of Sr. 
     The electronic band structure of LaB 3 C 3  also reveals the influence of the additional electron. Unlike SrB 3 C 3 , which is one electron short of a band gap, LaB 3 C 3  shows an indirect gap of 1.29 eV in the ┌ → L direction. The top of the valence band is composed primarily of hybrid p states of B and C, while the conduction band is constituted primarily of d states from La. LaB 3 C 3  is comparable to B-substituted silicon frameworks such as K 7 B 7 Si 39  and LiBSi 2  with electron precise frameworks in which the metals act as electron donors, as in Zintl phases. 
     Summarizing the above, a second carbon-boron clathrate in the bipartite sodalite structure was prepared with La guest atoms. Unlike metallic SrB 3 C 3 , LaB 3 C 3  has a balanced electron count and is predicted to be a semiconductor. This balanced count is also likely responsible for the improved synthetic conditions, i.e., &gt;20% reduction in synthetic pressure. The observation of two clathrate structures indicates that many other type-VII clathrates may be possible with different guest atoms. The inventors suggest that entirely different clathrate structure types also may be possible with different carbon:boron ratios and cage types. Carbon-based clathrates thus represent a class of new materials with highly tunable properties. 
       FIG.  8    provides a powder XRD pattern of LaB 3 C 3  during decompression at 5 GPa (top), and after opening the cell in air (bottom). Bragg positions from cubic LaB 3 C 3  are marked by blue lines. The diffraction intensity was limited by the opening angle above 2q=14°. When exposed to air, LaB 3 C 3  degrades into unknown, presumably, hydrolysis products. Diffraction from unreacted LaB 3 C 3  can still be observed in air several minutes after the cell was opened. 
       FIG.  10    provides phonon dispersion for LaB 3 C 3  at 1 atmosphere, showing dynamic stability. 
       FIG.  11    provides stress-strain calculations for LaB 3 C 3  with GGA and LDA. The calculations indicate a high ideal strength with hardness in near 30 GPa. 
     According to one aspect of the invention, the starting precursors to form the carbon clathrate compounds are borides, carbides and pure carbon. 
     Both the La and the Sr-based carbon clathrates are recoverable to ambient conditions but must be preserved under an inert atmosphere (e.g., argon) due to degradation in air. 
     SrB 3 C 3 , is predicted to be a hole metal with a high superconducting transition temperature near 42 K. LaB 3 C 3  is predicted to be a semiconductor. The clathrate structures exhibit a high bulk modulus near 25 GPa and are predicted to exhibit high strength and hardness. 
     According to one aspect embodiment of the invention, the carbon-based clathrate structures are composed of entirely sp 3  hybridized carbon and boron, which results in diamond-like bonding. In an embodiment of the invention, different guest atoms may be substituted within the cages to create a new class of diamond-like materials with tunable properties. Many other guest atoms are calculated to be possible within the clathrate cages. 
     Different boron substitution schemes may result in entirely different clathrate structures. Other elemental substitution schemes for the host lattice may be possible, for example N instead of B. 
     EXAMPLES 
     The following is a more detailed description of the present invention with reference to working examples for prediction, synthesis of compounds and characterization of the same. The present invention is in no way limited to the following examples. 
     Calculations 
     Structure-searching calculations were performed using the CALYPSO structure prediction method based on the global minimization of free energy using ab initio total-energy calculations. This method was benchmarked using various systems, ranging from elements to binary and ternary compounds. Total energy calculations were performed in the framework of density functional theory within the Perdew-Burke-Ernzerhof generalized gradient approximation as implemented in the VASP (Vienna Ab Initio Simulation Package) code. The projector-augmented wave (PAW) method (50) was adopted with the PAW potentials taken from the VASP library where 4s 2 4p 6 5s 2 , 2s 2 2p 1  and 2s 2 2p 2  are treated as valence electrons for Sr, B and C atoms, respectively. The use of a plane-wave kinetic energy cutoff of 520 eV and dense k-point sampling, adopted here, were shown to give excellent convergence of total energies. Electronic charges were calculated using a Bader charge analysis scheme using a 600 × 600 × 600 Fast Fourier Transform grid. Phonon dispersion calculations were performed to determine the dynamical stability of the predicted structures by using the finite displacement approach, as implemented in the Phonopy code. Electron-phonon coupling calculations for superconducting properties of stable phases were performed using density-functional perturbation theory (DFPT) with the Quantum-ESPRESSO package. To study interatomic interactions, crystal orbital Hamilton population (COHP) analysis was performed using the LOBSTER package. 
     Synthesis 
     Strontium metal was purified by sublimation (950° C., dynamic vacuum) onto Mo foil from a graphite crucible; graphite powder was pre-treated at 950° C. under dynamic vacuum for 16 h to remove adsorbed species. Strontium carbide (SrC 2 ) was prepared by heating 2.780 g Sr (31.7 mmol) with 0.371 g graphite powder (30.9 mmol) in a capped graphite crucible under dynamic vacuum at 825° C. for 16 h. Excess Sr was subsequently sublimed at 950° C. and the resulting, pale grey SrC 2  powder (1.400 g, 81% yield) isolated from the crucible under Ar. Powder XRD showed a small SrO impurity. SrB 6  (EPSI Metals, 99.5 %) and glassy carbon (Sigma Aldrich, 99.95%) were purchased commercially and used without further purification. Binary (SrB 6  + 3SrC 2 ) and ternary (SrB 6  + SrC 2  + 4C) mixtures were prepared under an inert Ar atmosphere, sealed, then milled vigorously using Si 3 N 4  media at 600 rpm for one-minute cycles over ~12 hours. The milled powders were removed from the media in an Ar glovebox and thin plates (10 µm) were pressed between two diamond anvils with 1 mm culets, then loaded into DAC sample chambers utilizing 100 - 300 pm culets and Re gaskets. Ne or alumina plates served as the pressure transmitting media and thermal insulation from the diamond anvils. After being compressed to the target pressure between 50-150 GPa, samples were heated to -2500 K using the double-sided laser heating systems at HPCAT or GSECARS. 
     X-Ray Diffraction 
     In situ X-ray diffraction patterns were collected at the Advanced Photon Source, Sector 16, HPCAT using a monochromatic wavelength of 0.4066 Å and at Sector 13, GSECARS using a monochromatic wavelength of 0.3344 Å. The X-ray beam was focused on the sample and scattered X-rays were detected using a PILATUS 1M or MARCCD detector. The sample-to-detector distance and geometrical parameters were calibrated using CeO 2  and LaB 6  standards with the DIOPTAS software. Pressure was calibrated using the equations of state of Ne and and/or Al 2 O 3  and cross-checked using the SrO equation of state and ruby fluorescence in some samples. Rietveld refinements of XRD patterns were conducted using Powdercell and GSAS with EXPGUI. 
       FIG.  12    provides raw XRD patterns before and after heating. The starting material is a mixture of SrB 6  (Pm 3 m), SrC 2  (/4 /mmm) and glassy C in a 1:1:4 molar ratio, e.g., 2Sr:6C:6B. SrC 2  and glassy carbon appear to be amorphous to X-rays and are not distinguishable from crystalline SrB 6 . After heating to a maximum temperature of ∼2900 K, the stoichiometric 1Sr:3C:3B mixture was converted to nearly phase pure SrB 3 C 3  with trace residual SrB 6 . 
       FIG.  13    provides raw XRD patterns with increasing pressure and heating. The starting material is a mixture of SrB 6 , (pm 3 m) and SrC2 (/4 /mmm) in a 1:3 molar ratio. Above 14 GPa, SrC 2  transforms to the R 3 m structure. All Bragg peaks broaden significantly and diminish in intensity with pressure due to significant stress accumulation prior to heating. Above ca. 50 GPa, SrC 2  appears to become amorphous to X-rays; one broad feature is potentially attributable to the R 3 m phase. After heating at -150 GPa, many new Bragg peaks appear, which can be attributed to SrB 3 C 3 , SrBC and a high-pressure form of SrB 6 , as described below. 
       FIG.  14    demonstrates Sr-B-C phase identification at 155(6) GPa with Ne medium. The starting material is a mixture of SrB 6  (Pm 3 m) and SrC 2  (/4 /mmm) in a 1:3 molar ratio. The data indicate the reaction SrB 6  + 3SrC 2  → SrB 3 C 3  + 3SrBC (and some unconverted SrB 6  in a high-pressure phase), which is calculated to be the most energetically favorable reaction for these conditions. Experimental data are shown as black points connected by a thin black line. Each phase is labeled with a different symbol; tick marks below the pattern indicate allowed reflections for SrB 3 C 3  (Pm 3 m). Rietveld refinement was conducted using only a Gaussian peak width and scale parameter. All atomic positions were taken from DFT-optimized structures shown in Tables S1 and S2 below. Powder averaging statistics vary between phases. Experimentally refined lattice parameters are compared with DFT (PBE) optimizations at 155 GPa. The inset to the right shows quantitative intensity agreement for SrB 3 C 3  to the limiting resolution of d = 1.16 Å. The 2D cake image is presented at the top of the figure with arrows indicating prominent sharp lines for SrB 3 C 3 . 
       FIG.  15    shows Sr-B-C phase identification at 142(10) GPa with Al 2 O 3  medium. The starting material is a mixture of SrB 6  (pm 3 m) and SrC 2  (/4 /mmm) in a 1:3 molar ratio. The data indicate the reaction SrB 6  + 3SrC 2  → SrB 3 C 3  + 3SrBC (and some unconverted SrB 6  in a high-pressure phase), which is calculated to be the most energetically favorable reaction for these conditions. Experimental data are shown as black points connected by a thin black line. Each phase is labeled with a different symbol; tick marks below the pattern indicate allowed reflections for SrB 3  C 3  (Pm 3 m). Rietveld refinement was conducted using only a Gaussian peak width and scale parameter. All atomic positions were taken from DFT-optimized structures shown in Tables S1 and S2 below. Powder averaging statistics vary between phases. Experimentally refined lattice parameters are compared with DFT (PBE) optimizations at 140 GPa. The inset to the right shows quantitative intensity agreement for SrB 3 C 3  to the limiting resolution of d = 0.93 Å. The 2D cake image is presented at the top of the figure with arrows indicating prominent sharp lines for SrB 3 C 3 . Single-crystalline Al 2 O 3  peaks were masked, but some powder intensity was unavoidable due to the breaking of crystals during compression to 140 GPa. 
       FIG.  16    shows stability of SrB 3 C 3  and analysis of different possible stoichiometries. Relatedly,  FIG.  1 B  shows predicted formation enthalpies of various Sr-B-C compounds with respect to elemental decomposition at 50 GPa. Compounds with enthalpy data represented by red points are on the convex hull and thermodynamically stable against decomposition. Black points show the formation enthalpies of metastable structures found in the CALYPSO structure searches. Blue points above the convex hull indicate the formation enthalpies of 600 additional structures containing up to 112 atoms with compositions SrB x C y  (excluding SrB 3 C 3 ). These 600 structures were determined from searches wherein the Sr atoms were constrained to a tolerance within the experimentally determined cubic positions, while the C and B atoms were allowed to freely fluctuate. No composition was found to be more stable than SrB 3 C 3 .  FIG.  16    shows the calculated EoS of the 600 SrB x C y  compounds presented in  FIG.  1 B . The red line shows the EoS of SrB 3 C 3  clathrate, which agrees with experiment. The grey lines represent the EoS of different SrB x C y  compounds with calculated B0 &lt; 150 GPa, while the cyan lines indicate the structures with calculated B0 in the range of 210 to 230 GPa. After a careful analysis, we found that all the structures from the cyan group can be classified as type-VII clathrate structures (e.g., 4668 cages with different B:C compositions) whereas structures from the grey group are not cage-like structures and are more compressible. As the experimental EoS is located in the cyan group, and SrB 3 C 3  is the lowest energy structure determined, the synthesized product is only consistent with a clathrate structure that should be identical or similar to the structure of cubic SrB 3 C 3 . 
       FIG.  17    shows X-ray diffraction of SrB 3 C 3  at atmospheric pressure. After synthesis at ~57 GPa, the cell was decompressed, the neon gas was released and XRD patterns were collected after the cell was sealed at 1 atm. The refined lattice parameter of a = 4.868 Å compares well with DFT-GGA calculations where a = 4.88 Å. Red bars indicate calculated positions for the cubic clathrate structure. When left in open air, SrB 3 C 3  appears to degrade and is thus sensitive to moisture and/or oxygen. 
       FIG.  18    provides electronic band structure for SrB 3 C 3  at 200 GPa projected onto atomic orbitals represented by different colors. The width of each band is proportional to the weight of the corresponding orbital character. The dashed line indicates the Fermi energy. 
       FIG.  19 A  and  FIG.  19 B  provide energetic stabilities as a function of pressure.  FIG.  19 A  shows calculated enthalpies (H) per SrC 2  unit as a function of pressure for various structures with respect to /4 /mmm structure.  FIG.  19 B  shows calculated enthalpies (H) per SrB 6  unit as a function of pressure for various structures with respect to Pm 3 m structure. 
       FIG.  22    provides an EDX map of recovered sample showing homogeneous elemental distribution in heated region (25×25 µm 2 ) with average composition La 1.00±0.07 B 2.95±0. 06 C 4.20±0.09 O 0.46±0.07.  The carbon content is biased by adventitious sources, also observed on pure standards. The color composite map is layered over the SEM image and individual elements from the same region are shown below. 
     Table 1 provides single-crystal diffraction and refinement paramters obtained from LaB 3 C 3  crystal synthesized at 57 Gpa and &gt;3000 K. The analysis confirms the bipartite sodalite framework and primitive cubic space group.  
     
       
         
          TABLE 1
           
               
               
             
               
                 Structural Information for LaB 3 C 3  from high-pressure single-crystal X-ray diffraction. 
               
               
                 Empirical formula 
                 LaB 3 C 3 
 
               
             
            
               
                 Formula weight 
                 207.37 
               
               
                 T/K 
                 293(2) 
               
               
                 Crystal system 
                 cubic 
               
               
                 Space group 
                 Pm 3 m 
               
               
                 a /Å 
                 4.6712(16) 
               
               
                 V/Å 3 
 
                 101.93(10) 
               
               
                 Z 
                 2 
               
               
                 P calc /g·cm -3 
 
                 6.757 
               
               
                 µ /mm -1 
 
                 20.514 
               
               
                 F(000) 
                 180.0 
               
               
                 Crystal size/mm 3 
 
                 0.005 × 0.005 × 0.005 
               
               
                 Radiation 
                 synchrotron (λ = 0.40663) 
               
               
                 2Θ range for data collection/° 
                 9.988 to 29.844 
               
               
                 Index ranges 
                 -4 ≤ h ≤ 5, 
               
               
                 -4 ≤ k ≤ 4, 
               
               
                 -4 ≤ / ≤ 3 
               
               
                 Reflections collected 
                 81 
               
               
                 Independent reflections 
                 21 [R int  = 0.0452, R sigma  = 0.0280] 
               
               
                 Data/restraints/parameters 
                 21/0/4 
               
               
                 Goodness-of-fit on F 2 
 
                 1.109 
               
               
                 Final R indexes [l&gt;=2σ (l)] 
                 R 1  = 0.0290, wR 2  = 0.0633 
               
               
                 Final R indexes [all data] 
                 R 1  = 0.0297, wR 2  = 0.0658 
               
               
                 Largest diff. peak/hole / eÅ -3 
 
                 +0.79/-0.77 
               
            
           
         
       
     
     As the present invention may be embodied in several forms without departing from the spirit or essential characteristics thereof, it will be understood that the invention is not limited by the details of the foregoing description, unless otherwise specified, but rather should be construed broadly within its spirit and scope as defined in the appended claims, and therefore all changes and modifications that fall within the metes and bounds of the claims. Accordingly, the invention is defined by the appended claims wherein: