The electronic structure of the zero-gap two-dimensional graphene has a charge neutrality point exactly at the Fermi level that limits the practical application of this material. There are several ways to modify the Fermi-level-region of graphene, e.g. adsorption of graphene on different substrates or different molecules on its surface. In all cases the so-called dispersion or van der Waals interactions can play a crucial role in the mechanism, which describes the modification of electronic structure of graphene. The adsorption of water on graphene is not very accurately reproduced in the standard density functional theory (DFT) calculations and highly-accurate quantum-chemical treatments are required. A possibility to apply wavefunction-based methods to extended systems is the use of local correlation schemes. The adsorption energies obtained in the present work by means of CCSD(T) are much higher in magnitude than the values calculated with standard DFT functional although they agree that physisorption is observed. The obtained results are compared with the values available in the literature for binding of water on the graphene-like substrates. / Dieser Beitrag ist mit Zustimmung des Rechteinhabers aufgrund einer (DFG-geförderten) Allianz- bzw. Nationallizenz frei zugänglich.

The electronic and geometrical structure of single difluoro-bora-1,3,5,7-tetraphenyl-aza-dipyrromethene (aza-BODIPY) molecules adsorbed on the Au(111) surface is investigated by low temperature scanning tunneling microscopy and spectroscopy in conjunction with ab initio density functional theory simulations of the density of states and of the interaction with the substrate. Our DFT calculations indicate that the aza-BODIPY molecule forms a chemical bond with the Au(111) substrate, with distortion of the molecular geometry and significant charge transfer between the molecule and the substrate. Nevertheless, most likely due to the low corrugation of the Au(111) surface, diffusion of the molecule is observed for applied bias in excess of 1 V. / Dieser Beitrag ist mit Zustimmung des Rechteinhabers aufgrund einer (DFG-geförderten) Allianz- bzw. Nationallizenz frei zugänglich.

The non-adiabatic quantum molecular dynamics (NA-QMD) method couples self-consistently classical nuclear motion with time-dependent density functional theory (TDDFT) in basis expansion for the electron dynamics. It has become a versatile approach to study the dynamics of atoms, molecules and clusters in a wide range of scenarios. This work presents applications of the NA-QMD method to important benchmark systems and its systematic extension to include quantum effects in the nuclear motion. Regarding the first objective, a complete study of the strong-field ionization and dissociation dynamics of nature’s simplest molecule H2+ is performed. By including all electronic and nuclear degrees of freedom and all reaction channels, molecular rotation is shown to play an important role in the ionization process. In addition, strong orientation effects in the energy deposition process of the Buckminster fullerene C60 in short intense laser pulses are surprisingly found in full dimensional calculations. Their consequences on the subsequent nuclear relaxation dynamics shed new light on available experimental data and future experiments are proposed to confirm the detailed predictions. Regarding the second objective, the NA-QMD formalism is basically extended to take electron-nuclear correlations into account. This extension is achieved by means of a trajectory surface hopping scheme in the adiabatic Kohn-Sham framework. First studied examples from collision physics and photochemistry illustrate the relevance and importance of quantum effects in the nuclear dynamics.

The interplay between different self-organized electronically ordered states and their relation to unconventional electronic properties like superconductivity constitutes one of the most exciting challenges of modern condensed matter physics. In the present thesis this issue is thoroughly investigated for the prototypical layered material 1T-TaS2 both experimentally and theoretically. At first the static charge density wave order in 1T-TaS2 is investigated as a function of pressure and temperature by means of X-ray diffraction. These data indeed reveal that the superconductivity in this material coexists with an inhomogeneous charge density wave on a macroscopic scale in real space. This result is fundamentally different from a previously proposed separation of superconducting and insulating regions in real space. Furthermore, the X-ray diffraction data uncover the important role of interlayer correlations in 1T-TaS2. Based on the detailed insights into the charge density wave structure obtained by the X-ray diffraction experiments, density functional theory models are deduced in order to describe the electronic structure of 1T-TaS2 in the second part of this thesis. As opposed to most previous studies, these calculations take the three-dimensional character of the charge density wave into account. Indeed the electronic structure calculations uncover complex orbital textures, which are interwoven with the charge density wave order and cause dramatic differences in the electronic structure depending on the alignment of the orbitals between neighboring layers. Furthermore, it is demonstrated that these orbital-mediated effects provide a route to drive semiconductor-to-metal transitions with technologically pertinent gaps and on ultrafast timescales. These results are particularly relevant for the ongoing development of novel, miniaturized and ultrafast devices based on layered transition metal dichalcogenides. The discovery of orbital textures also helps to explain a number of long-standing puzzles concerning the electronic self-organization in 1T-TaS2 : the ultrafast response to optical excitations, the high sensitivity to pressure as well as a mysterious commensurate phase that is commonly thought to be a special phase a so-called “Mott phase” and that is not found in any other isostructural modification.

