Abstract: Topological insulators (TIs) have attracted a lot of interest in recent years due to their topologically protected surface states, as well as exotic proximity-induced phenomena and device applications for TI heterostructures. Since the first experimental studies of TIs, angle-resolved photoemission spectra (ARPES) showed that the electronic structure of the topological surface states significantly changes as a function of time after cleavage. The origin and underlying mechanism of this ageing effect are still under debate, despite its importance. Here we investigate the evolution of the surface Dirac cone for Bi2Se3 films upon asymmetric potassium (K) adsorption, using density-functional theory and a tight-binding model. We find that the K adatoms induce short-ranged downward band bending within 2-3 nm from the surface, due to charge transfer from the adatoms to the TI. These findings are in contrast to earlier proposals in the literature, that propose a long-ranged downward band bending up to 15 nm from the surface. Furthermore, as the charge transfer increases, we find that a new Dirac cone, localized slightly deeper into the TI than the original one, appears at the K-adsorbed surface, originating from strong Rashba-split conduction-band states. Our results suggest possible reinterpretations of experiments because the new Dirac cone might have been observed in ARPES measurements instead of the original one that appears immediately after cleavage. Our findings are consistent with ARPES data and provide insight into building TI-heterostructure devices by varying the band-bending potential or film thickness.

Abstract: Ultrafine ZnO nanocrystals with a thickness down to 0.25 nm are grown by a metalorganic chemical vapor deposition method. Electronic band structures and native point defects of ZnO nanocrystals are studied by a combination of scanning tunneling microscopy/spectroscopy and first-principles density functional theory calculations. Below a critical thickness of nm ZnO adopts a graphitic-like structure and exhibits a wide band gap similar to its wurtzite counterpart. The hexagonal wurtzite structure, with a well-developed band gap evident from scanning tunneling spectroscopy, is established for a thickness starting from similar to 1.4 nm. With further increase of the thickness to 2 nm, V-O-V-Zn defect pairs are easily produced in ZnO nanocrystals due to the self-compensation effect in highly doped semiconductors.

Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

Abstract: Using first-principles density functional theory, we investigate the dynamical properties of the roomtemperature P21/n and ground-state P21/c phases of WO3. As a preliminary step, we assess the validity of various standard and hybrid functionals, concluding that the best description is achieved with the B1-WC hybrid functional while a reliable description can also be provided using the standard LDA functional. We also carefully rediscuss the structure and energetics of all experimentally observed and a few hypothetical metastable phases in order to provide deeper insight into the unusual sequence of phase transition of WO3 with temperature. Then, we provide a comprehensive theoretical study of the lattice dynamical properties of the P21/n and P21/c phases, reporting zone-center phonons, infrared and Raman spectra, as well as the full phonon dispersion curves, which attest to the dynamical stability of both phases. We carefully discuss the spectra, explaining the physical origin of their main features and evolution from one phase to another. We reveal a systematic connection between the dynamical and structural properties of WO3, highlighting that the number of peaks in the high-frequency range of the Raman spectrum appears as a fingerprint of the number of antipolar distortions that are present in the structure and a practical way to discriminate between the different phases.

Abstract: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.

Abstract: We calculate the quasiparticle energy spectrum of SnO2 within the GW approximation, properly taking into account the contribution of core levels to the energy corrections. The calculated fundamental gap is of 3.85 eV. We propose that the difference with respect to the experimental optical gap (3.6 eV) is due to excitonic effects in the latter. We further consider the effect applied on uniaxial pressure along the c-axis. Compared to GW, the effect of pressure on the quasiparticle energies and band gap is underestimated by the local-density approximation. The quasiparticle effective masses, however, appear to be well described by the latter.

Abstract: We consider a junction between two topological insulators and calculate the properties of the interface states with an effective low-energy Hamiltonian for topological insulators with a single cone on the surface. This system bears a close resemblance to bilayer graphene, as both result from the hybridization of Dirac cones. We find gapless interface states not only when the helicity directions of the topological surface states are oppositely oriented, but they can also exist if they are equally oriented. Furthermore, we find that the existence of the interface states can be understood from the closing of the bulk gap when the helicity changes orientation. Recently superluminal tachyonic excitations were also claimed to exist at the interface between topological insulators. However, here we show that these interface states do not exist. DOI: 10.1103/PhysRevB.87.115113

