Gj Gert Jan Kramer

The influence of electronic structure on hydrogen absorption in palladium alloys

The influence of the electronic structure and the lattice constant on hydrogen absorption in bulk Pd3M1 (M = Cd, Ag, Au, Pd, Cu, Ni, Pt, Pb, Sn, Fe, Rh, Ru) has been studied by density-functional calculations. We have assumed face-centred cubic structure for all the alloys, and hydrogen has been placed in the octahedral site surrounded by six Pd atoms. We have calculated the absorption energy of hydrogen in the alloys, and found that the influence of the electronic structure is much more important than that of the lattice constant. The results demonstrate that Miedemas empirical rule is also satisfied in this system, i.e., the higher the binding energy of the host alloy, the less stable the hydride. We have also calculated the detailed electronic structures of the alloys and their hydrides. We found that more stable hydrogen absorption is correlated with the hydrogen 1s electrons, palladium s electrons, palladium s-like electrons and the palladium d electrons moving higher in energy towards the Fermi level. The two latter relations have previously been described for bulk systems and surfaces respectively, while this study is apparently the first to point out the correlation between the position of the hydrogen band and the stability of the hydride, i.e., the deeper the hydrogen band, the less stable the hydride.

Predictions of melting, crystallization, and local atomic arrangements of aluminum clusters using a reactive force field

A parametrized reactive force field model for aluminum ReaxFF(Al) has been developed based on density functional theory (DFT) data. A comparison has been made between DFT and ReaxFF(Al) outputs to ascertain whether ReaxFF(Al) is properly parametrized and to check if the output of the latter has correlation with DFT results. Further checks include comparing the equations of state of condensed phases of Al as calculated from DFT and ReaxFF(Al). There is a good match between the two results, again showing that ReaxFF(Al) is correctly parametrized as per the DFT input. Simulated annealing has been performed on aluminum clusters Al(n) using ReaxFF(Al) to find the stable isomers of the clusters. A plot of stability function versus cluster size shows the existence of highly stable clusters (magic clusters). Quantum mechanically these magic clusters arise due to the complete filling of the orbital shells. However, since force fields do not care about electrons but work on the assumption of validity of Born-Oppenheimer approximation, the magic clusters are therefore correlated with high structural symmetry. There is a rapid decline in surface energy contribution due to the triangulated nature of the surface atoms leading to higher coordination number. The bulk binding energy is computed to be 76.8 kcal/mol. This gives confidence in the suitability of ReaxFF for studying and understanding the underlying dynamics in aluminum clusters. In the quantification of the growth of cluster it is seen that as the size of the clusters increase there is preference for the coexistence of fcc/hcp orders at the expense of simple icosahedral ordering, although there is some contribution from distorted icosahedral ordering. It is found that even for aluminum clusters with 512 atoms distorted icosahedral ordering exists. For clusters with N>/=256 atoms fcc ordering dominates, which implies that at this point we are already on the threshold of bulklike bonding.

Parametrization of a reactive force field for aluminum hydride

A reactive force field, REAXFF, for aluminum hydride has been developed based on density functional theory (DFT) derived data. REAXFF(AlH(3)) is used to study the dynamics governing hydrogen desorption in AlH(3). During the abstraction process of surface molecular hydrogen charge transfer is found to be well described by REAXFF(AlH(3)). Results on heat of desorption versus cluster size show that there is a strong dependence of the heat of desorption on the particle size, which implies that nanostructuring enhances desorption process. In the gas phase, it was observed that small alane clusters agglomerated into a bigger cluster. After agglomeration molecular hydrogen was desorbed from the structure. This thermodynamically driven spontaneous agglomeration followed by desorption of molecular hydrogen provides a mechanism on how mobile alane clusters can facilitate the mass transport of aluminum atoms during the thermal decomposition of NaAlH(4).

Modeling the sorption dynamics of NaH using a reactive force field

