arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 9 Jun 2006

Transcription

1 First-principles study on the intermediate compounds of LiBH 4 Nobuko Ohba, Kazutoshi Miwa, Masakazu Aoki, Tatsuo Noritake, and Shin-ichi Towata Toyota Central Research & Development Laboratories, Inc., Nagakute, Aichi , Japan arxiv:cond-mat/66228v1 [cond-mat.mtrl-sci] 9 Jun 26 Yuko Nakamori, and Shin-ichi Orimo Institute for Materials Research, Tohoku University, Sendai , Japan Andreas Züttel Physics Department, University of Fribourg, Perolles, CH-17 Fribourg, Switzerland (Dated: February 6, 28) We report the results of the first-principles calculation on the intermediate compounds of LiBH 4. The stability of LiB 3H 8 and Li 2B nh n(n = 5 12) has been examined with the ultrasoft pseudopotential method based on the density functional theory. Theoretical prediction has suggested that monoclinic Li 2B 12H 12 is the most stable among the candidate materials. We propose the following hydriding/dehydriding process of LiBH 4 via this intermediate compound : LiBH Li2B12H LiH H2 LiH + B +3 H2. The hydrogen content and enthalpy of the first reaction are estimated 2 to be 1 mass% and 56 kj/mol H 2, respectively, and those of the second reaction are 4 mass% and 125 kj/mol H 2. They are in good agreement with experimental results of the thermal desorption spectra of LiBH 4. Our calculation has predicted that the bending modes for the Γ-phonon frequencies of monoclinic Li 2B 12H 12 are lower than that of LiBH 4, while stretching modes are higher. These results are very useful for the experimental search and identification of possible intermediate compounds. I. INTRODUCTION Hydrogen is the most promising environmentally clean energy carrier to replace fossil fuels. The use of hydrogenbased energy in practical applications such as fuel cell vehicles, however, requires the development of safe and efficient hydrogen storage technology. Complex hydrides including the light metal lithium (Li) have sufficient gravimetric hydrogen storage capacity, and many researches and development about the lithium complex hydride such as LiBH 4 1,2,3 and LiNH 2 4,5,6 have been done recently. Particularly, LiBH 4 can desorb about 14 mass% of hydrogen by the following thermal decomposition : LiBH 4 LiH + B H 2. (1) However, the experimental value of enthalpy of this reaction is 69 kj/mol H 2 1 and LiBH 4 is too stable to release hydrogen at ambient condition. High temperature and high pressure are needed for rehydriding reaction, so its reversibility becomes a problem for practical use. Züttel et al. 1 investigated the hydrogen desorption from LiBH 4 in details, and reported that the thermal desorption spectra of LiBH 4 mixed with SiO 2 -powder exhibited three hydrogen desorption peaks. These peaks were observed around 5 K, 55 K and 6 K, and the products corresponded to LiBH 3.6, LiBH 3 and LiBH 2, respectively. The compounds in quotes were nominal compositions estimated from amount of desorbed hydrogen. This is a strong evidence for the existence of the intermediate compounds, and the hydrogen desorption reaction takes place in at least two steps not a single step like in Eq.(1). Although low temperature release of hydrogen and improvement of reversibility can be expected by use of the intermediate compound, little attention has been paid yet. In this study, the stability of the intermediate compounds of LiBH 4 has been investigated theoretically and we clarify a hydriding/dehydriding process of LiBH 4. The thermal desorption experiment observed a structural transition around 38 K and the melting of the compound at 55 K. The stability of the solid phases at absolute zero temperature is mainly discussed here. In Sec. II, we describe the details of the computational method and the possible intermediate compounds of LiBH 4. Section III reports the calculated results on the stability of some candidate compounds. And then, we discuss the hydriding/dehydriding reaction via the intermediate compounds. Furthermore, the electronic structures and Γ-phonon frequencies on the most stable compound are studied. We also investigate the stability of the borane complex anions, and then compare with the solid compounds. II. COMPUTATIONAL DETAILS First-principles calculations have been performed by the ultrasoft pseudopotential method 7 based on the density functional theory. 8 The generalized gradient approximation (GGA) formula 9 is applied to the exchangecorrelation energy. The interaction between the ion cores and electrons is described by the ultrasoft pseudopotential 7, and the norm-conservation constraint 1 is imposed on Li for the calculation efficiency improvement. The scalarrelativistic all-electron calculations are first carried out on the free atoms(ions). We chose 2s and 2p states for

