Abstract: When using hexagonal boron-nitride (hBN) as a substrate for graphene, the resulting moire pattern creates secondary Dirac points. By encapsulating a multilayer graphene within aligned hBN sheets the controlled moire stacking may offer even richer benefits. Using advanced tight-binding simulations on atomistically-relaxed heterostructures, here we show that the gap at the secondary Dirac point can be opened in selected moire-stacking configurations, and is independent of any additional vertical gating of the heterostructure. On the other hand, gating can broadly tune the gap at the principal Dirac point, and may thereby strongly compress the first moire mini-band in width against the moire-induced gap at the secondary Dirac point. We reveal that in hBN-encapsulated bilayer graphene this novel mechanism can lead to isolated bands flatter than 10 meV under moderate gating, hence presenting a convenient pathway towards electronically-controlled strongly-correlated states on demand.

Abstract: Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. In this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

Abstract: Large efforts have been made trying to understand the origin of the high catalytic activity of dealloyed nanoporous gold as a green catalyst for the selective promotion of chemical reactions at low temperatures. Residual silver, left in the sample after dealloying of a gold-silver alloy, has been shown to have a strong influence on the activity of the catalyst. But the question of how the silver is distributed within the porous structure has not finally been answered yet. We show by quantitative energy dispersive X-ray tomography measurements that silver forms clusters that are distributed irregularly, both on the surface and inside the ligaments building up the porous structure. Furthermore, we find that the role of the residual silver is ambiguous. Whereas CO oxidation is supported by more residual silver, methanol oxidation to methyl formate is hindered. Structural characterisation reveals larger ligaments and pores for decreasing residual silver concentration.

Abstract: MFI zeolite nanosheets with thickness of 2 and 8 nm were synthesized, transformed into bifunctional catalysts by loading with platinum and tested in n-decane isomerization and hydrocracking. Detailed analysis of skeletal isomers and hydrocracked products revealed that the MFI nanosheets display transition-state shape-selectivity similar to bulk MFI zeolite crystals. The suppressed formation of bulky skeletal isomers and C5 cracking products are observed both in the nanosheets and the bulk crystals grown in three dimensions. This is typical for restricted transition-state shape-selectivity, characteristic for the MFI type pores. It is a first clear example of transition-state shape-selectivity inside a zeolitic nanosheet. Owing to the short diffusion path across the sheets, expression of diffusion-based discrimination of reaction products in the MFI nanosheets was limited. The 2-methylnonane formation among monobranched C10 isomers and 2,7-dimethyloctane among dibranched C10 isomers, which in MFI zeolite are favored by product diffusion, was much less favored on the nanosheets compared to the reference bulk ZSM-5 material.

Abstract: Black phosphorus (BP) has recently triggered an unprecedented interest in the 2D community. However, many of its unique properties are not exploited and the well-known environmental vulnerability is not conquered. Herein, a type-I mixed-dimensional (0D-1D) van der Waals heterojunction is developed, where three-atomic-layer BP quantum dots (QDs) are assembled on a single ZnO nanorod (NR). By adjusting the indium (In) content in ZnO NRs, the degree and even the direction of surface charge transfer doping within the heterojunction can be tuned, which result in selective Raman scattering enhancements between ZnO and BP. The maximal enhancement factor is determined as 4340 for BP QDs with sub-ppm level. Furthermore, an unexpected long-term ambient stability (more than six months) of BP QDs is revealed, which is ascribed to the electron doping from ZnO:In NRs. The first demonstration of selective Raman enhancements between two inorganic semiconductors as well as the improved stability of BP shed light on this emerging 2D material.

Abstract: We present extensive calculations of the optical and plasmonic properties of a graphene sheet carrying a dc current. By calculating analytically the density-density response function of current-carrying states at finite temperature, we demonstrate that an applied dc current modifies the Pauli blocking mechanism and that absorption acquires a birefringent character with respect to the angle between the in-plane light polarization and current flow. Employing the random phase approximation at finite temperature, we show that graphene plasmons display a degree of non-reciprocity and collimation that can be tuned with the applied current. We discuss the possibility to measure these effects.