We have parametrized a reactive force field for NaH, ReaxFF(NaH), against a training set of ab initio derived data. To ascertain that ReaxFF(NaH) is properly parametrized, a comparison between ab initio heats of formation of small representative NaH clusters with ReaxFF(NaH) was done. The results and trend of ReaxFF(NaH) are found to be consistent with ab initio values. Further validation includes comparing the equations of state of condensed phases of Na and NaH as calculated from ab initio and ReaxFF(NaH). There is a good match between the two results, showing that ReaxFF(NaH) is correctly parametrized by the ab initio training set. ReaxFF(NaH) has been used to study the dynamics of hydrogen desorption in NaH particles. We find that ReaxFF(NaH) properly describes the surface molecular hydrogen charge transfer during the abstraction process. Results on heat of desorption versus cluster size shows that there is a strong dependence on the heat of desorption on the particle size, which implies that nanostructuring enhances desorption process. To gain more insight into the structural transformations of NaH during thermal decomposition, we performed a heating run in a molecular dynamics simulation. These runs exhibit a series of drops in potential energy, associated with cluster fragmentation and desorption of molecular hydrogen. This is consistent with experimental evidence that NaH dissociates at its melting point into smaller fragments.

Atomistic mechanisms for the (1×1)⇌hex surface phase transformations of Pt(100)

The atomistic mechanisms and dynamics of the (1 x 1) --> hex surface reconstruction of the clean Pt(100) and its lifting upon CO adsorption have been studied by means of both density functional theory and molecular dynamics simulations based on the modified embedded-atom method. It was found that during the surface reconstruction from the square (1 x 1) to the approximately 20% more dense (quasi)hexagonal phase, the required extra atoms are extracted from the second layer, leaving highly mobile subsurface vacancies, even in the presence of adatom islands or steps. These vacancies will, by surface diffusion, coalesce to form steps. In contrast, during the reverse process--where the hex surface reconstruction is lifted upon adsorption of CO--the roughly 20% excess surface atoms are in a collective manner ejected to form chains of adatoms. In turn, these adatoms coalesce into islands and steps by surface diffusion. A result of the two totally different mechanisms is that the reconstruction process can be concluded not to be the reverse of the deconstruction process.

Multiscale modeling of interaction of alane clusters on Al"111… surfaces: A reactive force field and infrared absorption spectroscopy approach

We have used reactive force field (ReaxFF) to investigate the mechanism of interaction of alanes on Al(111) surface. Our simulations show that, on the Al(111) surface, alanes oligomerize into larger alanes. In addition, from our simulations, adsorption of atomic hydrogen on Al(111) surface leads to the formation of alanes via H-induced etching of aluminum atoms from the surface. The alanes then agglomerate at the step edges forming stringlike conformations. The identification of these stringlike intermediates as a precursor to the bulk hydride phase allows us to explain the loss of resolution in surface IR experiments with increasing hydrogen coverage on single crystal Al(111) surface. This is in excellent agreement with the experimental works of Go et al. [E. Go, K. Thuermer, and J. E. Reutt-Robey, Surf. Sci. 437, 377 (1999)]. The mobility of alanes molecules has been studied using molecular dynamics and it is found that the migration energy barrier of Al(2)H(6) is 2.99 kcal/mol while the prefactor is D(0)=2.82 x 10(-3) cm(2)/s. We further investigated the interaction between an alane and an aluminum vacancy using classical molecular dynamics simulations. We found that a vacancy acts as a trap for alane, and eventually fractionates/annihilates it. These results show that ReaxFF can be used, in conjunction with ab initio methods, to study complex reactions on surfaces at both ambient and elevated temperature conditions.

Modeling of Hydrogen Storage Materials: A Reactive Force Field for NaH

Parameterization of a reactive force field for NaH is done using ab initio derived data. The parameterized force field(ReaxFFNaH) is used to study the dynamics governing hydrogen desorption in NaH. During the abstraction process of surface molecular hydrogen charge transfer is found to be well described by the parameterized force field. To gain more insight into the mechanism governing structural transformation of NaH during thermal decomposition a heating run in a molecular dynamics simulation is done. The result shows that a clear signature of hydrogen desorption is the fall in potential energy surface during heating.