Adhesion energy of Cu atoms on the MgO(001) surface

Magnetic nature and hyperfine interactions of transition metal atoms adsorbed on ultrathin insulating films: a challenge for DFT

The magnetic ground state and the hyperfine coupling parameters of some first-row transition metal (TM) atoms (Ti, Cr, Mn, Fe, Co, and Ni) adsorbed on ultrathin insulating oxide films are studied by means of DFT calculations. The results obtained using GGA, screened hybrid, and GGA+U functionals are compared for TMs adsorbed on free-standing MgO(100). Then, the case of adsorption on MgO mono- and bilayers supported on Ag(100) is studied using GGA+U. Along with the problematic aspects inherent to the calculation of hyperfine coupling constants, a critical dependence on the magnetic state and electron configuration of the TM is reported, which implies a real challenge for the state-of-the-art DFT methods. In the cases where all functionals considered provide a coherent magnetic and electron configuration, however, the calculated hyperfine parameters do not depend significantly on the choice of the functional. In this respect, the role of the metal support in the hyperfine coupling constants is highly system-dependent and becomes crucial in all cases where the support modifies the oxidation state of the adatom, induces a change in the bonding site or simply induces a rearrangement of the orbital energy diagram. This has important implications for the modelling of single TM atoms deposited on insulating ultrathin films supported on metals for application in quantum technologies or as memory devices.

Efficient and accurate description of adsorption in zeolites

Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.

Physico-Chemical Insights into Gas-Phase and Oxide-Supported Sub-Nanometre AuCu Clusters

Catalysis by AuCu nanoclusters is a promising scientific field. However, our fundamental understanding of the underlying mechanisms of mixing in AuCu clusters at the sub-nanometre scale and their physico-chemical properties in both the gas-phase and on oxide supports is limited. We have identified the global minima of gas-phase and MgO(100)-supported AuCu clusters with 3–10 atoms using the Mexican Enhanced Genetic Algorithm coupled with density functional theory. Au and Cu adatoms and supported dimers have been also simulated at the same level of theory. The most stable composition, as calculated from mixing and binding energies, is obtained when the Cu proportion is close to 50%. The structures of the most stable free AuCu clusters exhibit Cu-core/Au-shell segregation. On the MgO surface however, there is a preference for Cu atoms to lie at the cluster-substrate interface. Due to the interplay between the number of interfacial Cu atoms and surface-induced cluster rearrangement, on the MgO surface 3D structures become more stable than 2D structures. The O-site of MgO surface is found to be the most favourable adsorption site for both metals. All dimers favour vertical (V) configurations on the surface and their adsorption energies are in the order: AuCu < CuCu < AuAu < AuCu (where the underlined atom is bound to the O-site). For both adatoms and AuCu dimers, adsorption via Cu is more favourable than Au-adsorbed configurations, but, this disagrees with the ordering for the pure dimers due to a combination of electron transfer and the metal-on-top effect. Binding energy (and second difference) and HOMO-LUMO gap calculations show that even-atom (even-electron) clusters are more stable than the neighbouring odd-atom (odd- electron) clusters, which is expected for closed- and open-shell systems. Supporting AuCu clusters on the MgO(100) surface decreases the charge transfer between Au and Cu atoms calculated in free clusters. The results of this study may serve as a foundation for designing better AuCu catalysts.

Modeling doped and defective oxides in catalysis with density functional theory methods: room for improvements

Due to the well-known problem of the self-interaction, standard density functional theory (DFT) methods tend to produce delocalized holes and electrons in defective oxide materials even when there is ample experimental evidence of a strong localization. For late transition metal compounds or rare earth oxides, this results in the incorrect description of the electronic structure of the system (e.g., magnetic insulators are predicted to be metallic). Practical ways to correct this deficiency are based on the use of hybrid functionals or of the DFT+U approach. In this way, most of the limitations related to the self-interaction are removed, and the electronic structure is properly described. What is less clear is to what extent hybrid functionals, DFT+U approaches, or standard DFT functionals can properly describe the strength of the chemical bonds at the surface of an oxide. This is a crucial question if one is interested in the catalytic properties of oxide surfaces. Oxidation reactions often involve oxygen detachment from the surface and incorporation into an organic substrate. Oxides are doped with heteroatoms to create defects and facilitate oxygen removal from the surface, with formation of oxygen vacancies. Do standard DFT calculations provide a good binding energy of the missing oxygen despite the failure in giving the right electronic structure? Can hybrid functionals or the DFT+U approach provide a simple yet reliable way to get accurate reaction enthalpies and energy barriers? In this essay, we discuss these problems by analyzing some case histories and the relatively scarce data existing in the literature. The conclusion is that while modern electronic structure methods accurately reproduce and predict a wide range of electronic, optical, and magnetic properties of oxides, the description of the strength of chemical bonds still needs considerable improvements.

Structure of Ag Clusters Grown on Fs-Defect Sites of an MgO(1 0 0) Surface

