Interstitial hydrogen in Laves phases – local electronic structure modifications from first-principles

Jana Radaković*, Katarina Batalović, Ivan Mađarević and Jelena Belošević-Čavor
VINČA Institute of Nuclear Sciences, University of Belgrade, PO Box 522, 11001 Belgrade, Serbia. E-mail: janar@vinca.rs; Tel: +381 113408 601

Received 22nd August 2014 , Accepted 17th October 2014

First published on 17th October 2014


Abstract

Understanding the microscopic aspect of the hydride formation process provides an insight into the experimentally observed properties of prospective hydrogen storage materials. In this paper, we have studied the local structural and electronic modifications induced by hydrogen absorption in cubic C15 Laves phases AB2 (A = Zr; B = Cr, Mn, Ni), as well as the stability of the formed hydrides, by means of density functional theory (DFT). To address the effect of hydrogen absorbed in one of three tetrahedral sites (96g, 32e, and 8b) on the electronic structure of its surrounding atoms, we have calculated the electric field gradient (EFG) on the position of Cr, Mn, and Ni in pure and hydrogenated compounds. EFG is associated with the hydrogen site-preference, and formation enthalpies of ZrB2H hydrides are used to examine their formation feasibility. Obtained enthalpies reveal that ZrMn2H and ZrNi2H are both unstable regardless of the occupied site, and the only attainable hydride is ZrCr2H, with comparable occupational probability of sites 96g and 32e. EFG results indicate that a hydrogen distribution within the crystal depends on the level of induced electronic structure modifications; i.e., the hydrogen site-preference is governed by the condition of minimal divergence of the electronic charge from its initial distribution.


1. Introduction

Laves phases are a large class of intermetallic compounds that exhibit reasonably high hydrogen storage capacity and favourable absorption/desorption kinetics around room temperature and ambient pressures. Consequently, they are being extensively investigated as potential hydrogen storage materials, or electrodes for Ni/MH batteries;1–8 however, optimization of their performance and practical properties is still necessary for such applications. Formation enthalpies and binding energies of hydrides, and changes in the ground state structure that occur after hydrogenation are constantly being reported for various types of compounds, but linking the practical features of the material to its microscopic structure still remains a challenging issue. For example, subtle modifications in the electronic structure of hydrogen storage materials, induced by the absorbed hydrogen, can stabilize or destabilize initial compounds. Evaluating the cause of such behaviour contributes to better understanding of the hydride formation process and experimentally observed properties of hydrides. Consequently, it is of paramount importance to thoroughly investigate the effect that hydrogen has on the electronic structure of prospective hydrogen storage materials, and this study represents a contribution in that regard. It deals with the local electronic structure modifications that occur in cubic Laves phase compounds ZrCr2, ZrMn2, and ZrNi2, due to the formation of corresponding monohydrides ZrCr2H, ZrMn2H, and ZrNi2H.

Due to the experimentally observed storage capability of ZrCr2 and ZrMn2, as well as of the alloys from Zr–Ni phase system; all selected compounds have been previously studied. So far, the research focus has been based on their electrochemical performance and tuning of thermodynamic and kinetic properties by partial substitution of host atoms with various transition metals (e.g., V, Cr, Mn, Fe, Ni) and rare earths (La, Ce). In that regard, their hydrogen absorption capacity and absorption/desorption kinetics have been studied in the case of both pure and multicomponent systems.1,2,9,10 Pure and hydrogenated ZrCr2 is easily synthesized in cubic C15 structure, with the maximum storage capacity of 3 hydrogen atoms per formula unit. Additionally, pressure–composition–temperature (PCT) measurements revealed its satisfactory thermodynamic properties regarding the needed electrochemical stability.2,9 ZrMn2 is found to readily crystallize in the hexagonal C14 structure, and is able to store up to 3.6 hydrogen atoms per formula unit;10,11 however, several theoretical studies indicate the possibility of existing structural degeneracy at low temperatures, and forming of cubic C15 ZrMn2 compound and ZrMn2H0.5 hydride.1,7 Furthermore, optimization of the synthesis process of compounds in the Zr–Ni phase system led to the potential discovery of C15 Laves phase ZrNi2. Recently, Kumar et al.12 studied the effect of various concentrations of Ni on the ZrMn2−xNix phase transformation from hexagonal C14 ZrMn2 to cubic C15 ZrNi2 compound, and revealed the existence of cubic Laves phase ZrNi2. Acknowledging the previous studies, we have decided to examine the possibility of ZrMn2H and ZrNi2H formation in cubic configuration by calculating their heats of formation.

