Abstract: Graphene’s remarkable electronic properties have fuelled the vision of a graphene-based platform for lighter, faster and smarter electronics and computing applications. One of the challenges is to devise ways to tailor graphene’s electronic properties and to control its charge carriers. Here we show that a single-atom vacancy in graphene can stably host a local charge and that this charge can be gradually built up by applying voltage pulses with the tip of a scanning tunnelling microscope. The response of the conduction electrons in graphene to the local charge is monitored with scanning tunnelling and Landau level spectroscopy, and compared to numerical simulations. As the charge is increased, its interaction with the conduction electrons undergoes a transition into a supercritical regime where itinerant electrons are trapped in a sequence of quasi-bound states which resemble an artificial atom. The quasi-bound electron states are detected by a strong enhancement of the density of states within a disc centred on the vacancy site which is surrounded by halo of hole states. We further show that the quasi-bound states at the vacancy site are gate tunable and that the trapping mechanism can be turned on and off, providing a mechanism to control and guide electrons in graphene.

Abstract: Different allotropes of phosphorene are possible of which black and blue phosphorus are the most stable. While blue phosphorus has isotropic properties, black phosphorus is strongly anisotropic in its electronic and optical properties due to its anisotropic crystal structure. In this work, we systematically investigated the lattice thermal properties of black and blue phosphorene by using first-principles calculations based on the quasiharmonic approximation approach. Similar to the optoelectronic and electronic properties, we predict that black phosphorene has highly anisotropic thermal properties, in contrast to the blue phase. The linear thermal expansion coefficients along the zigzag and armchair direction differ up to 20% in black phosphorene. The armchair direction of black phosphorene is more expandable as compared to the zigzag direction and the biaxial expansion of blue phosphorene under finite temperature. Our comparative analysis reveals that the inclusion of finite-temperature effects makes the blue phase thermodynamically more stable over the black phase above 135 K.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: Water confined between two graphene layers with a separation of a few A forms a layered two-dimensional ice structure. Using large scale molecular dynamics simulations with the adoptable ReaxFF interatomic potential we found that flat monolayer ice with a rhombic-square structure nucleates between the graphene layers which is nonpolar and nonferroelectric. We provide different energetic considerations and H-bonding results that explain the interlayer and intralayer properties of two-dimensional ice. The controversial AA stacking found experimentally [Algara-Siller et al., Nature (London) 519, 443 (2015)] is consistent with our minimum-energy crystal structure of bilayer ice. Furthermore, we predict that an odd number of layers of ice has the same lattice structure as monolayer ice, while an even number of ice layers exhibits the square ice AA stacking of bilayer ice.

Abstract: The electron states in axially symmetric quantum wires are computed by means of the effective-mass Schrodinger equation, which is written in cylindrical coordinates phi, rho, and z. We show that a direct discretization of the Schrodinger equation by central finite differences leads to a non-symmetric Hamiltonian matrix. Because diagonalization of such matrices is more complex it is advantageous to transform it in a symmetric form. This can be done by the Liouville-like transformation proposed by Rizea et al. (2008), which replaces the wave function psi(rho) with the function F(rho) = psi(rho)root rho and transforms the Hamiltonian accordingly. Even though a symmetric Hamiltonian matrix is produced by this procedure, the computed wave functions are found to be inaccurate near the origin, and the accuracy of the energy levels is not very high. In order to improve on this, we devised a finite-difference scheme which discretizes the Schrodinger equation in the first step, and then applies the Liouville-like transformation to the difference equation. Such a procedure gives a symmetric Hamiltonian matrix, resulting in an accuracy comparable to the one obtained with the finite element method. The superior efficiency of the new finite-difference scheme (FDM) is demonstrated for a few p-dependent one-dimensional potentials which are usually employed to model the electron states in free-standing and core shell quantum wires. The new scheme is compared with the other FDM schemes for solving the effective-mass Schrodinger equation, and is found to deliver energy levels with much smaller numerical error for all the analyzed potentials. It also gives more accurate results than the scheme of Rizea et al., except for the ground state of an infinite rectangular potential in freestanding quantum wires. Moreover, the PT symmetry is invoked to explain similarities and differences between the considered FDM schemes. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: Employing density functional theory-based methods, we investigate monolayer and bilayer structures of hexagonal SnS2, which is a recently synthesized monolayer metal dichalcogenide. Comparison of the 1H and 1T phases of monolayer SnS2 confirms the ground state to be the 1T phase. In its bilayer structure we examine different stacking configurations of the two layers. It is found that the interlayer coupling in bilayer SnS2 is weaker than that of typical transition-metal dichalcogenides so that alternative stacking orders have similar structural parameters and they are separated with low energy barriers. A possible signature of the stacking order in the SnS2 bilayer has been sought in the calculated absorbance and reflectivity spectra. We also study the effects of the external electric field, charging, and loading pressure on the characteristic properties of bilayer SnS2. It is found that (i) the electric field increases the coupling between the layers at its preferred stacking order, so the barrier height increases, (ii) the bang gap value can be tuned by the external E field and under sufficient E field, the bilayer SnS2 can become a semimetal, (iii) the most favorable stacking order can be switched by charging, and (iv) a loading pressure exceeding 3 GPa changes the stacking order. The E-field tunable band gap and easily tunable stacking sequence of SnS2 layers make this 2D crystal structure a good candidate for field effect transistor and nanoscale lubricant applications.

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: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.

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.

Abstract: Within a tight binding approach we investigate the energy levels of hexagonal and triangular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We study AA- and AB-(Bernal) stacked BLG QDs and obtain the energy levels in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). Our results show that the size dependence of the energy levels is different from that of monolayer graphene QDs. The energy spectrum of AB-stacked BLG QDs with zigzag edges exhibits edge states which spread out into the opened energy gap in the presence of a perpendicular electric field. We found that the behavior of these edges states is different for the hexagonal and triangular geometries. In the case of AA-stacked BLG QDs, the electron and hole energy levels cross each other in both cases of armchair and zigzag edges as the dot size or the applied bias increases.

Abstract: Using first principles calculations we investigate the structural and electronic properties of interfaces between fluorographane and MoS2. Unsymmetrical functionalization of graphene with H and F results in an intrinsic dipole moment perpendicular to the plane of the buckled graphene skeleton. Depending on the orientation of this dipole moment, the electronic properties of a physically absorbed MoS2 monolayer can be switched from n-to p-type or vice versa. We show that one can realize vanishing n-type/p-type Schottky barrier heights when contacting MoS2 to fluorographane. By applying a perpendicular electric field, the size of the Schottky barrier and the degree of doping can be tuned. Our calculations indicate that a fluorographane monolayer is a promising candidate for bipolar doping of MoS2, which is vital in the design of novel technological applications based on two-dimensional materials.

Abstract: Using density functional theory, we study the electronic properties of a graphene-silicene bilayer (GSB). A single layer of silicene binds to the graphene layer with adhesion energy of about 25 meV/atom. This adhesion energy between the two layers follows accurately the well-known -1/z(2) dispersion energy as found between two infinite parallel plates. In small flakes of GSB with hydrogenated edges, negative charge is transferred from the graphene layer to the silicene layer, producing a permanent and a switchable polar bilayer, while in an infinite GSB, the negative charge is transferred from the silicene layer to the graphene layer. The graphene-silicene bilayer is a good candidate for a nanocapacitor with piezoelectric capabilities. We found that the permanent dipole of the bilayer can be tuned by an external perpendicular electric field.

Abstract: Using the tight-binding approach, we investigate confined states in two different hybrid monolayer-bilayer systems: (i) a hexagonal monolayer area surrounded by bilayer graphene in the presence of a perpendicularly applied electric field and (ii) a hexagonal bilayer graphene dot surrounded by monolayer graphene. The dependence of the energy levels on dot size and external magnetic field is calculated. We find that the energy spectrum for quantum dots with zigzag edges consists of states inside the gap which range from dot-localized states, edge states, to mixed states coexisting together, whereas for dots with armchair edges, only dot-localized states are observed.

Abstract: The Majorana zero modes and persistent spin current in mesoscopic d-wave-superconducting loops with spin-orbit (SO) interaction are investigated by numerically solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. For some appropriate strength of the SO coupling, Majorana zero-energy states and sharp jumps of the spin-polarized currents can be observed when the highest energy levels cross the Fermi energy in the spectrum, leading to spin currents with opposite chirality flowing near the inner and outer edges of the sample. When the threaded magnetic flux turns on, four flux-dependent patterns of the persistent spin current with step-like features show up, accompanied by Majorana edge modes at flux values where the energy gap closes. Moreover, the Majorana zero mode is highly influenced by the direction of the Zeeman field. A finite in-plane field can lead to the gap opening since the inversion symmetry is broken. Remarkably, multiple Majorana zero-energy states occur in the presence of an out-of-plane field h(z), and the number of steps in the spin current evolution can be effectively tuned by the field strength due to the shift of Majorana zero modes. Finally, when the loop sample contains surface indentation defects, zero-energy modes can always show up in the presence of an appropriate h(z). Interestingly, multiple Majorana states may be present in the system with a corner defect even if h(z) = 0.