A theoretical study on the photoluminescence of SrTiO3

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Ab-initio determination of x-ray absorption near edge structure (xanes) spectra in an ultrasoft and norm conserving pseudopotentials scheme

X-ray absorption spectroscopy (XAS) is a powerful means of investigation of structural and electronic properties in condensed -matter physics. Analysis of the near edge part of the XAS spectrum, the so – called X-ray Absorption Near Edge Structure (XANES), can typically provide the following information on the photoexcited atom: - Oxidation state and coordination environment. - Speciation of transition metal compounds. - Conduction band DOS projected on the excited atomic species (PDOS). Analysis of XANES spectra is greatly aided by simulations; in the most common scheme the multiple scattering framework is used with the muffin tin approximation for the scattering potential and the spectral simulation is based on a hypothetical, reference structure. This approach has the advantage of requiring relatively little computing power but in many cases the assumed structure is quite different from the actual system measured and the muffin tin approximation is not adequate for low symmetry structures or highly directional bonds. It is therefore very interesting and justified to develop alternative methods. In one approach, the spectral simulation is based on atomic coordinates obtained from a DFT (Density Functional Theory) optimized structure. In another approach, which is the object of this thesis, the XANES spectrum is calculated directly based on an ab – initio DFT calculation of the atomic and electronic structure. This method takes full advantage of the real many-electron final wavefunction that can be computed with DFT algorithms that include a core-hole in the absorbing atom to compute the final cross section. To calculate the many-electron final wavefunction the Projector Augmented Wave method (PAW) is used. In this scheme, the absorption cross section is written in function of several contributions as the many-electrons function of the finale state; it is calculated starting from pseudo-wavefunction and performing a reconstruction of the real-wavefunction by using a transform operator which contains some parameters, called partial waves and projector waves. The aim of my thesis is to apply and test the PAW methodology to the calculation of the XANES cross section. I have focused on iron and silicon structures and on some biological molecules target (myoglobin and cytochrome c). Finally other inorganic and biological systems could be taken into account for future applications of this methodology, which could become an important improvement with respect to the multiscattering approach.

Thiophene-containing organic semiconducting heteroacenes for electronic applications

