Search

We study theoretically electron interference in a Mach-Zehnder-like geometry formed by four zigzag graphene nanoribbons (ZGNRs) arranged in parallel pairs, one on top of the other, such that they form intersection angles of 60º. Depending on the interribbon separation, each intersection can be tuned to act either as an electron beam splitter or as a mirror, enabling tuneable circuitry with interfering pathways. Based on the mean-field Hubbard model and Green's function techniques, we evaluate the electron transport properties of such 8-terminal devices and identify pairs of terminals that are subject to self-interference. We further show that the scattering matrix formalism in the approximation of independent scattering at the four individual junctions provides accurate results as compared with the Green's function description, allowing for a simple interpretation of the interference process between two dominant pathways. This enables us to characterize the device sensitivity to phase shifts from an external magnetic flux according to the Aharonov-Bohm effect as well as from small geometric variations in the two path lengths. The proposed devices could find applications as magnetic field sensors and as detectors of phase shifts induced by local scatterers on the different segments, such as adsorbates, impurities or defects. The setup could also be used to create and study quantum entanglement.

Spin-flip excitations of iron porphyrin molecules on Au(111) are investigated with a low-temperature scanning tunneling microscope. The molecules adopt two distinct adsorption configurations on the surface that exhibit different magnetic anisotropy energies. Density functional theory calculations show that the different structures and excitation energies reflect unlike occupations of the Fe 3d levels. We demonstrate that the magnetic anisotropy energy can be controlled by changing the adsorption site, the orientation, or the tip-molecule distance.

The intermolecular interaction among phthalocyanine molecules deposited on Ag(111) has been investigated at submonolayer coverages. By means of infrared absorption spectroscopy and spot-profile analysis low-energy electron diffraction, unambiguous evidence for an attractive interaction between neighboring titanyl-phthalocyanine (TiOPc) molecules is found for dilute layers, which is in contrast to the intermolecular repulsion reported for the copper-phthalocyanine (CuPc)/Ag(111) system. Accordingly, the formation of dense two-dimensional islands of TiOPc molecules with a well-defined long-range order (commensurate phase) is observed upon cooling to temperatures below 75 K, while a disordered arrangement is retained for adsorbed CuPc. Using density functional theory (DFT) calculations as a guide, these differences have been traced to the competition of attractive van der Waals and repulsive electrostatic interactions. Specifically, the vertically oriented molecular dipoles are substantially smaller for adsorbed TiOPc as compared to CuPc due to a compensation of the oppositely oriented dipole moments arising from the Pauli push-back effect on the one hand and the internal molecular dipole associated with the axial Ti–O group on the other. A detailed account and comparison of the molecule–substrate interaction, molecular binding geometries, and intermolecular interactions of TiOPc and CuPc on the Ag(111) surface are provided in a theoretical analysis based on DFT., We gratefully acknowledge support from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation), Project ID 223848855-SFB 1083 “Structure and Dynamics of Internal Interfaces”. L.F. acknowledges financial support from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement MagicFACE no. 797109. D.S.-P. and A.X.B.-R. acknowledge the Spanish Agencia Estatal de Investigación (grant nos. MAT2016-78293-C6-4-R and PID2019-107338RB-C66) and Dep. Educación of the Basque Government and UPV/EHU (grant no. IT1246-19) for support.

The synthesis of novel organic prototypes combining different functionalities is key to achieve operational elements for applications in organic electronics. Here we set the stage towards individually addressable magneto-optical transducers by the on-surface synthesis of optically active manganese-phthalocyanine derivatives (MnPc) obtained directly on a metallic substrate. We created these 2D nanostructures under ultra-high vacuum conditions with atomic precision starting from a simple phthalonitrile precursor with reversible photo-induced reactivity in solution. These precursors maintain their integrity after powder sublimation and coordinate with the Mn ions into tetrameric complexes and then transform into MnPcs on Ag(111) after a cyclotetramerization reaction. Using scanning tunnelling microscopy and spectroscopy together with DFT calculations, we identify the isomeric configuration of two bi-stable structures and show that it is possible to switch them reversibly by mechanical manipulation. Moreover, the robust magnetic moment brought by the central Mn ion provides a feasible pathway towards magneto-optical transducer fabrication. This work should trigger further research confirming such magneto-optical effects in MnPcs both on surfaces and in liquid environments., We gratefully acknowledge financial support from the Spanish Ministries of Economy, Industry and Competitiveness (MINECO, grants MAT2016-78293-C6-4-R and MAT2016-78293-C6-6) and of Science and Innovation (MICINN, grant no. PID2019-107338RB-C62, PID2019-107338RB-C64, PID2019-107338RB-C66/AEI/10.13039/501100011033 and RED2018-102833-T), the regional Governments of Aragon (E13-20R and E12-20R), the Xunta de Galicia (Centro de Investigación de Galicia acreditation 2019–2022, ED431G 2019/03), the Dept. of Education of the Basque Government and UPV/EHU (grant no. IT1246-19), the Spanish Research Agency (AEI), the European Regional Development Fund (ERDF) under the program Interreg V-A España-Francia-Andorra (contract no. EFA194/16 TNSI) and the European Union through Horizon 2020 (FET-Open project SPRING grant. no. 863098)., We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

The excitation of low-energy electron–hole pairs is one of the most relevant processes in the gas–surface interaction. An efficient tool to account for these excitations in simulations of atomic and molecular dynamics at surfaces is the so-called local density friction approximation (LDFA). The LDFA is based on a strong approximation that simplifies the dynamics of the electronic system: a local friction coefficient is defined using the value of the electronic density for the unperturbed system at each point of the dynamics. In this work, we apply real-time time-dependent density functional theory to the problem of the electronic friction of a negative point charge colliding with spherical jellium metal clusters. Our non-adiabatic, parameter-free results provide a benchmark for the widely used LDFA approximation and allow the discussion of various processes relevant to the electronic response of the system in the presence of the projectile., This work has been supported in part by the Basque Departamento de Educación, Universidades e Investigación, the University of the Basque Country UPV/EHU (Grant No. IT1246-19), the Spanish Ministerio de Ciencia e Innovación (MICIN) Grants No. PID2019-107396GB-I00/AEI/10.13039/501100011033 and PID2019-107338RB-C66/AEI/10.13039/501100011033, as well as the Severo Ochoa Program for Centres of Excellence (CEX-2018-000867-S). NK acknowledges the funding from Spanish MICIN through grant PID2019-107338RB-C61/AEI/10.13039/501100011033, as well as a María de Maeztu award to Nanogune, Grant CEX2020-001038-M funded by MICIN/AEI/10.13039/501100011033.

