Search

Graphyne (GY) and graphdiyne (GDY)-based materials represent an intriguing class of two-dimensional (2D) carbon-rich networks with tunable structures and properties surpassing those of graphene. However, the challenge of fabricating atomically well-defined crystalline GY/GDY-based systems largely hinders detailed electronic structure characterizations. Here, we report the emergence of an unconventional band structure in mesoscopically regular (~1 {\mu}m) metallated GDY sheets featuring a honeycomb lattice on Ag(111) substrates. Employing complementary scanning tunnelling and angle-resolved photoemission spectroscopies, electronic band formation with a gap of 2.5 eV is rigorously determined in agreement with real-space electronic characteristics. Extensive density functional theory calculations corroborate our observations as well as recent theoretical predictions that doubly degenerate frontier molecular orbitals on a honeycomb lattice give rise to flat, Dirac and Kagome bands close to Fermi level. These results illustrate the tremendous potential of engineering novel band structures via molecular orbital and lattice symmetries in atomically precise 2D carbon scaffolds.

The production of high-quality hexagonal boron nitride (h-BN) is essential for the ultimate performance of two-dimensional (2D) materials-based devices, since it is the key 2D encapsulation material. Here, a decisive guideline is reported for fabricating high-quality h-BN on transition metals: It is crucial to exclude carbon from h-BN related process. Otherwise carbon prevails over boron and nitrogen due to its larger binding energy, thereupon forming graphene on metals after high-temperature annealing. We demonstrate the surface reaction-assisted conversion from h-BN to graphene with high-temperature treatments. The pyrolysis temperature Tp is an important quality indicator for h-BN/metals. When the temperature is lower than Tp, the quality of h-BN layer is improved upon annealing. While the annealing temperature is above Tp, in case of carbon-free conditions, the h-BN disintegrates and nitrogen desorbs from the surface more easily than boron, eventually leading to clean metal surfaces. However, once the h-BN layer is exposed to carbon, graphene forms on Pt(111) in the high-temperature regime. This not only provides an indispensable principle (avoid carbon) for fabricating high-quality h-BN materials on transition metals, but also offers a straightforward method for the surface reaction-assisted conversion from h-BN to graphene on Pt(111)., Comment: Accepted by SMALL on October 7, 2022

Ionic liquids (IL) are promising electrolytes for electrochemical applications due to their remarkable stability and high charge density. Molecular dynamics simulations are essential for better understanding the complex phenomena occurring at the electrode-IL interface. In this work, we have studied the interface between graphene and 1-ethyl-3-methyl-imidazolium tetrafluoroborate IL, using density functional theory-based molecular dynamics simulations at variable surface charge densities. We have disassembled the electrical double layer potential drop into two main components: one involving atomic charges and the other - dipoles. The latter component arises due to the electronic polarisation of the surface and is related to concepts hotly debated in the literature, such as the Thomas-Fermi screening length, effective surface charge plane, and quantum capacitance.

Phosphorene is a new two-dimensional material composed of a single or few atomic layers of black phosphorus. Phosphorene has both an intrinsic tunable direct band gap and high carrier mobility values, which make it suitable for a large variety of optical and electronic devices. However, the synthesis of single-layer phosphorene is a major challenge. The standard procedure to obtain phosphorene is by exfoliation. More recently, the epitaxial growth of single-layer phosphorene on Au(111) has been investigated by molecular beam epitaxy and the obtained structure has been described as a blue-phosphorene sheet. In the present study, large areas of high-quality monolayer phosphorene, with a band gap value at least equal to 0.8 eV, have been synthesized on Au(111). Our experimental investigations, coupled with DFT calculations, give evidence of two distinct phases of blue phosphorene on Au(111), instead of one as previously reported, and their atomic structures have been determined., Comment: This paper reports on the epitaxial synthesis of blue phosphorene

The adsorption of molecules on surfaces affects the surface dipole and thus changes in the work function may be expected. The effect in change of work function is particularly strong if charge between substrate and adsorbate is involved. Here we report the deposition of a strong electron acceptor molecule, tetrafluorotetracyanoquinodimethane C$_{12}$F$_4$N$_4$ (F$_{4}$TCNQ) on a monolayer of hexagonal boron nitride nanomesh ($h$-BN on Rh(111)). The work function of the F$_{4}$TCNQ/$h$-BN/Rh system increases upon increasing molecular coverage. The magnitude of the effect indicates electron transfer from the substrate to the F$_{4}$TCNQ molecules. Density functional theory calculations confirm the work function shift and predict doubly charged F$_{4}$TCNQ$^{2-}$ in the nanomesh pores, where the $h$-BN is closest to the Rh substrate, and to have the largest binding energy there. The preferred adsorption in the pores is conjectured from a series of ultraviolet photoelectron spectroscopy data, where the $\sigma$ bands in the pores are first attenuated. Scanning tunneling microscopy measurements indicate that F$_{4}$TCNQ molecules on the nanomesh are mobile at room temperature, as "hopping" between neighboring pores is observed.