Abstract: Since the advent of graphene, other 2D materials have garnered interest; notably the single element materials silicene, germanene, and stanene. Weinvestigate the ballistic current-voltage (I-V) characteristics of armchair silicene and stanene armchair nanoribbons (AXNRs with X = Si, Sn) using a combination of density functional theory and non-equilibrium Green's functions. The impact of out-of-plane electric field and in-plane uniaxial strain on the ribbon geometries, electronic structure, and (I-V)s are considered and contrasted with graphene. Since silicene and stanene are sp(2)/sp(3) buckled layers, the electronic structure can be tuned by an electric field that breaks the sublattice symmetry, an effect absent in graphene. This decreases the current by similar to 50% for Sn, since it has the largest buckling. Uniaxial straining of the ballistic channel affects the AXNR electronic structure in multiple ways: it changes the bandgap and associated effective carrier mass, and creates a local buckling distortion at the lead-channel interface which induces a interface dipole. Due to the increasing sp(3) hybridization character with increasing element mass, large reconstructions rectify the strained systems, an effect absent in sp(2) bonded graphene. This results in a smaller strain effect on the current: a decrease of 20% for Sn at 15% tensile strain compared to a similar to 75% decrease for C.

Abstract: <script type='text/javascript'>document.write(unpmarked('Quantum electrodynamics predicts that heavy atoms (Z \u003E Z(c) approximate to 170) will undergo the process of atomic collapse where electrons sink into the positron continuum and a new family of so-called collapsing states emerges. The relativistic electrons in graphene exhibit the same physics but at a much lower critical charge (Z(c) approximate to 1) which has made it possible to confirm this phenomenon experimentally. However, there exist conflicting predictions on the effect of a magnetic field on atomic collapse. These theoretical predictions are based on the continuum Dirac-Weyl equation, which does not have an exact analytical solution for the interplay of a supercritical Coulomb potential and the magnetic field. Approximative solutions have been proposed, but because the two effects compete on similar energy scales, the theoretical treatment varies depending on the regime which is being considered. These limitations are overcome here by starting from a tight-binding approach and computing exact numerical results. By avoiding special limit cases, we found a smooth evolution between the different regimes. We predict that the atomic collapse effect persists even after the magnetic field is activated and that the critical charge remains unchanged. We show that the atomic collapse regime is characterized: (1) by a series of Landau level anticrossings and (2) by the absence of root B scaling of the Landau levels with regard to magnetic field strength.'));

Abstract: Using the non-equilibrium Green's function formalism, we study carrier transport through imperfect two-dimensional (2D) topological insulator (TI) ribbons. In particular, we investigate the effect of vacancy defects on the carrier transport in 2D TI ribbons with hexagonal lattice structure. To account for the random distribution of the vacancy defects, we present a statistical study of varying defect densities by stochastically sampling different defect configurations. We demonstrate that the topological edge states of TI ribbons are fairly robust against a high concentration (up to 2%) of defects. At very high defect densities, we observe an increased inter-edge interaction, mediated by the localisation of the edge states within the bulk region. This effect causes significant back-scattering of the, otherwise protected, edge-states at very high defect concentrations (>2%), resulting in a loss of conduction through the TI ribbon. We discuss how this coherent vacancy scattering can be used to our advantage for the development of TI-based transistors. We find that there is an optimal concentration of vacancies yielding an ON-OFF current ratio of up to two orders of magnitude. Finally, we investigate the importance of spin-orbit coupling on the robustness of the edge states in the TI ribbon and show that increased spin-orbit coupling could further increase the ON-OFF ratio.

Abstract: Heterostructures of low-dimensional semiconducting materials, such as transition metal dichalcogenides (MX2), are promising building blocks for future electronic and optoelectronic devices. The patterning of one MX2 material on top of another one is challenging due to their structural similarity. This prevents an intrinsic etch stop when conventional anisotropic dry etching processes are used. An alternative approach consist in a two-step process, where a sacrificial silicon layer is pre-patterned with a low damage plasma process, stopping on the underlying MoS2 film. The pre-patterned layer is used as sacrificial template for the formation of the top WS2 film. This study describes the optimization of a cyclic Ar/Cl-2 atomic layer etch process applied to etch silicon on top of MoS2, with minimal damage, followed by a selective conversion of the patterned Si into WS2. The impact of the Si atomic layer etch towards the MoS2 is evaluated: in the ion energy range used for this study, MoS2 removal occurs in the over-etch step over 1-2 layers, leading to the appearance of MoOx but without significant lattice distortions to the remaining layers. The combination of Si atomic layer etch, on top of MoS2, and subsequent Si-to-WS2 selective conversion, allows to create a WS2/MoS2 heterostructure, with clear Raman signals and horizontal lattice alignment. These results demonstrate a scalable, transfer free method to achieve horizontally individually patterned heterostacks and open the route towards wafer-level processing of 2D materials.

