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: Based on extensive Langevin Dynamics simulations we investigate the structural properties of a two-dimensional ensemble of magnetic rods with a peapod-like morphology, i.e, rods consisting of aligned single dipolar beads. Self-assembled configurations are studied for different directions of the dipole with respect to the rod axis. We found that with increasing misalignment of the dipole from the rod axis, the smaller the packing fraction at which the percolation transition is found. For the same density, the system exhibits different aggregation states for different misalignment. We also study the stability of the percolated structures with respect to temperature, which is found to be affected by the microstructure of the assembly of rods.

Abstract: Starting from the simplified linear combination of atomic orbitals method in combination with first-principles calculations, we construct a tight-binding (TB) model in the two-centre approximation for borophene and hydrogenated borophene (borophane). The Slater and Koster approach is applied to calculate the TB Hamiltonian of these systems. We obtain expressions for the Hamiltonian and overlap matrix elements between different orbitals for the different atoms and present the SK coefficients in a nonorthogonal basis set. An anisotropic Dirac cone is found in the band structure of borophane. We derive a Dirac low-energy Hamiltonian and compare the Fermi velocities with that of graphene.

Abstract: In a recent paper [Phys. Rev. B 89, 155403 (2014)], the authors investigated the spectrum of a Coulomb impurity in graphene in the presence of magnetic and electric fields using the coupled series expansion approach. In the first part of their paper, they investigated how Coulomb impurity states collapse in the presence of a perpendicular magnetic field. We argue that the obtained spectrum does not give information about the atomic collapse and that their interpretation of the spectrum regarding atomic collapse is not correct. We also argue that the obtained results are only valid up to the dimensionless charge vertical bar alpha vertical bar = 0.5 and, to obtain correct results for alpha > 0.5, a proper regularization of the Coulomb interaction is required. Here we present the correct numerical results for the spectrum for arbitrary values of alpha.

Abstract: Starting from the shell structure in atoms and the significant correlation within electron pairs, we distinguish the exchange-correlation effects between two electrons of opposite spins occupying the same orbital from the average correlation among many electrons in a crystal. In the periodic potential of the crystal with lattice constant larger than the effective Bohr radius of the valence electrons, these correlated electron pairs can form a metastable energy band above the corresponding single-electron band separated by an energy gap. In order to determine if these metastable electron pairs can be stabilized, we calculate the many-electron exchange-correlation renormalization and the polaron correction to the two-band system with single electrons and electron pairs. We find that the electron-phonon interaction is essential to counterbalance the Coulomb repulsion and to stabilize the electron pairs. The interplay of the electron-electron and electron-phonon interactions, manifested in the exchange-correlation energies, polaron effects, and screening, is responsible for the formation of electron pairs (bipolarons) that are located on the Fermi surface of the single-electron band.

Abstract: We present a theoretical investigation of the plasmon-polariton modes in grating coupled graphene inside a Fabry-Perot cavity. The cavity or photon modes of the device are determined by the Finite Difference Time Domain (FDTD) simulations and the corresponding plasmon-polariton modes are obtained by applying a many-body self-consistent field theory. We find that in such a device structure, the electric field strength of the incident electromagnetic (EM) field can be significantly enhanced near the edges of the grating strips. Thus, the strong coupling between the EM field and the plasmons in graphene can be achieved and the features of the plasmon-polariton oscillations in the structure can be observed. It is found that the frequencies of the plasmon-polariton modes are in the terahertz (THz) bandwidth and depend sensitively on electron density which can be tuned by applying a gate voltage. Moreover, the coupling between the cavity photons and the plasmons in graphene can be further enhanced by increasing the filling factor of the device. This work can help us to gain an in-depth understanding of the THz plasmonic properties of graphene-based structures.