There is enormous recent interest in weak, van der Waals-type (vdW) interactions due to their fundamental relevance for two-dimensional materials and the so-called vdW heterostructures. Tackling this problem using computer simulation is very challenging due to the non-trivial, non-local nature of these interactions. We benchmark different treatments of London dispersion forces within the density functional theory (DFT) framework on hexagonal boron nitride or graphene monolayers on Ir(111) by comparing the calculated geometries to a comprehensive set of experimental data. The geometry of these systems crucially depends on the interplay between vdW interactions and wave function hybridisation, making them excellent test cases for vdW-treated DFT. Our results show strong variations in the calculated atomic geometry. While some of the approximations reproduce the experimental structure, this is rather based on \textit{a posteriori} comparison with the ``target results''. General predictive power in vdW-treated DFT is not achieved yet and might require new approaches., Comment: More or less published version, with Supporting Material

Hexagonal boron nitride (h-BN) is a prominent member in the growing family of two-dimensional materials with potential applications ranging from being an atomically smooth support for other 2D materials to templating growth of molecular layers. We have studied the structure of monolayer h-BN grown by chemical vapour deposition on Ir(111) by low-temperature scanning tunneling microscopy (STM) and spectroscopy (STS) experiments and state-of-the-art density functional theory (DFT) calculations. The lattice-mismatch between the h-BN and Ir(111) surface results in the formation of a moir\'e superstructure with a periodicity of $\sim$29 \AA\ and a corrugation of $\sim$0.4 \AA. By measuring the field emission resonances above the h-BN layer, we find a modulation of the work function within the moir\'e unit cell of $\sim$0.5 eV. DFT simulations for a 13-on-12 h-BN/Ir(111) unit cell confirm our experimental findings and allow us to relate the change in the work function to the subtle changes in the interaction between boron and nitrogen atoms and the underlying substrate atoms within the moir\'e unit cell. Hexagonal boron nitride on Ir(111) combines weak topographic corrugation with a strong work function modulation over the moir\'e unit cell. This makes h-BN/Ir(111) a potential substrate for electronically modulated thin film and hetero-sandwich structures., Comment: 8 pages, 5 figures

An interaction potential including chloride anion polarization effects, constructed from first-principles calculations, is used to examine the structure and transport properties of a series of chloroaluminate melts. A particular emphasis was given to the study of the equimolar mixture of aluminium chloride with 1-ethyl-3-methylimidazolium chloride, which forms a room temperature ionic liquid EMI-AlCl 4. The structure yielded by the classical simulations performed within the framework of the polarizable ion model is compared to the results obtained from entirely electronic structure-based simulations: An excellent agreement between the two flavors of molecular dynamics is observed. When changing the organic cation EMI+ by an inorganic cation with a smaller ionic radius (Li+, Na+, K+), the chloroaluminate speciation becomes more complex, with the formation of Al2Cl 7- in small amounts. The calculated transport properties (diffusion coefficients, electrical conductivity and viscosity) of EMI-AlCl4 are in good agreement with experimental data.

We report on experimental evidence for the formation of a two dimensional Si/Au(110) surface alloy. In this study, we have used a combination of scanning tunneling microscopy, low energy electron diffraction, Auger electron spectroscopy and ab initio calculations based on density functional theory. A highly ordered and stable Si-Au surface alloy is observed subsequent to growth of a sub-monolayer of silicon on an Au(110) substrate kept above the eutectic temperature.

