Researcher profile

Stefan Blügel

Stefan Blügel contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
41works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2024arXiv

Machine learning inspired models for Hall effects in non-collinear magnets

The anomalous Hall effect has been front and center in solid state research and material science for over a century now, and the complex transport phenomena in nontrivial magnetic textures have gained an increasing amount of attention, both in theoretical and experimental studies. However, a clear path forward to capturing the influence of magnetization dynamics on anomalous Hall effect even in smallest frustrated magnets or spatially extended magnetic textures is still intensively sought after. In this work, we present an expansion of the anomalous Hall tensor into symmetrically invariant objects, encoding the magnetic configuration up to arbitrary power of spin. We show that these symmetric invariants can be utilized in conjunction with advanced regularization techniques in order to build models for the electric transport in magnetic textures which are, on one hand, complete with respect to the point group symmetry of the underlying lattice, and on the other hand, depend on a minimal number of order parameters only. Here, using a four-band tight-binding model on a honeycomb lattice, we demonstrate that the developed method can be used to address the importance and properties of higher-order contributions to transverse transport. The efficiency and breadth enabled by this method provides an ideal systematic approach to tackle the inherent complexity of response properties of noncollinear magnets, paving the way to the exploration of electric transport in intrinsically frustrated magnets as well as large-scale magnetic textures.

preprint2023arXiv

Construction of Wannier functions from the spectral moments of correlated electron systems

When the first four spectral moments are considered, spectral features missing in standard Kohn-Sham (KS) density-functional theory (DFT), such as upper and lower Hubbard bands, as well as spectral satellite peaks, can be described, and the bandwidths can be corrected. Therefore, we have devised a \textit{moment-functional based spectral density functional theory} (MFbSDFT) recently. However, many computational tools in theoretical solid state physics, such as the construction of maximally localized Wannier functions (MLWFs), have been developed for KS-DFT and require modifications if they are supposed to be used in MFbSDFT. Here, we show how generalized Wannier functions may be constructed from the first four spectral moment matrices. We call these functions \textit{maximally localized spectral moment Wannier functions} (MLSMWFs). We demonstrate how MLSMWFs may be used to compute the anomalous Hall effect (AHE) in fcc Ni by Wannier interpolation. More generally, MLSMWFs may be computed from the first $2P$ moments ($P=1,2,3,\dots$). Using more than 4 moments opens the perspective of reproducing all spectral features accurately in MFbSDFT.

preprint2022arXiv

Bulk and surface electronic structure of Bi$_4$Te$_3$ from $GW$ calculations and photoemission experiments

We present a combined theoretical and experimental study of the electronic structure of stoichiometric Bi$_4$Te$_3$, a natural superlattice of alternating Bi$_2$Te$_3$ quintuple layers and Bi bilayers. In contrast to the related semiconducting compounds Bi$_2$Te$_3$ and Bi$_1$Te$_1$, density functional theory predicts Bi$_4$Te$_3$ to be a semimetal. In this work, we compute the quasiparticle electronic structure of Bi$_4$Te$_3$ in the framework of the $GW$ approximation within many-body perturbation theory. The quasiparticle corrections are found to modify the dispersion of the valence and conduction bands in the vicinity of the Fermi energy, leading to the opening of a small indirect band gap. Based on the analysis of the eigenstates, Bi$_4$Te$_3$ is classified as a dual topological insulator with bulk topological invariants $\mathbb{Z}_2$ (1;111) and magnetic mirror Chern number $n_M=1$. The bulk $GW$ results are used to build a Wannier-functions based tight-binding Hamiltonian that is further applied to study the electronic properties of the (111) surface. The comparison with our angle-resolved photoemission measurements shows excellent agreement between the computed and measured surface states and indicates the dual topological nature of Bi$_4$Te$_3$.

preprint2022arXiv

Chiral standing spin waves in skyrmion lattice

This work studies the resonance excitations of the three-dimensional skyrmions lattice in the finite thickness plate of an isotropic chiral magnet using spin dynamics simulations. We found that the absorption spectra and resonance modes differ from those predicted by the two-dimensional model and the model of the unconfined bulk crystal. The features observed on the spectra can be explained by the formation of chiral standing spin waves, which, contrary to conventional standing spin waves, are characterized by the helical profile of dynamic magnetization of fixed chirality defined by the Dzyaloshinskii-Moriya interaction. In this case, the dynamic susceptibility becomes a function of the plate thickness, which gives rise to an interesting effect that manifests itself in periodical fading of the intensity of corresponding modes and makes excitation of these modes impossible at specific thicknesses.

preprint2022arXiv

Construction of the spectral function from non-commuting spectral moment matrices

