Vacancy Ordering in Ultrathin Copper Oxide Films on Cu(111).

Oxide thin films grown on metal surfaces have wide applications in catalysis and beyond owing to their unique surface structures compared to their bulk counterparts. Despite extensive studies, the atomic structures of copper surface oxides on Cu(111), commonly referred to as "44" and "29", have remained elusive. In this work, we demonstrated an approach for the structural determination of oxide surfaces using element-specific scanning tunneling microscopy (STM) imaging enhanced by functionalized tips. This approach enabled us to resolve the atomic structures of "44" and "29" surface oxides, which were further corroborated by noncontact atomic force microscopy (nc-AFM) measurements and Monte Carlo (MC) simulations. The stoichiometry of the "44" and "29" frameworks was identified as Cu23O16 and Cu16O11, respectively. Contrary to the conventional hypothesis, we observed ordered Cu vacancies within the "44" structure manifesting as peanut-shaped cavities in the hexagonal lattice. Similarly, a combination of Cu and O vacancies within the "29" structure leads to bean-shaped cavities within the pentagonal lattice. Our study has thus resolved the decades-long controversy on the atomic structures of "44" and "29" surface oxides, advancing our understanding of copper oxidation processes and introducing a robust framework for the analysis of complex oxide surfaces.

Atomic-scale identification of nitrogen dopants in graphene on Ir(111) and Ru(0001).

Nitrogen (N) doped graphene materials have been synthesized using the sole precursor adenine on the Ir(111) and Ru(0001) surfaces. X-ray photoelectron spectroscopy and scanning tunneling microscopy (STM) have been used to characterize the obtained N-doped graphene materials. Several graphitic and pyridinic N dopants have been identified on the atomic scale by combining STM measurements and STM simulations based on density functional theory calculations.

Changing the Interaction of a Single-Molecule Magnetic Moment with a Superconductor.

The exchange interaction of a brominated Co-porphyrin molecule with the Cooper pair condensate of Pb(111) is modified by reducing the Co-surface separation. The stepwise dehalogenation and dephenylation change the Co adsorption height by a few picometers. Only the residual Co-porphine core exhibits a Yu-Shiba-Rusinov bound state with low binding energy in the Bardeen-Cooper-Schrieffer energy gap. Accompanying density functional calculations reveal that the Co dz2 orbital carries the molecular magnetic moment and is responsible for the intragap state. The calculated spatial evolution of the Yu-Shiba-Rusinov wave function is compatible with the experimentally observed oscillatory attenuation of the electron-hole asymmetry with increasing lateral distance from the magnetic porphine center.

Tuning Ferromagnetism in a Single Layer of Fe above Room Temperature.

The crystallographic and magnetic properties of an Fe monolayer (ML) grown on 2 ML Au/W(110) substrate are studied with spin-polarized low-energy electron microscopy, density functional theory, and relativistic screened Korringa-Kohn-Rostoker calculations. The single layer of iron atoms possesses hexagonal symmetry and reveals a ferromagnetic order at room temperature. We experimentally demonstrate the possibility of tuning the Curie temperature and the magnitude of magnetization of the Fe monolayer by capping with Au. Taking into account several structural models, the calculation results mostly show ferromagnetic states with enhanced magnetic moments of Fe atoms compared to their bulk value and a further increase in their value after covering with Au. The theoretically calculated Curie temperatures are in fair agreement with those obtained in the experiments. The calculations, furthermore, found evidence for the presence of frustrated isotropic Fe-Fe exchange interactions, and a discussion of the structural effects on the magnetic properties is provided herein.

On-surface synthesis of 2D COFs via molecular assembly directed photocycloadditions: a first-principles investigation.

Over the past decades, the rational synthesis of two-dimensional covalent organic framework (2D COFs) monolayer via on-surface chemistry has been widely explored. Herein, we propose the [2 + 2] photocycloaddition as a novel strategy for large-scale fabrication of COFs from theoretical perspective. Thanks to the symmetry forbidden of thermal [2 + 2] cycloaddition, the molecular precursors carrying vinyl groups will not chemically interact with each other during thermal annealing, which is essential to achieve molecular assembly. The subsequent photocycloaddition of these precursors may produce large-scale 2D COFs at low temperatures, in which the symmetry of molecular assembly remains unchanged. Our results show that 2D COFs can be produced via [2 + 2] photocycloadditions directed from self-assembled precursors, in which alkylbenzene molecules with vinyl groups on side chains exhibit appropriate intermolecular distances. By performing high-throughput calculations, several promising molecular precursors are proposed to achieve large-scale 2D COFs. This work provides an applicable strategy for the large-scale synthesis of 2D carbon materials.

