The aim of this thesis is to investigate the effects of additives, reactive elements and impurities, on the lifetime of thermal barrier coatings. The thesis consists of a number of studies on interface adhesion, impurity diffusion, grain boundary sliding and cleavage processes and their impact on the mechanical behaviour of grain boundaries. The effects of additives and impurity on interface adhesion were elaborated by using total energy calculations, electron localization and density of states, and by looking into the atomic separations. The results of these calculations allow the assessment of atomic level contributions to changes in the adhesive trend. Formation of new bonds across the interface is determined to improve the adhesion in reactive element(RE)-doped structures. Breaking of the cross interface bonds and sulfur(S)-oxygen(O) repulsion is found responsible for the decreased adhesion after S segregation. Interstitial and vacancy mediated S diffusion and the effects of Hf and Pt on the diffusion rate of S in bulk NiAl are studied. Hf is shown to reduce the diffusion rate, and the preferred diffusion mechanism of S and the influence of Pt are revealed to be temperature dependent. Finally, the effects of reactive elements on alumina grain boundary strength are studied. Reactive elements are shown to improve both the sliding and cleavage resistance, and the analysis of atomic separations suggest an increased ductility after the addition of quadrivalent Hf and Zr to the alumina grain boundaries.

The maximum theoretical energy efficiency of fuel cells is much larger than those of the steam-power-turbine cycles that are currently used for generating electrical power. Similarly, direct hydrocarbon fuel cells, DHFCs, can theoretically be much more efficient than hydrogen fuel cells. Unfortunately the current densities (overall reaction rates) of DHFCs are substantially smaller than those of hydrogen fuel cells. The problem is that the exchange current density (catalytic reaction rate) is orders of magnitude smaller for DHFCs. Other work at the University of Ottawa has been directed toward the development of polymer electrolytes for DHFCs that operate above the boiling point of water, making corrosion rates much slower so that precious metal catalysts are not required. Propane (liquefied petroleum gas, LPG) was the hydrocarbon chosen for this research partly because infrastructure for its transportation and storage in rural areas already exists. In this work nickel based catalysts, an inexpensive replacement for the platinum based catalysts used in conventional fuel cells, were examined using density functional theory, DFT. The heats of propane adsorption for 3d metals, when plotted as a function of the number of 3d electrons in the metal atom, had the shape of a volcano plot, with the value for nickel being the peak value of the volcano plot. Also the C-H bond of the central carbon atom was longer for propane adsorbed on nickel than when adsorbed on any of the other metals, suggesting that the species adsorbed on nickel was less likely to desorb than those on other metals. The selectivity of the propyl radical reaction was examined. It was found that propyl radicals

