This thesis is concerned with the muon-spin relaxation (musr) technique and its application in the study of a number of molecular magnetic systems that may be driven through a quantum phase transition at low temperatures through the application of a magnetic field or hydrostatic pressure. Musr is a highly sensitive probe of magnetism, but its utility can be severely limited by the lack of knowledge of the muon implantation site and the extent to which the muon perturbs its host. In a system of ionic fluorides, where partial information about the muon site is experimentally available, I demonstrate systematically that these problems can be addressed accurately using electronic-structure calculations. The F--$mu$--F complex formed by muons in many fluorides can be understood as an exotic molecule-in-a-crystal defect with a zero-point energy larger than that of any naturally-occurring triatomic molecule. I demonstrate the interesting possibility of controlling the magnetic dimensionality in a molecular magnet using applied pressure. musr and high-field magnetisation experiments under applied pressure on the coordination polymer CuF$_2$(H$_2$O)$_2$(pyrazine) show a transition from a quasi-two-dimensional to a quasi-one-dimensional antiferromagnetic phase. Density-functional theory calculations and calculations of the dipolar anisotropy complement the experiments. I describe how subtle differences in chemical composition can lead to starkly different structural and magnetic properties. [Cu(pyz)(H$_2$O)(gly)$_2$](ClO$_4$)$_2$ may be considered an antiferromagnetic chain that orders below 50 mK while the related compound [Cu(pyz)(gly)](ClO$_4$) is formed from Cu$^{2+}$ dimers and remains disordered down to 30 mK in zero field, but displays a field-temperature phase diagram consistent with the Bose-Einstein condensation of triplons. I also describe musr measurements on the strong-leg spin ladder DIMPY and on the molecular nanomagnets Cr$_8$Cd and Cr$_8$Mn which highlight some of the remaining challenges for longitudinal-field musr experiments.

When molecules are excited by photons or energetic particles, they will cool through the emission of photons, electrons, or by fragmenting. Such processes are often thermal as they occur after the excitation energy has been redistributed across all degrees-of-freedom in the system. Collisions with atoms or ions may also lead to ultrafast fragmentation in Rutherford-like scattering processes, where one or several atoms can literally be knocked out of the molecule by the incoming projectile before the energy can be completely redistributed. The resulting fragmentation pathways can in such knockout processes be very different from those in thermal processes. This thesis covers extensive studies of collisions between ions/atoms and isolated Polycyclic Aromatic Hydrocarbon (PAH) molecules, isolated fullerene molecules, or clusters of these. The high stabilities and distinct fragmentation channels make these types of molecules excellent test cases for characterizing knockout-driven fragmentation and the reactions that these processes can lead to. I will present experimental measurements for a wide range of energies and compare them with my own molecular dynamics simulations and quantum chemical calculations. In this thesis, I present an in-depth study of the role of knockout in the energetic processing of molecules and clusters. The competition between knockout and thermally driven fragmentation is discussed in detail. Knockout-driven fragmentation is shown to result in exotic fragments that are far more reactive than the intact parent molecules or fragments from thermal processes. When such reactive species are formed within molecular clusters efficient molecular growth can take place on sub-picosecond timescales. The cluster environments are crucial here because they protect the newly formed molecules by absorbing excess energy. This is a possible pathway for the growth of large PAHs, fullerenes, and similar carbonaceous complexes found in, for instance, the interstellar medium. / <p>At the time of the doctoral defense, the following paper was unpublished and had a status as follows: Paper 1: Submitted.</p>