Abstract: Transition metal dichalcogenides, intercalated with transition metals, are studied for their potential applications as dilute magnetic semiconductors. We investigate the magnetic properties of WSe2 doped with third-row transition metals (Co, Cr, Fe, Mn, Ti and V). Using density functional theory in combination with Monte Carlo simulations, we obtain an estimate of the Curie or Neel temperature. We find that the magnetic ordering is highly dependent on the dopant type. While Ti and Cr-doped WSe2 have a ferromagnetic ground state, V, Mn, Fe and Co-doped WSe2 are antiferromagnetic in their ground state. For Fe doped WSe2, we find a high Curie-temperature of 327 K. In the case of V-doped WSe2, we find that there are two distinct magnetic phase transitions, originating from a frustrated in-plane antiferromagnetic exchange interaction and a ferromagnetic out-of-plane interaction. We calculate the formation energy and reveal that, in contrast to earlier reports, the formation energy is positive for the intercalated systems studied here. We also show that in the presence of W-vacancies, it becomes favorable for Ti, Fe, and Co to intercalate in WSe2.

Abstract: Among the large variety of two-dimensional (2D) materials discovered to date, elemental monolayers that host superconductivity are very rare. Using ab initio calculations we show that recently synthesized gallium monolayers, coined gallenene, are intrinsically superconducting through electron-phonon coupling. We reveal that Ga-100 gallenene, a planar monolayer isostructural with graphene, is the structurally simplest 2D superconductor to date, furthermore hosting topological edge states due to its honeycomb structure. Our anisotropic Eliashberg calculations show distinctly three-gap superconductivity in Ga-100, in contrast to the alternative buckled Ga-010 gallenene which presents a single anisotropic superconducting gap. Strikingly, the critical temperature (T ( c )) of gallenene is in the range of 7-10 K, exceeding the T ( c ) of bulk gallium from which it is exfoliated. Finally we explore chemical functionalization of gallenene with hydrogen, and report induced multigap superconductivity with an enhanced T ( c ) in the resulting gallenane compound.

Abstract: By using first-principles simulations, we investigate the interaction of a 2D silicon layer with two classes of chalcogenide-layered compounds, namely MoX2 and GaX (X = S, Se, Te). A rather weak (van der Waals) interaction between the silicene layers and the chalcogenide layers is predicted. We found that the buckling of the silicene layer is correlated to the lattice mismatch between the silicene layer and the MoX2 or GaX template. The electronic properties of silicene on these different templates largely depend on the buckling of the silicene layer: highly buckled silicene on MoS2 is predicted to be metallic, while low buckled silicene on GaS and GaSe is predicted to be semi-metallic, with preserved Dirac cones at the K points. These results indicate new routes for artificially engineering silicene nanosheets, providing tailored electronic properties of this 2D layer on non-metallic substrates. These non-metallic templates also open the way to the possible integration of silicene in future nanoelectronic devices.

Abstract: Motivated by the recent synthesis of single layer TiSe2, we used state-of-the-art density functional theory calculations, to investigate the structural and electronic properties of zigzag and armchair-edged nanoribbons (NRs) of this material. Our analysis reveals that, differing from ribbons of other ultra-thin materials such as graphene, TiSe2 NRs have some distinctive properties. The electronic band gap of the NRs decreases exponentially with the width and vanishes for ribbons wider than 20 angstrom. For ultranarrow zigzag-edged NRs we find odd-even oscillations in the band gap width, although their band structures show similar features. Moreover, our detailed magnetic-ground-state analysis reveals that zigzag and armchair edged ribbons have non-magnetic ground states. Passivating the dangling bonds with hydrogen at the edges of the structures influences the band dispersion. Our results shed light on the characteristic properties of T phase NRs of similar crystal structures.