Electronic and Magnetic Properties of Building Blocks of Mn and Fe Atomic Chains on Nb(110).

We present results for the electronic and magnetic structure of Mn and Fe clusters on Nb(110) surface, focusing on building blocks of atomic chains as possible realizations of topological superconductivity. The magnetic ground states of the atomic dimers and most of the monatomic chains are determined by the nearest-neighbor isotropic interaction. To gain physical insight, the dependence on the crystallographic direction as well as on the atomic coordination number is analyzed via an orbital decomposition of this isotropic interaction based on the spin-cluster expansion and the difference in the local density of states between ferromagnetic and antiferromagnetic configurations. A spin-spiral ground state is obtained for Fe chains along the [11¯0] direction as a consequence of the frustration of the isotropic interactions. Here, a flat spin-spiral dispersion relation is identified, which can stabilize spin spirals with various wave vectors together with the magnetic anisotropy. This may lead to the observation of spin spirals of different wave vectors and chiralities in longer chains instead of a unique ground state.

Spin-orbit coupling induced splitting of Yu-Shiba-Rusinov states in antiferromagnetic dimers.

Magnetic atoms coupled to the Cooper pairs of a superconductor induce Yu-Shiba-Rusinov states (in short Shiba states). In the presence of sufficiently strong spin-orbit coupling, the bands formed by hybridization of the Shiba states in ensembles of such atoms can support low-dimensional topological superconductivity with Majorana bound states localized on the ensembles' edges. Yet, the role of spin-orbit coupling for the hybridization of Shiba states in dimers of magnetic atoms, the building blocks for such systems, is largely unexplored. Here, we reveal the evolution of hybridized multi-orbital Shiba states from a single Mn adatom to artificially constructed ferromagnetically and antiferromagnetically coupled Mn dimers placed on a Nb(110) surface. Upon dimer formation, the atomic Shiba orbitals split for both types of magnetic alignment. Our theoretical calculations attribute the unexpected splitting in antiferromagnetic dimers to spin-orbit coupling and broken inversion symmetry at the surface. Our observations point out the relevance of previously unconsidered factors on the formation of Shiba bands and their topological classification.

Polymorphic expressions of ultrathin oxidic layers of Mo on Au(111).

Ultrathin oxidic layers of Mo (i.e. O/Mo) on the Au(111) substrate are investigated using first-principles density-functional theory calculations. Various polymorphic structural models of these O/Mo layers are proposed and compared with previous experimental results - covering both spectroscopic and microscopic approaches of characterization. We find that, through the control of metal-oxygen coordination in these ultrathin oxidic O/Mo films on Au(111), the oxidation state of Mo atoms in the O/Mo layers can be modulated and reduced without intentional creation of oxygen vacancies. This is also assisted by a charge transfer mechanism from the Au substrate to these oxidic films, providing a direct means to tune the surface electronic properties of ultrathinoxide films on metal substrates.

Coordination-Controlled C-C Coupling Products via ortho-Site C-H Activation.

The coordination-restricted ortho-site C-H bond activation and dehydrogenative homocoupling of 4,4'-(1,3-phenylene)dipyridine (1,3-BPyB) and 4,4'-(1,4-phenylene)dipyridine (1,4-BPyB) on different metal surfaces were studied by a combination of scanning tunneling microscopy, noncontact atomic force microscopy, and density functional theory calculations. The coupling products on Cu(111) exhibited certain configurations subject to the spatial restriction of robust two-fold Cu-N coordination bonds. Compared to the V-shaped 1,3-BPyB, the straight backbone of 1,4-BPyB helped to further reduce the variety of reactive products. By utilizing the three-fold coordination of Fe atoms with 1,4-BPyB molecules on Au(111), a large-scale network containing single products was constructed. Our results offer a promising protocol for controllable on-surface synthesis with the aid of robust coordination interactions.

