# Gilbert damping in non-collinear magnetic systems

## Abstract

The modification of the magnetization dissipation or Gilbert damping caused by an inhomogeneous magnetic structure and expressed in terms of a wave vector dependent tensor is investigated by means of linear response theory. A corresponding expression for in terms of the electronic Green function has been developed giving in particular the leading contributions to the Gilbert damping linear and quadratic in . Numerical results for realistic systems are presented that have been obtained by implementing the scheme within the framework of the fully relativistic KKR (Korringa-Kohn-Rostoker) band structure method. Using the multilayered system (Cu/FeCo/Pt) as an example for systems without inversion symmetry we demonstrate the occurrence of non-vanishing linear contributions. For the alloy system bcc FeCo having inversion symmetry, on the other hand, only the quadratic contribution is non-zero. As it is shown, this quadratic contribution does not vanish even if the spin-orbit coupling is suppressed, i.e. it is a direct consequence of the non-collinear spin configuration.

###### pacs:

71.15.-m,71.55.Ak, 75.30.Ds^{1}

## I Introduction

The magnetization dissipation in magnetic materials is conventionally characterized by means of the Gilbert damping (GD) tensor that enters the Landau-Lifshitz-Gilbert (LLG) equation Gilbert (2004). This positive-definite second-rank tensor depends in general on the magnetization direction. It is well established that in the case of spatially uniformly magnetized ferromagnetic (FM) metals two regimes of slow magnetization dynamics can be distinguished, which are governed by different mechanisms of dissipation Kambersky (1970); Fähnle and Steiauf (2006); Gilmore et al. (2007): a conductivity-like behaviour occuring in the limiting case of ordered compounds that may be connected to the Fermi breathing mechanism and a resistivity-like behaviour shown by materials with appreciable structural, chemical or temperature induced disorder and connected to a spin-flip scattering mechanism. An important issue is that both mechanisms are determined by the spin-orbit coupling in the system (see e.g. Mankovsky et al. (2013); Kambersky (1970); Gilmore et al. (2007)). During the last years, it was demonstrated by various authors that first-principles calculations for the GD parameter for collinear ferromagntic materials allow to cover both regimes without use of any phenomenological parameters. In fact, in spite of the differences concerning the formulation for the damping parameter and the corresponding implementaion Starikov et al. (2010); Ebert et al. (2011); Turek et al. (2015), the numerical results are in general in rather good agreement with each other as well as with experiment.

In the case of a pronounced non-collinear magnetic texture, e.g. in the case of domain walls or topologically nontrivial magnetic configurations like skyrmions, the description of the magnetization dissipation assuming a spatial-invariant tensor is incomplete, and a non-local character of GD tensor in such systems has to be taken into account Tserkovnyak et al. (2009); Hals et al. (2015); Hals and Brataas (2014). This implies that the dissipative torque on the magnetization should be represented by the expression of the following general form Umetsu et al. (2015):

(1) |

In the case of a magnetic texture varying slowly in space, however, an expansion of the damping parameter in terms of the magnetization density and its gradients Hals and Brataas (2014) is nevertheless appropriate:

where the first term stands for the conventional isotropic GD and the second term is associated with the magneto-crystalline anisotropy (MCA). The third so-called chiral term is non-vanishing in non-centrosymmetric systems. The important role of this contribution to the damping was demonstrated experimentally when investigating the field-driven domain wall (DW) motion in asymmetric Pt/Co/Pt trilayers Jue et al. (2015).

As an alternative to the expansion in Eq. (I) one can discuss the Fourier transform of the damping parameter characterizing inhomogeneous magnetic systems, which enter the spin dynamics equation

(3) |

In this formulation the term linear in is the first chiral term appearing in the expansion of in powers of . Furthermore, it is important to note that it is directly connected to the term in Eq. (I).

