Bonding states of hydrogen for supported Ti clusters on pristine and defective graphene

To determine whether graphene-supported Ti clusters can synergistically store hydrogen through Kubas and spillover effect, we systematically investigate the growth pattern of Ti_n (n = 1–10) clusters on pristine and defective graphene and analyze multiple types of bonding states of hydrogen in detail. For pristine graphene, the most stable Ti_n (n = 1–10) clusters are the quasi-planar structure except for supported Ti₄ and Ti₅. The Ti dissociation energies and the binding energy of Ti_n clusters gradually increase with increasing n, which indicates that larger Ti_n clusters tend to form. Efficient spillover will occur on single-site Ti Catalysts at low hydrogen concentration due to the lower hydrogen spillover energy barrier (3.05 eV), while the energy barrier of hydrogen migration from Tin (n = 2–7) clusters to graphene on the cluster is 5.34–6.82 eV. The Ti: H ratio is a maximum of 1:8 for the single-site Ti catalyst, while decreases with the Tin cluster increases. Therefore, the pristine graphene-supported Ti nanoclusters are more suitable as substrates for hydrogen adsorbent rather than spillover. The introduced defects make Ti_n clusters have three-dimensional conical configurations from n = 4. Ti₃ and Ti₆ are the most stable clusters. Moreover, the migration energy barriers of H atoms on them decrease from 6.54 eV to 6.82 eV–4.32 eV and 3.42 eV, respectively. Our results explain recent experimental phenomena [Appl Phys Lett 2015; 106: 083,901. ACS Energy Lett 2022; 7 (7): 2297–2303] in depth at the molecular level.     

Effect of defects and dopants in graphene on hydrogen interaction in graphene-supported NaAlH4

Carbon-based materials have attracted great attention over the past few years due to their role as a support for sodium alanate improving the kinetics of H₂ release/uptake. Herein, we used graphene with defects and various dopants to simulate the carbon materials and performed a periodic density functional theory study on the impact of the modifications in the graphene substrates on the hydrogen interaction in, and hydrogen desorption from, the highly dispersed sodium alanate. Our results showed that the impact of various defects and dopants can be categorized in groups: (i) Pristine graphene and pentagon–heptagon (5–7) pairs defective graphene, as well as N- and S-doped graphene substrates show a weak interaction with the supported sodium alanate cluster, as reflected in the geometry change of the supported cluster and charge transfer between the supported cluster and the substrate. These defects and dopants do not promote H₂ formation and desorption. (ii) Carbon vacancies, as well as B and Cl dopants, cause instantaneous H₂ formation in supported NaAlH₄ upon relaxation. (iii) O-, P-, F- and OH-doped graphene substrates led to the formation of a meta-stable di-hydrogen state with a H–H distance of ∼0.96 Å. There is an activation barrier between the meta-stable di-hydrogen state and the most stable state with H₂ being formed. Furthermore, our results with the optB88-vdW functional show that van der Waals interaction strengthens the binding of the cluster on the substrates by 0.9–1.4 eV over the PBE results but does not alter the relative stability of the system.     

Lithium clusters on graphene surface and their ability to adsorb hydrogen molecules

This research describes the theoretical study of the adsorption of lithium clusters on graphene and the ability to capture hydrogen molecules. The results of the studied structures showed that the [Li₁C₅₄H₁₈]⁺ system is capable of accepting three hydrogen molecules showing adsorption energies of 0.12 eV. On the other hand, it is important to note that in [Li_nC₅₄H₁₈] n = 2–6 systems, the lithium atoms that do not interact with the graphene surface, they can adsorb up to four hydrogen molecules. The [Li₆C₅₄H₁₈]4H₂ system presented a higher adsorption energy value of 0.31 eV. Finally, the Li–H₂ interactions were characterized by a NBO analysis, which showed that hydrogen atoms are the donors and lithium atoms are the acceptors.     

Theoretical study of hydrogen adsorption on nitrogen doped graphene decorated with palladium clusters

