Hydrogen Storage in Boron Nitride and Carbon Nanomaterials

Boron nitride (BN) nanomaterials were synthesized from LaB6 and Pd/boron powder, and the hydrogen storage was investigated by differential thermogravimetric analysis, which showed possibility of hydrogen storage of 1–3 wt%. The hydrogen gas storage in BN and carbon (C) clusters was also investigated by molecular orbital calculations, which indicated possible hydrogen storage of 6.5 and 4.9 wt%, respectively. Chemisorption calculation was also carried out for B24N24 cluster with changing endohedral elements in BN cluster to compare the bonding energy at nitrogen and boron, which showed that Li is a suitable element for hydrogenation to the BN cluster. The BN cluster materials would store H2 molecule easier than carbon fullerene materials, and its stability for high temperature would be good. Molecular dynamics calculations showed that a H2 molecule remains stable in a C60 cage at 298 K and 0.1 MPa, and that pressures over 5 MPa are needed to store H2 molecules in the C60 cage.


Introduction
Hydrogen is a carrier with high energy density, and forms only water and heat.On the other hand, fossil fuels generate toxic fuels, such as COx, NOx and SOx.Therefore, clean hydrogen energy is expected as substitute of fossil fuel in the 21st century, and gas storage ability more than 6.5 wt% is needed for car application according to the US Department of Energy.Although LaNi5H6 is already used as an H2 gas storage materials, the ability is only 1 wt% because of the large atomic number of La and Ni [1].On the other hand, fullerene-like materials, which consist of light elements, such as boron, carbon and nitrogen, would store more H2 gas compared to the metal hydrides.Various works have been reported on hydrogen storage ability of carbon nanotubes, fullerenes and nanomaterials [1].It was reported that multi-walled carbon nanotubes could absorb hydrogen from 1 wt% up to 4.6 wt%.These results indicate carbon nanotubes might be a possible candidate of hydrogen storage materials although the evaluation of hydrogen storage measurements is necessary.Many studies on H2 gas storage in carbon (C) nanomaterials have been reported [2][3][4][5][6][7].Boron nitride (BN) nanomaterials are also expected in prospective application because they provide good stability at high temperatures with high electronic insulation in air [8].Hydrogen storage in BN nanomaterials has also been studied recently [9][10][11][12][13][14][15][16].It is difficult to absorb hydrogen by physisorption, both inside nanotubes and at the interstitial channels of nanotubes at room temperature, which is due to too weak bonding of physisorbed hydrogen to reach large storage at room temperature.Therefore, it is considered that chemisorption is a necessary requirement for hydrogen storage at room temperature.It was reported that chemisorption ratio observed at room temperature increases with increasing alkali/carbon rate.Therefore, if the energy of chemisorption can be lowered, hydrogen storage ability of carbon and BN nanotubes would be improved.
The purpose of the present study is twofold.The first is to prepare BN nanotubes, nanocapsules and nanocages, and to investigate hydrogen gas storage by thermogravimetry/differential thermogravimetric analysis (TG/DTA) and first principle calculation.LaB6 and Pd were selected in order to take advantage of their catalytic effect to produce the BN nanomaterials.La has shown excellent catalytic properties for producing a large number of single-walled carbon nanotubes and enlarging their diameter, and Pd is also expected to act as a hydrogen storage material.Although gas storage of hydrogen and argon in carbon nanotubes has been reported, carbon nanotubes are oxidized at 600 °C in air.On the other hand, BN starts to be oxidized into boron oxide and nitrogen ~900 °C in air, which indicates excellent heat resistance compared to carbon materials for gas storage.The differences between BN and C nanomaterials are summarized as shown in Table 1.To understand the formation of BN nanostructures, high-resolution electron microscopy (HREM) were carried out, which are very powerful methods for atomic structure analysis [17,18].For the BN nanomaterials, hydrogen gas storage measurements were carried out using TG/DTA.The second purpose of the present work is to investigate H2 gas storage ability of BN and C fullerene-like materials by theoretical calculation.Although huge amount of calculation is required to calculate nanotubes, it is considered that H2 molecules enter from the cap of nanotubes.Energy of chemisorption and stable hydrogen position inside the clusters were investigated as cap structures of nanotubes.Influence of endohedral element on H2 gas storage ability of BN fullerene-like materials was also investigated by theoretical calculations.Energy of chemisorption was investigated for cage clusters as a cap structure of nanotubes, and B24N24 clusters, reported by mass spectrometry [8,19], were selected for the storage material.Li, K and Na elements were also selected for the doping atoms because these elements are of high electropositive character.The present study will give us a guideline for hydrogenation of BN nanomaterials as hydrogen storage materials.

