Highly accurate quantum chemistry : spin-orbit splittings via multireference coupled-cluster methods and applications in heavy-atom main-group chemistry

Diese Dissertation demonstriert und verbessert die Vorhersagekraft der Coupled-Cluster-Theorie im Hinblick auf die hochgenaue Berechnung von Moleküleigenschaften. Die Demonstration erfolgt mittels Extrapolations- und Additivitätstechniken in der Single-Referenz-Coupled-Cluster-Theorie, mit deren Hilfe die Existenz und Struktur von bisher unbekannten Molekülen mit schweren Hauptgruppenelementen vorhergesagt wird. Vor allem am Beispiel von cyclischem SiS_2, einem dreiatomigen Molekül mit 16 Valenzelektronen, wird deutlich, dass die Vorhersagekraft der Theorie sich heutzutage auf Augenhöhe mit dem Experiment befindet: Theoretische Überlegungen initiierten eine experimentelle Suche nach diesem Molekül, was schließlich zu dessen Detektion und Charakterisierung mittels Rotationsspektroskopie führte. Die Vorhersagekraft der Coupled-Cluster-Theorie wird verbessert, indem eine Multireferenz-Coupled-Cluster-Methode für die Berechnung von Spin-Bahn-Aufspaltungen erster Ordnung in 2^Pi-Zuständen entwickelt wird. Der Fokus hierbei liegt auf Mukherjee's Variante der Multireferenz-Coupled-Cluster-Theorie, aber prinzipiell ist das vorgeschlagene Berechnungsschema auf alle Varianten anwendbar. Die erwünschte Genauigkeit beträgt 10 cm^-1. Sie wird mit der neuen Methode erreicht, wenn Ein- und Zweielektroneneffekte und bei schweren Elementen auch skalarrelativistische Effekte berücksichtigt werden. Die Methode eignet sich daher in Kombination mit Coupled-Cluster-basierten Extrapolations-und Additivitätsschemata dafür, hochgenaue thermochemische Daten zu berechnen.

Quantum effects and dynamics in hydrogen-bonded systems: a first-principles approach to spectroscopic experiments

Computer simulations have become an important tool in physics. Especially systems in the solid state have been investigated extensively with the help of modern computational methods. This thesis focuses on the simulation of hydrogen-bonded systems, using quantum chemical methods combined with molecular dynamics (MD) simulations. MD simulations are carried out for investigating the energetics and structure of a system under conditions that include physical parameters such as temperature and pressure. Ab initio quantum chemical methods have proven to be capable of predicting spectroscopic quantities. The combination of these two features still represents a methodological challenge. Furthermore, conventional MD simulations consider the nuclei as classical particles. Not only motional effects, but also the quantum nature of the nuclei are expected to influence the properties of a molecular system. This work aims at a more realistic description of properties that are accessible via NMR experiments. With the help of the path integral formalism the quantum nature of the nuclei has been incorporated and its influence on the NMR parameters explored. The effect on both the NMR chemical shift and the Nuclear Quadrupole Coupling Constants (NQCC) is presented for intra- and intermolecular hydrogen bonds. The second part of this thesis presents the computation of electric field gradients within the Gaussian and Augmented Plane Waves (GAPW) framework, that allows for all-electron calculations in periodic systems. This recent development improves the accuracy of many calculations compared to the pseudopotential approximation, which treats the core electrons as part of an effective potential. In combination with MD simulations of water, the NMR longitudinal relaxation times for 17O and 2H have been obtained. The results show a considerable agreement with the experiment. Finally, an implementation of the calculation of the stress tensor into the quantum chemical program suite CP2K is presented. This enables MD simulations under constant pressure conditions, which is demonstrated with a series of liquid water simulations, that sheds light on the influence of the exchange-correlation functional used on the density of the simulated liquid.

Higher-order molecular properties and excitation energies in single-reference and multireference coupled-cluster theory