Spin-hosting graphene nanostructures are promising metal-free systems for elementary quantum spintronic devices. Conventionally, spins are protected from quenching by electronic band gaps, which also hinder electronic access to their quantum state. Here, we present a narrow graphene nanoribbon substitutionally doped with boron heteroatoms that combines a metallic character with the presence of localized spin 1/2 states in its interior. The ribbon was fabricated by on-surface synthesis on a Au(111) substrate. Transport measurements through ribbons suspended between the tip and the sample of a scanning tunneling microscope revealed their ballistic behavior, characteristic of metallic nanowires. Conductance spectra show fingerprints of localized spin states in the form of Kondo resonances and inelastic tunneling excitations. Density functional theory rationalizes the metallic character of the graphene nanoribbon due to the partial depopulation of the valence band induced by the boron atoms. The transferred charge builds localized magnetic moments around the boron atoms. The orthogonal symmetry of the spin-hosting state’s and the valence band’s wave functions protects them from mixing, maintaining the spin states localized. The combination of ballistic transport and spin localization into a single graphene nanoribbon offers the perspective of electronically addressing and controlling carbon spins in real device architectures., We gratefully acknowledge financial support from Grants PID2019-107338RB-C61, PID2019-107338RB-C62, PID2019-107338RB-C66, PID2019-110037GB-I00, and PCI2019-111933-2 and the Maria de Maeztu Units of Excellence Program CEX2020-001038-M funded by MCIN/AEI/10.13039/501100011033, the European Regional Development Fund, the European Union (EU) H2020 program through the FET Open project SPRING (Grant Agreement No. 863098), the Xunta de Galicia (Centro de Investigación de Galicia accreditation 2019–2022, ED431G 2019/03), the Dpto. Educación Gobierno Vasco (Grant Nos. PIBA-2020-1-0014, IT1246-19, and IT-1569-22) and the Programa Red Guipuzcoana de Ciencia, Tecnología e Innovación 2021 (Grant No. 2021-CIEN-000070-01. Gipuzkoa Next). The authors acknowledge the financial support received from the IKUR Strategy under the collaboration agreement between Ikerbasque Foundation and DIPC on behalf of the Department of Education of the Basque Government.

A prerequisite to characterize magnetic materials is the capability to describe systems containing unpaired electrons. In this study, we benchmark the one-shot GW (G0W0) on top of different unrestricted mean-field solutions for open-shell molecules using Dunning's correlation-consistent basis sets expanded in terms of Gaussian functions. We find that the G0W0 correction to hybrid functionals provides reasonably accurate results for the ionization energies of open-shell systems when compared to those obtained from high-level ab initio methods. Moreover, the quality of the G0W0 exchange–correlation approximation is evaluated by the discrepancy between the ionization energy of the neutral molecules and the electron affinity of the corresponding cations. Furthermore, we assess the capability of the GW to reproduce the correct energy ordering of molecular spin–orbitals. To such an aim, we thoroughly discuss three open-shell molecules CN, NH2, and O2, for which approximate functionals fail to correctly capture the single-electron spectrum. Particularly, we demonstrate that the overestimation of the exchange energy in the studied spin–orbitals is reduced by the GW dynamic correlation term, restoring the molecular orbital ordering. Interestingly, we find that deviations of the exchange and correlation energies, in comparison with our ab initio reference, can be very different for molecular orbitals with different symmetry, e.g. σ and π-type orbitals., Authors acknowledge support from Spanish Agencia Estatal de Investigación (Grant Nos. PID2019-107338RB-C66, PID2019-109555GB-I00 and RTC-2016-5681-7), and the Eusko Jaurlaritza and UPV/EHU (Grant Nos. PIBA19-0004 and IT1246-19). DSP also acknowledges support from the European Union (EU) through Horizon 2020 (FET-Open project SPRING Grant No. 863098).

By means of spin-resolved density functional theory calculations using both atomic orbitals and plane-wave basis codes, we study the electronic and magnetic ground state of single-layer NbSe

By means of spin-resolved density functional theory calculations using both atomic orbitals and plane-wave basis codes, we study the electronic and magnetic ground state of single-layer NbSe2. We find that, for all the functionals considered, the most stable solution in this two-dimensional (2D) superconductor is the ferrimagnetic ground state with a magnetic moment of 1.09 μB at the Nb atoms and of 0.05 μB at the Se atoms pointing in the opposite direction. Our calculations show that the ferrimagnetic state precludes the development of charge density wave (CDW) order and their coexistence in the single-layer limit, unless graphene is considered as a substrate. The spin-resolved calculated density of states (DOS), a key fingerprint of the electronic and magnetic structure of a material, unambiguously reproduces the experimental DOS measured by scanning tunneling spectroscopy in single-layer NbSe2. Our work sets magnetism into play in this prototypical correlated 2D material, which is crucial to understand the formation mechanisms of 2D superconductivity and CDW order., We are indebted to Spanish Ministry of Science and Innovation for financial support through Grant FIS2015-64886-C5-5-P and MAT2017-83553-P. MMU acknowledges support by the Spanish MINECO under Grant No. MAT2017-88377-C2-1-R and by the ERC Starting Grant LINKSPM (Grant 758558). EB acknowledges support from the Scientific and Technological Research Council of Turkey (TUBİTAK) under Grant No. 1059B141801387 and DSP from Spanish Ministry of Science and Innovation under Grant PID2019-107338RB-C66.

