Abstract: The physical properties of uranium dioxide vary greatly with stoichiometry. Oxidation towards hyperstoichiometric UO2 – UO2+x – might be encountered at various stages of the nuclear fuel cycle if oxidative conditions are met; the impact of stoichiometry changes upon physical properties should therefore be properly assessed to ensure safe and reliable operations. These physical properties are intimately linked to the arrangement of atomic defects in the crystalline structure. The evolution of the defect concentration with environmental parameters – oxygen partial pressure and temperature – were evaluated by means of a point defect model where the reaction energies are derived from atomic-scale simulations. To this end, various configurations and net charge states of oxygen interstitial clusters in UO2 have been calculated. Various methodologies have been tested to determine the optimum cluster configurations and a rigid lattice approach turned out to be the most useful strategy to optimize defect configuration structures. Ultimately, results from the point defect model were discussed and compared to experimental measurements of stoichiometry dependence on oxygen partial pressure and temperature.

Abstract: Although recent studies have advanced the understanding of pyrolyzed
Fe−N−C materials as oxygen reduction reaction (ORR) catalysts, the atomic and
electronic structures of the active sites and their detailed reaction mechanisms still remain unknown. Here, based on first-principles density functional theory (DFT) computations, we discuss the electronic structures of three FeN4 catalytic centers with different local topologies of the surrounding C atoms with a focus on unraveling the mechanism of their ORR activity in acidic electrolytes. Our study brings back a forgotten, synthesized pyridinic Fe−N coordinate to the community’s attention, demonstrating that this catalyst can exhibit excellent activity for promoting direct four-electron ORR through the addition of a fifth ligand such as −NH2, −OH, and −SO4. We also identify sites with good stability properties through the combined use of our DFT calculations and Mössbauer spectroscopy data.

Abstract: The atomic scattering amplitudes of the various atoms of the systems Ga1−xInxAs, GaAs1−xNx and InAs1−xNx are calculated using the density functional theory (DFT) approach. The scattering amplitudes of N, Ga, As and In in the model systems are compared with the frequently used Doyle and Turner values. Deviation from the latter values is found for small scattering vectors (s<0.3Å−1) and for these scattering vectors dependence on the orientation of the scattering vector and the chemical environment is reported. We suggest a parametrization of these modified scattering amplitudes (MASAs) for small scattering vectors (s<1.0Å−1). The MASAs are exploited within zero pressure classical Metropolis Monte Carlo (MC), finite temperature calculations to investigate the effect of quantum well size on the electron {002} structure factor (SF) of Ga1−xInxAs quantum wells.

Abstract: We calculated the temperature dependence of the Debye-Waller factors for a variety of group IV, III-V and II-VI semiconductors from 0.1 to 1000 K. The approach used to fit the temperature dependence is described and resulting fit parameters are tabulated for each material. The Debye-Waller factors are deduced from generalized phonon densities of states which were derived from first principles using the WIEN2k and the ABINIT codes.

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: Aberration-corrected transmission electron microscopy, electron energy-loss spectroscopy, and density functional theory (DFT) calculations are used to solve several key questions about the surface structure, the particle morphology, and the distribution and nature of nitrogen impurities in detonation nanodiamond (DND) cleaned by a recently developed ozone treatment. All microscopy and spectroscopy measurements are performed at a lowered acceleration voltage (80/120kV), allowing prolonged and detailed experiments to be carried out while minimizing the risk of knock-on damage or surface graphitization of the nanodiamond. High-resolution TEM (HRTEM) demonstrates the stability of even the smallest nanodiamonds under electron illumination at low voltage and is used to image the surface structure of pristine DND. High resolution electron energy-loss spectroscopy (EELS) measurements on the fine structure of the carbon K-edge of nanodiamond demonstrate that the typical * pre-peak in fact consists of three sub-peaks that arise from the presence of, amongst others, minimal fullerene-like reconstructions at the nanoparticle surfaces and deviations from perfect sp(3) coordination at defects in the nanodiamonds. Spatially resolved EELS experiments evidence the presence of nitrogen within the core of DND particles. The nitrogen is present throughout the whole diamond core, and can be enriched at defect regions. By comparing the fine structure of the experimental nitrogen K-edge with calculated energy-loss near-edge structure (ELNES) spectra from DFT, the embedded nitrogen is most likely related to small amounts of single substitutional and/or A-center nitrogen, combined with larger nitrogen clusters.

