Abstract: We study the influence of confinement on the atomic collapse due to a Coulomb impurity placed at the center of a graphene quantum dot of radius R. We apply the zigzag or infinite-mass boundary condition and consider both a point-size and a finite-size impurity. As a function of the impurity strength Za, the energy spectra are discrete. In the case of the zigzag boundary condition, the degenerate (with respect to the angular momentum m) zero-energy levels are pulled down in energy as Z alpha increases, and they remain below epsilon = – Z alpha. Our results show that the energy levels exhibit a 1/R dependence in the subcritical regime [Z alpha < |km + 1/2|, k = 1 (-1) for the K (K') valley]. In the supercritical regime (Z alpha > |km + 1/2|) we find a qualitatively very different behavior where the levels decrease as a function of R in a nonmonotonic manner. While the valley symmetry is preserved in the presence of the impurity, we find that the impurity breaks electron-hole symmetry. We further study the energy spectrum of zigzag quantum dots in gapped graphene. Our results show that as the gap increases, the lowest electron states are pushed into the gap by the impurity.

Abstract: We report on work function measurements on graphene, which is exfoliated over a predefined array of wells in silicon oxide, by Kelvin probe force microscopy operating in a vacuum. The obtained graphene sealed microchambers can support large pressure differences, providing controllable stretching of the nearly impermeable graphene membranes. These measurements allow detecting variations of the work function induced by the mechanical stresses in the suspended graphene where the work function varies linearly with the strain and changes by 62 +/- 2 meV for 1 percent of strain. Our related ab initio calculations result in a work function variation that is a factor of 1.4 larger than the experimental value. The limited discrepancy between the theory and the experiment can be accounted for by a charge transfer from the unstrained to the strained graphene regions. Published by AIP Publishing.

Abstract: Optical and electronic properties of black phosphorus strongly depend on the number of layers and type of stacking. Using first-principles calculations within the framework of density functional theory, we investigate the electronic properties of bilayer black phosphorus with an interlayer twist angle of 90 degrees. These calculations are complemented with a simple (k) over right arrow . (p) over right arrow model which is able to capture most of the low energy features and is valid for arbitrary twist angles. The electronic spectrum of 90 degrees twisted bilayer black phosphorus is found to be x-y isotropic in contrast to the monolayer. However x-y anisotropy, and a partial return to monolayer-like behavior, particularly in the valence band, can be induced by an external out-of-plane electric field. Moreover, the preferred hole effective mass can be rotated by 90 degrees simply by changing the direction of the applied electric field. In particular, a +0.4 (-0.4) V angstrom(1) out-of-plane electric field results in a similar to 60% increase in the hole effective mass along the y (x) axis and enhances the m(y)*/m(x)* (m(x)*/m(y)*) ratio as much as by a factor of 40. Our DFT and (k) over right arrow . (p) over right arrow simulations clearly indicate that the twist angle in combination with an appropriate gate voltage is a novel way to tune the electronic and optical properties of bilayer phosphorus and it gives us a new degree of freedom to engineer the properties of black phosphorus based devices.

Abstract: Two-dimensional group-V elemental materials have attracted widespread attention due to their nonzero band gap while displaying high electron mobility. Using first-principles calculations, we propose a series of new elemental bilayers with group-V elements (Bi, Sb, As). Our study reveals the dynamical stability of four-, six-, and eight-atom ring structures, demonstrating their possible coexistence in such bilayer systems. The proposed structures for Sb and As are large-gap semiconductors that are potentially interesting for applications in future nanodevices. The Bi structures have nontrivial topological properties with a direct nontrivial band gap. The nontrivial gap is shown to arise from a band inversion at the Brillouin zone center due to the strong intrinsic spin-orbit coupling in Bi atoms. Moreover, we demonstrate the possibility of tuning the properties of these materials by enhancing the ratio of six-atom rings to four-and eight-atom rings, which results in wider nontrivial band gaps and lower formation energies.

Abstract: Recently proposed accurate correlation energies are used to determine the phase diagram of strongly coupled electron-hole graphene bilayers. The control parameters of the phase diagram are the charge carrier density and the insulating barrier thickness separating the bilayers. In addition to the electron-hole superfluid phase we find two new inhomogeneous ground states, a one dimensional charge density wave phase and a coupled electron-hole Wigner crystal. The elementary crystal structure of bilayer graphene plays no role in generating these new quantum phases, which are completely determined by the electrons and holes interacting through the Coulomb interaction. The experimental parameters for the new phases lie within attainable ranges and therefore coupled electron-hole bilayer graphene presents itself as an experimental system where novel emergent many-body phases can be realized.

Abstract: Starting from a Hamiltonian with anharmonic coupling between in-plane acoustic displacements and outof-plane (flexural) modes, we derived coupled equations of motion for in-plane displacements correlations and flexural mode density fluctuations. Linear response theory and time-dependent thermal Green's functions techniques are applied in order to obtain different response functions. As external perturbations we allow for stresses and thermal heat sources. The displacement correlations are described by a Dyson equation where the flexural density distribution enters as an additional perturbation. The flexural density distribution satisfies a kinetic equation where the in-plane lattice displacements act as a perturbation. In the hydrodynamic limit this system of coupled equations is at the basis of a unified description of elastic and thermal phenomena, such as isothermal versus adiabatic sound motion and thermal conductivity versus second sound. The general theory is formulated in view of application to graphene, two-dimensional h-BN, and 2H-transition metal dichalcogenides and oxides.

Abstract: We propose a tight-binding (TB) model, that includes spin-orbit coupling (SOC), to describe the electronic properties of methyl-substituted germanane (GeCH3). This model gives an electronic spectrum in agreement with first principle results close to the Fermi level. Using the Z(2) formalism, we show that a topological phase transition from a normal insulator (NI) to a quantum spin Hall (QSH) phase occurs at 11.6% biaxial tensile strain. The sensitivity of the electronic properties of this system on strain, in particular its transition to the topological insulating phase, makes it very attractive for applications in strain sensors and other microelectronic applications.

Abstract: The electron energy spectrum in monolayer-bilayer-monolayer and in bilayer-monolayer-bilayer graphene structures is investigated and the effects of a perpendicular magnetic field and electric bias are studied. Different types of monolayer-bilayer interfaces are considered as zigzag (ZZ) or armchair (AC) junctions which modify considerably the bulk Landau levels (LLs) when the spectra are plotted as a function of the center coordinate of the cyclotron orbit. Far away from the two interfaces, one obtains the well-known LLs for extended monolayer or bilayer graphene. The LL structure changes significantly at the two interfaces or junctions where the valley degeneracy is lifted for both types of junctions, especially when the distance between them is approximately equal to the magnetic length. Varying the nonuniform bias and the width of this junction-to-junction region in either structure strongly influence the resulting spectra. Significant differences exist between ZZ and AC junctions in both structures. The densities of states (DOSs) for unbiased structures are symmetric in energy whereas those for biased structures are asymmetric. An external bias creates interface LLs in the gaps between the LLs of the unbiased system in which the DOS can be quite small. Such a pattern of LLs can be probed by scanning tunneling microscopy.

Abstract: The energy spectrum and eigenstates of single-layer black phosphorus nanoribbons in the presence of a perpendicular magnetic field and an in-plane transverse electric field are investigated by means of a tight-binding method, and the effect of different types of edges is examined analytically. A description based on a continuum model is proposed using an expansion of the tight-binding model in the long-wavelength limit. Thewave functions corresponding to the flatband part of the spectrum are obtained analytically and are shown to agree well with the numerical results from the tight-binding method for both narrow (10 nm) and wide (100 nm) nanoribbons. Analytical expressions for the critical magnetic field at which Landau levels are formed and the ranges of wave numbers in the dispersionless flatband segments in the energy spectra are derived. We examine the evolution of the Landau levels when an in-plane lateral electric field is applied, and we determine analytically how the edge states shift withmagnetic field. For wider nanoribbons, the conductance is shown to have a characteristic staircase shape in combined magnetic and electric fields. Some of the stairs in zigzag and skewed armchair nanoribbons originate from edge states that are found in the band gap.

Abstract: <script type='text/javascript'>document.write(unpmarked('Two-dimensional materials are known to harbour properties very different from those of their bulk counterparts. Recent years have seen the rise of atomically thin superconductors, with a caveat that superconductivity is strongly depleted unless enhanced by specific substrates, intercalants or adatoms. Surprisingly, the role in superconductivity of electronic states originating from simple free surfaces of two-dimensional materials has remained elusive to date. Here, based on first-principles calculations, anisotropic Eliashberg theory, and angle-resolved photoemission spectroscopy (ARPES), we show that surface states in few-monolayer MgB2 make a major contribution to the superconducting gap spectrum and density of states, clearly distinct from the widely known, bulk-like sigma-and pi-gaps. As a proof of principle, we predict and measure the gap opening on the magnesium-based surface band up to a critical temperature as high as similar to 30 K for merely six monolayers thick MgB2. These findings establish free surfaces as an unavoidable ingredient in understanding and further tailoring of superconductivity in atomically thin materials.'));

Abstract: <script type='text/javascript'>document.write(unpmarked('The motion of Abrikosov vortices in type-II superconductors results in a finite resistance in the presence of an applied electric current. Elimination or reduction of the resistance via immobilization of vortices is the \u0022holy grail\u0022 of superconductivity research. Common wisdom dictates that an increase in the magnetic field escalates the loss of energy since the number of vortices increases. Here we show that this is no longer true if the magnetic field and the current are applied parallel to each other. Our experimental studies on the resistive behavior of a superconducting Mo0.79Ge0.21 nanostrip reveal the emergence of a dissipative state with increasing magnetic field, followed by a pronounced resistance drop, signifying a reentrance to the superconducting state. Large-scale simulations of the 3D time-dependent Ginzburg-Landau model indicate that the intermediate resistive state is due to an unwinding of twisted vortices. When the magnetic field increases, this instability is suppressed due to a better accommodation of the vortex lattice to the pinning configuration. Our findings show that magnetic field and geometrical confinement can suppress the dissipation induced by vortex motion and thus radically improve the performance of superconducting materials.'));

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: Two-dimensional crystals with strong interactions between layers has attracted increasing attention in recent years in a variety of fields. In particular, the growth of a single layer of oxide materials (e.g., MgO, CaO, NiO, and MnO) over metallic substrates were found to display different physical properties than their bulk. In this study, we report on the physical properties of a single layer of metallic oxide materials and compare their properties with their bulk and other two-dimensional (2D) crystals. We found that the planar structure of metallic monoxides are unstable whereas the buckled structures are thermodynamically stable. Also, the 2D-MnO and NiO exhibit different magnetic (ferromagnetic) and optical properties than their bulk, whereas band-gap energy and linear stiffness are found to be decreasing from NiO to MgO. Our findings provide insight into oxide thin-film technology applications.

Abstract: The interaction of monolayer graphene with specific substrates may break its sublattice symmetry and results in unidirectional chiral states with opposite group velocities in the different Dirac cones (Zarenia et al 2012 Phys. Rev. B 86 085451). Taking advantage of this feature, we propose a valley filter based on a transversal mass kink for low energy electrons in graphene, which is obtained by assuming a defect region in the substrate that provides a change in the sign of the substrate-induced mass and thus creates a non-biased channel, perpendicular to the kink, for electron motion. By solving the time-dependent Schrodinger equation for the tight-binding Hamiltonian, we investigate the time evolution of a Gaussian wave packet propagating through such a system and obtain the transport properties of this graphene-based substrate-induced quantum point contact. Our results demonstrate that efficient valley filtering can be obtained, provided: (i) the electron energy is sufficiently low, i.e. with electrons belonging mostly to the lowest sub-band of the channel, and (ii) the channel length (width) is sufficiently long (narrow). Moreover, even though the transmission probabilities for each valley are significantly affected by impurities and defects in the channel region, the valley polarization in this system is shown to be robust against their presence.