Coupled-cluster (CC) theory is one of the most successful approaches in high-accuracy quantum chemistry. The present thesis makes a number of contributions to the determination of molecular properties and excitation energies within the CC framework. The multireference CC (MRCC) method proposed by Mukherjee and coworkers (Mk-MRCC) has been benchmarked within the singles and doubles approximation (Mk-MRCCSD) for molecular equilibrium structures. It is demonstrated that Mk-MRCCSD yields reliable results for multireference cases where single-reference CC methods fail. At the same time, the present work also illustrates that Mk-MRCC still suffers from a number of theoretical problems and sometimes gives rise to results of unsatisfactory accuracy. To determine polarizability tensors and excitation spectra in the MRCC framework, the Mk-MRCC linear-response function has been derived together with the corresponding linear-response equations. Pilot applications show that Mk-MRCC linear-response theory suffers from a severe problem when applied to the calculation of dynamic properties and excitation energies: The Mk-MRCC sufficiency conditions give rise to a redundancy in the Mk-MRCC Jacobian matrix, which entails an artificial splitting of certain excited states. This finding has established a new paradigm in MRCC theory, namely that a convincing method should not only yield accurate energies, but ought to allow for the reliable calculation of dynamic properties as well. In the context of single-reference CC theory, an analytic expression for the dipole Hessian matrix, a third-order quantity relevant to infrared spectroscopy, has been derived and implemented within the CC singles and doubles approximation. The advantages of analytic derivatives over numerical differentiation schemes are demonstrated in some pilot applications.

Prediction of new crystal structures under extreme conditions

One of the most important challenges in chemistry and material science is the connection between the contents of a compound and its chemical and physical properties. In solids, these are greatly influenced by the crystal structure.rnrnThe prediction of hitherto unknown crystal structures with regard to external conditions like pressure and temperature is therefore one of the most important goals to achieve in theoretical chemistry. The stable structure of a compound is the global minimum of the potential energy surface, which is the high dimensional representation of the enthalpy of the investigated system with respect to its structural parameters. The fact that the complexity of the problem grows exponentially with the system size is the reason why it can only be solved via heuristic strategies.rnrnImprovements to the artificial bee colony method, where the local exploration of the potential energy surface is done by a high number of independent walkers, are developed and implemented. This results in an improved communication scheme between these walkers. This directs the search towards the most promising areas of the potential energy surface.rnrnThe minima hopping method uses short molecular dynamics simulations at elevated temperatures to direct the structure search from one local minimum of the potential energy surface to the next. A modification, where the local information around each minimum is extracted and used in an optimization of the search direction, is developed and implemented. Our method uses this local information to increase the probability of finding new, lower local minima. This leads to an enhanced performance in the global optimization algorithm.rnrnHydrogen is a highly relevant system, due to the possibility of finding a metallic phase and even superconductor with a high critical temperature. An application of a structure prediction method on SiH12 finds stable crystal structures in this material. Additionally, it becomes metallic at relatively low pressures.

Neue Materialien für die Magnetoelektronik: Heusler- und Halb-Heusler-Phasen

