We present a first-principles computational study of cation–Se Σ3 (112) grain boundaries in CuGaSe2. We discuss the structure of these grain boundaries, as well as the effect of native defects and Na impurities on their electronic properties. The formation energies show that the defects will tend to form preferentially at the grain boundaries, rather than in the grain interiors. We find that in Ga-rich growth conditions Cu vacancies as well as Ga at Cu and Cu at Ga antisites are mainly responsible for having the equilibrium Fermi level pinned toward the middle of the gap, resulting in carrier depletion. The Na at Cu impurity in its +1 charge state contributes to this. In Ga-poor growth conditions, on the other hand, the formation energies of Cu vacancies and Ga at Cu antisites are comparatively too high for any significant influence on carrier density or on the equilibrium Fermi level position. Thus, under these conditions, the Cu at Ga antisites give rise to a p-type grain boundary. Also, their formation energy is lower than the formation energy of Na at Cu impurities. Thus, the latter will fail to act as a hole barrier preventing recombination at the grain boundary, in contrast to what occurs in CuInSe2 grain boundaries. We also discuss the effect of the defects on the electronic properties of bulk CuGaSe2, which we assume reflect the properties of the grain interiors.

We present a first-principles computational study of CO and OH adsorption on non-polar ZnO (101¯0) surfaces doped with indium. The calculations were performed using a model ZnO slab. The position of the In dopants was varied from deep bulk-like layers to the surface layers. It was established that the preferential location of the In atoms is at the surface by examining the dependence of the defect formation energy as well as the surface energy on In location. The adsorption sites on the surface of ZnO and the energy of adsorption of CO molecules and OH-species were determined in connection to In doping. It was found that OH has higher bonding energy to the surface than CO. The presence of In atoms at the surface of ZnO is favorable for CO adsorption, resulting in an elongation of the C-O bond and in charge transfer to the surface. The effect of CO and OH adsorption on the electronic and conduction properties of surfaces was assessed. We conclude that In-doped ZnO surfaces should present a higher electronic response upon adsorption of CO.

Effect of active sites of nanocrystalline tin(IV), tungsten(VI), and zinc-indium oxides on sensing behavior to CO and ammonia is analyzed.The oxides were synthesized by aqueous deposition routes, and subsequently modified by catalytic (Pd, Ru) and/or adsorptive (In, SO 4 ) additives.Surface sites and gas-solid interaction were investigated by HRTEM, EDX-mapping, FTIR, probe molecule methods, sensing test and DRIFT spectroscopy.CO sensitivity of ZnO(In) increased with the dopant amount, in agreement with DFT-predicted catalytic activity of In(III) in CO oxidation.The largest catalytic effect favoring CO detection down to 1 ppm at room temperature was shown for Pd-modified sensors.Gas-solid interaction route in such systems was controlled by CO chemisorption on Pdclusters irrespective of the supporting matrix.The sensitivity of SnO 2 /SO 4 and WO 3 to ammonia was improved by increasing surface acidity.The largest NH 3 sensitivity was observed for the materials with modified by Ru clusters, which provide catalytic oxidation of target molecules.Surface acidity in the latter case should not be too strong for not to deactivate the adsorbed ammonia species.

A combined experimental and first-principles study is performed to study the origin of conductivity in ZnO:Al nanoparticles synthesized under controlled conditions via a reflux route using benzylamine as a solvent.

