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: The calculation of defect levels in semiconductors within a density functional theory approach suffers greatly from the band gap problem. We propose a band gap correction scheme that is based on the separation of energy differences in electron addition and relaxation energies. We show that it can predict defect levels with a reasonable accuracy, particularly in the case of defects with conduction band character, and yet is simple and computationally economical. We apply this method to ZnO doped with group III elements (Al, Ga, In). As expected from experiment, the results indicate that Zn substitutional doping is preferred over interstitial doping in Al, Ga, and In-doped ZnO, under both zinc-rich and oxygen-rich conditions. Further, all three dopants act as shallow donors, with the +1 charge state having the most advantageous formation energy. Also, doping effects on the electronic structure of ZnO are sufficiently mild so as to affect little the fundamental band gap and lowest conduction bands dispersion, which secures their n-type transparent conducting behavior. A comparison with the extrapolation method based on LDA+U calculations and with the HeydScuseriaErnzerhof hybrid functional (HSE) shows the reliability of the proposed scheme in predicting the thermodynamic transition levels in shallow donor systems.

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: It is well known that, in bulk, the solution of the Bogoliubovde Gennes equations is the same whether or not the HartreeFock term is included. Here the HartreeFock potential is position independent and so gives the same contribution to both the single-electron energies and the Fermi level (the chemical potential). Thus, the single-electron energies measured from the Fermi level (they control the solution) stay the same. This is not the case for nanostructured superconductors, where quantum confinement breaks the translational symmetry and results in a position-dependent HartreeFock potential. In this case its contribution to the single-electron energies depends on the relevant quantum numbers. We numerically solved the Bogoliubovde Gennes equations with the HartreeFock term for a clean superconducting nanocylinder and found a shift of the curve representing the thickness-dependent oscillations of the critical superconducting temperature to larger diameters.

Abstract: We present an investigation of the orientations and positions of tubular fullerene molecules (C90, ..., C200) encapsulated in single-walled carbon nanotubes (SWCNT), a series of so-called fullerene nanopeapods. We find that increasing the tube radius leads to the following succession of energetically stable regimes: (1) lying molecules positioned on the tube's long axis; (2) tilted molecules on the tube's long axis; and (3) lying molecules shifted away from the tube's long axis. As opposed to C70 and C80 molecules encapsulated in a SWCNT, standing orientations do not develop. Our results are relevant for the possible application of molecular-orientation-dependent electronic properties of fullerene nanopeapods, and also for the interpretation of future experiments on double-walled carbon nanotube formation by annealing fullerene peapod systems.

Abstract: Apart from the most studied tin-oxide compounds, SnO and SnO2, intermediate states have been claimed to exist for more than a hundred years. In addition to the known homologous series (Seko et al., Phys. Rev. Lett. 100, 045702 (2008)), we here predict the existence of several new compounds with an O concentration between 50 % (SnO) and 67 % (SnO2). All these intermediate compounds are constructed from removing one or more (101) oxygen layers of SnO2. Since the van der Waals (vdW) interaction is known to be important for the Sn-Sn interlayer distances, we use a vdW-corrected functional, and compare these results with results obtained with PBE and hybrid functionals. We present the electronic properties of the intermediate structures and we observe a decrease of the band gap when (i) the O concentration increases and (ii) more SnO-like units are present for a given concentration. The contribution of the different atoms to the valence and conduction band is also investigated.

Abstract: Chalcopyrite semiconductors are of considerable interest for application as absorber layers in thin-film photovoltaic cells. When growing films of these compounds, however, they are often found to contain CuAu-like domains, a metastable phase of chalcopyrite. It has been reported that for CuInS2, the presence of the CuAu-like phase improves the short circuit current of the chalcopyrite-based photovoltaic cell. We investigate the thermodynamic stability of both phases for a selected list of I-III-VI2 materials using a first-principles density functional theory approach. For the CuIn-VI2 compounds, the difference in formation energy between the chalcopyrite and CuAu-like phase is found to be close to 2 meV per atom, indicating a high likelihood of the presence of CuAu-like domains. Next, we calculate the spectroscopic limited maximum efficiency (SLME) of the CuAu-like phase and compare the results with those of the corresponding chalcopyrite phase. We identify several candidates with a high efficiency, such as CuAu-like CuInS2, for which we obtain an SLME of 29% at a thickness of 500 nm. We observe that the SLME can have values above the Shockley-Queisser (SQ) limit, and show that this can occur because the SQ limit assumes the absorptivity to be a step function, thus overestimating the radiative recombination in the detailed balance approach. This means that it is possible to find higher theoretical efficiencies within this framework simply by calculating the J-V characteristic with an absorption spectrum. Finally, we expand our SLME analysis to indirect band gap absorbers by studying silicon, and find that the SLME quickly overestimates the reverse saturation current of indirect band gap materials, drastically lowering their calculated efficiency.

Abstract: We investigate the role of transition metal atoms of group V-b (V, Nb, Ta) and VI-b (Cr, Mo, W) as n- or p-type dopants in anatase TiO$_2$ using thermodynamic
principles and density functional theory with the Heyd-Scuseria-Ernzerhof HSE06 hybrid functional. The HSE06 functional provides a realistic value for the band gap, which ensures a correct classification of dopants as shallow or deep donors or acceptors. Defect formation energies and thermodynamic transition levels are calculated taking into account the constraints imposed by the stability of TiO$_2$ and the solubility limit of the impurities.
Nb, Ta, W and Mo are identified as shallow donors. Although W provides two electrons, Nb and Ta show a considerable lower formation energy, in particular under O-poor conditions. Mo donates in principle one electron, but under specific conditions can turn into a double donor. V impurities are deep donors and Cr
shows up as an amphoteric defect, thereby acting as an electron trapping center in n-type TiO$_2$ especially under O-rich conditions. A comparison with the available experimental data yields excellent agreement.

