Engineering of metal oxide interfaces for renewable energy applications

Diminishing non-renewable energy resources and planet-wide de-pollution on our planet are among the major problems which mankind faces into the future. To solve these problems, renewable energy sources such as readily available and inexhaustible sunlight will have to be used. There are however no readily available photocatalysts that are photocatalytically active under visible light; it is well established that the band gap of the prototypical photocatalyst, titanium dioxide, is the UV region with the consequence that only 4% of sun light is utilized. For this reason, this PhD project focused on developing new materials, based on titanium dioxide, which can be used in visible light activated photocatalytic hydrogen production and destruction of pollutant molecules. The main goal of this project is to use simulations based on first principles to engineer and understand rationally, materials based on modifying TiO2 that will have the following properties: (1) a suitable band gap in order to increase the efficiency of visible light absorption, with a gap around 2 – 2.5 eV considered optimum. (2). The second key aspect in the photocatalytic process is electron and hole separation after photoexcitation, which enable oxidation/reduction reactions necessary to i.e. decompose pollutants. (3) Enhanced activity over unmodified TiO2. In this thesis I present results on new materials based on modifying TiO2 with supported metal oxide nanoclusters, from two classes, namely: transition metal oxides (Ti, Ni, Cu) and p-block metal oxides (Sn, Pb, Bi). We find that the deposited metal oxide nanoclusters are stable at rutile and anatase TiO2 surfaces and present an analysis of changes to the band gap of TiO2, identifying those modifiers that can change the band gap to the desirable range and the origin of this. A successful collaboration with experimental researchers in Japan confirms many of the simulation results where the origin of improved visible light photocatalytic activity of oxide nanocluster-modified TiO2 is now well understood. The work presented in this thesis, creates a road map for the design of materials with desired photocatalytic properties and contributes to better understanding these properties which are of great application in renewable energy utilization.

Computational modeling of defects in nanoscale device materials

This PhD thesis concerns the computational modeling of the electronic and atomic structure of point defects in technologically relevant materials. Identifying the atomistic origin of defects observed in the electrical characteristics of electronic devices has been a long-term goal of first-principles methods. First principles simulations are performed in this thesis, consisting of density functional theory (DFT) supplemented with many body perturbation theory (MBPT) methods, of native defects in bulk and slab models of In0.53Ga0.47As. The latter consist of (100) - oriented surfaces passivated with A12O3. Our results indicate that the experimentally extracted midgap interface state density (Dit) peaks are not the result of defects directly at the semiconductor/oxide interface, but originate from defects in a more bulk-like chemical environment. This conclusion is reached by considering the energy of charge transition levels for defects at the interface as a function of distance from the oxide. Our work provides insight into the types of defects responsible for the observed departure from ideal electrical behaviour in III-V metal-oxidesemiconductor (MOS) capacitors. In addition, the formation energetics and electron scattering properties of point defects in carbon nanotubes (CNTs) are studied using DFT in conjunction with Green’s function based techniques. The latter are applied to evaluate the low-temperature, low-bias Landauer conductance spectrum from which mesoscopic transport properties such as the elastic mean free path and localization length of technologically relevant CNT sizes can be estimated from computationally tractable CNT models. Our calculations show that at CNT diameters pertinent to interconnect applications, the 555777 divacancy defect results in increased scattering and hence higher electrical resistance for electron transport near the Fermi level.

Nonlinear piezoelectricity in electroelastic energy harvesters: Modeling and experimental identification

We propose and experimentally validate a first-principles based model for the nonlinear piezoelectric response of an electroelastic energy harvester. The analysis herein highlights the importance of modeling inherent piezoelectric nonlinearities that are not limited to higher order elastic effects but also include nonlinear coupling to a power harvesting circuit. Furthermore, a nonlinear damping mechanism is shown to accurately restrict the amplitude and bandwidth of the frequency response. The linear piezoelectric modeling framework widely accepted for theoretical investigations is demonstrated to be a weak presumption for near-resonant excitation amplitudes as low as 0.5 g in a prefabricated bimorph whose oscillation amplitudes remain geometrically linear for the full range of experimental tests performed (never exceeding 0.25% of the cantilever overhang length). Nonlinear coefficients are identified via a nonlinear least-squares optimization algorithm that utilizes an approximate analytic solution obtained by the method of harmonic balance. For lead zirconate titanate (PZT-5H), we obtained a fourth order elastic tensor component of c1111p =-3.6673× 1017 N/m2 and a fourth order electroelastic tensor value of e3111 =1.7212× 108 m/V. © 2010 American Institute of Physics.

Conductive junctions with parallel graphene sheets.