Abstract: Due to Klein tunneling, electrostatic confinement of electrons in graphene is not possible. This hinders the use of graphene for quantum dot applications. Only through quasi-bound states with finite lifetime has one achieved to confine charge carriers. Here, we propose that bilayer graphene with a local region of decoupled graphene layers is able to generate bound states under the application of an electrostatic gate. The discrete energy levels in such a quantum blister correspond to localized electron and hole states in the top and bottom layers. We find that this layer localization and the energy spectrum itself are tunable by a global electrostatic gate and that the latter also coincides with the electronic modes in a graphene disk. Curiously, states with energy close to the continuum exist primarily in the classically forbidden region outside the domain defining the blister. The results are robust against variations in size and shape of the blister which shows that it is a versatile system to achieve tunable electrostatic confinement in graphene. Published by AIP Publishing.

Abstract: There is an increasing need to understand interfaces between two-dimensional materials to realize an energy efficient boundary with low contact resistance and small heat dissipation. In this respect, we investigated the impact of charge and substitutional atom doping on the electronic transport properties of the hybrid metallic-semiconducting lateral junctions, formed between metallic (1T and 1T(d)) and semiconducting (1H) phases of MoS2 by means of first-principles and non-equilibrium Green function formalism based calculations. Our results clearly revealed the strong influence of the type of interface and crystallographic orientation of the metallic phase on the transport properties of these systems. The Schottky barrier height, which is the dominant mechanism for contact resistance, was found to be as large as 0.63 eV and 1.19 eV for holes and electrons, respectively. We found that armchair interfaces are more conductive as compared to zigzag termination due to the presence of the metallic Mo zigzag chains that are directed along the transport direction. In order to manipulate these barrier heights we investigated the influence of electron doping of the metallic part (i.e. 1T(d) -MoS2). We observed that the Fermi level of the hybrid system moves towards the conduction band of semiconducting 1H-MoS2 due to filling of 4d-orbital of metallic MoS2, and thus the Schottky barrier for electrons decreases considerably. Besides electron doping, we also investigated the effect of substitutional doping of metallic MoS2 by replacing Mo atoms with either Re or Ta. Due to its valency, Re (Ta) behaves as a donor (acceptor) and reduces the Schottky barrier for electrons (holes). Since Re and Ta based transition metal dichalcogenides crystallize in either the 1T(d) or 1T phase, substitutional doping with these atom favors the stabilization of the 1T(d) phase of MoS2. Co-doping of hybrid structure results in an electronic structure, which facilities easy dissociation of excitons created in the 1H part.

Abstract: We study topological phase transitions and topological quantum field effect transistor in monolayer molybdenum disulfide (MoS2) using a two-band Hamiltonian model. Without considering the quadratic (q(2)) diagonal term in the Hamiltonian, we show that the phase diagram includes quantum anomalous Hall effect, quantum spin Hall effect, and spin quantum anomalous Hall effect regions such that the topological Kirchhoff law is satisfied in the plane. By considering the q(2) diagonal term and including one valley, it is shown that MoS2 has a non-trivial topology, and the valley Chern number is non-zero for each spin. We show that the wave function is (is not) localized at the edges when the q(2) diagonal term is added (deleted) to (from) the spin-valley Dirac mass equation. We calculate the quantum conductance of zigzag MoS2 nanoribbons by using the nonequilibrium Green function method and show how this device works as a field effect topological quantum transistor.

Abstract: Recently, two stable allotropes (germancite and stancite) for the group IV elements (Ge and Sn) with a staggered layered dumbell structure were proposed to be three-dimensional (3D) topological Dirac semimetals [Phys. Rev. B 93, 241117 (2016)]. A pair of Dirac points is on the rotation axis away from the time-reversal invariant momentum, and the stability of the 3D bulk Dirac points is protected by the C-3 rotation symmetry. Here, we use the first principles calculations to investigate GexSny alloys which share the same rhombohedral crystal structure with the space group of D-3d(6). Six GexSny alloys are predicted to be energetically and dynamically stable, where (x, y) = (8, 6) and (6, 8) and the alpha and beta phases of (10, 4) and (4, 10). Our results demonstrate that all the six GexSny alloys are topological Dirac semimetals. The different nontrivial surface states and surface Fermi arcs are identified. Our work will substantially enrich the family of 3D Dirac semimetals which are within the reach of experimental realization. Published by AIP Publishing.