This article is part of the JCP Special Topic on Electronic Structure Software., A review of the present status, recent enhancements, and applicability of the SIESTA program is presented. Since its debut in the mid-1990s, SIESTA’s flexibility, efficiency, and free distribution have given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of SIESTA combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a realspace grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here, we describe the more recent implementations on top of that core scheme, which include full spin–orbit interaction, non-repeated and multiple-contact ballistic electron transport, density functional theory (DFT)+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as WANNIER90 and the second-principles modeling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering SIESTA runs, and various post-processing utilities. SIESTA has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low-level libraries, as well as data standards and support for them, particularly the PSeudopotential Markup Language definition and library for transferable pseudopotentials, and the interface to the ELectronic Structure Infrastructure library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments., Siesta development was historically supported by different Spanish National Plan projects (Project Nos. MEC-DGES-PB95-0202, MCyT-BFM2000-1312, MEC-BFM2003-03372, FIS2006-12117, FIS2009-12721, FIS2012-37549, FIS2015-64886-P, and RTC-2016-5681-7), the latter one together with Simune Atomistics Ltd. We are thankful for financial support from the Spanish Ministry of Science, Innovation and Universities through Grant No. PGC2018-096955-B. We acknowledge the Severo Ochoa Center of Excellence Program [Grant Nos. SEV-2015-0496 (ICMAB) and SEV-2017-0706 (ICN2)], the GenCat (Grant No. 2017SGR1506), and the European Union MaX Center of Excellence (EU-H2020 Grant No. 824143). P.G.-F. acknowledges support from Ramón y Cajal (Grant No. RyC-2013-12515). J.I.C. acknowledges Grant No. RTI2018-097895-B-C41. R.C. acknowledges the European Union’s Horizon 2020 Research and Innovation Program under Marie Skłodoswka-Curie Grant Agreement No. 665919. D.S.P, P.K., and P.B. acknowledge Grant No. MAT2016-78293-C6, FET-Open No. 863098, and UPV-EHU Grant No. IT1246-19. V. W. Yu was supported by a MolSSI Fellowship (U.S. NSF Award No. 1547580), and V.B. and V.W.Y. were supported by the ELSI Development by the NSF (Award No. 1450280).

Combining on-surface synthetic methods with the power of scanning tunneling microscopy to characterize novel materials at the single molecule level, we show how to steer the reactivity of one anthracene-based precursor towards different product nanostructures. Whereas using a Au(111) surface with three-fold symmetry results in the dominant formation of a starphene derivative, the two-fold symmetry of a reconstructed Au(110) surface allows the selective growth of non-benzenoid linear conjugated polymers. We further assess the electronic properties of each of the observed product structures via tunneling spectroscopy and DFT calculations, altogether advancing the synthesis and characterization of molecular structures of notable scientific interest that have been only scarcely investigated to date, as applies both to starphenes and to non-benzenoid conjugated polymers., We acknowledge funding from the European Union's Horizon 2020 programme (Grant Agreement No. 635919 and 863098), from the Spanish Agencia Estatal de Investigación (Grant No. MAT2016-78293-C6, PID2019-107338RB-C63, PID2019-107338RB-C66 and FIS2017-83780-P), Xunta de Galicia (Centro singular de investigación de Galicia, accreditation 2019–2022, ED431G 2009/03), Dep. Educación of the Basque Government and UPV/EHU (Grant No. IT-756-13), and the European Regional Development Fund (ERDF).

The catalytic oxidation of CO on transition metals, such as Pt, is commonly viewed as a sharp transition from the CO-inhibited surface to the active metal, covered with O. However, we find that minor amounts of O are present in the CO-poisoned layer that explain why, surprisingly, CO desorbs at stepped and flat Pt crystal planes at once, regardless of the reaction conditions. Using near-ambient pressure X-ray photoemission and a curved Pt(111) crystal we probe the chemical composition at surfaces with variable step density during the CO oxidation reaction. Analysis of C and O core levels across the curved crystal reveals that, right before light-off, subsurface O builds up within (111) terraces. This is key to trigger the simultaneous ignition of the catalytic reaction at different Pt surfaces: a CO-Pt-O complex is formed that equals the CO chemisorption energy at terraces and steps, leading to the abrupt desorption of poisoning CO from all crystal facets at the same temperature., We acknowledge financial support from the Spanish Ministry of Science and Innovation (Grants MAT-2017-88374-P, PID2019-107338RB-C63, PID2019-107338RB-C66), the Basque Government (Grants IT-1255-19, IT-1246-19), the Knut and Alice Wallenberg foundation (DNR KAW 2015.0058 „Atomistic Design of new Catalysts“) and the Swedish Research Council (DNR 349-2011-6491 „Catalysis on the atomic scale“).

For the layered transition metal dichalcogenide 1T-TaS2, we establish through a unique experimental approach and density functional theory, how ultrafast charge transfer in 1T-TaS2 takes on isotropic three-dimensional character or anisotropic two-dimensional character, depending on the commensurability of the charge density wave phases of 1T-TaS2. The X-ray spectroscopic core-hole-clock method prepares selectively in- and out-of-plane polarized sulfur 3p orbital occupation with respect to the 1T-TaS2 planes and monitors sub-femtosecond wave packet delocalization. Despite being a prototypical two-dimensional material, isotropic three-dimensional charge transfer is found in the commensurate charge density wave phase (CCDW), indicating strong coupling between layers. In contrast, anisotropic two-dimensional charge transfer occurs for the nearly commensurate phase (NCDW). In direct comparison, theory shows that interlayer interaction in the CCDW phase - not layer stacking variations - causes isotropic three-dimensional charge transfer. This is presumably a general mechanism for phase transitions and tailored properties of dichalcogenides with charge density waves., AF acknowledges funding by the Helmholtz Virtual Institute VI 419 and the ERC Advanced Investigator Grant No. 669531 EDAX. NM acknowledges funding by the ERC Advanced Investigator Grant No. 321319 ESUX and support from the Carl Tryggers foundation for Scientifc Research CTS. DS-P and MM acknowledge funding from EU-FP7 project THINFACE (grant No. 607232), Spanish MINECO (grant No. MAT2016-78293-C6-4-R), and Basque Dep.de Educación and UPV/EHU (grant No. IT-756-13).