Heusler- und Halb-Heusler-Phasen konnten als Verbindungsklasse identifiziert werden, neue Materialien für die Magnetoelektronik bereitzustellen. Auf Basis eines theoretischen Modells konnte das Auftreten eines hohen MR-Effekts auf bestimmte Valenzelektronenzahlen (VEK) präzisiert werden. Dazu muß ein dreiteiliger 'Fingerabdruck' in der Bandstruktur (Sattelpunkt, Spindichtewelle, lokales magnetisches Moment) vorliegen. Es existieren eine Vielzahl von halbmetallischen Ferromagneten in dieser Verbindungsklasse. Die variable Valenzelektronenkonzentrationen, die sich auch aus den hochsymmetrischen Strukturen ergeben, erlauben eine gute Dotierbarkeit der Phasen. Dadurch ist die Möglichkeit gegeben, die Phase exakt mit der geforderten VEK zu synthetisieren. Curietemperaturen > 500 K sind in Hinblick auf die technische Anwendbarkeit notwendig. PdMnTe hat nicht die Voraussetzungen für einen PMR-Effekt in der Bandstrukturrechnung, doch die Nähe zu einer halbmetallischen Zustandsdichte resultiert in einen negativen CMR-Effekt unterhalb des magnetischen Übergangs von MR0 = 18 % bei 4 K. Die Zusammenhänge von Probenpräparation zum magnetischen Sättigungsmoment konnten an Co2CrAl aufgedeckt werden. Die unter Anwendung des vorgestellten Modells synthetisierte Heusler-Phase Co2Cr0.6Fe0.4Al (VEK = 27.8) weisen den erwarteten MR-Effekt auf. Der gemessene PMR-Effekt ist größer als bei den GMR-Systemen (bei geringerer Feldempfindlichkeit) und anderen granularen Materialien wie CrO2 bei 295 K. Co2Cr0.6Fe0.4Al zeigt bei 295 K einen hohen negativen Magnetowiderstand von 30 % bei einem Sättigungsfeld von 0.2 Tesla. Durch die Beimischung von Oxiden und Polymeren sind Komposit-Materialien entwickelt worden, die MR0-Effekte von bis zu 88% mit einer verbesserten Feldempfindlichkeit von 0.1 Tesla bei Al2O3 und 0.05 Tesla bei den Oberflächenbeschichtungen zeigen.

Tuning the structural and magnetic properties of Sr 2 FeReO 6 by substituting Fe and Re with valence invariant metals

In the course of this work the effect of metal substitution on the structural and magnetic properties of the double perovskites Sr2MM’O6 (M = Fe, substituted by Cr, Zn and Ga; M’ = Re, substituted by Sb) was explored by means of X-ray diffraction, magnetic measurements, band structure calculations, Mößbauer spectroscopy and conductivity measurements. The focus of this study was the determination of (i) the kind and structural boundary conditions of the magnetic interaction between the M and M’ cations and (ii) the conditions for the principal application of double perovskites as spintronic materials by means of the band model approach. Strong correlations between the electronic, structural and magnetic properties have been found during the study of the double perovskites Sr2Fe1-xMxReO6 (0 < x < 1, M = Zn, Cr). The interplay between van Hove-singularity and Fermi level plays a crucial role for the magnetic properties. Substitution of Fe by Cr in Sr2FeReO6 leads to a non-monotonic behaviour of the saturation magnetization (MS) and an enhancement for substitution levels up to 10 %. The Curie temperatures (TC) monotonically increase from 401 to 616 K. In contrast, Zn substitution leads to a continuous decrease of MS and TC. The diamagnetic dilution of the Fe-sublattice by Zn leads to a transition from an itinerant ferrimagnetic to a localized ferromagnetic material. Thus, Zn substitution inhibits the long-range ferromagnetic interaction within the Fe-sublattice and preserves the long-range ferromagnetic interaction within the Re-sublattice. Superimposed on the electronic effects is the structural influence which can be explained by size effects modelled by the tolerance factor t. In the case of Cr substitution, a tetragonal – cubic transformation for x > 0.4 is observed. For Zn substituted samples the tetragonal distortion linearly increases with increasing Zn content. In order to elucidate the nature of the magnetic interaction between the M and M’ cations, Fe and Re were substituted by the valence invariant main group metals Ga and Sb, respectively. X-ray diffraction reveals Sr2FeRe1-xSbxO6 (0 < x < 0.9) to crystallize without antisite disorder in the tetragonal distorted perovskite structure (space group I4/mmm). The ferrimagnetic behaviour of the parent compound Sr2FeReO6 changes to antiferromagnetic upon Sb substitution as determined by magnetic susceptibility measurements. Samples up to a doping level of 0.3 are ferrimagnetic, while Sb contents higher than 0.6 result in an overall antiferromagnetic behaviour. 57Fe Mößbauer results show a coexistence of ferri- and antiferromagnetic clusters within the same perovskite-type crystal structure in the Sb substitution range 0.3 < x < 0.8, whereas Sr2FeReO6 and Sr2FeRe0.9Sb0.1O6 are “purely” ferrimagnetic and Sr2FeRe0.1Sb0.9O6 contains antiferromagnetically ordered Fe sites only. Consequently, a replacement of the Re atoms by a nonmagnetic main group element such as Sb blocks the double exchange pathways Fe–O–Re(Sb)–O–Fe along the crystallographic axis of the perovskite unit cell and destroys the itinerant magnetism of the parent compound. The structural and magnetic characterization of Sr2Fe1-xGaxReO6 (0 < x < 0.7) exhibit a Ga/Re antisite disorder which is unexpected because the parent compound Sr2FeReO6 shows no Fe/Re antisite disorder. This antisite disorder strongly depends on the Ga content of the sample. Although the X-ray data do not hint at a phase separation, sample inhomogeneities caused by a demixing are observed by a combination of magnetic characterization and Mößbauer spectroscopy. The 57Fe Mößbauer data suggest the formation of two types of clusters, ferrimagnetic Fe- and paramagnetic Ga-based ones. Below 20 % Ga content, Ga statistically dilutes the Fe–O–Re–O–Fe double exchange pathways. Cluster formation begins at x = 0.2, for 0.2 < x < 0.4 the paramagnetic Ga-based clusters do not contain any Fe. Fe containing Ga-based clusters which can be detected by Mößbauer spectroscopy firstly appear for x = 0.4.