2 2 both Li and B pseudopotential as the reference states with the cutoff radii of 2 a.u. (Li) and 1.5 a.u. (B), respectively. A single projector function is used for each angular momentum component. The 3d state is treated as the local part of pseudopotential. The hydrogen pseudopotential is constructed from 1s, 2s and 2p states with the cutoff radii of 1.1 a.u.(s) and 1.2 a.u. (p). We use double projector functions for the s component and a single projector function for the p component. For all pseudopotentials, the pseudo-wave functions and the pseudo-charge-augmentation functions are optimized by a method similar to that proposed by Rappe et al. 11 Also, the partial core correction 12 is taken into account for Li and B pseudopotentials. In the solid-state calculations, the pseudo-wave functions are expanded by plane waves with a cutoff energy equal to 15 hartrees. The cutoff energy for the charge density and potential is set to be 12 hartrees. The integral over the Brillouin zone is approximated by the summation on the k-grid elements of which the edge lengths closer to the target value of.15å 1 as possible. We confirmed that these calculation conditions gave a good convergence of energy within.2 ev/atom. The preconditioned conjugate-gradient technique is employed to minimize the Kohn-Sham energy functional. A procedure based on the iterative diagonalization scheme 13 and the Broyden charge mixing method 14 is adopted in this study. Optimization of crystal structures is performed till the atomic forces and the macroscopic stresses become less than hartree/bohr and.1 GPa, respectively. During the structural optimization process, the partial occupation numbers near the Fermi level are determined by the Fermi-Dirac distribution function with k B T = hartrees. The Helmholtz free-energy functional, 15 including the entropy term, is minimized instead of the Kohn-Sham energy functional. Then, the improved tetrahedron method 16 is used in order to minimize the Kohn-Sham energy functional in the optimized structure. The dynamical matrix is calculated by the force-constant method 17 to obtain the Γ phonon frequencies. The atomic displacement is set to be.2 Å. The further details of calculation are described in Refs. 2, 13 and references therein. As the candidates of the intermediate compounds, the existing alkali metal-b-h materials are used. For example, it is well known that boron (B) and hydrogen (H) form inorganic compounds called borane and the compounds of CsB 3 H 8, 18 K 2 B 6 H 6, 19 and K 2 B 12 H 12 2 are reported. The space group of crystal structure for the compound CsB 3 H 8 is Ama2 (No.4) and cation Cs + and anion [B 3 H 8 ] are arranged same as the NaCl-type structure with an orthorhombic distortion. The crystal structures of K 2 B 6 H 6 and K 2 B 12 H 12 are best described as an anti-caf 2 -type arrangement with K + cation in the center of a tetrahedron formed by [B n H n ] 2 dianions. We first assumed that LiB 3 H 8, Li 2 B 6 H 6 and Li 2 B 12 H 12 compounds, which had the same crystal structures as existing CsB 3 H 8, K 2 B 6 H 6 and K 2 B 12 H 12, respectively. The stability of these candidates was evaluated using the firstprinciples calculations. Since the existence of a series of closo-type dianions [B n H n ] 2 (n = 5 12) is also well known, our calculation have been expanded to Li 2 B n H n (n = 5 12) compounds. Although Li atom and [B n H n ] cluster are arranged in the anti-caf 2 -type structure and the unit cell parameters are supposed to be a facecentered-cubic symmetry (a = b = c, α = β = δ = 9 ) as starting points for the structural optimization, the output relaxed compounds have different symmetry depending on the structure of [B n H n ] clusters. The alkali metal salt of monoclinic K 2 B 12 (OH) with closo- [B 12 (OH) 12 ] 2 dianions has been reported, and monoclinic Li 2 B 12 H 12 with similar crystal structure is also examined. III. RESULTS AND DISCUSSIONS A. Stability of candidate compounds : LiB xh y Table I shows the calculated results on the structural parameters and the cohesive energies of LiB 3 H 8 and Li 2 B n H n (n = 5 12). We denote the cubic Li 2 B 12 H 12 based on K 2 B 12 H 12 with type-1 and monoclinic Li 2 B 12 H 12 based on K 2 B 12 (OH) 12 with type- 2, respectively. Compared with the cohesive energy of two type of Li 2 B 12 H 12, the value of monoclinic Li 2 B 12 H 12 (type-2) is larger than that of type-1. Therefore, the type-2 Li 2 B 12 H 12 is easy to form. After the next paragraph, the only result concerning monoclinic Li 2 B 12 H 12 (type-2) is shown. The enthalpy of formation for hydriding reactions including LiBH 4 from LiH (NaCl-type), α-b (rhombohedral), and H 2 molecule are given in Table II, where the zero-point energy corrections are not taken into consideration. They are provided using calculated cohesive energies of ev/atom for LiH, ev/atom for α-b, ev/atom for H 2 molecule, and ev/atom for LiBH 4 (orthorhombic Pnma symmetry), respectively. The enthalpies of formation for Li 2 B n H n (n = 1, 11, 12) are more negative than that for LiBH 4. Therefore these compounds have a great potential for generating in hydriding reactions from LiH and α-b as the intermediate phase of LiBH 4. In order to understand the stability of each compound intuitively, we represent the enthalpy of formation corresponding to the following reaction with hydrogen of δ mole in Fig. 1 : (1 δ)[lih+b]+δh 2 1 δ x LiB xh y + (1 δ)(x 1) LiH, x (2) where δ = (y 1)/(2x + y 1) shows the mole fraction of H 2 and x 1. From Fig. 1, monoclinic Li 2 B 12 H 12 is formed as the intermediate compound of the hydriding/dehydriding reaction of LiBH 4, because the actual reaction goes along the lowest state of enthalpy. Therefore, following hydriding/dehydriding process is proposed