The establishment of conductive graphene-molecule-graphene junction is investigated through first-principles electronic structure calculations and quantum transport calculations. The junction consists of a conjugated molecule connecting two parallel graphene sheets. The effects of molecular electronic states, structure relaxation, and molecule-graphene contact on the conductance of the junction are explored. A conductance as large as 0.38 conductance quantum is found achievable with an appropriately oriented dithiophene bridge. This work elucidates the designing principles of promising nanoelectronic devices based on conductive graphene-molecule-graphene junctions.

Uncovering compounds by synergy of cluster expansion and high-throughput methods.

Predicting from first-principles calculations whether mixed metallic elements phase-separate or form ordered structures is a major challenge of current materials research. It can be partially addressed in cases where experiments suggest the underlying lattice is conserved, using cluster expansion (CE) and a variety of exhaustive evaluation or genetic search algorithms. Evolutionary algorithms have been recently introduced to search for stable off-lattice structures at fixed mixture compositions. The general off-lattice problem is still unsolved. We present an integrated approach of CE and high-throughput ab initio calculations (HT) applicable to the full range of compositions in binary systems where the constituent elements or the intermediate ordered structures have different lattice types. The HT method replaces the search algorithms by direct calculation of a moderate number of naturally occurring prototypes representing all crystal systems and guides CE calculations of derivative structures. This synergy achieves the precision of the CE and the guiding strengths of the HT. Its application to poorly characterized binary Hf systems, believed to be phase-separating, defines three classes of alloys where CE and HT complement each other to uncover new ordered structures.

Modeling non-harmonic behavior of materials from experimental inelastic neutron scattering and thermal expansion measurements

Based on thermodynamic principles, we derive expressions quantifying the non-harmonic vibrational behavior of materials, which are rigorous yet easily evaluated from experimentally available data for the thermal expansion coefficient and the phonon density of states. These experimentally- derived quantities are valuable to benchmark first-principles theoretical predictions of harmonic and non-harmonic thermal behaviors using perturbation theory, ab initio molecular-dynamics, or Monte-Carlo simulations. We illustrate this analysis by computing the harmonic, dilational, and anharmonic contributions to the entropy, internal energy, and free energy of elemental aluminum and the ordered compound FeSi over a wide range of temperature. Results agree well with previous data in the literature and provide an efficient approach to estimate anharmonic effects in materials.

Stability of Sr adatom model structures for SrTiO3(001) surface reconstructions

We report results of first-principles calculations on the thermodynamic stability of different Sr adatom structures that have been proposed to explain some of the observed reconstructions of the (001) surface of strontium titanate (Kubo and Nozoye 2003 Surf Sci. 542 177). From surface free energy calculations, a phase diagram is constructed indicating the range of conditions over which each structure is most stable. These results are compared with Kubo and Nozoye's experimental observations. It is concluded that low Sr adatom coverage structures can only be explained if the surface is far from equilibrium. Intermediate coverage structures are stable only if the surface is in or very nearly in equilibrium with the strontium oxide.

Bismuth embrittlement of copper is an atomic size effect

Embrittlement by the segregation of impurity elements to grain boundaries is one of a small number of phenomena that can lead to metallurgical failure by fast fracture(1). Here we settle a question that has been debated for over a hundred years(2): how can minute traces of bismuth in copper cause this ductile metal to fail in a brittle manner? Three hypotheses for Bi embrittlement of Cu exist: two assign an electronic effect to either a strengthening(3) or weakening(4) of bonds, the third postulates a simple atomic size effect(5). Here we report first principles quantum mechanical calculations that allow us to reject the electronic hypotheses, while supporting a size effect. We show that upon segregation to the grain boundary, the large Bi atoms weaken the interatomic bonding by pushing apart the Cu atoms at the interface. The resolution of the mechanism underlying grain boundary weakening should be relevant for all cases of embrittlement by oversize impurities.

A stabilization mechanism of zirconia based on oxygen vacancies only

The microscopic mechanism leading to stabilization of cubic and tetragonal forms of zirconia (ZrO2) is analyzed by means of a self-consistent tight-binding model. Using this model, energies and structures of zirconia containing different vacancy concentrations are calculated, equivalent in concentration to the charge compensating vacancies associated with dissolved yttria (Y2O3) in the tetragonal and cubic phase fields (3.2 and 14.4% mol, respectively). The model is shown to predict the large relaxations around an oxygen vacancy, and the clustering of vacancies along the 111 directions, in good agreement with experiments and first principles calculations. The vacancies alone are shown to explain the stabilization of cubic zirconia, and the mechanism is analyzed. (C) 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

