Borophene Synthesis on Au(111) - ACS Nano (ACS Publications)

Mar 7, 2019 - ... J. Mannix , Brandon L. Fisher , Boris I. Yakobson , Mark C Hersam , and Nathan ... Field-Effect Transistors for Large-Area Flexible ...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF TEXAS DALLAS

Article

Borophene Synthesis on Au(111) Brian Kiraly, Xiaolong Liu, Luqing Wang, Zhuhua Zhang, Andrew J. Mannix, Brandon L. Fisher, Boris I. Yakobson, Mark C Hersam, and Nathan P Guisinger ACS Nano, Just Accepted Manuscript • Publication Date (Web): 07 Mar 2019 Downloaded from http://pubs.acs.org on March 7, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Borophene Synthesis on Au(111) AUTHORS Brian Kiraly1,2†, Xiaolong Liu3†, Luqing Wang4†, Zhuhua Zhang5†, Andrew J. Mannix1,2, Brandon L. Fisher1, Boris I. Yakobson4*, Mark C. Hersam2,3,6*, and Nathan P. Guisinger1* AFFILIATIONS 1Center for Nanoscale Materials, Argonne National Laboratory, 9700 South Cass Avenue, Building 440, Argonne, IL 60439, USA 2Department of Materials Science and Engineering, Northwestern University, 2220 Campus Drive, Evanston, IL 60208, USA 3Applied Physics Graduate Program, Northwestern University, 2220 Campus Drive, Evanston, IL 60208, USA 4Department of Materials Science and NanoEngineering and Department of Chemistry, Rice University, Houston, Texas 77005, USA 5State Key Laboratory of Mechanics and Control of Mechanical Structures and Institute of Nano Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China 6Department of Chemistry, Northwestern University, 2220 Campus Drive, Evanston, IL 60208, USA †

Equal Contribution: Kiraly, Liu, Wang, and Zhang are co-first authors Authors: [email protected]; [email protected]; [email protected] *Corresponding

1 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT Borophene (the first two-dimensional (2D) allotrope of boron) is emerging as a groundbreaking system for boron-based chemistry and, more broadly, the field of low-dimensional materials. Exploration of the phase space for growth is critical because borophene is a synthetic 2D material that does not have a bulk layered counterpart and thus cannot be isolated via exfoliation methods. Herein, we report synthesis of borophene on Au(111) substrates. Unlike previously studied growth on Ag substrates, boron diffuses into Au at elevated temperatures and segregates to the surface to form borophene islands as the substrate cools. These observations are supported by ab-initio modeling of interstitial boron diffusion into the Au lattice. Borophene synthesis also modifies the surface reconstruction of the Au(111) substrate, resulting in a trigonal network that templates growth at low coverage. This initial growth is comprised of discrete borophene nanoclusters, whose shape and size are consistent with theoretical predictions. As the concentration of boron increases, nano-templating breaks down and larger borophene islands are observed. Spectroscopic measurements reveal that borophene grown on Au(111) possesses metallic electronic structure, suggesting potential applications in 2D plasmonics, superconductivity, interconnects, electrodes, and transparent conductors.

Keywords: borophene, two-dimensional materials, Au(111), scanning tunneling microscopy, synthesis, allotrope

2 ACS Paragon Plus Environment

Page 2 of 23

Page 3 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

As the lightest valence-three element, boron possesses chemical properties that tend to favor complex bonding arrangements.1-8 From computational modeling to experimental synthesis, the recent discovery of two-dimensional (2D) boron polymorphs, collectively referred to as borophene,7 has further elaborated the richness of boron chemistry and ushered in a class of low-dimensional materials.916

Borophene is an elemental 2D metal with superlative theoretical properties

including a tensile strength comparable to graphene and a relatively high superconducting transition temperature.9,17 The initial breakthrough for borophene synthesis occurred using Ag(111) substrates, although theoretical modeling predicts stable low-energy structures on additional metal surfaces (e.g., Au and Cu). 11 Here, we present ultra-high vacuum (UHV) synthesis and characterization of borophene on Au(111). Borophene growth begins when boron from the bulk segregates to surface, while the sample is cooled. Growth on Au(111) is in stark contrast to boron limited to the surface on Ag(111). 9,11 Furthermore, the energy minimization and strain relief of the Au(111) surface results in nanoscale templating of borophene nucleation and growth. Increasing boron coverage results in a breakdown of the nanoscale templating mechanism, allowing the borophene to evolve from small ordered islands to larger sheets. The Au(111) surface is among the few low-energy metal surfaces that undergoes a reconstruction to minimize surface energy. This reconstruction is characterized by a 22 x √3 structure consisting of surface atoms located at facecentered cubic (fcc), hexagonal close packing (hcp), and bridge-type locations with respect to the bulk lattice.18,19 The reconstruction results in a strain-mediated 3 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

