CARBON-BASED CLATHRATE COMPOUNDS

Abstract

The present invention provides carbon-based clathrate compounds, including a carbon-based clathrate compound that includes a clathrate lattice with atoms of at least one element selected from the group consisting of carbon and boron as a host cage structure and guest atoms encapsulated within the clathrate lattice. The guest atoms can constitute single elements or multiple distinct elements. In one embodiment, the invention provides a carbon-based clathrate compound of the formula NdB.sub.3C.sub.3. In another embodiment, the invention provides a carbon-based clathrate-like compound of the formula CaB.sub.3C.sub.3. In another embodiment, the invention provides a carbon-based clathrate compound of the formula Na.sub.xSr.sub.(1-x)B.sub.3C.sub.3.

Claims

1. A carbon-based clathrate compound comprising: (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; and (ii) guest atoms encapsulated within the clathrate lattice wherein the guest atoms are selected from a group consisting of neodymium, scandium, yttrium, cerium, praseodymium, promethium, samarium, europium, erbium, thulium, ytterbium, lutetium, and lanthanum.

2. A carbon-based clathrate compound comprising: (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; and (ii) at least two different types of guest atoms encapsulated within the clathrate lattice.

3. The carbon-based clathrate compound of claim 2, wherein the at least two guest atoms are selected from the group consisting of a first guest atom pair sodium and strontium and a second guest atom pair calcium and strontium.

4. A compound of the formula CaB.sub.3C.sub.3.

5. The compound of claim 4, wherein the compound is a carbon-based clathrate-like structure.

6. The compound of claim 5, wherein the structure of the carbon-based clathrate-like structure is hexagonal.

7. The compound of claim 5 comprising a clathrate lattice-like structure with atoms selected from the group consisting of carbon and boron as a host environment structure.

8. The compound of claim 5, wherein the clathrate lattice-like structure is formed of sp.sup.3 hybridized carbon and boron.

9. The compound of claim 5, wherein calcium is a guest atom encapsulated within the clathrate lattice-like structure.

10. The compound of claim 8, wherein the clathrate lattice-like structure comprises star-shaped channels down the crystallographic c-axis.

11. A compound of the formula Na.sub.xSr.sub.(1-x)B.sub.3C.sub.3.

12. The compound of claim 11, wherein the compound is a carbon-based clathrate.

13. The compound of claim 12, wherein the structure of the carbon-based clathrate is cubic bipartite socialite.

14. The compound of claim 12 comprising a clathrate lattice with atoms selected from the group consisting of carbon and boron as a host cage structure.

15. The compound of claim 12, wherein the clathrate lattice is formed of sp.sup.3 hybridized carbon and boron.

16. The compound of claim 12, wherein sodium is a first guest atom encapsulated within the clathrate lattice and strontium is a second guest atom encapsulated within the clathrate lattice.

17. The compound of claim 15, wherein the clathrate lattice comprises cages, each cage comprising 24 atoms with six four-sided faces and eight six-sided faces.

18. The compound of claim 16, wherein the clathrate lattice comprises cages, each cage comprising 24 atoms with six four-sided faces and eight six-sided faces.

19. The compound of claim 18, wherein the first guest atom is encapsulated within the group of six four-sided faces and the second guest atom is encapsulated within the group of eight six-sided faces.

20. A compound of the formula Ca.sub.xSr.sub.(1-x)B.sub.3C.sub.3.

21. The compound of claim 20, wherein the compound is a carbon-based clathrate.

22. The compound of claim 21, wherein the structure of the carbon-based clathrate is cubic bipartite socialite.

23. The compound of claim 21 comprising a clathrate lattice with atoms selected from the group consisting of carbon and boron as a host cage structure.

24. The compound of claim 21, wherein the clathrate lattice is formed of sp.sup.3 hybridized carbon and boron.

25. The compound of claim 21, wherein calcium is a first guest atom encapsulated within the clathrate lattice and strontium is a second guest atom encapsulated within the clathrate lattice.

26. The compound of claim 24, wherein the clathrate lattice comprises cages, each cage comprising 24 atoms with six four-sided faces and eight six-sided faces.

27. The compound of claim 5, wherein the clathrate lattice comprises cages, each cage comprising 24 atoms with six four-sided faces and eight six-sided faces.

28. The compound of claim 27, wherein the first guest atom is encapsulated within the group of six four-sided faces and the second guest atom is encapsulated within the group of eight six-sided faces.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0022] 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.

[0023] The invention is more fully described by reference to the following detailed description and the accompanying drawings wherein:

[0024] FIG. 1A is a ternary phase diagram of a SrBC system at 50 GPa. Circles (green) represent thermodynamically stable compounds while squares (orange) represent additional metastable compositions used in the search.

[0025] FIG. 1B is a ternary convex hull for the SrBC 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.

[0026] FIG. 2 provides the structure of an SrB.sub.3C.sub.3 clathrate. The structure also may represent the LaB.sub.3C.sub.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 (4668). Different color cages are used to emphasize the stacking of cages that tile 3D space.

[0027] FIG. 3A provides experimental x-ray diffraction data collected at 57 (3) GPa with Rietveld refinement of the SrB.sub.3C.sub.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.sub.3C.sub.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.sub.3C.sub.3 phase. The inset shows a magnified view at high angle with sharp SrB.sub.3C.sub.3 peaks to a limiting resolution of 0.75 .

[0028] FIG. 3B provides experimental EoS and calculated EoS data for the SrB.sub.3C.sub.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.

[0029] FIG. 4A provides electronic properties of SrB.sub.3C.sub.3 and in particular a two-dimensional electron localized function (ELF) for the compound.

[0030] FIG. 4B provides electronic band structures for SrB.sub.3C.sub.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.sub.3C.sub.3 at 0 GPa are shown to the right. The Fermi energy is set to 0 eV (dashed line).

[0031] FIG. 5 provides an image showing in situ synthesis of SrB.sub.3C.sub.3 with laser heating in a diamond anvil cell.

[0032] FIG. 6A shows phonon dispersion, phonon density of states and electron-phonon integral of SrB.sub.3C.sub.3 at 0 GPa.

[0033] FIG. 6B shows Tc as a function of pressure for SrB.sub.3C.sub.3.

[0034] FIG. 7 (A-D) show calculated stress-strain relations of SrB.sub.3C.sub.3. FIG. 7A and FIG. 7B show calculated tensile stress-strain relations for SrB.sub.3C.sub.3 with GGA (FIG. 7A) and LDA (FIG. 7B). FIG. 7C and FIG. 7D show calculated shear stress-strain relations for SrB.sub.3C.sub.3 in the (110) easy cleavage plane with GGA (FIG. 7C) and LDA (FIG. 7D).

[0035] FIG. 8 provides an XRD pattern showing LaB.sub.3C.sub.3 recovered to low pressure. The Powder XRD pattern of LaB.sub.3C.sub.3 is shown during decompression at 5 GPa (top), and after opening the cell in air (bottom). Bragg positions from cubic LaB.sub.3C.sub.3 are marked by blue lines. The diffraction intensity was limited by the opening angle above 20=14. When exposed to air, LaB.sub.3C.sub.3 degrades into unknown, presumably, hydrolysis products. Diffraction from unreacted LaB.sub.3C.sub.3 can still be observed in air several minutes after the cell was opened.

[0036] FIG. 9A shows computed stability of LaB.sub.3C.sub.3 via different reaction pathways indicating exothermic formation from pure elements and from stable binaries above ca. 20 GPa.

[0037] FIG. 9B shows Rietveld refinement of LaB.sub.3C.sub.3 at 57 (2) GPa with a=4.66307 (8) , wRp-Bknd=0.042, .sup.2=3.8. The sample is nearly phase pure with a trace impurity of unreacted LaB.sub.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.sub.6 and Ne; ballmilled LaC.sub.2 and glassy carbon appear amorphous.

[0038] FIG. 9C provides unit cell volume of LaB.sub.3C.sub.3 as a function of pressure. Experimental data were fitted using a third-order Birch-Murnaghan EoS with B.sub.0=254 (9) GPa and B.sub.0=3.6 (4). DFT results yield B.sub.0=248 GPa with B.sub.0=3.8 (PBEGW) and B.sub.0=268 GPa with B.sub.0=3.9 (LDA).

[0039] FIG. 10 shows phonon dispersion for LaB.sub.3C.sub.3 at 1 atm showing dynamic stability.

[0040] FIG. 11 provides stress-strain calculations for LaB.sub.3C.sub.3 with GGA and LDA., indicating a high ideal strength with hardness near 30 GPa.

[0041] FIG. 12 provides raw XRD patterns of SrB.sub.3C.sub.3 before and after heating.

[0042] FIG. 13 provides raw XRD patterns of SrB.sub.3C.sub.3 with increasing pressure and heating.

[0043] FIG. 14 provides SrBC phase identification at 155 (6) GPa with Ne medium.

[0044] FIG. 15 provides SrBC phase identification at 142 (10) GPa with Al.sub.2O.sub.3 medium.

[0045] FIG. 16 presents data showing the stability of SrB.sub.3C.sub.3 and analysis of possible stoichiometries.

[0046] FIG. 17 shows X-ray diffraction of SrB.sub.3C.sub.3 at atmospheric pressure.

[0047] FIG. 18 provides electronic band structure for SrB.sub.3C.sub.3 at 200 GPa projected onto atomic orbitals.

[0048] FIG. 19A and FIG. 19B provide energetic stabilities as a function of pressure, with FIG. 19A showing calculated enthalpies per SrC.sub.2 unit as a function of pressure for various structures with respect to 14/mmm structure, and FIG. 19B showing calculated enthalpies per SrB.sub.6 unit as a function of pressure for various structures with respect to Pm3m structure.

[0049] FIG. 20 provides the structure of LaB.sub.3C.sub.3 viewed along different crystallographic directions. C and B atoms make up the framework, while La is trapped in the cages. The single 4668 cage has one unique B-C distance of 1.74 at atmospheric pressure.

[0050] FIG. 21A shows the difference in charge density for LaB.sub.3C.sub.3 (crystal density minus isolated atomic density) viewed in the (100) plane at 0 GPa.

[0051] FIG. 21B shows the electronic band structure (HSE) of LaB.sub.3C.sub.3 at 0 GPa with density of states. Contributions from different orbitals are indicated.

[0052] FIG. 22 provides an EDX map of recovered sample showing showing homogeneous elemental distribution in heated region (2525 m.sup.2) with average composition La.sub.1.000.07B.sub.2.950.06C.sub.4.200.09O.sub.0.460.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.

[0053] FIG. 23 shows the computed phonon dispersion relations for various MB.sub.3C.sub.3 compounds exhibiting dynamic stability at 1 atm.

[0054] FIG. 24 shows formation enthalpies of MB.sub.3C.sub.3 compounds versus pressure, showing calculated enthalpies per MB.sub.3C.sub.3 unit as a function of pressure for various clathrates in the Pm3n structure type.

[0055] FIG. 25A shows single-crystal refinement parameters for NdB.sub.3C.sub.3 at 50 GPa.

[0056] FIG. 25B shows the electronic band structures for NdB.sub.3C.sub.3 at 50 GPa projected onto atomic orbitals, where the width of each band is proportional to the weight of the corresponding orbital character. The Fermi energy is set to 0 eV (dashed line).

[0057] FIG. 25C shows the crystal structure of NdB.sub.3C.sub.3 where black, purple and orange spheres represent carbon, boron and neodymium atoms, respectively. The cubic structure is composed of face-sharing boron-carbon cages that encapsulate Nd atoms. Each cage contains 24 atoms with six four-sided faces and eight six-sided faces (4668).

[0058] FIG. 26A shows Powder X-ray diffraction patterns for binary Type-VII clathrates produced by CaB.sub.6+SrC.sub.2+4 C.fwdarw._CaSrB.sub.6C.sub.6. Peaks with asterisks are from unreacted precursor.

[0059] FIG. 26B shows Powder X-ray diffraction patterns for binary Type-VII clathrates produced by 2SrB.sub.6+Na.sub.2C.sub.2+10 C.fwdarw._2NaSrB.sub.6C.sub.6 (right). Peaks with asterisks are from unreacted precursor.

[0060] FIG. 27A shows a sample geometry for in situ electrical transport with HPHT synthesis of MB.sub.3C.sub.3 compounds.

[0061] FIG. 27B shows an image of in situ synthesis of SrB.sub.3C.sub.3 with laser heating in a diamond anvil cell.

[0062] FIG. 27C shows experimental data showing x-ray diffraction data for SrB.sub.3C.sub.3 before and after heating.

[0063] FIG. 28A shows experimental data showing the resistance drop during cooling of a first sample of SrB.sub.3C.sub.3 to a base temperature of 5 K at 54 GPa. The heating trace (green overlay) follows the same path with no detectable hysteresis. The inset shows the determination of the Tc onset by extrapolation of linear regions.

[0064] FIG. 28B shows experimental data showing the normalized resistance for several cooling runs at different pressures for the first sample of SrB.sub.3C.sub.3 and a second sample of SrB.sub.3C.sub.3. The inset shows the pressure dependence of experimental T.sub.c onset values SrB.sub.3C.sub.3 (points) compared with calculations using different values of * (lines).

[0065] FIG. 28C shows experimental data showing the resistance of a third sample of SrB.sub.3C.sub.3 at 40 GPa as a function of temperature and magnetic field.

[0066] FIG. 28D shows experimental data showing the upper critical field, .sub.0H.sub.c2(T), for the transition onset of SrB.sub.3C.sub.3 (black points) compared with the WHH model (=0, blue curve). A linear fit (red dashed line) is shown to guide the eye.

[0067] FIG. 29A shows the two-dimensional structure of the compound RbSrB.sub.6C.sub.6.

[0068] FIG. 29B shows the computed phonon dispersion, phonon density of states and electron-phonon integral of the compound RbSrB.sub.6C.sub.6 with dynamic stability at 1 atm.

[0069] FIG. 30A shows the resistance curve versus temperature of the compound Na.sub.xSr.sub.yB.sub.3C.sub.3

[0070] FIG. 30B shows the T.sub.c versus P curve of the compound SrB.sub.3C.sub.3 versus the T.sub.c versus P curve of a Na-doped Type-VII clathrate.

[0071] FIG. 31A shows the two-dimensional structure of the clathrate-like compound h-CaB.sub.3C.sub.3.

[0072] FIG. 31B shows the three-dimensional structure of the clathrate-like compound h-CaB.sub.3C.sub.3 along the c-axis.

[0073] FIG. 31C shows electronic band structure and density of states data for the compound h-CaB.sub.3C.sub.3. The Fermi energy is set to 0 eV (dashed line).

DETAILED DESCRIPTION OF THE INVENTION

[0074] Note: For purposes of this application, U.S. application Ser. No. 16/810,388 filed Mar. 5, 2020 and U.S. application Ser. No. 17/951,939 filed Sep. 23, 2022, both of which further claim priority to U.S. Provisional Application No. 62/814,024 filed Mar. 5, 2019 and all three of which are incorporated herein by reference, will be referred to collectively as the Prior Applications.

[0075] Carbon clathrate is an impressive 3D sp.sup.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.sup.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 previously been reported due to tremendous challenges associated with their synthesis.

[0076] 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 previously. 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.

[0077] 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.sub.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 previously, although three-dimensional polymers of C.sub.60 have been suggested to resemble carbon clathrate-like structures.

[0078] 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 previously 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.sup.3 carbon materials to a new class with tunable properties.

[0079] The inventors conducted an extensive computational search in the SrBC system (including pure elements, binaries and ternaries from SrB.sub.xC.sub.y with 0x, y6) 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 SrBC system, several high-pressure compounds were determined to be thermodynamically stable with respect to elemental mixtures (see FIG. 1A and FIG. 1B). At 50 GPa, the hexagonal P6.sub.3/mmc and -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.sub.5C.sub.2, SrC, SrC.sub.2, Sr.sub.5C.sub.2, SrB, SrB.sub.2, SrB.sub.4, and SrB.sub.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.

[0080] Two stable ternary compounds were predicted at 50 GPa. The first, hexagonal SrBC (space group P6.sub.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 Pm3n)) with composition 2Sr@B.sub.6C.sub.6 (SrB.sub.3C.sub.3) and has the type-VII clathrate structure, known for the clathrate hydrate HPF.sub.6.Math.6H.sub.2O and related to inorganic compounds such as BaPd.sub.2P.sub.4 through a small structural distortion. The topology of SrB.sub.3C.sub.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.