2. Details of calculations

2.1. Crystal structure of studied compounds and calculation methodology

The cubic AB2 Laves phase crystallizes in MgCu2 structure with Fd[3 with combining macron]m (O7h) space group, with atom A located on 8a cubic site and atom B on 16d noncubic Wyckoff site. During the absorption process, hydrogen can occupy three interstitial positions – 96g, 32e, or 8b. Unit cell of this structure with emphasized interstitial tetrahedra is presented in Fig. 1; for the sake of clarity some atoms are omitted from this figure.
image file: c4ra09082a-f1.tif
Fig. 1 Cubic C15 Laves phase unit cell, with indicated 96g (A2B2), 32e (AB3), and 8b (B4) tetrahedral interstitials. Some atoms are omitted for clarity.

Interstitial position 96g is a tetrahedral site formed by two A and two B atoms, A and B being the host atoms of AB2 Laves phase. For simplicity, this site and respective hydride will be addressed as A2B2. In analogy to previous position, 32e site, which is constructed by one A and three B atoms, will be referred to as AB3, and finally, 8b as B4.

Depending on the size of the interstitial site, and type of A and B atoms, hydrogen preferably occupies A2B2 and/or AB3 position, although, in some compounds B4 site can also be occupied, but with considerably lower probability.1,3,13 Considering that the foremost interest of this study is to observe the changes in electronic structure of atoms, in various tetrahedral arrangements, in the presence of hydrogen, and consequently provide an explanation for the hydrogen occupational site-preference, we have neglected the occupation probability and equally studied all three interstitial positions.

The used methodology assumed subsequent introduction of hydrogen into the A2B2, AB3, or B4 tetrahedral site, in order to construct AB2H hydride. It was expected that due to different hydrogen's first coordination spheres (see Fig. 1.) the local crystal structure and electronic properties of host metals would be differently modified. Because the electric field gradient (EFG) directly reflects both geometrical and electronic aspects, induced modifications were studied by calculating it on the atomic site of Cr, Mn, and Ni, before and after hydrogenation.

EFG originates from the charge density deviation from spherical symmetry in the vicinity of observed nucleus, and it is extremely sensitive even to trivial changes in the electronic structure. Mathematically, it is defined as the second partial derivative of the electrostatic potential near nucleus (V([r with combining right harpoon above (vector)]), [r with combining right harpoon above (vector)] → 0) with respect to spatial coordinates Vzz = ∂2V(0)/∂xixj.14,15 It is designated as a traceless, symmetric second rank tensor, diagonal in the principal-axis coordinate system, and it is always present at the atomic site with the noncubic point group symmetry. According to a convention, its three nonzero components are ordered as |Vzz| ≥ |Vyy| ≥ |Vxx|; however, it is completely described by its largest component, comparable with the experimental values. Furthermore, in full-potential formalism, used in present calculations, total EFG (Vtotii) equals a sum of the valence EFG (Vvalii), which is a linear combination of partial valence electronic charge contributions image file: c4ra09082a-t1.tif, and the lattice EFG (Vlatii).

The heats of formation were calculated by defining the energy difference between the total energy of studied hydrides and the total energy of the reactants; details can be found in our previous work.16 Our calculations reflect the situation at 0 K, and do not incorporate the vibrational frequencies associated with zero point energy, since in various types of Laves phases they are significantly smaller then ΔH(AB2H). In addition, their values in different interstitial sites are mutually comparable; therefore, they do not disrupt the calculated trend.17

2.2. Computational details

The calculations were performed using a full-potential (linearized) augmented plane-wave ((L)APW) + local orbitals (lo) method, as implemented in WIEN2k,14,18 a program package based on density functional theory (DFT).19 WIEN2k describes the charge density without any shape restriction in the entire region of the crystal. Therefore, the charge density and the Coulomb potential are calculated very accurately in close proximity of the nucleus, as well as in the interstitial region between the atoms. Consequently, quantities like EFG are calculated very precisely.14 The radii of the muffin-tin spheres, centred on the atomic nucleus, (Rmt) were set to 2.2 Bohr for Zr, 1.95 Bohr for Ni, Mn, and Cr, and 0.95 Bohr for H. The number of basis functions used in the interstitial region was determined by RmtKmax parameter, which was set to 4.5; Kmax signified the magnitude of the largest plane wave vector, while Rmt denoted the smallest MT radius in the unit cell. Inside the MT spheres wave functions were expanded in the spherical harmonics up to lmax = 10. The cutoff energy that separated the core and valence states was −6 Ry for ZrCr2 and ZrNi2, and −7 Ry for ZrMn2. The exchange and correlation effects were included within the generalized gradient approximation (GGA), using the Perdew–Burke–Ernzerhof (PBE) scheme.20 Brillouin zone integration was performed via the tetrahedron method, and depending on the occupied interstitial position, a mesh of 275 to 475 k points was used in the irreducible Brillouin zone. After testing the effect of spin polarization on studied values it was established that it has minor significance, as a result we have decided to neglect it. Self-consistency was achieved by demanding that the difference of integrated charge in successive iterations was less than 10−5e. Relaxation of all atomic positions was performed until forces acting on all atoms were less than 1 mRy per Bohr.