Abstract: Topological insulators (TIs) are renowned for their exotic topological surface states (TSSs) that reside in the top atomic layers, and hence, detailed knowledge of the surface top atomic layers is of utmost importance. Here we present the remarkable morphology changes of Bi2Te3 surfaces, which have been freshly cleaved in air, upon subsequent systematic annealing in ultrahigh vacuum and the resulting effects on the local and area-averaging electronic properties of the surface states, which are investigated by combining scanning tunneling microscopy (STM), scanning tunneling spectroscopy (STS), and Auger electron spectroscopy (AES) experiments with density functional theory (DFT) calculations. Our findings demonstrate that the annealing induces the formation of a Bi bilayer atop the Bi2Te3 surface. The adlayer results in n-type doping, and the atomic defects act as scattering centers of the TSS electrons. We also investigated the annealing-induced Bi bilayer surface on Bi2Te3 via voltage-dependent quasi-particle-interference (QPI) mapping of the surface local density of states and via comparison with the calculated constant-energy contours and QPI patterns. We observed closed hexagonal patterns in the Fourier transform of real-space QPI maps with secondary outer spikes. DFT calculations attribute these complex QPI patterns to the appearance of a “second” cone due to the surface charge transfer between the Bi bilayer and the Bi2Te3. Annealing in ultrahigh vacuum offers a facile route for tuning of the topological properties and may yield similar results for other topological materials.

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: Thermal diffuse scattered electrons significantly contribute to high-resolution transmission electron microscopy images. Their intensity adds to the background and is peaked at positions of atomic columns. In this paper we suggest an approximation to simulate intensity of thermal diffuse scattered electrons in plane-wave illumination transmission electron microscopy using an emission-potential multislice algorithm which is computationally less intensive than the frozen lattice approximation or the mutual intensity approach. Intensity patterns are computed for Au and InSb for different crystal orientations. These results are compared with intensities from the frozen lattice approximation based on uncorrelated vibration of atoms as well as with the frozen phonon approximation for Au. The frozen phonon method uses a detailed phonon model based on force constants we computed by a density functional theory approach. The comparison shows that our suggested emission-potential method is in close agreement with both the frozen lattice and the frozen phonon approximations.

Abstract: We present the quasiparticle band structure of ZnO in its zincblende (ZB) and rocksalt (RS) phases at the Γ point, calculated within the GW approximation. The effect of the pd hybridization on the quasiparticle corrections to the band gap is discussed. We compare three systems, ZB-ZnO which shows strong pd hybridization and has a direct band gap, RS-ZnO which is also hybridized but includes inversion symmetry and therefore has an indirect band gap, and ZB-ZnS which shows a weaker hybridization due to a change of the chemical species from oxygen to sulfur. The quasiparticle corrections are calculated with different numbers of valence electrons in the Zn pseudopotential. We find that the Zn20 + pseudopotential is essential for the adequate treatment of the exchange interaction in the self-energy. The calculated GW band gaps are 2.47 eV and 4.27 eV respectively, for the ZB and RS phases. The ZB-ZnO band gap is underestimated compared to the experimental value of 3.27 by ~ 0.8 eV. The RS-ZnO band gap compares well with the experimental value of 4.5 eV. The underestimation for ZB-ZnO is correlated with the strong pd hybridization. The GW band gap for ZnS is 3.57 eV, compared to the experimental value of 3.8 eV.

Abstract: 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 Cu2-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 currentvoltage
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 Cu2-II-SnS4 with II ¼ Cd, Hg, and Zn have a higher theoretical efficiency
compared with the materials currently used as absorber layer.

Abstract: Layered Li-rich oxides, demonstrating both cationic and anionic redox chemistry being used as positive electrodes for Li-ion batteries,have raised interest due to their high specific discharge capacities exceeding 250 mAh/g. However, irreversible structural transformations triggered by anionic redox chemistry result in pronounced voltagefade (i.e., lowering the specific energy by a gradual decay of discharge potential) upon extended galvanostatic cycling. Activating or suppressing oxygen anionic redox through structural stabilization induced by redox-inactivecation substitution is a well-known strategy. However, less emphasishas been put on the correlation between substitution degree and theactivation/suppression of the anionic redox. In this work, Ti4+-substituted Li2MnO3 was synthesizedvia a facile solution-gel method. Ti4+ is selected as adopant as it contains no partially filled d-orbitals. Our study revealedthat the layered “honeycomb-ordered” C2/m structure is preserved when increasing the Ticontent to x = 0.2 in the Li2Mn1-x Ti (x) O-3 solidsolution, as shown by electron diffraction and aberration-correctedscanning transmission electron microscopy. Galvanostatic cycling hintsat a delayed oxygen release, due to an improved reversibility of theanionic redox, during the first 10 charge-discharge cyclesfor the x = 0.2 composition compared to the parentmaterial (x = 0), followed by pronounced oxygen redoxactivity afterward. The latter originates from a low activation energybarrier toward O-O dimer formation and Mn migration in Li2Mn0.8Ti0.2O3, as deducedfrom first-principles molecular dynamics (MD) simulations for the“charged” state. Upon lowering the Ti substitution to x = 0.05, the structural stability was drastically improvedbased on our MD analysis, stressing the importance of carefully optimizingthe substitution degree to achieve the best electrochemical performance.

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.