By applying a gauge field theory, the origin of the non-collinear corrections to the GD can be ascribed to the emergent electromagnetic field created in the time-dependent magnetic texture Zhang and Zhang (2009); Tatara and Nakabayashi (2014). Such an emergent electromagnetic field gives rise to a spin current whose divergence characterizes the change of the angular momentum in the system. This allows to discuss the impact of non-collinearity on the GD via a spin-pumping formulation Tserkovnyak and Mecklenburg (2008); Zhang and Zhang (2009); Tserkovnyak et al. (2009). Some details of the physics behind this effect depend on the specific properties of the material considered. Accordingly, different models for magnetisation dissipation were discussed in the literature Tserkovnyak et al. (2009); Zhang and Zhang (2009); Zang et al. (2011); Umetsu et al. (2015); Chiba et al. (2015); Akosa et al. (2016). Non-centrosymmetric two-dimensional systems for which the Rashba-like spin-orbit coupling plays an important role have received special interest in this context. They have been discussed in particular by Akosa et al. Akosa et al. (2016), in order to explain the origin of chiral GD in the presence of a chiral magnetic structure.

The fourth term on the r.h.s. of Eq. (I) corresponds to a quadratic term of an expansion of with respect to . It was investigated for bulk systems with non-magnetic Hankiewicz et al. (2008) and magnetic Tserkovnyak et al. (2009) impurity atoms, for which the authors have shown on the basis of model consideration that it can give a significant correction to the homogeneous GD in the case of weak metallic ferromagnets. In striking contrast to the uniform part of the GD this contribution does not require a non-vanishing spin-orbit interaction.

To our knowledge, only very few ab-initio investigations on the Gilbert damping in non-collinear magnetic systems along the lines sketched above have been reported so far in the literature. Yuan et al. Yuan et al. (2014) calculated the in-plane and out-of-plane damping parameters in terms of the scattering matrix for permalloy in the presence of Néel and Bloch domain walls. Freimuth et al. Freimuth et al. (2017), discuss the properties of a -dependent Gilbert damping calculated for the one-dimensional Rashba model in the presence of the Néel-type non-collinear magnetic exchange field, demonstrating different GD for left-handed and right-handed DWs. Here we extend the formalism developed before to deal with the GD in ferromagnets Starikov et al. (2010), to get access to non-collinear system. The formalism based on linear response theory allows to expand the GD parameters with respect to a modulation of the magnetization expressed in terms of a wave vector . Corresponding numerical results will be presented and discussed.

## Ii Gilbert damping for non-collinear magnetization

In the following we focus on the intrinsic contribution to the Gilbert damping, excluding spin current induced magnetization dissipation which occurs in the presence of an external electric field. For the considerations on the magnetization dissipation an adiabatic variation of the magnetization in the time and space domain is assumed. Moreover, it is assumed that the magnitude of the local magnetic moments is unchanged during a change of the magnetization, i.e. the exchange field should be strong enough to separate transverse and longitudinal parts of the magnetic susceptibility. With these restrictions, the non-local Gilbert damping can be determined in terms of the spin susceptibility tensor

(4) |

where is the - and -dependent spin operator and reduced units have been used (). With this, the Fourier transformation of the real-space Gilbert damping can be represented by the expression Qian and Vignale (2002); Hankiewicz et al. (2007)

(5) |

Here is the gyromagnetic ratio, is the equilibrium magnetization and is the volume of the system. In order to avoid the calculation of the dynamical magnetic susceptibility tensor , which is the Fourier transformed of the real space susceptibility , it is convenient to represent in Eq. (5), in terms of a correlation function of time derivatives of . As corresponds to the torque , that may include non-dissipative and dissipative parts, one may consider instead the torque-torque correlation function Mori and Kawasaki (1962); Mahan (1970); Hasegawa (1971); Hankiewicz et al. (2007).

Assuming the magnetization direction parallel to one obtains the expression for the Gilbert damping

(6) |

where is the transverse Levi-Civita tensor. This implies the following relationship of the tensor elements with the elements of the torque-torque correlation tensor : and Hankiewicz et al. (2007).

Using Kubo’s linear response theory in the Matsubara representation and taking into account the translational symmetry of a solid the torque-torque correlation function can be expressed by (see, e.g. Mahan (2000)):

(7) | |||||

where is the Matsubara Green function and indicates a configurational average required in the presence of any disorder (chemical, structural or magnetic) in the system. Using a Lehman representation for the Green function Mahan (2000)

(8) |

with the retarded Green function and using the relation

for the sum over the Matsubara poles in Eq. (7), the torque-torque correlation function is obtained as:

(9) |

Perfoming finally the analytical continuation one arrives at the expression

(10) | |||||

