Search

We extensively characterize the electronic structure of ultranarrow graphene nanoribbons (GNRs) with armchair edges and zigzag termini that have five carbon atoms across their width (5-AGNRs), as synthesized on Au(111). Scanning tunneling spectroscopy measurements on the ribbons, recorded on both the metallic substrate and a decoupling NaCl layer, show well-defined dispersive bands and in-gap states. In combination with theoretical calculations, we show how these in-gap states are topological in nature and localized at the zigzag termini of the nanoribbons. In addition to rationalizing the driving force behind the topological class selection of 5-AGNRs, we also uncover the length-dependent behavior of these end states which transition from singly occupied spin-split states to a closed-shell form as the ribbons become shorter. Finally, we demonstrate the magnetic character of the end states via transport experiments in a model two-terminal device structure in which the ribbons are suspended between the scanning probe and the substrate that both act as leads., We acknowledge funding from the European Union’s Horizon 2020 programme (Grant Agreement Nos. 635919 and 863098 from ERC and FET Open projects, respectively), from the Spanish MINECO (Grant Nos. FIS2017-83780-P and MAT2016-78293-C6), and from the University of the Basque Country (Grant IT1246-19). D.G.O. thanks the Alexander von Humboldt Foundation for supporting his research stay at the MPI, and Klaus Kern for hosting him.

Scattering of electrons by localized spins is the ultimate process enabling detection and control of the magnetic state of a spin-doped material. At the molecular scale, scattering is mediated by the orbitals hosting the spin. Here we report the selective excitation of a molecular spin by tunneling through different molecular orbitals. Spatially resolved tunneling spectra on iron-porphyrins reveal that the inelastic spin excitation extends beyond the iron site, changing shape and symmetry along the molecule. Combining density functional theory simulations with a phenomenological scattering model, we show that the extension and lineshape of the inelastic signal are due to excitation pathways assisted by different frontier orbitals. By selecting the intramolecular site for electron injection, the relative weight of iron and pyrrole orbitals in the tunneling process is modified. Thus, the excitation mechanism, reflected by its spectral lineshape, depends on the degree of localization and energy alignment of the chosen molecular orbital., This work has been funded by the COST 15128 Molecular Spintronics project, by Marie Curie IF ARTE, by the Spanish Ministerio de Economía y Competitividad (MINECO) through the cooperative grant No. MAT2016-78293 and grant No. FIS2016-75862-P, and by the Basque Government (Dep. Industry, Grant PI-2015-1-42, Dep. Education, Grant PI-2016-1-27 and Grant IT-756-13), the EU project PAMS (610446), and the European Regional Development Fund (ERDF).

The theoretical modeling of metal-organic interfaces represents a formidable challenge, especially considering the delicate balance of various interaction mechanisms and the large size of the involved molecular species. In the present study, the energies of interface states, which are known to display a high sensitivity to the adsorption geometry and electronic structure of the deposited molecular species, have been used to test the suitability and reliability of current theoretical approaches. Two well-ordered overlayer structures (relaxed and compressed monolayers) of 1,4,5,8-naphthalene-tetracarboxylic acid dianhydride (NTCDA) on Ag(111) have been investigated using two-photon photoemission to derive precise interface-state energies for these closely related systems. The experimental values are reproduced by our density-functional theory (DFT) calculations with two approaches to treat dispersion interactions (semi-empirical correction DFT-D3 and parametrized functional optB88) and basis set approaches (localized numerical atomic orbitals, plane waves) with remarkable accuracy. Our results underline the trustworthiness and some of the limitations of current DFT-based methods regarding the description of geometric and electronic properties of metal-organic interfaces., This work is a project of the SFB 1083 “Structure and Dynamics of Internal Interfaces” funded by the Deutsche Forschungsgemeinschaft DFG (Germany). D.S.-P. acknowledges support from the Spanish MINECO (Spain), MAT2013-46593-C6-2-P.