[0081] The SrB.sub.3C.sub.3 clathrate framework (FIG. 2) is composed of a single truncated octahedral cage with six four-sided faces and eight six-sided faces (4668). 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.sub.12 is metastable by 0.667 eV/atom, while 2Sr@B6C.sub.6 lies on the convex hull.

[0082] DFT calculations show that at 50 GPa SrB.sub.3C.sub.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.sub.6, SrC.sub.2 and glassy C targeting the stoichiometric reaction SrB.sub.6+SrC.sub.2+4C.fwdarw.2SrB.sub.3C.sub.3. Additional experiments were conducted using only mixtures of binary compounds where the most energetically favorable reaction was calculated as SrB.sub.6+3SrC.sub.2.fwdarw.SrB.sub.3C.sub.3+3SrBC. Mixtures of the powders were compressed in Ne or Al.sub.2O.sub.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.

[0083] The starting SrB.sub.6 possesses the LaB.sub.6 (Pm3m) structure, whereas SrC.sub.2 takes on the acetylide structure of CaC.sub.2 (l4/mmm). When compressed at room temperature, SrB.sub.6 remains in the starting cubic phase and SrC.sub.2 transforms to the R3m structure (BaC.sub.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.sub.3C.sub.3 clathrate (Pm3n) with a calculated lattice parameter of a=4.593 at the same pressure.

[0084] The calculated X-ray diffraction pattern of the SrB.sub.3C.sub.3 clathrate is compared with experimental scattering data in FIG. 3A and FIG. 3B. 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.sub.3C.sub.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.sub.3C.sub.3 clathrate.

[0085] 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.sub.3C.sub.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 SrSr distances near 3.8 , which are >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.sub.3C.sub.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.

[0086] The SrB.sub.3C.sub.3 framework exhibits strong covalent bonding between sp.sup.3-hybridized B and C atoms, and weak interactions with the Sr guest. This strong sp.sup.3-hybridized covalent framework guarantees a high value for the bulk modulus (B0=249 (3) GPa) and incompressibility of SrB.sub.3C.sub.3 clathrate. Based on electron count, SrB.sub.3C.sub.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.sub.3B.sub.3].sup.3. SrB.sub.3C.sub.3 is one electron per formula unit short of this magic (insulator) electron count. Indeed, the band structure (FIG. 4B) shows thisa good gap for one electron more than SrB.sub.3C.sub.3. See also FIG. 4A.