for the imaginary part of the correlation function with . Accordingly one gets for the diagonal elements of Gilbert damping tensor the expression

(11) | |||||

where the index of the torque operator is related to the index according to Eq. 6, and the auxiliary functions

(12) |

expressed in terms of the retarded and advanced Green functions, and , respectively.

To account properly for the impact of spin-orbit coupling when dealing with Eqs. (11) and (12) a description of the electronic structure based on the fully relativistic Dirac formalism is used. Working within the framework of local spin density formalism (LSDA) this implies for the Hamiltonian the form Ebert (2000):

(13) |

Here and are the standard Dirac matrices, denotes the vector of relativistic Pauli matrices, is the relativistic momentum operator Rose (1961) and the functions and are the spin-averaged and spin-dependent parts, respectively, of the LSDA potential MacDonald and Vosko (1979) with giving the orientation of the magnetisation.

With the Dirac Hamiltonian given by Eq. (13), the torque operator may be written as . Furthermore, the Green functions entering Eqs. (11) and (12) are determined using the spin-polarized relativistic version of multiple scattering theory Ebert (2000); Ebert et al. (2016) with the real space representation of the retarded Green function given by:

(14) | |||||

Here refer to atomic cells centered at sites and , respectively, where is a function centered at the corresponding lattice vector . The four-component wave functions () are regular (irregular) solutions to the single-site Dirac equation labeled by the combined quantum numbers , with and being the spin-orbit and magnetic quantum numbers Rose (1961). Finally, is the so-called scattering path operator that transfers an electronic wave coming in at site into a wave going out from site with all possible intermediate scattering events accounted for.

Using matrix notation with respect to , this leads to the following expression for the auxilary damping parameters in Eq. (12):

(15) |

In the case of a uniform magnetization, i.e. for one obviously gets an expression for the Gilbert damping tensor as it was worked out before Ebert et al. (2011). Assuming small wave vectors, the term ) can be expanded w.r.t. leading to the series

(16) |

that results in a corresponding expansion for the Gilbert damping:

(17) |

with the following expansion coefficients:

(18) | |||||

(19) | |||||

(20) |

and with the g-factor in terms of the spin and orbital moments, and , respectively, and the total magnetic moment . The numerically cumbersome term in Eq. (20), that involves the second order derivative of the matrix of -dependent scattering path operator , can be reformulated by means of an integration by parts:

leading to the much more convenient expression:

(21) |

## Iii Results and discussions

The scheme presented above to deal with the Gilbert damping in non-collinear systems has been implemented within the SPR-KKR program package H. Ebert et al. (2017). To examine the importance of the chiral correction to the Gilbert damping a first application of Eq. (19) has been made for the multilayer system (Cu/FeCo/Pt) seen as a non-centrosymmetric model system. The calculated zero-order (uniform) GD parameter and the corresponding first-order (chiral) correction term for are plotted in Fig. 1 top and bottom, respectively, as a function of the Fe concentration .

Both terms, and , increase approaching the pure limits w.r.t. the FeCo alloy subsystem. In the case of the uniform parameter , this increase is associated with the dominating breathing Fermi-surface damping mechanism. This implies that the modification of the Fermi surface (FS) induced by the spin-orbit coupling (SOC) follows the magnetization direction that slowly varies with time. An additional contribution to the GD, having a similar origin, occurs for the non-centrosymmertic systems with helimagnetic structure. In this case, the features of the electronic structure governed by the lack of inversion symmetry result in a FS modification dependent on the helicity of the magnetic structure. This implies a chiral contribution to the GD which can be associated with the term proportional to the gradient of the magnetization. Obviously, this additional modification of the FS and the associated mechanism for the GD does not show up for a uniform ferromagnet. As is caused by the SOC one can expect that it vanishes for vanishing SOC. This was indeed demonstrated before Mankovsky et al. (2013). The same holds also for that is cased by SOC as well.

Another system considered is the ferromagnetic alloy system bcc FeCo. As this system has inversion symmetry the first-order term should vanish. This expectation could also be confirmed by calculations that account for the SOC. The next non-vanishing term of the expansion of the GD is the term . The corresponding second-order term is plotted in Fig. 2 (bottom) together with the zero-order term (top).