arXiv:1301.6489, We use time-dependent density functional theory to calculate the energy loss of an antiproton colliding with a small Al cluster previously excited. The velocity of the antiproton is such that non-linear effects in the electronic response of the Al cluster are relevant. We obtain that an antiproton penetrating an excited cluster transfers less energy to the cluster than an antiproton penetrating a ground state cluster. We quantify this difference and analyze it in terms of the cluster excitation spectrum. © 2013 Elsevier Ltd. All rights reserved., NEK acknowledges support from the CSIC JAE-predoc program, co-financed by the European Science Foundation. We also acknowledge the support of the Basque Departamento de Educación and the UPV/EHU (Grant No. IT-366-07), the Spanish Ministerio de Economía y Competitividad (Grant No. FIS2010-19609-CO2-02) and the ETORTEK program funded by the Basque Departamento de Industria and the Diputación Foral de Gipuzkoa.

We present here a comprehensive search for the structure of the Si(553)-Au reconstruction. More than two hundred different trial structures have been studied using first-principles density-functional calculations with the SIESTA code. An iterative procedure, with a step-by-step increase of the accuracy and computational cost of the calculations, was used to allow for the study of this large number of configurations. We have considered reconstructions restricted to the topmost bilayer and studied two types: i) "flat" surface-bilayer models, where atoms at the topmost bilayer present different coordinations and registries with the underlying bulk, and ii) nine different models based on the substitution of a silicon atom by a gold atom in different positions of a $\pi$-bonded chain reconstruction of the Si(553) surface. We have developed a compact notation that allows us to label and identify all these structures. This is very useful for the automatic generation of trial geometries and counting the number of inequivalent structures, i.e., structures having different bonding topologies. The most stable models are those that exhibit a honeycomb-chain structure at the step edge. One of our models (model f2) reproduces the main features of the room temperature photoemission and scanning-tunneling microscopy data. Thus we conclude that f2 structure is a good candidate for the high temperature structure of the Si(553)-Au surface., Comment: Submitted to Phys. Rev. B on December 10, 2007

Recent photoemission experiments on the Si(553)-Au reconstruction show a one-dimensional band with a peculiar ~1/4 filling. This band could provide an opportunity for observing large spin-charge separation if electron-electron interactions could be increased. To this end, it is necessary to understand in detail the origin of this surface band. A first step is the determination of the structure of the reconstruction. We present here a study of several structural models using first-principles density functional calculations. Our models are based on a plausible analogy with the similar and better known Si(557)-Au surface, and compared against the sole structure proposed to date for the Si(553)-Au system [Crain JN et al., 2004 Phys. Rev. B 69 125401 ]. Results for the energetics and the band structures are given. Lines for the future investigation are also sketched.

Scanning tunneling spectroscopy of p(2x2) Cs and Na ordered overlayers on Cu(111) reveals similar line widths of quasi two-dimensional quantum well states despite largely different binding energies. Detailed calculations show that 50% of the line widths are due to electron-phonon scattering while inelastic electron-electron scattering is negligible. A frequently ignored mechanism for ordered structures, i.e., enhanced elastic scattering due to Brillouin zone back folding, contributes the remaining width., Comment: 4 pages, 2 figures, 1 table

The electronic structure of an infinite poly-guanine - poly-cytosine DNA molecule in its dry A-helix structure is studied by means of density-functional calculations. An extensive study of 30 nucleic base pairs is performed to validate the method. The electronic energy bands of DNA close to the Fermi level are then analyzed in order to clarify the electron transport properties in this particularly simple DNA realization, probably the best suited candidate for conduction. The energy scale found for the relevant band widths, as compared with the energy fluctuations of vibrational or genetic-sequence origin, makes highly implausible the coherent transport of electrons in this system. The possibility of diffusive transport with sub-nanometer mean free paths is, however, still open. Information for model Hamiltonians for conduction is provided., Comment: 8 pages, 4 figures