[0087] SrB.sub.3C.sub.3 clathrate is likely the first member of a new class of strong and lightweight sp.sup.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.sup.3-bonded cages allows the carbon-based clathrates to possess diamond-like mechanical properties with a tunable electronic structure. While SrB.sub.3C.sub.3 is metallic, the electronic structure may be modified by substituting different guest atoms. For example, LaB.sub.3C.sub.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.sub.3C.sub.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.sub.3C.sub.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.

[0088] FIG. 5 provides an image showing in situ synthesis of SrB.sub.3C.sub.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.sub.3C.sub.3; the bottom portion is unconverted.

[0089] FIG. 6(A) and 6(B) demonstrate the superconducting transition temperature of SrB.sub.3C.sub.3. FIG. 6(A) shows phonon dispersion, and phonon density of states (PHDOS) and electron-phonon integral () of SrB.sub.3C.sub.3 at 0 GPa. FIG. 6(B) shows T.sub.c as a function of pressure. The insets show the evolution of and .sub.log with pressure. The superconducting transition temperature, T.sub.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.sub.c is 42 K at ambient pressure.

[0090] FIG. 7(A-D) show calculated stress-strain relations of SrB.sub.3C.sub.3. More specifically, FIGS. 7A and 7B show calculated tensile stress-strain relations for SrB.sub.3C.sub.3 with GGA (7A) and LDA (7B). FIGS. 7C and 7D show calculated shear stress-strain relations for SrB.sub.3C.sub.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.sub.3C.sub.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 <110> direction, which indicate that the <110> 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.sub.3C.sub.3. The lowest peak shear stress of 25 GPa (34 GPa with LDA) in the (110) [110] shear direction. These strength values place SrB.sub.3C.sub.3 as a very hard material with a hardness between 24-34 GPa.