A preconditioner to improve the convergence properties of Krylov subspace solvers is derived and analyzed in this work. This method is adapted to linear systems arising from a finite-element discretization of a phase-field crystal equation.

In this thesis, a series of computational studies based on density functional theory (DFT) and density functional tight-binding (DFTB) is presented to deeply understand experimental results on the synthesis of endohedral fullerenes and graphene/iron hybrids at atomic level. In the first part, a simple and efficient model is proposed to evaluate the strain experienced by clusters encapsulated in endohedral metallofullerenes (EMFs). Calculations for the sole cluster, either in the neutral or the charged state, cannot be used for this goal. However, when the effect of the carbon cage is mimicked by small organic π-systems (such as pentalene and sumanene), the cluster has sufficient freedom to adopt the optimal configuration, and therefore the energetic characteristics of the EMF-induced distortion of the cluster can be evaluated. Both nitride and sulfide clusters were found to be rather flexible. Hence, they can be encapsulated in carbon cages of different size and shape. For carbide M2C2 cluster the situation is more complex. The optimized cluster can adopt either butterfly or linear shapes, and these configurations have substantially different metal-metal distance. Whereas for Sc2C2 both structures are isoenergetic, linear form of the Y2C2 cluster is substantially less stable than the butterfly-shaped configuration. These results show that phenomenon of the “nanoscale fullerene compression” once proposed by Zhang et al. (J. AM. CHEM. SOC. (2012),134(20)) should be “nanoscale fullerene stretching”. Finally, the results also reveal that both Ti2S and Ti2C2 cluster are strained in corresponding EMF molecules, but the origin of the strain is opposite: C78-D3h(5) cage imposes too long Ti···Ti distance for the sulfide cluster and too short distance for the carbide cluster. In the second part of the thesis, possible fullerene geometries and electronic structures have been explored theoretically for the species detected in mass spectra of the Sc-EMF extract synthesized using CH4 as a reactive gas. Two most promising candidates, namely Sc4C@C80-Ih(7) and Sc4C3@C80-Ih(7), have been identified and further studied at the DFT level. For Sc4C@C80, the tetrahedral Sc4 cluster with the central μ4-C atom was found to be 10 kJ/mol more stable than the square cluster. For Sc4C3@C80, the calculation showed that the most stable is the Sc4C3 cluster in which the triangular C3 moiety is η3- and η2-coordinated to Sc atoms. Whereas Sc4C@C80 has rather small HOMO-LUMO gap and low ionization potential, the HOMO-LUMO gap of Sc4C3@C80 is substantially higher and exceeds that of Sc4C2@C80. In the third part, computational studies of structures and reactivity are described for a new type of EMFs with a heptagon that has been produced in the arc-discharge synthesis. DFT computations predict that LaSc2N@Cs(hept)-C80 is more stable than LaSc2N@D5h-C80, so the former should be synthesized in much higher yield than observed. This disagreement may be ascribed to the kinetic factors rather than thermodynamic stability. Because of prospective applications of this EMFs by introducing functional groups, the influence of the heptagon on the chemical properties have been further evaluated. Thermodynamically and kinetically preferred reaction sites are studied computationally for Prato and Bingel-Hirsch cycloaddition reactions. In both types of reactions the heptagon is not affected, and chemical reactivity is determined by the adjacent pentalene units. Thermodynamically controlled Prato addition is predicted to proceed regioselectively across the pentagon/pentagon edges, whereas the most reactive sites in kinetically-controlled Bingel-Hirsch reaction are the carbon atoms next to the pentagon/pentagon edge. Fourth, although various EMFs have been successfully synthesized and characterized, the formation mechanism is still not known in details, and hence control of the synthesis products is rather poor. Therefore, EMF self-assembly process in Sc/carbon vapor in the presence and absence of cooling gas (helium) and reactive gas (NH3 and CH4) is systematically investigated using quantum chemical molecular dynamics (QM/MD) simulations based on the DFTB potentials. The cooling gas effect is that the presence of He atoms accelerates formation of pentagons and hexagons and reduces the size of formed carbon cages in comparison to the analogous He-free simulations. As a result, the Sc/C/He system yields a large number of successful trajectories (i.e. leading to the Sc-EMFs) with more realistic cage-size distribution than the Sc/C system. Encapsulation of Sc atoms within the carbon cage was found to proceed via two parallel mechanisms. The main mechanism involves nucleation of the several hexagons and pentagons with Sc atoms already at the early stages of the carbon vapor condensation. In such proto-cages, both Sc–C σ-bonds and coordination bonds between Sc atoms and the π-system of the carbon network are present. Sc atoms are thus rather labile and can move along the carbon network, but the overall bonding is sufficiently strong to prevent dissociation even at high temperatures. Further growth of the carbon cage results in encapsulation of one or two Sc atoms within the forming fullerene. Another encapsulation mechanism is observed in rare cases. In this process, the closed cage is formed with Sc being a part of the carbon network, i.e. being bonded by three or four Sc–C σ-bonds. However, such intermediates are found to be unstable, and transform into the endohedral fullerenes within few picoseconds of annealing. In perfect agreement with experimental studies, extension of the simulation to Fe and Ti showed that Fe-EMFs are not formed at all, whereas Ti is prone to form Ti-EMFs with small cage sizes, including Ti@C28-Td and Ti@C30-C2v(3). The role of “reactive gas” in the EMF synthesis is revealed in dedicated simulations of the fullerene formation in the presence of several molecules of CH4 or NH3. When concentration of reactive gas is high, carbon vapor tends to form graphene flakes or other carbon species terminated by hydrogen atoms, whereas the yield of empty fullerenes is very low. Conversely, with additional metal atoms (Sc) and the same number of NH3 molecules, the yield of fullerenes constantly increase from 5 to 65% which is ascribed to the catalytic activity of metal atoms in the nucleation of carbon cages already at early stage. Moreover, due to the presence of hydrogen atoms from the reactive gas, the carbon cage formation requires much longer time, which provides sufficient reaction time to encapsulate 3 or 4 Sc atoms within one cage. It explains preferential formation of clusterfullerenes in experiments with reactive gas. At the same time, monometallofullerenes and dimetallofullerenes are the main products in absence of reactive gas. We also provide possible growth mechanisms of carbide and cyano-clusterfullerenes in details to elucidate how the intracluster goes into the cage. A possible growth mechanism of nitride clusterfullerenes has been proposed based on DFT results. In the last part, a free-standing crystalline single-atom thick layer of Fe has been studied theoretically. By investigating the energy difference, ΔE, between a suspended Fe monolayer and a nanoparticle using the equivalent number of Fe atoms, one can estimate that the largest stable membrane should be ca. 12 atoms wide or 3 × 3 nm2 which is in excellent agreement with the experimental observation. Otherwise, the possibility of C, O, N atoms embedded into the Fe membrane can been fully excluded by DFTB and DFT simulations, which agrees with electron energy loss spectroscopy (EELS) measurement. A significantly enhanced magnetic moment for single atom thick Fe membranes (3.08 μB) is predicted by DFT as compared to the bulk BCC Fe (2.1 μB), which originates from the 2D nature of the Fe membrane since the dz2 orbital is out-of-plane while the dxy orbital is in-plane.

