Source author record

Stefan Blügel

Stefan Blügel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.str-el cond-mat.other physics.comp-ph cond-mat.supr-con cond-mat.stat-mech math-ph math.MP nlin.PS physics.app-ph quant-ph

Catalog footprint

What is connected

125works

12topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

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.

preprint2016arXiv

Approximations to the exact exchange potential: KLI versus semilocal

In the search for an accurate and computationally efficient approximation to the exact exchange potential of Kohn-Sham density functional theory, we recently compared various semilocal exchange potentials to the exact one [F. Tran et al., Phys. Rev. B 91, 165121 (2015)]. It was concluded that the Becke-Johnson (BJ) potential is a very good starting point, but requires the use of empirical parameters to obtain good agreement with the exact exchange potential. In the present work, we extend the comparison by considering the Krieger-Li-Iafrate (KLI) approximation, which is a beyond-semilocal approximation. It is shown that overall the KLI and BJ-based potentials are the most reliable approximations to the exact exchange potential, however, sizeable differences, especially for the antiferromagnetic transition-metal oxides, can be obtained.

preprint2016arXiv

Chirality-driven orbital magnetic moments as a new probe for topological magnetic structures

When electrons are driven through unconventional magnetic structures, such as skyrmions, they experience emergent electromagnetic fields that originate several Hall effects. Independently, ground state emergent magnetic fields can also lead to orbital magnetism, even without the spin-orbit interaction. The close parallel between the geometric theories of the Hall effects and of the orbital magnetization raises the question: does a skyrmion display topological orbital magnetism? Here we first address the smallest systems with nonvanishing emergent magnetic field, trimers, characterizing the orbital magnetic properties from first-principles. Armed with this understanding, we study the orbital magnetism of skyrmions, and demonstrate that the contribution driven by the emergent magnetic field is topological. This means that the topological contribution to the orbital moment does not change under continous deformations of the magnetic structure. Furthermore, we use it to propose a new experimental protocol for the identification of topological magnetic structures, by soft x-ray spectroscopy.

preprint2016arXiv

Fermi surfaces, spin-mixing parameter, and colossal anisotropy of spin relaxation in transition metals from ab initio theory

The Fermi-surfaces and Elliott-Yafet spin-mixing parameter (EYP) of several elemental metals are studied by \emph{ab initio} calculations. We focus first on the anisotropy of the EYP as a function of the direction of the spin-quantization axis [Phys.~Rev.~Lett.\ \textbf{109}, 236603 (2012)]. We analyze in detail the origin of the gigantic anisotropy in $5d$ hcp metals as compared to $5d$ cubic metals by band-structure calculations and discuss the stability of our results against an applied magnetic field. We further present calculations of light (4$d$ and 3$d$) hcp crystals, where we find a huge increase of the EYP anisotropy, reaching colossal values as large as $6000\%$ in hcp Ti. We attribute these findings to the reduced strength of spin-orbit coupling, which promotes the anisotropic spin-flip hot loops at the Fermi surface. In order to conduct these investigations, we developed an adapted tetrahedron-based method for the precise calculation of Fermi surfaces of complicated shape and accurate Fermi-surface integrals within the full-potential relativistic Korringa-Kohn-Rostoker Green-function method.

preprint2016arXiv

Giant spin Nernst effect induced by resonant scattering at surfaces of metallic films

A new concept realizing giant spin Nernst effect in nonmagnetic metallic films is introduced. It is based on the idea of engineering an asymmetric energy dependence of the longitudinal and transverse electrical conductivities, as well as a pronounced energy dependence of the spin Hall angle in the vicinity of the Fermi level by the resonant impurity states at the Fermi level. We employ an analytical model and demonstrate the emergence of a giant spin Nernst effect in Ag(111) films using {\it ab-initio} calculations combined with the Boltzmann approach for transport properties arising from skew scattering off impurities.

preprint2016arXiv

Influence of complex disorder on skew-scattering Hall effects in $L1_0$-ordered FePt alloy

We show by first-principles calculations that the skew-scattering anomalous Hall and spin-Hall angles of L$1_0$-ordered FePt drastically depend on different types of disorder. A different sign of the AHE is obtained when slightly deviating from the stoichiometric ratio towards the Fe-rich side as compared to the Pt-rich side. For stoichiometric samples, short-range ordering of defects has a profound effect on the Hall angles and can change them by a factor of $2$ as compared to the case of uncorrelated disorder. This might explain the vast range of anomalous Hall angles measured in experiments, which undergo different preparation procedures and thus might differ in their crystallographic quality.

preprint2016arXiv

Interlayer exchange coupling -- a general scheme turning chiral magnets into magnetic multilayers carrying atomic-scale skyrmions

We report on a general principle using the interlayer exchange coupling to extend the regime of chiral magnetic films in which stable or metastable magnetic skyrmions can appear at zero magnetic field. We verify this concept on the basis of a first-principles model for a Mn monolayer on W(001) substrate, a prototype chiral magnet for which the atomic-scale magnetic texture is determined by the frustration of exchange interactions, impossible to unwind by laboratory magnetic fields. By means of \textit{ab initio} calculations for the Mn/W$_m$/Co$_n$/Pt/W(001) multilayer system we show that for certain thicknesses $m$ of the W spacer and $n$ of the Co reference layer, the effective field of the reference layer fully substitutes the required magnetic field for skyrmion formation.

preprint2016arXiv

Laser-induced torques in metallic ferromagnets

We study laser-induced torques in bcc Fe, hcp Co and $L1_0$ FePt based on first-principles electronic structure calculations and the Keldysh nonequilibrium formalism. We find that the torques have two contributions, one from the inverse Faraday effect (IFE) and one from the optical spin-transfer torque (OSTT). Depending on the ferromagnet at hand and on the quasiparticle broadening the two contributions may be of similar magnitude or one contribution may dominate over the other. Additionally, we determine the nonequilibrium spin polarization in order to investigate its relation to the torque. We find the torques and the perpendicular component of the nonequilibrium spin polarization to be odd in the helicity of the laser light, while the spin polarization that is induced parallel to the magnetization is helicity-independent. The parallel component of the nonequilibrium spin polarization is orders of magnitude larger than the perpendicular component. In the case of hcp Co we find good agreement between the calculated laser-induced torque and a recent experiment.

preprint2016arXiv

Microscopic theory of the residual surface resistivity of Rashba electrons

A microscopic expression of the residual electrical resistivity tensor is derived in linear response theory for Rashba electrons scattering at a magnetic impurity with cylindrical or non-cylindrical potential. The behavior of the longitudinal and transversal residual resistivity is computed for an Fe impurity at the Au(111) surface. We studied the evolution of the resistivity tensor elements as function of the Rashba spin-orbit strength and the magnetization direction of the impurity. We found that the absolute values of longitudinal resistivity reduces with increasing spin-orbit strength of the substrate and that the scattering of the conduction electrons at magnetic impurities with magnetic moments pointing in directions not perpendicular to the surface plane produce a planar Hall effect and an anisotropic magnetoresistance even if the impurity carries no spin-orbit interaction. Functional forms are provided describing the anisotropy of the planar Hall effect and the anisotropic magnetoresistance with respect to the direction of the impurity moment. In the limit of no spin-orbit interaction and a non-magnetic impurity of cylindrical symmetry, the expression of the residual resistivity of a two-dimensional electron-gas has the same simplicity and form as for the three-dimensional electron gas [J. Friedel, NuovoCimento Suppl. 7, 287 (1958)] and can also be expressed in terms of scattering phase shifts.

preprint2016arXiv

Report on imaging of individual skyrmions of MML systems made by different techniques

Deliverable 1.1 is a report on imaging of individual skyrmions of MML systems made by different techniques. This deliverable falls within work package 1 (WP1:Static equilibrium properties of interface induced skyrmions in magnetic films and multilayers) of the MAGicSky research proposal whose objectives are: the growth of magnetic multilayer systems with large interfacial chiral interaction, the determination of the amplitude of the Dzyaloshinskii--Moriya interactions (DMi), the detailed structural and magnetic characterization of interfaces, the magnetic imaging of individual skyrmions and the thermal stability and interaction of skyrmions with intrinsic and/or artificial defects. It should be emphasized that imaging magnetic skyrmionic structures forms the basis for all the planned subsequent material optimization, electrical transport experiments and devices preparation in the MAGicSky research program. In particular, imaging is paramount as a technique since it enables us to unequivocally demonstrate chiral magnetic skyrmions by directly imaging i) the effects of engineering multilayer samples with large DMis and ii) the skyrmionic structure and properties. In order to get complementary physical insight on the properties of magnetic skyrmions and to corroborate our research results, we utilise a plethora of key state-of-the-art advanced magnetic imaging techniques, namely: i) Spin-polarised Scanning Tunnelling Microscopy (SP-STM), Scanning Transmission X-ray Microscopy (STXM) Lorentz Transmission Electron Microscopy (Lorentz TEM) as well as, Magnetic Force Microscopy (MFM). These imaging experiments can therefore stand as the base for all the other objectives in WP1 as well as WP2 that will study Dynamic properties: nucleation, propagation and detection of skyrmions. Imaging has been aimed to test and get feedback from the theoretical and numerical investigations of skyrmion systems of WP3. Finally, for WP4: Towards skyrmion based devices for information storage and processing, some of the techniques detailed here in D1.1 will enable us to image and test the first devices and their functionality. In short, the MAGicSky consortium research results in D1.1 constitute a major milestone that acts as a stepping-stone to every other experimental effort and we have achieved it. The results relevant to the D1.1 deliverable were presented in 5 publications and several invited talks at international workshop and conferences:

preprint2016arXiv

Role of Berry phase theory for describing orbital magnetism: From magnetic heterostructures to topological orbital ferromagnets

We address the importance of the modern theory of orbital magnetization for spintronics. Based on an all-electron first-principles approach, we demonstrate that the predictive power of the routinely employed "atom-centered" approximation is limited to materials like elemental bulk ferromagnets, while the application of the modern theory of orbital magnetization is crucial in chemically or structurally inhomogeneous systems such as magnetic thin films, and materials exhibiting non-trivial topology in reciprocal and real space,~e.g.,~Chern insulators or non-collinear systems. We find that the modern theory is particularly crucial for describing magnetism in a class of materials that we suggest here $-$ topological orbital ferromagnets.

preprint2016arXiv

Spin-orbit torques and spin accumulation in FePt/Pt and Co/Cu thin films from first principles: the role of impurities