Computational challenges in atomic, molecular and optical physics

Six challenges are discussed. These are the laser-driven helium atom; the laser-driven hydrogen molecule and hydrogen molecular ion: electron scattering (with ionization) from one-electron atoms; the vibrational and rotational structure of molecules such as H-3(+) and water at their dissociation limits; laser- heated clusters; and quantum degeneracy and Bose-Einstein condensation. The first four concern fundamental few-body systems where use of high-performance computing (HPC) is currently making possible accurate modelling from first principles. This leads to reliable predictions and support for laboratory experiment as well as true understanding of the dynamics. Important aspects of these challenges addressable only via a terascale facility are set out. Such a facility makes the last two challenges in the above list meaningfully accessible for the first time, and the scientific interest together with the prospective role for HPC in these is emphasized.

Hydrogen bonding and collective proton modes in clusters and periodic layers of squaric acid: A density functional study

Hydrogen bonding in clusters and extended layers of squaric acid molecules has been investigated by density functional computations. Equilibrium geometries, harmonic vibrational frequencies, and energy barriers for proton transfer along hydrogen bonds have been determined using the Car-Parrinello method. The results provide crucial parameters for a first principles modeling of the potential energy surface, and highlight the role of collective modes in the low-energy proton dynamics. The importance of quantum effects in condensed squaric acid systems has been investigated, and shown to be negligible for the lowest-energy collective proton modes. This information provides a quantitative basis for improved atomistic models of the order-disorder and displacive transitions undergone by squaric acid crystals as a function of temperature and pressure. (C) 2001 American Institute of Physics.

Ferroelectricity and isotope effects in hydrogen-bonded KDP crystals

Based on an accurate first principles description of the energetics in H-bonded potassium-dihydrogen-phosphate crystals, we conduct a first study of nuclear quantum effects and of the changes brought about by deuteration. Tunneling is allowed only for clusters involving correlated protons and heavy ion displacements, the main effect of deuteration being a depletion of the proton probability density at the O-H-O bridge center, which in turn weakens its proton-mediated covalent bonding. The ensuing lattice expansion couples self-consistently with the proton off-centering, thus explaining both the giant isotope effect and its close connection with geometrical effects.

Interplay between proton ordering and ferroelectric polarization in H-bonded KDP-type crystals

The origin of ferroelectricity in KH2PO4 (KDP) is studied by first-principles electronic structure calculations. In the low-temperature phase, the collective off-centre ordering of the protons is accompanied by an electronic charge delocalization from the near and localization at the far oxygen within the O-H...O bonds. Electrostatic forces. then, push the K+ ions towards off-centre positions, and induce a macroscopic polarization. The analysis of the correlation between different geometrical and electronic quantities, in connection with experimental data. supports the idea that the role of tunnelling in isotopic effects is irrelevant. Instead, geometrical quantum effects appear to play a central role. (C) 2001 Elsevier Science B.V. All rights reserved.

Relevance of heterometallic binding energy for metal underpotential deposition

We present first-principles calculations for a number of metals adsorbed on several different metallic substrates. Some of these systems are very relevant in electrochemistry, especially in the field of underpotential deposition phenomena. The present studies reveal the existence of a relationship between the excess binding energy and the surface energy difference between substrate and adsorbate. Comparisons with experimental underpotential shifts show that excess binding energies are systematically underestimated. By analyzing experimental information on different systems, we conclude that this discrepancy between our vacuum calculations and experiments carried out in an electrolytic solution is likely to be due to anion adsorption and/or solvent effects.

Broken symmetry and pseudogaps in ropes of carbon nanotubes

Since the discovery of carbon nanotubes, it has been speculated that these materials should behave like nanoscale wires with unusual electronic properties and exceptional strength. Recently, 'ropes' of close-packed single-wall nanotubes have been synthesized in high yield. The tubes in these ropes are mainly of the (10,10) type3, which is predicted to be metallic. Experiments on individual nanotubes and ropes indicate that these systems indeed have transport properties that qualify them to be viewed as nanoscale quantum wires at low temperature. It has been expected that the close-packing of individual nanotubes into ropes does not change their electronic properties significantly. Here, however, we present first-principles calculations which show that a broken symmetry of the (10,10) tube caused by interactions between tubes in a rope induces a pseudogap of about 0.1 eV at the Fermi level. This pseudogap strongly modifies many of the fundamental electronic properties: we predict a semimetal-like temperature dependence of the electrical conductivity and a finite gap in the infrared absorption spectrum. The existence of both electron and hole charge carriers will lead to qualitatively different thermopower and Hall-effect behaviours from those expected for a normal metal.