herringbone pattern of Shockley partial dislocations with uniform spacing and domain orientations of 120°. Adsorption experiments reveal that this surface reconstruction influences the growth of weakly-interacting materials on Au(111), typically leading to nucleation at the corners of the herringbone reconstruction.2025

Strongly interacting adsorbates induce more significant modifications to the

herringbone reconstruction, including altering the symmetry to a highly ordered trigonal network.26-28 RESULTS AND DISCUSSION The pristine Au(111) surface, shown in the UHV scanning tunneling microscopy (STM) image of Fig. 1a, possesses the herringbone surface reconstruction that consists of Shockley partial dislocations running along the [112] directions that separate the fcc from hcp striped regions. Deposition of boron from an atomic source was studied for a series of substrate temperatures (see Supporting Fig. 1). Low temperature deposition results in boron clusters confined to the surface. At substrate temperatures of ~550 °C, boron clusters are no longer present and a transformation is observed from the conventional herringbone reconstruction to a trigonal network where small nanoscale borophene islands emerge, as seen in Fig. 1b (see also Supporting Fig. 1). Atomic-resolution imaging of the Au(111) surface confirms that the trigonal strain-relief network forms within the Au surface, while borophene islands with a calculated v1/12 structure nucleate at the nodes of the network. 11 The total dose of atomic boron needed to observe nanoscale borophene islands is an order of magnitude greater than what is needed to form a monolayer on Ag(111), which is the first indication that boron may be

4 ACS Paragon Plus Environment

Page 4 of 23

Page 5 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

dissolving into the bulk at these substrate temperatures. Figure 1(c) illustrates the computationally predicted atomic structure model for the borophene v1/12 on a Au(111) substrate. The trigonal network has a characteristic periodicity ranging from 5.5 nm to 8.0 nm and acts as a growth template for nanoscale borophene islands. Boron is not expected to form a stable alloy with Au, and we do not observe direct evidence of alloying. However, the trigonal network has been observed for clusters that are incorporated into the surface.26-28 To minimize surface energy, this reconstruction consists of alternating triangular regions of fcc and hcp stacked Au atoms,26-28 as indicated by the triangular dashed regions in Fig. 1b. The apparent embedded nature of the boron islands is likely responsible for the transition from the normal herringbone surface reconstruction to the trigonal network. The chemical properties of the boron-modified Au(111) surface are investigated by in situ X-ray photoelectron spectroscopy (XPS) and compared to that of a clean Au(111) surface, as shown in Fig. 2 and Supporting Fig. 2. The Au 4f7/2 core-level peak after room temperature boron deposition matches that of clean Au. However, the trigonal network samples show a 0.1 eV redshift (see Supporting Fig. 2). This shift is consistent with 1-2% strain-relief induced by the incorporation of boron into the Au(111) lattice and removal of Au atoms.29-31 Indeed, the transformation of the strain relief pattern from the herringbone reconstruction to the trigonal network indicates a change of strain status of the surface. The Au surface with boron particles also shows a clear peak in the B1s core-level (Fig. 2a), which is expected given limited solubility of boron in Au at

5 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