Using the Boltzmann formalism based on the first principles electronic structure and scattering rates, we investigate the current-induced spin accumulation and spin-orbit torques in FePt/Pt and Co/Cu bilayers in the presence of substitutional impurities. In FePt/Pt bilayers we consider the effect of intermixing of Fe and Pt atoms in the FePt layer, and find a crucial dependence of spin accumulation and spin-orbit torques on the details of the distribution of these defects. In Co/Cu bilayers we predict that the magnitude and sign of the spin-orbit torque and spin accumulation depend very sensitively on the type of the impurities used to dope the Cu substrate. Moreover, simultaneously with impurity-driven scattering we consider the effect of an additional constant quasiparticle broadening of the states at the Fermi surface to simulate phonon scattering at room temperature, and discover that even a small broadening of the order of 25 meV can drastically influence the magnitude of the considered effects. We explain our findings based on the analysis of the complex interplay of several competing Fermi surface contributions to the spin accumulation and spin-orbit torques in these structurally and chemically non-uniform systems.

preprint2016arXiv

Strong spin-orbit fields and Dyakonov-Perel spin dephasing in supported metallic films

Spin dephasing by the Dyakonov-Perel mechanism in metallic films deposited on insulating substrates is revealed, and quantitatively examined by means of density functional calculations combined with a kinetic equation. The surface-to-substrate asymmetry, probed by the metal wave functions in thin films, is found to produce strong spin-orbit fields and a fast Larmor precession, giving a dominant contribution to spin decay over the Elliott-Yafet spin relaxation up to a thickness of 70 nm. The spin dephasing is oscillatory in time with a rapid (sub-picosecond) initial decay. However, parts of the Fermi surface act as spin traps, causing a persistent tail signal lasting 1000 times longer than the initial decay time. It is also found that the decay depends on the direction of the initial spin polarization, resulting in a spin-dephasing anisotropy of 200% in the examined cases.

preprint2016arXiv

Sub-molecular modulation of a 4f driven Kondo resonance by surface-induced asymmetry

Coupling between a magnetic impurity and an external bath can give rise to many-body quantum phenomena, including Kondo and Hund's impurity states in metals, and Yu-Shiba-Rusinov states in superconductors. While advances have been made in probing the magnetic properties of d-shell impurities on surfaces, the confinement of f orbitals makes them difficult to access directly. Here we show that a 4f driven Kondo resonance can be modulated spatially by asymmetric coupling between a metallic surface and a molecule containing a 4f-like moment. Strong hybridisation of dysprosium double-decker phthalocyanine with Cu(001) induces Kondo screening of the central magnetic moment. Misalignment between the symmetry axes of the molecule and the surface induces asymmetry in the molecule's electronic structure, spatially mediating electronic access to the magnetic moment through the Kondo resonance. This work demonstrates the important role that molecular ligands play in mediating electronic and magnetic coupling and in accessing many-body quantum states.

preprint2016arXiv

Switching of chiral magnetic skyrmions by picosecond magnetic field pulses via transient topological states

Magnetic chiral skyrmions are vortex like spin structures that appear as stable or meta-stable states in magnetic materials due to the interplay between the symmetric and antisymmetric exchange interactions, applied magnetic field and/or uniaxial anisotropy. Their small size and internal stability make them prospective objects for data storage but for this, the controlled switching between skyrmion states of opposite polarity and topological charge is essential. Here we present a study of magnetic skyrmion switching by an applied magnetic field pulse based on a discrete model of classical spins and atomistic spin dynamics. We found a finite range of coupling parameters corresponding to the coexistence of two degenerate isolated skyrmions characterized by mutually inverted spin structures with opposite polarity and topological charge. We demonstrate how for a wide range of material parameters a short inclined magnetic field pulse can initiate the reliable switching between these states at GHz rates. Detailed analysis of the switching mechanism revealed the complex path of the system accompanied with the excitation of a chiral-achiral meron pair and the formation of an achiral skyrmion.

preprint2016arXiv

The inverse thermal spin-orbit torque and the relation of the Dzyaloshinskii-Moriya interaction to ground-state energy currents

Using the Kubo linear-response formalism we derive expressions to calculate the heat current generated by magnetization dynamics in magnets with broken inversion symmetry and spin-orbit interaction (SOI). The effect of producing heat currents by magnetization dynamics constitutes the Onsager reciprocal of the thermal spin-orbit torque (TSOT), i.e., the generation of torques on the magnetization due to temperature gradients. We find that the energy current driven by magnetization dynamics contains a contribution from the Dzyaloshinskii-Moriya interaction (DMI), which needs to be subtracted from the Kubo linear response of the energy current in order to extract the heat current. We show that the expressions of the DMI coefficient can be derived elegantly from the DMI energy current. Guided by formal analogies between the Berry phase theory of DMI on the one hand and the modern theory of orbital magnetization on the other hand we are led to an interpretation of the latter in terms of energy currents as well. Based on \textit{ab-initio} calculations we investigate the heat current driven by magnetization dynamics in Mn/W(001) magnetic bilayers. We predict that fast domain walls drive strong ITSOT heat currents.

preprint2015arXiv

Comparison between exact and semilocal exchange potentials: An all-electron study for solids

The exact-exchange (EXX) potential, which is obtained by solving the optimized-effective potential (OEP) equation, is compared to various approximate semilocal exchange potentials for a set of selected solids (C, Si, BN, MgO, Cu$_{2}$O, and NiO). This is done in the framework of the linearized augmented plane-wave method, which allows for a very accurate all-electron solution of electronic structure problems in solids. In order to assess the ability of the semilocal potentials to approximate the EXX-OEP, we considered the EXX total energy, electronic structure, electric-field gradient, and magnetic moment. An attempt to parameterize a semilocal exchange potential is also reported.

preprint2015arXiv

Direct and inverse spin-orbit torques

In collinear magnets lacking inversion symmetry application of electric currents induces torques on the magnetization and conversely magnetization dynamics induces electric currents. The two effects, which both rely on spin-orbit interaction (SOI), are reciprocal to each other and denoted direct spin-orbit torque (SOT) and inverse spin-orbit torque (ISOT), respectively. We derive expressions for SOT and ISOT within the Kubo linear response formalism. We show that expressions suitable for density-functional theory calculations can be derived either starting from a Kohn-Sham Hamiltonian with time-dependent exchange field or by expressing general susceptibilities in terms of the Kohn-Sham susceptibilities. For the case of magnetic bilayer systems we derive the general form of the ISOT current induced under ferromagnetic resonance. Using \textit{ab initio} calculations within density-functional theory we investigate SOT and ISOT in Co/Pt(111) magnetic bilayers. We determine the spatial distribution of spin and charge currents as well as torques in order to expose the mechanisms underlying SOT and ISOT and to highlight their reciprocity on the microscopic level. We find that the spin Hall effect is position-dependent close to interfaces.

preprint2015arXiv

Electronic phase transitions of bismuth under strain from relativistic self-consistent GW calculations

We present quasiparticle self-consistent GW (QSGW) calculations of semimetallic bulk Bi. We go beyond the conventional QSGW method by including the spin-orbit coupling throughout the self-consistency cycle. This approach improves the description of the electron and the hole pockets considerably with respect to standard density functional theory (DFT), leading to excellent agreement with experiment. We employ this relativistic QSGW approach to conduct a study of the semimetal-to-semiconductor and the trivial-to-topological transitions that Bi experiences under strain. DFT predicts that an unphysically large strain is needed for such transitions. We show, by means of the relativistic QSGW description of the electronic structure, that an in-plane tensile strain of only 0.3% and a compressive strain of 0.4% are sufficient to cause the semimetal-to-semiconductor and the trivial-to-topological phase transitions, respectively. Thus, the required strain moves into a regime that is likely to be realizable in experiment, which opens up the possibility to explore bulklike topological behavior of pure Bi.

preprint2015arXiv

Higher dimensional Wannier functions of multi-parameter Hamiltonians

When using Wannier functions to study the electronic structure of multi-parameter Hamiltonians $H^{(\boldsymbol k,\bf λ)}$ carrying a dependence on crystal momentum $\boldsymbol k$ and an additional periodic parameter $\bfλ$, one usually constructs several sets of Wannier functions for a set of values of $\bfλ$. We present the concept of higher dimensional Wannier functions (HDWFs), which provide a minimal and accurate description of the electronic structure of multi-parameter Hamiltonians based on a single set of HDWFs. The obstacle of non-orthogonality of Bloch functions at different $\bfλ$ is overcome by introducing an auxiliary real space, which is reciprocal to the parameter $\bfλ$. We derive a generalized interpolation scheme and emphasize the essential conceptual and computational simplifications in using the formalism, for instance, in the evaluation of linear response coefficients. We further implement the necessary machinery to construct HDWFs from ab initio within the full-potential linearized augmented plane-wave method (FLAPW). We apply our implementation to accurately interpolate the Hamiltonian of a one-dimensional magnetic chain of Mn atoms in two important cases of $\bfλ$: (i) the spin-spiral vector $\boldsymbol q$, and (ii) the direction of the ferromagnetic magnetization $\hat{\boldsymbol m}$. Using the generalized interpolation of the energy, we extract the corresponding values of magneto-crystalline anisotropy energy, Heisenberg exchange constants, and spin stiffness, which compare very well with the values obtained from direct first principles calculations. For toy models we demonstrate that the method of HDWFs can also be used in applications such as the virtual crystal approximation, ferroelectric polarization and spin torques.

preprint2015arXiv

Magnetic properties of 2D nickel nanostrips: structure dependent magnetic anomaly

We have investigated different geometries of two dimensional (2D) infinite length Ni nanowires of increasing width using spin density functional theory calculations. Our simulations demonstrate that the parallelogram motif is the most stable and structures that incorporate the parallelogram motif are more stable as compared to rectangular structures. The wires are conducting and the conductance channels increase with increasing width. The wires have a non-linear behavior in the ballistic anistropic magnetoresistance ratios with respect to the magnetization directions. All 2D nanowires as well as Ni (111) and Ni (100) monolayer investigated are ferromagnetic under the Stoner criterion and exhibit enhanced magnetic moments as compared to bulk Ni and the respective Ni monolayers. The Stoner parameter is seen to depend on the structure and the dimension of the Nws. The easy axis for all nickel nanowires under investigation is observed to be along the wire axis. The double rectangular nanowire exhibits a magnetic anomaly with a smaller magnetic moment when compared to Ni (100) monolayer and is the only structure with an easy axis perpendicular to the wire axis.

preprint2015arXiv

Oxygen-enabled control of Dzyaloshinskii-Moriya Interaction in ultra-thin magnetic films