3. Results and discussion

Lattice expansion, due to AB2H monohydride formation, is calculated with respect to pure compounds. The optimized lattice constants are presented in Table 1. According to obtained results lattice expansion depends on the type of occupied interstitial site and the chemical nature of host metals (Cr, Mn, or Ni). Relative to the initial volumes, the expansion ranges from 5% to 15% per hydrogen atom, which is in accordance with previous studies.1,3,13 Comparison of volume changes in series of monohydrides reveals nearly identical structural changes in all hydrogenated compounds. On the basis of this, it was difficult to evaluate to what extent hydrogen affects the studied structures. Therefore, we focused our research on observing the local crystal structure and electronic properties of atoms that surrounds the absorbed hydrogen.
Table 1 Calculated lattice parameters of AB2 and AB2H compounds
  ZrNi2 ZrMn2 ZrCr2
C15 6.991 6.995 7.134
B4 7.243 7.238 7.384
AB3 7.155 7.147 7.318
A2B2 7.136 6.958 7.302


3.1. Electric field gradients on the atomic site of Cr, Mn, and Ni in pure and hydrogenated compounds

EFG is considered to be a local property because it is determined by the electronic charge distribution of the selected atom, and by the spatial distribution of its neighbouring atoms. Subtle modifications of the local electronic structure induced by the absorbed hydrogen are sufficient to cause significant variations in the EFG value, and furthermore, to stabilize or destabilize a starting compound. Establishing the physical foundation of such correlation would be beneficial in providing the insight into the process of hydride formation, and experimentally observed properties of hydrogen storage materials.

Total EFGs on the atomic site of Cr, Mn, or Ni, the corresponding valence component, and partial orbital p (Vppzz) and d (Vddzz) contributions, have been calculated in pure and hydrogenated compounds and results are presented in Fig. 2.


image file: c4ra09082a-f2.tif
Fig. 2 Total EFG, valence EFG, and p and d orbital contributions calculated on the atomic site of Ni, Mn, or Cr, in pure and hydrogenated C15 Laves phases.

In all compounds total EFG is entirely determined by its valence component VtotzzVvalzz, that is, by the charge density within the MT sphere. In pure and hydrogenated compounds value and sign of the EFG depend on the type of B atom, while in hydrides it also depends on the atom's spatial arrangement in the tetrahedral site itself.

In series of pure compounds, stable ZrCr2 exhibits the lowest absolute EFG value, which indicates that compared to Mn or Ni a lesser charge deviation from the spherical symmetry occurs in the vicinity of Cr. In ZrCr2Hx and ZrMn2Hx (x = 0, 1) EFG exhibits an analogous increase of its absolute value throughout the series of tetrahedral sites. Established trend and its evident discrepancy in case of ZrNi2Hx can be understood in view of Cr, Mn, and Ni electronic configuration. However, to comprehend the importance of this, one should initially perceive that for all studied hydrides partial p and d terms predominantly contribute to the valence EFG in a form of linear combination VvalzzVppzz + Vddzz (see Fig. 2.). Atomic Cr and Mn have comparable electronic configuration (Cr: 4s1 3d5 4p0; Mn: 4s2 3d5 4p0), which differs from that of Ni (4s2 3d8). As a result of that, and their identical site symmetry, it can be expected that d and p charge in Cr and Mn are equivalently spatially distributed. Considering the spherical symmetry of s orbital, contribution from s states is negligible compared to Vppzz and Vddzz. Providing this explanation, it is relatively simple to understand the established trend and the lack of it in Ni-hydrides.