Abstract: We consider classical two-dimensional (2D) Coulomb clusters consisting of two species containing five particles with charge q1 and five with charge q2, respectively. Using Monte Carlo and molecular dynamics (MD) simulations, we investigated the ground state configurations as well as radial and angular displacements of particles as a function of temperature and their dependence on the ratio q = q2/q1. We found new configurations and a new multi-step melting behavior for q sufficiently different from the uniform charge limit q = 1.

Abstract: We evaluate the self-diffusion and transport diffusion of interacting particles in a discrete geometry consisting of a linear chain of cavities, with interactions within a cavity described by a free-energy function. Exact analytical expressions are obtained in the absence of correlations, showing that the self-diffusion can exceed the transport diffusion if the free-energy function is concave. The effect of correlations is elucidated by comparison with numerical results. Quantitative agreement is obtained with recent experimental data for diffusion in a nanoporous zeolitic imidazolate framework material, ZIF-8.

Abstract: We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self-and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self-and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

Abstract: We have performed a first-principles study of the p- and n-type conductivity in CuIn1−xGaxSe2 due to native point defects, based on the HSE06 hybrid functional. Band alignment shows that the band gap becomes larger with x due to the increasing conduction band minimum, rendering it hard to establish n-type conductivity in CuGaSe2. From the defect formation energies, we find that In/GaCu is a shallow donor, while VCu, VIn/Ga and CuIn/Ga act as shallow acceptors. Using the total charge neutrality of ionized defects and intrinsic charge carriers to determine the Fermi level, we show that under In-rich growth conditions InCu causes strongly n-type conductivity in CuInSe2. Under increasingly In-poor growth conditions, the conductivity type in CuInSe2 alters to p-type and compensation of the acceptors by InCu reduces, as also observed in photoluminescence experiments. In CuGaSe2, the native acceptors pin the Fermi level far away from the conduction band minimum, thus inhibiting n-type conductivity. On the other hand, CuGaSe2 shows strong p-type conductivity under a wide range of Ga-poor growth conditions. Maximal p-type conductivity in CuIn1−xGaxSe2 is reached under In/Ga-poor growth conditions, in agreement with charge concentration measurements on samples with In/Ga-poor stoichiometry, and is primarily due to the dominant acceptor CuIn/Ga.

Abstract: We introduce a modified version of the Hirshfeld charge analysis method and demonstrate its accurateness by calculating the charge transfer between the paramagnetic molecule NO2 and graphene. The charge transfer between paramagnetic molecules and a graphene layer as calculated with ab initio methods can crucially depend on the size of the supercell used in the calculation. This has important consequences for adsorption studies involving paramagnetic molecules such as NO2 physisorbed on graphene or on carbon nanotubes. © 2008 American Institute of Physics.

Abstract: We investigate scattering of the topological surface state of a three-dimensional time-reversal invariant topological insulator when graphene is deposited on the topological-insulator surface. Specifically, we consider the (111) surface of a Bi2Se3-like topological insulator. We present a low-energy model for the graphene-topological insulator heterostructure and we calculate the transmission probability at zigzag and armchair edges of the deposited graphene, and the conductance through graphene nanoribbon barriers, and show that its features can be understood from antiresonances in the transmission probability.

Abstract: We investigate the link between superconductivity and atomic vibrations in Nb3Sn films with a thickness ranging from 10 to 50 nm. The challenge of measuring the phonon density of states (PDOS) of these films has been tackled by employing the technique of nuclear inelastic scattering by Sn-119 isotopes to reveal the Sn-partial phonon density of states. With the support of ab initio calculations, we evaluate the effect of reduced film thickness on the PDOS. This approach allows us to estimate the changes in superconducting critical temperature T-c induced by phonon confinement, which turned out to be limited to a few tenths of K. The presented method is successful for the Nb3Sn system and paves the way for more systematic studies of the role of phonon confinement in Sn-containing superconductors.

Abstract: We investigated the effect of different stacking order of the four graphene layer system on the induced band gap when positively charged top and negatively charged back gates are applied to the system. A tight-binding approach within a self-consistent Hartree approximation is used to calculate the induced charges on the different graphene layers. We show that the electric field does not open an energy gap if the multilayer graphene system contains a trilayer part with the ABA Bernal stacking.