Intermetallic compounds are solid phases containing two or more metallic elements, whose crystal structure differs from that of its constituents [1]. The largest group among these compounds with more than 1400 binary and ternary representatives are the so-called Laves phases. The classification of an intermetallic compound as a Laves phase is solely based on the atomic configuration and the component ratio in the crystal structure. With the ideal composition AB2, the Laves phases crystallize in three closely related structure types which are named after their representatives, MgCu2 (C15, cubic), MgZn2 (C14, hexagonal) and MgNi2 (C36, hexagonal). Laves phases are built by almost all metals of the periodic system of the elements. A significant feature of many of these is the formation of broad homogeneity ranges by mutual substitution of atoms in combination with composition or temperature dependent phase transformations between the different Laves phase polytypes. Laves phases have received considerable attention in recent years as potential structural and functional materials. They combine high melting points with considerable creep resistance, high strength and fracture toughness and good corrosion and oxidation resistance. Some Laves phases like NbFe2 [2] or TaFe2 [3] show intriguing magnetic and electronic properties which provide a deeper inside into phenomenons like quantum criticality. Especially transition-metal based Laves phases like NbCr2 [4] and ZrCr2 [5] are promising candidates for the development of new high-temperature structural materials. The major drawback of the Laves phases, however, is their low-temperature brittleness. Many experimental and theoretical investigations have shown, that the low-temperature ductility can be improved by controlling the crystal structure with the help of phase transformations, by mechanical twinning or the addition of third elements. The addition of ternary alloying elements can alter the physical and electronic properties of the Laves phases and plays an important role in the composition dependent stability of the different polytypes. Some ternary Laves phases show an interesting phenomenon called site occupation reversal. It describes a composition dependent behavior of the alloying elements which prefer to occupy different crystallographic sites at different concentrations. The understanding of the point defect structure/mechanism and the site occupation of the alloying elements is thus of critical importance for the proper description of phase stability. The basis for the broad application of any metallic material is the knowledge of the corresponding phase diagram. The experimental determination of phase diagrams however, is tedious, time consuming and expensive work and the huge abundance of Laves phase makes this an impractical task. Thus, the time it takes to discover new advanced materials and to move them from the laboratory to the commercial market place is fairly long today. A cheap and fast enhancement for the development of new materials is the calculation of phase diagrams and physical properties using techniques like CALPHAD (CALculation of PHAse Diagrams) and DFT (Density Functional Theory). Very recently the Office of Science and Technology Policy of the United States White House announced to provide a budget of $100 million to launch the Materials Genome Initiative [6, 7]. The aim of this initiative is to provide the infrastructure and training needed to discover, develop, manufacture, and deploy advanced materials in a more expeditious and economical way [8]. One of the project’s three supporting legs is the calculation and prediction of crystal structures and physical properties using advanced Computational tools. "An early benchmark will be the ability to incorporate improved predictive modeling algorithms of materials behavior into existing product design tools. For example, the crystal structure and physical properties of the materials [. . . ]." [8]. Their computational tools of choice are the same as used in this work to predict crystal structures and site occupation factors. Contents of this work is the investigation of the substitutional disorder in binary and ternary Laves phases. This includes the experimental determination of the composition dependent stability of the Laves phase polytypes and the distribution of the substitution atoms in the crystal lattice of the respective phases, i.e., the site occupation factors (s.o.f.). For this purpose, detailed experimental studies on the two systems Cr–Co–Nb and Fe–Ta–V were performed and the Laves phase polytypes, their homogeneity ranges, the lattice parameters and the site occupation factors were determined. The experimental results are compared with the results obtained from quantum mechanical calculations. DFT is used to determine the composition dependent enthalpies of formation which serve as a measure for the stability of the different Laves phase polytypes. Additionally, the applicability of various approximations and their influence on the results has been checked. This study is thus also supposed to develop and improve the tools necessary for the calculation of phase stability and homogeneity ranges in ternary phases. Chapter two in the first part of this work describes the crystal structures of the Laves phases in detail with focus on the polytype stability, the site occupation and the c/a-ratio of hexagonal C14 Laves phases. Subsequently, the phase diagrams of the investigated systems and the occurring Laves phases are discussed. Chapter three briefly describes the experimental and theoretical methods used in this work. The last section of part one gives a detailed explanation of how the phase stability, the lattice parameters and the site occupation factors are calculated. The second part "Results and discussion" contains the discussion of the experimental and theoretical results for the intensively investigated systems Co–Cr–Nb (chapter five) and Fe–Ta–V (chapter six). Several other ternary C14 Laves phases and their site occupation behavior are studied in chapter seven. The thesis is concluded with a summary in chapter eight. Several additional information is contained in the appendix.

The aim of this thesis is to investigate the effects of additives, reactive elements and impurities, on the lifetime of thermal barrier coatings. The thesis consists of a number of studies on interface adhesion, impurity diffusion, grain boundary sliding and cleavage processes and their impact on the mechanical behaviour of grain boundaries. The effects of additives and impurity on interface adhesion were elaborated by using total energy calculations, electron localization and density of states, and by looking into the atomic separations. The results of these calculations allow the assessment of atomic level contributions to changes in the adhesive trend. Formation of new bonds across the interface is determined to improve the adhesion in reactive element(RE)-doped structures. Breaking of the cross interface bonds and sulfur(S)-oxygen(O) repulsion is found responsible for the decreased adhesion after S segregation. Interstitial and vacancy mediated S diffusion and the effects of Hf and Pt on the diffusion rate of S in bulk NiAl are studied. Hf is shown to reduce the diffusion rate, and the preferred diffusion mechanism of S and the influence of Pt are revealed to be temperature dependent. Finally, the effects of reactive elements on alumina grain boundary strength are studied. Reactive elements are shown to improve both the sliding and cleavage resistance, and the analysis of atomic separations suggest an increased ductility after the addition of quadrivalent Hf and Zr to the alumina grain boundaries.