[0091] 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.sub.3C.sub.3 is different in many respects from other carbon clathrate compounds discussed above.

[0092] The inventors report a carbon-based clathrate with composition 2La@B6C.sub.6 (LaB.sub.3C.sub.3). Like SrB.sub.3C.sub.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.sup.3+[B3C.sub.3].sup.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.

[0093] 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.sub.3C.sub.3, the inventors predicted that LaB.sub.3C.sub.3 becomes stable under high-pressure conditions. After conducting a computational survey of binary and ternary compounds in the LaBC system, the inventors estimate that LaB.sub.3C.sub.3 is the exothermic product of stable binary compounds and/or single components above about 20 GPa (FIG. 9A). Compared with calculations for SrB.sub.3C.sub.3, the pressure required for LaB.sub.3C.sub.3 is reduced by about 20 GPaa much more favorable synthetic condition that we tentatively attribute to the overall balanced electron count.

[0094] According to one aspect of the invention, LaB.sub.3C.sub.3 was synthesized according to the reaction LaB.sub.6+LaC.sub.2+4C.fwdarw.2LaB.sub.3C.sub.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. 9A), as expected for the clathrate (space group Pm3n). At 57 (2) GPa, the experimental lattice parameter a=4.66 is in good agreement with a DFT optimized LaB.sub.3C.sub.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. 9B). In addition, the stoichiometric conversion from the homogeneous precursor to a single phase strongly supports the 1La: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.sub.3C.sub.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.sub.3C.sub.3 synthesized near 50 GPa, this represents a significant pressure reduction and suggests promise for future stabilization efforts at even milder conditions.