Atomic and electronic structure of the Si(331)-(12 × 1) surface.

We report on the investigation of the atomic and electronic structures of a clean Si(331)-(12 × 1) surface using a first-principles approach with both plane wave and strictly localized basis sets. Starting from the surface structure proposed by Zhachuk and Teys [Phys. Rev. B 95, 041412(R) (2017)], we develop significant improvements to the atomic model and localized basis set which are critical for the correct description of the observed bias dependence of scanning tunneling microscopy (STM) images. The size mismatch between the Si pentamers from the surface model and those seen by STM is explained within the context of the Tersoff-Hamann model. The energy barriers that separate different Si(331) buckled configurations were estimated, showing that the surface structure is prone to dynamic buckling at room temperature. It is found that empty electronic states on Si(331) are essentially localized on the pentamers with interstitials and under-coordinated Si sp 2-like atoms between them, while filled electronic states are localized on under-coordinated Si sp 3-like atoms and dimers on trenches. The calculated electronic density of states exhibits two broad peaks in the fundamental bandgap of Si: one near the valence band top and the other near the conduction band bottom. The resulting surface bandgap of 0.58 eV is in an excellent agreement with spectroscopy studies.

On the Structure of Ultrathin FeO Films on Ag(111).

Ultrathin transition metal oxide films exhibit unique physical and chemical properties not observed for the corresponding bulk oxides. These properties, originating mainly from the limited thickness and the interaction with the support, make those films similar to other supported 2D materials with bulk counterparts, such as transition metal dichalcogenides. Ultrathin iron oxide (FeO) films, for example, were shown to exhibit unique electronic, catalytic and magnetic properties that depend on the type of the used support. Ag(111) has always been considered a promising substrate for FeO growth, as it has the same surface symmetry, only ~5% lattice mismatch, is considered to be weakly-interacting and relatively resistant to oxidation. The reports on the growth and structure of ultrathin FeO films on Ag(111) are scarce and often contradictory to each other. We attempted to shed more light on this system by growing the films using different preparation procedures and studying their structure using scanning tunneling microscopy (STM), low energy electron diffraction (LEED) and X-ray photoelectron spectroscopy (XPS). We observed the formation of a previously unreported Moiré superstructure with 45 Å periodicity, as well as other reconstructed and reconstruction-free surface species. The experimental results obtained by us and other groups indicate that the structure of FeO films on this particular support critically depends on the films' preparation conditions. We also performed density functional theory (DFT) calculations on the structure and properties of a conceptual reconstruction-free FeO film on Ag(111). The results indicate that such a film, if successfully grown, should exhibit tunable thickness-dependent properties, being substrate-influenced in the monolayer regime and free-standing-FeO-like when in the bilayer form.

Toward tailoring Majorana bound states in artificially constructed magnetic atom chains on elemental superconductors.

Realizing Majorana bound states (MBS) in condensed matter systems is a key challenge on the way toward topological quantum computing. As a promising platform, one-dimensional magnetic chains on conventional superconductors were theoretically predicted to host MBS at the chain ends. We demonstrate a novel approach to design of model-type atomic-scale systems for studying MBS using single-atom manipulation techniques. Our artificially constructed atomic Fe chains on a Re surface exhibit spin spiral states and a remarkable enhancement of the local density of states at zero energy being strongly localized at the chain ends. Moreover, the zero-energy modes at the chain ends are shown to emerge and become stabilized with increasing chain length. Tight-binding model calculations based on parameters obtained from ab initio calculations corroborate that the system resides in the topological phase. Our work opens new pathways to design MBS in atomic-scale hybrid structures as a basis for fault-tolerant topological quantum computing.

Inducing skyrmions in ultrathin Fe films by hydrogen exposure.