The LDA+U method is widely used to study the properties of realistic solids with strong electron correlations. One of its main shortcomings is that it does not provide direct access to the temperature dependence of material properties such as the Dzyaloshinskii-Moriya interaction, the anomalous Hall conductivity, and the spin-orbit torque. While the method of spectral moments allows us in principle to compute these quantities directly at finite temperatures, the standard two-pole approximation can be applied only to Hamiltonians that are effectively of single-band type. We do a first step to explore if the method of spectral moments may replace the LDA+U method in first-principles calculations of correlated solids with many bands in cases where the direct assessment of the temperature dependence of equilibrium and response functions is desired: The spectral moments of many-band Hamiltonians of correlated electrons do not commute and therefore they do not possess a system of common eigenvectors. We show that nevertheless the spectral function may be constructed from the spectral moments by solving a system of coupled non-linear equations. Additionally, we show how to compute the anomalous Hall conductivity of correlated electrons from this spectral function. We demonstrate the method for the Hubbard-Rashba model, where the standard two-pole approximation cannot be applied, because spin-orbit interaction (SOI)couples the spin-up and the spin-down bands. In the quest for new quantum states that arise from the combination of SOI and correlation effects, the Hartree-Fock approximation is frequently used to obtain a first approximation for the phase diagram. We propose that using the many-band generalization of the selfconsistent moment method instead of Hartree-Fock in such exploratory model calculations may improve the accuracy significantly, while keeping the computational burden low.

preprint2022arXiv

Diversity of states in a chiral magnet nanocylinder

The diversity of three-dimensional magnetic states in a FeGe nanocylinder is studied using micromagnetic simulations and off-axis electron holography in the transmission electron microscope. In particular, we report the observation of a dipole string -- a spin texture composed of two coupled Bloch points -- which becomes stable under geometrical confinement. Quantitative agreement is obtained between experimental and theoretical phase shift images by taking into account the presence of a damaged layer on the surface of the nanocylinder. The theoretical model is based on the assumption that the damaged surface layer, which results from focused ion beam milling during sample preparation, has similar magnetic properties to those of an amorphous FeGe alloy. The results highlight the importance of considering the magnetic properties of the surface layers of such nanoscale samples, which influence their magnetic states.

preprint2022arXiv

Driving spin chirality by electron dynamics in laser-excited antiferromagnets

Optical generation of complex spin textures is one of the most exciting challenges of modern spintronics. Here, we uncover a distinct physical mechanism for imprinting spin chirality into collinear magnets with short laser pulses. By simultaneously treating the laser-ignited evolution of electronic structure and magnetic order, we show that their intertwined dynamics can result in an emergence of quasi-stable chiral states. We find that laser-driven chirality does not require any auxiliary external fields or intrinsic spin-orbit interaction to exist, and it can survive on the time scale of nanoseconds even in the presence of thermal fluctuations, which makes the uncovered mechanism relevant for understanding various optical experiments on magnetic materials. Our findings open a new perspective at the interaction of complex chiral magnetism with light.

preprint2022arXiv

Engineering spin-orbit effects and Berry curvature by deposition of Eu on WSe$_{2}$

Motivated by recent progress in 2D spintronics, we present Eu deposited on a 1H-WSe$_2$ as a promising platform for engineering spin-orbit effects and Berry curvature. By first-principles calculations based on density functional theory, we show that Eu/WSe$_2$ exhibits intriguing properties such as high magnetic anisotropy, valley-dependent polarization of spin and orbital angular momenta, and their Rashba textures. These originate from magnetic and spin-orbit proximity effects at the interface and the interplay between localized $4f$ magnetic moments of Eu and mobile charge carriers of WSe$_2$. We find a pronounced anomalous Hall effect in the proposed system. Thus, we promote $4f$ rare-earth metals deposited on top of a transition-metal dichalcogenides as a promising platform for 2D spintronics.

preprint2022arXiv

Fermi-surface origin of skyrmion lattices in centrosymmetric rare-earth intermetallics

We show from first-principles that barrel-shaped structures within the Fermi surface of the centrosymmetric intermetallic compounds GdRu$_2$Si$_2$ and Gd$_2$PdSi$_3$ give rise to Fermi surface nesting, which determines the strength and sign of quasi-two-dimensional Ruderman-Kittel-Kasuya-Yosida pairwise exchange interactions between the Gd moments. This is the principal mechanism leading to their helical single-$q$ spin-spiral ground states, providing transition temperatures and magnetic periods in good agreement with experiment. Using atomistic spin-dynamic simulations, we draw a direct line between the subtleties of the three-dimensional Fermi surface topology and the stabilization of a square skyrmion lattice in GdRu$_2$Si$_2$ at applied magnetic fields as observed in experiment.

preprint2022arXiv

Is there a polaron signature in angle-resolved photoemission of CsPbBr3?

The formation of large polarons has been proposed as reason for the high defect tolerance, low mobility, low charge carrier trapping and low nonradiative recombination rates of lead halide perovskites. Recently, direct evidence for large-polaron formation has been reported from a 50% effective mass enhancement in angle-resolved photoemission of CsPbBr3 over theory for the orthorhombic structure. We present in-depth band dispersion measurements of CsPbBr3 and GW calculations which lead to almost identical effective masses at the valence band maximum of 0.203+/-0.016 m0 in experiment and 0.226 m0 in orthorhombic theory. We argue that the effective mass can be explained solely on the basis of electron-electron correlation and large-polaron formation cannot be concluded from photoemission data.

preprint2022arXiv

Link between Weyl-fermion chirality and spin texture

