Abstract: <script type='text/javascript'>document.write(unpmarked('We show that delaminations in bilayer graphene (BLG) with electrostatically induced interlayer symmetry can provide one with ballistic channels for electrons with energies inside the electrostatically induced BLG gap. These channels are formed by a combination of valley-polarized evanescent states propagating along the delamination edges (which persist in the presence of a strong magnetic field) and standing waves bouncing between them inside the delaminated region (in a strong magnetic field, these transform into Landau levels in the monolayers). For inverted stackings in BLGs on the left and right of the delamination (AB-2ML-BA or BA-2ML-AB, where 2ML indicates two decoupled monolayers of graphene), the lowest-energy ballistic channels are gapless, have linear dispersion, and appear to be weakly topologically protected. When BLG stackings on both sides of the delamination are the same (AB-2ML-AB or BA-2ML-BA), the lowest-energy ballistic channels are gapped, with a gap epsilon(g) scaling as epsilon(g) alpha W-1 with delamination width and epsilon(g) alpha delta(-1) with the on-layer energy difference in the delaminated part of the structure. Depending on the width, delaminations may also support several \u0022higher-energy\u0022 waveguide modes. Our results are based on both the analytical study of the wave matching of Dirac states and tight-binding model calculations, and we analyze in detail the dependence of the delamination spectrum on the electrostatic conditions in the structure, such as the vertical displacement field.'));

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 tight-binding approach, we investigate the energy spectrum of square, triangular, and hexagonal MoS2 quantum dots (QDs) in the presence of a perpendicular magnetic field. Novel edge states emerge in MoS2 QDs, which are distributed over the whole edge which we call ring states. The ring states are robust in the presence of spin-orbit coupling (SOC). The corresponding energy levels of the ring states oscillate as a function of the perpendicular magnetic field which are related to Aharonov-Bohm oscillations. Oscillations in the magnetic field dependence of the energy levels and the peaks in the magneto-optical spectrum emerge (disappear) as the ring states are formed (collapsed). The period and the amplitude of the oscillation decrease with the size of the MoS2 QDs.

Abstract: Defects are almost inevitably present in realistic materials and defective materials are expected to exhibit very different properties than their nondefective (perfect) counterparts. Here, using a combination of the tight-binding approach and the scattering matrix formalism, we investigate the electronic transport properties of defective phosphorene nanoribbons (PNRs) containing atomic vacancies. We find that for both armchair PNRs (APNRs) and zigzag PNRs (ZPNRs), single vacancies can create quasilocalized states, which can affect their conductance. With increasing vacancy concentration, three different transport regimes are identified: ballistic, diffusive, and Anderson localized ones. In particular, ZPNRs that are known to be metallic due to the presence of edge states become semiconducting: edge conductance vanishes and transport gap appears due to Anderson localization. Moreover, we find that for a fixed vacancy concentration, both APNRs and ZPNRs of narrower width and/or longer length are more sensitive to vacancy disorder than their wider and/or shorter counterparts, and that for the same ribbon length and width, ZPNRs are more sensitive to vacancy disorder than APNRs.