A nanographene formed by the fusion of 22 benzene rings has been prepared by combining an in-solution Pd-catalyzed cycloaddition reaction and on-surface Au-promoted cyclodehydrogenation. The structure and electronic properties of the resulting three-fold symmetric CH molecule have been characterized by scanning probe microscopy with atomic resolution and corroborated by theoretical modelling., This research was supported by the National Science Center, Poland (2017/26/E/ST3/00855), the EU project PAMS (610446), the Agencia Estatal de Investigacion (MAT2016-78293-C6-3-R, MAT2016-78293-C6-4-R, FIS2017-83780-P and CTQ2016-78157-R), the Xunta de Galicia (Centro singular de investigacion de Galicia, acreditacion 2016-2019, ED431G/09) and the European Regional Development Fund (ERDF). R. Z. acknowledges support received from the National Science Center, Poland (2017/24/T/ST5/00262).

We report the energy level alignment evolution of valence and conduction bands of armchair-oriented graphene nanoribbons (aGNR) as their band gap shrinks with increasing width. We use 4,4″-dibromo-para-terphenyl as the molecular precursor on Au(111) to form extended poly-para-phenylene nanowires, which can subsequently be fused sideways to form atomically precise aGNRs of varying widths. We measure the frontier bands by means of scanning tunneling spectroscopy, corroborating that the nanoribbon's band gap is inversely proportional to their width. Interestingly, valence bands are found to show Fermi level pinning as the band gap decreases below a threshold value around 1.7 eV. Such behavior is of critical importance to understand the properties of potential contacts in GNR-based devices. Our measurements further reveal a particularly interesting system for studying Fermi level pinning by modifying an adsorbate's band gap while maintaining an almost unchanged interface chemistry defined by substrate and adsorbate., The project leading to this publication has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 635919), from the European Union FP7 FET-ICT “Planar Atomic and Molecular Scale devices” (PAMS) project (Contract No. 610446), from the Spanish Ministry of Economy, Industry and Competitiveness (MINECO, Grant No. MAT2016-78293-C6), from the Basque Department of Education (PI-2016-1-0027), and from the University of Padova (Grant CPDA154322, Project AMNES).

We report the on-surface synthesis of 7-armchair graphene nanoribbons (7-AGNRs) substituted with nitrile (CN) functional groups. The CN groups are attached to the GNR backbone by modifying the 7-AGNR precursor. Whereas many of these groups survive the on-surface synthesis, the reaction process causes the cleavage of some CN from the ribbon backbone and the on-surface cycloisomerization of few nitriles onto pyridine rings. Scanning tunneling spectroscopy and density functional theory reveal that CN groups behave as very efficient n-dopants, significantly downshifting the bands of the ribbon and introducing deep impurity levels associated with the nitrogen electron lone pairs., This work was supported by FP7 FET-ICT “Planar Atomic and Molecular Scale devices” (PAMS) project (funded by the European Commission under Contract No. 610446), by the Agencia Estatal de Investigacion (Cooperative Grant No. MAT2016-78293 and Grant FIS2015-62538-ERC), the Basque Government (Dep. de Educacion and UPV/EHU, Grant No. IT-756-13, and Dep. Industry, Grant PI_2015_1_42), the Xunta de Galicia (Centro singular de investigacion de Galicia accreditation 2016−2019, ED431G/09), and the European Regional Development Fund (ERDF).

We analyze theoretically four-terminal electronic devices composed of two crossed graphene nanoribbons (GNRs) and show that they can function as beam splitters or mirrors. These features are identified for electrons in the low-energy region where a single valence or conduction band is present. Our modeling is based on p z orbital tight binding with Slater-Koster-type matrix elements fitted to accurately reproduce the low-energy bands from density functional theory calculations. We analyze systematically all devices that can be constructed with either zigzag or armchair GNRs in AA and AB stackings. From Green's function theory the elastic electron transport properties are quantified as a function of the ribbon width. We find that devices composed of relatively narrow zigzag GNRs and AA-stacked armchair GNRs are the most interesting candidates to realize electron beam splitters with a close to 50:50 ratio in the two outgoing terminals. Structures with wider ribbons instead provide electron mirrors, where the electron wave is mostly transferred into the outgoing terminal of the other ribbon, or electron filters where the scattering depends sensitively on the wavelength of the propagating electron. We also test the robustness of these transport properties against variations in the intersection angle, stacking pattern, lattice deformation (uniaxial strain), inter-GNR separation, and electrostatic potential differences between the layers. These generic features show that GNRs are interesting basic components to construct electronic quantum optical setups., This work was supported by the project Spanish Ministerio de Economía y Competitividad (MINECO) through the Grants no. FIS2017-83780-P (Graphene Nanostructures “GRANAS”) and no.MAT2016-78293-C6-4R, the Basque Departamento de Educación through the PhD fellowship no. PRE_2019_2_0218 (S.S.), the University of the Basque Country through the Grant no. IT1246-19, and the European Union (EU) through Horizon 2020 (FET-Open project SPRING Grant no. 863098).

Graphene nanoribbons (GNRs), low-dimensional platforms for carbon-based electronics, show the promising perspective to also incorporate spin polarization in their conjugated electron system. However, magnetism in GNRs is generally associated with localized states around zigzag edges, difficult to fabricate and with high reactivity. Here we demonstrate that magnetism can also be induced away from physical GNR zigzag edges through atomically precise engineering topological defects in its interior. A pair of substitutional boron atoms inserted in the carbon backbone breaks the conjugation of their topological bands and builds two spin-polarized boundary states around them. The spin state was detected in electrical transport measurements through boron-substituted GNRs suspended between the tip and the sample of a scanning tunneling microscope. First-principle simulations find that boron pairs induce a spin 1, which is modified by tuning the spacing between pairs. Our results demonstrate a route to embed spin chains in GNRs, turning them into basic elements of spintronic devices., We gratefully acknowledge financial support from Spanish Agencia Estatal de Investigación (AEI) (MAT2016-78293, PID2019-107338RB, FIS2017-83780-P, and the Maria de Maeztu Units of Excellence Programme MDM-2016-0618), from the European Union (EU) through Horizon 2020 (FET-Open project SPRING Grant. No. 863098), the Basque Departamento de Educación through the PhD fellowship No. PRE_2019_2_0218 (S.S.), the Xunta de Galicia (Centro de Investigación de Galicia accreditation 2019–2022, ED431G 2019/03), the University of the Basque Country (Grant IT1246-19), and the European Regional Development Fund (ERDF). I. P. also thanks Xunta de Galicia and European Union (European Social Fund, ESF) for the award of a predoctoral fellowship