The search for chiral magnetic textures in systems lacking spatial inversion symmetry has attracted a massive amount of interest in the recent years with the real space observation of novel exotic magnetic phases such as skyrmions lattices, but also domain walls and spin spirals with a defined chirality. The electrical control of these textures offers thrilling perspectives in terms of fast and robust ultrahigh density data manipulation. A powerful ingredient commonly used to stabilize chiral magnetic states is the so-called Dzyaloshinskii-Moriya interaction (DMI) arising from spin-orbit coupling in inversion asymmetric magnets. Such a large antisymmetric exchange has been obtained at interfaces between heavy metals and transition metal ferromagnets, resulting in spin spirals and nanoskyrmion lattices. Here, using relativistic first-principles calculations, we demonstrate that the magnitude and sign of DMI can be entirely controlled by tuning the oxygen coverage of the magnetic film, therefore enabling the smart design of chiral magnetism in ultra-thin films. We anticipate that these results extend to other electronegative ions and suggest the possibility of electrical tuning of exotic magnetic phases.

preprint2015arXiv

Quasi 2D electronic states with high spin-polarization in centrosymmetric MoS$_2$ bulk crystals

Time reversal dictates that nonmagnetic, centrosymmetric crystals cannot be spin-polarized as a whole. However, it has been recently shown that the electronic structure in these crystals can in fact show regions of high spin-polarization, as long as it is probed locally in real and in reciprocal space. In this article we present the first observation of this type of compensated polarization in MoS$_2$ bulk crystals. Using spin- and angle-resolved photoemission spectroscopy (ARPES) we directly observed a spin-polarization of more than 65% for distinct valleys in the electronic band structure. By additionally evaluating the probing depth of our method we find that these valence band states at the $\overline{\text{K}}$ point in the Brillouin zone are close to fully polarized for the individual atomic trilayers of MoS$_2$, which is confirmed by our density functional theory calculations. Furthermore, we show that this spin-layer locking leads to the observation of highly spin-polarized bands in ARPES since these states are almost completely confined within two dimensions. Our findings prove that these highly desired properties of MoS$_2$ can be accessed without thinning it down to the monolayer limit.

preprint2015arXiv

Realization of a vertical topological p-n junction in epitaxial $\mathrm{Sb_2Te_3 / Bi_2Te_3}$ heterostructures

3D topological insulators are a new state of quantum matter which exhibits both a bulk band structure with an insulating energy gap as well as metallic spin-polarized Dirac fermion states when interfaced with a topologically trivial material. There have been various attempts to tune the Dirac point to a desired energetic position for exploring its unusual quantum properties. Here we show a direct experimental proof by angle-resolved photoemission of the realization of a vertical topological p-n junction made of a heterostructure of two different binary 3D TI materials $\mathrm{Bi_2Te_3}$ and $\mathrm{Sb_2Te_3}$ epitaxially grown on Si(111). We demonstrate that the chemical potential is tunable by about 200 meV when decreasing the upper $\mathrm{Sb_2Te_3}$ layer thickness from 25 to 6 quintuple layers without applying any external bias. These results make it realistic to observe the topological exciton condensate and pave the way of exploring other exotic quantum phenomena in the near future.

preprint2015arXiv

Relativistic dynamical spin excitations of magnetic adatoms

We present a first-principles theory of dynamical spin excitations in the presence of spin-orbit coupling. The broken global spin rotational invariance leads to a new sum rule. We explore the competition between the magnetic anisotropy energy and the external magnetic field, as well as the role of electron-hole excitations, through calculations for 3$d$-metal adatoms on the Cu(111) surface. The spin excitation resonance energy and lifetime display non-trivial behavior, establishing the strong impact of relativistic effects. We legitimate the use of the Landau-Lifshitz-Gilbert equation down to the atomic limit, but with parameters that differ from a stationary theory.

preprint2015arXiv

Skew scattering in dilute ferromagnetic alloys

The challenging problem of skew scattering for Hall effects in dilute ferromagnetic alloys, with intertwined effects of spin-orbit coupling, magnetism and impurity scattering, is studied here from first principles. Our main aim is to identify chemical trends and work out simple rules for large skew scattering in terms of the impurity and host states at the Fermi surface, with particular emphasis on the interplay of the spin and anomalous Hall effects in one and the same system. The predicted trends are benchmarked by referring to three different \emph{ab initio} methods based on different approximations with respect to the electronic structure and transport properties.

preprint2014arXiv

Bulk Band Structure of Bi$_2$Te$_3$

The bulk band structure of Bi$_2$Te$_3$ has been determined by angle-resolved photoemission spectroscopy and compared to first-principles calculations. We have performed calculations using the local density approximation (LDA) of density functional theory and the one-shot $GW$ approximation within the all-electron full-potential linearized augmented-plane-wave (FLAPW) formalism, fully taking into account spin-orbit coupling. Quasiparticle effects produce significant changes in the band structure of \bite~when compared to LDA. Experimental and calculated results are compared in the spectral regions where distinct differences between the LDA and $GW$ results are present. Overall a superior agreement with $GW$ is found, highlighting the importance of many-body effects in the band structure of this family of topological insulators.

preprint2014arXiv

Dzyaloshinskii-Moriya interaction and chiral magnetism in 3$d$-5$d$ zig-zag chains: Tight-binding model and ab initio calculations

We investigate the chiral magnetic order in free-standing planar 3$d$-5$d$ bi-atomic metallic chains (3$d$: Fe, Co; 5$d$: Ir, Pt, Au) using first-principles calculations based on density functional theory. We found that the antisymmetric exchange interaction, commonly known as Dzyaloshinskii-Moriya interaction (DMI), contributes significantly to the energetics of the magnetic structure. We used the full-potential linearized augmented plane wave method and performed self-consistent calculations of homogeneous spin spirals, calculating the DMI by treating the effect of spin-orbit interaction (SOI) in the basis of the spin-spiral states in first-order perturbation theory. To gain insight into the DMI results of our ab initio calculations, we develop a minimal tight-binding model of three atoms and 4 orbitals that contains all essential features: the spin-canting between the magnetic $3d$ atoms, the spin-orbit interaction at the $5d$ atoms, and the structure inversion asymmetry facilitated by the triangular geometry. We found that spin-canting can lead to spin-orbit active eigenstates that split in energy due to the spin-orbit interaction at the $5d$ atom. We show that, the sign and strength of the hybridization, the bonding or antibonding character between $d$-orbitals of the magnetic and non-magnetic sites, the bandwidth and the energy difference between states occupied and unoccupied states of different spin projection determine the sign and strength of the DMI. The key features observed in the trimer model are also found in the first-principles results.

preprint2014arXiv

First-principles analysis of a homo-chiral cycloidal magnetic structure in a monolayer Cr on W(110)

The magnetic structure of a Cr monolayer on a W(110) substrate is investigated by means of first-principles calculations based on the noncollinear spin density functional theory (DFT). As magnetic ground state we find a long-period homochiral left-rotating spin spiral on-top of an atomic-scale anti-ferromagnetic order of nearest neighbor atoms. The rotation angle of the magnetic moment changes inhomogeneously from atom to atom across the spiral. We predict a propagation direction along the crystallographic [001] direction with a period length of 14.3 nm, which is in excellent agreement with a modulation of the local anti-ferromagnetic contrast observed in spin-polarized scanning tunneling microscope experiments by Santos et al. [New J. Phys. 10, 013005 (2008)]. We identify the Dzyaloshinskii-Moriya interaction (DMI) as origin of the homochiral magnetic structure, competing with the Heisenberg-type exchange interaction and magneto-crystalline anisotropy energy. From DFT calculations we extract parameters for a micromagnetic model and thereby determine a considerable inhomogeneity of the spin spiral, increasing the period length by 6% compared to homogeneous spin spirals. The results are compared to the behavior of a Mn and Fe monolayer and Fe doublelayer on a W(110) substrate.

preprint2014arXiv

Quantum well states and amplified spin-dependent Friedel oscillations in thin films

Electrons mediate many of the interactions between atoms in a solid. Their propagation in a material determines its thermal, electrical, optical, magnetic and transport properties. Therefore, the constant energy contours characterizing the electrons, in particular the Fermi surface, have a prime impact on the behavior of materials. If anisotropic, the contours induce strong directional dependence at the nanoscale in the Friedel oscillations surrounding impurities. Here we report on giant anisotropic charge density oscillations focused along specific directions with strong spin-filtering after scattering at an oxygen impurity embedded in the surface of a ferromagnetic thin film of Fe grown on W(001). Utilizing density functional theory, we demonstrate that by changing the thickness of the Fe films, we control quantum well states confined to two dimensions that manifest as multiple flat energy contours, impinging and tuning the strength of the induced charge oscillations which allow to detect the oxygen impurity at large distances ($\approx$ 50nm).

preprint2014arXiv

Spin relaxation and spin Hall transport in 5d transition-metal ultrathin films

The spin relaxation induced by the Elliott-Yafet mechanism and the extrinsic spin Hall conductivity due to the skew-scattering are investigated in 5d transition-metal ultrathin films with self-adatom impurities as scatterers. The values of the Elliott-Yafet parameter and of the spin-flip relaxation rate reveal a correlation with each other that is in agreement with the Elliott approximation. At 10-layer thickness, the spin-flip relaxation time in 5d transition-metal films is quantitatively reported about few hundred nanoseconds at atomic percent which is one and two orders of magnitude shorter than that in Au and Cu thin films, respectively. The anisotropy effect of the Elliott-Yafet parameter and of the spin-flip relaxation rate with respect to the direction of the spin-quantization axis in relation to the crystallographic axes is also analyzed. We find that the anisotropy of the spin-flip relaxation rate is enhanced due to the Rashba surface states on the Fermi surface, reaching values as high as 97% in 10-layer Hf(0001) film or 71% in 10-layer W(110) film. Finally, the spin Hall conductivity as well as the spin Hall angle due to the skew-scattering off self-adatom impurities are calculated using the Boltzmann approach. Our calculations employ a relativistic version of the first-principles full-potential Korringa-Kohn-Rostoker Green function method.

preprint2014arXiv

Spin Spirals in Surface Alloys on Ru(0001): A First-principles Study

We have used ab initio density functional theory to compute the magnetic ground states of the surface alloy systems FeAu$_2$/Ru(0001) and MnAu$_2$/Ru(0001). For both systems, we find that the lowest energy magnetic configuration corresponds to a left-rotating spin spiral, in which the sense of rotation is determined by the Dzyaloshinskii-Moriya interaction. These spirals are lower in energy than the ferromagnetic configuration by 3-4 meV per nm$^2$. We also find that FeAu$_2$/Ru(0001) has a significantly high magnetic anisotropy energy, of the order 1 meV per Fe atom. By comparing with the corresponding freestanding alloy monolayers, we find that the presence of the Ru substrate plays a significant role in determining the magnetic properties of the surface alloy systems.

preprint2014arXiv

Spin transport and spin-caloric effects in (Cr,Zn)Te half-metallic nanostructures: Effect of spin disorder at elevated temperatures from first principles