Experimental Procedures and Calculation Methods
Mixture powder compacts made of boron particles (99%, 40 mm, Nilaco Corp., Tokyo, Japan), LaB6 particles (99%, 1 mm, Wako Pure Chemical Industries Ltd., Osaka, Japan) and Pd particles (99.5%, 0.1 mm, Nilaco Corp., Tokyo, Japan), with the size of 3 mm in height and 30 mm in diameter were produced by pressing powder at 10 MPa.The atomic ratios of metal (M) to boron (B) were in the range of 1:1-1:10.The green compacts were set on a copper mold in an electric-arc furnace, which was evacuated down to 2.0 × 10 −3 Pa.After introducing a mixed gas of Ar 0.025 MPa and N2 0.025 MPa, arc-melting was applied to the samples at an accelerating voltage of 200 V and an arc current of 125 A for 2 s.Arc-melting was performed with a vacuum arc-melting furnace (NEV-AD03, Nissin Giken, Saitama, Japan), and the white/gray BN nanomaterial powders were collected from surface of the bulk.Samples for HREM observation were prepared by dispersing the materials on holy carbon grids.HREM observation was performed with a 300 kV electron microscope (JEM-3000F, JEOL, Tokyo, Japan).To confirm the formation of BN fullerene materials, energy dispersive x-ray spectroscopy analysis was performed using the EDAX system with a probe size of ~10 nm.In order to measure hydrogen gas storage in BN nanomaterials, BN nanomaterials were extracted from the obtained powder by a supersonic dispersing method based on the Stokes equation using ethanol.The Stokes equation is expressed as follows: where η is the viscosity of liquid; σ is the density of particles; ρ is the density of liquid; g is the gravitational acceleration; d is the diameter of particles; t is the subsidence time; and h is the height of liquid.Since there is a big difference in the size and density of the produced BN nanocapsules/nanotubes and other powders, it is believed that this method would be effective for separation of BN nanomaterials [20].After separation of BN nanomaterials, hydrogen storage was measured by TG/DTA at temperatures in the range of 20-300 °C in H2 atmosphere [21].
To search the optimized structure of B99N99 and B36N36 with H2 molecules, semi-empirical molecular orbital calculations (parameterized model revision 3, PM3) were performed.The energies of B36N36 with H2 molecules were calculated by first principle single point energy calculation using Gaussian.In the calculation, 3-21G was used as ground function with Hartree-Fock level.B24N24, B36N36, B60N60 and C60 were selected for cluster calculations.To investigate the optimized structures, semi-empirical molecular orbital calculations (parameterized model revision 5, PM5) were performed by using MOPAC.Eigenvector following (EF) method was used for geometry optimization, and charge is 0. The default self-consistent method was restricted Hartree-Fock (HF).Multi-electron configuration interaction (MECI) was used in order to prevent repulsion between electrons becoming excessive.Chemisorption calculations of hydrogen atoms were performed for B24N24, B36N36 and C60 by PM5 calculations [22,23].
M@B24N24 (M@: element encaged in B24N24) was selected for cluster calculations.To investigate the optimized structures and the chemisorption calculations of hydrogen atoms, molecular orbital calculations were performed by using PM5 in MOPAC.The EF method was used for geometry optimization.The default self-consistent method was restricted as HF, and the MECI was used.In the calculation, chemisorption of hydrogen was performed on boron and nitrogen positions outside the cage [24].Energy levels and densities of states (DOS) for B24N24 and Li@B24N24 were also calculated by the first principles calculation with discrete variational (DV)-Xα method.
Molecular dynamics calculations were carried out to confirm the stability of H2 molecules into C60 at 298 K and 0.1 MPa [25].Number, temperature, and pressure (NTP) ensembles were used in the calculation as follows: number of atoms (N), temperature (T) and pressure (P) are constant.Conditions of H2 gas storage were also calculated.Number, temperature, and enthalpy (NPH) ensembles were used in the calculation as follows: number of atoms (N), pressure (P) and enthalpy (H) are constant.These molecular dynamics were calculated by organic potential.

