Controlling the spin of metal atoms adsorbed on oxide surfaces: Ni on regular and defective sites of 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.

Towards understanding performance differences between approximate density functionals for spin states of iron complexes

Density functional theory has been widely used to investigate the structural and electronic properties of heme-containing proteins such as cytochrome P450. Nevertheless, recent studies have shown that approximate exchange-correlation energy density functionals can incorrectly predict the stability order of spin states in, for instance, iron-containing pyridine and imidazole systems. This raises questions about the validity of earlier theoretical studies. In this work, we systematically investigate a few typical inorganic and organic iron-containing complexes and try to understand the performance difference of various density functionals. Two oxidation states of iron, Fe(II) and Fe(III), with different spin states and both adiabatic and vertical structures are considered. A different description of the outmost molecular orbital is found to play the crucial role. Local density and generalized gradient based functionals bias the lower spin state and produce a more localized frontier orbital that is higher in energy than the hybrid functionals. Energy component analysis has been performed, together with comparison of numerous structural and electronic properties. Implications of the present work to the theoretical study of heme-containing biological molecules and other spin-related systems are discussed.

Single d-Metal Atoms on Fs and Fs+ Defects of MgO(001): A Theoretical Study across the Periodic Table

Single d-metal atoms on oxygen defects F(s) and F(s+) of the MgO(001) surface were studied theoretically. We employed an accurate density functional method combined with cluster models, embedded in an elastic polarizable environment, and we applied two gradient-corrected exchange-correlation functionals. In this way, we quantified how 17 metal atoms from groups 6-11 of the periodic table (Cu, Ag, Au; Ni, Pd, Pt; Co, Rh, Ir; Fe, Ru, Os; Mn, Re; and Cr, Mo, W) interact with terrace sites of MgO. We found bonding with F(s) and F(s+) defects to be in general stronger than that with O2- sites, except for Mn-, Re-, and Fe/F(s) complexes. In M/F(s) systems, electron density is accumulated on the metal center in a notable fashion. The binding energy on both kinds of O defects increases from 3d- to 4d- to 5d-atoms of a given group, at variance with the binding energy trend established earlier for the M/O2- complexes, 4d < 3d < 5d. Regarding the evolution of the binding energy along a period, group 7 atoms are slightly destabilized compared to their group 6 congeners in both the F(s) and F(s+) complexes; for later transition elements, the binding energy increases gradually up to group 10 and finally decreases again in group 11, most strongly on the F(s) site. This trend is governed by the negative charge on the adsorbed atoms. We discuss implications for an experimental detection of metal atoms on oxide supports based on computed core-level energies.

A systematic density functional theory study of the electronic structure of bulk and (001) surface of transition-metals carbides

A systematic study of the bulk and surface geometrical and electronic properties of a series of transition-metal carbides (TMC with TM = Ti, V, Zr, Nb, Mo, Hf, Ta, and W) by first-principles methods is presented. It is shown that in these materials the chemical bonding is strongly covalent, the cohesive energies being directly related to the bonding-antibonding gap although the shift of the center of the C(2s) band related peak in the density of states with respect to diamond indicates that some metal to carbon charge transfer does also take place. The (001) face of these metal carbides exhibits a noticeable surface rumpling which grows along the series. It is shown that neglecting surface relaxation results in very large errors on the surface energy and work function. The surface formation induces a significant shift of electronic energy levels with respect to the corresponding values in the bulk. The extent and nature of the shift can be understood from simple bonding-antibonding arguments and is enhanced by the structural rippling of this surface.

The adsorption and dissociation of Cl2 on the MgO (001) surface with vacancies: Embedded cluster model study

The adsorption of Cl(2) at a low-coordinated oxygen site (edge or corner site) and vacancy site (terrace, edge, corner F, F(+), or F(2+) center) has been studied by the density functional method, in conjunction with the embedded cluster models. First, we have studied the adsorption of Cl(2) at the edge and corner oxygen sites and the results show that Cl(2), energetically, is inclined to adsorb at the corner oxygen site. Moreover, similar to the most advantageous adsorption mode for Cl(2) on the MgO (001) perfect surface, the most favorable adsorption occurs when Cl(2) approaches the corner oxygen site along the normal direction. A small amount of electrons are transferred from the substrate to the antibonding orbital of the adsorbate, leading to the Cl-Cl bond strength weakened a little. Regarding Cl(2) adsorption at the oxygen vacancy site (F, F(+), or F(2+) center), both large adsorption energies and rather much elongation of the Cl-Cl bond length have been obtained, in particular at the corner oxygen vacancy site, with concurrently large amounts of electrons transferred from the substrate to the antibonding orbital of Cl(2). It suggests, at the oxygen vacancy site, that Cl(2) prefers to dissociate into Cl subspecies. And the potential energy surface indicates that the dissociation process of molecular Cl(2) to atomic Cl is virtually barrierless.