An important contribution to the thermoelectric and spin-caloric transport properties in magnetic materials at elevated temperatures is the formation of a spin-disordered state due to local moment fluctuations. This effect has not been largely investigated so far. We focus on various magnetic nanostructures of CrTe in the form of thin layers or nanowires embedded in ZnTe matrix, motivated by the miniaturization of spintronics devices and by recent suggestions that magnetic nanostructures can lead to extraordinary thermoelectric effects due to quantum confinement. The electronic structure of the studied systems is calculated within the multiple scattering screened Korringa-Kohn-Rostoker Green function (KKR-GF) framework. The Monte Carlo method is used to simulate the magnetization in the temperature induced spin disorder. The transport properties are evaluated from the transmission probability obtained using the Baranger-Stone approach within the KKR-GF framework. We find qualitative and quantitative changes in the thermoelectric and spin-caloric coefficients when spin-disorder is included in the calculation. Furthermore, we show that substitutional impurities in CrTe nanowires could considerably enhance the Seebeck coefficient and the thermoelectric figure of merit.

preprint2014arXiv

Spin-caloric transport properties of cobalt nanostructures: spin disorder effects from first principles

The fundamental aspects of spin-dependent transport processes and their interplay with temperature gradients, as given by the spin Seebeck coefficient, are still largely unexplored and a multitude of contributing factors must be considered. We used density functional theory together with a Monte-Carlo-based statistical method to simulate simple nanostructures, such as Co nanowires and films embedded in a Cu host or in vacuum, and investigated the influence of spin-disorder scattering on electron transport at elevated temperatures. While we show that the spin-dependent scattering of electrons due to temperature induced disorder of the local magnetic moments contributes significantly to the resistance, thermoelectric and spin-caloric transport coefficients, we also conclude that the actual magnitude of these effects cannot be predicted, quantitatively or qualitatively, without such detailed calculations.

preprint2014arXiv

Spin-flip hot spots in ultrathin films of monovalent metals: Enhancement and anisotropy of the Elliott-Yafet parameter

In contrast to the long-known fact that spin-flip hot spots, i.e., special \vc{k}-points on the Fermi surface showing a high spin-mixing parameter, do not occur in the bulk of monovalent (noble and alkali) metals, we found them on the surface Brillouin-zone boundary of ultrathin films of these metals. Density-functional calculations within the Korringa-Kohn-Rostoker Green function method for ultrathin (001) oriented Cu, Ag, and Au films of 10-layer thickness show that the region around the hot spots can have a substantial contribution, e.g.\ 52\% in Au(001), to the integrated spin-mixing parameter, that could lead to a significant enhancement of the spin-relaxation rate or spin-Hall angle in thin films. Owing to the appearance of spin-flip hot-spots, a large anisotropy of the Elliott-Yafet parameter [50\% for Au(001)] is also found in these systems. The findings are important for spintronics applications in which noble-metals are frequently used and in which the dimensionality of the sample is reduced.

preprint2014arXiv

Spin-orbit torques in Co/Pt(111) and Mn/W(001) magnetic bilayers from first principles

An applied electric current through a space-inversion asymmetric magnet induces spin-orbit torques (SOTs) on the magnetic moments, which holds much promise for future memory devices. We discuss general Green's function expressions suitable to compute the linear-response SOT in disordered ferromagnets. The SOT can be decomposed into an even and an odd component with respect to magnetization reversal, where in the limit of vanishing disorder the even SOT is given by the constant Berry curvature of the occupied states, while the odd part exhibits a divergence with respect to disorder strength. Within this formalism, we perform first principles density-functional theory calculations of the SOT in Co/Pt(111) and Mn/W(001) magnetic bilayers. We find the even and odd torque components to be of comparable magnitude. Moreover, the odd torque depends strongly on an additional capping layer, while the even torque is less sensitive. We show that the even torque is nearly entirely mediated by spin currents in contrast to the odd torque, which can contain an important contribution not due to spin transfer. Our results are in agreement with experiments, showing that our linear response theory is well-suited for the description of SOTs in complex ferromagnets.

preprint2014arXiv

Spin-orbit torques in L1$_0$-FePt/Pt thin films driven by electrical and thermal currents

Using the linear response formalism for the spin-orbit torque (SOT) we compute from first principles the SOT in a system of two layers of L1$_0$-FePt(001) deposited on an fcc Pt(001) substrate of varying thickness. We find that at room temperature the values of the SOTs that are even and odd with respect to magnetization generally lie in the range of values measured and computed for Co/Pt bilayers. We also observe that the even SOT is much more robust with respect to changing the number of layers in the substrate, and as a function of energy it follows the general trend of the even SOT exerted by the spin Hall current in fcc Pt. The odd torque, on the other hand, is strongly affected by modification of the electronic structure for a specific energy window in the limit of very thin films. Moreover, taking the system at hand as an example, we compute the values of the thermal spin-orbit torque (T-SOT). We predict that the gradients of temperature which can be experimentally created in this type of systems will cause a detectable torque on the magnetization. We also underline the correlation between the even T-SOT and the spin Nernst effect, thus motivating a more intensive experimental effort aimed at observation of both phenomena.

preprint2013arXiv

Ab-initio calculation of the effective on-site Coulomb interaction parameters for half-metallic magnets

Correlation effects play an important role in the electronic structure of half-metallic (HM) magnets. In particular, they give rise to non-quasiparticle states above (or below) the Fermi energy at finite temperatures that reduce the spin polarization and, as a consequence, the efficiency of spintronics devices. Employing the constrained random-phase approximation (cRPA) within the full-potential linearized augmented-plane-wave (FLAPW) method using maximally localized Wannier functions, we calculate the strength of the effective on-site Coulomb interaction (Hubbard $U$ and Hund exchange $J$) between localized electrons in different classes of HM magnets considering: (i) \emph{sp}-electron ferromagnets in rock-salt structure, (ii) zincblende 3\emph{d} binary ferromagnets, as well as (iii) ferromagnetic and ferrimagnetic semi- and full-Heusler compounds.

preprint2013arXiv

Anomalous Hall effect in ferromagnets with Gaussian disorder

Using the Kubo formalism we derived expressions and implemented the method for calculating the anomalous Hall conductivity (AHC) in ferromagnets with short-range Gaussian disorder directly from first-principles electronic structure of the perfect crystal. We used this method to calculate the AHC in bcc Fe, fcc Co, L1o-FePd, L1o-FePt as well as thin bcc Fe(001) films. Within our approach we can transparently decompose the conductivity into intrinsic (IC), side jump (SJ) and intrinsic skew-scattering (ISK) contributions. The existence of ISK, which originates from asymmetric Mott scattering but is clearly distinguishable from conventional skew-scattering in that it converges to a finite value in clean limit, was pointed out by Sinitsyn et al. [Phys. Rev. B 75, 045315 (2007)]. Here, we collect all contributions to the AHC in ferromagnets which result in "scattering-independent" AHE in clean limit, and analyze their relative magnitude from first principles calculations. By comparing our results to existing experiments we show that the Gaussian disorder is well suited to model various types of disorder present in real materials, to some extent including the effect of temperature. In particular, we show that in addition to intrinsic and side-jump AHE, the intrinsic skew-scattering can be a major player in determining the magnitude of the AHE in ferromagnets.

preprint2013arXiv

Berry phase theory of Dzyaloshinskii-Moriya interaction and spin-orbit torques

Recent experiments on current-induced domain wall motion in chiral magnets suggest important contributions both from spin-orbit torques (SOTs) and from the Dzyaloshinskii-Moriya interaction (DMI). We derive a Berry phase expression for the DMI and show that within this Berry phase theory DMI and SOTs are intimately related, in a way formally analogous to the relation between orbital magnetization (OM) and anomalous Hall effect (AHE). We introduce the concept of the \textit{twist torque moment}, which probes the internal twist of wave packets in chiral magnets in a similar way like the orbital moment probes the wave packet's internal self rotation. We propose to interpret the Berry phase theory of DMI as a theory of \textit{spiralization} in analogy to the modern theory of OM. We show that the twist torque moment and the spiralization together give rise to a Berry phase governing the response of the SOT to thermal gradients, in analogy to the intrinsic anomalous Nernst effect. The Berry phase theory of DMI is computationally very efficient because it only needs the electronic structure of the collinear magnetic system as input. As an application of the formalism we compute the DMI in Pt/Co, Pt/Co/O and Pt/Co/Al magnetic trilayers and show that the DMI is highly anisotropic in these systems.

preprint2013arXiv

Elimination of the linearization error and improved basis-set convergence within the FLAPW method

We analyze in detail the error that arises from the linearization in linearized augmented-plane-wave (LAPW) basis functions around predetermined energies $E_l$ and show that it can lead to undesirable dependences of the calculated results on method-inherent parameters such as energy parameters $E_l$ and muffin-tin sphere radii. To overcome these dependences, we evaluate approaches that eliminate the linearization error systematically by adding local orbitals (LOs) to the basis set. We consider two kinds of LOs: (i) constructed from solutions $u_l(r,E)$ to the scalar-relativistic approximation of the radial Dirac equation with $E>E_l$ and (ii) constructed from second energy derivatives $\partial^2 u_l(r,E) / \partial E^2$ at $E=E_l$. We find that the latter eliminates the error most efficiently and yields the density functional answer to many electronic and materials properties with very high precision. Finally, we demonstrate that the so constructed LAPW+LO basis shows a more favorable convergence behavior than the conventional LAPW basis due to a better decoupling of muffin-tin and interstitial regions, similarly to the related APW+lo approach, which requires an extra set of LOs to reach the same total energy, though.

preprint2013arXiv

Engineering quantum anomalous Hall phases with orbital and spin degrees of freedom

Combining tight-binding models and first principles calculations, we investigate the quantum anomalous Hall (QAH) effect induced by intrinsic spin-orbit coupling (SOC) in buckled honeycomb lattice with sp orbitals in an external exchange field. Detailed analysis reveals that nontrivial topological properties can arise utilizing not only spin but also orbital degrees of freedom in the strong SOC limit, when the bands acquire non-zero Chern numbers upon undergoing the so-called orbital purification. As a prototype of a buckled honeycomb lattice with strong SOC we choose the Bi(111) bilayer, analyzing its topological properties in detail. In particular, we show the emergence of several QAH phases upon spin exchange of the Chern numbers as a function of SOC strength and magnitude of the exchange field. Interestingly, we observe that in one of such phases, namely, in the quantum spin Chern insulator phase, the quantized charge and spin Hall conductivities co-exist. We consider the possibility of tuning the SOC strength in Bi bilayer via alloying with isoelectronic Sb, and speculate that exotic properties could be expected in such an alloyed system owing to the competition of the topological properties of its constituents. Finally, we demonstrate that 3d dopants can be used to induce a sizeable exchange field in Bi(111) bilayer, resulting in non-trivial Chern insulator properties.