3 3 TABLE I: Structural parameters and cohesive energies (E coh ) of the candidate compounds : LiB xh y. It is denoted that the cubic Li 2B 12H 12 based on K 2B 12H 12 with type-1 and monoclinic Li 2B 12H 12 based on K 2B 12(OH) 12 with type-2, respectively. Compound Space group Unit cell parameters Atomic position E coh Atom Site x/a y/b z/c (ev/atom) LiB 3H 8 Ama2 a = Å Li 4b (No.4) b = Å B1 4b c = Å B2 8c H1 4b H2 4b H3 8c H4 8c H5 8c Li 2B 5H 5 R3m a = Å Li1 3a (No.16) c = Å Li2 3a.7513 (hexagonal axes) B1 3a.876 B2 3a.9295 B3 9b H1 3a.1595 H2 3a.8586 H3 9b Li 2B 6H 6 Fm 3m a = Å Li 8c (No.225) B 24e.1553 H 24e.354 Li 2B 7H 7 I2 a = Å Li1 2b (No.5) b = Å Li2 2b c = Å B1 2a.1564 β = B2 4c B3 4c B4 4c H1 2a.2792 H2 4c H3 4c H4 4c Li 2B 8H 8 I 42m a = Å Li 4d (No.121) c = Å B1 8i B2 8i H1 8i H2 8i Li 2B 9H 9 R3 a = 7.44 Å Li1 3a (No.146) c = Å Li2 3a.87 (hexagonal axes) B1 9b B2 9b B3 9b H1 9b H2 9b H3 9b Li 2B 1H 1 I422 a = Å Li 4d (No.97) c = Å B1 4e.1827 B2 16k H1 4e.2984 H2 16k Li 2B 11H 11 I2 a = 7.83 Å Li1 2b (No.5) b = Å Li2 2b c = 7.8 Å B1 2a.1264 β = 7.52 B2 4c B3 4c B4 4c B5 4c B6 4c H1 2a.2355 H2 4c H3 4c H4 4c H5 4c H6 4c Li 2B 12H 12 Fm 3 a = 1.83 Å Li 8c (type-1) (No.22) B 48h H 48h Li 2B 12H 12 P2 1/n a = Å Li 4e (type-2) (No.14) b = Å B1 4e c = Å B2 4e β = B3 4e B4 4e B5 4e B6 4e H1 4e H2 4e H3 4e H4 4e H5 4e H6 4e

