A Theoretical Study on the Surfaces of Zinc Oxide

Zinc oxide is an important wide bandgap n-type semiconductor with uses ranging from electronics to catalysis. The chemical and physical properties related to its surfaces are of fundamental interest and also key to the material’s design. In this Thesis, computational methods have been used to model the surfaces of ZnO. We report a detailed theoretical study on the four main low-index wurtzite ZnO surfaces. For nonpolar surfaces, we focus on the stability, atomic structure and electronic properties of both clean and defective surfaces. Our calculations explain why steps are common on the (10-10) surface, as seen in experiment. We calculate the ionisation potential which is in good agreement with experiment. The electronic band edges of the nonpolar surfaces are seen to behave differently, with a local rise of the VBM and CBM for (10-10) and (11-20), respectively. For ZnO polar surfaces, our results can explain why experimental findings reported have been varied and even contradictory at times. The calculated surface energies indicate on average a slightly higher stability of the (000-1) surface compared to the (0001) surface. Structurally, triangular and hexagonal patterns are seen among the stable structures but a high level of disorder is predicted. We also report new interatomic potentials (IP) for the Cu/ZnO system. Our IP can work as a fast and reliable method to filter low energy Cu/ZnO structures. Global optimisation calculations show a preference for planar Cu clusters over the (10-10) surface, with a strong interaction between the Cu and Zn species. Finally, we study the surface atomic configurations for the MoO3/Fe2O3 catalytic system. The lowest energy structure was used in the fitting of EXAFS parameters. Overall, our Thesis shows the great utility of theoretical calculations in the explanation of experimental findings in surface science

Disorder-induced electron and hole trapping in amorphous TiO2

Thin films of amorphous (a)-TiO2 are ubiquitous as photocatalysts, protective coatings, photoanodes and in memory application, where they are exposed to excess electrons and holes. We investigate trapping of excess electrons and holes in the bulk of pure amorphous titanium dioxide, a-TiO2, using hybrid density functional theory (h-DFT) calculations. Fifty 270-atom a-TiO2 structures were produced using classical molecular dynamics and their geometries fully optimised using h-DFT simulations. They have the density, distribution of atomic coordination numbers and radial pair-distribution functions in agreement with experiment. The calculated average a-TiO2 band gap is 3.25 eV with no states splitting into the band gap. Trapping of excess electrons and holes in a-TiO2 is predicted at precursor sites, such as elongated Ti-O bonds. Single electron and hole polarons have average trapping energies (ET) of -0.4 eV and -0.8 eV, respectively. We also identify several types of electron and hole bipolaron states and discuss their stability. These results can be used for understanding the mechanisms of photo-catalysis and improving the performance of electronic devices employing a-TiO2 films.Comment: 12 pages, 10 figures. This article has been submitted to Physical Review

Making amorphous ZnO: Theoretical predictions of its structure and stability

ZnO is a transparent semiconductor with optoelectronic, thermoelectric, and sensor applications, where using amorphous thin films presents great advantages. However, growing amorphous (a) films of pure ZnO proved challenging due to their rapid crystallization. We investigated the ability of bulk ZnO to form glass structures using well-tested interatomic potentials and a melt and quench procedure within isothermal-isobaric (NPT) ensemble. The geometries of some of the resulting structures were further optimized using density functional theory (DFT) calculations with the PBE functional. We demonstrate that cooling rates in melt and quench procedure equal or exceeding 100 Kps−1 lead to formation of stable amorphous structures. However, ZnO samples tend to crystallize at lower cooling rates. This result does not depend on the size of the periodic cell used in the calculations for cells containing more than 324 atoms. Using simulation cells with up to 768 000 atoms, we demonstrate that the expected average glass density is about 5.04 gcm−3 and the coordination numbers of Zn and O atoms are around 3.9. We calculate radial distribution functions and characterize the structures of amorphous ZnO samples. Using both the activation-relaxation technique and simulated annealing, we show that the obtained amorphous structures have low propensity to crystallization

