Abstract: We provide a comprehensive analysis of the prominent tight-binding (TB) models for transition metal dichalcogenides (TMDs) available in the literature. We inspect the construction of these TB models, discuss their parameterization used and conduct a thorough comparison of their effectiveness in capturing important electronic properties. Based on these insights, we propose a novel TB model for TMDs designed for enhanced computational efficiency. Utilizing MoS2 as a representative case, we explain why specific models offer a more accurate description. Our primary aim is to assist researchers in choosing the most appropriate TB model for their calculations on TMDs.

Abstract: We explore the flatness of conduction and valence bands of interlayer excitons in MoS2/WSe2 van der Waals heterobilayers, tuned by interlayer twist angle, pressure, and external electric field. We employ an efficient continuum model where the moire pattern from lattice mismatch and/or twisting is represented by an equivalent mesoscopic periodic potential. We demonstrate that the mismatch moire potential is too weak to produce significant flattening. Moreover, we draw attention to the fact that the quasi-particle effective masses around the Gamma-point and the band flattening are reduced with twisting. As an alternative approach, we show (i) that reducing the interlayer distance by uniform vertical pressure can significantly increase the effective mass of the moire hole, and (ii) that the moire depth and its band flattening effects are strongly enhanced by accessible electric gating fields perpendicular to the heterobilayer, with resulting electron and hole effective masses increased by more than an order of magnitude – leading to record-flat bands. These findings impose boundaries on the commonly generalized benefits of moire twistronics, while also revealing alternative feasible routes to achieve truly flat electron and hole bands to carry us to strongly correlated excitonic phenomena on demand.

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: We report the formation of nanobubbles on graphene with a radius of the order of 1 nm, using ultralow energy implantation of noble gas ions (He, Ne, Ar) into graphene grown on a Pt(111) surface. We show that the universal scaling of the aspect ratio, which has previously been established for larger bubbles, breaks down when the bubble radius approaches 1 nm, resulting in much larger aspect ratios. Moreover, we observe that the bubble stability and aspect ratio depend on the substrate onto which the graphene is grown (bubbles are stable for Pt but not for Cu) and trapped element. We interpret these dependencies in terms of the atomic compressibility of the noble gas as well as of the adhesion energies between graphene, the substrate, and trapped atoms.

Abstract: We examine from first principles the flexoelectric properties of phosphorene nanoribbons under mechanical bending along armchair and zigzag directions. In both cases we find that the radial polarization depends linearly on the strain gradient. The flexoelectricity along the armchair direction is over 40% larger than along the zigzag direction. The obtained flexoelectric coefficients of phosphorene are four orders of magnitude larger than those of graphene and comparable to transition metal dichalcogenides. Analysis of charge density shows that the flexoelectricity mainly arises from the pz orbitals of phosphorus atoms. The electron mobilities in bent phosphorene can be enhanced by over 60% along the armchair direction, which is significantly higher than previous reports of mobility tuned by uniaxial strain. Our results indicate phosphorene is a candidate for a two-dimensional material applicable in flexible-electronic devices.

Abstract: The strain fields of periodically buckled graphene induce a periodic pseudomagnetic field (PMF) that modifies the electronic band structure. From the geometry, amplitude, and period of the periodic pseudomagnetic field, we determine the necessary conditions to access the regime of correlated phases by examining the band flattening. As compared to twisted bilayer graphene the proposed system has the advantages that (1) only a single layer of graphene is needed, (2) one is not limited to hexagonal superlattices, and (3) narrower flat bandwidth and larger separation between flat bands can be induced. We, therefore, propose that periodically strained graphene single layers can become a platform for the exploration of exotic many-body phases.