Nanoclusters de metais de transição (TM) com dimensões de &#126;1,0 nm têm atraído grande interesse em diversas aplicações tecnológicas, como microeletrônica, óptica, dispositivos magnéticos e principalmente como nanocatalisadores. Nesse contexto, a morfologia desses sistemas é um fator fundamental para otimzação de performace nessas aplicações, tendo em vista a relação direta entre as propriedades físico-químicas e a estrutura atômica do nanocluster determinada pelas condições do ambiente (temperatura e moléculas adsorvidas). No entanto, diante das limitações relativas à precisão das técnicas de caracterização experimental disponíveis atualmente, sobretudo em dimensões nanométricas, o uso de cálculos computacionais através de métodos de primeiros princípios (ab initio) e baseados na teoria do funcional da densidade se torna indispensável. Neste trabalho, foram investigados os efeitos de ligantes e temperaturas nas propriedades estruturais, eletrônicas, estabilidade e meta-estabilidade de nanoclusters de TM. (i) Os efeitos de ligantes de (PH3)n e (SH2)n (saturação gradativa em n = 1, 6, 12 e 18) adsorvidos em nanoclusters unários de Pt55 e Au55 com estruturas icosaédricas (ICO), cuboctaédricas (CUB) e desordenadas de core reduzido (DRC) foram estudados a 0K. Em fase gasosa, as estruturas DRC com 7 e 9 átomos na região do core são 5,34 eV (Pt55) e 2,20 eV (Au55) mais estáveis que modelo ICO com simetria Ih com 13 átomos no core. Os resultados mostraram que existe uma forte compressão do core catiônico pela superfície aniônica induzida por interações de Coulomb (core+-superfície-), levando ao colapso e redução de estresse das estruturas simétricas a partir da redução do número de átomos na região do core. No entanto, a estabilidade da estrutura ICO aumenta com o aumento do número de moléculas adsorvidas, de modo que DRC e ICO se tornam energeticamente degeneradas em < 0,5 eV. Além disso, a adição de ligantes na superfície aniônica reduz a transferência de cargas entre as regiões de core+-superfície-, contribuindo para a redução das interações de Coulomb e, consequentemente, aliviando o estresse interno da estrutura ICO. Resultados similares foram obtidos utilizando ligantes de trifenilfosfina (PPh3), nos quais as longas cadeias carbônicas adicionam interações laterais entre os ligantes. (ii) Nanoclusters de binários de PtnCo(55 - n) a 0K mostraram alta estabilidade em todo intervalo de composições, indicando uma correlação direta entre estabilidade e distribuição homogênea de átomos de Pt e Co formando a nanoliga com estrutura ICO. No entanto, sob uma atmosfera saturada de ligantes de CO adsorvidos, a estabilidade da nanoliga diminui (aumenta) para composições com grandes quantidades de Pt (Co). As análises mostraram que sob os efeitos da adsorção as composições permanecem com configuração ICO, exceto para Pt42Co13 (estrutura core@shell do tipo Co@Pt em fase gasosa), onde as moléculas de CO deslocam os átomos de Co para a superfície, e, então, induzindo um processo de amorfização na formação de uma estrutura formada com átomos de Pt ocupando o core. Para investigar os efeitos de temperatura na estabilidade e meta-estabilidade de nanoclusters e clusters de Au utilizando dinâmica molecular de Born-Oppenheimer, foram utilizados nanoclusters de tamanho médio de Au25, Au38 e Au40 em fase gasosa e clusters de Au13 sob diferentes atmosferas saturadas de CO. (iii) Observou-se que em temperaturas de 300, 400, 500 e 600 K os nanoclusters exibem estruturas dinâmicas para as regiões de core-superfície, com core tetraédrico (4 átomos) ou T-bipiramidal (5 átomos) catiônicos fracamente ligados à superfície aniônica flexível. Uma abordagem estatística através de um algoritmo de redução de dimensionalidades para representação no espaço euclideano bidimensional, chamado de sketch-map, foi proposta como uma nova linguagem para analisar a superfície de energia livre (FES) dos sistemas calculada na aproximação de multistate Bennet acceptance-ratio. A FES foi utilizada como uma análise qualitativa das configurações estáveis e meta-estáveis sob os efeitos de temperaturas, indicando as regiões preferenciais do espaço configuracional de cada nanocluster investigado. (iv) Incluindo os efeitos de ligantes e temperaturas em clusters de Au13, foi observado que em altas temperaturas os clusters de Au13 em fase gasosa tendem a configurações mais abertas com menor comprimento de ligação com relação às estruturas tridimensionais. Por outro lado, a saturação do ambiente com ligantes de CO restringem a mobilidade dos átomos nos clusters de Au13, favorecendo a amostragem de estruturas tridimensionais mesmo em temperaturas relativamente altas. A análise de população de cargas mostrou uma alta concentração de elétrons nos átomos de O, de modo que o forte catiônico dos clusters de Au13 leva ao aumento da coordenação atômica, contribuindo para a estabilização das estruturas mais tridimensionais. / Transition metal nanoclusters (TM) with dimensions of &#126;1,0 nm have attracted great interest in various technological applications such as microelectronics, optics, magnetic devices and mainly as nanocatalysts. In this way, the morphology of these systems is a fundamental factor for optimization of performance in these applications, considering the direct relationship between the physicochemical properties and the atomic structure of the nanocluster determined by the ambient conditions (temperature and adsorbed molecules). In this sense, in view of the limitations on the precision of the experimental characterization techniques currently available, especially in nanometric dimensions, the use of computational calculations using first principles methods (ab initio) and based on the density functional theory is indispensable. Thus, in this work, the effects of binders and temperatures on the structural, electronic properties, stability and metastability of TM nanoclusters were investigated. (i) In this way, the (PH3)n and (SH2)n ligands effects with gradual saturation at n = 1, 6, 12 and 18 adsorbed onto nanoclusters Pt55 and Au55 with icosahedral (ICO), cuboctahedron (CUB), and distorted reduced core (DRC) structures were studied at 0K. In the gas phase, the DRC structures with 7 and 9 atoms in the core region are 5,34 eV (Pt55) and 2,20 eV (Au55) more stable than ICO model with symmetry Ih with 13 atoms in core. The results showed that there is a strong compression of the cationic core by the anionic surface induced by interactions of Coulomb (core+-surface-), leading to collapse and stress reduction of the symmetrical structures from the reduction of the number of atoms in the core region. However, the stability of the ICO structure increases with increasing Number of molecules adsorbed, so that DRC and ICO become approximately degenerated in energy in < 0,5 eV. In addition, ligands on the anionic surface reduces the transfer of charges between core+- surface-, contributing to the reduction of interactions of Coulomb and, consequently, alleviating the internal stress of the ICO. Similar result were obtained using triphenylphosphine (PPh3) as large ligands, in which the long carbon chains add side interactions between the ligands. (ii) PtnCo(55 - n) binares nanoclusters at 0K showed a high stability across the range of compositions, indicating a direct correlation between stability and homogeneous distribution of Pt and Co atoms forming the nanoalloys with ICO structure. However, under a saturated atmosphere of adsorbed CO ligands, the stability of the nanoaaloys decreases (increases) to compositions with large amounts of Pt (Co). Analyzes have shown that under the effects of adsorption the compositions remain with ICO configuration, except for Pt42Co13 (Co@Pt as core@shell in gas-phase), where Co molecules displace the Co atoms to the surface, and then inducing an amorphization process for a structure formed with Pt atoms occupying the core. To investigate the temperature effects on the stability and metastability of Au nanoclusters and clusters by using Born-Oppenheimer molecular dynamics, medium size nanoclusters of Au25, Au38, and Au40 were used in gas-phase and clusters of Au13 under different saturated atmospheres of CO. (iii) It has been observed that at temperatures of 300, 400, 500, and 600 K, the nanoclusters exhibit dynamic structures for the core-surface regions with cationic tetrahedral (4 atoms) or T-bipyramidal (5 atoms) cores weakly bound to the flexible anionic surface. A statistical approach using a dimensionality reduction algorithm for two-dimensional Euclidean space representation, called sketch-map, was proposed as a new language to analyze the free energy surface (FES) of the systems calculated on the approximation of multistate Bennet acceptance-ratio. The FES was used as a qualitative analysis of the stable and metastable configurations under the effect of temperature, indicating the preferred regions of the configurational space of each nanocluster investigated. (iv) By including the effects of ligands and temperatures on Au13 clusters, it was observed that at high temperatures the Au13 clusters in gas-phase tend to have more open configurations with less bond length relative to three-dimensional structures. On the other hand, the saturation of the environment with CO ligands restricts the mobility of the atoms in Au13 clusters, favoring the sampling of three-dimensional structures even at temperatures relatively high. The charge population analysis showed a high concentration of electrons in the O atoms, so that the strong cationic character of the Au13 clusters leads to increased atomic coordination, contributing to the stabilization of the more three-dimensional.