Morphology of Cu clusters supported on reconstructed polar ZnO (0001) and (0001̄) surfaces

Unbiased Monte Carlo procedures are applied to investigate the structure of Cu clusters of various sizes deposited over reconstructed polar ZnO surfaces. Four distinct reconstructed polar ZnO surfaces (two Zn terminated (0001) reconstructions and two O terminated (000[1 with combining macron]) reconstructions) were investigated, having previously been determined to be the most stable under typical conditions, as revealed by the grand canonical ensemble studies. Random sampling was performed considering ∼400 000 random initial structural configurations of Cu atoms over the ZnO surfaces, with each structure being optimised using interatomic potential techniques, and the most stable resultant structures being refined using a plane-wave DFT approach. The investigation reveals the key role of surface adatoms and vacancies arising from the reconstruction of the polar ZnO surface in determining the morphology of deposited Cu clusters. Strong Cu–Zn interactions play an essential role in Cu cluster growth, with reconstructed polar ZnO surfaces featuring sites with undercoordinated Zn surface atoms promoting highly localised three dimensional Cu cluster morphologies, whist reconstructions featuring undercoordinated O atoms tend to result in more planar Cu clusters, in order to maximise the favourable Cu–Zn interaction. This is the first study that evaluates the thermodynamically most stable Cu/ZnO structures using realistic reconstructed ZnO polar surfaces, and thus provides valuable insights into the factors affecting Cu cluster growth over ZnO surfaces, as well as model catalyst surfaces that can be utilised in future computational studies to explore catalytic activity for key processes such as CO2 and CO hydrogenation to methanol

Development of Interatomic Potentials for Supported Nanoparticles: The Cu/ZnO Case

We present a potential model that has been parametrized to reproduce accurately metal−metal oxide interactions of Cu clusters supported on ZnO. Copper deposited on the nonpolar (101̅0) ZnO surface is investigated using the new pairwise Cu–ZnO interatomic potentials including repulsive Born–Mayer Cu–O and attractive Morse Cu–Zn potentials. Parameters of these interactions have been determined by fitting to periodic supercell DFT data using different surface terminations and Cu cluster sizes. Results of interatomic potential-based simulations show a good agreement both structurally and energetically with DFT data, and thus provide an efficient filter of configurations during a search for low DFT energy structures. Upon examining the low energy configurations of Cu clusters on ZnO nonpolar surfaces for a range of cluster sizes, we discovered why Cu islands are commonly observed on step edges on the (101̅0) surface but are rarely seen on terraces

Inverse Perovskite Oxysilicides and Oxygermanides as Candidates for Nontoxic Infrared Semiconductor and Their Chemical Bonding Nature

We have synthesized inverse-perovskite-type oxysilicides and oxygermanides represented by R3SiO and R3GeO (R = Ca and Sr) and studied their characteristics in the search for nontoxic narrow band gap semiconductors. These compounds exhibit a sharp absorption edge around 0.9 eV and a luminescence peak in the same energy range. These results indicate that the obtained materials have a direct-band electronic structure, which was confirmed by hybrid DFT calculations. These materials, made from earth abundant and nontoxic elements and with a relatively light electron/hole effective mass, represent strong candidates for nontoxic optoelectronic devices in the infrared range

Real and virtual polymorphism of titanium selenide with robust interatomic potentials

The first successful pairwise potential for a layered material, TiSe2, has been parameterised to fit the experimental data, using a genetic algorithm as the optimisation tool for the parameters of the interatomic potential

Structure and Migration Mechanisms of Small Vacancy Clusters in Cu: A Combined EAM and DFT Study