Density functional theory was used to study the adsorption of hydrogen on small palladium clusters (Pd_n, n = 1–4) supported on pyridine-like nitrogen doped graphene. Charge transfer and strong binding (up to four times higher than binding energy of Pd cluster supported on graphene) between graphene–nitrogen and palladium clusters prevent detachment of clusters and leads to three types of adsorption states of hydrogen. The first type is a molecular hydrogen physisorbed, the second one is an activated state of H₂ without adsorption barriers where H–H bond is relaxed and the third type is dissociated state. In dissociated cases, we found barriers below 0.56 eV. This means the process might occur spontaneously at room temperature. We also show that metal–N–vacancy complexes are very stable and involve nitrogen sp² and p_z, carbon p_z and palladium d orbitals. Besides, these systems exhibit an interesting magnetic behavior.     

Revealing the interfacial phenomenon of metal-hydride formation and spillover effect on C48H16 sheet by platinum metal catalyst (Pt (n=1-4)) for hydrogen storage enhancement

Hydrogen is a worldwide green energy carrier, however due its low storage capacity, it has yet to be widely used as an energy carrier. Therefore, the quantum chemical method is being employed in this investigation for better understand the hydrogen storage behaviour on Pt _(n = 1-4) cluster decorated C₄₈H₁₆ sheet. The Pt_(n = 1-4) clusters are strongly bonded on the surface of C₄₈H₁₆ sheet with binding energies of ?3.06, ?4.56, ?3.37, and ?4.03 eV respectively, while the charge transfer from Pt_(n = 1-4) to C₄₈H₁₆ leaves an empty orbital in Pt atom, which will be crucial for H₂ adsorption. Initially, the molecular hydrogen is adsorbed on Pt_(n = 1-4) decorated C₄₈H₁₆ sheet through the Kubas interaction with adsorption energies of ?0.85, ?0.66, ?0.72, and ?0.57 eV respectively, while H–H bond is elongated due to the transfer of electron from σ (HH) orbital to unfilled d orbital of the Pt atom, resulting in a Kubas metal-dihydrogen complexes. Furthermore, the dissociative hydrogen atoms adsorbed on Pt_(n = 1-4) decorated C₄₈H₁₆ sheet have adsorption energies of ?1.14 eV, ?1.02 eV, ?0.95 eV, and ?1.08 eV, which are greater than the molecular hydrogen adsorption on Pt_(n = 1-4) cluster supported C₄₈H₁₆ sheet with lower activation energy of 0.007, 0.109, 0.046, and 0.081 eV respectively. To enhance the dissociative hydrogen adsorption energy, positive and negative external electric fields are applied in the charge transfer direction. Increasing the positive electric field makes H–H bond elongation and good adsorption, whereas increasing the negative electric field results H–H bond contraction and poor adsorption. Thus, by applying a sufficient electric field, the H₂ adsorption and desorption processes are can be easily tailored.     

Theoretical calculation of hydrogen desorption energies of calcium hydride clusters

Recently calcium hydride has attracted attention as a possible component in ternary complex hydrides such as Ca(AlH₄)₂, Ca₂SiH_x and quaternary complex hydrides of the type Li–B–Ca–H. Although in bulk form CaH₂ decomposes reversibly above 600° centigrade we were motivated to see whether calcium hydride in cluster form has properties suitable for hydrogen storage. We report here the results of DFT calculations using VASP^® package for clusters Ca_nH_2n with n = 1, 2, 3, 4, 6, 8, 10, 12, 14, 16, 20 to get the ground state geometries, energies, bond lengths, and desorption energies, after molecular dynamics optimization. The desorption energy vs. cluster size n curve showed that the desorption energy goes up sharply to ∼1.4 eV per H₂ for n up to 4, followed by a broad maximum of ∼1.8 eV per H₂ around n = 12–14, and then tapers off to a nearly constant value of 1.6 eV per H₂ approximating bulk behavior, which compares favorably with previously reported results. Comparison of these results with those of Mg_nH_2n shows that Ca_nH_2n has a lesser potential as a hydrogen storage medium.     