4 4 TABLE II: The enthalpies of formation for the hydriding reaction of various compounds LiB xh y including LiBH 4 from LiH, α-b, and H 2 molecule, where the zero-point energy corrections are not taken into consideration. Hydriding reaction Enthalpy of formation (kj/mol H 2) LiH + 3B + 7/2H 2 LiB 3H LiH + 5B + 3/2H 2 Li 2B 5H LiH + 6B + 2 H 2 Li 2B 6H LiH + 7B + 5/2H 2 Li 2B 7H LiH + 8B + 3 H 2 Li 2B 8H LiH + 9B + 7/2H 2 Li 2B 9H LiH + 1B + 4 H 2 Li 2B 1H LiH + 11B + 9/2H 2 Li 2B 11H LiH + 12B + 5 H 2 Li 2B 12H LiH + B + 3/2H 2 LiBH 4 75 Enthalpy [kj] LiH + B δ H Li 2 B 5 H 5 Li 2 B 8 H 8 Li 2 B 6 H 6Li2 Li 2 B 9 H 9 B 7 H 7 Li 2 B 1 H 1 Li 2 B 11 H 11 Li 2 B 12 H 12 LiB 3 H 8 LiBH 4 FIG. 1: Enthalpy of formation relative to (1 δ)[lih + B] and molecular H 2 of δ mole (reaction Eq.(2) in text), as a function of the mole fraction of H 2, δ. : LiBH Li 2B 12 H LiH H 2, (3) LiH + B H 2. (4) The enthalpy without zero-point energy effects and hydrogen content of the first reaction (Eq.(3)) are 56 kj/mol H 2 and 1 mass% and those of the second reaction (Eq.(4)) are 125 kj/mol H 2 and 4 mass%, respectively. These results agree with the experimental ones, which reported that 9 mass% of hydrogen are liberated in the hydrogen desorption peak in thermal desorption spectra of LiBH 4 mixed with SiO 2 -powder in Ref. 1. The enthalpy value of 56 kj/mol H 2 for the first reaction is lower than computational one of 75 kj/mol H 2 for direct dehydriding reaction of LiBH 4 (Eq.(1)). 2 So the hydrogen desorption and absorption can occur at low temperature and pressure comparatively by using this intermediate compound. We also calculated other possible crystal structures such as Cu 2 B 1 H 1 -type 22 and [Li(thp) 3 ] 2 [B 11 H 11 ]- type. 23 Moreover, we performed the calculation for the Γ-phonon frequencies of the candidate materials. If there was a soft-mode, we lowered crystal structure symmetry with moving the atoms along the direction of the soft-mode eigenvectors. The maximum gain of cohesive energy was.3 ev/atom for Li 2 B 6 H 6, however, no compound or crystal structure more stable than monoclinic Li 2 B 12 H 12 was found. The present calculation result on the energy of monoclinic Li 2 B 12 H 12 provides the upper limit value for the thermodynamic stability. In conclusion, the existence of the intermediate compound of LiBH 4 was predicted theoretically. Recently Kang et al. 24 have reported that LiBH and LiB are the intermediate phases of LiBH 4, and propose the dehydriding reaction of LiBH 4 through LiBH as follows : LiBH 4 LiBH H 2. (5) We also performed the first-principles calculation of the orthorhombic phase of Pnma LiBH. The reaction enthalpy for Eq.(5) is 1.3eV per LiBH 4 formula unit (84kJ/molH 2 ), which is in good agreement with reported value of 1.28eV per LiBH 4 formula unit in Ref.24. On the other hand, the reaction enthalpy of Eq.(3) is.63ev per LiBH 4 formula unit. Therefore, the dehydriding reaction via the intermediate compound Li 2 B 12 H 12 is energetically more preferable one. B. Intermediate compound : Li 2B 12H 12 Here we describe the fundamental properties, such as the crystal structure, the electronic structure, and the Γ phonon frequencies, on Li 2 B 12 H 12 which is expected as intermediate compound of LiBH 4. Optimized crystal structure model of monoclinic Li 2 B 12 H 12 is shown in Fig. 2. The bond lengths between B and H atom of Li 2 B 12 H 12 are Å, and they are very close to those for LiBH 4 ( Å) reported in Ref. 2. As for the B-B bond lengths of intra-icosahedron, the values of Å for Li 2 B 12 H 12 are comparable to the experimental ones for the α-rhombohedral boron (α-b) of Å, 25 too. Since boron crystal has the icosahedral B 12 cluster as a common structural component, the decomposition of [B 12 H 12 ] 2 anion into B 12 cluster and hydrogen molecule is easy to understand. Figure 3 shows the total and partial density of states (DOS) for Li 2 B 12 H 12, which denotes that it has the energy gap of 5.6 ev. Since there is little contribution of Li orbital for occupied states, Li 2 B 12 H 12 consists of Li + and [B 12 H 12 ] 2 ions. The orbitals of B and H hybridize each other and the feature of occupied states in DOS is analogous with the distribution of the computed energies of the bonding molecular orbitals in [B 12 H 12 ] 2. 26