[0095] Synchrotron X-ray diffraction patterns collected during decompression confirm the clathrate framework of LaB.sub.3C.sub.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 B0=254 (9) GPa and it's derivative B0=3.6 (4). The bulk modulus compares favorably with calculations based on DFT (FIG. 11C), noting that LDA significantly underestimates the volume for a given pressure. LaB.sub.3C.sub.3 is very incompressible with B0 in the range of technical ceramics used for armor, such as B4C (B0=243 (3) GPa), and the observed bulk modulus cannot be explained by any structures other than the clathrate. Similar to SrB.sub.3C.sub.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.

[0096] Like SrB.sub.3C.sub.3, LaB.sub.3C.sub.3 is predicted to crystallize in the cubic bipartite sodalite structure (FIG. 20). This structure type is known for the clathrate hydrate HPF.sub.6.Math.6H.sub.2O and is related to inorganic compounds such as BaPd.sub.2P.sub.4 through a small structural distortion, as well as superhydrides such as Ca/YH.sub.6 and CaYH.sub.12. The clathrate structure is comprised of a single truncated octahedral B.sub.12C.sub.12 cage with six four-sided faces and eight six-sided faces (4668). 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.

[0097] 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.sub.3C.sub.3, the observation of {120} reflections with significant intensity confirms the primitive space group Pm3n, and the ordered bipartite structural model produces the lowest reliability factor between experimental observations compared with other possible coloring schemes.

[0098] At 0 GPa, the LaB.sub.3C.sub.3 lattice is slightly expanded compared with that of SrB.sub.3C.sub.3 (a=4.934 vs. 4.868 ). This is somewhat surprising as the ionic radius of La.sup.3+ is expected to be smaller than that of Sr.sup.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.sub.3C.sub.3 (1.74 vs 1.72 SrB.sub.3C.sub.3) can be understood as additional charge transferred from guest to host. Indeed, charge density analysis (FIG. 21A and FIG. 21B) show that 1.82 e.sup. are transferred from La to the B/C host lattice, whereas 1.35 e.sup. are transferred for the case of Sr.

[0099] The electronic band structure of LaB.sub.3C.sub.3 also reveals the influence of the additional electron. Unlike SrB.sub.3C.sub.3, which is one electron short of a band gap, LaB.sub.3C.sub.3 shows an indirect gap of 1.29 eV in the .fwdarw.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.sub.3C.sub.3 is comparable to B-substituted silicon frameworks such as K.sub.7B.sub.7Si.sub.39 and LiBSi.sub.2 with electron precise frameworks in which the metals act as electron donors, as in Zintl phases.

[0100] Summarizing the above, a second carbon-boron clathrate in the bipartite sodalite structure was prepared with La guest atoms. Unlike metallic SrB.sub.3C.sub.3, LaB.sub.3C.sub.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., >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.

[0101] FIG. 8 provides a powder XRD pattern of LaB.sub.3C.sub.3 during decompression at 5 GPa (top), and after opening the cell in air (bottom). Bragg positions from cubic LaB.sub.3C.sub.3 are marked by blue lines. The diffraction intensity was limited by the opening angle above 2=14. When exposed to air, LaB.sub.3C.sub.3 degrades into unknown, presumably, hydrolysis products. Diffraction from unreacted LaB.sub.3C.sub.3 can still be observed in air several minutes after the cell was opened.

[0102] FIG. 10 provides phonon dispersion for LaB.sub.3C.sub.3 at 1 atmosphere, showing dynamic stability.

[0103] FIG. 11 provides stress-strain calculations for LaB.sub.3C.sub.3 with GGA and LDA. The calculations indicate a high ideal strength with hardness in near 30 GPa.

[0104] According to one aspect of the invention, the starting precursors to form the carbon clathrate compounds are borides, carbides and pure carbon.

[0105] 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.

[0106] SrB.sub.3C.sub.3 is predicted to be a hole metal with a high superconducting transition temperature near 42 K. LaB.sub.3C.sub.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.

[0107] According to one aspect embodiment of the invention, the carbon-based clathrate structures are composed of entirely sp.sup.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.

[0108] 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

[0109] 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

[0110] 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.sup.24p.sup.65s.sup.2, 2s.sup.22p.sup.1 and 2s.sup.22p.sup.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 600600600 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

[0111] 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.sub.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.sub.2 powder (1.400 g, 81% yield) isolated from the crucible under Ar. Powder XRD showed a small SrO impurity. SrB.sub.6 (EPSI Metals, 99.5%) and glassy carbon (Sigma Aldrich, 99.95%) were purchased commercially and used without further purification. Binary (SrB.sub.6+3SrC.sub.2) and ternary (SrB.sub.6+SrC.sub.2+4C) mixtures were prepared under an inert Ar atmosphere, sealed, then milled vigorously using Si.sub.3N.sub.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 m 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

[0112] 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.sub.2 and LaB.sub.6 standards with the DIOPTAS software. Pressure was calibrated using the equations of state of Ne and and/or Al.sub.2O.sub.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.

[0113] FIG. 12 provides raw XRD patterns before and after heating. The starting material is a mixture of SrB.sub.6 (Pm3m), SrC.sub.2 (l4/mmm) and glassy C in a 1:1:4 molar ratio, e.g., 2Sr:6C:6B. SrC.sub.2 and glassy carbon appear to be amorphous to X-rays and are not distinguishable from crystalline SrB.sub.6. After heating to a maximum temperature of 2900 K, the stoichiometric 1Sr:3C:3B mixture was converted to nearly phase pure SrB.sub.3C.sub.3 with trace residual SrB.sub.6.

[0114] FIG. 13 provides raw XRD patterns with increasing pressure and heating. The starting material is a mixture of SrB.sub.6 (Pm3m) and SrC.sub.2 (l4/mmm) in a 1:3 molar ratio. Above 14 GPa, SrC.sub.2 transforms to the R3m 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.sub.2 appears to become amorphous to X-rays; one broad feature is potentially attributable to the R3m phase. After heating at 150 GPa, many new Bragg peaks appear, which can be attributed to SrB.sub.3C.sub.3, SrBC and a high-pressure form of SrB.sub.6, as described below.