First principles study on hydrogen storage in yttrium doped graphyne: Role of acetylene linkage in enhancing hydrogen storage

Through Density Functional Theory Simulations we predict that a Ytrrium atom attached on graphyne surface can adsorb up to a maximum of 9 molecular hydrogens (H₂), with a uniform binding energy of ∼0.3 eV/H₂ and an average desorption temperature of around 400 K (ideal for fuel cell applications), leading to 10 wt% of hydrogen, substantially higher than the requirement by DoE. The higher hydrogen wt% in Y doped graphyne compared to Y doped Single Walled Carbon Nanotubes (SWNT) and graphene is due to the presence of sp hybridized C atoms (in the acetylene linkage) supplying additional in-plane p_x-p_y orbitals leading to π (π*) bonding (antibonding) states. Charge transfer from metal to carbon nanostructure results in a redistribution of s, p, d orbitals of the metal leading to a non - spin polarized ground state in Y doped graphyne, due to the presence of the acetylene linkage, whereas Y doped SWNT and graphene remain magnetic like the isolated metal atom. In the non-magnetic graphyne + Y system, the net charge transfer from Y to successive H₂ molecules is less than in magnetic Y + graphene and Y + SWNT systems, enabling Y + graphyne to store a larger number of H₂ molecules. Furthermore, our ab initio MD simulations show that the system is stable even at room temperature and there is no dissociation of H₂ molecules, enabling the system to achieve 100% desorption. So Y doped graphyne is found to be a promising hydrogen storage device with high wt%, 100% recyclability and desirable desorption temperature.     

Reversible hydrogen absorption in sodium intercalated fullerenes

The hydrogen absorption of sodium intercalated fullerenes (Na_xC₆₀) was determined and compared to pure fullerenes (C₆₀). Up to 3.5 mass% hydrogen can reversibly be absorbed in Na_xC₆₀ at 200 °C and a hydrogen pressure of 200 bar. The absorbed amount of hydrogen is significantly higher than for the case when only the sodium would be hydrogenated (∼1 mass% for x = 10). At 200 bar the onset of hydrogen absorption is observed at 150 °C. At a pressure of 1 bar hydrogen the major desorption starts at 250 °C and is completed at 300 °C (heating rate 1 °C min⁻¹). This absorption and desorption temperatures are significantly reduced compared to pure C₆₀, either due to a catalytic reaction of hydrogen on sodium or due to the negatively charged C₆₀. The hydrogen ab/desorption is accompanied by a partial de/reintercalation of sodium. A minor part of the hydrogen is ionically bonded in NaH and the major part is covalently bonded in C₆₀H_x. The sample can be fully dehydrogenated and no NaH is left after desorption. In contrast to C₆₀, where the fullerene cages for high hydrogen loadings are destroyed during the sorption process, the Na_xC₆₀ sample stays intact. The samples were investigated by X-ray, in-situ neutron powder diffraction and infrared spectroscopy. Na_xC₆₀ was synthesized by reacting sodium azide (NaN₃) with C₆₀ (molar ratio of Na:C₆₀ is 10:1).     

Structure of V-doped Pdn (n = 2–12) clusters and their ability for H2 dissociation

The structure and stability of V-doped Pd_n (n = 2–12) clusters as well as their ability for hydrogen dissociation are analyzed using a successive growth algorithm coupled with density functional theory (DFT) computations. From the structural point of view, the lowest energy structures of these clusters are three-dimensional with exohedral geometries for n = 2–7 whereas endohedral for n = 8 onward. From their second-order energy differences, Pd₄V and Pd₁₀V are found to be the most stable ones. Among the Pd_nV(H₂) complexes, Pd₆V2H possesses the highest stability, as it is supported by the chemisorption energy, the vertical ionization potential (VIP), and the vertical electron affinity (VEA), respectively. Most importantly, the hydrogen dissociation pathway on Pd_nV clusters with n = 3, 4 and 10–12 shows that these clusters are rigid and suitable to dissociate H₂ while for n = 5–9 the structure of the clusters changes. The H₂ dissociation process on Pd_nV clusters with n = 8, 10, and 11 carries out barrierless.     