5 5 1 TOTAL 5 FIG. 2: (Color online) Crystal structure model of monoclinic Li 2B 12H 12(type-2). Large, middle, and small spheres denote Li, B, and H atoms, respectively. DOS(states/eV unit-cell) 1 1 Li B s p s p The Γ-phonon mode frequencies of monoclinic Li 2 B 12 H 12 have been calculated so that the vibrational properties can be compared with experiments easily and directly. The phonon density of states is shown in Fig. 4. It is divided into three regions in the same case as LiBH 4. The first region is less than 3 cm 1 where the displacements of Li atoms are dominant. The second region is between 45 cm 1 and 11 cm 1 where the B-H bond of [B 12 H 12 ] 2 dianions vibrates with changing in angle between them (bending modes), and the third region is between 24 cm 1 and 26 cm 1 where the interatomic distance of B-H bond is changing along bond axis (stretching modes). The frequencies of bending modes are lower than LiBH 4 described in Ref. 2, while stretching modes are higher. The investigation using in situ Raman spectroscopy is effective for the confirmation of the short-range order or bonding of LiBH We examine the atomistic vibrations of LiBH 4 during heating by in situ Raman spectroscopy, and the identification of spectra modes originating from Li 2 B 12 H 12 is now in progress. 28 C. Stability of complex anions : [B nh n] 2 Finally, we consider the relation between the stability of the compounds Li 2 B n H n and those of complex anions [B n H n ] 2. The energies of the isolated [B n H n ] 2 are obtained using a face-centered-cubic supercell with a = 2 Å. The single Γ point is used for the k-point sampling. The energies for the charged systems are computed by adding uniform background charges and improved with Makov and Payne correction 29 for the interaction between periodic image charges. Figure 5 shows the comparison of the formation energies, E solid f and E complex f, corresponding to the following 5 H s p Energy(eV) FIG. 3: (Color online) The total and partial densities of states (DOS) for Li 2B 12H 12. The energy is measured in ev relative to the top of valence states. PDOS (arb. unit) Frequency (cm -1 ) FIG. 4: Phonon density of states for monoclinic Li 2B 12H 12. The contribution of the TO Γ-phonon modes indicated by vertical bars is only taken into account and the Gaussian broadening with a width of 3 cm 1 is used. reactions: and 2LiH + nb + n 2 H 2 Li 2 B n H n, (6) 2 nb + n 2 H 2 + 2e [B n H n ] 2, (7) where the energies are normalized by the number of B-H pairs for comparison purposes. We can find a fairly good