The injection of spin-polarized electrons across interfaces is central to many technologies, and hence, it is important to understand the main ingredients controlling it. Here, we demonstrate that the spin dependence of ultrafast electron transfer at Ar/Co(0001) and Ar/Fe(110) interfaces is rooted in the details of the spin-split surface band structures. The injection dynamics are particularly sensitive to the sizes (in reciprocal space) of projected electronic band gaps around Γ̅. Our ab initio calculations back that minority electrons are injected significantly faster than majority electrons in line with recently reported experimental injection times. A simple tunnelling model incorporating the spin-dependent gap sizes confirms that this ingredient is crucial to rationalize the experimental results., We acknowledge illuminating discussions with Prof. U. Heinzmann and Prof. P. Feulner, as well as financial support from the EU FP7 programme under Grant Agreement No. 607232 (THINFACE), the Spanish Ministerio de Ciencia, Innovacion y Universidades (Grant. No. MAT2016-78293-C6-4-R and Grant No. PID2019-107338RB-C66) and the Basque Dep. de Educacion and the UPV/EHU (Grant No. IT1246-19).

The electronic properties of graphene nanoribbons (GNRs) can be precisely tuned by chemical doping. Here we demonstrate that amino (NH2) functional groups attached at the edges of chiral GNRs (chGNRs) can efficiently gate the chGNRs and lead to the valence band (VB) depopulation on a metallic surface. The NH2-doped chGNRs are grown by on-surface synthesis on Au(111) using functionalized bianthracene precursors. Scanning tunneling spectroscopy resolves that the NH2 groups significantly upshift the bands of chGNRs, causing the Fermi level crossing of the VB onset of chGNRs. Through density functional theory simulations we confirm that the hole-doping behavior is due to an upward shift of the bands induced by the edge NH2 groups., We acknowledge financial support from (i) AEI/FEDER-EU through grant nos. MAT2016-78293-C6, FIS2017-83780-P AEI/FEDER, EU, the Maria de Maeztu unit of excellence MDM-2016-0618, and the Severo Ochoa program (ICN2) SEV-2017-0706; (ii) the European Research Council (grant agreement no. 635919); (iii) the Xunta de Galicia (Centro Singular De Investigacion de Galicia, accreditation 2016-2019, ED431G/09); (iv) the European FET-OPEN project SPRING (#863098); (v) the European Regional Development Fund (ERDF) under the program Interreg V-A España-Francia-Andorra (contract no. EFA 194/16 TNI); (vi) the CERCA Program/Generalitat de Catalunya; (vii) the Gobierno Vasco UPV/EHU (project IT1246-19); (viii) Basque Government Project (IT-1255-19); and (ix) Diputacion Foral De Gipuzkoa RED 2019-096.

Graphene nanoribbons (GNRs) and their derivatives attract growing attention due to their excellent electronic and magnetic properties as well as the fine-tuning of such properties that can be obtained by heteroatom substitution and/or edge morphology modification. Here, we introduce graphene nanoribbon derivatives—organometallic hybrids with gold atoms incorporated between the carbon skeleton and side Cl atoms. We show that narrow chlorinated 5-AGNROHs (armchair graphene nanoribbon organometallic hybrids) can be fabricated by on-surface polymerization with omission of the cyclodehydrogenation reaction by a proper choice of tailored molecular precursors. Finally, we describe a route to exchange chlorine atoms connected through gold atoms to the carbon skeleton by hydrogen atom treatment. This is achieved directly on the surface, resulting in perfect unsubstituted hydrogen-terminated GNRs. This will be beneficial in the molecule on-surface processing when the preparation of final unsubstituted hydrocarbon structure is desired., This work was supported by the National Science Center, Poland (2019/35/B/ST5/02666), and the EU Project PAMS (610446). A.G.L., D.S.P., and P.B. acknowledge the Spanish Agencia Estatal de Investigación (Grants MAT2016-78293-C6-4-R, PID2019-107338RB-C66, and FIS2017-83780-P), Dep. Educación of the Basque Government and UPV/EHU (Grant IT-756-13), and the European Union (EU) through Horizon 2020 (FET-Open project SPRING Grant 863098) for support.

Miniaturization of electronic circuits into the single-atom level requires novel approaches to characterize transport properties. Due to its unrivaled precision, scanning probe microscopy is regarded as the method of choice for local characterization of atoms and single molecules supported on surfaces. Here we investigate electronic transport along the anisotropic germanium (001) surface with the use of two-probe scanning tunneling spectroscopy and first-principles transport calculations. We introduce a method for the determination of the transconductance in our two-probe experimental setup and demonstrate how it captures energy-resolved information about electronic transport through the unoccupied surface states. The sequential opening of two transport channels within the quasi-one-dimensional Ge dimer rows in the surface gives rise to two distinct resonances in the transconductance spectroscopic signal, consistent with phase-coherence lengths of up to 50 nm and anisotropic electron propagation. Our work paves the way for the electronic transport characterization of quantum circuits engineered on surfaces., This work was supported by the FP7 FET-ICT “Planar Atomic and Molecular Scale devices” (PAMS) project (funded by the European Commission under contract number 610446), the Polish Ministry for Science and Higher Education from financial resources for science in 2013–2017 granted for an international co-financed project (contract number 2913/7.PR/2013/2) and the Spanish Ministerio de Economía y Competitividad (MINECO) (Grant Numbers MAT2016-78293-C6-4-R and FIS2017-83780-P). M.K. acknowledges financial support received from the Polish Ministry of Science and Higher Education, contract number 0341/IP3/2016/74.

San Sebastian (Donostia in the Basque language) is a beautiful, medium-sized European city of about 200,000 inhabitants located at the heart of the Basque Country, next to the Atlantic Ocean and the Pyrenees. Long known as a tourist and culinary destination, San Sebastian has recently become a city of science and technology, with a strong focus on nanoscience research and characterized by an interdisciplinary approach and the support of intense industrial activity in the region. The San Sebastian nanoscience community has contributed considerably to ACS Nano with more than 100 publications during the past decade, some of which have already had significant impact and are highlighted in this virtual issue. This activity has turned this small point at the Basque coastline into a worldwide “nano”-pole for a variety of focused areas of nanoscience research.

The use of atomic orbitals in Hedin's

The evolution of titanyl-phthalocyanine (TiOPc) thin films on Ag(111) has been investigated using IRAS (infrared absorption spectroscopy), SPA-LEED (low energy electron diffraction) and STM (scanning tunneling microscopy). In the (sub)monolayer regime various phases are observed that can be assigned to a 2D gas, a commensurate and a point-on-line phase. In all three phases the nonplanar TiOPc molecule is adsorbed on Ag(111) in an oxygen-up configuration with the molecular π-conjugated backbone oriented parallel to the surface. The commensurate phase reveals a high packing density, containing two molecules at inequivalent adsorption sites within the unit cell. Both molecules assume different azimuthal orientations with respect to the Ag(111) substrate which is ascribed to the interplay between the occupation of preferred sites as well as distinct azimuthal orientations and a short ranging Pauli-Type repulsion to avoid overlapping of neighboring molecules. At full saturation of the monolayer the latter interaction leads to a further azimuthal reorientation that enables an even higher packing density. Thereby, the plain commensurate order is lost and a point-on-line phase is formed. DFT calculations have been used to study different adsorption geometries of TiOPc on Ag(111). The most stable configurations among those with pointing up oxygen atoms (bridge+, bridgex, topx) seem to correspond to those identified experimentally. The calculated dependence of the electronic structure and molecular dipole on the adsorption site and configuration is found to be rather small., We gratefully acknowledge support from the Deutsche Forschungsgemeinschaft DFG (Germany) through the collaborative research center “Structure and Dynamics of Internal Interfaces’’ (SFB 1083, Projects A2, A3, and GP1). DSP acknowledges support from the Spanish MINECO (Grant No. MAT2013-46593-C6-2-P), the Basque Departamento de Educacion, and the UPV/EHU (Grant No. IT-756-13).

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.

The use of atomic orbitals in Hedin's GW approximation provides, in principle, an inexpensive alternative to plane-wave basis sets, especially when modeling large molecules. However, benchmarking of the algorithms and basis sets is essential for a careful balance between cost and accuracy. In this paper, we present an implementation of the GW approximation using numerical atomic orbitals and a pseudopotential treatment of core electrons. The combination of a contour deformation technique with a one-shot extraction of quasiparticle energies provides an efficient scheme for many applications. The performance of the implementation with respect to the basis set convergence and the effect of the use of pseudopotentials has been tested for the 117 closed-shell molecules from the G2/97 test set and 24 larger acceptor molecules from another recently proposed test set. Moreover, to demonstrate the potential of our method, we compute the thermally averaged GW density of states of a large photochromic compound by sampling ab initio molecular dynamics trajectories at different temperatures. The computed thermal line widths indicate approximately twice as large electron-phonon couplings with GW than with standard DFT-GGA calculations. This is further confirmed using frozen-phonon calculations.

For the layered transition metal dichalcogenide 1T-TaS2, we establish through a unique experimental approach and density functional theory, how ultrafast charge transfer in 1T-TaS2 takes on isotropic three-dimensional character or anisotropic two-dimensional character, depending on the commensurability of the charge density wave phases of 1T-TaS2. The X-ray spectroscopic core-hole-clock method prepares selectively in- and out-of-plane polarized sulfur 3p orbital occupation with respect to the 1T-TaS2 planes and monitors sub-femtosecond wave packet delocalization. Despite being a prototypical two-dimensional material, isotropic three-dimensional charge transfer is found in the commensurate charge density wave phase (CCDW), indicating strong coupling between layers. In contrast, anisotropic two-dimensional charge transfer occurs for the nearly commensurate phase (NCDW). In direct comparison, theory shows that interlayer interaction in the CCDW phase – not layer stacking variations – causes isotropic three-dimensional charge transfer. This is presumably a general mechanism for phase transitions and tailored properties of dichalcogenides with charge density waves.

Understanding the mechanisms involved in the covalent attachment of organic molecules to surfaces is a major challenge for nanotechnology and surface science. On the basis of classical organic chemistry mechanistic considerations, key issues such as selectivity and reactivity of the organic adsorbates could be rationalized and exploited for the design of molecular-scale circuits and devices. Here we use tris(benzocyclobutadieno)triphenylene, a singular Y-shaped hydrocarbon containing antiaromatic cyclobutadienoid rings, as a molecular probe to study the reaction of polycyclic conjugated molecules with atomic scale moieties, dangling-bond (DB) dimers on a hydrogen-passivated Ge(001):H surface. By combining molecular design, synthesis, scanning tunneling microscopy and spectroscopy (STM/STS) and computational modeling, we show that the attachment involves a concerted [4+2] cycloaddition reaction that is completely site-selective and fully reversible. This selectivity, governed by the bond alternation induced by the presence of the cyclobutadienoid rings, allows for the control of the orientation of the molecules with respect to the surface DB-patterning. We also demonstrate that by judicious modification of the electronic levels of the polycyclic benzenoid through substituents, the reaction barrier height can be modified. Finally, we show that after deliberate tip-induced covalent bond cleavage, adsorbed molecules can be used to fine tune the electronic states of the DB dimer. This power to engineer deliberately the bonding configuration and electronic properties opens new perspectives for creating prototypical nanoscale circuitry., This work was supported by the FP7 FET-ICT “Planar Atomic and Molecular Scale devices” (PAMS) project (funded by the European Commission under contract no. 610446), by the Polish Ministry for Science and Higher Education from financial resources for science in 2013–2017 granted for an international co-financed project (contract no. 2913/7.PR/2013/2) and by projects CTQ2016-78157-R, MAT2016-78293-C6-3R and MAT2016-78293-C6-4R (AEI/FEDER, UE). EG, DP and DP are grateful for financial support from the Xunta de Galicia (Centro singular de investigación de Galicia accreditation 2016–2019, ED431G/09) and the European Union (European Regional Development Fund-ERDF).

This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence., Dangling bond (DB) arrays on Si(001):H and Ge(001):H surfaces can be patterned with atomic precision and they exhibit complex and rich physics making them interesting from both technological and fundamental perspectives. But their complex behavior often makes scanning tunneling microscopy (STM) images difficult to interpret and simulate. Recently it was shown that low-temperature imaging of unoccupied states of an unpassivated dimer on Ge(001):H results in a symmetric butterfly-like STM pattern, despite that the equilibrium dimer configuration is expected to be a bistable, buckled geometry. Here, based on a thorough characterization of the low-bias switching events, we propose a new imaging model featuring a dynamical two-state rate equation. This model allows us to reproduce the features of the observed symmetric empty-state images which strongly corroborates the idea that the patterns arise due to fast switching events and provides insight into the relation between the tunneling current and switching rates. Our new imaging model is general and can be applied to other systems that exhibit rapid fluctuations during STM experiments., The work is funded by the European Commission under the FP7 FET-ICT “Planar Atomic and Molecular Scale devices” (PAMS) project (Contract No. 610446) and National Science Centre, Poland (2014/15/D/ST3/02975). ME, TF, and DSP also acknowledge support from the Spanish Ministerio de Economía y Competitividad (MINECO) (Grant No. MAT2013-46593-C6-2-P) and the Basque Dep. de Educación and the UPV/EHU (Grant No. IT-756-13). MK acknowledges financial support received from the Foundation for Polish Science (FNP). RZ acknowledges support received from KNOW (scholarship KNOW/44/SS/RZ/2015)., We acknowledge support by the CSIC Open Access Publication Initiative through its Unit of Information Resources for Research (URICI).

Hydrogen molecules can be trapped in the nanocavity formed by the tip of a scanning tunneling microscope and a metal or molecular surface, and produce sharp inelastic non-linearities in the tunneling spectra. Here, we study the interaction effects of hydrogen in a tunneling junction created over Manganese phthalocyanines molecules in two oxidation states. The effect of hydrogen in the tunneling spectra varies strongly depending on the molecular species, but its force spectrum is fairly independent on the molecular state. We find that in mild tunneling conditions hydrogen-induced forces are weakly attractive during a small range of tip sample distance. The van der Waals interaction shows a maximum value of 140 pN, which faintly depends on electrostatic variations along the surface. These results show that AFM can be employed to resolve the complex interaction landscape of a trapped hydrogen molecule and deduce fainter effects such as molecular deformations or dipolar fields., The research was financially supported by the Spanish Ministerio de Economia y Competitividad, MINECO (Grant FIS 2015-62538-ERC, MAT2016-78293-C6, and the Maria de Maeztu Units of Excellence Programme MDM-2016-0618), the Basque Government (Departamento de Industria, grant no. PI-2015-1-42, and Departamento de Educación and UPV/EHU, grant No. IT-756-13), and the European Regional Development Fund (ERDF).

Plasmonic gaps are known to produce nanoscale localization and enhancement of optical fields, providing small effective mode volumes of about a few hundred nm. Atomistic quantum calculations based on time-dependent density functional theory reveal the effect of subnanometric localization of electromagnetic fields due to the presence of atomic-scale features at the interfaces of plasmonic gaps. Using a classical model, we explain this as a nonresonant lightning rod effect at the atomic scale that produces an extra enhancement over that of the plasmonic background. The near-field distribution of atomic-scale hot spots around atomic features is robust against dynamical screening and spill-out effects and follows the potential landscape determined by the electron density around the atomic sites. A detailed comparison of the field distribution around atomic hot spots from full quantum atomistic calculations and from the local classical approach considering the geometrical profile of the atoms' electronic density validates the use of a classical framework to determine the effective mode volume in these extreme subnanometric optical cavities. This finding is of practical importance for the community of surface-enhanced molecular spectroscopy and quantum nanophotonics, as it provides an adequate description of the local electromagnetic fields around atomic-scale features with use of simplified classical methods., Financial support from Project No. FIS2016-80174-P and MAT2016-78293-C6-4-R of MINECO, grant of Consolidated Groups at UPV/EHU (IT-756-13) of the Basque Government, and project 70NANB15H32 from U.S. Department of Commerce, National Institute of Standards and Technology, is acknowledged. M.U. acknowledges support from the University of the Basque Country through a Ph.D. grant, as well as DIPC and CFM at the initial stages of this work. M.B. acknowledges support from the Departamento de Educacion of the Basque Government through a Ph.D. grant, as well as from Euskampus and the DIPC at the initial stages of this work. P.K. acknowledges financial support from the Fellows Gipuzkoa program of the Gipuzkoako Foru Aldundia through the FEDER funding scheme of the European Union, “Una manera de hacer Europa”.

We present an algorithm and its implementation to calculate the properties of electronic excitations in molecules and clusters from first principles, using time-dependent density functional theory (TDDFT). The algorithm assumes the use of some localized functions as a basis set to represent the spatial degrees of freedom. It relies on an iterative computation of the induced density according to the Dyson-like equation for the linear response function. The current implementation is built upon so-called numerical atomic orbitals. It is suitable for a wide variety of density functional theory (DFT) software. In this work, we demonstrate TDDFT calculations starting from preceding DFT runs with SIESTA, GPAW and PySCF packages, while a coupling with the other DFT packages such as Fireball and OpenMX is planned. The mentioned packages are capable of performing ab initio molecular dynamics simulations, and the speed of our TDDFT implementation makes feasible to perform a configuration average of the optical absorption spectra. Our code is written mostly in Python language allowing for a quick and compact implementation of most numerical methods and data-managing tasks with the help of NumPy/SciPy libraries and Python intrinsic constructs. Part of the code is written in C and Fortran to achieve a competitive speed in particular sections of the algorithm. Many parts of the current algorithm and implementation are useful in other ab initio methods for electronic excited states, such as Hedin¿s GW, Bethe¿Salpeter equation and DFT with hybrid functionals. Corresponding proof-of-principles implementations are already part of the code., PK and DSP acknowledge support from Spanish MINECO Grants MAT2016-78293-C6-4-R and RTC-2016-5681-7 (SIESTA-PRO). PKacknowledges financial support from the Fellows Gipuzkoa program of the Gipuzkoako Foru Aldundia through the FEDER fundings cheme of the European Union.MB acknowledges support from the Departemento de Educación of the Basque Government through aPhD grant, as well as from Euskampus and the DIPC at the initial stages of this work

This article is part of themed collection: 2015 Hot Papers in Nanoscale.-- arXiv:1412.7698v1.-- et al., We present a quantitative exploration, combining experiment and simulation, of the mechanical and electronic properties, as well as the modifications induced by an alkylthiolated coating, at the single nanoparticle (NP) level. We determined the response of the NPs to external pressure in a controlled manner using an atomic force microscope tip. We found a strong reduction in their Young's modulus, as compared to bulk gold, and a significant influence of strain on the electronic properties of the alkylthiolated NPs. Electron transport measurements of tiny molecular junctions (NP/alkylthiol/CAFM tip) show that the effective tunnelling barrier through the adsorbed monolayer strongly decreases by increasing the applied load, which translates in a remarkable and unprecedented increase in the tunnel current. These observations are successfully explained using simulations based on the finite element analysis (FEA) and first-principles calculations that permit one to consider the coupling between the mechanical response of the system and the electric dipole variations at the interface., P.L. is a Senior Research Associate from the Fund for Scientific Research of Belgium (F.R.S. – FNRS). K.S. has been supported by the NordPas-de Calais Council fund and the ANR project SAGE III-V (no. ANR11BS1001203) and S.D. by EU project I-ONE (FP7 no. 280772). The experiments were partly funded by the SINGLEMOL project supported by the Nord-Pas-de Calais council fund. We also acknowledge funding from the Basque Government (Grant No. IT-756-13) and the Spanish MINECO (Grant Nos. MAT2013-46593-C6-2-P and FIS2013-48286-C2-2-P).

This work was partially supported by the Basque Departamento de Educacion, Universidades e Investigacion, the University of the Basque Country UPV/EHU (Grant No. IT-756-13) and the Spanish Ministerio de Economia y Competitividad (Grants No. FIS2016-76471-P and No. MAT2016-78293-C6-4-R). We are indebted to the Brazilian agencies CAPES, CNPq, and FAPERGS for partial support of the experimental part of this work.

Control of electron transport through molecular devices is a fundamental step toward design of functional molecular electronics. In this respect, the application of the field-effect transistor principle to molecular junctions appears to be a desirable strategy. Here we study the possibility of mechanically controlling the molecular orbital alignment in self-assembled monolayers via the electrostatic fields originating from dipoles at the metal-molecule interfaces. More specifically, we analyze first-principles simulations of prototype alkanedithiolate and alkanediamine monolayer junctions between Au(111) electrodes as a function of inclination of the molecules in the film. We find that the molecular orbital alignment and hence the low-bias conductance of the junctions, sensitively depends on the interface dipole. The dipole change with molecular tilt is rationalized in terms of two electrostatic effects: (i) the reorientation of a dipole associated with the anchoring group and (ii) a dipole modification arising from charge redistribution due to the metal-molecule bond. The first effect, dominating for the thiolates, is the desired way to gate the junctions by tuning of the molecular tilt. However, the second effect, equally important for the amines, may hamper the mechanical control of the level alignment because it depends on other geometric details than the tilt angle. Our results thus suggest that mechanical gating by tilt is achievable with molecules and anchoring groups with strong intrinsic dipole moments and well-defined binding geometry rather than with interface dipoles associated with weak and flexible metal-molecule bonds., G.F. acknowledges support from the CSIC JAE-predoc program, cofinanced by the European Science Foundation. We also acknowledge the support of the Basque Departamento de Educación and the UPV/EHU (Grant No. IT-756-13), the Spanish Ministerio de Economía y Competitividad (Grant No. FIS2010-19609-CO2-00), and the ETORTEK program funded by the Basque Departamento de Industria and the Diputación Foral de Gipuzkoa.

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).