preprint2013arXiv

Exchange interactions and local moment fluctuation corrections in ferromagnets at finite temperatures based on non-collinear density-functional calculations

We explore the derivation of interatomic exchange interactions in ferromagnets within density-functional theory (DFT) and the mapping of DFT results onto a spin Hamiltonian. We delve into the problem of systems comprising atoms with strong spontaneous moments together with atoms with weak induced moments. All moments are considered as degrees of freedom, with the strong moments thermally fluctuating only in angle and the weak moments thermally fluctuating in angle and magnitude. We argue that a quadratic dependence of the energy on the weak local moments magnitude, which is a good approximation in many cases, allows for an elimination of the weak-moment degrees of freedom from the thermodynamic expressions in favor of a renormalization of the Heisenberg interactions among the strong moments. We show that the renormalization is valid at all temperatures accounting for the thermal fluctuations and resulting in temperature-independent renormalized interactions. These are shown to be the ones directly derived from total-energy DFT calculations by constraining the strong-moment directions, as is done e.g. in spin-spiral methods. We furthermore prove that within this framework the thermodynamics of the weak-moment subsystem, and in particular all correlation functions, can be derived as polynomials of the correlation functions of the strong-moment subsystem with coefficients that depend on the spin susceptibility and that can be calculated within DFT. These conclusions are rigorous under certain physical assumptions on the measure in the magnetic phase space. We implement the scheme in the full-potential linearized augmented plane wave method using the concept of spin-spiral states, considering applicable symmetry relations and the use of the magnetic force theorem. Our analytical results are corroborated by numerical calculations employing DFT and a Monte Carlo method.

preprint2013arXiv

Exchange parameters and adiabatic magnon energies from spin-spiral calculations

We present a method of extracting the exchange parameters of the classical Heisenberg model from first-principles calculations of spin-spiral total energies based on density functional theory. The exchange parameters of the transition-metal monoxides MnO and NiO are calculated and used to estimate magnetic properties such as transition temperatures and magnon energies. Furthermore we show how to relate the magnon energies directly to differences in spin-spiral total energies for systems containing an arbitrary number of magnetic sublattices. This provides a comparison between magnon energies using a finite number of exchange parameters and the infinite limit.

preprint2013arXiv

Hybridisation at the organic-metal interface: a surface-scientific analogue of Hückel's rule?

We demonstrate that cyclooctatetraene (COT) can be stabilised in different conformations when adsorbed on different noble-metal surfaces due to varying molecule-substrate interaction. While at first glance the behaviour seems to be in accordance with Hückel's rule, a theoretical analysis reveals no significant charge transfer. The driving mechanism for the conformational change is hybridisation at the organic-metal interface and does not necessitate any charge transfer.

preprint2013arXiv

Magnetism- and impurity-assisted chain creation in Ir and Pt break junctions

Applying the generalization of the model for chain formation in break-junctions [JPCM 24, 135501 (2012)], we study the effect of light impurities on the energetics and elongation properties of Pt and Ir chains. Our model enables us with a tool ideal for detailed analysis of impurity assisted chain formation, where zigzag bonds play an important role. In particular we focus on H (s-like) and O (p-like) impurities and assume, for simplicity, that the presence of impurity atoms in experiments results in ..M-X-M-X-... (M: metal, X: impurity) chain structure in between the metallic leads. Feeding our model with material-specific parameters from systematic full-potential first-principles calculations, we find that the presence of such impurities strongly affects the binding properties of the chains. We find that while both types of impurities enhance the probability of chains to be elongated, the s-like impurities lower the chain's stability. We also analyze the effect of magnetism and spin-orbit interaction on the growth properties of the chains.

preprint2013arXiv

Scattering-Independent Anomalous Nernst Effect in Ferromagnets

Using the full-potential linearized augmented plane-wave method within the density functional theory, we compute all contributions to the scattering independent part of the thermoelectric conductivity tensor, namely the intrinsic contribution and the side-jump contribution. For the ferromagnetic materials bcc Fe, hcp Co, fcc Ni and L1_0 ordered alloys FePd and FePt, our investigations of the energy and temperature dependence of the intrinsic and side-jump contributions show that they are both of equal importance. Overall, our calculations are able to correctly reproduce the order of magnitude and sign of the experimentally measured signal, suggesting that the scattering independent part plays an important role in the anomalous Nernst effect of ferromagnets.

preprint2013arXiv

Spin relaxation and the Elliott-Yafet parameter in W(001) ultrathin films: surface states, anisotropy and oscillation effects

Using first-principles methods based on density-functional theory we investigate the spin relaxation in W(001) ultrathin films. Within the framework of the Elliott-Yafet theory we calculate the spin mixing of the Bloch states and we explicitly consider spin-flip scattering off self-adatoms. We find an oscillatory behavior of the spin-mixing parameter and relaxation rate as a function of the film thickness, which we trace back tosurface-state properties. We also analyze the Rashba effect experienced by the surface states and discuss its influence on the spin relaxation. Finally we calculate the anisotropy of the spin-relaxation rate with respect to the polarization direction of the excited spin population relative to the crystallographic axes of thefilm. We find that the spin-relaxation rate can increase by as much as 47% when the spin polarization is directed out of plane, compared to the case when it is in plane. Our calculations are based on the multiple-scattering formalism of the Korringa-Kohn-Rostoker Green-function method.

preprint2013arXiv

Strong magnon softening in tetragonal FeCo compounds

Magnons play an important role in fast precessional magnetization reversal processes serving as a heat bath for dissipation of the Zeeman energy and thus being responsible for the relaxation of magnetization. Employing \emph{ab initio} many-body perturbation theory we studied the magnon spectra of the tetragonal FeCo compounds considering three different experimental $c/a$ ratios, $c/a=$1.13, 1.18, and 1.24 corresponding to FeCo grown on Pd, Ir, and Rh, respectively. We find that for all three cases the short-wave-length magnons are strongly damped and tetragonal distortion gives rise to a significant magnon softening. The magnon stiffness constant $D$ decreases almost by a factor of two from FeCo/Pd to FeCo/Rh. The combination of soft magnons together with the giant magnetic anisotropy energy suggests FeCo/Rh to be a promising material for perpendicular magnetic recording applications.

preprint2013arXiv

Symmetry and magnitude of spin-orbit torques in ferromagnetic heterostructures

Current-induced spin torques are of great interest to manipulate the orientation of nanomagnets without applying external magnetic fields. They find direct application in non-volatile data storage and logic devices, and provide insight into fundamental processes related to the interdependence between charge and spin transport. Recent demonstrations of magnetization switching induced by in-plane current injection in ferromagnetic heterostructures have drawn attention to a class of spin torques based on orbital-to-spin momentum transfer, which is alternative to pure spin transfer torque (STT) between noncollinear magnetic layers and amenable to more diversified device functions. Due to the limited number of studies, however, there is still no consensus on the symmetry, magnitude, and origin of spin-orbit torques (SOTs). Here we report on the quantitative vector measurement of SOTs in Pt/Co/AlO trilayers using harmonic analysis of the anomalous and planar Hall effects as a function of the applied current and magnetization direction. We provide an all-purpose scheme to measure the amplitude and direction of SOTs for any arbitrary orientation of the magnetization, including corrections due to the interplay of Hall and thermoelectric effects. Based on general space and time inversion symmetry arguments, we show that asymmetric heterostructures allow for two different SOTs having odd and even behavior with respect to magnetization reversal. Our results reveal a scenario that goes beyond established models of the Rashba and spin Hall contributions to SOTs. The even SOT is STT-like but stronger than expected from the spin Hall effect in Pt. The odd SOT is composed of a constant field-like term and an additional component, which is strongly anisotropic and does not correspond to a simple Rashba field.

preprint2012arXiv

Anisotropy of spin relaxation and transverse transport in metals

Using first principles methods we explore the anisotropy of the spin relaxation and transverse transport properties in bulk metals with respect to the direction of the spin quantization axis in paramagnets or of the spontaneous magnetization in ferromagnets. Owing to the presence of the spin-orbit interaction the orbital and spin character of the Bloch states depends sensitively on the orientation of the spins relative to the crystal axes. This leads to drastic changes in quantities which rely on interband mixing induced by the spin-orbit interaction. The anisotropy is particularly striking for quantities which exhibit spiky and irregular distribution in the Brillouin zone, such as the spin-mixing parameter or the Berry curvature of the electronic states. We demonstrate this for three cases: (i) the Elliott-Yafet spin-relaxation mechanism in paramagnets with structural inversion symmetry; (ii) the intrinsic anomalous Hall effect in ferromagnets; and (iii) the spin Hall effect in paramagnets. We discuss the consequences of the pronounced anisotropic behavior displayed by these properties for spin-polarized transport applications.

preprint2012arXiv

Anisotropy of spin relaxation in metals

The concept of anisotropy of spin relaxation in non-magnetic metals with respect to the spin direction of the injected electrons relative to the crystal orientation is introduced. The effect is related to an anisotropy of the Elliott-Yafet parameter, arising from a modulation of the decomposition of the spin-orbit Hamiltonian into spin-conserving and spin-flip terms as the spin quantization axis is varied. This anisotropy, reaching gigantic values for uniaxial transition-metals (e.g. 830% for hcp Hf) as density-functional calculations show, is related to extended "spin-flip hot areas" on the Fermi surface created by the proximity of extended sheets of the surface, or "spin-flip hot loops" at the Brillouin zone boundary, and has no theoretical upper limit. Possible ways of measuring the effect as well as consequences in application are briefly outlined.

preprint2012arXiv

Correlations in sequences of generalized eigenproblems arising in Density Functional Theory

Density Functional Theory (DFT) is one of the most used ab initio theoretical frameworks in materials science. It derives the ground state properties of a multi-atomic ensemble directly from the computation of its one-particle density \nr .In DFT-based simulations the solution is calculated through a chain of successive self-consistent cycles; in each cycle a series of coupled equations (Kohn-Sham) translates to a large number of generalized eigenvalue problems whose eigenpairs are the principal means for expressing \nr. A simulation ends when \nr\ has converged to the solution within the required numerical accuracy. This usually happens after several cycles, resulting in a process calling for the solution of many sequences of eigenproblems. In this paper, the authors report evidence showing unexpected correlations between adjacent eigenproblems within each sequence. By investigating the numerical properties of the sequences of generalized eigenproblems it is shown that the eigenvectors undergo an "evolution" process. At the same time it is shown that the Hamiltonian matrices exhibit a similar evolution and manifest a specific pattern in the information they carry. Correlation between eigenproblems within a sequence is of capital importance: information extracted from the simulation at one step of the sequence could be used to compute the solution at the next step. Although they are not explored in this work, the implications could be manifold: from increasing the performance of material simulations, to the development of an improved iterative solver, to modifying the mathematical foundations of the DFT computational paradigm in use, thus opening the way to the investigation of new materials.