Abstract: Two-dimensional atomic crystals can radically change their properties in response to external influences, such as substrate orientation or strain, forming materials with novel electronic structure(1-5). An example is the creation of weakly dispersive, 'flat' bands in bilayer graphene for certain 'magic' angles of twist between the orientations of the two layers(6). The quenched kinetic energy in these flat bands promotes electron-electron interactions and facilitates the emergence of strongly correlated phases, such as superconductivity and correlated insulators. However, the very accurate fine-tuning required to obtain the magic angle in twisted-bilayer graphene poses challenges to fabrication and scalability. Here we present an alternative route to creating flat bands that does not involve fine-tuning. Using scanning tunnelling microscopy and spectroscopy, together with numerical simulations, we demonstrate that graphene monolayers placed on an atomically flat substrate can be forced to undergo a buckling transition(7-9), resulting in a periodically modulated pseudo-magnetic field(10-14), which in turn creates a 'post-graphene' material with flat electronic bands. When we introduce the Fermi level into these flat bands using electrostatic doping, we observe a pseudogap-like depletion in the density of states, which signals the emergence of a correlated state(15-17). This buckling of two-dimensional crystals offers a strategy for creating other superlattice systems and, in particular, for exploring interaction phenomena characteristic of flat bands. Buckled monolayer graphene superlattices are found to provide an alternative to twisted bilayer graphene for the study of flat bands and correlated states in a carbon-based material.

Abstract: A periodic spatial modulation, as created by a moire pattern, has been extensively studied with the view to engineer and tune the properties of graphene. Graphene encapsulated by hexagonal boron nitride (hBN) when slightly misaligned with the top and bottom hBN layers experiences two interfering moire patterns, resulting in a so-called supermoire (SM). This leads to a lattice and electronic spectrum reconstruction. A geometrical construction of the nonrelaxed SM patterns allows us to indicate qualitatively the induced changes in the electronic properties and to locate the SM features in the density of states and in the conductivity. To emphasize the effect of lattice relaxation, we report band gaps at all Dirac-like points in the hole doped part of the reconstructed spectrum, which are expected to be enhanced when including interaction effects. Our result is able to distinguish effects due to lattice relaxation and due to the interfering SM and provides a clear picture on the origin of recently experimentally observed effects in such trilayer heterostuctures.

Abstract: The mechanical and electronic properties of a graphene membrane placed on top of a triangular superlattice of nanopillars are investigated. We use molecular dynamics simulations to access the deformation fields and the tight-binding approaches to calculate the electronic properties. Ripples form in the graphene layer that span across the unit cell, connecting neighboring pillars, in agreement with recent experiments. We find that the resulting pseudo-magnetic field (PMF) varies strongly across the unit cell. We investigate the dependence of PMF on unit cell boundary conditions, height of the pillars, and the strength of the van der Waals interaction between graphene and the substrate. We find direct correspondence with typical experiments on pillars, showing intrinsic “slack” in the graphene membrane. PMF values are confirmed by the local density of states calculations performed at different positions of the unit cell showing pseudo-Landau levels with varying spacings. Our findings regarding the relaxed membrane configuration and the induced strains are transferable to other flexible 2D membranes.

Abstract: The in-plane dc conductivity of twisted bilayer graphene is calculated using an expansion of the real-space Kubo-Bastin conductivity in terms of Chebyshev polynomials. We investigate within a tight-binding approach the transport properties as a function of rotation angle, applied perpendicular electric field, and vacancy disorder. We find that for high-angle twists, the two layers are effectively decoupled, and the minimum conductivity at the Dirac point corresponds to double the value observed in monolayer graphene. This remains valid even in the presence of vacancies, hinting that chiral symmetry is still preserved. On the contrary, for low twist angles, the conductivity at the Dirac point depends on the twist angle and is not protected in the presence of disorder. Furthermore, for low angles and in the presence of an applied electric field, we find that the chiral boundary states emerging between AB and BA regions contribute to the dc conductivity, despite the appearance of localized states in the AA regions. The results agree qualitatively with recent transport experiments in low-angle twisted bilayer graphene.

Abstract: <script type='text/javascript'>document.write(unpmarked('Spin-triplet chiral p-wave superconductivity is typically described by a two-component order parameter, and as such is prone to unique emergent effects when compared to the standard single-component superconductors. Here we present the equilibrium phase diagram for small mesoscopic chiral p-wave superconducting disks in the presence of magnetic field, obtained by solving the microscopic Bogoliubov-de Gennes equations self-consistently. In the ultrasmall limit, the cylindrically symmetric giant-vortex states form the ground state of the system. However, with increasing sample size, the cylindrical symmetry is broken as the two components of the order parameter segregate into domains, and the number of fragmented domain walls between them characterizes the resulting states. Such domain walls are topological defects unique for the p-wave order, and constitute a dominant phase in the mesoscopic regime. Moreover, we find two possible types of domain walls, identified by their chirality-dependent interaction with the edge states.'));

