Abstract

The relationship between microstructure and hydrogenation properties of the mixed metals has been investigated via different spectroscopic techniques and the density functional theory (DFT). FESEM and TEM analyses demonstrated the nano-grains of Mg2NiH4 and MgH2 on the hydrogenated microstructure of the adsorbents that were confirmed by using XPS analysis technique. SAED pattern of hydrogenated metals attributed the polycrystalline nature of mixed metals and ensured the hydrogenation to Mg2NiH4 and MgH2 compounds. Flower-like rough surface of mixed metals showed high hydrogenation capacity. The density functional theory (DFT) predicted hydrogenation properties; enthalpy and entropy changes of hydrogenated microstructure of MgH2 and Mg2NiH4 are −62.90 kJ/mol, −158 J/mol·K and −52.78 kJ/mol, −166 J/mol·K, respectively. The investigation corresponds to the hydrogen adsorption feasibility, reversible range hydrogenation thermodynamics, and hydrogen desorption energy of 54.72 kJ/mol. DFT predicted IR band for MgH2 and Mg2NiH4 attributed hydrogen saturation on metal surfaces.

1. Introduction

The global energy demand is increasing in parallel with the population growth, economic expansion, and increasing demand of mobility. Hydrogen is known as an abundant source of clean energy carrier. There are few ways to utilize the hydrogen as source of energy. Metal-hydrogen system demonstrated the best option for hydrogenated system due to the higher hydrogen density in the system than liquid or gas phases [1]. Higher hydrogen density, reversible thermodynamics, low desorption temperature, and low production cost are vital requirements for metal-hydrogen system. Magnesium based mixed oxides are potential candidate to meet the above requirements [2]. The drawbacks of hydrogenation in magnesium based alloys are [3] stability of the MgH2 phase, formation of a surface oxide(s), slow dissociation of hydrogen at the metal surface, and slow diffusion of hydrogen through MgH2 grain. The problems can be overcome by developing microstructure that reduces the diffusion path and creates enough active site for dissociation [2, 3]. Grain boundaries provide active nucleation sites in hydride formation and decomposition of the hydride phases. Reduction of grain size or crystalline size of the materials enhances the hydrogen diffusion that moves faster along the grain boundaries.

Different strategies have been investigated to achieve reversible range hydrogen adsorption thermodynamics and kinetics. Homogenous and ultrafine microstructure of adsorbent is one of the important techniques to enhance the hydrogenation capacity and to achieve favorable thermodynamics [4]. Second technique is considered as addition of catalytic elements or rear earth metals or transitional metals to adsorbent. Transition metal oxides are an important catalyst in Mg-based adsorbent to get high hydrogenation capacity [5]. Mg2Ni alloy demonstrates higher hydrogenation kinetics and desorption at lower temperature because of the effect of transitional metal such as nickel [6]. Coprecipitation under low supersaturation synthesis method produces homogenously metals dispersed hydrotalcite [7]. The materials show high catalytic activity for hydrogenation and flexible to form different microstructure by optimizing the preparation conditions and metals.

Many DFT based studies have been reported for magnesium based hydrogen storage adsorbents [8, 9]. The density functional theory (DFT) based QM (quantum mechanical) descriptors describe adsorption properties such as adsorption thermodynamics and kinetics, interaction energy, and adsorption and desorption energy that support to observe the hydrogenation mechanism and properties [10] of the materials. Noncovalent week interaction energies can be determined precisely using M05-2X functional [11].

Microstructure influences hydrogenation or dehydrogenation properties strongly [12]. It is a novel initiative to make a clarification of hydrogenation on reduced mixed metal oxides microstructure via experimental and theoretical way. In this study, hydrogenated microstructure is investigated through electron microscopic techniques and studied hydrogenation/dehydrogenation properties using the density functional theory (DFT) to establish a relation between hydrogenation and microstructure of the material.

2. Experimental Section

2.1. Experimentation

Hydrotalcite based mixed oxides containing Mg, Ni, and Al were prepared by using the conventional coprecipitation method under low supersaturation method as described elsewhere [13]. The synthesis pH of 10 was optimized to get the flower-like rough surface. Optimized synthesis temperature of 55°C and molar ratio of metals were used as Mg : Ni : Al = 2 : 1 : 1. The fresh developed materials were calcined at 650°C for two hours to convert the mixed metals oxides.

The hydrogen adsorption equilibrium of mixed oxide was measured at the conditions of temperature (30°C) and pressure (0–760 mmHg) using a micrometrics ASAP 2020C sorptometer. A 20 mg sample was degassed using both preparation port and sample port. It was reduced by H2 gas purging at 20 mL/min, with a temperature ramping of 10°C/min up to 850°C, which further cool down to 30°C under vacuum condition. 99.99% hydrogen gas was used to inlet the adsorption unit. Hydrogen adsorption isotherm analyses were carried out at 30°C with a pressure of 0–760 mmHg. Reversible adsorption isotherms were collected by performing a 30-minute evacuation at the experiment temperature following the total adsorption isotherm.