Synthesis of Boron Nitride Nanomaterials
For some metals, formation enthalpies with boron (H for B) and nitrogen (H for N) are indicated in Figure 1a,b, respectively.The data were from theoretical calculations [8].Difference of formation enthalpy (H for N-H for B) is also shown in Figure 1c.The difference of formation enthalpy (H for N-H for B) is very important for the formation of BN fullerene nanomaterials, because reactivity with nitrogen and boron is decided by this enthalpy.Basically, BN nanotubes are formed when rare earth metals are used as catalytic metals, such as Y, Zr, Nb, Hf, Ta, W and La.These elements have minus enthalpy, as shown in Figure 1c.It means that catalytic elements for synthesis of BN nanotubes should be selected from those with minus formation enthalpy (H for N-H for B).From the present guideline, Sc element could be a good catalytic element to form BN nanotubes.The detailed formation mechanism of the nanostructures with different metals should be investigated further.
A HREM image of a BN nanotube produced using LaB6 powder is shown in Figure 2a.In Figure 2a, the diameter of the five-layered BN nanotube is changing from bottom to top, and amorphous patches are observed mostly at the top.Low magnification images of BN nanostructures were reported [17,20].A HREM image of BN nanocage produced from LaB6/B powder is shown in Figure 2b, which indicates square-like shape, and four-membered rings of BN exist at the corner of the cage.The BN nanocage has a network-like structure, whose atomic arrangement is basically consistent with the B36N36 cluster structure [8].BN nanocapsules with Pd nanoparticles were also produced as shown in Figure 2c, and Pd nanoparticles were covered by a few BN layers.In order to confirm the formation of BN nanocapsules, EDX analysis was carried out, which showed the atomic ratio of B:N ≈ 1.

Hydrogen Storage in Boron Nitride Nanomaterials Studied by Thermogravimetry/Differential Thermogravimetric Analysis
DTA and TG curve of BN nanomaterials produced from LaB6 powder is shown in Figure 3.At a temperature around 70 °C, an increase of sample weight of 0.3 mg is observed.Weight change for this sample was almost reversible, which indicates the reversibility of hydrogen adsorption.It also suggests that the hydrogen atoms would be physically absorbed.For the samples of La:B = 1:6 and Pd:B = 1:4, weight increases of 3.2% and 1.6% were observed, respectively, as listed in Table 2.  BN nanostructures would have energy barriers for H2 molecules to pass through tetragonal and hexagonal rings.Figure 4 are structure models in which H2 molecule passes from hexagonal rings of B99N99 and B36N36.Single point energies were calculated with changing set point of H2 molecule from the center of the cage at intervals of 0.1 nm.There is energy barrier that is given for H2 molecules to pass through the hexagonal rings.The energy barrier of B36N36 hexagonal rings showed the smallest value of 14 eV in the present calculation, which is smaller than that of 27 eV at tetragonal rings.The DE of C60 hexagonal rings was also calculated to be 16 eV for comparison.This value is higher than that of B36N36 hexagonal rings, and the H2 molecule would pass from hexagonal rings of B36N36 easier than from hexagonal rings of C60.It is known that H2 molecules are adsorbed on walls of single-walled carbon nanotubes over 7 MPa as an experimental result.As a result of the comparison, H2 molecules enter into B36N36 from hexagonal rings easier than tetragonal rings B36N36 and hexagonal rings of C60.A formation mechanism of BN nanotubes and nanocapsules synthesized in the present work is described below.Metal and boron particles are melted by arc-melting, and during the solidification of the liquid into metal and/or boride nanoparticles, excess boron would react with nitrogen to form BN layers at the surface of the nanoparticles.Because of electrical insulation, BN fullerene materials are usually fabricated by arc-discharge method with specific conducting electrodes such as HfB2 and ZrB2.The present arc-melting method from mixed powder has two advantages for BN nanomaterial production.Since the powder becomes conducting by pressing, special electrodes are not needed.In Although gas storages of hydrogen and argon in carbon nanotubes have been reported, there are few reports for gas storage in BN fullerene materials and for calculations [21].Weight increase of the sample in TG measurements was observed as shown in Figure 3.It might be due to the hydrogen gas storage in the BN nanomaterials.Since there would be metal and boron nanoparticles in the separated BN nanomaterials even after the separation, further qualification and evaluation of the samples are needed for hydrogen storage.
Carbon fullerenes and boron nitride fullerenes are sublimed at 600 and 1000 °C, respectively.Boron nitride fullerenes would storage H2 molecules with smaller energy than carbon fullerenes, and would give good stability at high temperature.Boron nitride fullerene materials would be a better candidate for H2 storage materials.