Diese Studie verfolgt das Konzept der "Oligomer-Ansatz", die von Müllen angesprochen wurde et.al. vor etwa 10 Jahren. Der Schwerpunkt dieser Arbeit war die Synthese, Charakterisierung und Anwendung von halbleitenden konjugierten heteroacenes für organische Elektronik.rnZur weiteren Entwicklung der Familie von schwefelhaltigen Pentacene, zwei Moleküle (Benzo [1,2 - b :4,5-b '] bis [b] benzodithiophene und dithieno-[2,3-d: 2', 3'-d ']-benzo-[1,2-b :4,5-b'] dithiophene)rnfacilely wurden synthetisiert und charakterisiert durch eine Kombination verschiedener Methoden. Die beiden neue Moleküle weisen hervorragende ökologische Stabilität und angewendet OFETs Geräte als p-Kanal-Material. Die Vorversuche gaben Ladungsträgerbeweglichkeiten von 0,1 cm2 V-1 s-1 undrn1,6 cm2 V-1 s-1 bzw. aus den beiden Molekülen.rnAusgelöst durch die Frage "je länger desto besser?", Eine Reihe von neuen heteroheptacenes wurden synthetisiert und intensiv im Hinblick auf ihre feste Struktur, Selbst-assenbly auf der studierte Oberfläche, opto-elektro-Eigenschaften und Eigenschaften des Orbits Grenze. Einer derrnheteroheptacene Moleküle wurden als die aktiven Kanäle in OFET Geräten angewendet. Jedoch in Trotz der mehr verlängert Konjugationslänge die Geräte auf der Basis zeigten heptacenes viel schlimmer Ladungsträgerbeweglichkeiten als die heteropentacenes. Viele Faktoren können Festlegung der endgültigen Leistung der Produkte und der chemischen Struktur ist nur einer von ihnen.rnIn dieser Hinsicht scheint es, dass es auch sinnvoll, den Einfluss der Heteroatome Studie und Alkylsubstitution auf der soliden und elektronischen Strukturen. Daher mehr heteroheptacenes wurden synthetisiert. Abwechslungsreiches in der Anzahl und Art der heteroatomare Brücke,rndiese Oligoazene ausgestellt dramatisch anders feste Struktur und opto-elektronischernEigenschaften. Darüber hinaus wurde eine kombinierte DFT Berechnung der Molekülorbitale dieser heptacenes darauf hingewiesen, dass die Einführung von Stickstoff Brücken wird die π-Orbitale zu destabilisieren, während stabilisieren den Schwefel Brücken sowohl HOMO und LUMO Energien. Dies ist wichtig, wenn man will hoch π verlängert Oligoazene synthetisieren und dabei eine angemessene Stabilität.

Theoretical study of correlated systems using hybrid functionals

Diese Arbeit unterstreicht das Potential von Hybridfunktionalen (B3LYP) für die Untersuchung einer großen Bandbreite von Systemen. Durch die Einbeziehung der exakten Hartree-Fock Austauschenergie kann B3LYP für molekulare und kristalline Systeme eingesetzt werden. Zum Beispiel können stark korrelierte Systeme mit B3LYP erfolgreich erforscht werden. Die elektronische Struktur von PAHs wurde mit B3LYP Hybriddichtefunktionalen untersucht. Mit der ∆SCF-Methode wurden Elektronenbindungsenergien bestimmt, welche die mit UPS gewonnenen experimentellen Resultate bestätigen und ergänzen. Symmetrieeigenschaften der molekularen Orbitale wurden analysiert, um eine Zuordnung und Einschätzung der zugehörigen Signalstärke zu ermöglichen. Während σ-artige Orbitale nur schwer durch UPS-Messungen an dünnen Filmen detektiert werden können, bieten Rechnungen eine detaillierte Einsicht in die verborgenen Teile der Spektren.rnWeiterhin wurden π−π-Komplexe untersucht, welche von verschiedenen Donor- und Akzeptor-Molekülen gebildet werden. Die Moleküle basieren auf polyzyklischen, aromatischen Kohlenwasserstoffen. Für Ladungstransferkomplexe finden DFT Rechnungen ein Minimum in der Oberfläche der potentiellen Energie. Diese attraktive Wechselwirkung wird durch Coulombanziehung verursacht. Allerdings ist die Coulombanziehung nicht die stärkste Wechselwirkung in Ladungstransferkomplexen. Die Einbeziehung von van der Waals-Korrekturen verbessert den intermolekularen Abstand und die Bindungsenergie.rnEine Verkleinerung der intermolekularen Abstände führt zu einer großen Verschiebung der HOMO- und LUMO-Energie.rnAus der Klasse der kristallinen korrelierten Systeme wurden Rb4O6 und FeSe untersucht. Im Falle von Rb4O6 führen Ladungsordnung und Korrelationen zu einem isolierenden Grundzustand. Das hypothetische druckabhängige Phasendiagramm wurde untersucht. Eine Erhöhung des Drucks führt zu einer vergrößerten Bandlücke. Bei etwa 75 GPa wird die Bandbreite W größer als der Bandabstand U und das System nimmt einen homogen gemischt valenten Zustand mit teilweise besetzten π−π-Orbitalen an. Für Drücke ab 160 GPa wird W sehr viel größer als U und das System wird metallisch.rnIm Fall von FeSe finden wir eine korrelierte und isolierende Phase bei hohen Drücken, während das System bei niedrigen Drücken supraleitendes Verhalten zeigt. Die Berechnungen der Elektronenstruktur mit dem Hybridfunktional B3LYP führt zum korrekten halbleitenden Grundzustand in der NiAs- und MnP-Struktur von FeSe. Die Rolle der Korrelationen, der Stöchiometrie und der Nähe zum Magnetismus wird besprochen. Im Speziellen wird gezeigt, dass die Phase mit NiAs-Struktur starke lokale Korrelationen aufweist, was zu einem halbleitenden Zustand in einem weiten Druckbereich führt.

Large spin splitting in the conduction band of transition metal dichalcogenide monolayers

We study the conduction band spin splitting that arises in transition metal dichalcogenide (TMD) semiconductor monolayers such as MoS2, MoSe2, WS2, and WSe2 due to the combination of spin-orbit coupling and lack of inversion symmetry. Two types of calculation are done. First, density functional theory (DFT) calculations based on plane waves that yield large splittings, between 3 and 30 meV. Second, we derive a tight-binding model that permits to address the atomic origin of the splitting. The basis set of the model is provided by the maximally localized Wannier orbitals, obtained from the DFT calculation, and formed by 11 atomiclike orbitals corresponding to d and p orbitals of the transition metal (W, Mo) and chalcogenide (S, Se) atoms respectively. In the resulting Hamiltonian, we can independently change the atomic spin-orbit coupling constant of the two atomic species at the unit cell, which permits to analyze their contribution to the spin splitting at the high symmetry points. We find that—in contrast to the valence band—both atoms give comparable contributions to the conduction band splittings. Given that these materials are most often n-doped, our findings are important for developments in TMD spintronics.

Van der Waals-corrected density functional theory: benchmarking for hydrogen-nanotube and nanotube-nanotube interactions

Density functional theory (DFT) is a powerful approach to electronic structure calculations in extended systems, but suffers currently from inadequate incorporation of long-range dispersion, or Van der Waals (VdW) interactions. VdW-corrected DFT is tested for interactions involving molecular hydrogen, graphite, single-walled carbon nanotubes (SWCNTs), and SWCNT bundles. The energy correction, based on an empirical London dispersion term with a damping function at short range, allows a reasonable physisorption energy and equilibrium distance to be obtained for H-2 on a model graphite surface. The VdW-corrected DFT calculation for an (8, 8) nanotube bundle reproduces accurately the experimental lattice constant. For H-2 inside or outside an (8, 8) SWCNT, we find the binding energies are respectively higher and lower than that on a graphite surface, correctly predicting the well known curvature effect. We conclude that the VdW correction is a very effective method for implementing DFT calculations, allowing a reliable description of both short-range chemical bonding and long-range dispersive interactions. The method will find powerful applications in areas of SWCNT research where empirical potential functions either have not been developed, or do not capture the necessary range of both dispersion and bonding interactions.

Towards NIR emission: 2-(1H-tetrazol-1-yl)pyridine as versatile ligand for the synthesis of two new families of cationic and neutral Ir(III) complexes.

In the last decades, cyclometalated Ir(III) complexes have drawn a large interest for their unique properties: they are excellent triplet state emitters, thus the emission is phosphorescent in nature; typically high quantum yields and good stability make them good candidates for luminescent materials. Moreover, through an opportune choice of the ligands, it is possible to tune the emission along the whole visible spectra. Thanks to these interesting features, Ir(III) complexes have found different applications in several areas of applied science, from OLEDs to bioimaging. In particular, regarding the second application, a remarkable red-shift in the emission is required, in order to minimize the problem of the tissue penetration and the possible damages for the organisms. With the aim of synthesizing a new family of NIR emitting Ir(III) complexes, we envisaged the possibility to use for the first time 2-(1H-tetrazol-1-yl)pyridine as bidentate ligand able to provide the required red-shift of the emission of the final complexes. Exploiting the versatility of the ligand, I prepared two different families of heteroleptic Ir(III) complexes. In detail, in the first case the 2-(1H-tetrazol-1-yl)pyridine was used as bis-chelating N^N ligand, leading to cationic complexes, while in the second case it was used as cyclometalating C^N ligand, giving neutral complexes. The structures of the prepared molecules have been characterised by NMR spectroscopy and mass spectrometry. Moreover, the neutral complexes’ emissive properties have been measured: emission spectra have been recorded in solution at both room temperature and 77K, as well as in PMMA matrix. DFT calculation has then been performed and the obtained results have been compared to experimental ones.

Hydrogen spillover mechanism on a pd-doped mg surface as revealed by ab initio density functional calculation

The hydrogenation kinetics of Mg is slow, impeding its application for mobile hydrogen storage. We demonstrate by ab initio density functional theory (DFT) calculations that the reaction path can be greatly modified by adding transition metal catalysts. Contrasting with Ti doping, a Pd dopant will result in a very small activation barrier for both dissociation of molecular hydrogen and diffusion of atomic H on the Mg surface. This new computational finding supports for the first time by ab initio simulationthe proposed hydrogen spillover mechanism for rationalizing experimentally observed fast hydrogenation kinetics for Pd-capped Mg materials.

Correlation of Oxygen Storage Capacity and Structural Distortion in Transition-Metal-, Noble-Metal-, and Rare-Earth-Ion-Substituted CeO2 from First Principles Calculation

Oxygen storage/release (OSC) capacity is an important feature common to all three-way catalysts to combat harmful exhaust emissions. To understand the mechanism of improved OSC for doped CeO2, we undertook the structural investigation by X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), H-2-TPR (temperature-programmed hydrogen reduction) and density functional theoretical (DFT) calculations of transition-metal-, noble-metal-, and rare-earth (RE)-ion-substituted ceria. In this report, we present the relationship between the OSC and structural changes induced by the dopant ion in CeO2. Transition metal and noble metal ion substitution in ceria greatly enhances the reducibility of Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu, Pd, Pt, Ru), whereas rare-earth-ion-substituted Ce(1-x)A(x)O(2-delta) (A = La, Y) have very little effect in improving the OSC. Our simulated optimized structure shows deviation in cation oxygen bond length from ideal bond length of 2.34 angstrom (for CeO2). For example, our theoretical calculation for Ce28Mn4O62 structure shows that Mn-O bonds are in 4 + 2 coordination with average bond lengths of 2.0 and 3.06 angstrom respectively. Although the four short Mn-O bond lengths spans the bond distance region of Mn2O3, the other two Mn-O bonds are moved to longer distances. The dopant transition and noble metal ions also affects Ce coordination shell and results in the formation of longer Ce-O bonds as well. Thus longer cation oxygen bonds for both dopant and host ions results in enhanced synergistic reduction of the solid solution. With Pd ion substitution in Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu) further enhancement in OSC is observed in H-2-TPR. This effect is reflected in our model calculations by the presence of still longer bonds compared to the model without Pd ion doping. The synergistic effect is therefore due to enhanced reducibility of both dopant and host ion induced due to structural distortion of fluorite lattice in presence of dopant ion. For RE ions (RE = Y, La), our calculations show very little deviation of bonds lengths from ideal fluorite structure. The absence of longer Y-O/La-O and Ce-O bonds make the structure much less susceptible to reduction.

