Ab-initio study of magnons in linear and non-linear regimes

Abstract

Low energy excitations (< 1eV) are of utmost importance for understanding the elec- tronic, magnetic, and thermodynamical properties of any material. The collective excita- tions of spin in this energy regime are magnons, which are quantized spin-waves. Recent developments in the field of spin-dynamics have opened up the world of femtomagnetism, whereby the spin degree-of-freedom is controlled using ultrafast laser pulses. Thus manipu- lating magnons using femtomagnetism holds great promise for future technological devices operating on ultrafast timescales. In order to reach this goal, it is vital that we are able to accurately describe these magnon excitations in order to understand, and ultimately control them. The objective of this thesis is to study magnetic excitations using an ab-initio approach, namely Time Dependent Density Functional Theory (TDDFT). These studies are divided into two sections. In the first part these excitations are studied in the linear regime where a new approximation is derived. In the second part these studies are extended into the real-time, non-equilibrium, regime where the response of magnons to ultrafast laser pulses is investigated. TDDFT encapsulates the electron-electron interactions of the many-body system in the exchange-correlation (XC) functional, which for practical applications must be approxi- mated. Despite the plethora of approximations for the XC energy functional only a few have been used for the XC kernel. Out of these, only the Adiabatic Local Density Approx- imation kernel has been implemented and applied to study the magnetic excitations. The work presented first focuses on deriving the Generalized Gradient Approximation (GGA) kernel by climbing up the Jacob’s ladder of functionals. Then the performance of the GGA kernel is studied by calculating the magnon spectra for ferromagnets and Heuslers. Results show that the GGA kernel generally worsens the spin-excitation spectra by overestimating the magnon energies. However, at the Brillouin Zone (BZ) boundary, suppression of the spin wave intensity is consistent with experimental findings. In the second part of this thesis the power of real-time TDDFT is demonstrated by simulating multiple magnon modes in multi-sublattice alloys, where these modes can be vi- sualized in real space and their response to external pulses can be investigated in real-time. The results in the case of two magnetic sublattice alloy suggests the existence of element specific decoupled magnon modes along with the coupled modes. These decoupled modes in- dicate that the constituent elements in the alloy respond at different timescales. Using these studies three distinct ways of ultrafast laser pulse control of magnon modes is demonstrated: (1) element selective destruction of magnon modes, (2) creation of a transient non-collinear state by destruction of specific magnon modes, (3) renormalization of the optical magnon frequency, where we found a linear dependence between the laser intensity and the decrease of the magnon frequency.

Niederenergetische Anregungen (<1 eV) sind für das Verständnis der elektronischen, mag- netischen und thermodynamischen Eigenschaften eines Materials von größter Bedeutung. Die kollektiven Anregungen des Spins in diesem Energieregime sind quantisierte Spinwellen, die man Magnonen nennt. Jüngste Entwicklungen auf dem Gebiet der Spindynamik er- möglichen den Zugang zur Kontrolle des Femtomagnetismus, wobei der Spinfreiheitsgrad mit ultraschnellen Laserpulsen gesteuert wird. Die Manipulation von Magnonen mithilfe des Femtomagnetismus ist daher vielversprechend für zukünftige technologische Geräte, die in ultraschnellen Zeitskalen arbeiten. Ein genaueres Verständnis der Magnonenanregungen ist dafür allerdings erforderlich. Das Ziel dieser Arbeit ist es daher, magnetische Anregungen mit einem Ab-initio-Ansatz zu untersuchen, nämlich der zeitabhängigen Dichtefunktionaltheorie (TDDFT). Die Arbeit ist in zwei Abschnitte unterteilt. Im ersten Abschnitt werden Magnonenanregungen im linearen Bereich untersucht und eine neue Näherung abgeleitet. Im zweiten Abschnitt werden diese Untersuchungen auf das Echtzeit-Nichtgleichgewichtsregime ausgedehnt, in dem die Reaktion von Magnonen auf ultraschnelle Laserpulse untersucht wird. TDDFT beschreibt die interelektronischen Wechselwirkungen des Vielteilchensystems mithilfe des Austauschkorrelationsfunktionals (XC), das für praktische Anwendungen an- genähert werden muss. Trotz einer Fülle von vorhandenen Näherungen für das XC En- ergiefunktional wurden bisher nur wenige für den XC-Kernel verwendet. Von diesen wurde bisher nur der Kernel des lokalen Dichtefunktionals in der adiabatischen Näherung imple- mentiert und angewendet, um die magnetischen Anregungen zu untersuchen. Der erste Abschnitt dieser Arbeit konzentriert sich daher auf die Ableitung des GGA-Kernels (Gen- eralized Gradient Approximation). Anschließend wird die Genauigkeit des GGA-Kernels untersucht, indem die Magnonenspektren für Ferromagnete und Heuslersche Legierungen berechnet werden. Die Ergebnisse zeigen, dass der GGA-Kernel die Spinanregungsspektren im Allgemeinen verschlechtert, indem die Magnonenenergien überschätzt werden. An der Brillouin Zone (BZ)-Grenze hingegen, stimmt die berechnete Verminderung der Spinwellen- intensität mit den experimentellen Ergebnissen überein. Im zweiten Abschnitt dieser Arbeit wird die Genauigkeit der Echtzeit-TDDFT-Methode demonstriert. Mehrere Magnonenmoden in Legierungen mit mehreren Untergittern wer- den simuliert, wobei die Magnonenmoden im realen Raum graphisch dargestellt und die Reaktion auf externe Laserpulse in Echtzeit untersucht werden. Die Ergebnisse im Fall von zwei magnetischen Untergitterlegierungen legen die Existenz sowohl von elementspezifschen entkoppelten Magnonenmoden als aucg gekoppelter Moden nahe. Sie deuten darauf hin, dass die laserinduzierte Rückreaktion der Elemente in der Legierung auf unterschiedlichen Zeitskalen liegt. Unter Verwendung dieser Studien werden drei verschiedene Möglichkeiten der ultraschnellen Laserpulssteuerung von Magnonenmoden dargelegt: (1) Elementselektive Zerstörung von Magnonenmoden, (2) Erzeugung eines transienten nichtkollinearen Zustands durch Zerstörung spezifischer Magnonenmoden, (3) Renormierung der Magnonenfrequenz im optischen Bereich, wobei wir eine lineare Abhängigkeit zwischen der Laserintensität und der Abnahme der Magnonenfrequenz feststellten.