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: 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: The valley splitting of the conduction-band states in a thin silicon-on-insulator layer is investigated using the 30-band k . p theory. The system composed of a few nm thick Si layer embedded within thick SiO2 layers is analyzed. The valley split states are found to cross periodically with increasing quantum well width, and therefore the energy splitting is an oscillatory function of the quantum well width, with period determined by the wave vector K-0 of the conduction band minimum. Because the valley split states are classified by parity, the optical transition between the ground hole state and one of those valley split conduction band states is forbidden. The oscillations in the valley splitting energy decrease with electric field and with smoothing of the composition profile between the well and the barrier by diffusion of oxygen from the SiO2 layers to the Si quantum well. Such a smoothing also leads to a decrease of the interband transition matrix elements. The obtained results are well parametrized by the effective two-valley model, but are found to disagree from previous 30-band calculations. This discrepancy could be traced back to the fact that the basis for the numerical solution of the eigenproblem must be restricted to the first Brillouin zone in order to obtain quantitatively correct results for the valley splitting.

Abstract: The electronic band alignment in semiconductor heterostructures is a key factor for their use in electronic applications. The alignment problem has been intensively studied for bulk systems but is less well understood for low-dimensional heterostructures. In this work we investigate the alignment in two-dimensional lateral heterostructures. First-principles calculations are used to show that the electronic band offset depends crucially on the width and thickness of the heterostructure slab. The particular heterostructures under study consist of thin hydrogenated and fluorinated diamond slabs which are laterally joined together. Two different limits for the band offset are observed. For infinitely wide heterostructures the vacuum potential above the two materials is aligned leading to a large step potential within the heterostructure. For infinitely thick heterostructure slabs, on the other hand, there is no potential step in the heterostructure bulk, but a large potential step in the vacuum region above the heterojunction is observed. The band alignment in finite systems depends on the particular dimensions of the system. These observations are shown to result from an interface dipole at the heterojunction that tends to align the band structures.

Abstract: Fulfillment of the promise of transparent electronics has been hindered until now largely by the lack of semiconductors that can be doped p-type in a stable way, and that at the same time present high hole mobility and are highly transparent in the visible spectrum. Here, a high-throughput study based on first-principles methods reveals four oxides, namely X2SeO2, with X = La, Pr, Nd, and Gd, which are unique in that they exhibit excellent characteristics for transparent electronic device applications – i.e., a direct band gap larger than 3.1 eV, an average hole effective mass below the electron rest mass, and good p-type dopability. Furthermore, for La2SeO2 it is explicitly shown that Na impurities substituting La are shallow acceptors in moderate to strong anion-rich growth conditions, with low formation energy, and that they will not be compensated by anion vacancies VO or VSe.

Abstract: We perform first-principles calculations to investigate the preferred positions of B and P dopants, both neutral and in their preferred charge state, in Si and Si/SiO2 core-shell nanowires (NWs). In order to understand the observed trends in the formation energy, we isolate the different effects that determine these formation energies. By making the distinction between the unrelaxed and the relaxed formation energy, we separate the impact of the relaxation from that of the chemical environment. The unrelaxed formation energies are determined by three effects: (i) the effect of strain caused by size mismatch between the dopant and the host atoms, (ii) the local position of the band edges, and (iii) a screening effect. In the case of the SiNW (Si/SiO2 NW), these effects result in an increase of the formation energy away from the center (interface). The effect of relaxation depends on the relative size mismatch between the dopant and host atoms. A large size mismatch causes substantial relaxation that reduces the formation energy considerably, with the relaxation being more pronounced towards the edge of the wires. These effects explain the surface segregation of the B dopants in a SiNW, since the atomic relaxation induces a continuous drop of the formation energy towards the edge. However, for the P dopants, the formation energy starts to rise when moving from the center but drops to a minimum just next to the surface, indicating a different type of behavior. It also explains that the preferential location for B dopants in Si/SiO2 core-shell NWs is inside the oxide shell just next to the interface, whereas the P dopants prefer the positions next to the interface inside the Si core, which is in agreement with recent experiments. These preferred locations have an important impact on the electronic properties of these core-shell NWs. Our simulations indicate the possibility of hole gas formation when B segregates into the oxide shell.