It was briefly mentioned that the governing factor in the occupation process is the volume of the interstitial position. Accordingly, because it possesses the largest volume, A2B2 site will be preferably occupied and followed by the occupation of somewhat smaller AB3 tetrahedra.21 The smallest available B4 site usually remains unoccupied; however, negative formation enthalpies for several B4 hydrides have been reported.13 At this point, we will correlate the obtained EFG values with the hydrogen occupational site-preference in studied compounds, and add a presumption to the “site-volume–site-preference” rule. According to our results, absorption of hydrogen in A2B2 site only slightly modifies the EFG values from the ones calculated in corresponding pure compounds. On the other hand, considerably pronounced divergence of the electronic charge from the original distribution is observed when hydrogen occupies B4 tetrahedral site. On the basis of this, it can be argued that the interstitial positions of studied Laves phases are being occupied under the condition that the electronic charge of its surrounding atoms minimally diverges from the initial distribution.

3.2. Density of states (DOS) of studied compounds

Previous presumption can be supported by analysing the density of states (DOS) of hydrides and corresponding pure compounds, which is presented in Fig. 3. DOS is decomposed according to atoms of all studied structures, and given with respect to the Fermi level (EF). In order to effectively signify hydrogen electron states, on the graphs of all hydrides are provided insets that represent the rescaled DOS in the equivalent energy range. The lowest conduction band in DOS graphs, which consists of hydrogen 1 s states and overlapping host metal bands, is located in the energy range approximately from −5 eV to −11 eV (see insets in Fig. 3.). However, its relative position with respect to EF depends on the occupied interstitial site. In comparison to A2B2 and AB3, in B4 hydride lowest conduction band is shifted further away from the EF. Also, in all studied structures it is closest to the EF in A2B2 hydride. This issue is addressed considering that several studies indicate that hydrides are more stable if hydrogen electrons states shift higher in energy; i.e., closer to EF.22,23 Furthermore, the shape of DOS, which corresponds to the distribution of occupied m-decomposed states, also depends on the type of the occupied interstitial site. DOS shape of A2B2 hydride, and to some extent of AB3, is comparable with the DOS outline of pure compounds, while in B4 hydrides it is significantly altered. This infers that compared to the hydrogen absorption in A2B2 or AB3 site, energetically and spatially less favoured modification of the initial electronic structure occurs after its absorption in B4 tetrahedral site. This confirms the aforementioned presumption, according to which modifications of the original electronic charge distribution of atoms surrounding the absorbed hydrogen should be minimal.
image file: c4ra09082a-f3.tif
Fig. 3 Density of states of ZrCr2, ZrMn2, and ZrNi2, as well as of their monohydrides, decomposed according to atoms. Insets: DOS rescaled to 1 per eV per atom to represent corresponding hydrogen states.

Further observation can be made regarding the correlation between the Fermi level position and the stability of studied structures. In pure ZrCr2, Fermi level clearly separates the bonding from the antibonding states, while in ZrMn2 it is shifted towards the antibonding area. This agrees with the pronounced stability of ZrCr2, and infers that ZrMn2 compound does not crystallize in cubic C15 structure.

In addition, according to the position of EF, all ZrMn2H interstitial hydrides are likely to be unstable. Position of EF in DOS of ZrCr2H reveals the stability of A2B2 and AB3 hydrides, and suggests the instability of B4 hydride. DOS of ZrNi2 and ZrNi2H does not afford such distinct assumptions, which is why in order to gain insight into its practical features formation enthalpies of this hydride have to be addressed.

3.3. Heats of formation of studied hydrides

In Fig. 4 are presented calculated formation enthalpies of all studied hydrides. It should be stressed out that the conditions under which the studied compounds are used are estimated according to a change in Gibbs free energy, not merely a change in enthalpy. However, several studies determined similar entropy changes in the solid hydride formation reactions, justifying the use of enthalpy change for the estimation of hydrides stability.21
image file: c4ra09082a-f4.tif
Fig. 4 Calculated formation enthalpies of studied hydrides.