Topological semimetals have recently attracted great attention due to prospective applications governed by their peculiar Fermi surfaces. Weyl semimetals host chiral fermions that manifest as pairs of non-degenerate massless Weyl points in their electronic structure, giving rise to novel macroscopic quantum phenomena such as the chiral anomaly, an unusual magnetoresistance, and various kinds of Hall effects These properties enable the engineering of non-local electric transport devices, magnetic sensors and memories, and spintronics devices. Nevertheless, little is known about the underlying spin- and orbital-degrees of freedom of the electron wave functions in Weyl semimetals, that govern the electric transport. Here, we give evidence that the chirality of the Weyl points in the Type-II Weyl semimetal MoTe$_2$ is directly linked to the spin texture and orbital angular momentum of the electron wave functions. By means of state-of-the-art spin- and momentum-resolved photoemission spectroscopy the spin- and orbital texture in the Fermi surface is directly resolved. Supported by first-principles calculations, we examined the relationship between the topological chiral charge and spin texture, which significantly contributes to the understanding of the electronic structure in topological quantum materials.

preprint2022arXiv

Long-Range Orbital Magnetoelectric Torque in Ferromagnets

While it is often assumed that the orbital response is suppressed and short-ranged due to strong crystal field potential and orbital quenching, we show that the orbital magnetoelectric response can be remarkably long-ranged in ferromagnets. In a bilayer consisting of a nonmagnet and a ferromagnet, spin injection from the interface results in spin accumulation and torque in the ferromagnet, which rapidly oscillate and decay by spin dephasing. In contrast, we find that even when an external electric field is applied only on the nonmagnet, we find substantially long-ranged orbital magnetoelectric response in the FM, which can go far beyond the spin dephasing length. This unusual feature is attributed to nearly degenerate orbital characters imposed by the crystal symmetry, which form hotspots for the intrinsic orbital response. Because only the states near the hotspots contribute dominantly, the induced orbital angular momentum does not exhibit destructive interference among states with different momentum as in the case of the spin dephasing. This gives rise to a distinct type of orbital torque on the magnetization, increasing with the thickness of the ferromagnet. Such behavior may serve as critical long-sought evidence of orbital transport to be directly tested in experiments. Our findings open the possibility of using long-range orbital magnetoelectric effect in orbitronic device applications.

preprint2022arXiv

Magnetic domain walls of the van der Waals material Fe$_3$GeTe$_2$

Among two-dimensional materials, Fe$_3$GeTe$_2$ has come to occupy a very important place owing to its ferromagnetic nature with one of the highest Curie temperatures among known van der Waals materials and the potential for hosting skyrmions. In this combined experimental and theoretical work, we investigate the magnetic bubble domains as well as the microscopic domain wall profile using spin-polarized scanning tunneling microscopy in combination with atomistic spin-dynamics simulations performed with parameters from density functional theory calculations. We find a weak magneto-electric effect influencing the domain wall width by the electric field in the tunneling junction and determine the critical magnetic field for the collapse of the bubble domains. Our findings shed light on the origins of complex magnetism that Fe$_3$GeTe$_2$ exhibits.

preprint2022arXiv

Non-Collinear Spin Current for Switching of Chiral Magnetic Textures

We propose a concept of non-collinear spin current, whose spin polarization varies in space even in non-magnetic crystals. While it is commonly assumed that the spin polarization of the spin Hall current is uniform, asymmetric local crystal potential generally allows the spin polarization to be non-collinear in space. Based on microscopic considerations we demonstrate that such non-collinear spin Hall currents can be observed for example in layered Kagome Mn$_3$X (X = Ge, Sn) compounds. Moreover, by referring to atomistic spin dynamics simulations we show that non-collinear spin currents can be used to switch the chiral spin texture of Mn$_3$X in a deterministic way even in the absence of an external magnetic field. Our theoretical prediction can be readily tested in experiments, which will open a novel route toward electric control of complex spin structures in non-collinear antiferromagnets.

preprint2022arXiv

Proximity induced superconductivity in a topological insulator

Interfacing a topological insulator (TI) with an $s$-wave superconductor (SC) is a promising material platform that offers the possibility to realize a topological superconductor through which Majorana-based topologically protected qubits can be engineered. In our computational study of the prototypical SC/TI interface between Nb and Bi$_2$Te$_3$, we identify the benefits and possible bottlenecks of this potential Majorana material platform. Bringing Nb in contact with the TI film induces charge doping from the SC to the TI, which shifts the Fermi level into the TI conduction band. For thick TI films, this results in band bending leading to the population of trivial TI quantum-well states at the interface. In the superconducting state, we uncover that the topological surface state experiences a sizable superconducting gap-opening at the SC/TI interface, which is furthermore robust against fluctuations of the Fermi energy. We also show that the trivial interface state is only marginally proximitized, potentially obstructing the realization of Majorana-based qubits in this material platform.

preprint2022arXiv

Reply to 'Comment on "Proper and improper chiral magnetic interactions" '