Die vorliegende Promotionsarbeit “Density functional study of the electronic and magnetic properties of selected transition metal complexes” beschäftigt sich mit dem Zusammenhang zwischen strukturellen Merkmalen sowie elektronischen und magnetischen Eigenschaften von Einzelmolekül-Magneten. Im Wesentlichen konnte dabei gezeigt werden, dass die magnetischen Eigenschaften sowohl von strukturellen Merkmalen als auch von den elektronischen Eigenschaften bestimmt werden. Des Weiteren ergab sich, dass verschiedene Kenngrößen der magnetischen Eigenschaften (im speziellen der magnetische Grundzustand S sowie die magnetische Anisotropie D) miteinander korreliert sind. Dies ist im Besonderen für eine mögliche Anwendung von Einzelmolekül-Magneten im Bereich der Datenspeicherung von Bedeutung.

In der vorliegenden Arbeit wurde ausgedehntes Si_1−x Ge_x für unterschiedliche Zusammensetzungen 0 ≤ x ≤ 1 untersucht. Die Untersuchungen basierten auf der DFT, wobei das Programm QuantumATK 18.06 zum Einsatz kam. Für die Korrektur der Bandlücke wurden empirische Pseudopotential Projektor Shifts verwendet [34]. Für jede untersuchte Zusammensetzung wurden 500 zufällig generierte Konfigurationen der 64-atomigen Superzelle berechnet und statistisch ausgewertet. Nach der Optimierung der Struktur erfolgte die Auswertung der Bandlücke indem über äquivalente Pfade in der Brillouinzone gemittelt wurde. Zusätzlich wurden nach dieser Art auch kleine Anteile an C untersucht. Die Ergebnisse der Berechnungen zeigen für die Bildungsenergie der Mischstrukturen positive Werte mit einem Maximum bei mittleren Zusammensetzungen. Zur Stabilitätsuntersuchung der Legierungen wurde die Gibbs-Energie berechnet. Es ergeben sich negative Werte, was die Stabilität von SiGe bestätigt. Die berechnete Gitterkonstante der relaxierten Strukturen zeigt eine leichte Überschätzung der experimentellen Werte. Die ermittelten Bandlücken reproduzieren den Übergang von Si-artigen zu Ge-artigen Bandlücken bei x = 0.85. Die Werte der Bandlücke zeigen eine gute Übereinstimmung mit dem Experiment. Aus den statistischen Untersuchungen wird deutlich, dass sowohl Bildungsenergie als auch Bandlücke Variationen von 10 % und mehr aufweisen. Es zeigt sich dadurch ein nicht zu vernachlässigender Unterschied zwischen verschie denen Konfigurationen der Superzelle, die alle eine Legierung mit gleicher Zusammensetzung beschreiben. Wird in die Strukturen Kohlenstoff eingebracht, so vergrößern sich die Variationen mit steigendem C-Anteil. Für die betrachteten kleine C-Anteile zeigt sich eine Erhöhung der Bildungsenergie und einer Verkleinerung der Gitterkonstante und der Bandlücke. Es wird deutlich, dass bereits wenig C einen Einfluss auf die wichtigen Eigenschaften der Legierung hat und für genaue Simulationen berücksichtigt werden muss. Wie die Ergebnisse zeigen, spielt die spezielle Konfiguration von Strukturen im nm-Bereich eine wichtige Rolle. Aus diesem Grund wurde im zweiten Teil der Arbeit ein Ge-Profil nachgebildet, wie es in der Basis von HBTs vorkommt. Die Ergebnisse zeigen eine Verkleinerung der Bandlücke im SiGe-Bereich, welche im Wesentlichen durch zusätzliche Valenzzustände hervorgerufen wird. Diese Zustände sind in die z-Richtung lokalisiert. Die Leitungsbandkante bleibt von der SiGe-Region nahezu unbeeinflusst. Die Vergrößerung der SiGe-Region verkleinert die Bandlücke.:Abkürzungsverzeichnis - 5 1 Motivation - 6 2 Theoretische Grundlagen der Dichtefunktionaltheorie - 8 2.1 Quantenmechanische Vielteilchensysteme - 8 2.2 Hohenberg-Kohn-Theoreme - 9 2.3 Austausch-Korrelations-Funktional und Kohn-Sham-Gleichung - 10 3 Siliziumgermanium - 12 3.1 Kristallstruktur und Gitterkonstante - 12 3.2 Bandstruktur - 13 3.2.1 Bandstruktur von Si und Ge - 13 3.2.2 Bandlücke von SiGe - 14 3.2.3 Bandlücke von SiGe:C - 15 4 Modellierung und Methoden - 16 4.1 Modellzellen - 16 4.1.1 8-atomige konventionelle Einheitszelle - 16 4.1.2 64-atomige Superzelle - 17 4.2 Bildungsenergie und Stabilität von Legierungen - 20 4.2.1 Gibbs-Energie - 21 4.3 Faltung der Bandstruktur - 22 4.4 Korrektur und Ermittlung der Bandlücke - 24 4.4.1 Korrektur der Bandlücke - 24 4.4.2 Bestimmung der Bandlücke von ungeordneten Legierungen - 26 4.5 Berechnungsverfahren der Kristallstrukturen - 28 5 Ergebnisse und Auswertung - 29 5.1 Gitterkonstante - 29 5.2 Bildungsenergie und Änderung der Gibbs-Energie - 32 5.3 Bandlücke - 36 5.3.1 Leitungsbandminimum - 38 5.3.2 Bildungsenergie - 40 5.4 Bandstruktur - 42 6 Anwendung für die Basis von HBTs - 44 6.1 Modellierung - 45 6.2 Ergebnisse - 46 7 Zusammenfassung und Ausblick Literatur - 49 Danksagung - 53 Selbstständigkeitserklärung - 54