Abstract: Chiral p-wave superconductors in applied magnetic field can exhibit more complex topological defects than just conventional superconducting vortices, due to the two-component order parameter (OP) and the broken time-reversal symmetry. We investigate the electronic properties of those exotic states, some of which contain clusters of one-component vortices in chiral components of the OP and/or exhibit skyrmionic character in the relative OP space, all obtained as a self-consistent solution of the microscopic Bogoliubov-de Gennes equations. We reveal the link between the local density of states (LDOS) of the novel topological states and the behavior of the chiral domain wall between the OP components, enabling direct identification of those states in scanning tunneling microscopy. For example, a skyrmion always contains a closed chiral domain wall, which is found to be mapped exactly by zero-bias peaks in LDOS. Moreover, the LDOS exhibits electron-hole asymmetry, which is different from the LDOS of conventional vortex states with same vorticity. Finally, we present the magnetic field and temperature dependence of the properties of a skyrmion, indicating that this topological defect can be surprisingly large in size, and can be pinned by an artificially indented nonsuperconducting closed path in the sample. These features are expected to facilitate the experimental observation of skyrmionic states, thereby enabling experimental verification of chirality in emerging superconducting materials.

Abstract: Experiments have confirmed that double monolayer graphene does not generate finite-temperature electron-hole superfluidity, because of very strong screening of the pairing attraction. The linear dispersing energy bands in monolayer graphene block any attempt to reduce the strength of the screening. We propose a hybrid device with two sheets of monolayer graphene in a modulated periodic perpendicular magnetic field. The field preserves the isotropic Dirac cones of the original monolayers but reduces the slope of the cones, making the monolayer Fermi velocity v(F) smaller. We demonstrate that with current experimental techniques, the reduction in vF can weaken the screening sufficiently to allow electron-hole superfluidity at measurable temperatures.

Abstract: We have grown epitaxial Cr-doped V2O3 thin films with Cr concentrations between 0% and 20% on (0001)-Al2O3 by oxygen-assisted molecular beam epitaxy. For the highly doped samples (>3%), a regular and monotonous increase of the resistance with decreasing temperature is measured. Strikingly, in the low doping samples (between 1% and 3%), a collapse of the insulating state is observed with a reduction of the low temperature resistivity by up to 5 orders of magnitude. A vacuum annealing at high temperature of the films recovers the low temperature insulating state for doping levels below 3% and increases the room temperature resistivity towards the values of Cr-doped V2O3 single crystals. It is well-know that oxygen excess stabilizes a metallic state in V2O3 single crystals. Hence, we propose that Cr doping promotes oxygen excess in our films during deposition, leading to the collapse of the low temperature insulating state at low Cr concentrations. These results suggest that slightly Cr-doped V2O3 films can be interesting candidates for field effect devices. (C) 2015 AIP Publishing LLC.

Abstract: We numerically investigate the mixed pairing symmetries inmesoscopic superconducting loops in the presence of spin correlations by solving the Bogoliubov-de Gennes equations self-consistently. The spatial variations of the superconducting order parameters and the spontaneous magnetization are determined by the band structure. When the threaded magnetic flux turns on, the charge and spin currents both emerge and depict periodic evolution. In the case of a mesoscopic loop with dominant triplet p(x) +/- ip(y)-wave symmetry, a slight change of the chemical potential may lead to novel flux-dependent evolution patterns of the ground-state energy and the magnetization. The spin-polarized currents show pronounced quantum oscillations with fractional periods due to the appearance of energy jumps in flux, accompanied with a steplike feature of the enhanced spin current. Particularly, at some appropriate flux, the peaks of the zero-energy local density of states clearly indicate the occurrence of the odd-frequency pairing. In the case of a superconducting loop with dominant singlet d(x2-y2)-wave symmetry, the spatial profiles of the zero-energy local density of states and the magnetization show spin-dependent features on different sample diagonals. Moreover, the evolution of the flux-induced spin current always exhibits an hc/e periodicity.