The structure of AgN clusters (N=1-4, 6, 8, 10), both in the gas phase and grown on the MgO(1 0 0) surface containing Fs-defects, has been investigated by a density functional basin-hopping (DF-BH) approach. In analogy with what observed in the case of gold clusters, it is found that the presence of the defect implies a double frustration and a cylindrical invariance of the metal-surface interaction, causing small Ag clusters growing around the Fs defect to be highly fluxional. Nevertheless, two different structural crossovers are found to be induced by the metal-defect interaction for the adsorbed clusters such that: 1) planar structures prevail for N<or=4 (as in the gas phase); 2) noncrystalline (fivefold symmetric) structures, which are the lowest energy ones in the gas phase for medium sized AgN clusters (N>or=7), prevail for N=6 and N=8; 3) distorted face-centered cubic (fcc) structures grown pseudomorphically on the defected surface prevail for N=10. The transition from fivefold to fcc motifs is rationalized in terms of the double-frustration effect, which increases the bond strain of the noncrystalline structures. Detrapping energies from the defect were also calculated; the lowest energy pathway corresponds to the detachment of a dimer.

Adsorption of carbon on Pd clusters of nanometer size: A first-principles theoretical study

Adsorbed atomic C species can be formed in the course of surface reactions and commonly decorate metal catalysts. We studied computationally C adsorption on Pd nanoclusters using an all-electron scalar relativistic density functional method. The metal particles under investigation, Pd(55), Pd(79), Pd(85), Pd(116), Pd(140), and Pd(146), were chosen as fragments of bulk Pd in the form of three-dimensional octahedral or cuboctahedral crystallites, exposing (111) and (100) facets as well as edge sites. These cluster models are shown to yield size-converged adsorption energies. We examined which surface sites of these clusters are preferentially occupied by adsorbed C. According to calculations, surface C atoms form strongly adsorbed carbide species (with adsorption energies of more than 600 kJ mol(-1)) bearing a significant negative charge. Surface sites allowing high, fourfold coordination of carbon are overall favored. To avoid effects of adsorbate-adsorbate interaction in the cluster models for carbon species in the vicinity of cluster edges, we reduced the local symmetry of selected adsorption complexes on the nanoclusters by lowering the global symmetry of the nanocluster models from point group O(h) to D(4h). On (111) facets, threefold hollow sites in the center are energetically preferred; adsorbed C is calculated to be slightly less stable when displaced to the facet borders.

Optical properties of Cu nanoclusters supported on MgO(100)

The vertical transitions of Cu atoms, dimers, and tetramers deposited on the MgO surface have been investigated by means of ab initio calculations based either on complete active space second-order perturbation theory or on time-dependent density functional theory. Three adsorption sites have been considered as representative of the complexity of the MgO surface: regular sites at flat (100) terraces, extended defects such as monoatomic steps, and point defects such as neutral oxygen vacancies (F or color centers). The optical properties of the supported Cu clusters have been compared with those of the corresponding gas-phase units. Upon deposition a substantial modification of the energy levels of the supported cluster is induced by the Pauli repulsion with the substrate. This causes shifts in the optical transitions going from free to supported clusters. The changes in cluster geometry induced by the substrate have a much smaller effect on the optical absorption bands. On F centers the presence of filled impurity levels in the band gap of MgO results in a strong mixing with the empty levels of the Cu atoms and clusters with consequent deep changes in the optical properties of the color centers. The results allow to interpret electron energy loss spectra of Cu atoms deposited on MgO thin films.

Bond energies for molecules, clusters, and deposit systems

A new approach is suggested to the assignment of bond energies in molecules and clusters. It uses a shareholder principle for the redistribution of the shifts in atomic energies, which arise in a molecule, on the bonds. The scheme is directly suitable for semiempirical methods, where only one- and two-center terms occur. MSINDO calculations are performed to demonstrate the suitability of the approach for molecules and clusters. As an application the bonding in a deposit system is analyzed for the case of copper on magnesium oxide. It is found that copper atoms do not only bind to the preferred oxygen sites but also substantially to the magnesium sites. The copper-copper bonds are the strongest and will determine the structure of copper clusters on magnesium oxide surfaces.

Metal adsorption and adhesion energies on MgO(100)

The calormetically measured heats of adsorption of Cu, Ag, and Pb on MgO(100), previously measured in our group, are correlated with bulk properties of the metals and their sticking probabilities and film morphologies. The low-coverage heats of adsorption (when the metals are mainly in two-dimensional (2D) islands) are used to estimate metal-MgO(100) bond energies within a pairwise bond additivity model. These values correlate well with the observed initial sticking probabilities and saturation island densities of the metals. This supports a transient mobile precursor model for adsorption. The values also correlate with their bulk sublimation energies, which suggests that covalent metal-Mg bonding dominates the interaction at low coverage, probably due to very strong bonding at defects. The heats of adsorption integrated up to multilayer coverages provide the metal-MgO(100) adhesion energies and metal-MgO(100) bond energies for metals in 3D films. These values correlate with the sum of magnitudes of the metal's bulk sublimation energy plus the heat of formation of the bulk oxide of the metal per mole of metal atoms. This suggests that local chemical bonds, both metal-oxygen and covalent metal-Mg, dominate the interfacial bonding for 3D films.

Theoretical analysis of the growth mode for thin metallic films on oxide substrates

We show how the growth mode of a thin metallic film on an insulating substrate can be predicted theoretically by combining thermodynamic considerations with ab initio calculations for ordered metal/insulator interfaces at low coverage. Our approach is illustrated by calculations for Ag film deposited on an MgO substrate. Ab initio calculations predict high mobility of adsorbed Ag atoms on MgO, even at low temperatures, which greatly aids their aggregation.