The microstructure of hydrogenated mixed metals was analyzed using field emission scanning electron microscopy (FESEM) technique. CARL Zeiss Supra 55VP FESEM instrument equipped with the Oxford INCA 400 EDX microanalysis system with an operating voltage in the range of 0.1–30 kV was used. The system was vacuumed before analyses and the mixed metals were introduced in the sample grid.

The microstructure of hydrogenated adsorbent was investigated using a CARL Zeiss LIBRA 200FE, TEM using ultrasonically dispersed samples in isopropanol. Selected area electron diffraction (SAED) patterns were collected for hydrogenated single sheet. The diffraction rings of SAED pattern were indexed using the diameter of the rings and calculating the ratio of those diameters () for different phases.

The surface composition of the hydrogenated adsorbent was analysed by X-ray photoelectron spectroscopy (XPS) using a Perkin Elmer 1257 model, operating at a base pressure of 3.5 × 10−6 Pa at 300 K with a nonmonochromatized Al-Kα line at 1486.6 eV and a hemispherical sector analyzer.

2.2. Computational Methods

The hydrogenated clusters or hydrides were optimized (geometrically) using M05-2X functional [11] along with the 6-311+G(d,p) basis set [14]. All computations were performed using Gaussian 09W computational program package [15]. Frequency calculations were conducted on optimized geometries of clusters to ensure the minima on potential energy surface (PES). Cluster descriptors such as chemical potential [16], strength of charge transfer [17], and electrophilicity [18] were calculated by using finite approximation. Energy of interaction (Δ), hydrogenation energy (Δ), and average desorption energy (Δ) [19] per hydrogen molecule were calculated using different equation as stated by Salam et al. [10]. Thermodynamics parameters such as heat of reaction and free energies have been calculated using computed energies of reactants and products [20]; The IR spectrum of hydride clusters was generated using computational software (Gaussian 09W) via frequency calculation.

3. Results and Discussion

3.1. Experimental Analyses
3.1.1. Hydrogen Adsorption Mechanism

Isothermal adsorption analysis predicted the hydrogen storage capacity of 3.9 wt% hydrogen of the reduced mixed metals oxides containing magnesium, nickel, and aluminium [13]. A hydrogenation mechanism and causes of hydrogenation can be explained by pointing out different factors.

The surface of s-orbital electron and d-orbital electrons of atom plays a key role in the hydrogen dissociation on the metals surface. H2 molecules come to the metals surface in the hydrogenation process. The molecular orbitals of hydrogen (1s2) start to overlap on the metals surface s-orbital electrons. According to the Pauli Exclusion Principle, hydrogen s-orbital electron repels due to the fully occupied valence electron of metals (Mg: 1s22s22p63s2), resulting in the energy barrier. The metals with d-orbital electrons contribute in the dissociation of H2. H2 comes closer to the metal surface; the charge transfers from the H2 s-orbital to the d-orbital of metal (Ni (28): 1s22s22p63s23p63d84s2) surface and simultaneously a back-donation of charge from the surface d-orbital to the H2 antibonding state. Thus, the interaction between d-orbital of nickel and s-orbital of H2 enhances the H2 dissociation on the adsorbent surface with a small energy barrier. The main group element of aluminum (Al (13): 1s22s22p63s23p1) valence electron is three (3).

The electron affinity of Al is −44 kJ/mol and shows low affinity to H. Since the p-state of Al can accommodate electron more, it competes with H (hydrogen) for the valence electrons of Mg (Mg (12): 1s22s22p63s2). The electron affinities of magnesium and hydrogen are 21 kJ/mol and −74 kJ/mol. Therefore, hydrogen binds with weak interaction (lower bonding strength) with the Mg or Mg2Ni and forms reversible hydride of Mg2NiH4 and MgH2. Aluminum also works as cofactor to heat transfer quickly for better decomposition of hydrides from microstructure.

3.1.2. Microstructure of Mixed Metals Oxides