A aproximação da densidade local para a teoria do funcional da densidade tem, no passado, levado a resultados conflitantes para a dimerização de cadeias infinitas do trans-poliacetileno (trans-PA). Estes resultados variaram deste forte dimerização, próximo dos resultados experimentais, até fraca dimerização, ou até mesmo nenhuma dimerização. Desde que a aproximação da densidade local usualmente descreve transições de fase estruturais corretamente, esta situação insatisfatória no caso do trans-PA clama por uma investigação detalhada. Neste trabalho estudamos o problema descrevendo a molécula como um modelo de Hubbard-Peierls unidimensional. Estabelecemos uma teoria do funcional da densidade e construímos uma aproximação da densidade local para este modelo. Em acordo com os cálculos ab initio encontramos que este tipo de aproximação não descreve adequadamente a dimerização. Propomos então um Formalismo alternativo, baseado em funcionais da matriz-densidade. Neste formalismo a energia de troca-correlação é escrita como um funcional da densidade de carga (como no método tradicional) e do parâmetro de ordem para a transição. Desta maneira obtemos um tratamento aperfeiçoado para a fase dimerizada. Nossos resultados sugerem que uma descrição adequada de sistemas espacialmente correlacionados dentro da teoria do funcional da densidade requer uma nova classe de funcionais, que vai além da aproximação da densidade local por levar explicitamente em conta parâmetros de ordem espaciais de longo alcance. / The local-density approximation to density-functional theory has, in the past, lead to conflicting results for the dimerization of infinite trans-polyacetylene (trans-PA) chains. These results range from strong dimerization, close to the one observed experimentally, to weak dimerization, or even no dimerization at all. Since the local-density approximation usually describes structural phase transitions correctly, this unsatisfactory situation in the case of trans-PA calls for a detailed investigation. In this work we study the problem by describing the polyacetylene molecule as a one-dimensional Hubbard-Peierls model. We set up a density-functional theory for this model, and construct a local-density approximation. In agreement with the ab initio calculations we find that this type of approximation does not consistently describe the dimerization. We therefore propose an alternative formalism, based on density-matrix functionals. In this formalism the exchange-correlation energy is written as a functional of the charge density (as in the traditional method) and the order parameter for the dimerization transition. In this way we achieve an improved treatment of the dimerized phase. Our results suggest that a reliable description of spatially correlated systems within density-functional theory requires a new type of functionals, going beyond the local-density approximation by explicitly accounting for long-range spatial order.