Computational studies on the structure and stabilities of magnetic inter-metallic compounds

The purpose of this thesis is to further the understanding of the structural, electronic and magnetic properties of ternary inter-metallic compounds using density functional theory (DFT). Four main problems are addressed. First, a detailed analysis on the ternary Heusler compounds is made. It has long been known that many Heusler compounds ($X_2YZ$; $X$ and $Y$ transition elements, $Z$ main group element) exhibit interesting half-metallic and ferromagnetic properties. In order to understand these, the dependence of magnetic and electronic properties on the structural parameters, the type of exchange-correlation functional and electron-electron correlation was examined. It was found that almost all Co$_2YZ$ Heusler compounds exhibit half-metallic ferromagnetism. It is also observed that $X$ and $Y$ atoms mainly contribute to the total magnetic moment. The magnitude of the total magnetic moment is determined only indirectly by the nature of $Z$ atoms, and shows a trend consistent with Slater-Pauling behaviour in several classes of these compounds. In contrast to experiments, calculations give a non-integer value of the magnetic moment in certain Co$_2$-based Heusler compounds. To explain deviations of the calculated magnetic moment, the LDA+$U$ scheme was applied and it was found that the inclusion of electron-electron correlation beyond the LSDA and GGA is necessary to obtain theoretical description of some Heusler compounds that are half-metallic ferromagnets. The electronic structure and magnetic properties of substitutional series of the quaternary Heusler compound Co$_2$Mn$_{1-x}$Fe$_x$Si were investigated under LDA+$U$. The calculated band structure suggest that the most stable compound in a half-metallic state will occur at an intermediate Fe concentration. These calculated findings are qualitatively confirmed by experimental studies. Second, the effect of antisite disordering in the Co$_2$TiSn system was investigated theoretically as well as experimentally. Preservation of half-metallicity for Co$_2$TiSn was observed with moderate antisite disordering and experimental findings suggest that the Co and Ti antisites disorder amounts to approximately 10~% in the compound. Third, a systematic examination was carried out for band gaps and the nature (covalent or ionic) of bonding in semiconducting 8- and 18-electron or half-metallic ferromagnet half-Heusler compounds. It was found that the most appropriate description of these compounds from the viewpoint of electronic structures is one of a $YZ$ zinc blende lattice stuffed by the $X$ ion. Simple valence rules are obeyed for bonding in the 8- and 18-electron compounds. Fourth, hexagonal analogues of half-Heusler compounds have been searched. Three series of compounds were investigated: GdPdSb, GdAutextit{X} (textit{X} = Mn, Cd and In) and EuNiP. GdPdSb is suggested as a possible half-metallic weak ferromagnet at low temperature. GdAutextit{X} (textit{X} = Mn, Cd and In) and EuNiP were investigated because they exhibit interesting bonding, structural and magnetic properties. The results qualitatively confirm experimental studies on magnetic and structural behaviour in GdPdSb, GdAutextit{X} (textit{X} = Mn, Cd and In) and EuNiP compounds. ~

