Current ethanol oxidation catalysts in direct ethanol fuel cells (typically platinum-based) suffer from low conversion and are susceptible to CO poisoning. Therefore we determined to find viable alternative catalysts for ethanol oxidation based on iridium using density functional theory to model bimetallic alloy (111) surfaces. Iridium was alloyed with another transition metals M in an overlayer (one layer of metal M on top of bulk iridium) or subsurface configuration (M is inserted under the first layer of iridium). Complete oxidation of ethanol is limited by the breaking of strong C-C bonds, so any catalyst must lower the barriers for C-C bond breaking. We modeled the reaction CH+CO →CHCO.Segregation energies were calculated and the subsurface configuration was the most stable configuration in the vast majority of alloy cases. CO adsorption was also studied and a lower CO adsorption energy was found in many alloy cases compared to pure Pt (, providing encouraging results about the possibility of reducing CO poisoning. Activation energies were lowered for the vast majority of the alloys used in an underlayer structure, reinforcing our interest in the underlayer structures or “subsurfaceâ€� alloys. Finally, we found, based on the CO adsorption energies, activation energies of the C-C breakage reaction, and metal cost, three important catalyst descriptors, a number of promising catalysts for the ethanol oxidation reaction. The most interesting alloys all adopted the underlayer structure Ir/M/Ir. With M = Ta, Hf, Nb, V, Zr, they demonstrated enhanced reactivity and high CO tolerance, having the advantage of reducing the cost of the catalyst, potentially substituting expensive platinum group metals by more affordable components.

This work 1) draws insights on molten salt structure and properties directly from ab initio molecular dynamics (AIMD) simulations, 2) demonstrates the advantageous symbiosis of computational and experimental collaborations on molten salt research, and 3) simultaneously generates ab initio data sets for fitting an interatomic potential model for classical molecular dynamics (MD) simulations. This work discusses the motivations for AIMD simulations of molten salts, thermophysical properties and structural characteristics of interest, advanced methodologies for AIMD simulations, and several completed AIMD studies on molten salts. Of import are the methodological contributions of this work to AIMD simulations, primarily the radical increase in generalized gradient planewave energy cutoff used to more accurately model the electron distribution across a highly-polarizable molten salt. Cutoffs of up to 2500 Rydbergs are used in this work, but 2000 Rydbergs is found to be sufficient for most AIMD NpT modelling of molten fluorides. The equilibrium liquid density of eutectic FLiNaK as a function of temperature is found to agree with the experimental density reported by Chrenkova et al. to within 0.2%, and the equilibrium liquid density of eutectic FMgNaK is found to agree with experimental measurements reported herein to within 4%. Self-diffusion coefficients in FMgNaK are also considered, with applicability to other halide salts. Molybdenum, Cesium, iodide, nickel, hydrogen, oxide, and uranium complexation are examined. It is found that solvation strength can be qualitatively determine via AIMD simulations, and that poorly solvated solutes will minimize the surface area of interaction with the salt solution. Cesium in particular is shown to be volatile or retainable in FLiBe at 500, 650, and 800 ËšC based on complexation and validated experimentally. It is shown that the chemical potential of an anion varies between melts as influenced by the different cations present in each melt. Hence, attempts to use a common electrochemical reference reaction for different salt mixtures are at best an approximation.

Density functional theory (DFT) calculations of polyethylene (PE) HVDC cable insulation have been performed for systems containing four different chemical impurities: acetophenone, cumene, $\alpha$-methyl styrene and $\alpha$-cumyl alcohol. Systems were generated by molecular dynamics (MD) equilibration at four different temperatures relevant for cable insulation applications: 277 K, 293 K, 343 K and 363 K. With the goal of gaining better measure of variations in hole and electron traps energies, four initial configurations were also stochastically generated at each temperature, which yielded four different final configurations after equilibration. The counterpoise correction scheme was implemented for DFT calculations, by distributing ghost atoms thought any empty pockets of space in between the PE chains. The PBE functional was selected for DFT simulations. The resulting band gaps were in agreement with those of earlier GGA-based studies, and thus lower by 3 eV than empirical band gaps. For all impurities, the first HOMO state and the first two LUMO states were generally located on the impurity molecule, forming one hole trap and two electron traps, but certain configurations generated increased electron trap numbers, or eliminated hole traps. No dependence could be derived between temperature and trap depth for either electron or hole traps. Mean electron trap energies were largely in agreement with results from earlier studies, they were deepest for acetophenone, and they varied by as much as 0.6 eV between different configurations. Hole traps are universally shallow and vary by up to 0.7 eV between configurations, and are similar in depth for all impurities. Results suggest that electron trap depths correlate with the presence of molecular features such as oxygen atoms and conjugated double bonds. The dependence of trap depth on the spatial configuration of the impurity molecule suggests that results could be improved by more precise quantum mechanical treatment of the dynamics of the impurity. / Täthetsfunktionalteori (DFT) har använts för beräkningar av isolering till HVDC kablar som består av polyeten innehållande fyra olika kemiska orenheter: acetofenon, kumen, alfa-metylstyren och alfa-kumylalkohol. System att studera genererades genom molekylärdynamisk ekvilibrering vid fyra olika temperaturer relevanta för tillämning till kabelisolering: 277 K, 293 K, 343 K och 363 K. För att få ett mått på de variationer som existerar i energierna på hål- och elektronfällor genererades stokastiskt fyra initialkonfigurationer vid varje temperatur, vilket fyra olika konfigurationer efter relaxering. Motviktskorrigering implementerades för DFT-beräkningar, genom att fördela "spökatomer" i de tomrum som bildas mellan PE-kedjorna i den amorfa fasen. PBE-funktionalen användes för DFT-simuleringar. De resulterande bandgapen stämde överens med tidigare GGA-baserade studier, och var därmed runt 3 eV smalare än empiriskt uppmätta bandgap. För alla orenheter var det första HOMO-tillståndet och de två första LUMO-tillstånden i allmänhet placerade på orenheten, vilket resulterade i en hålfälla och två elektronfällor, men vissa konfigurationer gav upphov till fler elektronfällor, eller eliminerade hålfällorna. Inget samband kunde härledas mellan temperaturen och djupet på fällorna för vare sig elektron- eller hålfällor. Medelvärdet på elektronfällornas energier överensstämde till stor del med resultat från tidigare studier, energierna var högst för acetofenon, och de varierade med så mycket som 0,6 eV mellan olika konfigurationer. Hålfällorna var genomgående grunda, varierade med upp till 0,7 eV mellan olika konfigurationer, och hade likartat djup för alla orenheter. Resultaten indikerar att variationerna elektronfällornas medeldjup uppstår på grund av orenheternas olika molekylära uppbyggnad: förekomst av syreatomer och konjugerade dubbelbindningar i orenheterna leder till djupare elektronfällor. Det faktum att djupet på elektron- och hålfällor varierar mellan olika rumsliga konfigurationer av av orenheten och polyetenstrukturen ger en antydan om att resultaten kan komma att förbättras om dynamiken hos orenheten simuleras med mer exakta kvantmekanisk metoder.