Understanding the conditions which favor crystallisation or vitrification of liquids has been a long-standing scientific problem. Another connected, and not yet well understood question is the relationship between the glassy and the various possible crystalline forms a system may adopt. In this context, B2O3 is a puzzling case of study since i) it is one of the best glass-forming systems despite an apparent lack of low-pressure polymorphism ii) it vitrifies in a glassy form abnormally different from the only known crystalline phase at ambient pressure iii) it never crystallises from the melt unless pressure is applied, an intriguing behaviour known as the crystallisation anomaly. Here, by means of ab-initio calculations, we discover the existence of novel B2O3 crystalline polymorphs with structural properties similar to the glass and formation energies comparable to the known ambient crystal. The resulting configurational degeneracy drives the system vitrification at ambient pressure. The degeneracy is lifted under pressure, unveiling the origin of the crystallisation anomaly. This work reconciles the behaviour of B2O3 with that from other glassy systems and reaffirms the role played by polymorphism in a system's ability to vitrify. Some of the predicted crystals are cage-like materials entirely made of three-fold rings, opening new perspectives for the synthesis of boron-based nanoporous materials., Comment: 4 pages, 3 figures. This is the initial version as submitted to Nature Materials. Final version contains supplementary information. See doi:10.1038/nmat3416

We determine, by means of density functional theory, the stability and the structure of graphene nanoribbon (GNR) edges in presence of molecules such as oxygen, water, ammonia, and carbon dioxide. As in the case of hydrogen-terminated nanoribbons, we find that the most stable armchair and zigzag configurations are characterized by a non-metallic/non-magnetic nature, and are compatible with Clar's sextet rules, well known in organic chemistry. In particular, we predict that, at thermodynamic equilibrium, neutral GNRs in oxygen-rich atmosphere should preferentially be along the armchair direction, while water-saturated GNRs should present zigzag edges. Our results promise to be particularly useful to GNRs synthesis, since the most recent and advanced experimental routes are most effective in water and/or ammonia-containing solutions., Comment: accepted for publication in PRB

We show that Clar's theory of the aromatic sextet is a simple and powerful tool to predict the stability, the \pi-electron distribution, the geometry, the electronic/magnetic structure of graphene nanoribbons with different hydrogen edge terminations. We use density functional theory to obtain the equilibrium atomic positions, simulated scanning tunneling microscopy (STM) images, edge energies, band gaps, and edge-induced strains of graphene ribbons that we analyze in terms of Clar formulas. Based on their Clar representation, we propose a classification scheme for graphene ribbons that groups configurations with similar bond length alternations, STM patterns, and Raman spectra. Our simulations show how STM images and Raman spectra can be used to identify the type of edge termination.

We develop a first-principles scheme based on the continued fraction approach an d ultrasoft pseudopotentials to calculate K-edge X-ray absorption spectra in solids. The method allows for calculations of K-edge X-ray absorption spectra in transition metal and rare-earths compounds with substantially reduced cutoffs respect to the norm-conserving case. We validate the method by calculating Si and O K-edges in $\alpha$ quartz, Cu K-edge in Copper and in La$_2$CuO$_4$. For the case of Si and O edges in $\alpha$ quartz and in Copper we obtain a good agreement with experimental data. In the Cu K-edge spectra of La$_2$CuO$_4$, a material considered a real challeng e for density functional theory we attribute all the near-edge and far-edge peaks to single particle excitations., Comment: minor corrections, some references added, 9 pages, 5 figures, to appear on Phys. Rev. B

We compute the orbital magnetization in real materials by evaluating a recently discovered formula for periodic systems, within density functional theory. We obtain improved values of the orbital magnetization in the ferromagnetic metals Fe, Co, and Ni, by taking into account the contribution of the interstitial regions neglected so far in literature. We also use the orbital magnetization to compute the EPR $g$-tensor in molecules and solids. The present approach reproduces the $g$-tensor obtained by linear response (LR), when the spin-orbit can be treated as a perturbation. However, it can also be applied to radicals and defects with an orbital-degenerate ground-state or containing heavy atoms, that can not be properly described by LR., Comment: 6 pages, 6 tables; two new tables added in the auxiliary material

We determine the stability, the geometry, the electronic and magnetic structure of hydrogen-terminated graphene-nanoribbons edges as a function of the hydrogen content of the environment by means of density functional theory. Antiferromagnetic zigzag ribbons are stable only at extremely-low ultra-vacuum pressures. Under more standard conditions, the most stable structures are the mono- and di-hydrogenated armchair edges and a zigzag edge reconstruction with one di- and two mono-hydrogenated sites. At high hydrogen-concentration ``bulk'' graphene is not stable and spontaneously breaks to form ribbons, in analogy to the spontaneous breaking of graphene into small-width nanoribbons observed experimentally in solution. The stability and the existence of exotic edge electronic-states and/or magnetism is rationalized in terms of simple concepts from organic chemistry (Clar's rule), Comment: 4 pages, 3 figures, accepted for publication by Physical Review Letters