A Teoria do funcional da Densidade (DFT) é muito empregada no estudo da densidade eletrônica e energia do estado fundamental de sistemas interagente de muitos elétrons. Uma de suas desvantagens é que, apesar de formalmente exata, a DFT depende de aproximações no funcional de troca e correlação Exc[n]. Uma de suas vantagens é a possibilidade de trabalhar com a conexão adiabática, que permite conectar explicitamente o sistema interagente de interesse com uma versão não-interagente de mesma densidade. Baseado em seu escopo esta dissertação tem como objetivo tratar da DFT no modelo de Hubbard unidimensional. Uma das vantagens em trabalhar com o Hubbard é a existência de uma solução formalmente exata para a energia do modelo homogêneo e unidimensional via ansatz de Bethe, tal solução serve de base para os nossos cálculos e resultados. Outra vantagem do modelo de Hubbard é a existência de um gap que aproximações usuais na DFT falham em reproduzir. Um de nossos resultados é calcular o gap de duas formas. Uma através da derivada da solução via ansatz de Bethe e outro é tratar o Hamiltoniano pela conexão adiabática, calcular o gap pelo potencial químico, este último calculado pelo ansatz de Bethe. Ao final comparamos a precisão dos métodos, apenas para descobrir que o método via potencial químico é mais preciso. Por último usamos a conexão adiabática e a aproximação de Hartree-Fock para teorizar um funcional híbrido no modelo de Hubbard unidimensional. / The Density Functional Theory (DFT) is a widely used in the study of electronic density and energy of the ground-state of interacting systems consisting of many electrons. One of its disadvantages is that, although it is formally exact, the DFT depends on approximations of the exchange-correlation functional Exc[n]. One of its advantages consists of the possibility of working with the adiabatic connection, which allows a explicit connection between the interacting system of our interest and a non-interacting system that yields the same density. Based on its scope this dissertation aims to address the DFT in the one-dimensional Hubbard model. One of the Hubbards model advantage consists of the existence of a formally exact solution to the energy of the homogeneous one-dimensional model via Bethe ansatz, such a solution serves as the basis for our calculations and results. Anothe Hubbards model advantage is the existence of a gap that usual approximations in DFT fail to calculate. One of our results is to calculate the gap in two different ways. One through the derivative of the Bethe ansatz solution and the other is to address the Hamiltonian through the adiabatic connection, calculate the gap through the chemical potential, the latter calculated via Bethe ansatz. In the end we compare their precisions, only to find that the one via chemical potential is more precise. At last we use the adiabatic connection and the Hartree-Fock approximation to theorize a hybrid functional in the one-dimensional Hubbard model.