Epitaxial growth is an important technique for the fabrication of film structures with good crystalline quality, e.g., monoatomic overlayers, multilayers, compound materials, and ordered alloys. Such epitaxially grown films are technologically important materials with, e.g., adjustable electronic, magnetic, and optical properties. In case of coherent or pseudomorphic epitaxy, the overlayer adapts the in-plane lattice parameters of the substrate, i.e., the overlayer is strained to match the lattice parameters parallel to the substrate surface (in-plane directions). Simultaneously, a relaxation of the film dimension perpendicular to the substrate-film interface occurs (out-of-plane direction). Thus, coherent epitaxy provides a method to put phases under strain, and it can stabilise a metastable state of the film material, if the substrate lattice matches this metastable structure. Bulk-like properties in thick overlayers, which adopt the body-centred tetragonal (BCT) crystal structure and which grow coherently on a suitable substrate with quadratic surface symmetry, are modelled by the epitaxial Bain path (EBP) in this thesis. The knowledge of the EBP allows to study properties of the overlayer as function of the substrate lattice parameter. In particular, strain effects on the film material, magnetic order in the overlayer, and the existence of possible metastable states are investigated by means of density functional theory (DFT) in the local spin density approximation (LSDA), and in the singular case of uranium, employing the generalised gradient approximation (GGA). Note that a symmetry property of the BCT structure states, that it is identical to the body-centred cubic (BCC) structure or the face-centred cubic (FCC) structure for definite ratios of the tetragonal lattice parameters. Our definition of the EBP has two, previously not considered consequences for EBPs in general: an EBP can be discontinuous, and the high symmetry cubic structures (FCC and BCC) need not be points on the EBP. Both cases occurred for several elements considered in this thesis. If, however, a cubic structure is a point on the EBP, then a symmetry property guarantees that the total energy along the EBP, E(a), is stationary at this cubic structure. We computed the EBPs of all transition metals (TMs), the post TMs Zn, Cd, and Hg, the alkaline earth metals Ca, Sr, and Ba, the lanthanides La and Lu, and the actinide U (35 elements were treated in total). For each element but Zr, Hg, and U, there are exactly two structures whose energies are minima on the EBP, and which exhibit neither in-plane nor out-of-plane stresses; for Zr, Hg, and U there are three minima each. All other states on the EBP exhibit in-plane stresses because they are a strained form of the stress-free structures. The possibility of metastability of these particular, stress-free structures, i.e., stabilisation of these structures without bonding to the substrate, was investigated by stability conditions based on linear elasticity theory (except for U). We predict that ten FCC structures and three BCT structures not known from the respective phase diagrams may be metastable. We studied the properties of ferromagnetic (FM) states on the EBP for the elements Fe, Co, and Ni, and moreover predict, that Mn, Ru, Os, and U order ferromagnetically for certain states of the EBP. The latter three elements are paramagnetic in their ground states. The onset of ferromagnetism in Os and U is not accompanied by a simultaneously fulfilled Stoner criterion. According to our results, antiferromagnetic order (with moment sequences up-down or up-up-down-down on successive (001) planes) is never more stable than FM order on any EBP for any element investigated. On the basis of our comprehensive results for all TMs, we analysed trends across each of the three TM series and similarities among the three series. We demonstrate, that the type of the EBP (a classification of extrema of E(a) by symmetry into types) follows a characteristic trend across each of the three TM series. We discuss exceptions (Mn, Fe, and Zr) to this trend. Another trend, identical for the three series, is found for the BCT­-FCC structural energy difference as function of the d-band filling (evaluated for BCT structures that define extrema of E(a)), which follows a similar trend as the well studied BCC­-FCC structural energy difference. Clear similarities among the three periods of elements are also reflected in the bulk moduli and in the elastic constants of the cubic or tetragonal structures, that define the global and local minima of E(a). The mentioned similarities suggest, that many properties which are associated with the EBPs of TMs, can be attributed to the occupation of the d-band, which is the most dominant feature of the electronic structure of TMs.