Voids in face-centered cubic (fcc) metals are commonly assumed to form via the aggregation of vacancies; however, the mechanisms of vacancy clustering and diffusion are not fully understood. In this study, we use computational modeling to provide a detailed insight into the structures and formation energies of primary vacancy clusters, mechanisms and barriers for their migration in bulk copper, and how these properties are affected at simple grain boundaries. The calculations were carried out using embedded atom method (EAM) potentials and density functional theory (DFT) and employed the site-occupation disorder code (SOD), the activation relaxation technique nouveau (ARTn) and the knowledge led master code (KLMC). We investigate stable structures and migration paths and barriers for clusters of up to six vacancies. The migration of vacancy clusters occurs via hops of individual constituent vacancies with di-vacancies having a significantly smaller migration barrier than mono-vacancies and other clusters. This barrier is further reduced when di-vacancies interact with grain boundaries. This interaction leads to the formation of self-interstitial atoms and introduces significant changes into the boundary structure. Tetra-, penta-, and hexa-vacancy clusters exhibit increasingly complex migration paths and higher barriers than smaller clusters. Finally, a direct comparison with the DFT results shows that EAM can accurately describe the vacancy-induced relaxation effects in the Cu bulk and in grain boundaries. Significant discrepancies between the two methods were found in structures with a higher number of low-coordinated atoms, such as penta-vacancies and di-vacancy absortion by grain boundary. These results will be useful for modeling the mechanisms of diffusion of complex defect structures and provide further insights into the structural evolution of metal films under thermal and mechanical stress

Morphological Features and Band Bending at Nonpolar Surfaces of ZnO

We employ hybrid density functional calculations to analyze the structure and stability of the (101̅0) and (112̅0) ZnO surfaces, confirming the relative stability of the two surfaces. We then examine morphological features, including steps, dimer vacancies, and grooves, at the main nonpolar ZnO surface using density functional methods. Calculations explain why steps are common on the (101̅0) surface even at room temperature, as seen in experiment. The surface structure established has been used to obtain the definitive ionization potential and electron affinity of ZnO in good agreement with experiment. The band bending across the surface is analyzed by the decomposition of the density of states for each atomic layer. The upward surface band bending at the (101̅0) surface affects mostly the valence band by 0.32 eV, which results in the surface band gap closing by 0.31 eV; at the (112̅0) surface, the valence band remains flat and the conduction band bends up by 0.18 eV opening the surface band gap by 0.12 eV

Morphology of Cu clusters supported on reconstructed polar ZnO (0001) and (000[1]) surfaces†

Unbiased Monte Carlo procedures are applied to investigate the structure of Cu clusters of various sizes deposited over reconstructed polar ZnO surfaces. Four distinct reconstructed polar ZnO surfaces (two Zn terminated (0001) reconstructions and two O terminated (000[1 with combining macron]) reconstructions) were investigated, having previously been determined to be the most stable under typical conditions, as revealed by the grand canonical ensemble studies. Random sampling was performed considering ∼400 000 random initial structural configurations of Cu atoms over the ZnO surfaces, with each structure being optimised using interatomic potential techniques, and the most stable resultant structures being refined using a plane-wave DFT approach. The investigation reveals the key role of surface adatoms and vacancies arising from the reconstruction of the polar ZnO surface in determining the morphology of deposited Cu clusters. Strong Cu–Zn interactions play an essential role in Cu cluster growth, with reconstructed polar ZnO surfaces featuring sites with undercoordinated Zn surface atoms promoting highly localised three dimensional Cu cluster morphologies, whist reconstructions featuring undercoordinated O atoms tend to result in more planar Cu clusters, in order to maximise the favourable Cu–Zn interaction. This is the first study that evaluates the thermodynamically most stable Cu/ZnO structures using realistic reconstructed ZnO polar surfaces, and thus provides valuable insights into the factors affecting Cu cluster growth over ZnO surfaces, as well as model catalyst surfaces that can be utilised in future computational studies to explore catalytic activity for key processes such as CO2 and CO hydrogenation to methanol