Molecular Orbital Calculations of Hydrogen Storage in Boron Nitride and C Clusters
Figure 5 is a structural model of hydrogen atoms chemisorbed on boron and nitrogen for BN clusters and carbon clusters.Atoms bonded with hydrogen are moved outside from the clusters.Energies for hydrogen chemisorption on each position are summarized as Table 3. Hydrogen bonding with nitrogen is more stable than that with boron, and hydrogen bonding with carbon is more stable than C60.This result indicates the chemical modification of carbon fullerenes.When two hydrogen atoms were chemisorbed on carbon clusters, energies of carbon clusters increased.3. Chemisorption energy of hydrogen (H) atoms on C60 and B24N24.ΔE = (Heat of formation after hydrogen addition) − (Heat of formation before hydrogen addition).In the calculation, chemisorption of hydrogen was performed on outside of cage, and on boron, nitrogen and carbon of tetragonal and hexagonal rings.The B36N36 cluster has tetragonal and hexagonal BN rings, and there are four kinds of boron and nitrogen positions for the hydrogen chemisorption, as shown in Figure 6.For the stability calculations of H2 molecules in these clusters, 30 H2 molecules were also introduced in B24N24, B36N36, B60N60 and C60. Figure 5 is a structural model of hydrogen atoms chemisorbed on boron and nitrogen for B36N36.Energies for hydrogen chemisorption on each position are summarized as Table 3. Hydrogen bonding with nitrogen is more stable than that with boron because nitrogen atoms are more electrophilic compared to boron atoms.In addition, hydrogen bonding on tetragonal ring is more stable than that of hexagonal ring.Chemisoption of hydrogen with C60 reduced the energy.When two hydrogen atoms were chemisorbed on carbon clusters, energies of carbon clusters increased.To investigate stability of H2 molecules in clusters, energies were calculated for H2 molecules introduced in the center of clusters.The structural models are shown in Figure 7. Energies of C60, B24N24, B36N36 and B60N60 were calculated to be 0.58, 20.71, 20.93 and 20.82 eV/mol atom, respectively.This result indicates that B24N24, B36N36 and B60N60 with H2 molecules are more stable than C60 with H2 molecule, and that B36N36 is the most stable in BN clusters.Hydrogen storage (wt%) of B60N60 cluster was calculated, and structural models are shown in Figure 7.Other C60, B24N24 and B36N36 were calculated, as summarized in Table 4. From Table 4, C and BN cluster showed possibility of hydrogen storage of ~6.5 and ~4.9 wt%, respectively.

Cluster Number of H atoms Additional position of H Heat of formation (eV) ΔE (eV) Before addition After addition
Although hydrogen storage (wt%) of C60 is better than those of BN clusters, energy increase by hydrogen addition is higher for C60 (4.6-5.2 eV/H atom) compared to BN clusters (1.0-3.1 eV/H atom), which would be due to C-H interaction.This indicates that needed energy for hydrogen storage in BN clusters is lower compared to the C60.From Table 4, stability of H2 molecules in B24N24 and B36N36 seems to be higher than that of C60.Hydrogen atoms were calculated to be chemisorbed inside C60, as shown in Table 4.This indicates that inside of C60 also has good reactivity as well as outside of the cage.Chemisorption inside the C60 clusters may indicate that hydrogen capacity may gradually reduce under the adsorption-desorption cycles.On the other hand, BN clusters have no chemisorption inside the clusters, which indicates that the BN clusters would be a better candidate for stable adsorption-desorption cycles.C60 cluster shows the minimum energies in spite of the positive values.It is believed that p-electrons outside and inside of the cage would increase the energy.However, s-electrons in the 5-and 6-membered carbon rings would stabilize the cage structure.Hydrogen addition to the C60 decreased the energy, which agrees with the above description.The hydrogen storage of the present BN cluster was calculated to be ~4.9 wt%, and the measured values were 1-3 wt%.This might be due to the remained catalytic metals in the sample, which reduces the storage ratio.Although carbon nanotubes are oxidized at ~600 °C, BN nanomaterials are almost stable ~900 °C in air, which indicates BN fullerenes have higher thermal and chemical stability than those of carbon fullerenes.BN fullerenes with good thermal and chemical stability can store H2 molecules with less energy, and they have the same chemisorption energy and higher stability, compared to carbon clusters.BN fullerene materials would be better candidates for H2 storage materials.