In our previous Letter [Phys. Rev. B 103, L140408 (2021)], we presented a discussion of the fundamental physical properties of the interactions parameterizing atomistic spin models in connection to first-principles approaches that enable their calculation for a given material. This explained how some of those approaches can apparently lead to magnetic interactions that do not comply with the expected physical properties, such as Dzyaloshinskii-Moriya interactions which are non-chiral and independent of the spin-orbit interaction, and which we consequently termed `improper'. In the preceding Comment [Phys. Rev. B 105, 026401], the authors present arguments based on the distinction between global and local approaches to the mapping of the magnetic energy using first-principles calculations to support their proposed non-chiral Dzyaloshinskii-Moriya interactions and their dismissal of our distinction between `proper' and `improper' magnetic interactions. In this Reply, we identify the missing step in the local approach to the mapping and explain how ignoring this step leads to the identification of magnetic interactions which do not comply with established physical principles and that we have previously termed `improper'.

preprint2022arXiv

Spin-orbit enabled all-electrical readout of chiral spin-textures

Chirality and topology are intimately related fundamental concepts, which are heavily explored to establish spin-textures as potential magnetic bits in information technology. However, this ambition is inhibited since electrical reading of chiral attributes is highly non-trivial with conventional current perpendicular-to-plane (CPP) sensing devices. Here we demonstrate from extensive first-principles simulations and multiple scattering expansion the emergence of the chiral spin-mixing magnetoresistance (C-XMR) enabling highly efficient all-electrical readout of the chirality and helicity of respectively one- and two-dimensional magnetic states of matter. It is linear with spin-orbit coupling in contrast to the quadratic dependence associated with the newly unveiled non-local spin-mixing anisotropic MR (X-AMR). Such transport effects are systematised on various non-collinear magnetic states -- spin-spirals and skyrmions -- and compared to the uncovered spin-orbit-independent multi-site magnetoresistances. Owing to their simple implementation in readily available reading devices, the proposed magnetoresistances offer exciting and decisive ingredients to explore with all-electrical means the rich physics of topological and chiral magnetic objects.

preprint2021arXiv

$Z_2$ topology of bismuth

While first-principles calculations with different levels of sophistication predict a topologically trivial $Z_2$ state for bulk bismuth, some photoemission experiments show surface states consistent with the interpretation of bismuth being in a topologically non-trivial $Z_2$ state. We resolve this contradiction between theory and experiment by showing, based on quasiparticle self-consistent $GW$ calculations, that the experimental surface states interpreted as supporting a non-trivial phase are actually consistent with a trivial $Z_2$ invariant. We identify this contradiction as the result of a crosstalk effect arising from the extreme penetration depth of the surface states into the bulk of Bi. A film of Bi can be considered bulk-like only for thicknesses of about 1000 bilayers ($\approx$ 400 nm) and more.

preprint2021arXiv

Charge and spin photocurrents in the Rashba model

In metallic noncentrosymmetric crystals and at surfaces the response of spin currents and charge currents to applied electric fields contains contributions that are second order in the electric field, which are forbidden by symmetry in centrosymmetric systems. Thereby, photocurrents and spin photocurrents can be generated in inversion asymmetric metals by the application of femtosecond laser pulses. We study the laser-induced charge current in the ferromagnetic Rashba model with in-plane magnetization and find that this \textit{magnetic photogalvanic effect} can be tuned to be comparable in size to the laser-induced photocurrents measured experimentally in magnetic bilayer systems such as Co/Pt. Additionally, we show that femtosecond laser pulses excite strong spin currents in the nonmagnetic Rashba model when the Rashba parameter is large.

preprint2021arXiv

Density functional Bogoliubov-de Gennes analysis of superconducting Nb and Nb(110) surfaces

We report on the implementation of the Bogoliubov-de Gennes method into the JuKKR code [https: //jukkr.fz-juelich.de], an implementation of the relativistic all-electron, full-potential Korringa-Kohn-Rostoker Green function method, which allows a material-specific description of inhomogeneous super-conductors and heterostructures on the basis of density functional theory. We describe the formalism and report on calculations for the s-wave superconductor Nb, a potential component of the materials platform enabling the realization of the Majorana zero modes in the field of topological quantum computing. We compare the properties of the superconducting state both in the bulk and for (110) surfaces. We compare slab calculations for different thicknesses and comment on the importance of spin-orbit coupling, the effect of surface relaxations and the influence of a softening of phonon modes on the surface for the resulting superconducting gap.

preprint2021arXiv

Direct and inverse spin-orbit torques in antiferromagnetic and ferromagnetic FeRh/W(001)

We use \textit{ab-initio} calculations to investigate spin-orbit torques (SOTs) in FeRh(001) deposited on W(100). Since FeRh undergoes a ferromagnetic-antiferromagnetic phase transition close to room temperature, we consider both phases of FeRh. In the antiferromagnetic case we find that the effective magnetic field of the even torque is staggered and therefore ideal to induce magnetization dynamics or to switch the antiferromagnet (AFM). At the antiferromagnetic resonance the inverse SOT induces a current density, which can be determined from the SOT. In the ferromagnetic case our calculations predict both even and odd components of the SOT, which can also be used to describe the ac and dc currents induced at the ferromagnetic resonance. For comparison we compute the SOTs in the c($2\times 2$) AFM state of Fe/W(001).

preprint2021arXiv

Heisenberg representation of nonthermal ultrafast laser excitation of magnetic precessions

We derive the Heisenberg representation of the ultrafast inverse Faraday effect that provides the time evolution of magnetic vectors of a magnetic system during its interaction with a laser pulse. We obtain a time-dependent effective magnetic operator acting in the Hilbert space of the total angular momentum that describes a process of nonthermal excitation of magnetic precessions in an electronic system by a circularly polarized laser pulse. The magnetic operator separates the effect of the laser pulse on the magnetic system from other magnetic interactions. The effective magnetic operator provides the equations of motion of magnetic vectors during the excitation by the laser. We show that magnetization dynamics calculated with these equations is equivalent to magnetization dynamics calculated with the time-dependent Schrödinger equation, which takes into account the interaction of an electronic system with the electric field of light. We model and compare laser-induced precessions of magnetic sublattices of an easy-plane and an easy-axis antiferromagnetic systems. Using these models, we show how the ultrafast inverse Faraday effect induces a net magnetic moment in antiferromagnets and demonstrate that a crystal field environment and the exchange interaction play essential roles for laser-induced magnetization dynamics even during the action of a pump pulse. Using our approach, we show that light-induced precessions can start even during the action of the pump pulse with a duration several tens times shorter than the period of induced precessions and affect the position of magnetic vectors after the action of the pump pulse.

preprint2021arXiv

Laser-induced torques in spin spirals

We investigate laser-induced torques in magnetically non-collinear ferromagnets with a spin-spiral magnetic structure using \textit{ab-initio} calculations. Since spin-spirals may be used to approximate the magnetization gradients locally in domain walls and skyrmions, our method may be used to obtain the laser-induced torques in such objects from a multiscale approach. Employing the generalized Bloch-theorem we obtain the electronic structure computationally efficiently. We employ our method to assess the laser-induced torques in bcc Fe, hcp Co, and L$_{1}0$ FePt when a spin-spiral magnetic structure is imposed. We find that the laser-induced torques in these magnetically noncollinear systems may be orders of magnitude larger than those in the corresponding magnetically collinear systems and that they exist both for linearly and circularly polarized light. This result suggests that laser-induced torques driven by noncollinear magnetic order or by magnetic fluctuations may contribute significantly to processes in ultrafast magnetism.

preprint2021arXiv

Magnetic particle-antiparticle creation and annihilation

A fundamental property of particles and antiparticles, such as electrons and positrons, is their ability to annihilate one another. Similar behavior is predicted for magnetic solitons -- localized spin textures that can be distinguished by their topological index $Q$. Theoretically, magnetic topological solitons with opposite values of $Q$, such as skyrmions and their antiparticles -- antiskyrmions -- are expected to be able to merge continuously and to annihilate. However, experimental verification of such particle-antiparticle pair production and annihilation processes has been lacking. Here, we report the creation and annihilation of skyrmion-antiskyrmion pairs in an exceptionally thin film of the cubic chiral magnet B20-type FeGe observed using transmission electron microscopy. Our observations are highly reproducible and are fully consistent with micromagnetic simulations. Our findings provide a new platform for fundamental studies of particles and antiparticles based on magnetic solids and open new perspectives for practical applications of thin films of isotropic chiral magnets.

preprint2021arXiv

Short period magnetization texture of B20-MnGe explained by thermally fluctuating local moments

B20-type compounds, such as MnSi and FeGe, host helimagnetic and skyrmion phases at the mesoscale, which are canonically explained by the combination of ferromagnetic isotropic interactions with weaker chiral Dzyaloshinskii-Moriya ones. Mysteriously, MnGe evades this paradigm as it displays a noncollinear magnetic state at a much shorter nanometer scale. Here we show that the length scale and volume-dependent magnetic properties of MnGe stem from purely isotropic exchange interactions, generally obtained in the paramagnetic state. Our approach is validated by comparing MnGe with the canonical B20-helimagnet FeGe. The free energy of MnGe is calculated, from which we show how triple-q magnetic states can stabilize by adding higher-order interactions.

preprint2021arXiv

The AiiDA-Spirit plugin for automated spin-dynamics simulations and multi-scale modelling based on first-principles calculations

Landau-Lifshitz-Gilbert (LLG) spin-dynamics calculations based on the extended Heisenberg Hamiltonian is an important tool in computational materials science involving magnetic materials. LLG simulations allow to bridge the gap from expensive quantum mechanical calculations with small unit cells to large supercells where the collective behavior of millions of spins can be studied. In this work we present the AiiDA-Spirit plugin that connects the spin-dynamics code Spirit to the AiiDA framework. AiiDA provides a Python interface that facilitates performing high-throughput calculations while automatically augmenting the calculations with metadata describing the data provenance between calculations in a directed acyclic graph. The AiiDA-Spirit interface thus provides an easy way for high-throughput spin-dynamics calculations. The interface to the AiiDA infrastructure furthermore has the advantage that input parameters for the extended Heisenberg model can be extracted from high-throughput first-principles calculations including a proper treatment of the data provenance that ensures reproducibility of the calculation results in accordance to the FAIR principles. We describe the layout of the AiiDA-Spirit plugin and demonstrate its capabilities using selected examples for LLG spin-dynamics and Monte Carlo calculations. Furthermore, the integration with first-principles calculations through AiiDA is demonstrated at the example of $γ$-Fe, where the complex spin-spiral ground state is investigated.

preprint2020arXiv

Ab-initio Theory of Fourier-transformed Quasiparticle Interference Maps and Application to the Topological Insulator Bi$_2$Te$_3$

The quasiparticle interference (QPI) technique is a powerful tool that allows to uncover the structure and properties of electronic structure of a material combined with scattering properties of defects at surfaces. Recently this technique has been pivotal in proving the unique properties of the surface state of topological insulators which manifests itself in the absence of backscattering. In this work we derive a Green function based formalism for the ab initio computation of Fourier-transformed QPI images. We show the efficiency of our new implementation at the examples of QPI that forms around magnetic and non-magnetic defects at the Bi$_2$Te$_3$ surface. This method allows a deepened understanding of the scattering properties of topologically protected electrons off defects and can be a useful tool in the study of quantum materials in the future.

preprint2020arXiv

Atomistic perspective of long lifetimes of small skyrmions at room temperature

The current development to employ magnetic skyrmions in novel spintronic device designs has led to a demand for room temperature-stable skyrmions of ever smaller size. We present extensive studies on skyrmion stability in atomistic magnetic systems in two- and three-dimensional geometries. We show that for materials described by the same micromagnetic parameters, the variation of the atomistic exchange between different neighbors, the stacking order, and the number of layers of the atomic lattice can significantly influence the rate of the thermally activated decay of a skyrmion. These factors alone are important considerations, but it is shown that their combination can open up novel avenues of materials design in the search for sub-10nm skyrmions, as their lifetime can be extended by several orders of magnitude.

preprint2020arXiv

Effective Seiberg-Witten gauge theory of noncollinear magnetism

Smoothly varying magnetization textures such as domain walls, skyrmions or hopfions serve as promising candidates for the information bits of the future. Understanding their physical properties is both a major field of interest and a theoretical challenge, involving the physics on different length scales. Here, we apply the phase space formulation of quantum mechanics to magnetic insulators and metals in the limit of zero temperature to obtain a gradient expansion in terms of real-space derivatives of the magnetization. Our primary focus is the anomalous Hall effect in noncollinear magnets which serves as an important proxy in the detection of localized magnetic structures. We formulate the problem in the language of noncommutative fiber bundles and make the central finding that the semiclassical expansion of the density matrix and the Berry curvature is governed by a construction from string theory which is known as the Seiberg-Witten map. Originally discovered in the effective low-energy behavior of D-branes, this map now gives a geometrical underpinning to gradient expansion techniques in noncollinear magnets and offers a radically new perspective on their electronic properties.

preprint2020arXiv

Electric dipole moment as descriptor for interfacial Dzyaloshinskii-Moriya interaction

Chiral magnets are of fundamental interest and have important technological ramifications. The origin of chiral magnets lies in the Dzyaloshinskii-Moriya interaction (DMI), an interaction whose experimental and theoretical determination is laborious. We derive an expression that identifies the electric dipole moment as descriptor for the systematic design of chiral magnetic multilayers. Using density functional theory calculations, we determine the DMI of (111)-oriented metallic ferromagnetic $Z$/Co/Pt multilayers of ultrathin films. The non-magnetic layer $Z$ determines the DMI at the Co-Pt interface. The results validate the electric and magnetic dipole moments as excellent descriptors. We found a linear relation between the electric dipole moment of Pt, the Allen electronegativity of $Z$, and the contribution of Pt to the total DMI.

preprint2020arXiv

First-principles study of the electronic and magnetic properties of cubic GdCu compound

The structural, electronic, and magnetic properties of bulk GdCu (CsCl-type) are investigated using spin density functional theory, where highly localized $4f$ orbitals are treated within LDA+$U$ and GGA+$U$ methods. The calculated magnetic ground state of GdCu using collinear as well as spin spiral calculations exhibits a C-type antiferromagnetic configuration representing a spin spiral propagation vector $\mathbf{Q}=\frac{2π}{a}(\frac{1}{2},\frac{1}{2},0)$. The parameters of the effective Heisenberg Hamiltonian are evaluated from a self-consistent electronic structure and are used to determine the magnetic transition temperature. The estimated Néel temperature of the cubic GdCu using GGA+$U$ and LDA+$U$ density functionals within the mean field and random phase approximations are in good agreement with the experimentally measured values. In particular, the theoretical understanding of the experimentally observed core Gd $4f$ levels shifting in photoemission spectroscopy experiments is investigated in detail. By employing the self-consistent constrained random-phase approximation we determined the strength of the effective Coulomb interaction (Hubbard $U$) between localized $4f$ electrons. We find that, the shift of Gd-$4f$ states in GdCu with respect to bulk Gd within DFT+$U$ is sensitive to choice of lattice parameter. The calculations for $4f$-level shifts using DFT+$U$ methods as well as Hubbard-1 approximation are not consistent with the experimental findings.

preprint2020arXiv

Giant Anomalous Nernst Effect in Noncollinear Antiferromagnetic Mn-based Antiperovskite Nitrides

The anomalous Nernst effect (ANE) - the generation of a transverse electric voltage by a longitudinal heat current in conducting ferromagnets or antiferromagnets - is an appealing approach for thermoelectric power generation in spin caloritronics. The ANE in antiferromagnets is particularly convenient for the fabrication of highly efficient and densely integrated thermopiles as lateral configurations of thermoelectric modules increase the coverage of heat source without suffering from the stray fields that are intrinsic to ferromagnets. In this work, using first-principles calculations together with a group theory analysis, we systematically investigate the spin order-dependent ANE in noncollinear antiferromagnetic Mn-based antiperovskite nitrides Mn$_{3}X$N ($X$ = Ga, Zn, Ag, and Ni). The ANE in Mn$_{3}X$N is forbidden by symmetry in the R1 phase but amounts to its maximum value in the R3 phase. Among all Mn$_{3}X$N compounds, Mn$_{3}$NiN presents the most significant anomalous Nernst conductivity of 1.80 AK$^{-1}$m$^{-1}$ at 200 K, which can be further enhanced if strain, electric, or magnetic fields are applied. The ANE in Mn$_{3}$NiN, being one order of magnitude larger than that in the famous Mn$_{3}$Sn, is the largest one discovered in antiferromagnets so far. The giant ANE in Mn$_{3}$NiN originates from the sharp slope of the anomalous Hall conductivity at the Fermi energy, which can be understood well from the Mott relation. Our findings provide a novel host material for realizing antiferromagnetic spin caloritronics which promises exciting applications in energy conversion and information processing.

preprint2020arXiv

Material systems for FM-/AFM-coupled skyrmions in Co/Pt-based multilayers

By means of systematic first-principles calculations based on density functional theory we search for suitable materials that can host antiferromagnetically coupled skyrmions. We concentrate on fcc-stacked (111)-oriented metallic $Z$/Co/Pt ($Z=4d$ series: Y$-$Pd, the noble metals: Cu, Ag, Au, post noble metals: Zn and Cd) magnetic multilayers of films of monatomic thickness. We present quantitative trends of magnetic properties: Magnetic moments, interlayer exchange coupling, spin-stiffness, Dzyaloshinskii-Moriya interaction, magnetic anisotropy, and the critical temperature. We show that some of the $Z$ elements (Zn, Y, Zr, Nb, Tc, Ru, Rh, and Cd) can induce antiferromagnetic interlayer coupling between the magnetic Co layers, and that they influence the easy magnetization axis. Employing a multiscale approach, we transfer the micromagnetic parameters determined from $ab$ $initio$ to a micromagnetic energy functional and search for one-dimensional spin-spiral solutions and two-dimensional skyrmions. We determine the skyrmion radius by numerically solving the equation of the skyrmion profile. We found an analytical expression for the skyrmion radius that covers our numerical results and is valid for a large regime of micromagnetic parameters. Based on this expression we have proposed a model that allows to extrapolate from the $ab$ $initio$ results of monatomic films to multilayers with Co films consisting of several atomic layers containing $10\,$nm skyrmions. We found thickness regimes where tiny changes of the film thickness may alter the skyrmion radius by orders of magnitude. We estimated the skyrmion size as function of temperature and found that the size can easily double going from cryogenic to room temperature. We suggest promising material systems for ferromagnetically and antiferromagnetically coupled spin-spiral and skyrmion systems.

preprint2020arXiv

Plumbene on a magnetic substrate: a combined STM and DFT study

As heavy analog of graphene, plumbene is a two-dimensional material with strong spin-orbit coupling effects. Using scanning tunneling microscopy (STM), we observe that Pb forms a flat honeycomb lattice on an Fe monolayer on Ir(111). In contrast, without the Fe layer, a c(2x4) structure of Pb on Ir(111) is found. We use density functional theory (DFT) calculations to rationalize these findings and analyze the impact of the hybridization on the plumbene band structure. In the unoccupied states the splitting of the Dirac cone by spin-orbit interaction is clearly observed while in the occupied states of the freestanding plumbene we find a band inversion that leads to the formation of a topologically non-trivial gap. Exchange splitting as mediated by the strong hybridization with the Fe layer drives a quantum spin Hall to quantum anomalous Hall state transition.

preprint2020arXiv

Solution to the modified Helmholtz equation for arbitrary periodic charge densities

We present a general method for solving the modified Helmholtz equation without shape approximation for an arbitrary periodic charge distribution, whose solution is known as the Yukawa potential or the screened Coulomb potential. The method is an extension of Weinert's pseudo-charge method [M. Weinert, J. Math. Phys. 22, 2433 (1981)] for solving the Poisson equation for the same class of charge density distributions. The inherent differences between the Poisson and the modified Helmholtz equation are in their respective radial solutions. These are polynomial functions, for the Poisson equation, and modified spherical Bessel functions, for the modified Helmholtz equation. This leads to a definition of a modified pseudo-charge density and modified multipole moments. We have shown that Weinert's convergence analysis of an absolutely and uniformly convergent Fourier series of the pseudo-charge density is transferred to the modified pseudo-charge density. We conclude by illustrating the algorithmic changes necessary to turn an available implementation of the Poisson solver into a solver for the modified Helmholtz equation.

preprint2020arXiv

Theory of Current-Induced Angular Momentum Transfer Dynamics in Spin-Orbit Coupled Systems

Motivated by the importance of understanding competing mechanisms to current-induced spin-orbit torque in complex magnets, we develop a unified theory of current-induced spin-orbital coupled dynamics. The theory describes angular momentum transfer between different degrees of freedom in solids, e.g., the electron orbital and spin, the crystal lattice, and the magnetic order parameter. Based on the continuity equations for the spin and orbital angular momenta, we derive equations of motion that relate spin and orbital current fluxes and torques describing the transfer of angular momentum between different degrees of freedom. We then propose a classification scheme for the mechanisms of the current-induced torque in magnetic bilayers. Based on our first-principles implementation, we apply our formalism to two different magnetic bilayers, Fe/W(110) and Ni/W(110), which are chosen such that the orbital and spin Hall effects in W have opposite sign and the resulting spin- and orbital-mediated torques can compete with each other. We find that while the spin torque arising from the spin Hall effect of W is the dominant mechanism of the current-induced torque in Fe/W(110), the dominant mechanism in Ni/W(110) is the orbital torque originating in the orbital Hall effect of W. It leads to negative and positive effective spin Hall angles, respectively, which can be directly identified in experiments. This clearly demonstrates that our formalism is ideal for studying the angular momentum transfer dynamics in spin-orbit coupled systems as it goes beyond the "spin current picture" by naturally incorporating the spin and orbital degrees of freedom on an equal footing. Our calculations reveal that, in addition to the spin and orbital torque, other contributions such as the interfacial torque and self-induced anomalous torque within the ferromagnet are not negligible in both material systems.

preprint2020arXiv

Topological magneto-optical effects and their quantization in noncoplanar antiferromagnets

Reflecting the fundamental interactions of polarized light with magnetic matter, magneto-optical effects are well known since more than a century. The emergence of these phenomena is commonly attributed to the interplay between exchange splitting and spin-orbit coupling in the electronic structure of magnets. Using theoretical arguments, we demonstrate that topological magneto-optical effects can arise in noncoplanar antiferromagnets due to the finite scalar spin chirality, without any reference to exchange splitting or spin-orbit coupling. We propose spectral integrals of certain magneto-optical quantities that uncover the unique topological nature of the discovered effect. We also find that the Kerr and Faraday rotation angles can be quantized in insulating topological antiferromagnets in the low-frequency limit, owing to nontrivial global properties that manifest in quantum topological magneto-optical effects. Although the predicted topological and quantum topological magneto-optical effects are fundamentally distinct from conventional light-matter interactions, they can be measured by readily available experimental techniques.

preprint2020arXiv

Topological-chiral magnetic interactions driven by emergent orbital magnetism

Two hundred years ago, André-Marie Ampère discovered that electric loops in which currents of electrons are generated by a penetrating magnetic field can interact with each other. Here we show that Ampère's observation can be transferred to the quantum realm of interactions between triangular plaquettes of spins on a lattice, where the electrical currents at the atomic scale are associated with a peculiar type of the orbital motion of electrons in response to the non-coplanarity of neighbouring spins playing the role of a magnetic field. The resulting topological orbital moment underlies the relation of the orbital dynamics with the topology of the spin structure. We demonstrate that the interactions of the topological orbital moments with each other and with the spins of the underlying lattice give rise to a new class of magnetic interactions $-$ topological chiral interactions $-$ which can dominate over the celebrated Dzyaloshinskii-Moriya interaction, thus opening a path for the realization of new classes of chiral magnetic materials with three-dimensional magnetization textures such as magnetic hopfions.

preprint2019arXiv

Systematic derivation of realistic spin-models for beyond-Heisenberg solids

We present a systematic derivation of effective lattice spin Hamiltonians derived from a rotationally invariant multi-orbital Hubbard model including a term ensuring Hund's rule coupling. The Hamiltonians are derived down-folding the fermionic degrees of freedom of the Hubbard model into the proper low-energy spin sector using Löwdin partitioning, which will be outlined in detail for the case of two sites and two orbitals at each site. Correcting the ground state systematically up to fourth order in the hopping of electrons, we find for spin $S\geq 1$ a biquadratic, three-spin and four-spin interaction beyond the conventional Heisenberg term. Comparing the puzzling energy spectrum of the magnetic states for a single Fe monolayer on Ru(0001), obtained from density functional theory, with the spin Hamiltonians taken at the limit of classical spins, we show that the previously ignored three-spin interaction can be comparable in size to the conventional Heisenberg exchange.

preprint2019arXiv

The chiral Hall effect of magnetic skyrmions from a cyclic cohomology approach

We demonstrate the emergence of an anomalous Hall effect in chiral magnetic textures which is neither proportional to the net magnetization nor to the well-known emergent magnetic field that is responsible for the topological Hall effect. Instead, it appears already at linear order in the gradients of the magnetization texture and exists for one-dimensional magnetic textures such as domain walls and spin spirals. It receives a natural interpretation in the language of Alain Connes' noncommutative geometry. We show that this chiral Hall effect resembles the familiar topological Hall effect in essential properties while its phenomenology is distinctly different. Our findings make the re-interpretation of experimental data necessary, and offer an exciting twist in engineering the electrical transport through magnetic skyrmions.

preprint2019arXiv

Wannier90 as a community code: new features and applications

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials.