京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第20368号 / 工博第4305号 / 新制||工||1667(附属図書館) / 京都大学大学院工学研究科材料工学専攻 / (主査)教授 田中 功, 教授 乾 晴行, 教授 邑瀬 邦明 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

The main objective of this work is to gain basis for rational design of catalysts used in fuel cells for conversion of chemical energy stored in hydrogen molecules into electric energy, as well as photo-catalysts used for hydrogen production from water under solar irradiation. This objective is achieved by applying the first principles computational approach to reveal relationship among compositions of materials under consideration, their electronic structure and catalytic activity. A major part of the work is focused on electro-catalysts for hydrogen fuel cells. Platinum (Pt) is widely used in the electrodes of fuel cells due to its good catalytic properties. However, Pt is an expensive and scarce element, its catalytic activity is not optimal and also it suffers from CO poisoning at anode. Therefore the search for new catalytic materials is needed for large scale implementation of fuel cells. The main direction of search of more efficient electro-catalysts is based in the design in which an active element monoatomic layer (AE) is deposited on a metal substrate (MS) made of a cost-effective material. Two goals are achieved by doing this: on the one hand, the cost of the catalytic system is reduced by reducing the amount of the AE in the system and on the other hand the catalytic properties of the AE can be tuned through its interactions with the MS. In the first part of this work the Pd-based alloys and layered structures have been studied as promising electro-catalysts for the ORR on the fuel cell cathodes, more precisely Pd-Co alloys and Pd/M/Pd (M=Co,Fe). There exists a robust model linking the activity of a surface toward ORR to computable thermodynamic properties of the system and further to the binding energies iv of the ORR intermediates on the catalyst surface. A more challenging task is to find how to tune these binding energies through modification of the surface electronic structure that can be achieved by varying the surface composition and/or morphology. To resolve this challenge, the electronic structure, binding energies of intermediates and the ORR free energies have been calculated within the density functional theory (DFT) approximation. The results presented in this work show that in contrast to the widely accepted notion, the strain exerted by a substrate on AE hardly affects the surface activity toward ORR, while the hybridization of the electronic states of the AE-and MS-electronic states is the key factor controlling the catalytic properties of these systems. Next it is shown that the catalytic activity of the promising anode electrocatalysts, such as Pt/M, M=Au, Ru and Pd, is also determined by the AE-MS hybridization with a minor effect of the strain. Furthermore, we have shown that, if AE is weakly bound to the substrate (as it is for Pt/Au), surface reconstruction occurs. This leads to the breaking of the relation between the electronic structure of the clean surface and the reactivity of the sytem. Other kind of promising ORR catalysts is designed in the form of Ru nanoparticles modified by chalcogens. In this work, I present the results obtained for small Ru clusters and flat Ru facets modified with chalcogens (S, Se and Te). The O and OH binding energies are chosen as descriptors of the ORR. The results on the two systems are compared, concluding that large clusters with relative large flat facets have higher catalytic activity due to the absence of low coordinated and thus high reactive Ru atoms. Regarding the problem of the hydrogen production via photo-catalytic splitting of water, one of the challenges is tuning the band gap of the photo-anodes to optimal levels. Graphitic carbon nitride (g-C3N4) is a promising material to be used as a photo-anode, however, a v reduction of the band gap width by rational doping of the material would improve the efficiency significantly. This issue is addressed in the last chapter of this work. Two problems are considered: a) the stability of the doped system and b) the band gap width. To address the first problem the ab-initio thermodynamics approach has been used, finding that the substitution of C and N with the doping agent (B, C, N, O, Si and P) is thermodynamically preferred over the interstitial addition of dopant to the g-C3N4 structure. However, due to high kinetic energy barriers for the detachment of C and N atoms, involved in the substitution doping, the interstitial addition found to be kinetically more favorable. Since the density functional theory fails to reproduce the band gap of semiconductors correctly, the GW approximation was used to study the band gap of the system. The results indicate that the g-C3N4 system maintain its semiconductor character if doped with B, O and P under certain conditions, while reducing the band gap.

No description available.