Abstract: Based on first-principles calculations, we study the structural stability of the H and T phases of monolayer MoS2 upon Li doping. Our calculations demonstrate that it is possible to stabilize a distorted T phase of MoS2 over the H phase through adsorption of Li atoms on the MoS2 surface. Through molecular dynamics and phonon calculations, we show that the T phase of MoS2 is dynamically unstable and undergoes considerable distortions. The type of distortion depends on the concentration of adsorbed Li atoms and changes from zigzag-like to diamond-like when increasing the Li doping. There exists a substantial energy barrier to transform the stable H phase to the distorted T phases, which is considerably reduced by increasing the concentration of Li atoms. We show that it is necessary that the Li atoms adsorb on both sides of the MoS2 monolayer to reduce the barrier sufficiently. Two processes are examined that allow for such two-sided adsorption, namely, penetration through the MoS2 layer and diffusion over the MoS2 surface. We show that while there is only a small barrier of 0.24 eV for surface diffusion, the amount of energy needed to pass through a pure MoS2 layer is of the order of similar or equal to 2 eV. However, when the MoS2 layer is covered with Li atoms the amount of energy that Li atoms should gain to penetrate the layer is drastically reduced and penetration becomes feasible.

Abstract: A thorough investigation of three-dimensional kagome lattices of group IV elements is performed with first-principles calculations. The investigated kagome lattices of silicon and germanium are found to be of similar stability as the recently proposed carbon kagome lattice. Carbon and silicon kagome lattices are both direct-gap semiconductors but they have qualitatively different electronic band structures. While direct optical transitions between the valence and conduction bands are allowed in the carbon case, no such transitions can be observed for silicon. The kagome lattice of germanium exhibits semimetallic behavior but can be transformed into a semiconductor after compression.

Abstract: A first-principles study of the structural and electronic properties of carbon impurities in CuInSe2 and CuGaSe2 is presented. Carbon is present in organic molecules in the precursor solutions used in non-vacuum growth methods for CuInSe2 and CuGaSe2 based photovoltaic cells. These growth methods make more efficient use of material, time, and energy than traditional vacuum methods. The formation energies of several carbon impurities are calculated using the hybrid HSE06 functional. C Cu acts as a shallow donor, CIn and interstitial C yield deep donor levels in CuInSe2, while in CuGaSe2 CGa and interstitial C act as deep amphoteric defects. So, these defects reduce the majority carrier (hole) concentration in p-type CuInSe2 and CuGaSe2 by compensating the acceptor levels. The deep defects are likely to act as recombination centers for the photogenerated charge carriers and are thus detrimental for the performance of the photovoltaic cells. On the other hand, the formation energies of the carbon impurities are high, even under C-rich growth conditions. Thus, few C impurities will form in CuInSe2 and CuGaSe2 in thermodynamic equilibrium. However, the deposition of the precursor solution in non-vacuum growth methods presents conditions far from thermodynamic equilibrium. In this case, our calculations show that C impurities formed in non-equilibrium tend to segregate from CuInSe2 and CuGaSe2 by approaching thermodynamic equilibrium, e.g., via thorough annealing.

Abstract: We study the adsorption and desorption kinetics of interacting particles moving on a one-dimensional lattice. Confinement is introduced by limiting the number of particles on a lattice site. Adsorption and desorption are found to proceed at different rates, and are strongly influenced by the concentration-dependent transport diffusion. Analytical solutions for the transport and self-diffusion are given for systems of length 1 and 2 and for a zero-range process. In the last situation the self- and transport diffusion can be calculated analytically for any length.