This thesis work was performed at AB Sandvik Coromant and aimed to enhance the knowledge about the relationships between the performance of TiN and TiAlN-coated cutting tools in metal turning and their mechanical and chemical properties. Measurements of coating material properties and turning wear tests in annealed tool steel Sverker 21, stainless steel 316L, grey cast iron V314 and nodular cast iron SS0727 were performed. The cutting temperatures were estimated from FEM-simulations. To find the dominant wear mechanism and identify the properties that are most important for the resistance against that particular wear, a correlation analysis was performed together with a wear study using LOM, SEM and EDS. The results show that relations between cutting performance and mechanical properties and/or composition of the coatings can be established. The FEM-simulations suggested that the peak tool temperature was highest, ~750°C, for turning in 316L and lowest for turning in Sverker 21, ~300°C. Turning in cast iron resulted in temperatures around 500-550°C. A mechanism for the growth of the crater on inserts tested in stainless steel 316L is proposed. Wear due to thermo-mechanical load and adhesion are believed to be the dominating wear mechanisms. The performance of the tool showed a high correlation to the composition of the coatings, with a decreased tool life for higher Al-contents. The reason for this might lie in an increased brittleness of these coatings, accelerating formation of lateral cracks above the crater. Calculated ratios of bulk and shear modulus suggests an increased brittleness for higher Al-contents. A higher tendency to stick to the work piece material might also contribute to a decrease in tool life. An Increased Al-content could also drive the formation of c-AlN to h-AlN, causing even higher wear rates. The coatings with higher substrate bias showed an enhanced performance, even though the crack pattern was worsened for these variants. The reason for the enhanced performance seen for these variants might instead originate in an enhanced adhesion to the substrate. In the flank wear resistance test in Sverker 21 the Al-content proved to be important, with an improved performance for higher Al-contents. In contrast to the test in 316L, a change in bias or hardness had no effect on the performance in this test. Scratch patterns on the flank supports that an abrasive wear mechanism is present, but no correlation between hardness and tool life could be obtained. Either some other material property than hardness is of importance for the abrasive resistance in this test, or another wear mechanism, occurring simultaneously with abrasion, is the wear rate deciding. The second part of this thesis work was to evaluate the ability of a quantum mechanical computational method, density functional theory, to predict material properties. The method predicts the lattice parameters and bulk moduli in excellent agreement with experimental values. The method also well predicts other elastic properties, with results consistent with reference values. There seems to be a constant shift of about 50-100 GPa between the calculated elastic modulus and the experimentally measured values, probably originating in contributions from grain boundaries, texture, stresses and defects present in the real coatings, and possibly also in errors in the experimental method due to an influence from the substrate. The calculated hardness values did not follow the trend of an increased hardness for TiAlN compared to TiN, which is seen in experiments.

In this thesis, a number of clean-energy materials for hydrogen generation, hydrogen storage, and Li-ion battery energy storage applications have been investigated through state-of-the-art density functional theory. As an alternative fuel, hydrogen has been regarded as one of the promising clean energies with the advantage of abundance (generated through water splitting) and pollution-free emission if used in fuel cell systems. However, some key problems such as finding efficient ways to produce and store hydrogen have been hindering the realization of the hydrogen economy. Here from the scientific perspective, various materials including the nanostructures and the bulk hydrides have been examined in terms of their crystal and electronic structures, energetics, and different properties for hydrogen generation or hydrogen storage applications. In the study of chemisorbed graphene-based nanostructures, the N, O-N and N-N decorated ones are designed to work as promising electron mediators in Z-scheme photocatalytic hydrogen production. Graphene nanofibres (especially the helical type) are found to be good catalysts for hydrogen desorption from NaAlH4. The milestone nanomaterial, C60, is found to be able to significantly improve the hydrogen release from the (LiH+NH3) mixture. In addition, the energetics analysis of hydrazine borane and its derivative solid have revealed the underlying reasons for their excellent hydrogen storage properties.  As the other technical trend of replacing fossil fuels in electrical vehicles, the Li-ion battery technology for energy storage depends greatly on the development of electrode materials. In this thesis, the pure NiTiH and its various metal-doped hydrides have been studied as Li-ion battery anode materials. The Li-doped NiTiH is found to be the best candidate and the Fe, Mn, or Cr-doped material follows. / <p>QC 20130925</p>