Small Pd cluster adsorbed double vacancy defect graphene sheet for hydrogen storage: A first-principles study

Stability and electronic properties of small Pd_n clusters (n = 1–5), adsorbed on different types of double vacancy (DV) defect graphene sheets are thoroughly investigated by both density functional theory (DFT) and molecular dynamics (MD). Defect bridge sites of DV(555-777) defect graphene sheet are identified to be the most favorable for Pd₄ cluster adsorption. MD calculations, performed using a canonical ensemble, showed this system to be highly stable up to 800 K. Much better hybridization between C 2p and Pd 4d and 5s orbitals near Fermi level as well as higher charge transfer to graphene sheet was found to be the governing reason for enhanced stability of Pd₄ cluster on DV(555-777) defect site. Comparative analysis of H₂ storage on Pd₄ cluster adsorbed pristine and DV(555-777) defect graphene sheet showed, while adsorption energy/H₂ molecule for both cases lie well within desirable energy window for a hydrogen storage media, the later is much more efficient energetically as distorted in plane sp² hybridization reduces the saturations of C–C bonds in the defect regions, making more electron density available for bonding; which leads to higher net charge gain of Pd₄ cluster and higher charge sharing with H₂ molecule.     

Li-decorated double vacancy graphene for hydrogen storage application: A first principles study

Lithium decoration is an effective strategy for improving the hydrogen adsorption binding energy and the storage capacity in carbon nanostructures. Here, it is shown that Li-decorated double carbon vacancy graphene (DVG) can be used as an efficient hydrogen storage medium by means of Density Functional Theory (DFT) based calculations. The Li binding energy in DVG is 4.04 eV, which is much higher than that of pristine graphene. A maximum of four hydrogen molecules adsorb on Li decorated on one side of DVG and this leads to a gravimetric storage capacity of 3.89 wt% with an average adsorption binding energy of 0.23 eV/H₂. When Li is decorated on both sides of DVG, the gravimetric storage capacity reaches 7.26 wt% with a binding energy of 0.26 eV/H₂ which shows that desorption would take place at ambient conditions.     

Hydrogen storage in lithium-decorated benzene complexes

Based on first-principles calculations, we find Li-decorated benzene complexes are promising materials for high-capacity hydrogen storage. Lithium atoms in the complexes are always positively charged, and each one can bind at most four H₂ molecules by a polarization mechanism. Therefore, a hydrogen uptake of 8.6 wt% and 14.8 wt% can be achieved in isolated C₆H₆–Li and Li–C₆H₆–Li complexes, respectively. The binding energy in the two cases is 0.22 eV/H₂ and 0.29 eV/H₂, respectively, suitable for reversible hydrogen storage. Various dimers may form, but the hydrogen storage capacity remains high or uninfluenced. This study provides not only a promising hydrogen storage medium but also comprehensions to other hydrogen storage materials containing six-carbon rings.     

Hydrogen flux of BCC and FCC PdCuAg membranes

First-principles calculations have been used to study the effects of Ag addition on adsorption and dissociation of H₂ on BCC and FCC PdCu surfaces as well as hydrogen diffusion and recombinative hydrogen desorption through the PdCu membranes. It is found that the Ag addition makes it energetically difficult for the adsorption of H₂ on PdCu surfaces and hydrogen diffusion through PdCu, while could help the recombinative desorption of H atoms from both BCC (110) and FCC (111) surfaces of PdCu. Moreover, substitution of Ag for Pd or Cu would impede or improve the dissociation of H₂ on PdCu surface. Calculations also reveal that the overall hydrogen flux of BCC Pd₈Cu₈ (Pd₈Cu₇Ag) membranes is determined by the recombinative desorption and diffusion when the membrane thickness is smaller and bigger than 10.31 (4.73) μm, respectively. In addition, hydrogen diffusion is the dominant step of hydrogen permeation of FCC PdCu and PdCuAg as well as BCC Pd₇Cu₈Ag membranes. The present results not only agree well with experimental observations in the literature, but also deepen the understanding of the effect of Ag alloying on hydrogen permeation through PdCu membranes.     