preprint2012arXiv

Electronic properties of graphene nanoribbons under gate electric fields

Quantum-dot states in graphene nanoribbons (GNR) were calculated using density-functional theory, considering the effect of the electric field of gate electrodes. The field is parallel to the GNR plane and was generated by an inhomogeneous charge sheet placed atop the ribbon. Varying the electric field allowed to observe the development of the GNR states and the formation of localized, quantum-dot-like states in the band gap. The calculation has been performed for armchair GNRs and for armchair ribbons with a zigzag section. For the armchair GNR a static dielectric constant of ε approx. 4 could be determined.

preprint2012arXiv

GW study of the half-metallic Heusler compounds Co2MnSi and Co2FeSi

Quasiparticle spectra of potentially half-metallic Co2MnSi and Co2FeSi Heusler compounds have been calculated within the one-shot GW approximation in an all-electron framework without adjustable parameters. For Co2FeSi the many-body corrections are crucial: a pseudogap opens and good agreement of the magnetic moment with experiment is obtained. Otherwise, however, the changes with respect to the density-functional-theory starting point are moderate. For both cases we find that photoemission and x-ray absorption spectra are well described by the calculations. By comparison with the GW density of states, we conclude that the Kohn-Sham eigenvalue spectrum provides a reasonable approximation for the quasiparticle spectrum of the Heusler compounds considered in this work.

preprint2012arXiv

Information transfer by vector spin chirality in finite magnetic chains

Vector spin chirality is one of the fundamental characteristics of complex magnets. For a one-dimensional spin-spiral state it can be interpreted as the handedness, or rotational sense of the spiral. Here, using spin-polarized scanning tunneling microscopy, we demonstrate the occurrence of an atomic-scale spin-spiral in finite individual bi-atomic Fe chains on the (5x1)-Ir(001) surface. We show that the broken inversion symmetry at the surface promotes one direction of the vector spin chirality, leading to a unique rotational sense of the spiral in all chains. Correspondingly, changes in the spin direction of one chain end can be probed tens of nanometers away, suggesting a new way of transmitting information about the state of magnetic objects on the nanoscale.

preprint2012arXiv

Magnetic adatom induced skyrmion-like spin texture in surface electron waves

When a foreign atom is placed on a surface of a metal, the surrounding sea of electrons responds screening the additional charge leading to oscillations or ripples. On surfaces, those electrons are sometimes confined to two-dimensional surface states, whose spin-degeneracy is lifted due to the Rashba effect arising from the spin-orbit interaction of electrons and the inversion asymmetric environment. It is believed that at least for a single adatom scanning tunneling microscopy measurements are insensitive to the Rashba splitting i.e. no signatures in the charge oscillations will be observed. Resting on scattering theory, we demonstrate that, if magnetic, one single adatom is enough to visualize the presence of the Rashba effect in terms of an induced spin-magnetization of the surrounding electrons exhibiting a twisted spin texture described as superposition of two skyrmionic waves of opposite chirality.

preprint2012arXiv

Parity effect in ground state localization of antiferromagnetic chains coupled to a ferromagnet

We investigate the ground states of antiferromagnetic Mn nanochains on Ni(110) by spin-polarized scanning tunneling microscopy in combination with theory. While the ferrimagnetic linear trimer experimentally shows the predicted collinear classical ground state, no magnetic contrast was observed for dimers and tetramers where non-collinear structures were expected based on ab-initio theory. This striking observation can be explained by zero-point energy motion for even numbered chains derived within a classical equation of motion leading to non classical ground states. Thus, depending on the parity of the chain length, the system shows a classical or a quantum behavior.

preprint2012arXiv

Precise response functions in all-electron methods: Application to the optimized-effective-potential approach

The optimized-effective-potential (OEP) method is a special technique to construct local Kohn-Sham potentials from general orbital-dependent energy functionals. In a recent publication [M. Betzinger, C. Friedrich, S. Blügel, A. Görling, Phys. Rev. B 83, 045105 (2011)] we showed that uneconomically large basis sets were required to obtain a smooth local potential without spurious oscillations within the full-potential linearized augmented-plane-wave method (FLAPW). This could be attributed to the slow convergence behavior of the density response function. In this paper, we derive an incomplete-basis-set correction for the response, which consists of two terms: (1) a correction that is formally similar to the Pulay correction in atomic-force calculations and (2) a numerically more important basis response term originating from the potential dependence of the basis functions. The basis response term is constructed from the solutions of radial Sternheimer equations in the muffin-tin spheres. With these corrections the local potential converges at much smaller basis sets, at much fewer states, and its construction becomes numerically very stable. We analyze the improvements for rock-salt ScN and report results for BN, AlN, and GaN, as well as the perovskites CaTiO3, SrTiO3, and BaTiO3. The incomplete-basis-set correction can be applied to other electronic-structure methods with potential-dependent basis sets and opens the perspective to investigate a broad spectrum of problems in theoretical solid-state physics that involve response functions.

preprint2012arXiv

Strength of the Effective Coulomb Interaction at Metal and Insulator Surfaces

The effective on-site Coulomb interaction (Hubbard $U$) between localized electrons at crystal surfaces is expected to be enhanced due to the reduced coordination number and reduced subsequent screening. By means of first principles calculations employing the constrained random-phase approximation (cRPA) we show that this is indeed the case for simple metals and insulators but not necessarily for transition metals and insulators that exhibit pronounced surface states. In the latter case, the screening contribution from surface states as well as the influence of the band narrowing increases the electron polarization to such an extent as to overcompensate the decrease resulting from the reduced effective screening volume. The Hubbard $U$ parameter is thus substantially reduced in some cases, e.g., by around 30% for the (100) surface of bcc Cr.

preprint2012arXiv

Topological phases of Bi(111) bilayer in an external exchange field

Using first principles methods, we investigate topological phase transitions as a function of exchange field in a Bi(111) bilayer. Evaluation of the spin Chern number for different magnitudes of the exchange field reveals that when the time reversal symmetry is broken by a small exchange field, the system enters the time-reversal broken topological insulator phase, introduced by Yang {\it et al.} in Phys. Rev. Lett. 107, 066602 (2011). After a metallic phase in the intermediate region, the quantum anomalous Hall phase with non-zero Chern number emerges at a sufficiently large exchange field. We analyze the phase diagram from the viewpoint of the evolution of the electronic structure, edge states and transport properties, and demonstrate that different topological phases can be distinguished by the spin-polarization of the edge states as well as spin or charge transverse conductivity.

preprint2011arXiv

Ab Initio Theory of Scattering-Independent Anomalous Hall Effect

We report on first-principles calculations of the side-jump contribution to the anomalous Hall conductivity (AHC) directly from the electronic structure of a perfect crystal. We implemented our approach for a short-range scattering disorder model within the density functional theory and computed the full scattering-independent AHC in elemental bcc Fe, hcp Co, fcc Ni, and L1o FePd and FePt alloys. The full AHC thus calculated agrees systematically with experiment to a degree unattainable so far, correctly capturing the previously missing elements of side-jump contributions, hence paving the way to a truly predictive theory of the anomalous Hall effect and turning it from a characterization tool to a probing tool of multi-band complex electronic band structures.

preprint2011arXiv

Anisotropic intrinsic anomalous Hall effect in ordered 3dPt alloys

By performing first principles calculations we investigate the intrinsic anomalous Hall conductivity (AHC) and its anisotropy in ordered L1o FePt, CoPt and NiPt ferromagnets, and their intermediate alloys. We demonstrate that the AHC in this family of compounds depends strongly on the direction of the magnetization in the crystal. We predict that such pronounced orientational dependence in combination with the general decreasing trend of the AHC when going from FePt to NiPt leads to a sign change of the AHC upon rotating the magnetization direction in the crystal of CoPt alloy. We also suggest that for a range of concentration x in Co(x)Ni(1-x)Pt alloy it is possible to achieve a complete quenching of the anomalous Hall current for a certain direction of the magnetization in the crystal. By analyzing the spin-resolved AHC in 3dPt alloys we endeavor to relate the overall trend of the AHC in these compounds to the changes in their densities of d-states around the Fermi energy upon varying the atomic number. Moreover, we show the generality of the phenomenon of anisotropic anomalous Hall effect by demonstrating its occurrence within the three-band tight-binding model.

preprint2011arXiv

Effective Coulomb interaction in transition metals from constrained random-phase approximation

The effective on-site Coulomb interaction (Hubbard $U$) between localized \textit{d} electrons in 3\textit{d}, 4\textit{d}, and 5\textit{d} transition metals is calculated employing a new parameter-free realization of the constrained random-phase approximation using Wannier functions within the full-potential linearized augmented-plane-wave method. The $U$ values lie between 1.5 and 5.7 eV and depend on the crystal structure, spin polarization, \textit{d} electron number, and \textit{d} orbital filling. On the basis of the calculated $U$ parameters, we discuss the strength of the electronic correlations and the instability of the paramagnetic state towards the ferromagnetic one for 3\textit{d} metals.

preprint2011arXiv

Electrically tunable quantum anomalous Hall effect in 5d transition-metal adatoms on graphene

The combination of the unique properties of graphene with spin polarization and magnetism for the design of new spintronic concepts and devices has been hampered by the small Coulomb interaction and the tiny spin-orbit coupling of carbon in pristine graphene. Such device concepts would take advantage of the control of the spin degree of freedom utilizing the widely available electric fields in electronics or of topological transport mechanisms such as the conjectured quantum anomalous Hall effect. Here we show, using first-principles methods, that 5d transition-metal (TM) adatoms deposited on graphene display remarkable magnetic properties. All considered TM adatoms possess significant spin moments with colossal magnetocrystalline anisotropy energies as large as 50 meV per TM atom. We reveal that the magneto-electric response of deposited TM atoms is extremely strong and in some cases offers even the possibility to switch the spontaneous magnetization direction by a moderate external electric field. We predict that an electrically tunable quantum anomalous Hall effect can be observed in this type of hybrid materials.

preprint2011arXiv

Ferromagnetic spin coupling of 2p-impurities in band insulators stabilized by intersite Coulomb interaction: Nitrogen-doped MgO

For a nitrogen dimer in insulating MgO, a ferromagnetic coupling between spin-polarized $2p$-holes is revealed by calculations based on the density functional theory amended by an on-site Coulomb interaction and corroborated by the Hubbard model. It is shown that the ferromagnetic coupling is facilitated by a T-shaped orbital arrangement of the $2p$-holes, which is in its turn controlled by an intersite Coulomb interaction due to the directionality of the $p$-orbitals. We thus conjecture that this interaction is an important ingredient of ferromagnetism in band insulators with $2p$ dopants.

