The shortcomings of semi-local and hybrid functionals: what we can learn from surface science studies

A study of the adsorption of CO on late 4d and $5d$ transition metal (111) surfaces (Ru, Rh, Pd, Ag, Os, Ir, and Pt) considering atop and hollow site adsorption is presented. The applied functionals include the gradient corrected PBE and BLYP functional, and the corresponding hybrid Hartree-Fock density functionals HSE and B3LYP. We find that PBE based hybrid functionals (specifically HSE) yield, with the exception of Pt, the correct site order on all considered metals, but they also considerably overestimate the adsorption energies compared to experiment. On the other hand, the semi-local BLYP functional and the corresponding hybrid functional B3LYP yield very satisfactory adsorption energies and the correct adsorption site for all surfaces. We are thus faced with a Procrustean problem: the B3LYP and BLYP functionals seem to be the overall best choice for describing adsorption on metal surfaces, but they simultaneously fail to account well for the properties of the metal, vastly overestimating the equilibrium volume and underestimating the atomization energies. Setting out from these observations, general conclusions are drawn on the relative merits and drawbacks of various semi-local and hybrid functionals. The discussion includes a revised version of the PBE functional specifically optimized for bulk properties and surface energies (PBEsol), a revised version of the PBE functional specifically optimized to predict accurate adsorption energies (rPBE), as well as the aforementioned BLYP functional. We conclude that no semi-local functional is capable to describe all aspects properly, and including non-local exchange also only improves some, but worsens other properties.Comment: 12 pages, 6 figures; to be published in New Journal of Physic

Modeling the sorption dynamics of NaH using a reactive force field

We have parametrized a reactive force field for NaH, ReaxFFNaH, against a training set of ab initio derived data. To ascertain that ReaxFFNaH is properly parametrized, a comparison between ab initio heats of formation of small representative NaH clusters with ReaxFFNaH was done. The results and trend of ReaxFFNaH 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 ReaxFFNaH. There is a good match between the two results, showing that ReaxFFNaH is correctly parametrized by the ab initio training set. ReaxFFNaH has been used to study the dynamics of hydrogen desorption in NaH particles. We find that ReaxFFNaH 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

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

A parametrized reactive force field model for aluminum ReaxFFAl has been developed based on density functional theory (DFT) data. A comparison has been made between DFT and ReaxFFAl outputs to ascertain whether ReaxFFAl 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 ReaxFFAl. There is a good match between the two results, again showing that ReaxFFAl is correctly parametrized as per the DFT input. Simulated annealing has been performed on aluminum clusters Aln using ReaxFFAl 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

The role of interfacial energy in zeolite synthesis

The thermodn. factors playing a role in the crystn. of Si-rich zeolites are discussed. The Ostwald ripening model for nucleation and the Ostwald law of successive transformations make it possible to describe the kinetic phenomena in terms of the relative stability of the intermediate phases. The contribution of the interfacial free energy between the large internal surface area in the micropores and the intracryst. liquor to the chem. potential of the zeolitic material plays a crucial role. A model for the elementary steps occurring in pentasil zeolite formation is propose