6 6 Ef complex (ev/b-h pair) n=1 n=7 n=11 n=6 n=9 n=12 (type-2) n=12 (type-1) n=8 n= E f solid (ev/b-h pair) FIG. 5: Comparison of the formation energies, Ef solid and E complex f, corresponding to the reactions Eq.(6) and Eq.(7) in text, respectively, where the energies are normalized by the number of B-H pairs. correlation between both energies. This is probably due to the fact that the electrostatic interaction between Li + and [B n H n ] 2 is not sensitive to n and a large part of this energy cancels out with that of LiH in Eq. (6). Among the closo-type dianions considered here, [B 12 H 12 ] 2 is the most stable one. LiBH 4 desorbs hydrogen at temperatures above the melting point. The latent heat of fusion has been reported 3 to be.78 ev/formula-unit for LiBH 4 and similar values are expected for Li 2 B n H n, which is considerably smaller than the energy difference between [B 12 H 12 ] 2 and other close-type dianions. The stability of the [B 12 H 12 ] 2 anion supports our main conclusion in this study, that is, the intermediate phase suggested by the experiment is Li 2 B 12 H 12. We have examined the stability of LiB 3 H 8 and Li 2 B n H n (n = 5 12), which are the possible intermediate compounds of LiBH 4, by the first-principles calculation. Our computational results for enthalpy of the hydriding reactions provide that monoclinic Li 2 B 12 H 12 is the most stable one among the candidates. The following hydriding/dehydriding process of LiBH 4 is proposed : LiBH 4 1/12Li 2 B 12 H /6LiH +13/12H 2 LiH + B + 3/2H 2. The hydrogen content of the first and the second reaction are 1 mass% and 4 mass%, respectively, which agree well with the thermal desorption spectra (TDS) experiment on LiBH 4. 1 The heat of formation without zero point energy corrections for the first reaction, which is estimated from the solid state LiBH 4 with Pnma symmetry (not liquid LiBH 4 ), is 56 kj/molh 2. This value is lower than that for the direct reaction (LiBH 4 LiH + B H 2). Therefore, low temperature release of hydrogen can be expected by use of this intermediate compound. We have calculated the electronic structure and the Γ- phonon frequencies of monoclinic Li 2 B 12 H 12. This compound has the energy gap of 5.6 ev and consists of Li + and [B 12 H 12 ] 2 ions. From the phonon density of states, it is predicted that its bending modes have lower frequencies than that of LiBH 4, 2 while stretching modes are higher. The identification of the experimental Raman spectra modes originating from Li 2 B 12 H 12 is now in progress. The stability of closo borane complex anions [B n H n ] 2 was also examined. We found a fairly good correlation between the formation energies of the solid phases and the isolated dianions. This result supports the validity of the intermediate compound indicated in the TDS experiment being Li 2 B 12 H 12. There are various kinds of borane, such as nido-type [B n H n+4 ] and arachno-type [B n H n+6 ], in addition to closo-type borane [B n H n+2 ]. 31 The nido and arachno borane are derived from closo borane by removing one and two vertices, respectively. These generate the salts in the ionized state as well as closo-type dianions. These alkali metal salts are also the candidates of the intermediate compound of LiBH 4. We will study the details of the hydriding/dehydriding process for LiBH 4 including the stability of these materials in future. IV. SUMMARY Acknowledgments We wish to thank M. Matsumoto, R. Jinnouchi and S. Hyodo for helpful discussions. This work was partially supported by the New Energy and Industrial Technology Development Organization(NEDO), International Joint Research under the Development for Safe Utilization and Infrastructure of Hydrogen Project (24-25). 1 A. Züttel, P. Wenger, S. Rentsch, P. Sudan, Ph. Mauron, and Ch. Emmenegger, J. Power Sources 118, 1 (23). 2 K. Miwa, N. Ohba, S. Towata, Y. Nakamori, and S. Orimo, Phys. Rev. B 69, (24). 3 S. Orimo, Y. Nakamori, and A. Züttel, Mater. Sci. Eng. B 18, 51 (24). 4 P. Chen, Z. Xiong, J. Luo, J. Lin, and K. L. Tan, Nature 42,32 (22). 5 Y. Nakamori and S. Orimo, Mater. Sci. Eng. B 18, 48 (24). 6 K. Miwa, N. Ohba, S. Towata, Y. Nakamori, and S. Orimo, Phys. Rev. B 71, (25). 7 D. Vanderbilt, Phys. Rev. B 41, R7892 (199); K. Laasonen, A. Pasquarello, R. Car, C. Lee, and D. Vanderbilt,