Effects of Endohedral Element in B24N24 Clusters on Hydrogenation
Figure 8 is a structural model of hydrogen atoms chemisorbed on nitrogen position for M@B24N24.Energies for hydrogen chemisorption on boron and nitrogen positions are summarized as Table 5. Heats of formation (eV) of M@B24N24 by hydrogenation are indicated in Figure 9.In Table 5 and Figure 9, "None" means that metal catalyst was not encaged in B24N24, and "BH" and "NH" means that hydrogen atom was chemisorbed on boron and nitrogen, respectively.Figure 9 shows that hydrogen bonding with nitrogen is more stable than that with boron because nitrogen atoms are more electrophilic compared to boron atoms.Figure 9 also indicates that energies of chemisorption on M@B24N24 are much lower than that of B24N24.From this result, Li atom works as a good endohedral element for hydrogen chemisorption.Metal catalysts in the present work have been reported to generate hydrides such as LiH because of strong interaction between hydrogen and metal atoms.In the present work, it is clarified that Li doping and nitrogen positions are suitable for hydrogenation for the B24N24 clusters.Bond lengths of B-H and N-H for M@B24N24 are summarized in Table 6.B-H and N-H distance decreased by doping element in the BN cluster.Bond lengths of N-H are shorter than that of B-H for these BN clusters, and the bond length of N-H for Li@B24N24 is the shortest.The endohedral atoms appear to decrease the repulsion energy between the electrons of the hydrogen atom and the p-electrons of B24N24.Figure 8d is a structural model of M@B24N24H24, which indicates hydrogenation on all nitrogen positions for M@B24N24, and the hydrogen storage capacity is summarized in Table 7. Li is also a good element for hydrogen storage capacity because Li is the lightest element of these.Energy levels and density of states for B24N24 and Li@B24N24 are shown in Figure 10.Fermi levels in energy level diagrams and DOS diagrams correspond to 0 eV.B24N24 and Li@B24N24 show energy gaps of 4.8873 and 0.0247 eV between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), respectively.This means that electron of Li element transferred to the B24N24 cage, and electronic state of BN cluster would change from semiconductor to metallic property by Li doping in BN clusters.Although it was reported that BN nanotubes were produced by lithium vapor, synthesis of Li@BN has not been reported yet.BN fullerenes have high thermal and chemical stability, and M@BN fullerenes have lower energy of chemisorption compared to the present work.Since the BN clusters were reported to be doped with metal elements, M@BN clusters would be produced.BN fullerene materials with endohedral element such as Li would be better candidates for H2 storage materials.

Molecular Dynamics Calculations of Hydrogen Storage in C60 Clusters
C60 included a H2 molecule kept in stable state at T = 298 K and P = 0.1 MPa.This unit cell is shown in Figure 11.This model was calculated with N = 62 atoms, T = 298 K and P = 0.1 MPa.Although the H2 molecule vibrated in the C60 cage, it was not discharged from the cage.Some H2 molecules were stored in the C60 cage when the pressure was 5 MPa.This unit cell is shown in Figure 12, which is composed of 32 C60 (fcc) with 288 H2 molecules (fcc).This model was calculated with N = 2496 atoms and P = 5 MPa.H2 molecules pass through the hexagonal rings of the C60 cage at 0.5 ps.After introducing H2 molecules into the C60 cage at 2.5 ps, they are stored and stable in C60. Figure 13 shows a schematic model of a single H2 molecule stored in a hexagonal ring of C60.These results indicate that fullerene-like materials can store H2 gas in cage at T = 298 K and P = 0.1 MPa, and some H2 molecules were stored in the C60 cage when the pressure was greater than 5 MPa.It is known that H2 molecules are adsorbed on the walls of single-walled carbon nanotubes over 7 MPa as an experimental result [26].Since a C60 cluster has a large curvature and H2 molecules are very small compared to C60, it is considered that adsorption and storage of H2 gas occur at the same time.Detailed theoretical calculation on hydrogen storage in BN nanomaterials has also been performed [27][28][29][30][31].