Methods for melting temperature calculation

Melting temperature calculation has important applications in the theoretical study of phase diagrams and computational materials screenings. In this thesis, we present two new methods, i.e., the improved Widom's particle insertion method and the small-cell coexistence method, which we developed in order to capture melting temperatures both accurately and quickly.

We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

We propose the small-cell coexistence method based on the statistical analysis of small-size coexistence MD simulations. It eliminates the risk of a metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated melting points. The method converges to the exact result in the limit of a large system size. An accuracy within 100 K in melting temperature is usually achieved when the simulation contains more than 100 atoms. DFT examples of Tantalum, high-pressure Sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the melting temperature is a design criterion.

We present in detail two examples of refractory materials. First, we demonstrate how key material properties that provide guidance in the design of refractory materials can be accurately determined via ab initio thermodynamic calculations in conjunction with experimental techniques based on synchrotron X-ray diffraction and thermal analysis under laser-heated aerodynamic levitation. The properties considered include melting point, heat of fusion, heat capacity, thermal expansion coefficients, thermal stability, and sublattice disordering, as illustrated in a motivating example of lanthanum zirconate (La₂Zr₂O₇). The close agreement with experiment in the known but structurally complex compound La₂Zr₂O₇ provides good indication that the computation methods described can be used within a computational screening framework to identify novel refractory materials. Second, we report an extensive investigation into the melting temperatures of the Hf-C and Hf-Ta-C systems using ab initio calculations. With melting points above 4000 K, hafnium carbide (HfC) and tantalum carbide (TaC) are among the most refractory binary compounds known to date. Their mixture, with a general formula Ta_xHf_1-xC_y, is known to have a melting point of 4215 K at the composition Ta₄HfC₅, which has long been considered as the highest melting temperature for any solid. Very few measurements of melting point in tantalum and hafnium carbides have been documented, because of the obvious experimental difficulties at extreme temperatures. The investigation lets us identify three major chemical factors that contribute to the high melting temperatures. Based on these three factors, we propose and explore a new class of materials, which, according to our ab initio calculations, may possess even higher melting temperatures than Ta-Hf-C. This example also demonstrates the feasibility of materials screening and discovery via ab initio calculations for the optimization of "higher-level" properties whose determination requires extensive sampling of atomic configuration space.