This study seeks to evaluate the temperature effects on the electronic properties of thermoelectric materials, using first principles Density Functional Theory (DFT) calculations by incorporating the temperature effects on structural properties of the material. Using the electronic properties attained, the charge carrier scattering relaxation times were determined. The effect of interface between PbTe and SrTe on the charge carrier mobility was studied by finding out the relative alignment of energy bands at the semiconductor heterojunction. The crystal shape of the SrTe precipitates in the PbTe host matrix was evaluated from the interface energies using the Wulffman construction. We also attempted to develop a relation between the interface energies and electronic band alignment for different interface orientations. In this research, we incorporated the temperature effects on the structural properties of PbTe to get the temperature dependence of electronic properties like energy bandgap and effective masses of charge carriers. We used the values of bandgap and effective masses to determine the charge carrier scattering relaxation time at different temperatures which is used in evaluating the transport properties of thermoelectric materials like the Seebeck coefficient and electrical conductivity. / Master of Science

This thesis is a theoretical study of configurational disorder in icosahedral boron-rich solids, in particular boron carbide, including also the development of a methodological framework for treating configurational disorder in such materials, namely superatom-special quasirandom structure (SA-SQS). In terms of its practical implementations, the SA-SQS method is demonstrated to be capable of efficiently modeling configurational disorder in icosahedral boron-rich solids, whiles the thermodynamic stability as well as the properties of the configurationally disordered icosahedral boron-rich solids, modeled from the SA-SQS method, can be directly investigated, using the density functional theory (DFT). In case of boron carbide, especially B4C and B13C2 compositions, the SA-SQS method is used for modeling configurational disorder, arising from a high concentration of low-energy B/C substitutional defects. The results, obtained from the DFT-based calculations, demonstrate that configurational disorder of B and C atoms in boron carbide is not only thermodynamically favored at high temperature, but it also plays an important role in altering the properties of boron carbide − for example, restoration of higher rhombohedral symmetry of B4C, a metal-to-nonmetal transition and a drastic increase in the elastic moduli of B13C2. The configurational disorder can also explain large discrepancies, regarding the proper- ties of boron carbide, between experiments and previous theoretical calculations, having been a long standing controversial issue in the field of icosahedral boron- rich solids, as the calculated properties of the disordered boron carbides are found to be in qualitatively good agreement with those, observed in experiments. In order to investigate the configurational evolution of B4C as a function of temperature, beyond the SA-SQS level, a brute-force cluster-expansion method in combination with Monte Carlo simulations is implemented. The results demonstrate that configurational disorder in B4C indeed essentially takes place within the icosahedra in a way that justifies the focus on lowenergy defect patterns of the superatom picture. The investigation of the thermodynamic stability of icosahedral carbon-rich boron carbides beyond the believed solubility limit of carbon (20 at.% C) demonstrates that, apart from B4C generally addressed in the literature, B2.5C represented by B10Cp2(CC) is predicted to be thermodynamically stable with respect to B4C as well as pure boron and carbon under high pressure, ranging between 40 and 67 GPa, and also at elevated temperature. B2.5C is expected to be metastable at ambient pressure, as indicated by its dynamical and mechanical stabilities at 0 GPa. A possible synthesis route of B2.5C and a fingerprint for its characterization from the simulations of x-ray powder diffraction pattern are suggested. Besides modeling configurational disorder in boron carbide, the SA-SQS method also opens up for theoretical studies of new alloys between different icosahedral boron-rich solids − for example, (B6O)1−x(B13C2)x and B12(As1−xPx)2. As for the pseudo-binary (B6O)1−x(B13C2)x alloy, it is predicted to display a miscibility gap resulting in B6O-rich and either ordered or disordered B13C2-rich domains for intermediate global compositions at all temperatures up to melting points of the materials. However, some intermixing of B6O and B13C2 to form solid solutions is also predicted at high temperature. A noticeable mutual solubility of icosahedral B12As2 and B12P2 in each other to form B12(As1−xPx)2 disordered alloy is predicted even at room temperature, and a complete closure of a pseudo-binary miscibility gap is achieved at around 900 K. Apart from B12(As1−xPx)2, the thermodynamic stability of other compounds and alloys in the ternary B-As-P system is also investigated. For the binary B-As system, zincblende BAs is found to be thermodynamically unstable with respect to icosahedral B12As2 and gray arsenic at 0 K and increasingly so at higher temperature, indicating that BAs may merely exist as a metastable phase. This is in contrast to the binary B-P system, in which zinc-blende BP and icosahedral B12P2 are both predicted to be stable. Owing to the instability of BAs with respect to B12As2 and gray arsenic, only a tiny amount of BAs is predicted to be able to dissolve in BP to form BAs1−xPx disordered alloy at elevated temperature. For example, less than 5% BAs can dissolve in BP at 1000 K. As for the binary As-P system, As1−xPx disordered alloys are predicted at elevated temperature − for example, a disordered solid solution of up to ∼75% As in black phosphorus as well as a small solubility of ∼1% P in gray arsenic at 750 K, together with the presence of miscibility gaps. The thermodynamic stability of three different compositions of α-rhombohedral boron-like boron subnitride, having been proposed so far in the literature, is investigated. Those are, B6N, B13N2, and B38N6, represented respectively by B12(N-N), B12(NBN), and [B12(N-N)]0.33[B12(NBN)]0.67. It is found that, out of these sub- nitrides, only B38N6 is thermodynamically stable from 0 GPa up to ∼7.5 GPa, depending on the temperature, and is thus concluded as a stable composition of α-rhombohedral boron-like boron subnitride.