preprint2011arXiv

First-principles calculation of electronic excitations in solids with SPEX

We describe the software package SPEX, which allows first-principles calculations of quasiparticle and collective electronic excitations in solids using techniques from many-body perturbation theory. The implementation is based on the full-potential linearized augmented-plane-wave (FLAPW) method, which treats core and valence electrons on an equal footing and can be applied to a wide range of materials, including transition metals and rare earths. After a discussion of essential features that contribute to the high numerical efficiency of the code, we present illustrative results for quasiparticle band structures calculated within the GW approximation for the electronic self-energy, electron-energy-loss spectra with inter- and intraband transitions as well as local-field effects, and spin-wave spectra of itinerant ferromagnets. In all cases the inclusion of many-body correlation terms leads to very good quantitative agreement with experimental spectroscopies.

preprint2011arXiv

Graphene on Ir(111): Physisorption with chemical modulation

The non-local van der Waals density functional (vdW-DF) approach is applied to calculate the binding of graphene to Ir(111). The precise agreement of the calculated mean height h = 3.41 Å of the C atoms with their mean height h = (3.38 \pm 0.04) Å as measured by the X-ray standing wave (XSW) technique provides a benchmark for the applicability of the non-local functional. We find bonding of graphene to Ir(111) to be due to the van der Waals interaction with an antibonding average contribution from chemical interaction. Despite its globally repulsive character, in certain areas of the large graphene moiré unit cell charge accumulation between Ir substrate and grapheneC atoms is observed, signaling a weak covalent bond formation.

preprint2011arXiv

Local exact exchange potentials within the all-electron FLAPW method and a comparison with pseudopotential results

We present a general numerical approach to construct local Kohn-Sham potentials from orbital-dependent functionals within the all-electron full-potential linearized augmented-plane-wave (FLAPW) method, in which core and valence electrons are treated on an equal footing. As a practical example, we present a treatment of the orbital-dependent exact-exchange (EXX) energy and potential. A formulation in terms of a mixed product basis, which is constructed from products of LAPW basis functions, enables a solution of the optimized-effective-potential (OEP) equation with standard numerical algebraic tools and without shape approximations for the resulting potential. We find that the mixed product and LAPW basis sets must be properly balanced to obtain smooth and converged EXX potentials without spurious oscillations. The construction and convergence of the exchange potential is analyzed in detail for diamond. Our all-electron results for C, Si, SiC, Ge, GaAs semiconductors as well as Ne and Ar noble-gas solids are in very favorable agreement with plane-wave pseudopotential calculations. This confirms the adequacy of the pseudopotential approximation in the context of the EXX-OEP formalism and clarifies a previous contradiction between FLAPW and pseudopotential results.

preprint2011arXiv

Modelling impurity-assisted chain creation in noble-metal break junctions

In this work we present the generalization of the model for chain formation in break-junctions, introduced by Thiess {\it et al.} [\onlinecite{Alex08}], to zigzag transition-metal chains with $s$ and $p$ impurities. We apply this extended model to study the producibility trends for noble-metal chains with impurities, often present in break junction experiments, namely, Cu, Ag and Au chains with H, C, O and N adatoms. Providing the material-specific parameters for our model from systematic full-potential linearized augmented plane-wave first-principles calculations, we find that the presence of such impurities crucially affects the binding properties of the noble-metal chains. We reveal that both, the impurity-induced bond strengthening and the formation of zigzag bonds, can lead to a significantly enhanced probability for chain formation in break junctions.

preprint2011arXiv

The HSE hybrid functional within the FLAPW method and its application to GdN

We present an implementation of the Heyd-Scuseria-Ernzerhof (HSE) hybrid functional within the full-potential linearized augmented-plane-wave (FLAPW) method. Pivotal to the HSE functional is the screened electron-electron interaction, which we separate into the bare Coulomb interaction and the remainder, a slowly varying function in real space. Both terms give rise to exchange potentials, which sum up to the screened nonlocal exchange potential of HSE. We evaluate the former with the help of an auxiliary basis, defined in such a way that the bare Coulomb matrix becomes sparse. The latter is computed in reciprocal space, exploiting its fast convergence behavior in reciprocal space. This approach is general and can be applied to a whole class of screened hybrid functionals. We obtain excellent agreement of band gaps and lattice constants for prototypical semiconductors and insulators with electronic-structure calculations using plane-wave or Gaussian basis sets. We apply the HSE hybrid functional to examine the ground-state properties of rocksalt GdN, which have been controversially discussed in literature. Our results indicate that there is a half-metal to insulator transition occurring between the theoretically optimized lattice constant at 0 K and the experimental lattice constant at room temperature. Overall, we attain good agreement with experimental data for band transitions, magnetic moments, and the Curie temperature.

preprint2011arXiv

The role of spin-flip transitions in the anomalous Hall effect of FePt alloy

We carry out first-principles calculations which demonstrate the importance of the non-spin-conserving part of the spin-orbit interaction for the intrinsic anomalous Hall conductivity of ordered FePt alloys. The impact of this interaction is strongly reduced if Pt is replaced by the lighter isoelectronic element Pd. An analysis of the interband transitions responsible for the anomalous velocity reveals that spin-flip transitions occur not only at avoided band crossings near the Fermi level, but also between well-separated pairs of bands with similar dispersions. We also predict a strong anisotropy in the anomalous Hall conductivity of FePt caused entirely by low-frequency spin-flip transitions.

preprint2010arXiv

Anisotropic spin Hall effect from first principles

We report on first principles calculations of the anisotropy of the intrinsic spin Hall conductivity (SHC) in nonmagnetic hcp metals and in antiferromagnetic Cr. For most of the metals of this study we find large anisotropies. We derive the general relation between the SHC vector and the direction of spin polarization and discuss its consequences for hcp metals. Especially, it is predicted that for systems where the SHC changes sign due to the anisotropy the spin Hall effect may be tuned such that the spin polarization is parallel either to the electric field or to the spin current.

preprint2010arXiv

Complex magnetic behavior and high spin polarization in Fe(3-x)Mn(x)Si alloys

Fe3Si is a ferromagnetic material with possible applications in magnetic tunnel junctions. When doped with Mn, the material shows a complex magnetic behavior, as suggested by older experiments. We employed the Korringa-Kohn-Rostoker (KKR) Green function method within density-functional theory (DFT) in order to study the alloy Fe(3-x)Mn(x)Si, with 0 < x < 1. Chemical disorder is described within the coherent potential approximation (CPA). In agreement with experiment, we find that the Mn atoms align ferromagnetically to the Fe atoms, and that the magnetization and Curie temperature drop with increasing Mn-concentration $x$. The calculated spin polarization P at the Fermi level varies strongly with x, from P=-0.3 at x=0 (ordered Fe3Si) through P=0 at x=0.28, to P=+1 for x>0.75; i.e., at high Mn concentrations the system is half-metallic. We discuss the origin of the trends of magnetic moments, exchange interactions, Curie temperature and the spin polarization.

preprint2010arXiv

Efficient implementation of the GW approximation within the all-electron FLAPW method

We present an implementation of the GW approximation for the electronic self-energy within the full-potential linearized augmented-plane-wave (FLAPW) method. The algorithm uses an all-electron mixed product basis for the representation of response matrices and related quantities. This basis is derived from the FLAPW basis and is exact for wave-function products. The singularity of the bare and screened interaction potentials gives rise to a numerically important self-energy contribution, which we treat analytically to achieve good convergence with respect to the k-point sampling. As numerical realizations of the GW approximation typically suffer from the high computational expense required for the evaluation of the nonlocal and frequency-dependent self-energy, we demonstrate how the algorithm can be made very efficient by exploiting spatial and time-reversal symmetry as well as by applying an optimization of the mixed product basis that retains only the numerically important contributions of the electron-electron interaction. Furthermore, we demonstrate that one can employ an extrapolar approximation for high-lying states to reduce the number of empty states that must be taken into account explicitly in the construction of the polarization function and the self-energy. We show convergence tests, CPU timings, and results for prototype semiconductors and insulators as well as ferromagnetic nickel.

preprint2010arXiv

Hybrid functionals within the all-electron FLAPW method: implementation and applications of PBE0

We present an efficient implementation of the PBE0 hybrid functional within the full-potential linearized augmented-plane-wave (FLAPW) method. The Hartree-Fock exchange term, which is a central ingredient of hybrid functionals, gives rise to a computationally expensive nonlocal potential in the one-particle Schroedinger equation. The matrix elements of this exchange potential are calculated with the help of an auxiliary basis that is constructed from products of FLAPW basis functions. By representing the Coulomb interaction in this basis the nonlocal exchange term becomes a Brillouin-zone (BZ) sum over vector-matrix-vector products. We show that the Coulomb matrix can be made sparse by a suitable unitary transformation of the auxiliary basis, which accelerates the computation of the vector-matrix-vector products considerably. Additionally, we exploit spatial and time-reversal symmetry to identify the nonvanishing exchange matrix elements in advance and to restrict the k summations for the nonlocal potential to an irreducible set of k points. Favorable convergence of the self-consistent-field cycle is achieved by a nested density-only and density-matrix iteration scheme. We discuss the convergence with respect to the parameters of our numerical scheme and show results for a variety of semiconductors and insulators, including oxide materials, where the PBE0 hybrid functional improves the band gaps and the description of localized states in comparison with the PBE functional. Furthermore, we find that in contrast to conventional local exchange-correlation functionals ferromagnetic EuO is correctly predicted to be a semiconductor.

preprint2010arXiv

Real-space electronic-structure calculations with full-potential all-electron precision for transition-metals

We have developed an efficient computational scheme utilizing the real-space finite-difference formalism and the projector augmented-wave (PAW) method to perform precise first-principles electronic-structure simulations based on the density functional theory for systems containing transition metals with a modest computational effort. By combining the advantages of the time-saving double-grid technique and the Fourier filtering procedure for the projectors of pseudopotentials, we can overcome the egg box effect in the computations even for first-row elements and transition metals, which is a problem of the real-space finite-difference formalism. In order to demonstrate the potential power in terms of precision and applicability of the present scheme, we have carried out simulations to examine several bulk properties and structural energy differences between different bulk phases of transition metals, and have obtained excellent agreement with the results of other precise first-principles methods such as a plane wave based PAW method and an all-electron full-potential linearized augmented plane wave (FLAPW) method.

preprint2010arXiv

Spin-polarization of platinum (111) induced by the proximity to cobalt nanostripes

