Strong electron correlations in electronic structure theory are purely quantum effects arising as a result of degeneracies in molecules and materials, and exhibit significantly different yet interesting characters than do weak correlations. Although weak correlations have recently been able to be described very efficiently and accurately within single particle pictures, less known are good prescriptions for treating strong correlations efficiently. Brute-force calculations of strong correlations in wave function theories tend to be very computationally-intensive, and are usually limited to small molecules for applications. Breaking symmetry in a mean-field approximation is an efficient alternative to acquire strong correlations with, in many cases, qualitatively accurate results. The symmetry broken in quantum chemistry has been traditionally of spin, in so-called unrestricted methods, which typically break spatial symmetry as a consequence, and vice versa, in most situations. In this work, we present a novel approach to accurately describing strong correlations with a mean-field cost by means of Hartree- Fock-Bogoliubov (HFB) theory. We are inspired by the number-symmetry-breaking in HFB, which, with an attractive particle interaction, accounts for strong correlations, while maintaining spin and spatial symmetry. We show that this attractive interaction must be restricted to the chemically-relevant orbitals in an active space to obtain physically meaningful results. With such constraints, our constrained pairing mean-field theory (CPMFT) can accurately describe potential energy curves of various strongly-correlated molecular systems, by cleanly separating strong and weak correlations. To achieve the correct dissociation limits in hetero-atomic molecules, we have modified our CPMFT functional by adding asymptotic constraints. We also include weak correlations by combining CPMFT with density functional theory for chemically accurate results, and reveal the connection between CPMFT and traditional unrestricted methods. The similarity between CPMFT and unrestricted methods leads us to the idea of constrained active space unrestricted mean-field approaches. Motivated by CPMFT, we partially retrieve spin-symmetry that has been fully broken in unrestricted methods. We allow symmetry breaking only in an active space. This constrained unrestricted Hartree-Fock (CUHF) is an interpolation between two extrema: the fully broken-symmetry solution and the symmetry preserved solution. This thesis defines the theory behind and reports the results of CUHF. We first show that, if an active space is chosen to include only open-shell electrons, CUHF reduces to restricted open-shell Hartree-Fock (ROHF), and such CUHF proves in many ways significantly

Orientadores: Gaston Eduardo Barberis, Pascoal José Giglio Pagliuso / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Física Gleb Wataghin / Made available in DSpace on 2018-08-17T22:20:28Z (GMT). No. of bitstreams: 1 CalderonFilho_CesarJose_M.pdf: 21461723 bytes, checksum: 3089cf622e6e81e41b460dc218211be3 (MD5) Previous issue date: 2011 / Resumo: Na física da materia condensada, o estudo de sistemas de eletrons fortemente correlacionados é, com certeza, um dos problemas mais interessantes tanto do ponto de vista experimental como teórico, e são estes materiais que tem sido utilizados recentemente em aplicações tecnológicas. Destes compostos, os multiferroicos apresentam um conjunto de propriedades físicas muito rico. Estes materiais apresentam pelo menos duas das seguintes correlações de longo alcance: (anti)ferromagnetismo, ferroelasticidade e ferroeletricidade. Porém, as transições não precisam ser necessariamente correlacionadas, mas quando são, estas ocorrem simultaneamente, e o efeito magnetoelétrico pode ser induzido por campo. Neste trabalho, foram desenvolvidos cálculos numéricos que simulam o acoplamento magnetoelétrico presente nos multiferróicos minimizando a energia através da técnica de Monte Carlo. Foram desenvolvidos dois modelos muito simples. O primeiro modelo acopla uma rede de Ising 2D com spin 1/2 com uma rede de dipolos elétricos tambem 2D; este acoplamento e tal que a mudança de direção de um dado spin reorienta uma dada componente perpendicular do dipolo elétrico vizinho a este mesmo spin. Assim, para este primeiro modelo, as transições de fase das redes elétrica e magnetica ocorrem na mesma temperatura, sendo o hamiltoniano dependente de três parâmetros. Para o segundo modelo, foram utilizadas novamente duas redes, uma rede de Ising 2D com spin 1/2, e uma rede elétrica que se comporta da mesma maneira que uma rede de Ising 2D. Neste caso, o acoplamento entre o spin e o dipolo eletrico ocorre através de um sistema de dois níveis, gerando a possibilidade de temperaturas de transição independentes para as duas redes. Este segundo modelo tambem depende de três parâmetros / Abstract: In condensed matter physics, the study of strongly correlated electron systems is certainly one of the most interesting problems both from the experimental and the theoretical points of view, also these materials recently being used in technological applications. Among these compounds, the multiferroics show a very rich set of physical properties. These materials have at least two of the following long-range correlations: (anti)ferromagnetism, ferroelasticity and ferroelectricity. However, the transitions need not necessarily to be correlated, but when it happens, they occur simultaneously, and the magnetoelectric effect can be induced by field. In this work, numerical calculations have been developed to simulate the magnetoelectric coupling present in the multiferroics minimizing the energy through Monte Carlo technique. Two simple models have been developed. The first model couples a spin 1/2 2D Ising magnetic lattice with to a 2D lattice of classic electric dipoles; this coupling is such that the change in the spin direction reorients a perpendicular component of the electric dipole neighbor of this same spin. Therefore, for this first model, the phase transitions of the magnetic and electric lattices occur at the same temperature, and the Hamiltonian is dependent of three parameters. For the second model, two lattices have been used again, a 2D Ising lattice for the magnetic system and an electric lattice that also behaves as a 2D Ising lattice. In this case, the coupling between the spin and the electric dipole occurs through a two-level system, generating the possibility of the independent transition temperatures for the two systems. This second model also contains three independent parameters / Mestrado / Física da Matéria Condensada / Mestre em Física