Os métodos de primeiros princípios dentro da Teoria do Funcional Densidade têm se desenvolvido bastante, devido aos avanços computacionais ocorridos nas últimas décadas. Por outro lado, a Termodinâmica e a Mecânica Estatística têm representado um papel importante na compreenssão da Física de Materiais, em especial no estudo dos materiais intermetálicos ordenados. A aliança entre estas duas abordagens têm se tornado cada vez mais factível, uma vez que os cálculos de estrutura eletrônica de primeiros princípios são hoje capazes de proporcionar resultados extremamente precisos para energias de formação de compostos estequiométricos. É dentro deste contexto que neste trabalho investigamos o diagrama de fases composição-temperatura do sistema Fe-Al-Mo, na estrutura cúbica de corpo centrado utilizando o método Full-Potential Linear Augmented Plane Wave (FP-LAPW) aliado ao Método Variacional de Clusters (CVM) na aproximação do tetraedro irregular. Através do método FP-LAPW determinamos a energia total de configurações cristalinas do sistema Fe-Mo-Al cúbico de corpo centrado. Estes valores são utilizados como parâmetros de entrada do CVM para determinação do potencial termodinâmico do sistema em suas diferentes fases e os correspondentes equilíbrio entre estas fases em função da composição e da temperatura, ou seja, o diagrama de fases do sistema. Embora o Fe-Al tenha grande interesse tecnológico, o seu comportamento magnético é bastante complicado e normalmente mal descrito pelos métodos de cálculos usuais. Este trabalho procura enfrentar este problema com ferramentas \"estados de arte\" e apontam para as limitações inerentes ao procedimento geral aqui utilizado. / First-principles methods within the framework of Density Functional Theory, have been developed sufficiently, due to the computational advances occurrences in the last decades. On the other hand, the Thermodynamics and the Statistical Mechanics have represented an important role in the understanding of the Physics of Materials, in special in the study to ordered intermetallic compounds. The alliance between these two approaches has become each more feasible with the time due to the fact that electronic structure calculations of first principles are today to provide almost exact results for the formations energies of stoichiometric compounds. It is inside of this context that in this work we investigate the composition-temperature phases diagram of the Fe-Al-Mo system, in the body centered cubic structure (BCC), using the Full-Potential Liner Augmented Plane Wave (FP-LAPW) method allied to the Cluster Variation Method (CVM) in the irregular tetrahedron approximation. Through the FP-LAPW method we determined the total energy of crystalline configurations of BCC Fe-Mo-Al system. These values are used as input parameters of the CVM for the determination of the thermodynamic potential of the system in its different phases and corresponding to the equilibria between these phases as a function of the composition and the temperature, that is the phase diagram of the system. Although the Fe-Al has great technological interest, its magnetic behavior is complicated and badly described by the methods used in the actual ab-initio calculations. This work faces this problem with state-of-art tools and points to the inherent limitations of the general procedure used here.