Pentacene molecules have recently been observed to form a well-ordered monolayer on the (110) surface of rutile TiO, with the molecules adsorbed lying flat, head to tail. With the geometry favorable for direct optical excitation and given its ordered character, this interface seems to provide an intriguing model to study charge-transfer excitations where the optically excited electrons and holes reside on different sides of the organic-inorganic interface. In this work, we theoretically investigate the structural and electronic properties of this system by means of ab initio calculations and compute its excitonic absorption spectrum. Molecular states appear in the band gap of the clean TiO surface, which enables charge-transfer excitations directly from the molecular HOMO to the TiO conduction band. The calculated optical spectrum shows a strong polarization dependence and displays excitonic resonances corresponding to the charge-transfer states, which could stimulate new experimental work on the optical response of this interface., This work was financed by the Deutsche Forschungsgemeinschaft (DFG) through the SFB 1083 project. MPL, DSP and PK also acknowledges support from Spanish MINECO (Grants No. MAT2013-46593-C6-2-P and MAT2016-78293-C6-4-R) and PK also acknowledges financial support from Fellows Guipuzcoa program from Gipuzkoako Foru Aldundia through the FEDER funding scheme of the European Union.

Graphene nanoribbons (GNRs) are promising components in future nanoelectronics due to the large mobility of graphene electrons and their tunable electronic band gap in combination with recent experimental developments of on-surface chemistry strategies for their growth. Here, we explore a prototype 4-terminal semiconducting device formed by two crossed armchair GNRs (AGNRs) using state-of-the-art first-principles transport methods. We analyze in detail the roles of intersection angle, stacking order, inter-GNR separation, GNR width, and finite voltages on the transport characteristics. Interestingly, when the AGNRs intersect at θ=60°, electrons injected from one terminal can be split into two outgoing waves with a tunable ratio around 50% and with almost negligible back-reflection. The split electron wave is found to propagate partly straight across the intersection region in one ribbon and partly in one direction of the other ribbon, i.e., in analogy with an optical beam splitter. Our simulations further identify realistic conditions for which this semiconducting device can act as a mechanically controllable electronic beam splitter with possible applications in carbon-based quantum electronic circuits and electron optics. We rationalize our findings with a simple model suggesting that electronic beam splitters can generally be realized with crossed GNRs., The authors acknowledge financial support from FP7 FET-ICT “Planar Atomic and Molecular Scale devices” (PAMS) project (funded by the European Commission under Contract No. 610446), the Spanish Ministerio de Economia y Competitividad (MINECO) (Grant No. MAT2013-46593-C6-2-P), the Basque Departamento de Educacion, and the UPV/EHU (Grant No. IT-756-13). N.P. acknowledges financial support from No. EU H2020 Project No. 676598, “MaX: Materials at the eXascale” Center of Excellence in Supercomputing Applications.