Neste trabalho utilizamos métodos computacionais ab initio, baseados na Teoria do Funcional da Densidade (DFT), para simular em nível atomístico propriedades estruturais, eletrônicas e magnéticas de nanoestruturas de grafeno. Estudamos nanoflocos de grafeno (GNFs) em estado pristino e GNFs com defeitos intrínsecos (monovacância, divacância e Stone-Wales). Escolhemos GNFs com diferentes terminações e formas, e estudamos também empilhamentos duplos - biflocos - em diferentes composições. Empregamos dois enfoques diferentes de DFT, a aproximação de gradiente generalizado simples no nível teórico de Perdew-Burke-Ernzerhof (PBE), e PBE híbrida (PBEh), incorporando uma fração de troca de Hartree-Fock. Todos os cálculos foram realizados através do código all-electron AIMS, incluindo correções de van der Waals. Nossos GNFs foram escolhidos com simetrias específicas: D2h, D3h e D6h, e com diferentes bordas, armchair (AC), zigue-zague (ZZ) e misturas das duas. Os flocos hexagonais D6h apresentam um gap de energia e nao apresentam spin, enquanto flocos perfeitos com bordas zigue-zague e mistas apresentam spin intrínseco. Esse spin não nulo é devido à diferença no numero de átomos entre uma e outra subrede do grafeno (Liebs imbalance). Defeitos em materiais de carbono sao frequentes, e tem sido estudados experimental e teoricamente. Aqui, estudamos a monovaçancia, através de modelos de cluster e supercélulas, e obtemos para esse defeito o momento magnético de = 2B, (B ´e o magneton de Bohr). Mostramos que as diferenças entre resultados anteriores são oriundas do erro de auto-interação presente na DFT simples, amenizado através do uso de PBEh. Através da mesma metodologia estudamos a interação de nanoestruturas de grafeno com moléculas de água, focalizando em propriedades estruturais. A grafite é um material hidrofóbico, mas a nanoestrutura poderia favorecer a interação com a água. Obtemos que pequenos agregados de água são adsorvidos na superfície de GNFs e biflocos, entretanto a inclusão desses agregados na região interna dos biflocos é altamente desfavorável. Assim podemos esperar que essas nanoestruturas empilhadas sejam também hidrofóbicas. / In this work, computational ab initio methods based on density functional theory (DFT) are used to simulate on an atomistic level the structural, electronic and magnetic properties of graphene nanostructures. We study pristine graphene nanoflakes (GNFs), and GNFs with intrinsic defects (monovacancy, divacancy, Stone-Wales). We design GNFs with different terminations and shapes and also studied stacked forms -biflakes- in different compositions. We employed two DFT approaches, plain generalized gradient approximation in the Perdew-Burke-Ernzerhof (PBE) level of the theory, and hybrid PBE (PBEh) incorporating a fraction of Hartree-Fock exchange. All calculations were performed with the all-electron code AIMS, including van der Waals corrections. Our GNFs were chosen from three symmetry groups: D2h, D3h and D6h, and with different edges, armchair (AC), zigzag (ZZ) and a mixture of both. Our chosen D6h- hexagonal flakes present an energy gap and no spin, while perfect trigonal zigzag and mixed edges GNFs have an intrinsic spin. This non-zero spin is due to the graphene sublattice imbalance (Liebs imbalance). Defects are common in carbon materials, and have been experimentally and theoretically studied in graphene. Here, the single vacancy in graphene was studied, by cluster and supercell approaches, finding that the vacancy induces a magnetic moment = 2B (Bohr magneton). We show that conflicting results for the magnetic moment coming from theoretical studies come from the self-interaction error present in plain PBE, cured through the use of PBEh. Using the same methodology we studied the interaction of carbon nanostructures with water molecules, focusing on structural properties. Graphite is a hydrophobic material but nanostructuring could favor the interaction with water. We obtained that small water groups are adsorbed on the surface of GNFs and biflakes, however the inclusion of these groups in the internal region of biflakes is highly unfavorable, thus we can expect these stacked nanostructures to be also hydrophobic.