A DFT study on isomorphously substituted MCM-22 zeolite

The density functional theory has been used to study the isomorphously substituted MCM-22 zeolite for the first time. The effect of the basis sets on the calculation results is discussed in details. Data of several index properties for characterizing the relative acidity of T-MCM-22 (T = B, Al, Ga, and Fe), including proton affinity, bond length and bond angle, OH stretching frequency, and charge on the acidic proton, show that the acidity of T-MCM-22 increases in the sequence of B-MCM-22 < Fe-MCM-22 < Ga-MCM-22 < Al-MCM-22. After making a correction, the calculated OH stretching frequencies for Al-MCM-22 and Fe-MCM-22 show a reasonable agreement with the experimental data. On the basis of an equilibrium structure of the B-MCM-22 zeolite, the effect of the B element in the synthesis of the Ti-MCM-22 is also discussed. The adding of the B element during the synthesis of the Ti-MCM-22 can decrease greatly the Ti substitution energy because of the forming of a structure quite similar to the terminal silanol group. The results can provide some constructively information for zeolite synthesis.

Different electronic structures and spectroscopic properties of cationic [M(ppy)2(N--N)]+ (M = Rh, Ir; N--N = Hcmbpy, H2dcbpy), a DFT study.

We report a comparative quantum-chemical study of the electronic structures and optical properties of a series of cationic complexes [M(ppy)(2)(N--N)](+) (N--N = Hcmbpy = 4-carboxy-4'-methylbpy, M = Rh (Rh1), Ir(Ir2); N--N = H(2)dcbpy = 4,4'-dicarboxy-bpy, M = Rh (Rh3) and Ir (Ir4)). The theoretical calculation reveals that the increased number of -COOH groups on the bpy ligand can decrease the energy levels of LUMO more than HOMO and narrow down the HOMO-LUMO energy gaps, which results in the red-shifted of the lowest-lying absorption and phosphorescent spectra. The lowest-lying singlet absorptions were categorized as d(M,M = Rh or Ir) + pi(ppy) -->pi*(bpy) with MLCT and LLCT characters.