<p>The computational success of density functional theory relies on the construction of suitable approximations to the exchange-correlation energy functional. Use of functional approximations depending explicitly upon the density alone appear unable to address all aspects of many-body interactions, such as the fundamental constraint that the ground state energy is a piecewise linear function of the total number of electrons, and the ability to model nonlocal effects. Functionals depending explicitly upon occupied and unoccupied Kohn–Sham orbitals are considered necessary to address these and other issues. This dissertation considers certain issues relevant to the successful implementation of explicitly orbital-dependent functionals through the optimized effective potential (OEP) approach, as well as extending the potential functional formalism that provides the formal basis for the OEP approach to systems in the presence of noncollinear magnetic fields.</p><p>The self-consistent implementation of orbital-dependent energy functionals is correctly done through the optimized effective potential approach—minimization of the ground state energy with respect to the Kohn–Sham potential that generates the set of orbitals employed in the energy evaluation. The focus on the potential can be problematic in finite basis set approaches as determining the exchange-correlation potential in this manner is an inverse problem which, depending upon the combination of orbital and potential basis sets employed, is often ill-posed. The ill-posed nature manifests itself as nonphysical exchange-correlation potentials and total energies. We address the problem of determining meaningful exchange-correlation potentials for arbitrary combinations of orbital and potential basis sets through an L-curve regularization approach based on biasing towards smooth potentials in the energy minimization. This approach generates physically reasonable potentials for any combination of basis sets as shown by comparisons with grid-based OEP calculations on atoms, and through direct comparison with DFT calculations employing functionals not depending on orbitals for which OEP can also be performed. This work ensures that the OEP methodology can be considered a viable many-body computational methodology.</p><p>A separate issue of our OEP implementation is that it can suffer from a lack of size-extensivity—the total energy of a system of infinitely separated monomers may not scale linearly with the total number of monomers depending upon how we construct the Kohn–Sham potential. Typically, a fixed reference potential is employed to aid in the convergence of a finite basis set expansion of the Kohn–Sham potential. This reference potential can be utilized to ensure other desirable properties of the resulting potential. In particular, it can enforce the correct asymptotic behavior. The Fermi–Amaldi potential is often used for this purpose but suffers from size-nonextensivity owing to the explicit dependence of the potential on the total number of electrons. This error is examined and shown to be rather small and rapidly approaches a limiting linear behavior. A size-extensive reference potential with the correct asymptotic behavior is suggested and examined.</p><p>We also consider a formal aspect of the potential-based approach that provides the underlying justification of the OEP methodology. The potential functional formalism of Yang, Ayers, and Wu is extended to include systems in the presence of noncollinear magnetic fields. In doing so, a solution to the nonuniqueness issue associated with mapping between potentials and wave functions in such systems is provided, and a computational implementation of the OEP in noncollinear systems is suggested.</p><p>Finally, as an example of an issue for which orbital-dependent functionals seem necessary to obtain a correct description, we consider the ground state structures of C_{4<italic>N</italic> + 2} rings which are believed to exhibit a geometric transition from angle-alternation (<italic>N</italic> ≤ 2) to bond-alternation (<italic>N</italic> > 2). So far, no published DFT approach has been able to reproduce this behavior owing to the tendency of common density functional approximations to bias towards delocalized electron densities. Calculations are presented with the rCAM-B3LYP exchange-correlation functional that correctly predict the structural evolution of this system. This is rationalized in terms of the recently proposed delocalization error for which rCAM-B3LYP explicitly attempts to address.</p> / Dissertation