The bottom panel shows in addition results for that have been obtained by calculations with the SOC suppressed. As one notes the results for the full SOC and for SOC suppressed are very close to each other. The small difference between the curves for that reason have to be ascribed to the hybridization of the spin-up and spin-down subsystems due to SOC. As discussed in the literature Hankiewicz et al. (2008); Tserkovnyak et al. (2009); Zang et al. (2011) a non-collinear magnetic texture has a corresponding consequence but a much stronger impact here. In contrast to the GD in uniform FM systems where SOC is required to break the total spin conservation in the system, is associated with the spin-pumping effect that can be ascribed to an emergent electric field created in the non-uniform magnetic system. In this case magnetic dissipation occurs due to the misalignment of the electron spin following the dynamic magnetic profile and the magnetization orientation at each atomic site, leading to the dephasing of electron spins Tserkovnyak and Mecklenburg (2008)

## Iv Summary

To summarize, expressions for corrections to the GD of homogeneous systems were derived which are expected to contribute in the case of non-collinear magnetic systems. The expression for the GD parameter seen as a function of the wave vector is expanded in powers of . In the limit of weakly varying magnetic textures, this leads to the standard uniform term, , and the first- and second-order corrections, and , respectively. Model calculations confirmed that a non-vanishing value for can be expected for systems without inversion symmetry. In addition, SOC has been identified as the major source for this term. The second-order term, on the other hand, may also show up for systems with inversion symmetry. In this case it was demonstrated by numerical work, that SOC plays only a minor role for , while the non-collinearity of the magnetization plays the central role.

## V Acknowledgement

Financial support by the DFG via SFB 1277 (Emergente relativistische Effekte in der Kondensierten Materie) is gratefully acknowledged.

### Footnotes

- preprint: APS/123-QED

### References