Abstract: We report on a first-principles study of the structural and electronic properties of a Sigma3 (112) grain boundary model in CuInSe2. The study focuses on a coherent, stoichiometry preserving, cation–Se terminated grain boundary, addressing the properties of the grain boundary as such, as well as the effect
of well known defects in CuInSe2. We show that in spite of its apparent simplicity, such a grain boundary exhibits a very rich phenomenology, providing an explanation for several of the experimentally observed properties of grain boundaries in CuInSe2 thin films. In particular, we show that the combined effect of Cu vacancies and cation antisites can result in the observed Cu depletion with no In enrichment at the grain boundaries. Furthermore, Cu vacancies are unlikely to produce a hole barrier at the grain boundaries, but Na may indeed have such an effect. We find that Na-on-Cu defects will tend to form abundantly at
the grain boundaries, and can provide a mechanism for the carrier depletion and/or type inversion experimentally reported.

Abstract: When slow ions incident on a surface are neutralized, the excess potential energy is passed on to an electron inside the surface, leading to emission of secondary electrons. The microscopic description of this process, as
well as the calculation of the secondary electron yield, is a challenging problem due to its complexity as well
as its sensitivity to surface properties. One of the first quantitative descriptions was articulated in the 1950s by
Hagstrum, who based his calculation on a parametrization of the density of states of the material. In this paper, we
present a model for calculating the secondary electron yield, derived from Hagstrum’s initial approach. We use
first-principles density functional theory calculations to acquire the necessary input and introduce the concept of
electron cascades to Hagstrum’s model in order to improve the calculated spectra, as well as remove its reliance
on fitting parameters. We apply our model to He+ and Ne+ ions incident on Ge(111) and Si(111) and obtain
yield spectra that match closely to the experimental results of Hagstrum.

Abstract: Lithium-rich layered oxides such as Li₂MnO₃have shown great potential as cathodes in Li-ion batteries, mainly because of their large capacities. However, these materials still suffer from structural degradation as the battery is cycled, reducing the average voltage and capacity of the cell. The voltage fade is believed to be related to the migration of transition metals into the lithium layer, linked to the formation of O-O dimers with a short bond length, which in turn is driven by the presence of oxygen holes due to the participation of oxygen in the redox process. We investigate the formation of O-O dimers for partially charged O1-Li₂MnO₃using a first-principles density functional theory approach by calculating the reaction energy and kinetic barriers for dimer formation. Next, we perform similar calculations for partially charged O1-Li₂IrO₃, a Li-rich material for which the voltage fade was not observed during cycling. When we compare the stability of the oxygen framework, we conclude that the formation of O-O dimers is both thermodynamically and kinetically viable for O1-Li_0.5MnO₃. For O1-Li_0.5IrO₃, we observe that the oxygen lattice is much more stable, either returning to its original state when perturbed, or resulting in a structure with an O-O dimer that is much higher in energy. This can be explained by the mixed redox process for Li₂IrO₃, which is also shown from the calculated magnetic moments. The lack of O-O dimer formation in O1-Li_0.5IrO₃provides valuable insight as to why Li₂IrO₃does not demonstrate a voltage fade as the battery is cycled, which can be used to design Li-rich battery cathodes with an improved cycling performance.

Abstract: We consider an interface between two strong time-reversal invariant topological insulators having surface states with opposite spin chirality, or equivalently, opposite mirror Chern number. We show that such an interface supports gapless modes that are protected by mirror symmetry. The interface states are investigated with a continuum model for the Bi2Se3 class of topological insulators that takes into account terms up to third order in the crystal momentum, which ensures that the model has the correct symmetry. The model parameters are obtained from ab initio calculations. Finally, we consider the effect of rotational mismatch at the interface, which breaks the mirror symmetry and opens a gap in the interface spectrum.

Abstract: The strain fields of periodically buckled graphene induce a periodic pseudomagnetic field (PMF) that modifies the electronic band structure. From the geometry, amplitude, and period of the periodic pseudomagnetic field, we determine the necessary conditions to access the regime of correlated phases by examining the band flattening. As compared to twisted bilayer graphene the proposed system has the advantages that (1) only a single layer of graphene is needed, (2) one is not limited to hexagonal superlattices, and (3) narrower flat bandwidth and larger separation between flat bands can be induced. We, therefore, propose that periodically strained graphene single layers can become a platform for the exploration of exotic many-body phases.

Abstract: The moire pattern observed in stacked noncommensurate crystal lattices, such as heterobilayers of transition metal dichalcogenides, produces a periodic modulation of their band gap. Excitons subjected to this potential landscape exhibit a band structure that gives rise to a quasiparticle dubbed the moire exciton. In the case of MoS2/WSe2 heterobilayers, the moire trapping potential has honeycomb symmetry and, consequently, the moire exciton band structure is the same as that of a Dirac-Weyl fermion, whose mass can be further tuned down to zero with a perpendicularly applied field. Here we show that, analogously to other Dirac-like particles, the moire exciton exhibits a trembling motion, also known as Zitterbewegung, whose long timescales are compatible with current experimental techniques for exciton dynamics. This promotes the study of the dynamics of moire excitons in van der Waals heterostructures as an advantageous solid-state platform to probe Zitterbewegung, broadly tunable by gating and interlayer twist angle.

Abstract: A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasiparticles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low-energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition-metal dichalcogenides. However, the method is easily adapted for other cases such as electrons in few-layer black phosphorus or any Schrodinger quasiparticles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in mono-, bi-, and trilayer graphene, where specular-and retro-reflection cases are observed for electrons scattered by a steplike superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.