[0115] FIG. 14 demonstrates SrBC phase identification at 155 (6) GPa with Ne medium. The starting material is a mixture of SrB.sub.6 (Pm3m) and SrC.sub.2 (l4/mmm) in a 1:3 molar ratio. The data indicate the reaction SrB.sub.6+3SrC.sub.2.fwdarw.SrB.sub.3C.sub.3+3SrBC (and some unconverted SrB.sub.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.sub.3C.sub.3 (Pm3m). 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.sub.3C.sub.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.sub.3C.sub.3.

[0116] FIG. 15 shows SrBC phase identification at 142 (10) GPa with Al.sub.2O.sub.3 medium. The starting material is a mixture of SrB.sub.6 (Pm3m) and SrC.sub.2 (l4/mmm) in a 1:3 molar ratio. The data indicate the reaction SrB.sub.6+3SrC.sub.2.fwdarw.SrB.sub.3C.sub.3+3SrBC (and some unconverted SrB.sub.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.sub.3C.sub.3 (Pm3m). 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.sub.3C.sub.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.sub.3C.sub.3. Single-crystalline Al.sub.2O.sub.3 peaks were masked, but some powder intensity was unavoidable due to the breaking of crystals during compression to 140 GPa.

[0117] FIG. 16 shows stability of SrB.sub.3C.sub.3 and analysis of different possible stoichiometries. Relatedly, FIG. 1B shows predicted formation enthalpies of various SrBC 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.sub.xC.sub.y (excluding SrB.sub.3C.sub.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.sub.3C.sub.3. FIG. 16 shows the calculated EoS of the 600 SrB.sub.xC.sub.y compounds presented in FIG. 1B. The red line shows the EoS of SrB.sub.3C.sub.3 clathrate, which agrees with experiment. The grey lines represent the EoS of different SrB.sub.xC.sub.y compounds with calculated B0<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.sub.3C.sub.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.sub.3C.sub.3.

[0118] FIG. 17 shows X-ray diffraction of SrB.sub.3C.sub.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.sub.3C.sub.3 appears to degrade and is thus sensitive to moisture and/or oxygen.

[0119] FIG. 18 provides electronic band structure for SrB.sub.3C.sub.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.

[0120] FIG. 19A and FIG. 19B provide energetic stabilities as a function of pressure. FIG. 19A shows calculated enthalpies (H) per SrC.sub.2 unit as a function of pressure for various structures with respect to 14/mmm structure. FIG. 19B shows calculated enthalpies (H) per SrB.sub.6 unit as a function of pressure for various structures with respect to Pm3m structure.

[0121] FIG. 22 provides an EDX map of recovered sample showing homogeneous elemental distribution in heated region (2525 m2) with average composition La.sub.1.000.07B.sub.2.950.06C.sub.4.200.09O.sub.0.460.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.

[0122] Table 1 provides single-crystal diffraction and refinement parameters obtained from LaB.sub.3C.sub.3 crystal synthesized at 57 GPa and >3000 K. The analysis confirms the bipartite sodalite framework and primitive cubic space group.

TABLE-US-00001 TABLE 1 Structural Information for LaB.sub.3C.sub.3 from high- pressure single-crystal X-ray diffraction. Empirical formula LaB.sub.3C.sub.3 Formula weight 207.37 T/K 293(2) Crystal system cubic Space group Pm3n a/ 4.6712(16) V/.sup.3 101.93(10) Z 2 p.sub.calc/g .Math. cm.sup.3 6.757 /mm.sup.1 20.514 F(000) 180.0 Crystal size/mm.sup.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 l 3 Reflections collected 81 Independent reflections 21 [R.sub.int = 0.0452, R.sub.sigma = 0.0280] Data/restraints/parameters 21/0/4 Goodness-of-fit on F.sup.2 1.109 Final R indexes [I >= 2.sub. (I)] R.sub.1 = 0.0290, wR.sub.2 = 0.0633 Final R indexes [all data] R.sub.1 = 0.0297, wR.sub.2 = 0.0658 Largest diff. peak/hole/e.sup.3 +0.79/0.77

[0123] With respect to a clathrate compound with composition LaB.sub.3C.sub.3, work related to this compound was published in 2021 (Strobel, Angew. Chem. 2021,https://doi.org/10.1002/ange.202012821). LaB.sub.3C.sub.3 is a semiconductor, opposed to metallic/superconducting SrB.sub.3C.sub.3. Given that LaB.sub.3C.sub.3 is isostructural with SrB.sub.3C.sub.3, it is likely that many additional compositions with different guest atoms are possible, as well as different boron/carbon distribution schemes and entirely different clathrate structure types. The present application demonstrates that a large number of compositions are indeed possible, and further demonstrates the formation and properties of new structure types.

[0124] In addition to the SrB.sub.3C.sub.3 and LaB.sub.3C.sub.3 clathrates, the inventors also computationally predicted a related clathrate with composition ScB.sub.3C.sub.3. This clathrate has a different structure than Sr/LaB.sub.3C.sub.3 and crystallizes in the noncentrosymmetric space group Ama2. ScB.sub.3C.sub.3 is predicted to be a ferroelectric material with a critical temperature near 370 K. This work was published in 2020 (Zhu, Phys. Rev. Lett., 2020, https://doi.org/10.1103/PhysRevLett.125.127601).

[0125] Given the stability of SrB.sub.3C.sub.3 and LaB.sub.3C.sub.3, as demonstrated above and in the Prior Applications, the inventors performed extensive density functional theory (DFT) calculations to examine the energetic and dynamic stabilities of related compounds with various guest atom substitutions. The computational and experimental methods used to investigate these compounds were similar to those described in the Prior Applications, except where otherwise disclosed in this application.

[0126] After structural optimization at ambient pressure, it was determined that BaB.sub.3C.sub.3, Y B3C.sub.3, CeB.sub.3C.sub.3, PrB.sub.3C.sub.3, NdB.sub.3C.sub.3, PmB.sub.3C.sub.3, SmB.sub.3C.sub.33, EuB.sub.3C.sub.3, ErB.sub.3C.sub.3, TmB.sub.3C.sub.3, YbB.sub.3C.sub.3, and LuB.sub.3C.sub.3 are all dynamically stable as shown in FIG. 23. The compounds show negative formation enthalpies at pressures as low as 8 GPa as shown in FIG. 24. The phonon dispersion relations and formation enthalpies were calculated with methods similar to those above (e.g. in FIG. 6A) and were recently published in the supporting material of Geng et al. J. Am. Chem. Soc, 2023, https://doi.org/10.1021/jacs.2c10089. These results strongly support the claim that Type-VII B-C clathrates in the bipartite sodalite structure (cubic MB.sub.3C.sub.3) are possible for many different guest metals, and may provide a variety of useful compounds with varying and tunable properties. The inventors have recently also synthesized a Type-VII CaB.sub.3C.sub.3 clathrate-like compound. Thus, Type-VII clathrates represent a large class of new materials with highly tunable properties depending on the specific guest atoms.

[0127] Experimentally, a new MB.sub.3C.sub.3 clathrate was verified for Nd guest atoms, showing that the theorized compounds can be prepared and analyzed in the physical world. Single-crystal refinement parameters for NdB.sub.3C.sub.3 are shown in FIG. 25A, with the spin-polarized density of states shown in FIG. 25B. Single-crystal X-ray diffraction measurements show the cubic bipartite sodalite lattice shown in FIG. 25C. Nd atoms within the structure were calculated to exhibit ferromagnetic ordering, indicating that Type-VII B-C clathrates may also serve as novel magnetic materials.

[0128] These findings clearly support the claim that MB.sub.3C.sub.3 clathrates represent a broad class of novel materials with tunable properties. Depending on the metal guest atom these materials may be metals/superconductors, semiconductors and/or ferroelectric/ferromagnetic materials.

[0129] In addition to MB.sub.3C.sub.3 clathrates with one type of guest atom in the cages, the inventors have demonstrated that different metal atoms may occupy the cages in the form of binary guest clathrates, i.e. M1M2B6C.sub.6. Starting from precursor mixtures of metal borides and metal carbides with different metals, i.e., M1B6+M2C.sub.2 the high-pressure synthetic process results in the formation of binary guest clathrates. FIG. 26A shows powder X-ray diffraction patterns for one such compound, namely CaSrB.sub.6C.sub.6 (Ca.sub.0.5Sr.sub.0.5B.sub.3C.sub.3) obtained after high-pressure laser heating of CaB.sub.6+SrC.sub.2+4 C. FIG. 26B shows powder X-ray diffraction patterns for another such compound, namely NaSrB.sub.6C.sub.6 (Na.sub.xSr.sub.(1-x)B.sub.3C.sub.3) obtained after high-pressure laser heating of 2SrB.sub.6+Na.sub.2C.sub.2+10 C. Based on these results, a large, combinatorial number of M1.sub.xM2.sub.yB.sub.6C.sub.6 (or similar) compounds may be produced following this approach. The below discussion of FIGS. 27A, 27B, and 27C demonstrates a technique for the generation and examination of such compounds.

[0130] As mentioned above, MB.sub.3C.sub.3 compounds may serve as moderately high-temperature superconductors depending on the specific guest composition. In 2017, the inventors predicted that SrB.sub.3C.sub.3 exhibits superconductivity at 1 atm with T.sub.c42 K (u*=0.1). (Li Zhu, Hanyu Liu, R. E. Cohen, Roald Hoffmann, Timothy A. Strobel, Prediction of a thermodynamically stable carbon-based clathrate, arXiv: 1708.03483, https://doi.org/10.48550/arXiv.1708.03483) In June 2017 the inventors also predicted that BaB.sub.3C.sub.3 exhibits superconductivity with T.sub.c45 K (*=0.1). Recently, the inventors have validated the predictions for SrB.sub.3C.sub.3 experimentally and have shown T.sub.c22 K at 23 GPa. The inventors developed a novel experimental technique that allows for the measurement of in situ electrical transport in high-pressure, high-temperature diamond anvil cells (DAC) with in situ synthesis at 50 GPa and temperatures >3000 K. This was achieved by producing an insulating DAC gasket equipped with Pt electrical contacts that rest on top of a sapphire crystal that thermally insulates the diamond anvils during high-temperature laser heating.

[0131] FIG. 26A shows a sample geometry for in situ electrical transport with HPHT synthesis of SrB.sub.3C.sub.3. The 1Sr:3B:3C precursor is pressed between diamond anvils, separated by sapphire single crystals that serve as thermal insulation. Four Pt electrical probes sit on the bottom sapphire crystal and are secured within the insulating gasket. FIG. 27B shows photomicrographs of the sample chamber before loading (with view through the top diamond), and after loading sample (view through the bottom diamond). FIG. 27C shows representative XRD patterns before and after heating. The heated region shows 60 wt. % SrB.sub.3C.sub.3 with unconverted cubic SrB.sub.6. FIG. 27D shows an X-ray radiograph of the sample chamber with four Pt wires. Outline of the sapphire is indicated by the white dashed line. The dashed black square contains the laser-heated region between the Pt contacts. The zoom window shows an x-ray-diffraction intensity map for SrB.sub.3C.sub.3 and Pt showing the synthesized sample between the contacts.

[0132] After in situ HPHT synthesis, SrB.sub.3C.sub.3 samples within DACs were placed within open-flow cryostats to measure temperature-dependent electrical transport. Due to resistive precursor present at the sample contacts, samples initially show semiconductor behavior upon cooling. After cooling below 20 K, samples show sharp, non-hysteretic, drops in resistance, consistent with a superconducting transition under percolative transport as shown in FIG. 28A. The onset temperature of the resistance drop follows the calculated pressure dependence for T. of SrB.sub.3C.sub.3 with a Coulomb pseudopotential (*) value of 0.15 as shown in FIG. 28B.

[0133] Superconductivity was verified by measuring electrical transport under applied magnetic fields using a non-magnetic DAC. With increasing magnetic field strength, the onset temperature of the resistance drop is suppressed up to a critical field H.sub.c2 (T), where superconductivity is destroyed. Resistance at a pressure of 40 GPa as a function of temperature and magnetic field is shown in FIG. 28C. Complementary differential resistance measurements also confirm superconductivity, showing a peak a zero bias when measured below T.sub.c and H.sub.c2. FIG. 28D shows .sub.0H.sub.c2 (T) for the transition onset (black points) compared with the WHH model (=0). A linear fit is shown to guide the eye. The experimental work demonstrating these results was recently published in Phys. Rev. Research 5, 013012 (2023) https://journals.aps.org/prresearch/abstract/10.1103/PhysRevResearch.5.013012 (also posted on arXiv at https://arxiv.org/ftp/arxiv/papers/1708/1708.03483.pdf).

[0134] In addition to simple (single metal guest) Type-VII clathrates such as SrB.sub.3C.sub.3 and BaB.sub.3C.sub.3, the inventors predict superconductivity in binary Type-VII clathrates. In this case, different ratios of metal atoms containing different valences can be used to precisely tune the density of states at the Fermi level and significantly enhance T.sub.c. For example, in September 2018, the inventors predicted that RbSrB.sub.6C.sub.6 (FIG. 29A, based on a 221 supercell of SrB.sub.3C.sub.3) is dynamically stable at 1 atm and possessed T.sub.c of 50 K, which is enhanced by hole doping, shifting the Fermi level closer to a peak in the electronic density of states. As above, this can also be denoted as Rb.sub.0.5Sr.sub.0.5B.sub.3C.sub.3. The phonon dispersion of RbSrB.sub.6C.sub.6 at 1 atm showing dynamic stability is shown in FIG. 29B.

[0135] Following the inventors' predictions on hole-doped, binary Type-VII clathrates and positive synthesis experiments (See FIG. 26), the inventors performed electrical transport measurements on Na.sub.xSr.sub.(1-x)B.sub.6C.sub.6 clathrate samples produced at 60 GPa and 3000 K. Upon cooling, the inventors observed a sharp drop in the electrical resistance similar to SrB.sub.3C.sub.3, but the onset temperature is shifted nearly +20 K as shown in FIG. 30A. This result suggests an enhancement in the superconducting transition temperature and is consistent with theoretical predictions. As shown in FIG. 30B, the T.sub.c vs P curve of an Na-doped compound is similar to that of the compound SrB.sub.3C.sub.3, only shifted by nearly +20K.

[0136] The inventors recently published predictions of high-T.sub.c superconductors for other binary M1M2B.sub.6C.sub.6 clathrates including novel compositions such at KPbB.sub.6C.sub.6 (Geng et al. J. Am. Chem. Soc, 2023, https://doi.org/10.1021/jacs.2c10089). Over the past year other groups have also made predictions of high-T.sub.c, binary Type-VII clathrates based on the inventors' initial work on SrB.sub.3C.sub.3 (See: Ting-Ting Gai, Peng-Jie Guo, Huan-Cheng Yang, Yan Gao, Miao Gao, and Zhong-Yi Lu, Van Hove Singularity induced phonon-mediated superconductivity above 77 K in hole-doped SrB.sub.3C.sub.3, Phys. Rev. B 105, 224514, Published 27 Jun. 2022.

[0137] https://journals.aps.org/prb/abstract/10.1103/PhysRevB.105.224514. See also: Simone Di Cataldo, Shadi Qulaghasi, Giovanni B. Bachelet, and Lilia Boeri, High-Tc superconductivity in doped boron-carbon clathrates, Phys. Rev. B 105, 064516, Published 25 Feb. 2022, https://journals.aps.org/prb/abstract/10.1103/PhysRevB.105.064516).

[0138] According to one aspect of the invention, the inventors have predicted and synthesized the carbon-based clathrate-like structure with Ca as the guest atom. The clathrate-like compound h-CaB.sub.3C.sub.3 is different in many respects from other carbon clathrate compounds discussed above. Note that the inventors have also produced the cubic bipartite sodalite clathrate with the same compositions, but distinct structure.

[0139] The inventors report having synthesized a new compound with hexagonal symmetry in the space group P-62m and composition CaB.sub.3C.sub.3. While the composition is the same as cubic Type-VII clathrate, the structure type is entirely new as shown in FIG. 31A. Hexagonal CaB.sub.3C.sub.3 is clathrate-like but does not contain closed polyhedral cages. It contains star-like channels down the crystallographic c-axis, which potentially indicate mobility of Ca.sup.2+ ions for possible divalent battery applications as shown in FIG. 31B. It is also metallic with appreciable density of states at the Fermi level, indicating potential for superconductivity. As shown by the experimental data in FIG. 31C, with channels down the c-axis, the electronic band structure and density of states show that the compound is metallic.

[0140] 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: