Abstract: The intrinsic field effect, the change in surface conductance with an applied transverse electric field, of prototypal strongly correlated VO(2) has remained elusive. Here we report its measurement enabled by epitaxial VO(2) and atomic layer deposited high-kappa dielectrics. Oxygen migration, joule heating, and the linked field-induced phase transition are precluded. The field effect can be understood in terms of field-induced carriers with densities up to approximately 5x10(13) cm(-2) which are trongly localized, as shown by their low, thermally activated mobility ( approximately 1x10(-3) cm(2)/V s at 300 K). These carriers show behavior consistent with that of Holstein polarons and strongly impact the (opto)electronics of VO(2).

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: Using density-functional theory in combination with the nonequilibrium Green's function formalism, we study the effect of iodide/chloride and iodide/bromide mixing on the electronic transport in lead based organometallic perovskite CH3NH3PbI3, which is known to be an effective tool to tune the electronic and optical properties of such materials. We found that depending on the level and position of the halide mixing, the electronic transport can be increased by more than a factor of 4 for a given voltage biasing. The largest current is observed for small concentration of bromide substitutions located at the equatorial sites. However, full halide substitution has a negative effect on the transport properties of this material: the current drops by an order of magnitude for both CH3NH3PbCl3 and CH3NH3PbBr3 samples. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: Simulations are carried out based on the dynamical mean-field theory (DMFT) in order to investigate the properties of correlated thin films for various values of the chemical potential, temperature, interaction strength, and applied transverse electric field. Application of a sufficiently strong field to a thin film at half filling leads to the appearance of conducting regions near the surfaces of the film, whereas in doped slabs the application of a field leads to a conductivity enhancement on one side of the film and a gradual transition to the insulating state on the opposite side. In addition to the inhomogeneous DMFT, a local density approximation (LDA) is considered in which the particle density n, quasiparticle residue Z, and spectral weight at the Fermi level A(ω=0) of each layer are approximated by a homogeneous bulk environment. A systematic comparison between the two approaches reveals that the less expensive LDA results are in good agreement with the DMFT approach, except close to the metal-to-insulator transition points and in the layers immediately at the film surfaces. LDA values for n are overall more reliable than those for Z and A(ω=0). The hysteretic behavior (memory effect) characteristic of the bulk doping driven Mott transition persists in the slab.

Abstract: Using the tight-binding approach, we calculate the electronic structure of triangular and hexagonal MoS2 quantum dots. Due to the orbital asymmetry we show that it is possible to form quantum dots with the same shape but having different electronic properties. The electronic states of triangular and hexagonal quantum dots are explored, as well as the local and total density of states and the convergence towards the bulk spectrum with dot size is investigated. Our calculations show that: (1) edge states appear in the band gap, (2) that there are a larger number of electronic states in the conduction band as compared to the valence band, and (3) the relative number of edge states decreases with increasing dot size.

Abstract: Often real samples of graphene consist of islands of both monolayer and bilayer graphene. Bound states in such hybrid quantum dots are investigated for (i) a circular single-layer graphene quantum dot surrounded by an infinite bilayer graphene sheet and (ii) a circular bilayer graphene quantum dot surrounded by an infinite single-layer graphene. Using the continuum model and applying zigzag boundary conditions at the single-layer-bilayer graphene interface, we obtain analytical results for the energy levels and the corresponding wave spinors. Their dependence on perpendicular magnetic and electric fields are studied for both types of quantum dots. The energy levels exhibit characteristics of interface states, and we find anticrossings and closing of the energy gap in the presence of a bias potential.

Abstract: Fluorinated graphene is one of the most important derivatives of graphene and has been found to have great potential in optoelectronic and photonic nanodevices. However, the stability of F atoms on fluorinated graphene under different conditions, which is essential to maintain the desired properties of fluorinated graphene, is still unclear. In this work, we investigate the diffusion of F atoms on pristine graphene, graphene with defects, and at graphene/fluorographene interfaces by using density functional theory calculations. We find that an isolated F atom diffuses easily on graphene, but those F atoms can be localized by inducing vacancies or absorbates in graphene and by creating graphene/fluorographene interfaces, which would strengthen the binding energy of F atoms on graphene and increase the diffusion energy barrier of F atoms remarkably.

Abstract: Stable monolayer materials based on existing, well known and stable two-dimensional crystal fluorographene and molybdenum disulfide are predicted to exhibit a huge magnetocrystalline anisotropy when functionalized with adsorbed transition metal atoms at vacant sides. Ab initio calculations within the density-functional theory formalism were performed to investigate the adsorption of the transitional metals in a single S (or F) vacancy of monolayer molybdenum disulfide (or fluorographene). We found strong bonding of the transitional metal atoms to the vacant sites with binding energies ranging from 2.5 to 5.2 eV. Our calculations revealed that these systems with adsorbed metal atoms exhibit a magnetic anisotropy, specifically the structures including Os and Ir show a giant magnetocrystalline anisotropy energy of 31-101 meV. Our results demonstrate the possibility of obtaining stable monolayer materials with huge magnetocrystalline anisotropy based on preexisting, well known and stable two-dimensional crystals: fluorographene and molybdenum disulfide. We believe that the results obtained here are useful not only for deeper understanding of the origin of magnetocrystalline anisotropy but also for the design of monolayer optoelectronic devices with novel functionalities.

Abstract: Nitrogenated holey graphene (NHG) is a recently synthesized two-dimensional material. In this paper the structural and electronic properties of heterostructures of graphene and NHG are investigated using first-principles and tight-binding calculations. Due to the lattice mismatch between NHG and graphene, the formation of a moire pattern is preferred in the graphene/NHG heterostructure, instead of a lattice-coherent structure. In moire-patterned graphene/NHG, the band gap opening at the K point is negligible, and the linear band dispersion of graphene survives. Applying an electric field modifies the coupling strength between the two atomic layers. The Fermi velocity upsilon(F) is reduced as compared to the one of pristine graphene, and its magnitude depends on the twist angle theta between graphene and NHG: For theta = 0 degrees, upsilon(F) is 30% of that of graphene, and it increases rapidly to a value of 80% with increasing theta. The heterostructure exhibits electron-hole asymmetry in upsilon(F), which is large for small theta. In NHG encapsulated between two graphene layers, a “Dirac ring” appears around the K point. Its presence is robust with respect to the relative stacking of the two graphene layers. These findings can be useful for future applications of graphene/NHG heterostructures.

Abstract: We investigate how the interplay of quantum confinement and level broadening caused by disorder affects superconducting correlations in ultra-small metallic grains. We use the electron-phonon interaction-induced electron mass renormalization and the reduced static-path approximation of the BCS formalism to calculate the critical temperature as a function of the grain size. We show how the strong electron-impurity scattering additionally smears the peak structure in the electronic density of states of a metallic grain and imposes additional limits on the critical temperature under strong quantum confinement.

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.