A escolha do catalisador no processo de reforma a vapor de etanol é de suma importância na determinação da eficiência de produção de hidrogênio, utilizado em células combustíveis. Superfícies de metais de transição (MT) são normalmente utilizadas como catalisadores, porém resultados experimentais têm apresentado que clusters MT suportados em óxidos são mais eficientes e mais seletivos que seus análogos macroscópicos. Assim, para se compreender o mecanismo de interação de etanol e água, e também do radical hidroxila (simples caracterização da ligação O-H), com clusters mágicos MT13 (MT = Ni, Cu, Pd, Ag, Pt e Au), foram utilizados cálculos ab initio baseados na teoria do funcional da densidade com o funcional proposto por Perdew-Burke-Ernzerhof (PBE) e a aproximação de Tkatchenko-Scheffler para as correções de van der Waals (vdW). Os orbitais de Kohn-Sham foram expandidos através de orbitais numéricos centrados nos átomos (NAOs), que estão implementados no pacote Fritz Haber Institute - ab initio molecular simulations (FHI-aims). Foi encontrado que clusters 3d e 4d apresentam configurações mais compactas, enquanto os sistemas de Pt13 e Au13 adotam geometrias mais abertas. Em particular, Au13 possui uma estrutura planar como configuração mais estável, enquanto o seu primeiro isômero mais energético possui uma estrutura 3D. Para os cálculos de adsorção foram selecionados os clusters de menor energia (LOW), e os clusters icosaédricos (ICO). Observou-se que a hidroxila (OH) não possui sítio preferencial ou dependência geométrica durante adsorção, enquanto água e etanol se ligam preferencialmente através do oxigênio no sítio top. Para o caso do etanol também é possível se encontrar isômeros com energias mais altas, na ordem de kT, que se ligam via hidrogênio do grupo CH e que são prováveis a temperatura ambiente. OH possui um elétron desemparelhado, o que favorece a formação de ligações químicas entre a molécula e MT13, de forma que a energia de adsorção varia entre -4.11 e -2.94 eV, condizente com a transferência de carga do cluster para a molécula. Quando se compara a adsorção de água (etanol) sobre superfícies MT(111) e MT13, nota-se que a energia de adsorção, que segue a ordem 3d &gt; 4d &gt; 5d para MT(111), é maior em relação aos clusters e varia entre -255 meV (-317 meV) &lt; -Ead &lt; -670 meV (-837 meV); porém a magnitude da interação ocorre na ordem 3d &gt; 4d, mas 4d &lt; 5d. A quebra da tendência pode ser explicada através dos efeitos de tamanho, que potencializam a atividade catalítica dos clusters de platina e ouro. É possível se observar ainda que as moléculas interagem mais fortemente com os MT13 dos metais com orbitais d semi-preenchidos, uma vez que esses proporcionam um maior rearranjo da densidade eletrônica. Ocorre uma pequena transferência de carga das moléculas de água e etanol para os clusters, apontando para a fraca interação. Cálculos vibracionais suportam as modificações estruturais e energéticas que ocorrem no sistema, onde &nu;O-H aumenta para a hidroxila e diminui para água e etanol. / The choice of the catalyst in steam reforming of ethanol process has a central hole in determining the hydrogen production efficiency, used in fuel cells. Transition-metal (TM) surfaces are commonly used as catalysts, but experimental results have shown that TM clusters supported on oxides are more efficient and more selective than their macroscopic analogues. Therefore, to understand the mechanism of interaction between ethanol and water, and also the hydroxyl radical (simple characterization for the O-H bonding) with TM13 magic clusters (TM = Ni, Cu, Pd, Ag, Pt and Au), we performed ab initio calculations based on density functional theory within the functional proposed by Perdew-Burke-Ernzerhof (PBE) and Tkatchenko-Scheffler approach for van der Waals corrections (vdW). The Kohn-Sham orbitals were expanded by numerical atom-centered orbitals (NAOs), which are implemented at the Fritz Haber Institute - ab initio molecular simulations package (FHI-AIMS). We found that clusters 3d and 4d prefer more compact configurations, while Pt13 and Au13 systems adopt more open geometries. In particular case of Au13, our lowest energy configuration is a planar structure, while the first high energy isomer is a 3D configuration. For adsorption calculations, we selected the lowest energy configurations (LOW), and the icosahedral clusters (ICO). We observed that the hydroxyl (OH) has no adsorption site preference and also it does not show a geometric trend, while water and ethanol preferentially bind via oxygen at the top site. In the case of ethanol, it is also possible to find isomers with higher energies, which bind via hydrogen from CH group, and they are probable at room temperature. OH has an electron unpaired, which favors the formation of chemical bonds between the molecule and TM13, so that the adsorption energy varies from -4.11 to -2.94 eV, consistent with charge transfer from cluster to the molecule. When we compare the water and ethanol adsorption on surfaces TM(111) and on TM13, we note that the adsorption energy, that follows the order 3d &gt; 4d &gt; 5d for MT(111), is higher compared to the clusters energy. For water (ethanol), this property varies between -255 meV (-317 meV) &lt; -Ead &lt; -670 meV (-837 meV); however, the magnitude of the interaction occurs in the order 3d &gt; 4d, but also 4d &lt; 5d. The breaking of the trend can be explained by size effects, that enhances the catalytic activity of platinum and gold clusters. It is also observed that the stronger interaction between molecules and TM13 with d-orbitais partly occupied, since these metals provide a higher rearrangement of the electron density. There is a small charge transference from water and ethanol molecules to the clusters, pointing to the weak interaction. Vibrational calculations support the structural and energetic changes that occur in the system, where &nu;O-H increases for the hydroxyl and decreases to water and ethanol.