The structural and mechanical properties of s-I methane hydrate have been investigated by first principles calculations. For the first time, the fully elastic constant tensor of s-I methane hydrate is obtained entirely ab-initio. The calculated lattice parameter, bulk modulus, and elastic constants were found to be in good agreement with experimental data at ambient pressure. The Young modulus, Poisson ratio and bulk sound velocities are estimated from the calculated elastic constants and compared with wave speed measurements available.

Orientador: Antonio Carlos Dias Angelo / Banca: Hamilton Brandão Varela de Albuquerque / Banca: Nelson Henrique Morgon / O Programa de Pós-Graduação em Ciência e Tecnologia de Materiais, PosMat, tem caráter institucional e integra as atividades de pesquisa em materiais de diversos campi da Unesp / Resumo: A utilização das fases intermetálicas ordenadas como eletrocatalisadores em células a combustível já pode ser considerada como uma solução iminente para os problemas que envolvem a eficiência e as questões econoômicas. Para assegurar que as propriedades geométricas e eletrônicas destes materiais sejam realmente as almejadas para atender a todas as exigências na eletrocatálise das reações de oxidação das moléculas é necessário um estudo aprofundado de caracterização das fases intermetálicas, comumente realizado por criteriosas técnicas experimentais. Entretanto, experimentalmente, a caracterização destes materiais não fornece informações precisas que permitam correlacionar as propriedades dos materiais com o seu desempenho frente a uma dada reação eletrostática. Desta forma, uma estratégia metodológica para se obter um conhecimento mais adequado no estudo das fases intermetálicas é a utilização de métodos computacionais, baseados na Teoria do Funcional da Densidade (DFT). A metodologia empregada neste trabalho aborda uma sistemática para a otimização das propriedades geométricas através da minimização da energia total do sistema, bem como uma avaliação da estrutura eletrônica para estes materiais por meio de projeções sobre os orbitais atômicos na densidade de estados e de mapas de densidade de carga. O processo de otimização é feito por cálculos de campo auto-consistente sucessivos que variam o parâmetro de rede até encontrar uma estrutura que possua energia mínima, este processo pode ser realizado de duas formas, manual e automaticamente pelo código computacional, os resultados obtidos mostram que ambas as formas possuem a mesma precisão, levando a valores quase idênticos e que permitem reproduzir bem os cristais para os materiais estudados. A análise comparativa entre os dados cristalográficos da literatura e os resultados... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: The use of ordered intermetallic phases as electrocatalysts in fuel cells can now be regarded as an imminent solution for the problems concerning the efficiency of the device and for economic issue. To ensure that the geometric and electronic properties of these materials are actually suitable for the requirements in the electrocatalysis of melecules oxidation reactions need a meticulous characterization of the intermetallic phases, in general done by standard experimental techniques. However, the characterization of these materials performed solely experimentally does not provide accurate information to enable correlation of the properties of the materials with their performance against a given electrocatalytic reaction. Thus, a methodological strategy for obtaining a better knowledge in the study of ordered intermetallic phases is the use of computational methods, based on the Density Functional Theory. The methodology used in in thius research presents a sytematic optimization of the geometric properties by minimizing the total energy of the system and an evaluation of the electronic structure for these materials by means of the density of states projected onto atomic orbital and charge density maps. The optimization process is done by successively self-consistent field calculations that very the lattice parameter to find a structure that has a minimum energy, this process can be accomplished in two ways, manually and automatically by the computer code, the results show that both forms have the same precision, leading to almost identical identical values and allow to reproduce well the crystals of the studied materials. A comparative analysis of the crystallographic data from the literature and the results presented here show very small errors (in the order of 2-3% for most of the materials), which can be attributed exclusively to the various mathematical approaches applied... (Complete abstract click electronic access below) / Mestre

