Mar 22, 2016 - caused by the doping atoms extends over both shells due to a hybridization of their electronic states, yielding a reduction of the ... senting a drawback in the design of an electronic device, the creation of mobility gaps as a result
Oct 12, 2011 - like structure involving a high density of defects, which complicates the investigation of nitrogen configurations and nitrogen-induced .... Therefore, it is instructive to compare the STS measurement with optical absorption experiment
(GGA) to exchange-correlation potential in the PBE form  was used. The full-lattice optimization for both ... bands with mainly Fe 3d character; these bands form electron pockets centered at M and hole pockets centered at Ð. ....  A.L. Ivanov
cleansing that takes place at interfaces between graphene, hBN and transition metal dichalcogenides. Surface ... achieving ultra-high electronic quality is the self-cleansing process previously reported for graphene on hBN11 and .... graphene's DoS.1
Dec 5, 2016 - as opened, i.e. without the cap which contains the pentagonal defects. The second ... Besides DOS, one more quantity is defined - the local density of states (LDOS). It is DOS ...... Defects in Bilayer Graphene [Internet]. 2013.
Mar 19, 2006 - Department of Physics, University of Rijeka, Omladinska 14, 51000 ..... was done on the cluster of PCs at the University Computing Center ...
Mott state, commensurate Wigner crystal, mixed site-bond 4KF CDW. The order parameteris is ~ cos(2Ï+xÏ/a), its phase Ï is locked by the commensurability energy HU. ~-Ucos(2Ï-2Î±). The potential U comes from the CO, the phase centre shift Î± is th
mean free path at room temperature of â 10Ã is found. The resistivity increases nearly linearly with temperature to 1000K in spite of such a short mean free path that resistivity saturation would be expected. .... ferromagnetic transition exceeds
Apr 27, 2015 - applications. Discovered at the outset of the second world War by F. Bloch and E.M. Purcell, the nuclear magnetic ..... to establish that the Na sites are in an ordered pattern with three sites per unit cell for this Na .... Kittel C.,
the p-d exchange coupling between the local magnetic moments of Mn ions and spins of ... experimental data on the thermoelectric power in GaMnAs at high ...
therefore not surprising that graphene-based electronics holds great promise for future applications. Apart from the ... configuration is a highly flexible platform for investigating the electronic properties of FLG. In these ...... that these struct
Oct 1, 2009 - methods, however, possibly at the cost of a decreased computational accuracy. The aim of this paper is ... Such scaling relations may be helpful when modeling on-going experiments on graphene antidot ..... calculations (TB) were carried
Nov 26, 2009 - Section III presents and discusses some details of the electronic band structure calculations that were used in the subsequent sections. Particular attention is paid to the effects of s-o interaction mentioned above and also to the opt
Jan 5, 2018 - and orbital- and spin-projected density of states. Halogen adsorptions can create the non-magnetic, ferromagnetic or anti-ferromagnetic metals, being mainly determined by concentrations and edge structures. The number of holes per an un
1996 Royal Society Typescript. Printed in Great Britain. 1. TEX Paper. Page 2. 2. S. Sachdev. 2. General considerations at zero temperature. We will phrase our discussion in terms of the following Hamiltonian. H = â â i
compounds and calls for further investigation in a more systematic way. Unfortunately, the determination of the ..... Acknowledgments: We thank the computer centers CINES and [email protected] in Montpellier for their support. ....  B. C. Sales, B. C. Chako
properties. Pristine ByNRs adopt a magnetic ground state, anti-ferromagnetic (AFM) or ferromagnetic (FM) depending on the ribbon width, while pristine BxNRs are non-magnetic (NM). Upon ... From graphene to GNRs, the change is not merely morphological
... inertness, endurance, strong oxidizing power, non-toxicity and lower production cost. ....  D.S. Seo, J.K. Lee, H. Kim, J. Crystal Growth 229, 428 (2001).
from bulk sensitive soft X-ray photoemission that CaAgAs belongs to the family of nodal line semimetals . We study the ... CaAgAs single crystals were grown by the flux method with bismuth as an external flux. Ca, Ag, As, and Bi ... crystals of C
Because of comparable electronegativity with silicon, transition metals tend to form strong covalent bond with it. In realizing phonon-mediated .... nature of bonding among the ions within a compound. These factors, in turn, determine the ... The uni
Nov 27, 2015 - the isolated molecule (3 ÂµB for VPc, 2 ÂµB for CrPc). As. U increases from 4 to 8 eV, the MMs get closer to the isolated molecule value. For Mn ...
systems by the methods of electronic structure theory of solids. The basic tool ... corrections due to strong electronic correlation effects appear in the functional in its simplest HartreeâFock form known as the ..... multiple scattering technique
Apr 20, 2013 - strength [37, 38], or optical signature . In the following we review the electronic properties of bilayer graphene. Section II is an overview of the electronic tight-binding Hamiltonian and resulting band structure describing the l
arXiv:1806.02385v1 [cond-mat.mes-hall] 6 Jun 2018
Electronic Properties of Substitutionally Boron-doped Graphene Nanoribbons on a Au(111) Surface Eduard Carbonell-Sanrom`a,† Aran Garcia-Lekue,∗,‡,¶ Martina Corso,∗,§,‡ Guillaume Vasseur,§ Pedro Brandimarte,‡ Jorge Lobo-Checa,k,⊥ Dimas G. de Oteyza,†,‡,¶ Jingcheng Li,† Shigeki Kawai,# Shohei Saito,@ Shigehiro Yamaguchi,4 J. Enrique Ortega,§ Daniel S´anchez-Portal,§ and Jose Ignacio Pascual∗,†,¶ †CIC nanoGUNE, 20018 Donostia-San Sebasti´an, Spain ‡Donostia International Physics Center, 20018 Donostia-San Sebasti´an, Spain ¶Ikerbasque, Basque Foundation for Science, 48013 Bilbao, Spain §Centro de F´ısica de Materiales CFM/MPC (CSIC-UPV/EHU), 20018 Donostia-San Sebasti´an, Spain kInstituto de Ciencia de Materiales de Arag´on (CSIC-UZ), E-50009 Zaragoza, Spain ⊥Dpto. F´ısica de la Materia Condensada, Univ. de Zaragoza, E-50009 Zaragoza, Spain #International Center for Materials Nanoarchitectonics, National Institute for Materials Science, 1-1, Namiki, Tsukuba, Ibaraki 305-0044, Japan. @Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan 4Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan E-mail: [email protected]; [email protected]; [email protected]
ing their stronger interaction with the metal. We find two boron-rich flat bands emerging as impurity states inside the GNR band gap, one of them particularly broadened after its hybridization with the gold surface states. In addition, the boron atoms shift the conduction and valence bands of the pristine GNR away from the gap edge, and leave unaffected the bands above and below, which become the new frontier bands and have negligible boron character. This is due to the selective mixing of boron states with GNR bands according to their symmetry. Our results depict that the GNRs band structure can be tuned by modifying the separation between di-boron moieties.
June 8, 2018
Abstract High quality graphene nanoribbons (GNRs) grown by on-surface synthesis strategies with atomic precision can be controllably doped by inserting heteroatoms or chemical groups in the molecular precursors. Here, we study the electronic structure of armchair GNRs substitutionally doped with diboron moieties at the center, through a combination of scanning tunneling spectroscopy, angle-resolved photoemission, and density functional theory simulations. Boron atoms appear with a small displacement towards the surface signal1
7-AGNR band gap and their inter-band distance depends on the separation between 2B moieties.
Graphene nanoribbons (GNRs) have emerged as a low dimensional platform for nanoelectronics with promising applications due to the predictive control on their functionality. 1 Most GNRs are predicted to be one-dimensional semiconductors, with a band-gap that depends on their width and orientation. 2,3 However, due to their low dimensions, their electronic properties are extremely sensitive to atomic-scale defects and inhomogeneities appearing during their fabrication. On-surface chemistry strategies allow producing atomically precise nanoribbons with perfectly well-defined widths and edges, using organic precursors and directing their reaction over a metal surface. 4–6 This method succeeded in fabricating nanoribbons with controlled size 4,7–11 and orientation, 12,13 in doping them with chemical groups, 14–20 as well as in producing hybrid structures with enhanced functionality. 21–23 Here we study the electronic structure of 7 atoms wide armchair GNRs (7-AGNRs) substitutionally doped with di-boron moieties at the center of the carbon backbone (2B-7-AGNRs, Scheme 1), adsorbed on a Au(111) surface. In our study, we combine scanning tunneling spectroscopy (STS), angle-resolved photoemission (ARPES), and density functional theory (DFT) simulations to unveil the effect of inserting boron atoms in the electronic structure of the pristine ribbon. The synthesis of these ribbons was previously reported in Refs. 15,16, and a recent report proposed that their bands should exhibit a different topological configuration than bare 7-AGNRs. 24 We find that the di-boron impurities modify the electronic structure of the 7-AGNR in two ways. First, the boron states interact with both valence and conduction bands (VB and CB) of the pristine 7-AGNR, while the VB-1 and CB+1 bands remain unaffected. This causes a band-inversion of the frontier bands. Second, two impurity (boron-rich) flat bands appear above and below the Fermi level (EF ). The unoccupied one, already predicted in Refs. 15,16, is found here as a broad band, strongly mixed with Au(111) states. Both impurity bands lie well inside the
Methods SPM methods. Our STS measurements were performed in an ultra-high vacuum system, composed of a chamber for the sample preparation and an independent chamber hosting a custom-made low temperature (4.8 K) Scanning Tunneling Microscope (STM). We measured the differential conductance using a lock-in demodulation method. Angle-resolved photoemission measurements. ARPES measurements were performed in a separate UHV system equipped with a room temperature STM and a low-energy electron diffraction (LEED) system used to check the sample quality. ARPES data were taken using a Phoibos 150 SPECS high-resolution hemispherical electron analyzer on the sample cooled down to 150 K. He-I (hν = 21.2 eV) radiation was provided by a high intensity UVS300 SPECS discharge lamp coupled to a TMM302 SPECS monochromator. Binding energies values given here are referred to EF . Because ARPES is an ensemble-averaging technique, a measurement of the electronic dispersion along the ribbons requires uniaxially aligned GNRs. Therefore, we used the stepped Au(322) substrate instead of the Au(111), which consists of narrow terraces with the Au(111) orientation with parallel steps separated by ∼1.2 nm.
Scheme 1: Chemical structure of the precursor 9,10-bis(10-bromoanthracen-9-yl)-9,10dihydro-9,10-diboraanthracene (1), and the two-step reaction pathway undergone on the Au(111) surface: the Ulmann-like C-C coupling into extended polymers (2), and the cyclodehydrogenation into the flat di-borondoped graphene nanoribbon 2B-7-AGNRs (3).
tion of the (111) surface of Au. 30 This extended, optimized basis set describes adequately the decay of the metal electron density into vacuum, and improves the description of the Au(111) surface state. A 3 × 1 × 1 Monkhorst-Pack mesh was used for the k-point sampling of the three-dimensional Brillouin zone, where the 3 k-points are taken along the direction of the ribbon. A cutoff of 300 Ry was used for the realspace grid integrations. For the free standing 2B-7-AGNR, we performed DFT calculations (including geometrical optimization) using the same simulation parameters and the same unit cell as for the adsorbed ribbon. This allows us to directly compare the band structures obtained in the isolated and the Au-supported cases. All the DFT results presented in this article are spin-restrictred simmulations, and are qualitatively identical to similar calculations including spin as an additional degree of freedom. So, even for the free ribbon, the system does not show any net spin-polarization according ot our DFT simulations.
This surface is a template for the growth of GNRs uniaxially aligned along the narrow terraces. 25 The dense array of steps have the additional effect of lowering the work function of the Au(111) surface by 0.25 eV. 25 Sample preparation. On both the Au(111) and the Au(322) surfaces, the boron-doped graphene nanoribbons 2B-7-AGNRs, were obtained by following a two-step on-surface reaction of 1 (Scheme 1) precursors, as previously reported in Refs. 4,15,16. First, a clean Au(111) single crystal surface was precovered with precursor 1 and annealed to 200 ◦ C for ∼ 5 minutes to activate their Ullmann-like coupling and polymerization into structure 2. Further annealing to 400 ◦ C (∼3 minutes) induced a cyclodehydrogenation reaction, leading to additional C-C bonds between neighbor anthracene units and the creation of the flat GNR backbone 3. As shown in Scheme 1, the two boron atoms of the precursor appear at the center of the ribbon, a pair every 3-unit cells of the pristine 7-AGNR ribbon, amounting to a total ∼5% of the total C atoms of the ribbon. Theoretical methods. Ab initio calculations were carried out using DFT, as implemented in the SIESTA code. 26 For the description of the 2B-7-AGNR adsorbed on gold, we used a supercell composed of a slab containing 3 layers of the Au(111) surface and the nanoribbon on top. We employed the stacking geometry suggested √ in Ref.√27, where 2 × 2 graphene and Au(111) 3 × 3 R30◦ unit cells are directly matched. Following Ref. 28, the bottom Au surface was passivated with hydrogen atoms to quench the Au(111)’s surface state of the bottom surface, avoiding spurious effects due to the interaction of the surface states on both sides of the slab. Moreover, to avoid interactions between periodic surfaces, we considered a vacuum region of more than 20 ˚ A. The GNR and the top two Au layers were fully relaxed until forces were < 0.01 eV/˚ A. Dispersion interactions were taken into account using the nonlocal optB88-vdW functional. 29 The basis set consisted of double-ζ plus polarization (DZP) orbitals for C, B, H, and bulk Au atoms. Au atoms on the topmost layer were treated with a DZP basis set optimized for the descrip-
Results and Discussion Figure 1a shows an overview of the resulting ribbons on the Au(111) surface. As reported in Refs. 15,16, the on-surface reaction of precursor 1 is complete, and produces extended ribbons with a small internal contrast in STM images. A close-up view on a ribbon segment (Fig. 1b) highlights the distinctive topographic contrast of 2B-7-AGNRs, where borylated regions appear as slight depressions along the ribbon backbone. To enhance the resolution of the ribbon we functionalized the STM tip with a CO molecule 31 and mapped the tunneling current in a constant height image at very low bias. 32 This method achieves high resolution images of the ring structure of the graphenoid backbone (Fig. 1c), which are comparable with constant-height AFM images (Fig. 1d). Interestingly, while the boron-regions are imaged as dark areas in AFM, they can be resolved in COSTM, but with a lower contrast (red arrows in Fig. 1c). Although both imaging methods have some sensitivity to electronic states, they 3
1.0 2.0 3.0
0.0 -0.93 eV
Binding Energy (eV)
Sample Bias (V)
NR -AG 2B-7 ) 2 2 (3 Au
1.6 2.0 0.8
Binding Energy (eV)
dI/dV @ V=-0.6 V (arb. units)
-1 Sample Bias (V)
Boron segment Carbon segment Au (111)
dI/dV (arb. units)
Intensity (arb. units)
Figure 2: (a) dI/dV spectra acquired at the edge of a borylated (red) and of a carbon (blue) segment (spectrum in gray is measured over the Au(111) as a reference). A dI/dV step at ∼1.6 eV is attributed to an unoccupied band. We associate a broad feature between 0.2 and 1 V to the distorted surface state hybridized with an impurity band (see text and Fig. 3). The inset shows a constant current (closed feedback) dI/dV spectra acquired at the edges of the ribbon, resolving a weak onset of conductance at ∼-0.6 V, which is attributed to the valence band. (Open feedback conditions Vs =-1.1 V, It =0.4 nA, Vrms =18 mV, f =751.9 Hz; Inset: closed feedback conditions It =0.5 nA,Vrms =12 mV, f =769 Hz) (b) Constant current dI/dV map measured at the bias value of the dI/dV onset, at V=-0.6 V. Arrows mark the position of the borylated units. (c) Angle-resolved photoemission spectral maps of the bare Au(322) stepped surface and of 2B-7-AGNRs grown on a Au(322), where the ribbons align parallel to the steps. On the right, ARPES spectra at the maximum of the 2B-7-AGNRs band and of the clean surface are plotted.
Figure 1: (a) Overview of the 2B-7-AGNR on Au(111) (Vs =1 V It =170 pA). (b) Small scale STM image of a 2B-7-AGNR where the lower topography regions correspond to the borylated areas (Vs =-150 mV It =410 pA). (c-d) High resolution dI/dV (c) and AFM (d) images taken with a CO functionalized tip in constant height mode. Red arrows highlight the borylated segments. suggest that the borylated regions are lower 15 and, consequently, hint towards a stronger hybridization of boron atoms with the gold surface. Valence band determination. We performed STS measurements on 2B-7-AGNR to determine the impact of the boron substitution on the 7-AGNR electronic structure. Figure 2a shows the dI/dV spectra obtained at the edges of the borylated (red) and carbon (blue) regions of the 2B-7-AGNRs and the reference spectrum of the clean Au(111) surface (grey). The most prominent feature is the presence of
a strong onset at ∼1.6 V, similar to the case of the conduction bands in pristine 7-AGNR. 33,34 Thus, we attribute this onset to the conduction 4
band of the borylated ribbon CBB . The identical alignment of this band than in the pristine case suggests that no significant charge transfer with the gold surface is to be expected in the borylated ribbon. Above this onset the conductance in the carbon segments reaches a plateau, while on the borylated segments we detect a lager density of states at higher energies. In contrast to pristine 7-AGNRs, the onset of the valence band (VB) is not easily distinguishable in dI/dV spectra of 2B-7-AGNRs. 15,16 However, we found evidence of the presence of occupied frontier bands from dI/dV spectra in constant-current (closed feedback) mode. The spectra (inset of Fig. 2a) shows an onset of dI/dV signal at ∼-0.6 V, which extends to bias values below -1 V, both at the carbon and the borylated regions. The dI/dV signal associated to this band appears stronger at the edges and also shows some higher signal over the boron sites (Fig. 2b). However, this is far from representing the band’s real shape due to spurious effects regarding the different extension of the GNR orbitals into the vacuum. 34 As we shall show later, the probable origin of the weak valence signal in dI/dV spectra as compared to the pristine version of 7-AGNR is a band order reversal with respect to the 7-AGNR band structure. This reversal is induced by boron states that selectively mix with the carbon πsystem, as we detailed in Ref. 17. To explore the region of occupied bands of the borylated nanoribbons we also performed angle-resolved photoelectron spectroscopy of 2B-7-AGNRs on Au(322). The phoemission results resolve a clear dispersing band with an effective mass of 0.12 me and its onset at a binding energy of 0.93 eV (Fig. 2c). We note that the work function of this surface is 0.25 eV lower than that of Au(111). 25 Hence, if we consider that for weakly interacting adsorbates the bands shift according to changes in the workfunction, these results would correspond to a band’s onset energy of 0.68 eV on the Au(111) surface. This value is close to the value observed in STS, although the low effective mass of the bands electrons is intriguing. As we shall show later, DFT finds that this band corresponds instead to the VB-1 of the pristine rib-
bon, which is however not visible in STS due to its nodal-plane structure. 9,34,35 Impurity states of 2B-AGNR. The dI/dV spectra in Fig. 2a also show the presence of a broad conductance signal, between EF and CBB , which is absent in the electronic structure of the pristine 7-AGNR. To explore in more detail the origin of this broad feature, we present in Fig. 3b and 3c two stack plots of dI/dV spectra (spectral map) taken across the borylated and the carbon sections of the 2B-7-AGNR, respectively (Figure 3a). The spectral maps highlight the presence of a broad feature above 0.3 V, with maximum intensity at around 1 V, i.e. within the CBB -VBB band-gap, strongly localized at the center of the ribbon, and extending along the ribbon’s axis. Its energy and spatial distribution suggest that it corresponds to a new unoccupied impurity band (UIB) originated by the boron substitution. A similar feature was denominated acceptor band in Ref. 16. Interestingly, the spectral maps of Fig. 3 show an evolution of features attributed to the Au(111) surface state (SS) in the proximity of the ribbon (red arrows in Figures 3b and c). The SS signatures gradually shift and merge with the UIB. This effect suggests a sizable hybridization of boron-rich states with the SS, which is consistent with the broad line width of the UIB. Constant height dI/dV maps at the UIB’s maximum signal (1 V) and the CBB ’s onset show distinct modulations, both of which detail the absence of a nodal plane along the AGNR axis, and thus evidences the even symmetry of these bands across the ribbon (Fig. 3d and 3e). Furthermore, the edges of the ribbon are depleted of states at the borylated segments in both bands, whereas the carbon segments accumulate states, specially for the UIB-derived band. In spite of the hybridization with the SS, the dI/dV maps agree fairly well with the simulated local density of states (LDOS) maps for the unoccupied bands of an ribbon adsorbed on the Au(111) surface (shown in Fig. 4d-4f). Density Functional Theory results. A better understanding of the electronic structure of the boron doped nanoribbons and their interaction with the underlying metal substrate 5
a x (nm)
2.5 2.0 1.5 dI/dV (arb. units)
0 1 Sample bias (V)
0 1 Sample bias (V)
dI/dV (arb. untis)
dI/dV (arb. untis)
Figure 3: (a) STM topography of a 2B-7-AGNR. (Vs =-150 mV It =410 pA) dI/dV spectral maps taken across (b) a boron segment and (c) a carbon segment of the 2B-7-AGNR (dashed lines in (a)), including dI/dV plots at the edge and at the center of the ribbon. The UIB is localized at the center of the ribbon, while the CBB appears stronger at the edges. The Au(111) SS is visible outside the ribbon at -500mV and shifts upward underneath (following the red arrows). Constant height dI/dV maps of (d) the UIB and (e) the CBB , measured by approaching a CO-functionalized tip by 0.8 ˚ A and 1.2 ˚ A, respectively, from the open feedback conditions Vs = 20 mV, It = 100 pA. Red arrows mark the position of the borylated regions (Vrms =14 mV, f =760 Hz). is provided by ab initio simulations. In particular, we performed DFT calculations for 2B-7AGNRs, both free standing and adsorbed on a Au(111) substrate. In the latter case, special care was taken to ensure a correct description of the Shockley state of the metal surface (more details in the methods section above). The optimized geometry of the 2B-7-AGNR/Au(111) system is shown in Figure 4a. The presence of the boron pairs induces a significant deformation of the nanoribbon: the borylated regions appear ∼0.3 ˚ A closer to the metal surface, in close agreement with the results from Ref. 15. This deformation is absent in the free ribbon, and hence, reflects a sizable nanoribbonmetal hybridization. In the following we de-
scribe the main ingredients that characterize the band structure of the 2B-7-AGNRs. In-gap impurity bands: Figure 4b shows the band structure of the free-standing nanoribbon projected onto carbon and boron orbitals, thus reflecting the contribution from different species to each band. Notice that this band structure has been computed using a unit cell that contains two 2B-defects in order to facilitate the comparison with the adsorbed case (Fig. 4a). The most characteristic feature is a pair of bands with large boron contribution, located close to the Fermi level (∼ ±0.3 eV in Fig. 4b). They correspond to an occupied impurity-induced band (OIB) and an unoccupied impurity-induced band (UIB). The pres-
Figure 4: (a) Side view of a 2B-7-AGNR on a Au(111) slab after geometry optimization (dashed red line indicates the obtained profile using the atoms in the GNR backbone, where the borylated region are 0.3 ˚ A closer to the metallic surface). A top view of the relaxed structure is shown on the right. (b) Band structure of a free-standing 2B-7-AGNR along the Γ−X direction (Xdirection as indicated in (a)). Left and right panels present the carbon and boron character of the bands, respectively. The radius of the black (red) circles represent the C (B) contribution of each electronic state. Two clear boron heavy bands appear at ±0.3 eV with respect to EF , which for the free standing nanoribbon is arbitrarily defined to be in the middle of the gap between them. For visualization purposes, the C and B contributions have been multiplied by 1/3 and 3, respectively. (c) Same as in b but for 2B-7-AGNR on top of the gold slab, with the B contributions multiplied by 10. (d-f) Constant LDOS isosurfaces (0.0004 e− /Bohr3 ) of the CBB , CBB +1, and UIB bands for the Au-supported 2B-7-AGNR. The considered energy range of integration for each isosurface is indicated on the right side of panel c. In the caption of Fig. 5 refer the labeling scheme of the bands to those of the pristine ribbon. identified in the high dilution limit in Ref. 17 These two impurity bands can be rationalized by considering that boron states of neighboring 2B sites hybridize with each other and form these flat bands. For very low impurity concentration in a free-standing ribbon (e.g. large separation between 2B dopants, see Supporting Information in Ref. 17), two boron states lie very close to EF with a shape that resembles the impurity bands of the present case. In Fig. 5a we show the evolution of the band structure as the borylated moieties are approached. The two
ence of these two bands at EF was already identified at low doping concentration in Ref. 17 (also shown in Fig. 5). Intriguingly, very little B content is present in the bands immediately below and above these impurity bands, i.e. in the delocalized frontier bands of the 2B7-AGNR, whose origin will be discussed later in more detail. We also find an additional band with large B contribution at ∼ -1.9 eV. As in the case of both in-gap impurity states, this band is rather flat, denoting its large degree of localization at the 2B sites. This band was also 7
Figure 5: (a) Band structure of a free-standing 7-AGNR (left) and of 2B-7-AGNRs with decreasing doping concentrations (from left-to-right) reflected by the increase of the unit cell size containing a single 2B impurity. The size of the red shaded area around a given band is proportional to its boron character. The labels of the 7-AGNR frontier bands VB-1/VB/CB/CB+1, are assigned to the borylated ribbons according to their wave-function’s amplitude isosurfaces, shown in (b), such that CBB =CB+1, CBB +1=CB, VBB =VB-1, and VBB -1=VB, thus defining the reversal of frontier bands explained in the text. The occupied/unoccupied mid-gap impurity-induced bands are labeled as OIB/UIB. The Fermi energy is arbitrarily assigned in the middle of the gap between the highest occupied and lowest unoccupied bands, and all the band structures are aligned with respect to CB+1. (b) Representation of the unit cell (top) and the real part of the wave-functions at Γ corresponding to the frontier and impurity-induced bands marked in (a). The color indicate the phase of the wave function, whereas the surfaces are constant amplitude (0.01) isosurfaces. boron states found at EF for the diluted case (n=11) gradually increase their energy spacing as the boron moieties become closer than 7 unit cells of the pristine 7-AGNR, and develop a weak dispersion. Figure 5b shows the evolution of the wave functions of the localized boron states. As the 2B sites get closer the impurities’ wave functions overlap and become the new impurity bands (OIB and UIB) of the cur-
rent (n=3) 2B-7-AGNR. From these results, we foresee that the OIB-UIB energy spacing can be tuned by adjusting the length of the carbon segments used for the on-surface reaction. Interestingly, the wave-functions for the OIB and UIB of the 2B-7-AGNR (obtained at the Γ point, shown in the second column of Figure 5 b) exhibit an even symmetry across the ribbon for both impurity states, 17 in agreement with 8
dI/dV maps of Fig. 3. At the Γ point (Fig. 5) the OIB has a nodal plane in the middle of the 2B-defect, while UIB presents even parity with respect to the 2B-defect center. This contrasts with the apparent nodal plane observed in the experimental map of the UIB (Fig. 3d). The situation, however, is reversed as we move towards the Brillouin zone boundary (not shown in the figure). Hence, a direct comparison of experimental STS maps with wave-function isosurfaces of periodic systems at a particular k point can lead to errors in the band assignment. We thus plot in Fig. 4(d-f) the computed LDOS maps of the unoccupied bands for the 2B-7AGNR on the Au(111) surface, which can be more directly correlated with the dI/dV maps of Fig. 3. The even symmetry across the ribbon is maintained for the B-rich bands, while the low signal between the two B atoms of the UIB is now clearly reproduced by the simulations. Band-reversal at the band-gap edges: Surprisingly, the original VB of the pristine 7-AGNR is absent in the band structure of the 2B-7-AGNR in Fig. 4b. Instead, the first carbon-rich ribbon state is a dispersing band, slightly below the OIB, which has the characteristic wave-function (plotted at Γ in Fig. 5b) of the VB-1 of the pristine ribbon. This band is not affected by the presence of the boron impurities because its has odd symmetry across the ribbon axis, 17 and presents a clear nodal plane along the axis, where the 2B moieties lie. The negligible interaction of VB-1 with the boron impurities is evident from Fig. 5b, where the corresponding wave-function is identical to that of the pristine ribbon and not affected by the change in the separation between the boron moieties. Additionally, this band has many nodal-planes along the ribbon, which explains why it is difficult to detect in normal STS measurements, as it is the case for the VB-1 of the pristine ribbons. 34 In contrast, the original VB is known to interact strongly with boron states because both have even symmetry across the ribbon. As previously reported, 17 this band is strongly scattered by 2B moieties and confined in the carbon 7-AGNR segments between them, appearing as quantized flat bands at lower energies.
The quantized VB of the pristine ribbon appears in Fig. 5b, for the n=11 case, split in several modes with large B content. The VB band onset shifts down significantly as the distance between 2B moieties (i.e. the confined region) becomes smaller. For the n=3 case, the first VB derived QWS appears in Fig. 4b as a quite flat band at ∼ -1.3 eV with a moderate B content that indicates its hybridization with the 2B impurities. The results in Fig. 5 also confirm that a similar reversal of frontier bands due to symmetryselective band mixing applies to the conduction bands. The original CB of 7-AGNR is even across the ribbon and, thus, interacts with the 2B states, in contrast to the pristine CB+1, which is odd. However, the interaction is weaker in this case, and the pristine CB moves to lie slightly above the CB+1, which becomes the new frontier band. In contrast to the case of valence bands, the pristine CB was in this case the “hidden” band in STS spectra because of the numerous nodal planes along the ribbon, 34 and the one usually detected was the CB+1, which has odd symmetry across and is not affected by the boron. This explain that in both pristine and borylated ribbons, the conduction band detected by STM, the CB+1, lies at the same position. Adsorption on the Au(111) surface: When adsorbed on the Au(111) slab, the main ingredients of the electronic structure of the freestanding 2B-7-AGNRs can be still identified. One of the most significant changes though, is the mixing of the impurity bands with Au states. Figure 4c shows that UIB undergoes a larger distortion and mixture with Au states upon adsorption, resulting in a broad band of boron-rich states that accounts for its spectral line shape in our experiment. In contrast, the dispersion of the carbon-rich bands and the OIB can still be recognized in the band plots, signaling their weaker hybridization with the states at the surface. We observe a small mixing between the originally orthogonal VB-1 and OIB, probably due to their proximity and the mediation of the metal. Since the OIB is a rather flat band, we conclude that the band resolved by ARPES in Fig.2c corre9
sponds to the VB-1, while the one measured by STM in Fig. 2a is probably the OIB, because it appears modulated with the periodicity of the boron subunits.
search Council (grant agreement no. 635919), and the European Regional Development Fund.
(1) Celis, A.; Nair, M. N.; Taleb-Ibrahimi, A.; Conrad, E. H.; Berger, C.; De Heer, W. A.; Tejeda, A. Graphene nanoribbons: Fabrication, properties and devices. J. Phys. D. Appl. Phys. 2016, 49, 143001.
In conclusion, we have studied the structural and electronic effects of a boron-substitution on armchair graphene nanoribbons. The resulting structures exhibit a high boron-dopant concentration and a well-defined periodicity. The most prominent change on the electronic structure of 2B-7-AGNRs is the introduction of two in-gap impurity bands, which arise from the overlap of nearby borylated moieties, and mix strongly with Au(111) surface states resulting in hybrid metal-ribbon states. The 2B dopants also cause a reversal of the carbon-rich bands at the gap edges, compared to the pristine 7-AGNR ribbon. This is caused by the symmetry selective mixing of bands with the states of the diboron defect, which leads to quantum confinement of original VB and CB and their shift away from the gap. The new valence and conductance bands have a different symmetry both across and along the ribbon. Indeed, they correspond to the VB-1 and CB+1 of the pristine 7-AGNR. This is the reason for the previously reported 15,16 disappearance of the VB from STS spectra. Our results further define basic rules for tuning the band structure of the ribbons: varying the distance between 2B groups can be employed to tune the gap between impurity bands; rotation of the 2B moiety with respect to the ribbon axis changes the symmetry of the carbon bands that they mix with, thus affecting the character of the frontier transport bands.
(2) Yang, L.; Park, C.-H. H.; Son, Y.-W. W.; Cohen, M. L.; Louie, S. G. Quasiparticle energies and band gaps in graphene nanoribbons. Phys. Rev. Lett. 2007, 99, 6–9. (3) Wakabayashi, K.; Sasaki, K.-i.; Nakanishi, T.; Enoki, T. Electronic states of graphene nanoribbons and analytical solutions. Sci. Technol. Adv. Mater. 2010, 11, 054504. (4) Cai, J.; Ruffieux, P.; Jaafar, R.; Bieri, M.; Braun, T.; Blankenburg, S.; Muoth, M.; Seitsonen, A. P.; Saleh, M.; Feng, X.; Mullen, K.; Fasel, R. Atomically precise bottom-up fabrication of graphene nanoribbons. Nature 2010, 466, 470–473. (5) Talirz, L.; Ruffieux, P.; Fasel, R. OnSurface Synthesis of Atomically Precise Graphene Nanoribbons. Adv. Mater. 2016, 6222–6231. (6) Corso, M.; Carbonell-Sanrom`a, E.; de Oteyza, D. In On-surface synthesis II ; de Oteyza, D., Rogero, C., Eds.; Springer: Switzerland, 2018; Chapter Bottomup Fabrication of Atomically Precise Graphene Nanoribbons, pp 113–152.
Acknowledgement We acknowledge financial support from the Basque Government (Departamento de Industria, grant no. PI2015-1-42, and Departamento de Educaci´on and UPV/EHU, grant No. IT-756-13), and the Spanish Ministerio de Econom´ıa y Competitividad (Grants No. MAT2016-78293-C6, FIS2015-62538-ERC, FIS2017-83780-P), and the Maria de Maeztu Units of Excellence Programme MDM-2016-0618), the European Re-
(8) Chen, Y.-C.; de Oteyza, D. G.; Pedramrazi, Z.; Chen, C.; Fischer, F. R.; Crommie, M. F. Tuning the Band Gap of Graphene Nanoribbons Synthesized from Molecular Precursors. ACS Nano 2013, 7, 6123–6128.
N = 13 Sulfur-Doped Graphene Nanoribbons. J. Phys. Chem. C 2016, 120, 2684– 2687. (15) Kawai, S.; Saito, S.; Osumi, S.; Yamaguchi, S.; Foster, A. S.; Spijker, P.; Meyer, E. Atomically controlled substitutional boron-doping of graphene nanoribbons. Nat. Commun. 2015, 6, 8098.
(9) Talirz, L. et al. On-Surface Synthesis and Characterization of 9-Atom Wide Armchair Graphene Nanoribbons. ACS Nano 2017, 1380–1388.
(16) Cloke, R. R.; Marangoni, T.; Nguyen, G. D.; Joshi, T.; Rizzo, D. J.; Bronner, C.; Cao, T.; Louie, S. G.; Crommie, M. F.; Fischer, F. R. SiteSpecific Substitutional Boron Doping of Semiconducting Armchair Graphene Nanoribbons. J. Am. Chem. Soc. 2015, 137, 8872–8875.
(10) Abdurakhmanova, N.; Amsharov, N.; Stepanow, S.; Jansen, M.; Kern, K.; Amsharov, K. Synthesis of wide atomically precise graphene nanoribbons from para-oligophenylene based molecular precursor. Carbon 2014, 77, 1187 – 1190.
(17) Carbonell-Sanrom`a, E.; Brandimarte, P.; Balog, R.; Corso, M.; Kawai, S.; GarciaLekue, A.; Saito, S.; Yamaguchi, S.; Meyer, E.; S´anchez-Portal, D.; Pascual, J. I. Quantum Dots Embedded in Graphene Nanoribbons by Chemical Substitution. Nano Lett. 2017, 17, 50–56.
(11) Merino-D´ıez, N.; Garcia-Lekue, A.; Carbonell-Sanrom`a, E.; Li, J.; Corso, M.; Colazzo, L.; Sedona, F.; S´anchezPortal, D.; Pascual, J. I.; De Oteyza, D. G. Width-Dependent Band Gap in Armchair Graphene Nanoribbons Reveals Fermi Level Pinning on Au(111). ACS Nano 2017, 11, 11661–11668.
(18) Carbonell-Sanrom`a, E.; Hieulle, J.; VilasVarela, M.; Brandimarte, P.; Iraola, M.; Barrag´an, A.; Li, J.; Abadia, M.; Corso, M.; S´anchez-Portal, D.; Pe˜ na, D.; Pascual, J. I. Doping of Graphene Nanoribbons via Functional Group Edge Modification. ACS Nano 2017, 11, 7355– 7361.
(12) de Oteyza, D. G.; Garc´ıa-Lekue, A.; VilasVarela, M.; Merino-D´ıez, N.; CarbonellSanrom`a, E.; Corso, M.; Vasseur, G.; Rogero, C.; Guiti´an, E.; Pascual, J. I.; Ortega, J. E.; Wakayama, Y.; Pe˜ na, D. Substrate-Independent Growth of Atomically Precise Chiral Graphene Nanoribbons. ACS Nano 2016, 10, 9000–9008.
(19) Durr, R. A.; Haberer, D.; Lee, Y.L.; Blackwell, R.; Kalayjian, A. M.; Marangoni, T.; Ihm, J.; Louie, S. G.; Fischer, F. R. Orbitally Matched EdgeDoping in Graphene Nanoribbons. J. Am. Chem. Soc. 2018, 140, 807–813.
Dispersive Hybrid Bands at an OrganicMetal Interface. Phys. Rev. Lett. 2008, 100, 156805.
(22) Chen, Y.-C.; Cao, T.; Chen, C.; Pedramrazi, Z.; Haberer, D.; de Oteyza, D. G.; Fischer, F. R.; Louie, S. G.; Crommie, M. F. Molecular bandgap engineering of bottom-up synthesized graphene nanoribbon heterojunctions. Nat. Nanotechnol. 2015, 10, 156–160.
(29) Klimeˇs, J.; Bowler, D. R.; Michaelides, A. Chemical accuracy for the van der Waals density functional. J. Phys. Cond. Matt. 2010, 22, 022201. (30) Garc´ıa-Gil, S.; Garc´ıa, A.; Lorente, N.; Ordej´on, P. Optimal strictly localized basis sets for noble metal surfaces. Phys. Rev. B 2009, 79, 075441.
(23) Li, J.; Merino-D´ıez, N.; CarbonellSanrom`a, E.; Vilas-Varela, M.; de Oteyza, D. G.; Pe˜ na, D.; Corso, M.; Pascual, J. I. Survival of spin state in magnetic porphyrins contacted by graphene nanoribbons. Sci. Adv. 2018, 4, eaaq0582.
(31) Gross, L.; Moll, N.; Mohn, F.; Curioni, A.; Meyer, G.; Hanke, F.; Persson, M. HighResolution Molecular Orbital Imaging Using a p -Wave STM Tip. Phys. Rev. Lett. 2011, 107, 086101.
(24) Cao, T.; Zhao, F.; Louie, S. G. Topological Phases in Graphene Nanoribbons: Junction States, Spin Centers, and Quantum Spin Chains. Phys. Rev. Lett. 2017, 119, 076401.
(32) Hieulle, J.; Carbonell-Sanrom`a, E.; VilasVarela, M.; Garcia-Lekue, A.; Guiti´an, E.; Pe˜ na, D.; Pascual, J. I. On-Surface Route for Producing Planar Nanographenes with Azulene Moieties. Nano Lett. 2018, 18, 418–423.
(25) Merino-D´ıez, N.; Li, J.; Garcia-Lekue, A.; Vasseur, G.; Vilas-Varela, M.; CarbonellSanrom`a, E.; Corso, M.; Ortega, J. E.; Pe˜ na, D.; Pascual, J. I.; de Oteyza, D. G. Unraveling the Electronic Structure of Narrow Atomically Precise Chiral Graphene Nanoribbons. J. Phys. Chem. Lett. 2018, 9, 25–30.
(33) Ruffieux, P.; Cai, J.; Plumb, N. C.; Patthey, L.; Prezzi, D.; Ferretti, A.; Molinari, E.; Feng, X.; M¨ ullen, K.; Pignedoli, C. A.; Fasel, R. Electronic Structure of Atomically Precise Graphene Nanoribbons. ACS Nano 2012, 6, 6930– 6935. (34) S¨ode, H.; Talirz, L.; Gr¨oning, O.; Pignedoli, C. A.; Berger, R.; Feng, X.; M¨ ullen, K.; Fasel, R.; Ruffieux, P. Electronic band dispersion of graphene nanoribbons via Fourier-transformed scanning tunneling spectroscopy. Phys. Rev. B 2015, 91, 045429.
(26) Soler, J. M.; Artacho, E.; Gale, J. D.; Garc´ıa, A.; Junquera, J.; Ordej´on, P.; S´anchez-Portal, D. The SIESTA method for ab initio order- N materials simulation. J. Phys. Condens. Matter 2002, 14, 2745. (27) Giovannetti, G.; Khomyakov, P. A.; Brocks, G.; Karpan, V. M.; van den Brink, J.; Kelly, P. J. Doping Graphene with Metal Contacts. Phys. Rev. Lett. 2008, 101, 026803.
(35) Senkovskiy, B. V.; Usachov, D. Y.; Fedorov, A. V.; Haberer, D.; Ehlen, N.; Fischer, F. R.; Gr¨ uneis, A. Finding the hidden valence band of N=7 armchair graphene nanoribbons with angle-resolved photoemission spectroscopy. 2D Mater. 2018, 5, 035007.
(28) Gonzalez-Lakunza, N.; Fern´andezTorrente, I.; Franke, K. J.; Lorente, N.; Arnau, A.; Pascual, J. I. Formation of