A study on the hydrogen storage performance of graphene–Pd(T)–graphene structure

The 1–6 H₂ molecule adsorption energy and electronic properties of sandwich graphene–Pd(T)–Graphene (G–Pd(T)–G) structure were studied by the first-principle analysis. The binding energies, adsorption energies, and adsorption distances of Pd atoms-modified single-layer graphene and bilayer graphene with H₂ molecules at B, H, T adsorption sites were calculated. In bilayer graphene, the adsorption properties at T sites were found to be more stable than those at B and H sites. The binding energy of Pd atoms (4.16 eV) on bilayer graphene was higher than the experimental cohesion energy of Pd atoms (3.89 eV), and this phenomenon eliminated the impact of metal clusters on adsorption properties. It was found that three H₂ molecules were stably adsorbed on the G–Pd(T)–G structure with an average adsorption energy of 0.22 eV. Therefore, it can be speculated that G–Pd(T)–G is an excellent hydrogen storage material.     

Kinetic improvement of H2 absorption and desorption properties in Mg/MgH2 by using niobium ethoxide as additive

The hydrogen absorption and desorption properties of a MgH₂ – 1 mol.% Nb(V) ethoxide mixture are reported. The material was prepared by hand mixing the additive with previously ball-milled MgH₂. Nb ethoxide reacts with MgH₂ during heating, releasing C₂H₆ and H₂, and producing MgO and Nb or Nb hydride. Hydriding and dehydriding are greatly enhanced by the use of the alkoxide. At 250 °C the material with Nb takes up 1.8 wt% in 30 s compared with 0.1 wt% of pure Mg, and releases 4.2 wt% in 30 min, whereas MgH₂ without Nb does not appreciably desorb hydrogen. The absorption and desorption activation energies are reduced from 153 kJ/mol H₂ to 94 kJ/mol H₂, and from 176 kJ/mol H₂ to 75 kJ/mol H₂, respectively. The hydrogen sorption properties remain stable after 10 cycles at 300 °C. The kinetic improvement is attributed to the fine distribution of amorphous/nanometric NbH_x achieved by the dispersion of the liquid additive.     

Hydrogen release mechanisms of MgH2 over NiN4-embedded graphene nanosheet: First-principles calculations

In this study, first-principles calculations were performed to investigate the catalytic effect of NiN4-G on the dehydrogenation of MgH₂. Side-on MgH₂ can be adsorbed stably on the NiN4-G monolayer and is preferentially adsorbed on the NiN4 site compared with the graphene site. The hydrogen desorption process, in which MgH₂ dissociated to the Mg atom on the NiN4 site or graphene site and an H₂ molecule in the vacuum, should overcome lower barriers than pure MgH₂. This is because the corresponding Mg–H bond is weakened owing to the electron transfer between the Mg atom and the substrate. The hydrogen desorption enthalpies of the (MgH₂)₅ cluster on the NiN4 active and graphene sites are significantly smaller (0.11 eV and 1.50 eV, respectively) when H2+H2 is released from the cluster compared with those of the undoped MgH₂ cluster (2.48 eV). Therefore, the NiN4-G monolayer can provide the double effect of the NiN4 active and graphene sites on improving the dehydrogenation performance of MgH₂.     

First-principles based modeling of hydrogen permeation through Pd–Cu alloys