Abstract: Bi2Se3 is a three-dimensional topological insulator which has been extensively studied because it has a single Dirac cone on the surface, inside a relatively large bulk band gap. However, the effect of two-dimensional topological insulator Bi bilayers on the properties of Bi2Se3 and vice versa, has not been explored much. Bi bilayers are often present between the quintuple layers of Bi2Se3, since (Bi2)n(Bi2Se3)m form stable ground-state structures. Moreover, Bi2Se3 is a good substrate for growing ultrathin Bi bilayers. By first-principles techniques, we first show that there is no preferable surface termination by either Bi or Se. Next, we investigate the electronic structure of Bi bilayers on top of, or inside a Bi2Se3 slab. If the Bi bilayers are on top, we observe a charge transfer to the quintuple layers that increases the binding energy of the surface Dirac cones. The extra states, originating from the Bi bilayers, were declared to form a topological Dirac cone, but here we show that these are ordinary Rashba-split states. This result, together with the appearance of a new Dirac cone that is localized slightly deeper, might necessitate the reinterpretation of several experimental results. When the Bi bilayers are located inside the Bi2Se3 slab, they tend to split the slab into two topological insulators with clear surface states. Interface states can also be observed, but an energy gap persists because of strong coupling between the neighboring quintuple layers and the Bi bilayers.

Abstract: We present a high-throughput computing scheme based on density functional theory (DFT) to generate a class of oxides and screen them with the aim of identifying those that might be electronically appropriate for transparent conducting oxide (TCO) applications. The screening criteria used are a minimum band gap to ensure sufficient transparency, a band edge alignment consistent with easy n- or p-type dopability, and a minimum thermodynamic phase stability to be experimentally synthesizable. Following this scheme we screened 23 binary and 1518 ternary bixbyite oxides in order to identify promising candidates, which can then be a subject of an in-depth study. The results for the known TCOs are in good agreement with the reported data in the literature. We suggest a list of several new potential TCOs, including both n- and p-type compounds.

Abstract: Positrons can be trapped in localized states at the surface of a material, and thus quite selectively interact with core or valence surface electrons. Hence, advanced surface positron spectroscopy techniques can present the ideal tools to study a topological insulator, where surface states play a fundamental role. We analyze the problem of a positron at a TI surface, assuming that it is a weakly physisorbed positronium (Ps) atom. To determine if the surface of interest in a material can sustain such a physisorption, an accurate description of the underlying van der Waals (vdW) interaction is essential. We have developed a first-principles parameterfree method, based on the density functional theory, to extract key parameters determining the vdW interaction potential between a Ps atom and the surface of a given material. The method has been successfully applied to quartz and preliminary results on Bi2Te2Se indicate the existence of a positron surface state. We discuss the robustness of our predictions versus the most relevant approximations involved in our approach.

Abstract: By considering a cylindrically symmetric generalization of a plane wave, the first-order Born approximation of screened Coulomb scattering unfolds two new dimensions in the scattering problem: transverse momentum and orbital angular momentum of the incoming beam. In this paper, the elastic Coulomb scattering amplitude is calculated analytically for incoming Bessel beams. This reveals novel features occurring for wide-angle scattering and quantitative insights for small-angle vortex scattering. The result successfully generalizes the well-known Rutherford formula, incorporating transverse and orbital angular momentum into the formalism.

Abstract: The electronic and optical properties of zincblende nanowires are investigated in the presence of a uniform magnetic field directed along the [001] growth direction within the k . p method. We focus our numerical study on core-shell nanowires consisting of the III-V materials GaAs, AlxGa1-xAs and AlyGa1-y/0.51In0.49P. Nanowires with electrons confined in the core exhibit a Fock-Darwin-like spectrum, whereas nanowires with electrons confined in the shell show Aharonov-Bohm oscillations. Thus, by properly choosing the core and the shell materials of the nanowire, the optical properties in a magnetic field can be tuned in very different ways.