Figure 1 shows the flower-like rough microstructure of the mixed metals oxides. At pH = 10, cations Al3+, Ni2+, and Mg2+ precipitated simultaneously as Mg(OH)2, Ni(OH)2, and Al(OH)3 that immediately converted to the Mg(Ni)Al-CO3 structure. The flower-like particle of material can be considered as composed of hexagonal nanosheets which are dependent on isoelectric point of metal hydroxides. The nanosheet or nanoflakes were observed via high resolution FESEM image (not shown here). The nanosheets seem to be flexible and can be wrinkled during precipitation at experimental conditions that help to form flower-like particle. If pH = IEP (isoelectric point), the MgAl-CO3 nanoparticle is in net zero point charge and the formation of hydrotalcite becomes slow along the 001 plane. In the initial stage, the interfacial energy between the nuclei and the support is low and metal hydroxide nuclei are grown within their 001 plane. Since the reaction is continuous (during synthesis), the subsequent nucleation is grown preferentially on the surface of the previously grown particles. Thus, the growth of the primary nanoparticles turns to flower-like morphology. After calcination, the materials retained its shape. The reason is that the trivalent cations of Al3+ substitute the divalent cations of Mg2+, Ni2+ to the brucite sheets and converted to periclase-like Mg-Al-O or Mg(Ni)Al-O solid solution. Thus, the brucite-like shape of the crystal remained in the resulting Mg/Ni-Al-O structure. Different synthesis conditions demonstrate different microstructure for this metals combination.

EDX analysis demonstrated the elemental composition of microstructure of the mixed metals oxides in Figure 2. Since the material is in mixed oxides form, it shows high amount of oxygen (30 wt%). Magnesium, nickel, and aluminium maintain the ratio (based on wt%) of 2 : 1 : 1.

3.1.3. Microstructure of Hydrogenated Mixed Metals

Microstructure of the hydrogenated mixed metals is shown in Figure 3. Hydrogen makes bond with every flake of the flower-like particle, thus forming very small grains. Therefore, well-defined grains attribute the faster hydrogenation and high hydrogen capacity of Mg2NiH4 [6]. The size range of hydrogenated grains is in 20–40 nm and synthesized adsorbent form 100–180 nm sizes of flower-like particles in optimized experimental conditions. The rough surfaces of flower-like particles offer active sites to attach hydrogen to each flake and enhance the hydrogenation capacity. The hydride nucleation takes place at the end of flakes including grain boundaries. The hydrogenation rate and hydrogen capacity of material increase rapidly with increasing active site and for the homogenous distribution of catalytic metals of nickel on microstructure. The surface with hexagonal particle demonstrates less hydrogenation in same conditions [13] which is due to the shortage of active site of materials. The β-MgH2 is formed at ambient temperature as a rutile structure as observed by Zachariasen et al., [21].

The TEM image of the adsorbent (Figure 4) shows the flower-like particle which is also observed by FESEM image in Figure 1(a). The newly formed irregular shape of grain adjacent to the flakes of the particle is due to the formation of different phases of hydride at ambient conditions. The adsorbent shows the well-shaped diffusion path and grain boundaries. SAED (selected area electron diffraction) pattern of hydrogenated metals is shown in Figure 5 which attributes the crystal planes of (110), (202), (220), and (330) and corresponding the presence of MgH2, Mg2NiH4, and MgO in crystal. Similar planes of SAED pattern were observed by Zou et al., [22]. Different metals’ combination creates a different microstructure due to the structural change and balance band of the adsorbent that influence the hydrogenation on material surface.

3.1.4. Hydrogenated Phases of Adsorbent

Hydrogenated microstructure was analyzed by XPS techniques to confirm the hydride phases and for further studies. Binding energy of hydrides is higher than elemental species, because metals loss their electron in forming hydrides [23]. X-ray photoelectron spectra of Mg (1s), Mg (2p), and Ni (2p) are shown in previous investigation [6]. Hydride phases of MgH2 and Mg2NiH4 can be observed by measuring their higher binding energies as compared to corresponding elements. Combination of magnesium and nickel modifies the electronic structure to Mg2NiH4 that enhances the hydrogenation. Mg and Ni change the structure that reduces the formation of Mg(OH)2 and increases hydrogenation [24].

3.2. Microstructure Analysis Using DFT

The study highlights the hydrogenation properties of hydrogenated microstructure or hydrogenated phases of MgH2 and Mg2NiH4 that were confirmed by experimental analyses in this study. The geometry of clusters has been optimized successfully and shown in Figure 6. To observe the variation of hydrogenation properties, Mg-(H2)2 is included along with two hydride clusters of MgH2 and Mg2NiH4. Frequency analyses predict their minima on potential energy surface (PES) and free from imaginary frequencies (Figures 7 and 8). The electronic properties of clusters and feasibility of their application are discussed next.

Total electronic energy and thermodynamical descriptors, hardness, electronegativity, and electrophilicity of hydrides are calculated and shown in Table 1. The increasing trend of electronic energy of Mg-(H2)n based on gradual number of captured hydrogen corresponds hydrogenation on metals surface. Chemical hardness determines the clusters strength to donate/accept the charges and the chemical potential represents the electron shifting trend from systems or clusters. Minimum electrophilicity principle (MEP) and maximum hardness principle (MHP) attribute the Mg2NiH4 hydride which is more stable than Mg-(H2)n [25].