The solubility and diffusivity of hydrogen in disordered fcc Pd_1−xCu_x alloys are investigated using a combination of first-principles calculations, a composition-dependent local cluster expansion (CDLCE) technique, and kinetic Monte Carlo simulations. We demonstrate that a linear CDCLE model can accurately describe interstitial H in fcc Pd_1−xCu_x alloys over the entire composition range (0 ≤ x ≤ 1) with accuracy comparable to that of direct first-principles calculations. Our predicted H solubility and permeability results are in reasonable agreement with experimental measurements. The proposed model is quite general and can be employed to rapidly and accurately screen a large number of alloy compositions for potential membrane applications. Extension to ternary or higher-order alloy systems should be straightforward. Our study also highlights the significant effect of local lattice relaxations on H energetics in size-mismatched disordered alloys, which has been largely overlooked in the literature.     

ZrO2@Nb2CTx composite as the efficient catalyst for Mg/MgH2 based reversible hydrogen storage material

While Mg/MgH₂ system has a high hydrogen storage capacity, its sluggish hydrogen desorption rate has hindered practical applications. Herein, we report that the hydrogen absorption and desorption kinetics of Mg/MgH₂ system can be significantly improved by using the synergetic effect between Nb₂CT_x MXene and ZrO₂. The catalyst of Nb₂CT_x MXene loading with ZrO₂ (ZrO₂@Nb₂CT_x) is successfully synthesized, and the dehydrogenation activation energy of MgH₂ becomes as low as 60.0 kJ/mol H₂ when ZrO₂@Nb₂CT_x is used as the catalyst, which is far smaller than the case of ZrO₂ (94.8 kJ/mol H₂) and Nb₂CT_x MXene (125.6 kJ/mol H₂). With the addition of ZrO₂@Nb₂CT_x catalyst, MgH₂ can release about 6.24 wt.% and 5.69 wt.% of hydrogen within 150 s at 300 °C and within 900 s even at 240 °C, respectively. Moreover, it realizes hydrogen absorption at room temperature, which can uptake 2.98 wt.% of hydrogen within 1800 s. The catalytic mechanism analysis demonstrates that the in-situ formed nanocomposites can weaken the Mg–H bonding and provide more hydrogen diffusion channels, enabling the dissociation and recombination of hydrogen under milder reaction conditions.     

Ti decorated heterocyclic rings for hydrogen storage

This study uses first-principles calculations to investigate and compare the hydrogen storage properties of Ti doped benzene (C₆H₆Ti) and Ti doped borazine (B₃N₃H₆Ti) complexes. C₆H₆Ti and B₃N₃H₆Ti complex each can adsorb four H₂ molecules, but the former has a 0.11 wt% higher H₂ uptake capacity than the latter. Ti atoms bind to C₆H₆ more strongly than B₃N₃H₆. The hydrogen adsorption energies with Gibbs free energy correction for C₆H₆Ti and B₃N₃H₆Ti complexes are 0.17 and 0.45 eV, respectively, indicating reversible hydrogen adsorption. The hydrogen adsorption properties of C₆H₆Ti have also been studied after boron (B) and nitrogen (N) atom substitutions. Several B and N substituted structures between C₆H₆Ti and B₃N₃H₆Ti with different boron and nitrogen concentration and at different positions were considered. Initially, one boron and one nitrogen atom is substituted for two carbon atoms of benzene at three different positions and three different structures are obtained. Seven structures are possible when four carbon atoms of benzene are replaced by two boron and two nitrogen atoms at different positions. The hydrogen storage capacity of the C₆H₆Ti complex increases as boron and nitrogen atom concentrations increases. The positions of substituted boron and nitrogen atoms have less impact on H₂ uptake capacity for the same B and N concentration. The position and concentration of B and N affects the H₂ adsorption energy as well as the temperature and pressure range for thermodynamically favorable H₂ adsorption. The H₂ desorption temperature for all the complexes is found to be higher than 250 K indicates the stronger binding of H₂ molecules with these complexes.