Os métodos de primeiros princípios dentro da Teoria do Funcional Densidade têm se desenvolvido bastante, devido aos avanços computacionais ocorridos nas últimas décadas. Por outro lado, a Termodinâmica e a Mecânica Estatística têm representado um papel importante na compreenssão da Física de Materiais, em especial no estudo dos materiais intermetálicos ordenados. A aliança entre estas duas abordagens têm se tornado cada vez mais factível, uma vez que os cálculos de estrutura eletrônica de primeiros princípios são hoje capazes de proporcionar resultados extremamente precisos para energias de formação de compostos estequiométricos. É dentro deste contexto que neste trabalho investigamos o diagrama de fases composição-temperatura do sistema Fe-Al-Mo, na estrutura cúbica de corpo centrado utilizando o método Full-Potential Linear Augmented Plane Wave (FP-LAPW) aliado ao Método Variacional de Clusters (CVM) na aproximação do tetraedro irregular. Através do método FP-LAPW determinamos a energia total de configurações cristalinas do sistema Fe-Mo-Al cúbico de corpo centrado. Estes valores são utilizados como parâmetros de entrada do CVM para determinação do potencial termodinâmico do sistema em suas diferentes fases e os correspondentes equilíbrio entre estas fases em função da composição e da temperatura, ou seja, o diagrama de fases do sistema. Embora o Fe-Al tenha grande interesse tecnológico, o seu comportamento magnético é bastante complicado e normalmente mal descrito pelos métodos de cálculos usuais. Este trabalho procura enfrentar este problema com ferramentas \"estados de arte\" e apontam para as limitações inerentes ao procedimento geral aqui utilizado. / First-principles methods within the framework of Density Functional Theory, have been developed sufficiently, due to the computational advances occurrences in the last decades. On the other hand, the Thermodynamics and the Statistical Mechanics have represented an important role in the understanding of the Physics of Materials, in special in the study to ordered intermetallic compounds. The alliance between these two approaches has become each more feasible with the time due to the fact that electronic structure calculations of first principles are today to provide almost exact results for the formations energies of stoichiometric compounds. It is inside of this context that in this work we investigate the composition-temperature phases diagram of the Fe-Al-Mo system, in the body centered cubic structure (BCC), using the Full-Potential Liner Augmented Plane Wave (FP-LAPW) method allied to the Cluster Variation Method (CVM) in the irregular tetrahedron approximation. Through the FP-LAPW method we determined the total energy of crystalline configurations of BCC Fe-Mo-Al system. These values are used as input parameters of the CVM for the determination of the thermodynamic potential of the system in its different phases and corresponding to the equilibria between these phases as a function of the composition and the temperature, that is the phase diagram of the system. Although the Fe-Al has great technological interest, its magnetic behavior is complicated and badly described by the methods used in the actual ab-initio calculations. This work faces this problem with state-of-art tools and points to the inherent limitations of the general procedure used here.

<p>The solid oxide fuel cell (SOFC) technology has been under research since thelate 1950s, and most of the research has been on designs utilizing yttria stabilized zirconia (YSZ) as the electrolyte of choice. However, the SOFC technology has the major drawback of requiring high operation temperatures (up to 1000 degrees Celcius), so research of alternative materials have come into interest that would possibly require a lower working temperature without any significant loss of conductivity.One such material of interest for the electrolyte is compounds of ceriumdioxide (ceria). Ceria is well known for its ability to release oxygen by formingoxygen vacancies under oxygen-poor conditions, which increases its oxygen ionconductivity, and works at a lower temperature than the YSZ compounds whenproperly doped. Conversely, ceria is also able to absorb oxygen under oxygen-rich conditions, and those two abilities make it a very good material to use in catalytic converters for reduction of carbon monoxide and nitrogen oxide emission. The ability for the oxygen ions to easily relocate inbetween the different lattice sites is likely the key property of oxygen ion transportation in ceria. Also, in oxygen-rich conditions, the absorbed oxygen atom is assumed to join the structure at either the roomy octrahedral sites, or the vacant tetrahedral sites. Following that, the oxygen atom may relocate to other vacant locations, given it can overcome a possible potential barrier.</p><p>This thesis studies how those interstitial oxygen vacancies (defects) affect theenergy profile of ceria-based supercells by first principles calculations. The system is modeled within the density functional theory (DFT) with aid of (extended) local density approximation (LDA+U) using the software VASP. Furthermore, it is studied how those vacancies affect neighbouring oxygen atoms, and wether or not it is energetically benificial for the neighbouring atoms to readjust their positions closer or further away from the vacancy. The purpose of this thesis is to analyze wether or not it is theoretically possible that interstitial oxygen vacancies may cause neighbouring oxygen atoms to naturally relocate to the octahedral site in ceria, and how this affects the overall energy profile of the material.</p>