The IR studies on hydrogenated clusters provide two significant properties. IR spectrum of hydrogenated clusters shows that they are free from imaginary frequency indicating that the minima on potential energy surface (PES). Secondly, it confirms the hydrogen saturation on metal surface. IR peaks in the frequency region of 400–2000 cm−1 confirm the hydrogen saturation on metal surface of Mg and Ni containing clusters. Linear structure of MgH2 was produced in the method of UM052X/6-311+G(d,p) calculation and a strong Mg-H antisymmetric stretching frequency mode at 1637 cm−1 and a bending frequency mode at 460 cm−1 is observed. The inclusion of ground-state Mg (3s1, S) into hydrogen experiences large energy barrier. Therefore, Mg atom exited to the electronic state of (3s3p, 1P) to overcome the energy barrier and insert into hydrogen. The vibration between 1000 and 1200 cm−1 is because of the bridge of Mg-H-Mg bond stretching. Outcome of calculation suggests that stronger interaction for H-Mg-H with hydrogen on metal surface is due to the strong H-H stretching modes induced by MgH2. The band at 1645 cm−1 is symmetrical stretching mode of H-Ni-H. Intermetallic ternary hydride Mg2Ni-(H2)2 (=Mg2NiH4) displayed the vibrational band (Figure 8) at 222, 528, 650, 7685, 1100, and 2030 cm−1 which includes the common band of Ni-(H2)n; here, n = 1 or 2 [26].

Different energy parameters have been investigated for Mg-(H2)n and Mg2NiH4 and shown in Table 2. Interaction energies between hydrogen and metal cluster increase with gradual number of captured hydrogen. Hydrogenation and interaction energies are increasing based on gradual addition of hydrogen. These phenomena correspond to the favorable hydrogenation on metals. Table 2 shows that adsorption energy values are higher than 50 kJ/mol and indicate chemisorption range hydrogenation (strongly bonded). According to the calculated data, (Table 2) interaction energy corresponds to the fact that the interaction between metals and hydrogen is electrostatic range. Magnesium shows higher affinity with hydrogen than nickel. Addition of nickel decreases strength of bonding between magnesium and hydrogen. The average desorption energies (Table 2) of the hydrides attribute that desorption may occur at higher condition (temperature, pressure) for the less adsorbed hydrogen and maximum captured hydrogen decompose at low temperature. Predicted hydrogen desoprtion energies (or decomposition energy) of MgH2 and Mg2NiH4 are 79.31 kJ/mol and 54.72 kJ/mol which is one of the important concerns for reversible hydrogenated materials and those are close to the calculated desorption energy value of 76 kJ/mol (for MgH2) by Swart et al., [27]. Reaction electrophilicity values of hydrides are decreasing upon increasing number of hydrogen captured. This phenomenon is very favorable for high hydrogenation [28].

The thermodynamic descriptors of hydride phases such as change of Gibb’s free energies (Δ), enthalpy change (Δ), and entropy change (Δ) were computed to ensure the feasibility of hydrogenation on metals surfaces. The negative (−) values of Gibb’s free energies (Δ), shown in Table 3, ensure the favorable hydrogenation phenomenon at room temperature and indicate the impetuous behavior of hydrogenation. DFT predicted values of hydrogenation enthalpy and entropy changes of hydrogenated clusters of MgH2 and Mg2NiH4 are −62.90 kJ/mol, −158 J/mol·K and −52.78 kJ/mol, −166 J/mol·K, respectively, that are close to reversible hydrogen storage materials thermodynamics values of −40 kJ/mol and −130 J/mol·K [6]. Negative enthalpy (Δ) of hydrides indicates that the hydrogenation is exothermic. Decreasing trend of entropy corresponds to the limitation of hydrogen movement on surface due to the less degrees of freedom of hydrogen. The reaction constant values are in the order of ( × 1013) that attribute to the slow hydrogenation on this materials surface.

4. Conclusion

The density functional theory (DFT) study has explained the mechanism of the hydrogenation and the hydrogenation properties of microstructure. Variation of the thermodynamics and energy descriptor values upon increasing captured hydrogen molecules correspond to the promising activities of the hydride microstructure and their feasible hydrogenation properties. Experiment and DFT based investigation demonstrated similar IR frequencies, which confirm the hydride formation of MgH2 and Mg2NiH4. Isothermal adsorption analysis predicted high hydrogenation capacity which is higher than the stoichiometric capacity of Mg2NiH4 (3.6 wt% hydrogen). Polycrystalline and flower-like interactive rough microstructure is an important material’s properties that can be formed in an optimal condition to enhance the hydrogenation.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors acknowledge gratefully the financial support of this study by the FRGS Grant 158-200-092, Malaysia, and the computational facility provided by University of Malaya.