In this study, we used the Density functional theory (DFT) to obtain the relationship between mechanical property and electronic property of Boron nitride nanotubes (BNNTs) under the uni-axial strain. Moreover, we also investigated one CO molecule adsorbed on the BNNTs under the uni-axial strain. We also use the molecular dynamics to introduce the mechanical property and dynamic behavior of (8,8)BNNT under the uni-axial strain. There were three parts in this study: The first part: The effect of uni-axial strain on the electronic properties of (5,5) and (8,0)boron nitride nanotubes were obtained by DFT calculation. We used the HOMO-LUMO Gap¡Bbond angle¡Bbond length and radial buckling to analyze the electronic properties and mechanical properties. The stress-strain profiles indicated that different BNNTs types displayed very similar mechanical properties, but there were variations in HOMO-LUMO gaps at different strains, indicating that the electronic properties of BNNTs not only depend on uni-axial strain, but on BNNT type. In addition, the variations in nanotube geometries, partial density of states (PDOS) and charges of boron and nitride atoms were also discussed for (8,0) and (5,5) BNNTs at different strains. The second part: The DFT was used to investigate electronic properties of CO molecule adsorbed on BNNT under the uni-axial strain. The stress-strain profiles indicated that the CO molecule adsorption on BNNT leaded only to a local mechanical deformation. The strength of BNNT could not be affected when the CO molecule adsorbed on that. Moreover, we obtained that the charge of CO will slightly transfer to the adsorbed atom of BNNT when strain increased. Hence, the adsorption energy increased slightly under the uni-axial strain. The third part: The molecular dynamics simulations were performed to investigate deformation behaviors of (8,8)BN nanotubes under axial tensile strains at 300k. Variations with the tensile strain in the axial stress, bond lengths, bond angles, radial buckling, and slip vectors were all examined. The axial, radial, and tangential components of the slip vector were also employed to monitor, respectively, the local elongation, necking, and twisting deformation near the failure of the nanotube. The components of the slip vector grew rapidly and abruptly after the failure strain, especially for the axial component. This implies that the local elongation dominates the failure of the loaded BN nanotube and finally results in a chain-like tensile failure mode.

In this study, we employed density functional theory (DFT) and molecular dynamics (MD) to investigate the mechanical and electronic properties of one-dimensional zinc oxide nanostructure. This study can be arranged into two parts: In part I: We investigated the mechanical and electronic properties of one-dimensional zinc oxide nanostructure under axial mechanical deformations by density functional theory. In this case, we could find both the highest occupied molecular orbital and the lowest unoccupied molecular orbital gap (HOMO-LUMO gap) and value of radial buckling will decrease linearly with the increase of axial strain. The changes of bond lengths and bond angles show the variation of nanostructure dependence to the increase of axial strain. This study also used partial density of state (PDOS), bond order (BO) and deformation density to analyse the differences of the electronic properties between the zinc oxide nanotubes under axial strain. In part II: This study, which employed molecular dynamics combines Buckingham and Core-Shell potentials, shows the different physical parameters, such as yield stress, young¡¦s modulus and slip vector to research the mechanical behavior and variation of structure of nanotube under axial strain.

The structural and electronic properties of small tungsten nanoparticles Wn (n=2-16) were investigated by density functional theory (DFT) calculation. For the W10 nanoparticle, ten lowest-energy structures were first obtained by basin-hopping method (BH) and ten by big-bang method (BB) with the tight-binding many-body potential for bulk tungsten material. These fifty structures were further optimized by the DFT calculation in order to find the better parameters of tight-binding potential adquately for W nanoparticles. With these modified parameters of tight-binding potentials, several lowest-energy W nanoparticles of different sizes can be obtained by BH and BB methods and then further refined by DFT calculation. According to the values of binding energy and second-order energy difference, it reveals that the structure W12 has a relatively higher stability than those of other sizes. The vertical ionization potential (VIP), adiabatic electron affinity (AEA) and HOMO-LUMO Gap are also discussed for W nanoparticles of different sizes.