We measured a spin polarization above a Pt (111) surface in the vicinity of a Co nanostripe by spin-polarized scanning tunneling spectroscopy. The spin polarization is exponentially decaying away from the Pt/Co interface and is detectable at distances larger than 1 nm. By performing self-consistent ab-initio calculations of the electronic-structure for a related model system we reveal the interplay between the induced magnetic moments within the Pt surface and the spin-resolved electronic density of states above the surface.

preprint2010arXiv

Theory of real space imaging of Fermi surfaces

A scanning tunneling microscope can be used to visualize in real space Fermi surfaces with buried impurities far below substrates acting as local probes. A theory describing this feature is developed based on the stationary phase approximation. It is demonstrated how a Fermi surface of a material acts as a mirror focusing electrons that scatter at hidden impurities.

preprint2010arXiv

Thermally activated magnetization reversal in monoatomic magnetic chains on surfaces studied by classical atomistic spin-dynamics simulations

We analyze the spontaneous magnetization reversal of supported monoatomic chains of finite length due to thermal fluctuations via atomistic spin-dynamics simulations. Our approach is based on the integration of the Landau-Lifshitz equation of motion of a classical spin Hamiltonian at the presence of stochastic forces. The associated magnetization lifetime is found to obey an Arrhenius law with an activation barrier equal to the domain wall energy in the chain. For chains longer than one domain-wall width, the reversal is initiated by nucleation of a reversed magnetization domain primarily at the chain edge followed by a subsequent propagation of the domain wall to the other edge in a random-walk fashion. This results in a linear dependence of the lifetime on the chain length, if the magnetization correlation length is not exceeded. We studied chains of uniaxial and tri-axial anisotropy and found that a tri-axial anisotropy leads to a reduction of the magnetization lifetime due to a higher reversal attempt rate, even though the activation barrier is not changed.

preprint2010arXiv

Wannier-function approach to spin excitations in solids

We present a computational scheme to study spin excitations in magnetic materials from first principles. The central quantity is the transverse spin susceptibility, from which the complete excitation spectrum, including single-particle spin-flip Stoner excitations and collective spin-wave modes, can be obtained. The susceptibility is derived from many-body perturbation theory and includes dynamic correlation through a summation over ladder diagrams that describe the coupling of electrons and holes with opposite spins. In contrast to earlier studies, we do not use a model potential with adjustable parameters for the electron-hole interaction but employ the random-phase approximation. To reduce the numerical cost for the calculation of the four-point scattering matrix we perform a projection onto maximally localized Wannier functions, which allows us to truncate the matrix efficiently by exploiting the short spatial range of electronic correlation in the partially filled d or f orbitals. Our implementation is based on the FLAPW method. Starting from a ground-state calculation within the LSDA, we first analyze the matrix elements of the screened Coulomb potential in the Wannier basis for the 3d transition-metal series. In particular, we discuss the differences between a constrained nonmagnetic and a proper spin-polarized treatment for the ferromagnets Fe, Co, and Ni. The spectrum of single-particle and collective spin excitations in fcc Ni is then studied in detail. The calculated spin-wave dispersion is in good overall agreement with experimental data and contains both an acoustic and an optical branch for intermediate wave vectors along the [100] direction. In addition, we find evidence for a similar double-peak structure in the spectral function along the [111] direction.

preprint2009arXiv

Exchange coupling in transition-metal nano-clusters on Cu(001) and Cu(111) surfaces

We present results of density-functional calculations on the magnetic properties of Cr, Mn, Fe and Co nano-clusters (1 to 9 atoms large) supported on Cu(001) and Cu(111). The inter-atomic exchange coupling is found to depend on competing mechanisms, namely ferromagnetic double exchange and antiferromagnetic kinetic exchange. Hybridization-induced broadening of the resonances is shown to be important for the coupling strength. The cluster shape is found to weaken the coupling via a mechanism that comprises the different orientation of the atomic d-orbitals and the strength of nearest-neighbour hopping. Especially in Fe clusters, a correlation of binding energy and exchange coupling is also revealed.

preprint2006arXiv

Changing the Magnetic Configurations of Nanoclusters Atom-by-Atom

The Korringa-Kohn-Rostoker Green (KKR) function method for non-collinear magnetic structures was applied on Mn and Cr ad-clusters deposited on the Ni(111) surface. By considering various dimers, trimers and tetramers, a large amount of collinear and non-collinear magnetic structures is obtained. Typically all compact clusters have very small total moments, while the more open structures exhibit sizeable total moments, which is a result of the complex frustration mechanism in these systems. Thus, as the motion of a single adatom changes the cluster structure from compact to open and vice versa, this can be considered as a magnetic switch, which via the local exchange field of the adatom allows to switch the cluster moment on and off, and which might be useful for future nanosize information storage.

preprint2006arXiv

Non-collinear magnetism of Cr nanostructures on Fe$_{3ML}$/Cu(001): first--principles and experimental investigations

A combined experimental, using X-ray Magnetic Circular Dichroism and theoretical investigation, using full-potential Korringa-Kohn-Rostoker (KKR) Green function method, is carried out to study the spin structure of small magnetic Cr adatom--clusters on the surface of 3 monolayers of $fcc$ Fe deposited on Cu(001). The exchange interaction between the different Cr adatoms as well as between the Cr atoms and the Fe atoms is of antiferromagnetic nature and of comparable magnitude, leading due to frustration to complex non-collinear magnetic configurations. The presence of non-collinear magnetic coupling obtained by {\it ab initio} calculations is confirmed by the experimental results.

preprint2005arXiv

Surface state scattering by adatoms on noble metals

When surface state electrons scatter at perturbations, such as magnetic or nonmagnetic adatoms or clusters on surfaces, an electronic resonance, localized at the adatom site, can develop below the bottom of the surface state band for both spin channels. In the case of adatoms, these states have been found very recently in scanning tunneling spectroscopy experiments\cite{limot,olsson} for the Cu(111) and Ag(111) surfaces. Motivated by these experiments, we carried out a systematic theoretical investigation of the electronic structure of these surface states in the presence of magnetic and non-magnetic atoms on Cu(111). We found that Ca and all 3$d$ adatoms lead to a split-off state at the bottom of the surface band which is, however, not seen for the $sp$ elements Ga and Ge. The situation is completely reversed if the impurities are embedded in the surface: Ga and Ge are able to produce a split-off state whereas the 3$d$ impurities do not. The resonance arises from the s-state of the impurities and is explained in terms of strength and interaction nature (attraction or repulsion) of the perturbing potential.

Stefan Blügel

What is connected

Connect this record

See the researcher in context

Building this map preview

125 published item(s)

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

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

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

Chiral standing spin waves in skyrmion lattice

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

Diversity of states in a chiral magnet nanocylinder

Driving spin chirality by electron dynamics in laser-excited antiferromagnets

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

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

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

Link between Weyl-fermion chirality and spin texture

Long-Range Orbital Magnetoelectric Torque in Ferromagnets

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

Non-Collinear Spin Current for Switching of Chiral Magnetic Textures

Proximity induced superconductivity in a topological insulator

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

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

$Z_2$ topology of bismuth

Charge and spin photocurrents in the Rashba model

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

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

Heisenberg representation of nonthermal ultrafast laser excitation of magnetic precessions

Laser-induced torques in spin spirals

Magnetic particle-antiparticle creation and annihilation

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

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

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

Atomistic perspective of long lifetimes of small skyrmions at room temperature

Effective Seiberg-Witten gauge theory of noncollinear magnetism

Electric dipole moment as descriptor for interfacial Dzyaloshinskii-Moriya interaction

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

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

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

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

Solution to the modified Helmholtz equation for arbitrary periodic charge densities

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

Topological magneto-optical effects and their quantization in noncoplanar antiferromagnets

Topological-chiral magnetic interactions driven by emergent orbital magnetism

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

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

Wannier90 as a community code: new features and applications

Approximations to the exact exchange potential: KLI versus semilocal

Chirality-driven orbital magnetic moments as a new probe for topological magnetic structures

Fermi surfaces, spin-mixing parameter, and colossal anisotropy of spin relaxation in transition metals from ab initio theory

Giant spin Nernst effect induced by resonant scattering at surfaces of metallic films

Influence of complex disorder on skew-scattering Hall effects in $L1_0$-ordered FePt alloy

Interlayer exchange coupling -- a general scheme turning chiral magnets into magnetic multilayers carrying atomic-scale skyrmions

Laser-induced torques in metallic ferromagnets

Microscopic theory of the residual surface resistivity of Rashba electrons

Report on imaging of individual skyrmions of MML systems made by different techniques

Role of Berry phase theory for describing orbital magnetism: From magnetic heterostructures to topological orbital ferromagnets

Spin-orbit torques and spin accumulation in FePt/Pt and Co/Cu thin films from first principles: the role of impurities

Strong spin-orbit fields and Dyakonov-Perel spin dephasing in supported metallic films

Sub-molecular modulation of a 4f driven Kondo resonance by surface-induced asymmetry

Switching of chiral magnetic skyrmions by picosecond magnetic field pulses via transient topological states

The inverse thermal spin-orbit torque and the relation of the Dzyaloshinskii-Moriya interaction to ground-state energy currents

Comparison between exact and semilocal exchange potentials: An all-electron study for solids

Direct and inverse spin-orbit torques

Electronic phase transitions of bismuth under strain from relativistic self-consistent GW calculations

Higher dimensional Wannier functions of multi-parameter Hamiltonians

Magnetic properties of 2D nickel nanostrips: structure dependent magnetic anomaly

Oxygen-enabled control of Dzyaloshinskii-Moriya Interaction in ultra-thin magnetic films

Quasi 2D electronic states with high spin-polarization in centrosymmetric MoS$_2$ bulk crystals

Realization of a vertical topological p-n junction in epitaxial $\mathrm{Sb_2Te_3 / Bi_2Te_3}$ heterostructures

Relativistic dynamical spin excitations of magnetic adatoms

Skew scattering in dilute ferromagnetic alloys

Bulk Band Structure of Bi$_2$Te$_3$

Dzyaloshinskii-Moriya interaction and chiral magnetism in 3$d$-5$d$ zig-zag chains: Tight-binding model and ab initio calculations

First-principles analysis of a homo-chiral cycloidal magnetic structure in a monolayer Cr on W(110)

Quantum well states and amplified spin-dependent Friedel oscillations in thin films

Spin relaxation and spin Hall transport in 5d transition-metal ultrathin films

Spin Spirals in Surface Alloys on Ru(0001): A First-principles Study

Spin transport and spin-caloric effects in (Cr,Zn)Te half-metallic nanostructures: Effect of spin disorder at elevated temperatures from first principles

Spin-caloric transport properties of cobalt nanostructures: spin disorder effects from first principles