7 7 ibid. 47, (1993). 8 P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964); W. Kohn and L. J. Sham, Phys. Rev. 14, A1133 (1965). 9 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996); 78, 1396 (E) (1997). 1 D. R. Hamann, M. Schlüter, and C. Chiang, Phys. Rev. Lett. 43, 1494 (1979). 11 A. M. Rappe, K. M. Rabe, E. Kaxiras, and J. D. Joannopoulos, Phys. Rev. B 41, 1227 (199). 12 S. G. Louie, S. Froyen, and M. L. Cohen, Phys. Rev. B 26, 1738 (1982). 13 A. Fukumoto and K. Miwa, Phys. Rev. B 55, (1997). 14 V. Eyret, J. Comput. Phys. 124, 271 (1996). 15 N. D. Mermin, Phys. Rev. 137, A1441 (1965). 16 P. E. Blochl, O. Jepsen, and O. K. Andersen, Phys. Rev. B 49, (1994). 17 K. Kunc and R. M. Martin, Phys. Rev. Lett. 48, 46 (1982). 18 H. J. Deiseroth, O. Sommer, H. Binder, K. Wolfer, and B. Frei, Z. anorg. Allg. Chem. 571, 21 (1989). 19 I.Y. Kuznetsov, D.M. Vinitskii, K.A. Solntsev, N.T. Kuznetsov, and L.A. Butman, Russian J. Inorg. Chem. 32, 183 (1987). 2 I. Tiritiris and T. Schleid, Z. anorg. Allg. Chem. 629, 139 (23). 21 T. Peymann, C. B. Knobler, S. I. Khan, and M. F. Hawthorne, J. Am. Chem. Soc. 123, 2182 (21). 22 R. D. Dobrott and W. N. Lipscomb, J. Chem. Phys. 37, 1779 (1962). 23 O. Volkov, W. Dirk, U. Englert, and P. Paetzold, Z. Anorg. Allg. Chem. 625, 1193 (1999). 24 J. K. Kang, S. Y. Kim, Y. S. Han, R. P. Muller, and W. A. Goddard III, Appl. Phys. Lett. 87, (25). 25 B. Morosin, A. W. Mullendore, D. Emin, and G. A. Slack, in Boron-Rich Solids, edited by D. Emin, T. L. Aselage, C. L. Beckel, I. A. Howard, and C. Wood, AIP Conf. Proc. No. 14 (AIP, New York, 1986), p R. B. King, I. Silaghi-Dumitrescu, and A. Lupan, Inorg. Chem. 44, 7819 (25). 27 S. Orimo, Y. Nakamori, G. Kitahara, K. Miwa, N. Ohba, S. Towata, and A. Züttel, J. Alloys Comp , 427 (25). 28 S. Orimo, Y. Nakamori, N. Ohba, K. Miwa, M. Aoki, S. Towata, and A. Züttel, Appl. Phys. Lett., to be published. 29 G. Makov and M. C. Payne, Phys. Rev. B 51, 414 (1995). 3 M. B. Smith and G. E. Bass Jr., J. Chem. Eng. Data 8, 342 (1963). 31 G. A. Olah, K. Wade, and R. E. Williams, ELEC- TRON DEFICIENT BORON AND CARBON CLUS- TERS, (JOHN WILEY & SONS, 1991).