- T. L. Gilbert, IEEE Transactions on Magnetics 40, 3443 (2004).
- V. Kambersky, Can. J. Phys. 48, 2906 (1970), URL http://www.nrcresearchpress.com/doi/abs/10.1139/p70-361.
- M. Fähnle and D. Steiauf, Phys. Rev. B 73, 184427 (2006).
- K. Gilmore, Y. U. Idzerda, and M. D. Stiles, Phys. Rev. Lett. 99, 027204 (2007), URL http://link.aps.org/doi/10.1103/PhysRevLett.99.027204.
- S. Mankovsky, D. Ködderitzsch, G. Woltersdorf, and H. Ebert, Phys. Rev. B 87, 014430 (2013), URL http://link.aps.org/doi/10.1103/PhysRevB.87.014430.
- A. A. Starikov, P. J. Kelly, A. Brataas, Y. Tserkovnyak, and G. E. W. Bauer, Phys. Rev. Lett. 105, 236601 (2010), URL http://link.aps.org/doi/10.1103/PhysRevLett.105.236601.
- H. Ebert, S. Mankovsky, D. Ködderitzsch, and P. J. Kelly, Phys. Rev. Lett. 107, 066603 (2011), eprint http://arxiv.org/abs/1102.4551v1, URL http://link.aps.org/doi/10.1103/PhysRevLett.107.066603.
- I. Turek, J. Kudrnovský, and V. Drchal, Phys. Rev. B 92, 214407 (2015), URL http://link.aps.org/doi/10.1103/PhysRevB.92.214407.
- Y. Tserkovnyak, E. M. Hankiewicz, and G. Vignale, Phys. Rev. B 79, 094415 (2009), URL https://link.aps.org/doi/10.1103/PhysRevB.79.094415.
- K. M. D. Hals, K. Flensberg, and M. S. Rudner, Phys. Rev. B 92, 094403 (2015), URL https://link.aps.org/doi/10.1103/PhysRevB.92.094403.
- K. M. D. Hals and A. Brataas, Phys. Rev. B 89, 064426 (2014), URL https://link.aps.org/doi/10.1103/PhysRevB.89.064426.
- N. Umetsu, D. Miura, and A. Sakuma, Phys. Rev. B 91, 174440 (2015), URL https://link.aps.org/doi/10.1103/PhysRevB.91.174440.
- E. Jue, C. K. Safeer, M. Drouard, A. Lopez, P. Balint, L. Buda-Prejbeanu, O. Boulle, S. Auffret, A. Schuhl, A. Manchon, et al., Nature Materials 15, 272 (2015), URL http://dx.doi.org/10.1038/nmat4518.
- S. Zhang and S. S.-L. Zhang, Phys. Rev. Lett. 102, 086601 (2009), URL https://link.aps.org/doi/10.1103/PhysRevLett.102.086601.
- G. Tatara and N. Nakabayashi, Journal of Applied Physics 115, 172609 (2014), eprint http://dx.doi.org/10.1063/1.4870919, URL http://dx.doi.org/10.1063/1.4870919.
- Y. Tserkovnyak and M. Mecklenburg, Phys. Rev. B 77, 134407 (2008), URL https://link.aps.org/doi/10.1103/PhysRevB.77.134407.
- J. Zang, M. Mostovoy, J. H. Han, and N. Nagaosa, Phys. Rev. Lett. 107, 136804 (2011), URL https://link.aps.org/doi/10.1103/PhysRevLett.107.136804.
- T. Chiba, G. E. W. Bauer, and S. Takahashi, Phys. Rev. B 92, 054407 (2015), URL https://link.aps.org/doi/10.1103/PhysRevB.92.054407.
- C. A. Akosa, I. M. Miron, G. Gaudin, and A. Manchon, Phys. Rev. B 93, 214429 (2016), URL https://link.aps.org/doi/10.1103/PhysRevB.93.214429.
- E. M. Hankiewicz, G. Vignale, and Y. Tserkovnyak, Phys. Rev. B 78, 020404 (2008), URL https://link.aps.org/doi/10.1103/PhysRevB.78.020404.
- Z. Yuan, K. M. D. Hals, Y. Liu, A. A. Starikov, A. Brataas, and P. J. Kelly, Phys. Rev. Lett. 113, 266603 (2014), URL https://link.aps.org/doi/10.1103/PhysRevLett.113.266603.
- F. Freimuth, S. Blügel, and Y. Mokrousov, Phys. Rev. B 96, 104418 (2017), URL https://link.aps.org/doi/10.1103/PhysRevB.96.104418.
- Z. Qian and G. Vignale, Phys. Rev. Lett. 88, 056404 (2002), URL http://link.aps.org/doi/10.1103/PhysRevLett.88.056404.
- E. M. Hankiewicz, G. Vignale, and Y. Tserkovnyak, Phys. Rev. B 75, 174434 (2007), URL https://link.aps.org/doi/10.1103/PhysRevB.75.174434.
- H. Mori and K. Kawasaki, Progr. Theor. Phys. 27, 529 (1962), URL http://ptp.ipap.jp/link?PTP/27/529/.
- G. Mahan, Journal of Physics and Chemistry of Solids 31, 1477 (1970), ISSN 0022-3697, URL http://www.sciencedirect.com/science/article/pii/0022369770900326.
- M. Hasegawa, Journal of the Physical Society of Japan 31, 649 (1971), eprint https://doi.org/10.1143/JPSJ.31.649, URL https://doi.org/10.1143/JPSJ.31.649.
- G. D. Mahan, Many-particle physics, Physics of Solids and Liquids (Springer, New York, 2000), ISBN ISBN 978-1-4757-5714-9, URL https://openlibrary.org/works/OL4474340W/Many-particle_physics.
- H. Ebert, in Electronic Structure and Physical Properties of Solids, edited by H. Dreyssé (Springer, Berlin, 2000), vol. 535 of Lecture Notes in Physics, p. 191.
- M. E. Rose, Relativistic Electron Theory (Wiley, New York, 1961), URL http://openlibrary.org/works/OL3517103W/Relativistic_electron_theory.
- A. H. MacDonald and S. H. Vosko, J. Phys. C: Solid State Phys. 12, 2977 (1979), URL http://iopscience.iop.org/0022-3719/12/15/007.
- H. Ebert, J. Braun, D. Ködderitzsch, and S. Mankovsky, Phys. Rev. B 93, 075145 (2016), URL http://link.aps.org/doi/10.1103/PhysRevB.93.075145.
- H. Ebert et al., The Munich SPR-KKR package, version 7.7, http://olymp.cup.uni-muenchen.de/ak/ebert/SPRKKR (2017), URL http://olymp.cup.uni-muenchen.de/ak/ebert/SPRKKR.