This thesis describes the neutron scattering measurements of magnetic excitations in spin-chains and ladders. The first part discusses an experimental investigation of the copper oxide family Sr₁₄Cu₂₄O₄₁ composed of edge-sharing chains and spin-ladders. The study of La₄Sr₁₀Cu₂₄O₄₁ comprises a slightly hole-doped chain and an undoped ladder structure where the chain can be modeled by a ferromagnetic nearest and an antiferromagnetic next-nearest neighbor coupling. The hole effects are apparent in gaps in the dispersion relation and can be described by a charge-density wave agreeing with the commensuration of the dispersion. Investigating the undoped ladder establishes the exchange constants including a cyclic exchange manifested by the two-magnon continuum and the suppression of the S = 1 bound mode. An orbital consideration provides an explanation for the exchanges including the different sizes of rung and leg coupling. The excitation spectrum of the doped ladder in Ca₂.₅Sr₁₁.₅Cu₂₄O₄₁ can be described by a direct comparison with the undoped ladder and the differences consisting of a higher energy mode and subgap scattering can be successfully modeled by the charge spectrum of the ladder calculated from the free electron model. The second part of the thesis investigates the alternating chain material Cu(NO₃)₂ · 2.5D2O and establishes the gapped one-magnon dispersion, the two-magnon continuum and for the first time the S =1 bound mode. Applying magnetic field drives the system through two critical field transitions, condensation of magnons into the ground state and saturation. The modes beyond saturation can be modeled by spin wave theory and the excitations at the first critical field follow Luttinger Liquid behavior. Additionally investigated are the temperature effects with the excitations being of a different nature but containing the signature of a strong correlated system. For an outlook the measurements including temperature and field are provided with further theoretical descriptions necessary.