The aim of this thesis was to find anomalies of elastic properties induced by topological changes of the Fermi surface. The latter are called "Lifshitz transitions". Lifshitz transitions are an interesting subject to study because a topological change of the Fermi surface results in a van Hove singularity of the density of states at the Fermi energy, which again induces an anomaly in the free energy and therefore yield anomalies of observable physical quantities. In all cases the question arose, if the corresponding van Hove singularities are large enough to cause anomalies in the elastic properties, which are measurable by nowadays experimental techniques and computable within the accuracy reachable in nowadays computer calculations. The calculations have been done with the Full-Potential nonorthogonal Local-Orbital minimum-basis band-structure code FPLO. To shift the van Hove singularities through the Fermi energy we used hydrostatic pressure, which is mimicked in the computations by decreasing the volume of the unit cell. The materials under consideration had been YCo5 and LaCo5 as examples for magnetic compounds and the element Osmium as an example for a non-magnetic material. All these materials exhibit hexagonal symmetry. In the case of YCo5 our calculations yield a first order Lifshitz transition. Here, an extraordinarily large peak in the spin-up part of the DOS, which is caused by a nearly dispersionless band in the hexagonal plane, crosses the Fermi level under a pressure of about 21 GPa. Thus, the spin-up 3d states become partly depopulated, which results in a drop of the total magnetic moment of 35%. Therefore the transition can be regarded as a transition from strong to weak ferromagnetism. Further, the transition results in a volume collapse of 1.4%. Though the volume collapse is isomorphic, it exhibits the following anisotropy: while the lattice constant in the hexagonal plane is almost smoothly contracting with increasing pressure, the lattice constant in c-direction collapses at the transition-pressure. This volume collapse has been verified in experiment. Analogous calculations have been performed for the compound LaCo5, which is isoelectronic to YCo5. Here as well we predict a first order Lifshitz transition, taking place at a pressure of about 23 GPa. The mechanism of the transition is the same than in YCo5. Again we find a volume collapse under pressure together with a decrease of the magnetic moment. The relative volume change amounts to 1.3%. Like in YCo5, the unit cell dimensions in the hexagonal plane are decreasing almost smoothly with pressure while in c-direction the lattice constant collapses at the transition-pressure. For LaCo5 there are no such experiments done so far to the best of our knowledge. For Osmium we found, that LDA reproduces the ground state volume very well. Furthermore, we could detect three Lifshitz transitions taking place at very high pressures of about 72 GPa, 81 GPa, and 122 GPa. At first, a hole ellipsoid appears at the Gamma-point (V=24.6Å^3, P=72 GPa), then a neck is created at the symmetry-line LH (V=24.2Å^3, P=81 GPa), and finally a hole ellipsoid appears at the L-point (V=23.2 Å^3, P=122 GPa). Due to a degeneracy in the band structure, the hole ellipsoid at the L-point appears at the same pressure when the necks, situated at the symmetry-lines LH merge at L. The corresponding van Hove singularities in the DOS are very tiny and thus no anomalies in the elastic properties could be detected. Furthermore, we showed that the kink in c/a at 25 GPa and at 27 GPa found by Occelli et al. [Occelli et al., Phys. Rev. Lett. 93, 095502 (2004)] and Ma et al. [Ma et al., Phys. Rev. B 72, 174103 (2005)], respectively, is not statistically significant and that (c/a)(P) can be fitted equally well by a smooth function as by piece-wise linear functions as proposed in these references. / Das Ziel dieser Arbeit war es, Anomalien in den elastischen Eigenschaften zu finden, die durch topologische Änderungen der Fermifläche - genannt "Lifschitz Übergänge" - hervorgerufen werden. Lifschitz Übergänge sind ein interessantes Forschungsgebiet, denn die topologische Änderung der Fermifläche führt zu einer van Hove Singularität in der Zustandsdichte, die an der Fermienergie liegt und eine Anomalie in der freien Energie hervorruft und deswegen zu Anomalien in beobachtbaren physikalischen Größen führt. In allen Fällen kam die Frage auf, ob die entsprechenden van Hove Singularitäten groß genug sind, um mit heutigen Methoden meßbare und berechenbare Anomalien in den elastischen Eigenschaften zu verursachen. Die Daten wurden mit dem Computerprogramm FPLO (Full-Potential nonorthogonal Local-Orbital minimum-basis band-structure scheme) berechnet. Um die van Hove Singularitäten durch die Fermienergie zu schieben, verkleinerten wir das Volumen der Einheitszelle, um hydrostatischen Druck zu simulieren. Als zu untersuchende Stoffe wurden YCo5 und LaCo5 als Beispiele für magnetische Verbindungen gewählt und Osmium als Beispiel für ein nicht magnetisches Element. Im Falle von YCo5 fanden wir einen Lifschitz Übergang erster Ordnung. Hier springt ein besonders großer Peak im Spin-auf Teil der Zustandsdichte unter einem Druck von ca. 21 GPa über die Fermienergie. Dadurch werden die Spin-auf 3d Zustände teilweise unbesetzt und das magnetische Moment verringert sich um 35%. Deswegen kann man den Übergang als einen Übergang von starkem Ferromagnetismus zu schwachem Ferromagnetismus bezeichnen. Das Volumen verkleinert sich hierbei um 1.4%. Obwohl dieser Volumenkollaps isomorph ist, zeigt er folgende Anisotropie: während die Gitterkonstante in der hexagonalen Ebene mit zunehmendem Druck mehr oder weniger glatt kontrahiert, kollabiert am Übergangsdruck die Gitterkonstante in c-Richtung. Dieser Volumenkollaps wurde vom Experiment verifiziert. Analoge Rechnungen wurden für die Verbindung LaCo5, die isoelektronisch zu YCo5 ist, durchgeführt. Auch hier sagen wir einen Lifschitz Übergang erster Ordnung voraus, der bei einem Druck von ca. 23 GPa stattfinden wird. Der Mechanismus dieses Übergangs ist der selbe wie in YCo5. Wiederum finden wir einen Volumenkollaps unter Druck zusammen mit einer Verringerung des magnetischen Moments. Die relative Volumenänderung beträgt hier 1.3%. Wie in YCo5 verläuft hier die Kontraktion der Gitterkonstante in der hexagonalen Ebene mehr oder weniger glatt, während die Gitterkonstante in c-Richtung am Übergang kollabiert. Für LaCo5 existieren unseres Wissens hierzu noch keine Experimente. Im Falle von Osmium fanden wir drei Lifschitz Übergänge bei sehr hohen Drücken von ca. 72 GPa, 81 GPa, und 122 GPa. Zuerst bildet sich ein Lochellipsoid am Gamma-Punkt (V=24.6Å^3, P=72 GPa), dann bildet sich ein Hals an der Symmetrielinie LH (V=24.2Å^3, P=81 GPa), und zum Schluß erscheint ein Lochellipsoid am L-Punkt (V=23.2 Å^3, P=122 GPa). Auf Grund einer Entartung in der Bandstruktur taucht das Lochellipsoid am L-Punkt an dem Druck auf, an dem sich auch die Hälse auf der Symmetrielinie LH bei L verbinden. Die entsprechenden van Hove Singularitäten in der Zustandsdichte sind jedoch extrem klein und deswegen können keine Anomalien in den elastischen Eigenschaften detektiert werden. Desweiteren zeigten wir, daß der Knick in c/a, den Occelli et al. [Occelli et al., Phys. Rev. Lett. 93, 095502 (2004)] bei 25 GPa und Ma et al. [Ma et al., Phys. Rev. B 72, 174103 (2005)] bei 27 GPa fanden, statistisch nicht relevant ist und daß (c/a)(P) genauso gut von einer glatten Funktion gefittet wird als von stückweise linearen Funktionen.