We have developed and implemented a self-consistent density functional method using standard norm-conserving pseudopotentials and a flexible, numerical LCAO basis set, which includes multiple-zeta and polarization orbitals. Exchange and correlation are treated with the local spin density or generalized gradient approximations. The basis functions and the electron density are projected on a real-space grid, in order to calculate the Hartree and exchange-correlation potentials and matrix elements, with a number of operations that scales linearly with the size of the system. We use a modified energy functional, whose minimization produces orthogonal wavefunctions and the same energy and density as the Kohn-Sham energy functional, without the need of an explicit orthogonalization. Additionally, using localized Wannier-like electron wavefunctions allows the computation time and memory, required to minimize the energy, to also scale linearly with the size of the system. Forces and stresses are also calculated efficiently and accurately, thus allowing structural relaxation and molecular dynamics simulations., Comment: 22 pages, 8 figures

A brief review of the SIESTA project is presented in the context of linear-scaling density-functional methods for electronic-structure calculations and molecular-dynamics simulations of systems with a large number of atoms. Applications of the method to different systems are reviewed, including carbon nanotubes, gold nanostructures, adsorbates on silicon surfaces, and nucleic acids. Also, progress in atomic-orbital bases adapted to linear-scaling methodology is presented., Comment: 5 pages + 4 figures, submitted to Proc. Int. Conf. Solid State Spectroscopy (Sept 1999)

We apply first principles calculations to study the opening of single-wall carbon nanotubes (SWNT's) by oxidation. We show that an oxygen rim can stabilize the edge of the open tube. The sublimation of CO$_2$ molecules from the rim with the subsequent closing of the tube changes from endothermic to exothermic as the tube radius increases, within the range of experimental feasible radii. We also obtain the energies for opening the tube at the cap and at the wall, the latter being significantly less favorable., Comment: 4 pages, 3 figures. Submitted to Phys. Rev. B

A study based on ab initio calculations is presented on the estructural, elastic, and vibrational properties of single-wall carbon nanotubes with different radii and chiralities. We use SIESTA, an implementation of pseudopotential-density-functional theory which allows calculations on systems with a large number of atoms per cell. Different quantities like bond distances, Young moduli, Poisson ratio and the frequencies of different phonon branches are monitored versus tube radius. The validity of expectations based on graphite is explored down to small radii, where some deviations appear related to the curvature effects. For the phonon spectra, the results are compared with the predictions of the simple zone-folding approximation. Except for the known defficiencies of this approximation in the low-frequency vibrational regions, it offers quite accurate results, even for relatively small radii., Comment: 13 pages, 7 figures, submitted to Phys. Rev. B (11 Nov. 98)

The projection of the eigenfunctions obtained in standard plane-wave first-principle calculations is used for analyzing atomic-orbital basis sets. The "spilling" defining the error in such a projection allows the evaluation of the quality of an atomic-orbital basis set for a given system and its systematic variational optimization. The same projection allows obtaining the band structure and the Hamiltonian matrix elements in the previously optimized atomic basis. The spilling is shown to correlate with the mean square error in the energy bands, indicating that the basis optimization via spilling minimization is qualitatively equivalent to the energy-minimization scheme, but involves a much smaller computational effort. The spilling minimization also allows the optimization of bases for uses different than total-energy calculations, like the description of the band gap in semiconductors. The method is applied to the characterization of finite-range pseudo-atomic orbitals {[}O. F. Sankey and D. J. Niklewski, Phys. Rev. B {\bf 40}, 3979 (1989){]} in comparison to infinite-range pseudo-atomic and Slater-type orbitals. The bases are evaluated and optimized for several zincblende semiconductors and for aluminum. The quality of the finite-range orbitals is found to be perfectly comparable to the others with the advantage of a limited range of interacions. A simple scheme is proposed to expand the basis without increasing the range. It is found that a double-$z$ basis substantially improves the basis performance on diamond, whereas $d$-polarization is required for Si and Al for similar results. Finally, the projection allows the chemical analysis of the plane-wave results via population analysis on the previously optimized atomic basis., Comment: 11 pages, RevTeX, plus 13 uuencoded, compressed, tared figures