Magnetic skyrmions are localized nanometer-sized spin configurations with particle-like properties, which are envisioned to be used as bits in next-generation information technology. An essential step toward future skyrmion-based applications is to engineer key magnetic parameters for developing and stabilizing individual magnetic skyrmions. Here we demonstrate the tuning of the non-collinear magnetic state of an Fe double layer on an Ir(111) substrate by loading the sample with atomic hydrogen. By using spin-polarized scanning tunneling microscopy, we discover that the hydrogenated system supports the formation of skyrmions in external magnetic fields, while the pristine Fe double layer does not. Based on ab initio calculations, we attribute this effect to the tuning of the Heisenberg exchange and the Dzyaloshinsky-Moriya interactions due to hydrogenation. In addition to interface engineering, hydrogenation of thin magnetic films offers a unique pathway to design and optimize the skyrmionic states in low-dimensional magnetic materials.

Atomistic Origins of Surface Defects in CH3NH3PbBr3 Perovskite and Their Electronic Structures.

The inherent instability of CH3NH3PbX3 remains a major technical barrier for the industrial applications of perovskite materials. Recently, the most stable surface structures of CH3NH3PbX3 have been successfully characterized by using density functional theory (DFT) calculations together with the high-resolution scanning tunneling microscopy (STM) results. The two coexisting phases of the perovskite surfaces have been ascribed to the alternate orientation of the methylammonium (MA) cations. Notably, similar surface defect images (a dark depression at the sites of X atoms) have been observed on surfaces produced with various experimental methods. As such, these defects are expected to be intrinsic to the perovskite crystals and may play an important role in the structural decomposition of perovskite materials. Understanding the nature of such defects should provide some useful information toward understanding the instability of perovskite materials. Thus, we investigate the chemical identity of the surface defects systematically with first-principles density functional theory calculations and STM simulations. The calculated STM images of the Br and Br-MA vacancies are both in good agreement with the experimental measurements. In vacuum conditions, the formation energy of Br-MA is 0.43 eV less than the Br vacancy. In the presence of solvation effects, however, the formation energy of a Br vacancy becomes 0.42 eV lower than the Br-MA vacancy. In addition, at the vacancy sites, the adsorption energies of water, oxygen, and acetonitrile molecules are significantly higher than those on the pristine surfaces. This clearly demonstrated that the structural decomposition of perovskites are much easier to start from these vacancy sites than the pristine surfaces. Combining DFT calculations and STM simulations, this work reveals the chemical identities of the intrinsic defects in the CH3NH3PbX3 perovskite crystals and their effects on the stability of perovskite materials.

Nanostructural adsorption of vanadium oxide on functionalized graphene: a DFT study.

Rutile-monoclinic phase transitions of vanadium oxide (VO2) nanocrystals adsorbed on graphene-based substrates are of current scientific interest, although their adsorption and growth mechanisms have not been investigated theoretically. In this study, we use density functional theory (DFT) calculations for determining the binding energies and predicting the corresponding directions of growth of VO2 nanostructures (rutile and M1-monoclinic) interacting with both pure graphene and functionalized graphene nanoribbons. Several adsorption sites of pure graphene including the top, bridge, and hollow sites are considered, while additional adsorption sites of functionalized graphene nanoribbons, epoxy, alcohol and carboxylate are investigated. Vanadium oxide nanostructures are found to favor physical adsorption on the hollow sites of pure graphene, while chemical adsorption is favored on the carboxylate sites of functionalized graphene nanoribbons (FGNRs). Charge density maps showed the electron distribution originating from the interaction between VO2 and graphene substrates, helping to understand the mechanism of charge transfer. Electronic local potentials showed vertical growth tendencies for rutile VO2, while M1-monoclinic VO2 showed horizontal growth tendencies. Partial density of states (PDOS) helped examine the electronic structure of metallic rutile VO2 binding to hollow and carboxylate sites of functionalized graphene. These results provide an improved understanding of the controlled and oriented growth of VO2 nanocrystals on graphene-based substrates which can enable various properties such as the metal-insulator transition (MIT) of VO2 in light regulation applications.

Microscopic origin of chiral shape induction in achiral crystals.