We report on a first-principles study of the structural and electronic properties of a Σ3 (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.

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/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 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.

In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the u̇l-component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.

Cu-based chalcogenides are promising materials for thin-film solar cells with more than 20% measured cell efficiency. Using first-principles calculations based on density functional theory, the optoelectronic properties of a group of Cu-based chalcogenides Cu-2-II-IV-VI4 is studied. They are then screened with the aim of identifying potential absorber materials for photovoltaic applications. The spectroscopic limited maximum efficiency (SLME) introduced by Yu and Zunger [Phys. Rev. Lett. 108, 068701 (2012)] is used as a metric for the screening. After constructing the current-voltage curve, the SLME is calculated from the maximum power output. The role of the nature of the band gap, direct or indirect, and also of the absorptivity of the studied materials on the maximum theoretical power conversion efficiency is studied. Our results show that Cu2II-GeSe4 with II = Cd and Hg, and Cu-2-II-SnS4 with II = Cd, Hg, and Zn have a higher theoretical efficiency compared with the materials currently used as absorber layer. Published by AIP Publishing.

The motivation for the present study is to be found in the recent experiments of Fuchizaki and Hamaya on GeI4. They observed a rapid increase in the melting temperature Tm in going from atmospheric pressure to p 2: 6 GPa. Tm was found to be largely independent of pressure above this value. In this paper, heuristic arguments are presented to support the idea that until some critical pressure, a crystalline phase of SnI4, CCl4 and GeI4 molecular solids melts into a low density liquid. However, at this critical pressure, a phase boundary intersects Tmdp, separating a low density liquid phase from a high density liquid. The new phase boundary is between the crystal and an amorphous molecular solid with increasing polymerisation as the pressure is increased.

We present the results of our study on p-type dopability of bixbyite V2O3 using the Heyd, Scuseria, and Ernzerhof hybrid functional (HSE06) within the density functional theory (DFT) formalism. We study vanadium and oxygen vacancies as intrinsic defects and substitutional Mg, Sc, and Y as extrinsic defects. We find that Mg substituting V acts as a shallow acceptor, and that oxygen vacancies are electrically neutral. Hence, we predict Mg-doped V2O3 to be a p-type conductor. Our results also show that vanadium vacancies are relatively shallow, with a binding energy of 0.14 eV, so that they might also lead to p-type conductivity.

Past research efforts aiming at obtaining stable p-type ZnO have been based on complexes involving nitrogen doping. A recent experiment by (J. G. Reynolds et al., Appl. Phys. Lett., 2013, 102, 152114) demonstrated a significant (∼10(18) cm(-3)) p-type behavior in N-doped ZnO films after appropriate annealing. The p-type conductivity was attributed to a VZn-NO-H shallow acceptor complex, formed by a Zn vacancy (VZn), N substituting O (NO), and H interstitial (Hi). We present here a first-principles hybrid functional study of this complex compared to the one without hydrogen. Our results confirm that the VZn-NO-H complex acts as an acceptor in ZnO. We find that H plays an important role, because it lowers the formation energy of the complex with respect to VZn-NO, a complex known to exhibit (unstable) p-type behavior. However, this additional H atom also occupies the hole level at the origin of the shallow behavior of VZn-NO, leaving only two states empty higher in the band gap and making the VZn-NO-H complex a deep acceptor. Therefore, we conclude that the cause of the observed p-type conductivity in experiment is not the presence of the VZn-NO-H complex, but probably the formation of the VZn-NO complex during the annealing process.

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. CCu 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.

ZnM2O4 (M = Co, Rh, Ir) spinels are considered as a class of potential p-type transparent conducting oxides (TCOs). We report the formation energy of acceptor-like defects using first principles calculations with an advanced hybrid exchange-correlation functional (HSE06) within density functional theory (DFT). Due to the discrepancies between the theoretically obtained band gaps with this hybrid functional and the - scattered - experimental results, we also perform GW calculations to support the validity of the description of these spinels with the HSE06 functional. The considered defects are the cation vacancy and antisite defects, which are supposed to be the leading source of disorder in the spinel structures. We also discuss the band alignments in these spinels. The calculated formation energies indicate that the antisite defects ZnM (Zn replacing M, M = Co, Rh, Ir) and VZn act as shallow acceptors in ZnCo2O4, ZnRh2O4 and ZnIr2O4, which explains the experimentally observed p-type conductivity in those systems. Moreover, our systematic study indicates that the ZnIr antisite defect has the lowest formation energy in the group and it corroborates the highest p-type conductivity reported for ZnIr2O4 among the group of ZnM2O4 spinels. To gain further insight into factors affecting the p-type conductivity, we have also investigated the formation of localized small polarons by calculating the self-trapping energy of the holes.

In order to account explicitly for the existence of long-periodic layered structures and the strong structural relaxations in the most common binary and ternary alloys of the Bi-Sb-Te-Se system, we have developed a one-dimensional cluster expansion (CE) based on first-principles electronic structure calculations, which accounts for the Bi and Sb bilayer formation. Excellent interlayer distances are obtained with a van der Waals density functional. It is shown that a CE solely based on pair interactions is sufficient to provide an accurate description of the ground-state energies of Bi-Sb-Te-Se binary and ternary systems without making the data set of ab initio calculated structures unreasonably large. For the binary alloys A(1-x)Q(x) (A = Sb, Bi; Q = Te, Se), a ternary CE yields an almost continuous series of (meta) stable structures consisting of consecutive A bilayers next to consecutive A(2)Q(3) for 0 < x < 0.6. For x > 0.6, the binary alloy segregates into pure Q and A(2)Q(3). The Bi-Sb system is described by a quaternary CE and is found to be an ideal solid solution stabilized by entropic effects at T not equal 0 K but with an ordered structure of alternating Bi and Sb layers for x = 0.5 at T = 0 K. A quintuple CE is used for the ternary Bi-Sb-Te system, where stable ternary layered compounds with an arbitrary stacking of Sb2Te3, Bi2Te3, and Te-Bi-Te-Sb-Te quintuple units are found, optionally separated by mixed Bi/Sb bilayers. Electronic properties of the stable compounds were studied taking spin-orbit coupling into account.

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 parameter-free 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.

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 (Bi-2)(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.