Electronic band structure of 1D π-d hybridized narrow-gap metal-organic polymers.

One-dimensional (1D) metal-organic (MO) nanowires are captivating from fundamental and technological perspectives due to their distinctive magnetic and electronic properties. The solvent-free synthesis of such nanomaterials on catalytic surfaces provides a unique approach for fabricating low-dimensional single-layer materials with atomic precision and low amount of defects. A detailed understanding of the electronic structure of MO polymers such as band gap and dispersive bands is critical for their prospective implementation into nanodevices such as spin sensors or field-effect transistors. Here, we have performed the on-surface reaction of quinoidal ligands with single cobalt atoms (Co-QDI) on a vicinal Au(788) surface in ultra-high vacuum. This procedure promotes the growth and uniaxial alignment of Co-QDI MO chains along the surface atomic steps, while permitting the mapping of their electronic properties with space-averaging angle-resolved photoemission spectroscopy. In the direction parallel to the principal chain axis, a well-defined 1D band structure with weakly dispersive and dispersive bands is observed, confirming a pronounced electron delocalization. Low-temperature scanning tunneling microscopy/spectroscopy delves into the atomically precise structure of the nanowires and elucidates their narrow bandgap. These findings are supported with GW0 band structure calculations showing that the observed electronic bands emanate from the efficient hybridization of Co(3d) and molecular orbitals. Our work paves the way towards a systematic search of similar 1D π-d hybridized MO chains with tunable electronic and magnetic properties defined by the transition or rare earth metal atom of choice.

Role of the Magnetic Anisotropy in Atomic-Spin Sensing of 1D Molecular Chains.

One-dimensional metal-organic chains often possess a complex magnetic structure susceptible to modification by alteration of their chemical composition. The possibility to tune their magnetic properties provides an interesting playground to explore quasi-particle interactions in low-dimensional systems. Despite the great effort invested so far, a detailed understanding of the interactions governing the electronic and magnetic properties of the low-dimensional systems is still incomplete. One of the reasons is the limited ability to characterize their magnetic properties at the atomic scale. Here, we provide a comprehensive study of the magnetic properties of metal-organic one-dimensional (1D) coordination polymers consisting of 2,5-diamino-1,4-benzoquinonediimine ligands coordinated with Co or Cr atoms synthesized under ultrahigh-vacuum conditions on a Au(111) surface. A combination of integral X-ray spectroscopy with local-probe inelastic electron tunneling spectroscopy corroborated by multiplet analysis, density functional theory, and inelastic electron tunneling simulations enables us to obtain essential information about their magnetic structures, including the spin magnitude and orientation at the magnetic atoms, as well as the magnetic anisotropy.

Magnetic Properties of Metal⁻Organic Coordination Networks Based on 3d Transition Metal Atoms.

The magnetic anisotropy and exchange coupling between spins localized at the positions of 3d transition metal atoms forming two-dimensional metal⁻organic coordination networks (MOCNs) grown on a Au(111) metal surface are studied. In particular, we consider MOCNs made of Ni or Mn metal centers linked by 7,7,8,8-tetracyanoquinodimethane (TCNQ) organic ligands, which form rectangular networks with 1:1 stoichiometry. Based on the analysis of X-ray magnetic circular dichroism (XMCD) data taken at T = 2.5 K, we find that Ni atoms in the Ni⁻TCNQ MOCNs are coupled ferromagnetically and do not show any significant magnetic anisotropy, while Mn atoms in the Mn⁻TCNQ MOCNs are coupled antiferromagnetically and do show a weak magnetic anisotropy with in-plane magnetization. We explain these observations using both a model Hamiltonian based on mean-field Weiss theory and density functional theory calculations that include spin⁻orbit coupling. Our main conclusion is that the antiferromagnetic coupling between Mn spins and the in-plane magnetization of the Mn spins can be explained by neglecting effects due to the presence of the Au(111) surface, while for Ni⁻TCNQ the metal surface plays a role in determining the absence of magnetic anisotropy in the system.

Intermode Coupling Drives the Irreversible Tautomerization in Porphycene on Copper(111) Induced by Scanning Tunnelling Microscopy.

In this contribution, we develop a nonadiabatic theory that explains, from first-principles, the recently reported irreversible trans → cis tautomerization of porphycene on Cu(111) induced by a scanning tunnelling microscope at finite bias. The inelastic contribution to the STM current is found to excite a large number of skeletal vibrational modes of the molecule, thereby inducing a deformation of the potential energy landscape along the hydrogen transfer coordinate. Above a threshold bias, the stability of the tautomers is reversed, which indirectly drives the reaction via intermode coupling. The proposed potential deformation term accounts effectively for the excitation of all internal vibrational modes without increasing the dimensionality of the problem. The model yields information about reaction rates, explains the reaction irreversibility at low temperatures, and accounts for the presence of resonant processes.

On the tautomerisation of porphycene on copper (111): Finding the subtle balance between van der Waals interactions and hybridisation.

We use density-functional theory (DFT) to analyse the interaction of trans- and cis-porphycene with Cu(111) and their interconversion by intramolecular H-transfer. This tautomerisation reaction is characterised by small values for the reaction energy and barrier, on the order of ∼0.1 eV, where the trans configuration is thermodynamically more stable upon adsorption according to the experiments [J. N. Ladenthin et al., ACS Nano 9, 7287-7295 (2015)]. To gain even a qualitatively correct description of this reaction at the DFT level, an accurate treatment of dispersion interactions and a careful choice of the exchange contribution are required in order to predict the subtle energetics. Analysis of the electronic structure shows that adsorption is contributed by a van der Waals (vdW) interaction, mainly responsible for stabilising the polyaromatic fragments, and by a significant charge redistribution localised between Cu and the unsaturated N atoms of the molecule central cavity. We find that different vdW functionals can produce qualitatively different electronic structures, while yielding small trans vs. cis energy differences. Unlike other functionals surveyed here, vdW-DF with PBE exchange satisfactorily reproduces not only the experimental energetics but also the scanning tunneling microscopy images. This gives us confidence that this functional achieves a reliable balance between the two mechanisms contributing to the adsorption of porphycene.

Manipulating interfacial hydrogens at palladium via STM.

In this contribution, we provide a detailed dynamical analysis of the interfacial hydrogen migration mediated by scanning tunneling microscopy (STM). Contributions from the STM-current and from the non-adiabatic couplings are taken into account using only first principle models. The slight asymmetry of the tunnelling rates with respect to the potential bias sign inferred from experimental observations is reproduced by weighting the contributions of the metal acceptor-donor states for the propagation of the impinging electrons. The quasi-thermal inelastic collision mechanism is treated perturbatively. The influence of hydrogen pre-coverage is also investigated using new potential energy surfaces obtained from periodic density functional theory calculations. Fully quantum dynamical simulations of the system evolution are performed by solving the Pauli master equation, providing insight into the reaction mechanism of STM manipulation of subsurface hydrogens. It is observed that the hydrogen impurity favors resurfacing over occupation of the bulk and subsurface sites whenever possible. The present simulations give strong indication that the experimentally observed protuberances after STM-excitation are due to hydrogen accumulating in the vicinity of the surface.

Ab initio molecular dynamics calculations on scattering of hyperthermal H atoms from Cu(111) and Au(111).

Energy loss from the translational motion of an atom or molecule impinging on a metal surface to the surface may determine whether the incident particle can trap on the surface, and whether it has enough energy left to react with another molecule present at the surface. Although this is relevant to heterogeneous catalysis, the relative extent to which energy loss of hot atoms takes place to phonons or electron-hole pair (ehp) excitation, and its dependence on the system's parameters, remain largely unknown. We address these questions for two systems that present an extreme case of the mass ratio of the incident atom to the surface atom, i.e., H + Cu(111) and H + Au(111), by presenting adiabatic ab initio molecular dynamics (AIMD) predictions of the energy loss and angular distributions for an incidence energy of 5 eV. The results are compared to the results of AIMDEFp calculations modeling energy loss to ehp excitation using an electronic friction ("EF") model applied to the AIMD trajectories, so that the energy loss to the electrons is calculated "post" ("p") the computation of the AIMD trajectory. The AIMD calculations predict average energy losses of 0.38 eV for Cu(111) and 0.13-0.14 eV for Au(111) for H-atoms that scatter from these surfaces without penetrating the surface. These energies closely correspond with energy losses predicted with Baule models, which is suggestive of structure scattering. The predicted adiabatic integral energy loss spectra (integrated over all final scattering angles) all display a lowest energy peak at an energy corresponding to approximately 80% of the average adiabatic energy loss for non-penetrative scattering. In the adiabatic limit, this suggests a way of determining the approximate average energy loss of non-penetratively scattered H-atoms from the integral energy loss spectrum of all scattered H-atoms. The AIMDEFp calculations predict that in each case the lowest energy loss peak should show additional energy loss in the range 0.2-0.3 eV due to ehp excitation, which should be possible to observe. The average non-adiabatic energy losses for non-penetrative scattering exceed the adiabatic losses to phonons by 0.9-1.0 eV. This suggests that for scattering of hyperthermal H-atoms from coinage metals the dominant energy dissipation channel should be to ehp excitation. These predictions can be tested by experiments that combine techniques for generating H-atom beams that are well resolved in translational energy and for detecting the scattered atoms with high energy-resolution.

Co nanodot arrays grown on a GdAu2 template: substrate/nanodot antiferromagnetic coupling.

Controlling anisotropy and exchange coupling in patterned magnetic nanostructures is the key for developing advanced magnetic storage and spintronic devices. We report on the antiferromagnetic interaction between a Co nanodot array and its supporting GdAu2 nanotemplate that induces large anisotropy values in individual Co nanodots. In clear contrast with nonmagnetic Au substrates, GdAu2 triggers an earlier switch from out-of-plane anisotropy in monatomic high dots to in-plane when the dot height becomes biatomic.

Determination of the surface structure of CeO2(111) by low-energy electron diffraction.

We determine the atomic structure of the (111) surface of an epitaxial ceria film using low-energy electron diffraction (LEED). The 3-fold-symmetric LEED patterns are consistent with a bulk-like termination of the (111) surface. By comparing the experimental dependence of diffraction intensity on electron energy (LEED-I(V) data) with simulations of dynamic scattering from different surface structures, we find that the CeO2(111) surface is terminated by a plane of oxygen atoms. We also find that the bond lengths in the top few surface layers of CeO2(111) are mostly undistorted from their bulk values, in general agreement with theoretical predictions. However, the topmost oxygen layer is further from the underlying cerium layer than the true bulk termination, an expansion that differs from theoretical predictions.