Structure and properties of intermetallic ternary rare earth compounds

The so called material science is an always growing field in modern research. For the development of new materials not only the experimental characterization but also theoretical calculation of the electronic structure plays an important role. A class of compounds that has attracted a great deal of attention in recent years is known as REME compounds. These compounds are often referred to with RE designating rare earth, actinide or an element from group 1 - 4, M representing a late transition metal from groups 8 - 12, and E belonging to groups 13 - 15. There are more than 2000 compounds with 1:1:1 stoichiometry belonging to this class of compounds and they offer a broad variety of different structure types. Although many REME compounds are know to exist, mainly only structure and magnetism has been determined for these compounds. In particular, in the field of electronic and transport properties relatively few efforts have been made. The main focus in this study is on compounds crystallizing in MgAgAs and LiGaGe structure. Both structures can only be found among 18 valence electron compounds. The f electrons are localized and therefor not count as valence electrons. A special focus here was also on the magnetoresistance effects and spintronic properties found among the REME compounds. An examination of the following compounds was made: GdAuE (E = In, Cd, Mg), GdPdSb, GdNiSb, REAuSn (RE = Gd, Er, Tm) and RENiBi (RE = Pr, Sm, Gd - Tm, Lu). The experimental results were compared with theoretic band structure calculations. The first half metallic ferromagnet with LiGaGe structure (GdPdSb) was found. All semiconducting REME compounds with MgAgAs structure show giant magnetoresistance (GMR) at low temperatures. The GMR is related to a metal-insulator transition, and the value of the GMR depends on the value of the spin-orbit coupling. Inhomogeneous DyNiBi samples show a small positive MR at low temperature that depends on the amount of metallic impurities. At higher fields the samples show a negative GMR. Inhomogeneous nonmagnetic LuNiBi samples show no negative GMR, but a large positive MR of 27.5% at room temperature, which is interesting for application.

Unconventional correlated systems

This work comprises three different types of unconventional correlated systems.rnChapters 3-5 of this work are about the open shell compounds Rb4O6 and Cs4O6. These mixed valent compounds contain oxygen in two different modifications: the closed-shell peroxide anion is nonmagnetic, whereas the hyperoxide anion contains an unpaired electrons in an antibonding pi*-orbital. Due to this electron magnetic ordering is rendered possible. In contrast to theoretical predictions, which suggested half-metallic ferromagnetism for Rb4O6,rndominating antiferromagnetic interactions were found in the experiment. Besidesrna symmetry reduction due to the mixed valency, strong electronic correlations of this highly molecular system determine its properties; it is a magnetically frustrated insulator. The corresponding Cs4O6 was found to show similar properties.rnChapters 6-9 of this work are about intermetallic Heusler superconductors. rnAll of these superconductors were rationally designed using the van Hove scenario as a working recipe. A saddle point in the energy dispersion curve of a solid leads to a van Hove singularity in the density of states. In the Ni-based and Pd-based Heusler superconductors presented in this work this sort of a valence instability occurs at the high-symmetry L point and coincides or nearly coincides with the Fermi level. The compounds escape the high density of states at the Fermi energy through a transition into the correlated superconducting state.rnChapter 10 of this work is about the tetragonally distorted ferrimagnetic DO22 phase of Mn3Ga. This hard-magnetic modification is technologically useful for spin torque transfer applications. The phase exhibits two different crystallographic sites that are occupied by Mn atoms and can thus be written as Mn2MnGa. The competition between the mainly itinerant moments of the Mn atoms at the Wyckoff position 4d and the localized moments of the Mn atoms at the Wyckoff position 2b leads to magnetic correlations. The antiferromagnetic orientation of these moments determines the compound to exhibit a resulting magnetic moment of approximately 1 µB per formula unit in a partially compensated ferrimagnetic configuration.

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.