Conclusions
HREM observation showed the formation of BN nanotubes, nanocapsules and nanocages, which were synthesized from mixtures of LaB6, Pd, and boron powder by using an arc melting method.Although samples produced with Pd include only BN nanocapsule structures, samples produced with LaB6 present BN nanocapsule, nanotube and nanocage structures.After separation of BN nanomaterials using ethanol, hydrogen storage was measured by TG/DTA, and the BN nanomaterials produced from LaB6 and Pd/boron powder showed possibility of hydrogen storage of ~3 wt%.Energies for BN and C clusters with hydrogen were investigated by molecular orbital calculations.Stabilities of these clusters were estimated from the energies, and possibility of H2 storage ability was predicted by these results.B24N24, B36N36, B60N60 and C60 clusters were selected as the tip structure of the nanotubes.Chemisorption calculation of hydrogen for BN clusters showed that hydrogen bonding with nitrogen atoms was more stable than that with boron atoms, and that hydrogen bonding with tetragonal ring is more stable than that with hexagonal ring.Energy increase by hydrogen addition to C60 is higher compared to BN clusters because of C-H interaction, which indicates that the BN clusters have higher stability with hydrogen atoms.BN cluster and C cluster showed possibility of H2 storage of ~4.9 wt% and ~6.5 wt%, respectively.H2 storage of C cluster is better than those of BN clusters.However, stability of H2 molecules in BN clusters might be higher than that of C clusters.Energies of hydrogenation for B24N24 were also investigated by molecular orbital calculations.Chemisorption calculation of hydrogen in the B24N24 clusters showed that hydrogen bonding with nitrogen atoms was more stable than that with boron atoms.Chemisorption calculations also indicate that endohedral elements decreased energies of hydrogenation and Li atom is suitable element for hydrogen chemisorption.Molecular dynamics calculations showed that a single H2 molecule remains in a stable state in a C60 cage at T = 298 K and P = 0.1 MPa.It is confirmed that pressure of over 5 MPa is required to store H2 molecules in a C60 cage.As a result of SPE calculations, H2 molecules enter from hexagonal rings of fullerene-like materials.The energy required of H2 discharge from fullerene-like materials is similar to that of H2 storage.The present study indicates that BN fullerene materials could be a one of the possible candidates as hydrogen gas storage materials.

Figure 1 .
Figure 1.(a) Formation enthalpy with boron (H for B) and (b) nitrogen (H for N); (c) difference of formation enthalpy (H for N-H for B).

Figure 4 .
Figure 4. Structure models that H2 molecule passes through hexagonal BN rings of (a,b) tip of B99N99 nanotube and (c) B36N36 cluster.
, commercial arc-melting furnaces can be used.These advantages indicate a simpler fabrication method compared to the ordinary arc-discharge methods.

Figure 6 .
Figure 6.Structure models for hydrogen chemisorption on B36N36: (a) structural model of B36N36; (b,c) structural models of hydrogen chemisorption on N and B positions of tetragonal BN ring, respectively; and (d,e) structural models of hydrogen chemisorbed on N and B positions of hexagonal BN ring, respectively.

Figure 9 .
Figure 9. Heats of formation of hydrogenation for M@B24N24 clusters by endohedral elements.

Figure 10 .
Figure 10.Energy levels diagrams and density of states of (a) B24N24 and (b) Li@B24N24.

Figure 11 .
Figure 11.Molecular dynamics calculation to confirm stability of H2 molecule into C60 at 298 K and 105 Pa.NTP ensembles and organic potential were used in the calculation.

Figure 12 .
Figure 12.Molecular dynamics calculation to find condition of H2 gas storage.Unit cell is composed from 32 C60 (fcc) with 288 H2 molecules (fcc).NPH ensembles and organic potential were used in the calculation.

Figure 13 .
Figure 13.Schematic model of H2 molecule stored in hexagonal ring in the order of (a-d).

Table 2 .
Atomic ratio of starting materials, produced structures and hydrogen storage.

Table 4 .
Energy of clusters with hydrogen.
* A C-C bond was broken.

Table 5 .
Chemisorption energies of a hydrogen atom on M@B24N24.

Table 6 .
Bond length of B-H and N-H on M@B24N24.

Table 7 .
Hydrogen storage capacity of chemisorption on nitrogen position of M@B24N24.