Complex-oxides have seen an enormous amount of attention in the realm of Condensed Matter Physics and Materials Science/Engineering over the last several decades. Their ability to host a wide variety of novel physical properties has even caused them to be exploited commercially as dielectric, metallic and magnetic materials. Indeed, since the discovery of high temperature superconductivity in the “Cuprates” in the late 1980’s there has been an explosion of activity involving complex-oxides. Further, as the experimental techniques and equipment for fabricating thin films and heterostructures of these materials has improved over the last several decades, the search for new and more exotic properties has intensified. These properties stem from the interfaces formed by depositing these materials onto one another. Whether it be interfacial strain induced by the mismatch between the crystal structures, modified exchange interactions, or some combination of these and other interactions, thin films and heterostuctures provide an invaluable tool the modern condensed matter community. Simply put, a “complex-oxide” is any compound that contains Oxygen and at least two other elements; or one atom in two different oxidation states. Transition Metal Oxides (TMO’s) are a subset of complex-oxides which are of particular interest because of their strong competition between their charge, spin and orbit degrees of freedom. As we progress down the periodic table from 3d to 4d to 5d transition metals, the crystal field, electron correlation and spin-orbit energies become more and more comparable. Therefore, TMO thin films and heterostructures are indispensable to the search for novel physical properties. KTaO3 (KTO) is a polar 5d TMO which has been investigated for its high-k dielectric properties. It is a band insulator with a cubic perovskite crystal structure which is isomorphic to SrTiO3 (STO). This is important because non-polar STO is famous for forming a highly mobile, 2-Dimensional Electron Gas (2DEG) at the hetero-interface with polar LaAlO3 (LAO) as a result of the so-called “polar catastrophe”. Here, I use this concept of polarity to ask an important question: “What happens at hetero-interfaces where two different polar complex oxides meet?” From this question we propose that a hetero-interface between two polar complex-oxides with opposite polarity (I-V/III-III) should be impossible because of the strong Coulomb repulsion between the adjacent layers. However, we find that despite this proposed conflict we are able to synthesize KTO thin films on (110) oriented GdScO3 (GSO) substrates and the conflict is avoided through atomic reconfiguration at the hetero-interface. SrRuO3 (SRO) is a 4d TMO, and an itinerant ferromagnet that is used extensively as an electrode material in capacitor and transistor geometries and proof-of-concept devices. However, in the thin film limit the ferromagnetic transition temperature, TC, and conductivity drop significantly and even become insulating and lose their ferromagnetic properties. Therefore, we ask “Are the transport properties of SRO thin films inherently inferior to single crystals, or is there a way to maintain and/or enhance the metallic properties in the thin film limit?” We have fabricated SRO thin films of various thickness on GSO substrates (tensile strain) and find that all of our samples have enhanced metallic properties and even match those of single crystals. Finally, we ask “Can these enhanced metallic properties in SRO thin films allow us to observe evidence of a topological phase without the complexity of off-stoichiometry and/or additional hetero-structural layers?” Recent reports of oxygen deficient EuO films as well as hetero-structures and superlattices of SRO mixed with SrIrO3 or La0.7Sr0.3MnO3 have suggested that a magnetic skyrmion phase may exist in these systems. By measuring the Hall resistivity, we are able to observer a topological Hall effect which is likely a result of a magnetic skyrmion. We find that of the THE exists in a narrow temperature range and the proposed magnetic skyrmions range in size from 20-120 nm. Therefore, the SRO/GSO system can provide a more viable means for investigating magnetic skyrmions and their fundamental interactions.

The central problem in the modern electronic structure theory is the calculation of cor- relation energy, possibly by an approach that would account for both static and dynamic correlation in an efficient, balanced and accurate way. In this thesis, I present a collection of methods that combine the effective treatment of dynamic correlation by the coupled cluster theory with density matrix renormalization group, a well-established technique for calculations of strongly correlated systems. The connection between them is achieved via the tailored coupled clusters (TCC) ansatz, which conveniently does not impose any ad- ditional computational costs. After the successful initial assessment, we developed more efficient implementations of these methods by employing the local approaches based on pair natural orbitals. This way, we extended the range of possible applications to larger systems with thousands of basis functions. To assess the accuracy of TCC as well as its local counterparts, we performed a variety of benchmark calculations ranging from small, yet challenging systems such as the nitrogen molecule or tetramethyleneethane diradical, to larger molecules like oxo-Mn(Salen) or Fe(II)-porphyrin model. 1