room temperature. However, the reduced peak height from the trigonal network samples indicate that the surface coverage of boron is much lower in these cases, which suggests that boron diffuses into the subsurface region at higher substrate temperatures. The absence of obvious peak shift of the B1s core level, albeit low in intensity, suggests minimum valence changes, in contrast to the case of boron oxides (peak shifts to 192 eV.9 Increasing the boron dose eventually results in the breakdown of the trigonal network and the formation of larger borophene islands that are embedded in the topmost layer of Au, as illustrated in Fig. 2b. The dashed hexagons in the insets of Fig. 2b illustrate the trigonal network. Borophene sheets increase in size with boron dose, leading to domains that extend at least several tens of nanometers at the highest doses in this study as seen in the bottom two panels of Fig. 2b. Both STM and XPS (Fig. 2 and Supporting Fig. 2) reveal that boron clusters form on the surface at low temperature, while boron appears to dissolve into the bulk at elevated temperatures. Thus, we suggest that the borophene growth results from surface segregation of boron from the bulk when the sample is cooled (illustrated in Fig. 3a). This surface segregation of boron to form borophene on Au(111) is contrary to the surface-only growth on Ag(111); these two schemes are analogous to the two primary growth modes observed for graphene on metal substrates.32 To gain further insight into the mechanisms underlying the apparent surface segregation and recessed growth modes of boron on Au(111), we performed density functional theory (DFT) calculations to quantify the adsorption energetics

6 ACS Paragon Plus Environment

Page 6 of 23

Page 7 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

of boron atoms on this surface. We do not compute all of the atom-state sequences, as in the nanoreactor diagram,33 but only key initial configurations. Figure 3c shows the resulting DFT predictions of the relaxed energies of boron atoms at different sites within and on top of Au(111). The calculations indicate that while both bulk solution and surface locations are energetically favorable, they are significantly higher in energy than the near-surface dissolution of boron atoms. In particular, the most stable configurations of isolated boron atoms is in subsurface sites in the topmost layer, with sites along Au(111) step edges being favored by ~0.3 eV over planar subsurface sites. These energetics differ significantly from those of boron on Ag(111).9,11 The surface strain is greater in Au, up to 5.5%,34 than in Ag so that inlayer growth is preferred on Au. Despite their similar chemical properties, the noble metal surfaces have strongly differing electronic densities near the Fermi level (Ef), which leads to the stabilization of distinct boron configurations. Importantly, subsurface boron undergo a less steric effect from surrounding Au and in turn can form a stronger Au-B interaction to adopt a significantly more energetically favorable adsorption site. Modeling of the boron diffusion energetics within the first Au(111) layer shows further consistency with the experimental observation of embedded islands. For example, Fig. 3c (right panel) illustrates the energetics of multi-atom boron clusters on the Au(111) surface. As shown with respect to the energy of a single subsurface boron atom, the calculations favor the formation of small boron clusters within the Au(111) subsurface. The increasing depth of the potential well for dimers and subsequently trimers indicates the preference for boron-boron bond formation

7 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

within this subsurface environment. Furthermore, the most stable configuration involves a Au atom being expelled out of the surface to accommodate the boron trimer cluster, attributed largely to further relief of surface stress on Au(111). The nucleated boron trimer then extends into 2D islands upon further docking of boron atoms and the removal of redundant Au atoms. Again, this behavior marks a significant difference from the 2D surface growth motif of borophene overlayers on Ag(111), 9,11 confirming our experimental observations. To experimentally confirm that boron is dissolving into the bulk, the distribution of boron in Au is further investigated via time-of-flight secondary ion mass spectrometry (ToF-SIMS). In addition to a trigonal network sample and a boron cluster sample where boron is deposited on Au held at room temperature, an as-received and a UHV cleaned Au samples are also explored as control samples. In Fig. 3d and Supporting Fig. 3, higher boron signals are seen in the SIMS spectra for the boron cluster sample compared to the trigonal network sample, suggesting a higher surface concentration of boron. This observation is in agreement with the XPS B1s level results. Depth-profiling of the four samples reveals that the boron concentration in the trigonal network sample only decays after ~30 nm into the bulk, whereas the boron signals decay more quickly in the three other cases once the surface is sputtered (Supporting Fig. 3). In agreement with the XPS results, the SIMS data suggest that boron diffuses deeply into the Au subsurface region at elevated temperatures. The UHV cleaned Au shows no obvious differences compared to as-received Au, indicating negligible effect of UHV annealing.

8 ACS Paragon Plus Environment

Page 8 of 23

Page 9 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

On larger borophene islands, STM faintly resolves a herringbone reconstruction in the Au directly below the 2D sheets. The zoomed-in image of a borophene island (Fig. 4a) reveals apparent atomic-scale periodicity, while herringbone stripes concurrently remain visible (black arrows) below the 2D sheet. The ability to resolve the underlying Au(111) herringbone reconstruction suggests that the boron islands are atomically thin. The STM image of Fig. 4b further reveals an apparent atomic-scale periodicity within the borophene sheet. However, this periodicity is distorted by a larger superstructure that is likely related to a Moiré pattern between borophene and Au(111). A fast Fourier transform (FFT) of the STM image in Fig. 4b shows a strong periodicity of ~0.66 nm as is illustrated in Fig. 4c. This ~0.66 nm periodicity is consistently observed in the extended borophene sheets. At the early stages of nucleation and growth, the borophene islands are templated at the nanoscale. As seen in the STM images in Fig. 1, the specific shape of the borophene islands is both compact and anisotropic. Closer inspection reveals that each island can be approximately decomposed into rhombohedral units, as highlighted in Figs. 4d and 4e. Furthermore, the coalescence of multiple rhombohedral building blocks can account for the shape of larger borophene islands. The size of an individual rhombohedral unit (Figs. 4d and 4e) ranges between 0.9 and 1.4 nm2 with an aspect ratio of 1.5 +/- 0.2, where we note that the accuracy of this measurement is limited by the gold adatoms that are often found at the edges of the borophene islands. A histogram generated from over 500 islands (Fig. 4e) reveals clear discretization of island area at integer values of

9 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

square nanometers. These rhombohedral units, both in terms of size and shape, are consistent with investigations of planar size-selected boron clusters.7, 35-38 The smaller borophene islands are often at the single nanometer length scale, where strong electronic confinement effects are likely to occur. Furthermore, their discrete rhombohedral building blocks may show quantum effects. To further understand the electronic properties of the borophene islands, a series of biasdependent images and concurrently acquired dI/dV tunneling conductance maps (proportional to the electronic density of states) are presented in Fig. 5a. The biasdependent conductance maps reveal that not only is the borophene electronically distinct from the Au but that the fcc (red dashed triangle) and hcp (yellow dashed triangle) regions are also distinct, where the contrast inverts between -0.2 V and 0.4 V sample bias. At positive sample bias, the Au surface state becomes pronounced as we observe the standing wave reflecting from the step-edge and interference patterns between the borophene islands, further confirming the embedded nature of the borophene islands. The borophene islands themselves exhibit electronic variation at the nanometer length scale (particularly at -0.2 V or +0.4 V), which could be related to electronic confinement effects. DFT calculations reveal that the previously predicted borophene structure on Au(111) is metallic, as illustrated by the band structure and density of states (DOS) in Fig. 5b.11 Scanning tunneling point spectra verify that the 2D borophene on Au(111) is indeed metallic as shown in the blue curve in Fig. 5c. It is noteworthy that the spectra for the large borophene island on Au(111) are similar to the spectra seen for borophene on Ag(111).9 In contrast, the gray curve is taken on the clean

10 ACS Paragon Plus Environment

Page 10 of 23

Page 11 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Au(111) surface and shows the characteristic Shockley surface state at approximately -0.52 eV.39 Spectra taken on the borophene nano-islands differ significantly from larger borophene sheets, with several additional spectral features. For example, in the black curve in Fig. 5c, a shoulder is observed at approximately +0.3 V in addition to a symmetric inflection at -0.3 V (indicated by the dashed lines). We believe that these shoulders are a result of additional electron confinement within the borophene clusters.

Electron confinement in

quantum dot structures has previously been shown to result in symmetric states around the Fermi level.40,41 CONCLUSIONS In conclusion, nanoscale borophene structures have been synthesized and characterized on Au(111) in UHV. The observed growth mode is consistent with boron dissolving into the bulk at high temperatures and segregating to the surface to form borophene as the sample is cooled. The initial stages of growth for low concentrations of boron result in nanoscale templating of the borophene islands across the surface, as the Au(111) reconstructs in a trigonal network to minimize energy. At low boron concentration, the borophene islands are comprised of 1 to 8 rhombohedral units with a characteristic area of 1 nm2, as observed in statistical analysis of the island size. The electronic properties of borophene on Au(111) is confirmed to be metallic, with the nanoscale borophene islands showing evidence of electron confinement. Overall, the nanoscale templated growth of borophene on Au(111) presents opportunities for studying and manipulating the electronic properties of 2D metals.

11 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

METHODS The Au(111) on mica substrates were commercially purchased from Keysight Technologies. The Au(111) surface was cleaned in UHV via repeated sputter/anneal cycles until atomically pristine. Characterization was performed using a variable temperature STM (Scienta Omicron) and a home-built UHV STM system adopting a Lyding-design microscope. 42,43 STM images and spectroscopy in the main text were gathered at a sample temperature of 55 K with electrochemically etched W tips. Borophene was grown via thermal deposition from a pure boron rod using electron-beam evaporators (Focus EFM and SPECS EBE-1). The substrate temperature during borophene growth was ~550C. The evaporation rate is ~ 0.02 ML/min (based on borophene growth on Ag(111) with the same power), and an deposition time ranging from 30 min to 3 hrs. DFT calculations were performed using ultrasoft pseudo-potentials for the core region and generalized gradient approximation of Perdew–Burke–Ernzerhof (PBE) functional for describing the exchange correlation, as implemented in VASP code.44,45 A kinetic energy cutoff of 400 eV is chosen for the plane-wave expansion. In all structures, the vacuum region between two adjacent periodic images is fixed to 10 Å to eliminate spurious interaction. The Brillouin-zone integration is densely sampled according to the supercell size, ensuring approximately the same k-point density among different-sized supercells. The k-point density is approximately 45 Å-1. For example, for a supercell with a size of 10 Å×8.65 Å, the k-point sampling is 5×5×1, while “1” is for the vacuum layer.The metal substrate consists of four atomic layers, with atoms in the bottom layer being

12 ACS Paragon Plus Environment

Page 12 of 23

Page 13 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

fixed to the bulk lattice. The positions of all other atoms are relaxed using the conjugate-gradient method until the force on each atom is less than 0.01 eV/Å. The minimum-energy path was mapped out using the climbing image-nudged elastic-band method. 46 XPS spectra were taken with an Omicron DAR 400M X-ray source, XM 500 X-ray monochromator, and EA 125 energy analyzer at a base pressure of ~4 × 1010

Torr. The XPS system is integrated with the home-built UHV STM system for in

situ measurements with an energy resolution of ~0.6 eV. ToF-SIMS spectra were acquired on a PHI TRIFT III ToF-SIMS (Physical Electronics) system with a liquid metal ion gun in the positive ion mode. The analysis chamber pressure was ~1×10-9 Torr.

ASSOCIATED CONTENT Supporting information contains boron dosing information as a function of temperature, additional XPS data, and ToF-SIMS data.

13 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACKNOWLEDGMENTS Use of the Center for Nanoscale Materials, an Office of Science user facility, was supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. MCH and XL acknowledge funding from the Office for Naval Research (N00014-17-1-2993). Work at Rice (ZZ and BIY) was supported by the DOE grant DE-SC0012547. COMPETING INTERESTS The authors declare no conflict of interest. DATA AVAILABILITY STATEMENT The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. CONTRIBUTIONS BK, AJM, MCH, and NPG conceived the project; BK, XL, and AJM carried out the experiment and data analysis. ZZ, LW, and BIY provided theoretical support and computational modeling. BLF provided engineering support for the experiment. NPG, BK, ZZ, XL, and MCH wrote and revised the manuscript. All authors contributed to the manuscript.

14 ACS Paragon Plus Environment

Page 14 of 23

Page 15 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

REFERENCES 1. Ogitsu, T.; Schwegler, E.; Galli, G. -rhombohedral Boron: at the Crossroads of the Chemistry of Boron and the Physics of Frustration. Chem. Rev. 2013,113, 3425–3449. 2. Douglas, B.; Ho, S.-M. Structure and Chemistry of Crystalline Solids (Springer Science & Business Media, New York, 2007). 3. Oganov, A. R.; Chen, J.; Gatti, C.; Ma, Y.; Ma, Y.; Glass, C. W.; Liu, Z.; Yu, T.; Kurakevych, O. O.; Solozhenko, V. L. Ionic High-pressure Form of Elemental Boron. Nature 2009, 457, 863–867. 4. Zhai, H.-J.; Kiran, B.; Li, J.; Wang L.-S. Hydrocarbon Analogues of Boron Clusters – Planarity, Aromaticity, and Antiaromaticity. Nat. Mater. 2003, 2, 827– 833. 5. Sergeeva, A. P.; Popov, I. A.; Piazza, Z. A.; Li, W. -L.; Romanescu, C.; Wang, L. -S.; Boldyrev, A. I. Understanding Boron Through Size-Selected Clusters: Structure, Chemical Bonding, and Fluxionality. Acc. Chem. Res. 2014, 47, 1349– 1358. 6. Zhai, H.-J.; Zhao, Y. -F.; Li, W. -L.; Chen, Q.; Bai, H.; Hu, H. -S.; Piazza, Z. A.; Tian, W. -J.; Lu, H. -G.; Wu, Y. -B.; Mu, Y. -W.; Wei, G. -F.; Liu, Z. -P.; Li, J.; Li, S. -D. Observation of an All-boron Fullerene.” Nat. Chem. 2014, 6, 727–731. 7. Piazza, Z. A.; Hu, H. -S.; Li, W. -L.; Zhao, Y. -F.; Li, J. Planar Hexagonal B36 as a Potential Basis for Extended Single-Atom Layer Boron Sheets. Nat. Commun. 2014, 5, 3113. 8. Yi, W. -C.; Liu, W.; Botana, J.; Zhao, L.; Liu, Z.; Liu, J. -Y.; Miao, M. -S. Honeycomb Boron Allotropes with Dirac Cones: A True Analogue to Graphene. J. Phys. Chem. Lett. 2017, 8, 2647. 9. Mannix, A. J.; Zhou, X. -F.; Kiraly, B.; Wood, J. D.; Alducin, D.; Myers, B. D.; Liu, X.; Fisher, B. L.; Santiago, U.; Guest, J. R.; Yacaman, M. J.; Ponce, A.; Oganov, A. R.; Hersam, M. C.; Guisinger, N. P.; Synthesis of Borophenes: Anisotropic, Two-Dimensional Boron Polymorphs. Science 2015, 350, 1513-1516. 10. Feng. B.; Zhang, J.; Zhong, Q.; Li, W.; Li, S.; Li, H.; Cheng, P.; Meng, S.; Chen, L.; Wu, K. Experimental Realization of Two-Dimensional Boron Sheets. Nat. Chem. 2016, 8, 563-568. 11. Zhang, Z.; Yang, Y.; Gao, G.; Yakobson, B. I. Two-dimensional Boron Monolayers Mediated by Metal Substrates. Angew. Chem., Int. Ed. 2015, 54, 13022-13026. 12. Zhang, Z.; Mannix, A. J.; Hu, Z.; Kiraly, B.; Guisinger, N. P.; Hersam, M. C.; Yakobson, B. Substrate-induced Nanoscale Undulations of Borophene on Silver. Nano Lett. 2016, 16, 6622-6627. 13. Liu, X.; Wei, Z.; Balla, I.; Mannix, A. J.; Guisinger, N. P.; Luijten, E.; Hersam, M. C. Self-Assembly of Electronically Abrupt Borophene/Organic Lateral Heterostructures. Sci. Adv. 2017, 3, e1602356. 14. Liu, Y.; Penev, E. S.; Yakobson, B. I. Probing the Synthesis of TwoDimensional Boron by First-Principles Computations. Angew. Chem., Int. Ed. 2013, 52, 3156-3159.

15 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

15. Mannix, A. J.; Kiraly, B.; Hersam, M. C.; Guisinger N. P. Synthesis and Chemistry of Elemental 2D Materials. Nat. Rev. Chem. 2017, 1, 0014. 16. Liu, X.; Zhang, Z.; Wang, L.; Yakobson, B. I.; Hersam, M. C. Intermixing and Periodic Self-Assembly of Borophene Line Defects. Nat. Mater. 2018, 17, 783-788. 17. Penev, E. S.; Kutana, A.; Yakobson, B. I. Can Two-dimensional Boron Superconduct? Nano Lett. 2016, 16, 2522-2526. 18. Barth, J. V.; Brune, H.; Ertl, G.; Behm, R. J. Scanning Tunneling Microscopy Observations on the Reconstructed Au(111) Surface: Atomic Structure, LongRange Superstructure, Rotational Domains, and Surface Defects. Phys. Rev. B 1990, 42, 9307-9318. 19. Wöll, Ch.; Chiang, S.; Wilson, R. J.; Lippel, P. H. Determination of Atom Positions at Stacking-Fault Dislocations on Au(111) by Scanning Tunneling Microscopy. Phys. Rev. B 1989, 39, 7988-7991. 20. Meyer, J. A.; Baikie, I. D.; Kopatzki, E.; Behm, R. J. Preferential Island Nucleation at the Elbows of the Au(111) Herringbone Reconstruction Through Place Exchange. Surf. Sci. 1996, 365, L647-L651. 21. Poirier G. E.; Pylant, E. D. The Self-Assembly Mechanism of Alkanethiols on Au(111). Science 1996, 272, 1145-1148. 22. Chambliss, D. D.; Wilson, R. J.; Chiang, S. Nucleation of Ordered Ni Island Arrays on the Au(111) by Surface-Lattice Dislocations. Phys. Rev. Lett. 1991, 66, 1721-1724. 23. Yokoyama, T.; Yokoyama, S.; Kamikado, T.; Okuno, Y.; Mashiko, S. Selective Assembly on a Surface of Supramolecular Aggregates with Controlled Size and Shape. Nature 2001, 413, 619-621. 24. Bartels, L. Tailoring Molecular Layers at Metal Surfaces. Nat. Chem. 2010, 2, 87-95. 25. Voigtländer, B.; Meyer, G.; Amer, N. M. Epitaxial Growth of Thin Magnetic Cobalt Films on Au(111) Studied by Scanning Tunneling Microscopy. Phys. Rev. B 1991, 44, 10354-10357. 26. Brune, H.; Giovannini, M.; Bromann, K.; Kern, K. Self-Organized Growth of Nanostructure Arrays on Strain-Relief Patterns. Nature 1998, 394, 451-453. 27. Brune, H.; Röder, H.; Boragno, C.; Kern, K. Strain Relief at Hexagonal-ClosePacked Interfaces. Phys. Rev. B 1994, 49, 2997-3000. 28. Corso, M.; Fernández, L.; Schiller, F.; Ortega, J. E. Au(111)-Based Nanotemplates by Gd Alloying. ACS Nano 2010, 4, 1603-1611. 29. Bagus, P. S.; Wieckowski, A.; Freund, H. The Contribution of Lattice Strain to Core-Level Binding Energy Shifts in Metal Nanoparticles: Generality and Origin of the Shifts. Comput. Theor. Chem. 2012, 987, 22-24. 30. Richter, B.; Kuhlenbeck, H.; Freund, H.-J.; Bagus, P. S. Cluster Core-Level Binding-Energy Shifts: The Role of Lattice Strain. Phys. Rev. Lett. 2004, 93, 026805. 31. Walter, A. L.; Schiller, F.; Corso, M.; Merte, L. R.; Bertram, F.; Lobo-Checa, J.; Shipilin, M.; Gustafson, J.; Lundgren, E.; Brión-Ríos, A. X.; Cabrera-Sanfelix, P.; Sánchez-Portal, D.; Ortega, J. E. X-ray Photoemission Analysis of Clean and Carbon Monoxide-Chemisorbed Platinum(111) Stepped Surfaces Using a Curved Crystal. Nat. Comm. 2015, 6, 8903. 16 ACS Paragon Plus Environment

Page 16 of 23

Page 17 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

32. Batzill, M. The Surface Science of Graphene: Metal Interfaces, CVD Synthesis, Nanoribbons, Chemical Modifications, and Defects. Surf. Sci. Rep. 2011, 67, 83115. 33. Artyukhov, V. I.; Liu, Y.; Yakobson, B. I. Equilibrium at the Edges and Atomistic Mechanisms of Graphene Growth. Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 1513615140. 34. Bolou, H.; Goyhenex, C. Local Strain Analysis of the Herringbone Reconstruction of Au(111) Through Atomistic Simulations. Phys. Rev. B 2002, 65, 045407. 35. Sergeeva, A. P.; Popov, I. A.; Piazza, Z. A.; Li, W. -L.; Romanescu, C.; Wang, L. -S.; Bolderyv, A. I. Understanding Boron Through Size-Selected Clusters: Structure, Chemical Bonding, and Fluxionality. Acc. Chem. Res. 2014, 47, 13491358. 36. Zhai, H.-J.; Kiran, B.; Li, J.; Wang, L.-S. Hydrocarbon Analogues of Boron Clusters — Planarity, Aromaticity and Antiaromaticity. Nat. Mater. 2003, 2, 827833. 37. Zhai, H.-J.; Zhao, Y. -F.; Li, W. -L.; Chen, Q.; Bai, H.; Hu, H. -S.; Piazza, Z. A.; Tian, W. -J.; Lu, H. -G.; Wu, Y. -B.; Mu, Y. -W.; Wei, G. -F.; Liu, Z. -P.; Li, J. Li, S. -D.; Wang, L. -S.; Observation of an All-boron Fullerene. Nat. Chem. 2014, 6, 727731. 38. Wang, L.-S. Photoelectron Spectroscopy of Size-Selected Boron Clusters: From Planar Structures to Borophenes and Borospherenes. Int. Rev. Phys. Chem. 2016, 35, 69-142. 39. Chen, W.; Madhavan, V.; Jamneala, T.; Crommie, M. F. Scanning Tunneling Microscopy Observation of an Electronic Superlattice at the Surface of Clean Gold. Phys. Rev. Lett. 1998, 80, 1469-1472. 40. Millo, O.; Katz, D.; Cao, Y.-W.; Banin, U. Scanning Tunneling Spectroscopy of InAs Nanocrystal Quantum Dots. Phys. Rev. B 2000, 61, 16773-16777. 41. Banin, U.; Cao, Y.-W.; Katz, D.; Millo, O. Identification of Atomic-like Electronic States in Indium Arsenide Nanocrystal Quantum Dots. Nature 1999, 400, 542-544. 42. Foley, E. T.; Yoder, N. L.; Guisinger, N. P.; Hersam, M. C. Cryogenic Variable Temperature Ultrahigh Vacuum Scanning Tunneling Microscope for Single Molecule Studies on Silicon Surfaces. Rev. Sci. Inst. 2004, 75, 5280-5287. 43. Brockenbrough, R.; Lyding, J. W. Inertial Tip Translator for a Scanning Tunneling Microscope. Rev. Sci. Inst. 1993, 64, 2225-2228. 44. Kresse, G.; Hafner, J. Ab initio Molecular-Dynamics Simulation of the LiquidMetal-Amorphous-Semiconductor Transition in Germanium. Phys. Rev. B 1994, 49 14251. 45. Kresse, G.; Futhermüller, J. Efficient Iterative Schemes for Ab Initio Totalenergy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169. 46. Henkelman, G. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901.

17 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

FIGURE 1

Figure 1. (a) STM image of a clean Au(111) surface that shows the herringbone reconstruction (V = -1.5 V, I = 100 pA). (b) Following boron deposition at 550 °C, the herringbone reconstruction is modified to a trigonal network, where nanoscale borophene islands (one highlighted by a white dashed line) emerge at the nodes resulting in templated growth across the surface (V = -0.1 V, I = 200 pA). (c) The atomic structure of the borophene v1/12 computationally modeled and predicted for Au(111) (left – top down, right – cross-sectional view).

18 ACS Paragon Plus Environment

Page 18 of 23

Page 19 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

FIGURE 2

Figure 2. (a) B1s core-level spectra for room temperature B deposition (orange), trigonal network with low-dose (green), trigonal network with high-dose (red), and clean Au(111) (blue). (b) Increasing boron dose results in the breakdown of the trigonal network and growth of larger borophene islands (top panels: V = -2 V, I = 100 pA; bottom panels: V= 3.4 V, I = 60 pA).

19 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

FIGURE 3

Figure 3: (a) Schematic illustration of borophene growth dynamics. At low substrate temperature, the boron is predominately on the surface and forms boron clusters. For higher temperatures, the boron dissolves into the bulk and then segregates to the surface upon cooling to form 2D borophene sheets. (b) Minimum energy path for boron diffusion on Au(111) (red), penetration from the fcc site on the surface into the subsurface (gray), and diffusion in the subsurface (blue). (c) DFT calculations of the free energies of single boron atoms at various locations on the Au(111) surface and in the Au(111) bulk, as well as boron dimer and trimer in the subsurface. The gold and boron atoms are displayed in yellow and red, respectively. (d) SIMS spectra show that boron is limited to the surface when deposited at room temperature, but is present ~30 nm into the subsurface region in the trigonal network sample that was prepared at high temperature.

20 ACS Paragon Plus Environment

Page 20 of 23

Page 21 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

FIGURE 4

Figure 4: (a) STM image of a larger borophene island that shows the herringbone reconstruction from the underlying Au(111) substrate and atomic-scale periodicity (V = 1 V, I = 100 pA). (b) Atomic-scale periodicity is apparent in larger borophene sheets (V = -0.4 V, I = 80 pA). (c) An FFT of the STM image in (b) reveals that the periodic length scale is ~0.66 nm. (d) STM image following the initial stages of growth of the templated borophene nanostructures (V = -0.4 V, I = 200 pA). (e) Histogram of the island sizes with observable peaks corresponding to integer values of the basic rhombohedral unit.

21 ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

FIGURE 5

Figure 5: (a) Series of bias-dependent STM images showing topography (top panels) and dI/dV maps (bottom panels) of the same region at the biases indicated above each panel. Red and yellow triangles represent fcc and hcp regions of the Au(111) surface, as noted in Fig. 1. The borophene islands appear as depressions in topography at all biases (I = 200 pA). (b) Calculated band structure and electronic density of states (DOS) for monolayer boron on Au(111) showing metallic characteristics. The red and grey lines represent electronic bands contributed by borophene and Au, respectively. The thickness represents the weight of the contribution. (c) dI/dV point spectra taken on (gray) Au(111) for reference, (red) large borophene sheets on Au(111) for comparison, and (dark blue) borophene nano-islands on Au(111). Dashed lines at +/- 300 mV highlight the boron island spectral features that are discussed in the text.

22 ACS Paragon Plus Environment

Page 22 of 23

Page 23 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

TABLE OF CONTENTS FIGURE

23 ACS Paragon Plus Environment