Neste trabalho estudamos o magnetismo no modelo de Hubbard em uma dimensão (1DMH) via Teoria do Funcional da Densidade (DFT). Dentro desse contexto, abordamos três tópicos: (i) Baseados numa recente proposta de generalização da DFT para magnetismo não colinear, estudamos numericamente por meio de diagonalização exata o comportamento magnético do 1DMH adaptado para o caso não colinear. (ii) Desenvolvemos e aplicamos o formalismo da DFT para magnetismo não colinear no espaço discreto e comparamos os resultados do cálculo do sistema Kohn-Sham (KS) com os obtidos no passo (i) através da diagonalização exata. Aqui testamos uma proposta de aproximação para o funcional de troca e correlação no estado de onda de densidade de spin (SDW). (iii) Construímos e testamos uma aproximação da densidade local (LDA) para o funcional de troca e correlação baseados na solução exata do modelo de Hubbard homogêneo, dada pelo Ansatz de Bethe (BA). Comparamos os resultados do cálculo numérico oriundos da BA-LDA com os provenientes da diagonalização exata, obtendo resultados satisfatórios tanto para o estado SDW como para vários sistemas não homogêneos (impurezas, super-redes, potencial binário e potencial degrau). O bom acordo entre a diagonalização exata e a DFT mostra que esta última pode ser aplicada com sucesso nos casos em que a diagonalização exata é inviável. Paralelamente, adaptamos o formalismo da DFT-SDW, usado no passo (i) para aplicação em cálculos ab initio de estrutura de bandas na aproximação Muffin- Tin. Objetivando, nesse úl¬timo contexto, desenvolver as equações matemáticas necessárias para a implementação numérica desse método, fornecendo assim o formalismo necessário para futuros cálculos ab initio / In this work we study the magnetism of the Hubbard Model in one dimension (1DMH) via the Density Functional Theory (DFT). In this context, we considered three different phys¬ical context topics: (i) based on a recent generalization proposed to DFT for non-collinear magnetism, we performed an exact numerical diagonalization in order to study the magnetic behavior of the 1DMH; (ii) we developed and applied the DFT formalism, using the Kohn-Sham (KS) scheme, for non-collinear magnetism in the discreet space and we compared the results with those obtained in the step (i) through exact diagonalization. In that approach we tested an approximation for exchange and correlation functional in the framework of the spin density wave state (SDW); (iii) finally, we built and tested a local density approximation (LDA) for the functional of exchange and correlation based on the exact solution of the homogeneous Hubbard model, given by the Bethe Ansatz (BA). We compared the results of BA-LDA obtained from the numerical calculation using the exact diagonalization, yielding satisfactory results for the SDW state, as well as several inhomogeneous systems (impurity, super-lattice, binary and step potential). The good agreement between exact diagonalization and DFT shows that this last can be applied with success in the cases the diagonalization exact is not viable. We also adapted the DFT-SDW formalism, used in the step (i), for applications in ab initio calculations of band structure in the Muffin-Tin approximation. Our goal, in this last context, is to develop the necessary mathematical equations for the numerical implementation of this method, thus provinding the necessary formalism for future ab initio calculations.