Le présent travail concerne l’avancement des approximations de l’énergie d’échange- corrélation (XC) de la théorie fonctionnelle de la densité (DFT) de Kohn-Sham (KS) basée sur l’approche du facteur de corrélation (CF). Le travail est organisé en trois parties où chaque partie est construite sur des modèles et méthodes précédents. La première partie du travail introduit une nouvelle condition physique à travers la déri- vation du développement en série du quatrième ordre du trou d’échange exact. La dérivation détaillée des formules requises est suivie d’une analyse approfondie qui montre que le terme de quatrième ordre peut ajouter des informations supplémentaires importantes qui sont par- ticulièrement pertinentes pour les molécules par rapport aux atomes. Sur la base de ces résultats, nous explorons les fonctionnelles d’échange qui dépendent du terme de quatrième ordre de l’expansion du trou d’échange. Nous constatons également que les développements d’ensembles de base gaussiens, fréquemment utilisés dans les codes de structure électronique, donnent des représentations insatisfaisantes du terme de quatrième ordre. La deuxième partie de ce travail porte sur la mise en œuvre de nouvelles versions du modèle CF initial [J. P. Precechtelova, H. Bahmann, M. Kaupp et M. Ernzerhof, J. Chem. Phys. 143, 144102 (2015)] dans lequel le trou XC est approximé. Étant donné que diverses contraintes satisfaites par le trou XC sont connues, des approximations peuvent être conçues pour éviter en grande partie des ajustements empiriques. Dans l’approche CF, le trou XC est écrit comme le produit d’un trou d’échange multiplié par un facteur de corrélation. Une contrainte importante satisfaite par le modèle CF est qu’il reproduit correctement l’éner- gie d’échange exacte dans la limite de haute densité. Ceci est réalisé en utilisant l’énergie d’échange exacte par particule comme variable d’entrée, c’est-à-dire que le modèle CF s’ap- puie sur l’échange exact. Des variations du modèle CF initial sont proposées qui assurent que la réponse exacte est obtenue dans la limite homogène. De plus, nous appliquons une correction à la profondeur du trou XC qui est conçue pour capturer une forte corrélation. Les fonctions d’échange-corrélation qui s’appuient sur un échange exact, comme les hybrides, échouent souvent pour les systèmes qui présentent une corrélation électronique importante. Malgré ce fait et malgré la réduction de l’empirisme à un seul paramètre dans CF, des énergies d’atomisation précises sont obtenues pour des composés de métaux de transition fortement corrélés. Le modèle CF montre des résultats significativement supérieurs aux fonctionnelles populaires comme Perdew-Burke-Ernzerhof (PBE), PBE hybride et Tao-Perdew-Staroverov- Scuseria (TPSS). La troisième partie du travail s’appuie sur les modèles CF précédents développés dans notre groupe et aborde l’erreur d’auto-interaction à un électron et introduit un modèle de facteur de corrélation modifié où f C (r, u) est construit tel qu’il se réduit à un dans les régions à un électron d’un système à plusieurs électrons. Ce trou XC avec une correction d’auto- interaction est ensuite utilisé pour générer la fonctionnelle énergie XC correspondante. La nouvelle fonctionnelle est évaluée en l’implémentant dans un programme KS et en calculant diverses propriétés moléculaires. Nous constatons que, dans l’ensemble, une amélioration significative est obtenue par rapport aux versions précédentes du modèle de facteur de cor- rélation. / The present work is concerned with the advancement of approximations to the exchangecorrelation (XC) energy of Kohn-Sham (KS) density functional theory (DFT) based on the correlation factor (CF) approach. The work is organized in three parts where each part is build upon previous models and methods. The first part of the work introduces a new physical condition through the derivation of the fourth-order series expansion of the exact exchange hole. The detailed derivation of the required formulas is followed by a thorough analysis that shows that the fourth-order term can add important additional information that is particularly relevant for molecules compared to atoms. Drawing on these findings, we explore exchange functionals that depend on the fourth-order term of the expansion of the exchange hole. We also find that Gaussian basis set expansions, frequently used in electronic structure codes, result in unsatisfactory representations of the fourth-order term. The second part of this work addresses the implementation of new versions of the initial CF model [J. P. Precechtelova, H. Bahmann, M. Kaupp, and M. Ernzerhof, J. Chem. Phys. 143, 144102 (2015)] in which the XC hole is approximated. Since various constraints satisfied by the XC hole are known, approximations to it can be designed which largely avoid empirical adjustments. In the CF approach, the XC-hole is written as a product of an exchange hole times a correlation factor. An important constraint satisfied by the CF model is that it correctly reproduces the exact exchange energy in the high density limit. This is achieved by employing the exact exchange-energy per particle as an input variable, i.e., the CF model builds on exact exchange. Variations of the initial CF model are proposed which ensure that the exact answer is obtained in the homogeneous limit. Furthermore, we apply a correction to the depth of the XC-hole that is designed to capture strong correlation. Exchangecorrelation functionals that build on exact exchange, such as hybrids, often fail for systems that exhibit sizeable electron correlation. Despite this fact and despite the reduction of empiricism to a single parameter within CF, accurate atomization energies are obtained for strongly-correlated transition metal compounds. The CF model significantly improves upon widely used functionals such as Perdew-Burke-Ernzerhof (PBE), PBE hybrid, and Tao-Perdew-Staroverov-Scuseria (TPSS) density functionals. The third part of the work builds on the previous CF models developed in our group and addresses the one-electron, self-interaction error and introduces a modified correlation factor model where fC(r, u) is constructed such that it reduces identically to one in oneelectron regions of a many-electron system. This self-interaction corrected XC-hole is then used to generate the corresponding XC-energy functional. The new functional is assessed by implementing it into a KS program and by calculating various molecular properties. We find that, overall, a significant improvement is obtained compared to previous versions of the correlation factor model.