1 Introduction
At present, hydrogen has been widely regarded as a promising carbon-free renewable energy carrier [
1-
3]. Immense efforts have been made worldwide to find feasible hydrogen storage materials with high gravimetric and volumetric densities. To get hydrogen uptake and release at ambient conditions [
4,
5], the binding energies of H
2 are estimated as 0.2–0.4 eV/H
2 [
6], which is intermediate between physisorption and chemisorption. Carbon nanotubes have ever attracted much attention as candidates for high-capacity hydrogen storage materials, due to the hollow tubular structure and large surface area. However, according to abundant investigations, pure carbon nanotubes exhibit low hydrogen storage potential [
7,
8]. Metal-decorated carbon nanostructures have been proposed to achieve prospective hydrogen storage capacity [
9-
11].
Boron nitride (BN) nanotubes [
12-
18] have many analogical aspects with carbon nanotubes. In addition, in terms of the chemical and thermal stability, BN nanostructures may be stable lightweight hydrogen storage media. However, theoretical studies also disclose that pristine BN nanomaterials cannot serve as practical hydrogen storage materials at ambient conditions in view of the low hydrogen adsorption energy and gravimetric density [
19,
20]. It was reported that BN nanotubes can take up 1.8–2.6 wt% hydrogen under
10 MPa at room temperature [
21] and that the collapsed BN nanotubes prepared through metal-catalyzed growth mechanism exhibited a higher hydrogen storage capacity up to 4.2 wt% at
10 MPa [
22]. The large hydrogen uptake in BN nanotubes was attributed primarily to chemisorption. The collapsed surfaces of BN nanotubes effectively increase the hydrogen adsorption capacity due to the increase of dangling bonds and specific surface area [
23]. The enhanced hydrogen physisorption and chemsorption may also be attributed to interstitial regions in realistic nanotube bundles [
24,
25]. Another reasonable possibility is that there remain metal cluster contaminants in BN nanotubes [
26-
28], and the enhanced hydrogen dissociation and adsorption may result from the strong interaction between metal atoms and hydrogen molecules. However, metal atoms tend to form clusters on the surface of nanomaterials due to the larger cohesive energy, and clustering reduces the hydrogen storage capacity and changes the hydrogen adsorption state [
29,
30]. To avoid the clustering problem, the interaction between metal atoms and the substrates should be enhanced by introducing chemically active sites; therefore, vacancy defects are effective to increase the metal binding characteristics.
Recently, much experimental progress has been made in the preparation of BN nanosheets [
31-
33], and extensive computations have been performed to explore their physical and chemical properties and potential applications [
34-
37]. As the analogue of graphene, BN monolayer nanosheets may also be applicable to hydrogen storage media. In this work, we performed detailed density functional theory (DFT) computations to gain an insight into the hydrogen adsorption in BN nanosheets.
2 Computational methods
All computations were carried out within the DFT framework by using the Vienna
ab initio simulation package (VASP) [
38]. The electron exchange correlation is treated within the generalized gradient approximation (GGA) with PW91 functional [
39]. Also, local density approximation (LDA) with CA functional [
40] was employed to evaluate the hydrogen molecule adsorption energies and compare with the GGA results, as previous studies showed that GGA usually underestimates the H
2 adsorption energy, whereas LDA overestimates the interaction [
24,
41-
43], and the real value may be located in the range between GGA and LDA results; therefore, GGA and LDA results are widely adopted to evaluate the range of the hydrogen adsorption energies before better methodologies appear. Meanwhile, the electron–ion interactions are modeled by the ultrasoft pseudopotentials [
44]. The cutoff energy for the plane wave-basis expansion was chosen to be 360 eV, and the convergence threshold was set as 10
-4 eV in energy and 10
-3 eV/Å in force. As for simulating two-dimensional BN sheets, a 4 × 4 supercell was used. The Monkhorst–Pack
k point 4 × 4 × 1 was adopted for the system. The positions of all the atoms in the supercell were fully relaxed during the geometry optimization.
The hydrogen adsorption energy, , is defined as , where , , and is the energy of the complex of BN nanosheets and hydrogen molecules, the separate BN nanosheets, and H2 molecule, respectively. n is the number of hydrogen molecules absorbed.
3 Results and discussion
3.1 Ca and Sc adsorption on pristine BN nanosheets
Since H
2 is weakly adsorbed onto BN nanosheet [
19], we attempted to introduce Ca and Sc atoms to BN nanosheets in order to increase the hydrogen adsorption energy (Fig. 1). The computed binding energies of Ca and Sc atom on pristine BN nanosheets, and the cohesive energies of bulk Ca and Sc metals are summarized in Table 1. After full relaxation, the distance between the Ca atom and the B or N atoms of the hexagon ring in pristine BN nanosheets is about 4.13 Å. The binding energy of Ca on pristine BN nanosheets is only 0.15 eV, which is even lower than that of Ca on pristine graphene [
45]. The cohesive energy of the Ca bulk metal is 1.94 eV/Ca atom; therefore, Ca atoms definitely tend to form clusters on the surface of pristine BN nanosheets. Subsequently, the binding state of Sc atoms on pristine BN nanosheets was also examined. The distance between the Sc atom and the B or N atoms of the hexagon ring is about 2.68 Å in the optimized structure, and the binding energy of Sc on pristine BN nanosheets is 1.84 eV, which is much larger than that of Ca on BN nanosheets, but it is still lower than the cohesive energy of the Sc bulk metal (4.13 eV/Sc atom). Similar to Ca atoms, it is inevitable for Sc atoms to agglomerate on the surface of pristine BN nanosheets. To avoid the clustering problem, the binding states between metal atoms and BN nanosheets should be enhanced by introducing defects.
3.2 Ca and Sc atom adsorption on BN nanosheets with vacancies
Next, we consider BN nanosheets with various vacancies, so we remove one boron or nitrogen atom from the BN nanosheet to create a vacancy defect. Figure 2 is a schematic representation of the boron vacancy (
and nitrogen vacancy (
defect in BN nanosheets. After full relaxation, the distances between the three nitrogen atoms around the boron vacancy is 2.69 Å, which is larger than that in the pristine BN nanostructures (2.51 Å). The results are consistent with previous experimental and theoretical reports that the N–N distances are
2.70 Å surrounding boron vacancy and
2.50 Å in perfect BN nanostructures [
46,
47]. There is a contraction in the distance between the boron atoms around the nitrogen vacancy (
as compared with the pristine BN nanosheets, and the distance is about 2.35 Å, which is shorter than that in the pristine BN nanosheets (2.52 Å). The structural modification is important for the structural stability of B vacancy. Owning to the larger distance and the stronger repulsion of nitrogen atoms around boron vacancy, it is difficult to form chemical bonds among the nitrogen atoms around boron vacancy and get a pentagonal ring.
When we attached one Ca atom onto the nitrogen vacancy of BN nanosheet, as presented in Fig. 3(a), the optimized distance between the Ca atom and the three boron atoms around the nitrogen vacancy is 3.29 Å, and the binding energy between Ca and BN nanosheet with is 0.95 eV, which is greatly enhanced as compared with that of pristine BN nanosheet. The increased binding energy is due to the larger charge transfer between metal atoms and BN substrate. The Hirshfeld charge analysis indicates that Ca carries a 0.34 e positive charge. However, the binding energy is still lower than the cohesive energy of the bulk Ca metal (1.94 eV), and Ca atoms persist in clustering on BN nanosheets with nitrogen vacancies. Similar results appeared when we introduced one Sc atom onto the nitrogen vacancy of BN nanosheet [Fig. 3(b)]. After full relaxation, the distance between the Sc atom and the three boron atoms around the nitrogen vacancy is 2.35 Å, and the binding energy is 2.35 eV, which is also lower than the cohesive energy of Sc bulk metal (4.13 eV). The Hirshfeld charge analysis indicates that 0.43 e charge is transferred from Sc atom to BN nanosheet with nitrogen vacancy. The smaller binding energy is unfavorable to exclude the formation of the clusters.
Then, we checked the effects of boron vacancies on the binding energies of metal atoms on BN nanosheets. When we attached one Ca atom onto the boron vacancy of BN nanosheet, as shown in Fig. 3(c), the boron vacancy retains its structure, and no apparent distortion occurs. The effect of Ca on the vacancy structure is inconsistent with the case of graphene. Adding one Ca atom to the vacancy of graphene results in the formation of a 5–9 carbon ring [
45]. After full relaxation, the distance between the Ca atom and the three nitrogen atoms around the boron vacancy is 2.17 Å. As listed in Table 1, the binding energy between Ca and the BN nanosheet with boron vacancy is 6.99 eV, which is increased so significantly as to exceed the cohesive energy of Ca bulk metal (1.94 eV). Ca atoms strongly bind to BN nanosheets due to charge transfer between metal atoms and BN substrate. We found that 0.73
e
charge is transferred from Ca atom to BN nanosheet through Hirshfeld charge analysis. Similar results were obtained when we introduced one Sc atom onto the boron vacancy of BN nanosheet [Fig. 3(d)]. In the optimized structure, the distance between the Sc atom and the three nitrogen atoms around the boron vacancy is 1.98 Å, and the binding energy between Ca and vacancy is 11.29 eV, which is far larger than the cohesive energy of Sc bulk metal (4.13 eV). The Hirshfeld charge analysis shows that Sc carries a 0.71
e
positive charge. As summarized in Table 1 and Fig. 4, lager binding energies are obtained in the cases of Ca and Sc adsorption in BN nanosheets with boron vacancies, accompanied by greater charge transfer. The larger binding energies are favorable to exclude the clustering problem.
The boron vacancies efficiently enhance the metal binding and dispersion. Recently, Iijima
et al. have found abundant vacancy defects in the single layer of hexagonal BN prepared by means of controlled energetic electron irradiation, and they successfully validated that the predominant atomic defects are boron vacancies [
47]. Since Ca and Sc atoms can bind tightly on BN nanosheets with
and avoid the clustering problems, metal-decorated BN nanosheets may be available experimentally and are promising hydrogen adsorption media.
The enhanced binding energies of transition metals at boron vacancy sites can be simply attributed to the more pronounced charge transfer; however, to get more insight, we have computed the density of states (DOS) for Sc-decorated pristine and defective BN nanosheets and compared with pristine and defective BN nanosheets. As shown in Fig. 5, pristine BN nanosheet is a typical semiconductor with a wide band gap. When an Sc atom is adsorbed on the BN nanosheet, the Fermi level is shifted into the conduction band, and new states appear around Fermi level. The partial DOS analysis demonstrates that these states are mainly from the Sc atom. The hybridization between Sc and BN nanosheet is found, but the 4s orbital of Sc atom is still occupied, which indicates that the hybridization between Sc and BN nanosheet is not strong. When a boron vacancy is introduced into BN nanosheet, an acceptor-like state appears near the Fermi level, which is mainly contributed by the N atoms around the vacancy site [ in Fig. 5(c)]. Once Sc atom is bound, the above acceptor-like state disappears, and both 4s and 3d orbitals of Sc are unoccupied, indicating a strong hybridization between Sc atom and BN nanosheet, which can further explain the large binding energy of Sc atom in this case.
3.3 Hydrogen adsorption in metal-decorated BN nanosheets with vacancies
First, we checked the hydrogen adsorption on BN nanosheets with vacancies. We placed one hydrogen molecule onto the boron vacancy. After full relaxation, the dissociation of hydrogen molecule into hydrogen atoms is detected in the simulation. As shown in Fig. 6(a), molecular hydrogen is chemisorbed to boron vacancy; strong chemical bonds form between hydrogen and nitrogen atoms around the B vacancy. The bond length of N–H is 1.02 Å, and the binding energy is 5.32 eV/H
2. The whole dissociation process of hydrogen molecule is exothermic when introducing defects to BN nanosheets, in agreement with a previous report on BN nanotubes [
48]. Similar results are obtained when we add hydrogen molecule onto the nitrogen vacancy, as depicted in Fig. 6(b). The hydrogen molecule is dissociated to form B–H bonds with the bond length of 1.20 Å, which is a little longer than that of N–H, and the binding energy is 0.88 eV/H
2. Chemisorption of hydrogen with large binding energies is unfavorable for the release of hydrogen at ambient conditions. In order to find the intermediate state between physisorption and chemisorption with a binding energy of 0.2–0.4 eV and achieve the reversible hydrogen uptake and release at ambient conditions, we then used Ca or Sc-decorated defective BN nanosheets to attach hydrogen molecules.
We investigated the interaction between Ca-decorated defective BN nanosheets and hydrogen molecules. As illustrated in Fig. 7, the equilibrium Ca–H bond length is about 2.55 Å. The H–H bond length is elongated to 0.76 Å (0.75 Å for isolated H
due to the interaction between Ca and H
2, and the binding energy is –0.1 eV/H
2 at the GGA level, which is smaller than those of many Ca decorated complexes [
49-
54]. The larger difference in electronegativity between nitrogen and boron leads to partially ionic bonding; therefore, the electronic structures of BN nanosheets and graphene differ apparently in spite of their isoelectronic nature. The polarized B–N bonds may account for the smaller binding energies. As more H
2 molecules approach Ca-decorated BN nanosheets with V
B, the average hydrogen adsorption energies, the distances between H
2 and Ca, and the H–H bond lengths change accordingly. The binding energy is slightly reduced from -0.10 eV to -0.07 eV at the GGA level, which may be due to the steric repulsion when the number of H
2 molecules increase. The H–H bond length is about 0.76 Å, and the average bond length between H
2 and Ca is about 2.66 Å. The adsorption energy is so weak that the hydrogen adsorption is in an unstable state at ambient conditions.
Next, we investigated the interaction between Sc-decorated defective BN nanosheets and hydrogen molecules. Hydrogen molecule is adsorbed with an elongated bond of 0.76 Å, and the bond length of Sc–H is 2.47 Å. The hydrogen adsorption energy is –0.19 eV/H
2 for GGA and –0.35 eV/H
2 for LDA. Such optimal molecular hydrogen binding energies make hydrogen adsorption and desorption feasible at ambient conditions, which is critical for practical applications. H
2 nondissociative binding to Sc-decorated BN nanosheets with V
B presents the adsorption energies significantly larger than those of typical vdW force but much smaller than those of metal hydrides. The larger adsorption energies result from the Kubas interaction [
55]. The
σ electrons of H
2 are transferred to the empty d orbitals of Sc atoms. At the same time, some filled d orbitals donate its electrons back to
antibonding orbital of H
2 to form back–donation bonds. The back–donation bond is crucial for H
2 adsorption to Sc; the bond length of H
2 is enlongated, and the degree of enlongation is related to the intensity of the back–donation bonds. If transition metals donate enough electrons to the
antibonding orbital, the H–H bond will be even dissociated.
As four H2 molecules are adsorbed onto Sc-decorated BN nanosheets, the optimized H–H bond length is 0.76 Å, and the average bond length between H2 and Sc is about 2.76 Å (Fig. 8). The binding energy is slightly reduced to –0.1 eV/H2 at GGA and –0.26 eV/H2 at LDA. When more hydrogen molecules are added onto Sc atoms, there are not enough electrons that can back donate to the antibonding orbital of H2 to destabilize H2, and the Kubas interaction is weakened correspondingly. The gravimetric density of hydrogen storage is estimated as 6 wt% if Sc atoms are placed on both sides of a 2×2 supercell of BN nanosheet with a B vacancy. The principles disclosed in this work may guide the future material design for hydrogen storage.
4 Conclusions
In summary, we have investigated the hydrogen adsorption in metal-decorated defective BN nanosheets by means of first-principles computations. The effects of vacancies are investigated on the binding of metal atoms to BN nanosheets. The binding energies of Ca and Sc on pristine BN nanosheets are far lower than the cohesive energies of the corresponding bulk metals; therefore, Ca and Sc atoms tend to form clusters on the surface of pristine BN nanosheets. Although the binding of metal atoms on nitrogen vacancies of BN nanosheets is enhanced, the binding energies are still smaller than the cohesive energies of the bulk metals, and then, Ca and Sc atoms persist in clustering on BN nanosheets with N vacancies. The situation changes in the case of BN nanosheets with B vacancies; the binding energies of Ca and Sc atoms on boron vacancies are much larger than the cohesive energies of the bulk Ca and Sc metals, and Ca and Sc atoms do not cluster on the BN nanosheets with B vacancies. Boron vacancies efficiently avoid the clustering problems due to charge transfer between metal atoms and BN substrates.
When hydrogen molecules are placed onto the vacancies of BN nanosheets, molecular hydrogen dissociates into hydrogen atoms to form chemical bonds with the atoms around the vacancy. However, as we place hydrogen molecules onto Ca or Sc decorated BN nanosheets with vacancies, the nondissociative hydrogen binding is achieved, and the intermediate hydrogen adsorption energy is obtained. Sc decorated BN nanosheets with boron vacancies demonstrate the optimal hydrogen adsorption characteristics in terms of the hydrogen adsorption energy up to –0.19–0.35 eV/H2. The Kubas interaction is responsible for the enhanced hydrogen adsorption.
Metal atoms can bind tightly on BN nanosheets with experimentally available boron vacancies and avoid clustering problems. Metal-decorated defective BN nanosheets are promising hydrogen storage media.
©Higher Education Press and Springer-Verlag Berlin Heidelberg 2011