Aminic nitrogen-bearing polydentate Schiff base compounds as corrosion inhibitors for iron in acidic media: A quantum chemical calculation

Quantum chemical calculations based on DFT method were performed on three polydentate Schiff base compounds (PSCs) used as corrosion inhibitors for iron in acid media to determine the relationship between the molecular structure of PSC and inhibition efficiency. The structural parameters, such as the frontier molecular orbital energy HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital), the charge distribution of the studied inhibitors, the absolute electronegativity (chi) values, and the fraction of electrons (Delta N) transfer from inhibitors to iron, were also calculated and correlated with inhibition efficiencies. The results showed that the inhibition efficiency of PSCs increased with the increase in E-HOMO and decrease in E-LUMO-E-HOMO; and the areas containing N atoms are most possible sites for bonding the metal iron surface by donating electrons to the metal. (C) 2007 Elsevier Ltd. All rights reserved.

Calculation of gamma spectra for positron annihilation on molecules

Calculations of gamma spectra for positron annihilation for a selection of molecules, including methane and its fluoro-substitutes, ethane, propane, butane and benzene are presented. The contribution to the ?-spectra from individual molecular orbitals is obtained from electron momentum distributions calculated using the density functional theory (DFT) based B3LYP/TZVP model. For positrons thermalised to room temperature, the calculation, in its simplest form, effectively treats the positron as a plane wave and gives positron annihilation ?-spectra linewidths that are broader (30-40%) than experiment, although the main chemical trends are reproduced. The main physical reason for this is the neglect of positron repulsion from the nuclei. We show that this effect can be incorporated through momentum-dependent correction factors, determined from positron-atom calculations, e.g., many-body perturbation theory. Inclusion of these factors in the calculation gives linewidths that are in improved agreement with experiment.