According to the formation enthalpies of all three studied cubic Laves phase compounds the only stable monohydride is ZrCr2H. ZrNi2H and ZrMn2H have positive formation enthalpies, which indicates that, regardless of the occupied interstitial site, all suggested hydrides are physically unattainable. When hydrogen is accommodated in A2B2 or AB3 interstitial site of ZrCr2, negative formation enthalpies are obtained, according to which both hydrides are stable. On the other hand, positive enthalpy of B4 hydride suggests that occupation of this site destabilizes the initial compound. Calculated enthalpies of A2B2 or AB3 hydrides are mutually comparable (−16.98 kJ mol−1 H−1 for AB3 and −18.17 kJ mol−1 H−1 for A2B2) and lead to an assumption that in ZrCr2 their occupation is competitive. To be able to directly compare our values with those from previous studies, experimental or computational, it would be beneficial to have the formation enthalpies for ZrCr2Hx, with x = 1. However, in lack of these results we can relate our values to the existing ones in respect to the range of available hydrogen concentrations. Previously studied concentrations ranged from 0.5 H per f.u. to ∼3 H per f.u., and published formation enthalpies listed from the experimental value of −18 kJ mol−1 H−1 for ZrCr2H3 (ref. 10 and 24) to theoretical of −33.5 kJ mol−1 H−1 for ZrCr2H0.5.1–3 Our value for A2B2 hydride agrees excellently with the experimental one for 3 hydrogen atoms per formula unit.

Positive formation enthalpies of all ZrNi2H and ZrMn2H indicate that regardless of the occupied interstitial site these hydrides are always unstable. The obtained values suggest difficulties in any prospective synthesis attempt; nevertheless, hydrogenation of potentially synthesized cubic C15 ZrMn2 or ZrNi2 may also result in the phase transformation into hexagonal C14 Laves phase, or even some other stable structure. Our results are in agreement with the experimental findings of Kumar and coworkers.12 On the other hand, present computational and previous experimental results contradict the first-principles calculations conducted on hypothetical ZrMn2H0.5,1 for which the authors obtained negative formation enthalpies. In addition, it is expected that the formation of hydrides in Zr–Ni system would be less probable with increasing Ni content, given that Ni, being the late transition metal, does not form stable binary hydride, while Zr does (ZrH2). Baring in mind that ZrNi2 forms a hydride of significantly decreased stability, with respect to the pronounced stability of ZrNi hydrides,25,26 further investigation of ZrNixHy alloy system (1 < x < 2) could present a prospective way of tailoring the desired stability of the Zr–Ni hydrides.

4. Conclusions

First-principles calculations of the electronic structure of AB2 cubic C15 Laves phases (A = Zr; B = Cr, Mn, Ni), and of their monohydrides, showed how the absorption of only one hydrogen atom substantially modifies the local electronic structure. Electronic charge distribution of hydrogen surrounding atoms was associated with the hydrogen formation process by calculating the electric field gradient (EFG) on Cr, Mn, or Ni, when hydrogen occupied one of three interstitial positions (A2B2, AB3, B4). Comparison of EFG values indicated that the hydrogen distribution within the crystal lattice is governed by the condition of minimal modification of the initial electronic charge distribution. Consequently, compared to the A2B2, or even AB3 site, absorption of hydrogen in the B4 interstitial site is not likely to occur in any of the studied compounds, as it considerably modifies the starting electronic structure of its surrounding atoms. This observation is in accordance with the previous results; in addition, it complements the former conclusion, according to which occupational site-preference of atomic hydrogen in cubic Laves phases is governed by the volume of interstitial tetrahedra (A2B2 → AB3 → B4). Furthermore, provided presumption was supported by the analysis of densities of states obtained before and after hydrogenation, which implied energetically less favoured modification of the electronic structure with respect to pure compounds in case of B4 or AB3 hydrides.

Calculated formation enthalpies indicated that of all studied hydrides only ZrCr2H is the stable one; with hydrogen occupying A2B2 or AB3 interstitial site. Both ZrMn2H and ZrNi2H were found to be unstable, regardless of the occupied interstitial site. Our values agree with experimental results for all hydrides; however, they contradict the previously calculated formation enthalpies of Mn-hydride, according to which this hydride in cubic C15 structure can be stable.

Presented results suggest the existence of correlation between the electronic structure of atoms surrounding the absorbed hydrogen and hydride formation process in the cubic Laves phases. However, to provide a deeper insight into such correlation, and additionally study provided conclusions, it would be beneficial if an extensive investigation is conducted on a broader class of prospective hydrogen storage materials.

Acknowledgements

The Ministry of Science and Technological Development of the Republic of Serbia provided the financial support for this study through the Project no. 171001. Authors would like to thank Dr Carmen Rangel from the LNEG institute in Lisbon, Portugal, for useful discussions during the meetings organized through Portugal/Serbia bilateral project (no.: 451-03-02328/2012-14/04).