Ligações de hidrogênio é um tema que tem despertado o interesse da comunidade científica desde o final do século XIX. Sua importância é enorme nos processos ligados à vida, como por exemplo na estabilização das estruturas de DNA e na manutenção da água em seu estado líquido. Várias metodologias teóricas foram desenvolvidas para o estudo de sistemas moleculares e das ligaçõas de hidrogênio, entre elas está o emprego de cálculos de teoria de perturbação de muitos corpos (MBPT). Uma alternativa aos cálculos moleculares com MBPT, que tem crescido em termos de aplicação e confiança, é o emprego da teoria do funcional da densidade (DFT). Nesta tese, calculamos propriedades de sistemas hidrogênio-ligados em clusters e líquidos, usando ambas as metodologias DFT e MBPT. Entre as propriedades consideradas estão constantes rotacionais, momentos de dipolo, energias de ligação, deslocamentos espectroscópicos quando da formação do complexo e espalhamento de luz. Parte desta tese é dedicada a salientar as diferenças entre as propriedades de um cluster otimizado e a estrutura de um líquido gerada por simulação de Monte Carlo. Comparamos os resultados obtidos para o complexo uréia-água nestas duas situações e reforçamos o fato de que líquido e cluster são situações físicas distintas, cujas propriedades também são diferentes. Os sistemas estudados foram HCN, CH IND. 3CN, HC IND. 3N, HC IND. 2NC, HCN...H IND. 2O, CH IND. 3CN...H IND.2O, (CH IND. 3) IND.3CCN...H IND. 2O e (NH IND. 2) IND. 2CO...H IND. 2O. Dos resultados obtidos nesta fase, verificamos que DFT é de fato uma alternativa completamente viável para a obtenção de propriedades de moléculas e biomoléculas hidrogênio-ligadas. / Hydrogen bonding is a topic of interest in the scientific community since the end of the XIX century. Its importance is enormous in processes related to life as, for example, the stabilization of DNA structures and the maintenance of water in its liquid state. Several theoretical methodologies were developed to study molecular systems and hydrogen bonds, among them is the use of many-body perturbation theory (MBPT). An alternative to MBPT, that has gained confidence, is the employment of the density functional theory (DFT). In the present thesis we calculate properties of hydrogen-bonded systems in clusters and liquids using both methodologies, DFT and MBPT. Among the properties considered are rotational constants, dipole moments, biding energies, spectroscopic shifts upon complex formation and light scattering. Part of this thesis is dedicated also to point out the difference between the properties of an optimized cluster and a liquid structure generated by Monte Carlo simulation. We compare the results obtained for the urea-water complex in these two situations and reinforce the fact that liquid and cluster are different physical situations, whose properties are also different. The systems studied were HCN, CH IND.3CN, HC IND.3N, HC IND.2NC, HCN H IND.2O and (NH IND.2)IND.2CO H IND.2O. from the results obtained in this thesis, we verify that DFT is indeed feasible to obtain properties of hydrogen bonded molecules and biomolecules.

Le déploiement constant de nouveaux réseaux de capteurs dans des endroits confinés des turbomachines ou difficilement accessibles nécessite l'apport d'énergie pour alimenter ces capteurs. De nombreuses recherches s'orientent donc sur une alimentation à demeure des capteurs afin de monitorer de nouveaux composants. Une voie possible d'alimentation (à demeure) est l'utilisation de modules thermoélectriques afin de convertir l'énergie thermique issue de la combustion dans les turbines. Parmi les matériaux thermoélectriques utilisables dans ce genre d'applications (au-delà de 700 °C), les siliciures semblent constituer la famille la plus prometteuse. Le sujet de cette thèse est l'étude de siliciures connus pour leurs résistances à l'environnement en température afin de réaliser des modules thermoélectriques pouvant fonctionner à des températures de l'ordre de 700 °C. À cette fin, différents représentants de la famille des siliciures ont été étudiés par calcul ab initio à l'aide de la théorie de la fonctionnelle de la densité. Cet outil a permis de calculer leurs propriétés thermoélectriques potentielles et a montré que les deux meilleurs candidats à des applications bon marché étaient les disiliciures de manganèse et de fer. Le calcul a également permis de montrer le caractère métallique de nombreux siliciures ternaires. Les siliciures prometteurs ont été élaborés par frittage et leur vieillissement ainsi que les coefficients de dilatation thermique ont été étudiés. Ces connaissances ont permis de mettre au point des modules thermoélectriques à base de siliciures résistants à 900 °C sous air. Bien que possédant des propriétés thermoélectriques modestes, leur résistance à l'oxydation à haute température permet d'envisager la fabrication de modules thermoélectriques / The constant deployment of new sensors networks in confined areas of turbomachines or difficult to access, requires the input of energy to power these sensors. Many researches are thus focused on a permanent power supply of sensors to monitor new components. One possible way of permanent power supply is the use of thermoelectric modules to convert the thermal energy from combustion into the turbines. Among the thermoelectric materials, silicides seem to be the most promising family for high temperature applications (above 700 °C). The subject of this thesis is the study of silicides known for their high temperature oxidation resistance in order to produce thermoelectric modules that can operate at optimal temperatures around 700 °C. To this end, different representatives of the silicide family have been studied by ab initio calculation using the density functional theory. This tool allowed to calculate their potential thermoelectric properties and showed that the two best candidates for cheap applications were the disilicides of manganese and iron. The calculation also showed the metallic character of many ternary silicides. The promising silicides were developed by sintering method and their aging as well as their coefficients of thermal expansion were studied. This knowledge has made it possible to develop thermoelectric modules based on silicides resistant to 900 °C in air. Although possessing modest thermoelectric properties, their resistance to high temperature oxidation makes it possible to envisage the manufacture of thermoelectric modules