Abstract: Massless Dirac fermions (MDFs) emerge as quasiparticles in various novel materials such as graphene and topological insulators, and they exhibit several intriguing properties, of which Veselago focusing is an outstanding example with a lot of possible applications. However, up to now Veselago focusing merely occurred in p-n junction devices based on the isotropic MDF, which lacks the tunability needed for realistic applications. Here, motivated by the emergence of novel Dirac materials, we investigate the propagation behaviors of anisotropic MDFs in such a p-n junction structure. By projecting the Hamiltonian of the anisotropic MDF to that of the isotropic MDF and deriving an exact analytical expression for the propagator, precise Veselago focusing is demonstrated without the need for mirror symmetry of the electron source and its focusing image. We show a tunable focusing position that can be used in a device to probe masked atom-scale defects. This study provides an innovative concept to realize Veselago focusing relevant for potential applications, and it paves the way for the design of novel electron optics devices by exploiting the anisotropic MDF.

Abstract: We use a variational technique to study the many-body phase diagram of electrons and holes in n-doped and p-doped monolayer transition metal dichalcogenides (TMDs). We find a total of four different phases. (i) A fully spin polarized and valley polarized ferromagnetic state. (ii) A state with no global spin polarization but with spin polarization in each valley separately, i.e., spin-valley locking. (iii) A state with spin polarization in one of the valleys and little to no spin polarization in the other valley. (iv) A paramagnetic state with no valley polarization. These phases are separated by first-order phase transitions and are determined by the particle density and the dielectric constant of the substrate. We find that in the presence of a perpendicular magnetic field the four different phases persist. In the case of n-doped MoS2, a fifth phase, which is completely valley polarized but not spin polarized, appears for magnetic fields larger than 7 T and for magnetic fields larger than 23 T completely replaces the second phase.

Abstract: Coupled 2D sheets of electrons and holes are predicted to support novel quantum phases. Two experiments of Coulomb drag in electron-hole (e-h) double bilayer graphene (DBLG) have reported an unexplained and puzzling sign reversal of the drag signal. However, we show that this effect is due to the multiband character of DBLG. Our multiband Fermi liquid theory produces excellent agreement and captures the key features of the experimental drag resistance for all temperatures. This demonstrates the importance of multiband effects in DBLG: they have a strong effect not only on superfluidity, but also on the drag.

Abstract: The spectrum of a Coulomb impurity in bilayer graphene is investigated as function of the strength of a perpendicular magnetic field for different values of the angular quantum number m and for different values of the gate voltage. We point out fundamental differences between the results from the two-band and four-band model. The supercritical instability and fall-to-center phenomena are investigated in the presence of a magnetic field. We find that in the four-band model the fall-to-center phenomenon occurs as in monolayer graphene, while this is not the case in the two-band model. We find that in a magnetic field the supercritical instability manifests itself as a series of anticrossings in the hole part of the spectrum for states coming from the low-energy band. However, we also find very distinct anticrossings in the electron part of the spectrum that continue into the hole part, which are related to the higher energy band of the four-band model. At these anticrossings, we find a very sharp peak in the probability density close to the impurity, reminiscent for the fall-to-center phenomenon. In this paper, these peculiar and interesting effects are studied for different magnetic field, interlayer coupling, and bias potential strengths.

Abstract: Edge functionalization of graphene nanoribbons with nitrogen atoms for various adatom configurations at armchair and zigzag edges are investigated. We provide comprehensive information on the electronic and magnetic properties and investigate the stability of the various systems. Two types of rippling of the nanoribbons, namely edge and bulk rippling depending on the sign of edge stress induced at the edge, are found. They are found to play the decisive role for the stability of the structures. We also propose a type of edge decoration in which every third nitrogen adatom at the zigzag edges is replaced by an oxygen atom. In this way, the electron count is compatible with a full aromatic structure, leading to additional stability and a disappearance of magnetism that is usually associated with zigzag nanoribbons.