Notes and references

  1. S. Hong and C. L. Fu, Phys. Rev. B: Condens. Matter Mater. Phys., 2002, 66, 094109 CrossRef.
  2. J. Bodega, J. F. Fernandez, F. Leardini, J. R. Ares and C. Sanchez, J. Phys. Chem. Solids, 2011, 72, 1334–1342 CrossRef CAS.
  3. S. B. Gesari, M. E. Pronsato, A. Visintin and A. Juan, J. Phys. Chem. C, 2010, 114, 16832–16836 CAS.
  4. A. Ormeci, F. Chu, J. M. Wills, T. E. Mitchell, R. C. Albers, D. J. Thoma and S. P. Chen, Phys. Rev. B: Condens. Matter Mater. Phys., 1996, 54, 12753 CrossRef CAS.
  5. A. Skripov, J. Combet, H. Grimm, R. Hempelmann and V. Kozhanov, J. Phys.: Condens. Matter, 2000, 12, 3313–3324 CrossRef CAS.
  6. B. A. Talagañis, M. R. Esquivel and G. A. Meyer, Int. J. Hydrogen Energy, 2009, 34, 2062–2068 CrossRef.
  7. X. Chen, W. Wolf, R. Podloucky, P. Rogl and M. Marsman, Phys. Rev. B: Condens. Matter Mater. Phys., 2005, 72, 054440 CrossRef.
  8. X. Chen, W. Wolf, R. Podloucky and P. Rogl, Phys. Rev. B: Condens. Matter Mater. Phys., 2005, 71, 174101 CrossRef.
  9. J.-M. Joubert, D. Sun, M. Latroche and A. Percheron-Guegan, J. Alloys Compd., 1997, 253–254, 564–569 CrossRef CAS.
  10. A. T. Pedziwiatr, R. S. Craig, W. E. Wallace and F. Pourarian, J. Solid State Chem., 1983, 46, 336–341 CrossRef CAS.
  11. F. Pourarian, V. K. Sinha, W. E. Wallace and H. Kevin Smith, J. Less-Common Met., 1982, 88, 451–458 CrossRef CAS.
  12. V. Kumar, D. Pukazhselvan and A. K. Tyagi, Mater. Renew. Sustain. Energy, 2013, 8, 12–20 CrossRef.
  13. J. Radaković, J. Belošević-Čavor and V. Koteski, Int. J. Hydrogen Energy, 2013, 38, 9229–9235 CrossRef.
  14. K. Schwarz, C. Ambrosch-Draxl and P. Blaha, Phys. Rev. B: Condens. Matter Mater. Phys., 1990, 42, 422051 CrossRef.
  15. K. Schwarz, J. Solid State Chem., 2003, 176, 319–328 CrossRef CAS.
  16. K. Ćirić, V. Koteski, D. Stojić, J. Radaković and V. Ivanovski, Int. J. Hydrogen Energy, 2010, 35, 3572–3577 CrossRef.
  17. F. Li, J. Zhao, D. Tian, H. Zhang, X. Ke and B. Johansson, J. Appl. Phys., 2009, 105, 043707 CrossRef.
  18. P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka and J. Luitz, WIEN2k An Augmented Wave Plus Local Orbitals Program for Calculating Crystal Properties. User's Guide, Vienna Universtiy of Technology, Austria, 2012 Search PubMed.
  19. W. Kohn, Rev. Mod. Phys., 1999, 71, 1253–1266 CrossRef CAS.
  20. J. P. Perdew, S. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865–3868 CrossRef CAS PubMed.
  21. A. Züttel, A. Borgschulte and L. Schlapbach, in Hydrogen as a Future Energy Carrier, ed. A. Züttel, Wiley-VCH Verlag and Co, Weinheim, 2007, ch. 6, p. 192 Search PubMed.
  22. R. Griessen and A. Driessen, Phys. Rev. B: Condens. Matter Mater. Phys., 1984, 30, 4372–4381 CrossRef CAS.
  23. X. Ke, G. J. Kramer and O. M. Løvvik, J. Phys.: Condens. Matter, 2004, 16, 6267–6277 CrossRef CAS.
  24. A. V. Krupenchenko, E. P. Magomedbekov and I. I. Vedernikova, J. Phys. Chem., 1990, 64, 2897–2908 CAS.
  25. S. F. Matar, Chem. Phys. Lett., 2009, 473, 61–65 CrossRef CAS.
  26. K. Batalović, V. Koteski and D. Stojić, J. Phys. Chem. C, 2013, 117, 26914–26920 Search PubMed.

This journal is © The Royal Society of Chemistry 2014
Click here to see how this site uses Cookies. View our privacy policy here.