Abstract: Structural, mechanical and electronic properties of two-dimensional single-layer hexagonal structures in the (111) crystal plane of IIIAs-ZnS systems (III = B, Ga and In) are studied by first-principles calculations based on density functional theory (DFT). Elastic and phonon dispersion relation display that 2D h-IIIAs systems (III = B, Ga and In) are both mechanical and dynamically stable. Electronic structures analysis show that the semiconducting nature of the 3D-IIIAs compounds is retained by their 2D single layer counterpart. Furthermore, density of states reveals the influence of sigma and pi bonding in the most stable geometry (planar or buckled) for 2D h-IIIAs systems. Calculations of elastic constants show that the Young's modulus, bulk modulus and shear modulus decrease for 2D h-IIIAs binary compounds as we move down on the group of elements of the periodic table. In addition, as the bond length between the neighboring cation-anion atoms increases, the 2D h-IIIAs binary compounds display less stiffness and more plasticity. Our findings can be used to understand the contribution of the r and p bonding in the most stable geometry (planar or buckled) for 2D h-IIIAs systems. Structural and electronic properties of h-IIIAs systems as a function of the number of layers have been also studied. It is shown that h-BAs keeps its planar geometry while both h-GAs and h-InAs retained their buckled ones obtained by their single layers. Bilayer h-IIIAs present the same bandgap nature of their counterpart in 3D. As the number of layers increase from 2 to 4, the bandgap width for layered h-IIIAs decreases until they become semimetal or metal. Interestingly, these results are different to those found for layered h-GaN. The results presented in this study for single and few-layer h-IIIAs structures could give some physical insights for further theoretical and experimental studies of 2D h-IIIV-like systems. (C) 2017 Elsevier B.V. All rights reserved.

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: The self rotation of a graphene flake over graphite is controlled by the size, initial misalignment and temperature. Using both ab initio calculations and molecular dynamics simulations, we investigate annealing effects on the self rotation of a graphene flake on a graphene substrate. The energy barriers for rotation and drift of a graphene flake over graphene is found to be smaller than 25 meV/atom which is comparable to thermal energy. We found that small flakes (of about similar to 4 nm) are more sensitive to temperature and initial misorientation angles than larger one (beyond 10 nm). The initial stacking configuration of the flake is found to be important for its dynamics and time evolution of misalignment. Large flakes, which are initially in the AA-or AB-stacking state with small misorientation angle, rotate and end up in the AB-stacking configuration. However small flakes can they stay in an incommensurate state specially when the initial misorientation angle is larger than 2 degrees. Our results are in agreement with recent experiments.

Abstract: By taking account of the electric-field-induced charge screening, a self-consistent calculation within the framework of the tight-binding approach is employed to obtain the electronic band structure of gated multilayer phosphorene and the charge densities on the different phosphorene layers. We find charge density and screening anomalies in single-gated multilayer phosphorene and electron-hole bilayers in dual-gated multilayer phosphorene. Due to the unique puckered lattice structure, both intralayer and interlayer charge screenings are important in gated multilayer phosphorene. We find that the electric-field tuning of the band structure of multilayer phosphorene is distinctively different in the presence and absence of charge screening. For instance, it is shown that the unscreened band gap of multilayer phosphorene decreases dramatically with increasing electric-field strength. However, in the presence of charge screening, the magnitude of this band-gap decrease is significantly reduced and the reduction depends strongly on the number of phosphorene layers. Our theoretical results of the band-gap tuning are compared with recent experiments and good agreement is found.

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: Landau levels (LLs) are modified by the Frohlich interaction which we investigate within the improved Wigner-Brillouin theory for energies both below and above the longitudinal-optical-continuum in monolayer MoS2.., WS2, MoSe2, and WSe2. Polaron corrections to the LLs are enhanced in monolayer MoS2 as compared to WS2. A series of levels are found at h omega(LO) + lh omega(c), and in addition, the Frohlich interaction lifts the degeneracy between the levels nh omega(c) and h omega(LO) + lh omega(c) resulting in an anticrossing. The screening effect due to the environment plays an important role in the polaron energy corrections, which are also affected by the effective thickness r(eff) parameter. The polaron anticrossing energy gap E-gap decreases with increasing effective thickness r(eff). Published by AIP Publishing.

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: The permeation of water between neighboring graphene oxide (GO) flakes, i.e., 2D nanochannels, are investigated using a simple model for the GO membrane. We simulate the hydrophilic behavior of nanocapillaries and study the effect of surface charge on the dynamical properties of water flow and the influence of Na+ and Cl- ions on water permeation. Our approach is based on extensive equilibrium molecular dynamics simulations to obtain a better understanding of water permeation through charged nanochannels in the presence of ions. We found significant change in the slippage dynamics of confined water such as a profound increase in viscosity/slip length with increasing charges over the surface. The slip length decreases one order of magnitude (i.e., 1/30) with increasing density of surface charge, while it increases by a factor of 2 with ion concentration. We found that commensurability induced by nanoconfinement plays an important role on the intrinsic dynamical properties of water.

Abstract: Thin superconducting films can exhibit negative magnetoresistance when an in-plane external magnetic field is aligned parallel with the transport current. We explain this effect as due to appearance of parallel vortices in the plain of the film at the first critical magnetic field H-c1 which leads to an enhancement of the superconducting properties and impedes the motion of the current induced perpendicular vortices. Our theoretical results are based on a numerical solution of the time-dependent and stationary 3D Ginzburg-Landau equations.