In biomineralization, inorganic materials are formed with remarkable control of the shape and morphology. Chirality, as present in the biomolecular world, is therefore also common for biominerals. Biomacromolecules, like proteins and polysaccharides, are in direct contact with the mineral phase and act as modifiers during nucleation and crystal growth. Owing to their homochirality--they exist only as one of two possible mirror-symmetric isomers--their handedness is often transferred into the macroscopic shape of the biomineral crystals, but the way in which handedness is transmitted into achiral materials is not yet understood at the atomic level. By using the submolecular resolution capability of scanning tunnelling microscopy, supported by photoelectron diffraction and density functional theory, we show how the chiral 'buckybowl' hemibuckminsterfullerene arranges copper surface atoms in its vicinity into a chiral morphology. We anticipate that such new insight will find its way into materials synthesis techniques.

Contrast stability and 'stripe' formation in scanning tunnelling microscopy imaging of highly oriented pyrolytic graphite: the role of STM-tip orientations.

Highly oriented pyrolytic graphite (HOPG) is an important substrate in many technological applications and is routinely used as a standard in Scanning Tunnelling Microscopy (STM) calibration, which makes the accurate interpretation of the HOPG STM contrast of great fundamental and applicative importance. We demonstrate by STM simulations based on electronic structure obtained from first principles that the relative local orientation of the STM-tip apex with respect to the HOPG substrate has a considerable effect on the HOPG STM contrast. Importantly for experimental STM analysis of HOPG, the simulations indicate that local tip-rotations maintaining a major contribution of the d(3z(2)-r(2)) tip-apex state to the STM current affect only the secondary features of the HOPG STM contrast resulting in 'stripe' formation and leaving the primary contrast unaltered. Conversely, tip-rotations leading to enhanced contributions from m ≠ 0 tip-apex electronic states can cause a triangular-hexagonal change in the primary contrast. We also report a comparison of two STM simulation models with experiments in terms of bias-voltage-dependent STM topography brightness correlations and discuss our findings for the HOPG(0 0 0 1) surface in combination with tungsten tip models of different sharpnesses and terminations.

Arbitrary tip orientation in STM simulations: 3D WKB theory and application to W(110).

We extend the orbital-dependent electron tunnelling model implemented within the three-dimensional (3D) Wentzel-Kramers-Brillouin (WKB) atom-superposition approach for simulating scanning tunnelling microscopy (STM) by including arbitrary tip orientations. The orientation of the tip is characterized by a local coordinate system centred on the tip apex atom obtained by a rotation with respect to the sample coordinate system. The rotation is described by the Euler angles. Applying our method, we highlight the role of the real-space shape of the electron orbitals involved in the tunnelling, and analyse the convergence and the orbital contributions of the tunnelling current above the W(110) surface depending on the orientation of a model tungsten tip. We also simulate STM images at constant-current condition, and find that their quality depends very much on the tip orientation. Some orientations result in protrusions on the images that do not occur above W atoms. The presence of such apparent atom positions makes it difficult to identify the exact position of surface atoms. It is suggested that this tip orientation effect should be considered in the evaluation of experimental STM images on other surfaces as well. The presented computationally efficient tunnelling model could prove to be useful for obtaining more information on the local tip geometry and orientation by comparing STM experiments to a large number of simulations with systematically varied tip orientations.

Role of hydrogen in giant spin polarization observed on magnetic nanostructures.

We demonstrate that the giant spin contrast observed by scanning tunneling microscopy for double-layer Coislands on Pt(111) is caused by adsorbates at the apex of the Cr-coated W tip. The most likely candidate, in ab initio simulations, is hydrogen. Here, the electron charge is highly polarized by the adjacent Cr layers. The hydrogen adsorption site is shown to change from hollow to on top due to the electric field at the tip apex, created by the tunnel voltage.

Site- and orientation-selective anchoring of a prototypical molecular building block.

The controlled anchoring of molecular building blocks on appropriate templates is a major prerequisite for the rational design and fabrication of supramolecular architectures on surfaces. We report on a particularly selective adsorption process of hexa-peri-hexabenzocoronene on Au(111), which leads to well-controlled adsorption position and orientation of the polycyclic aromatic hydrocarbons. Scanning tunneling microscopy reveals selective adsorption on monatomic steps in the fcc stacking regions with a specific orientation of 18 degrees between the molecular axis and the step normal. Ab initio calculations for various adsorption sites reveal the lowest total energy for adsorption on a kink site. Energy considerations and the excellent agreement between experimental and simulated images show that adsorption on kink sites is responsible for the specific adsorption angle.