Microscopic Theory for Coupled Atomistic Magnetization and Lattice Dynamics

A coupled atomistic spin and lattice dynamics approach is developed which merges the dynamics of these two degrees of freedom into a single set of coupled equations of motion. The underlying microscopic model comprises local exchange interactions between the electron spin and magnetic moment and the local couplings between the electronic charge and lattice displacements. An effective action for the spin and lattice variables is constructed in which the interactions among the spin and lattice components are determined by the underlying electronic structure. In this way, expressions are obtained for the electronically mediated couplings between the spin and lattice degrees of freedom, besides the well known inter-atomic force constants and spin-spin interactions. These former susceptibilities provide an atomistic ab initio description for the coupled spin and lattice dynamics. It is important to notice that this theory is strictly bilinear in the spin and lattice variables and provides a minimal model for the coupled dynamics of these subsystems and that the two subsystems are treated on the same footing. Questions concerning time-reversal and inversion symmetry are rigorously addressed and it is shown how these aspects are absorbed in the tensor structure of the interaction fields. By means of these results regarding the spin-lattice coupling, simple explanations of ionic dimerization in double anti-ferromagnetic materials, as well as, charge density waves induced by a non-uniform spin structure are given. In the final parts, a set of coupled equations of motion for the combined spin and lattice dynamics are constructed, which subsequently can be reduced to a form which is analogous to the Landau-Lifshitz-Gilbert equations for spin dynamics and damped driven mechanical oscillator for the ...Comment: 22 pages, including 7 pages of Appendix and references, 6 figure

Gilbert damping tensor within the breathing Fermi surface model: anisotropy and non-locality

In magnetization dynamics, the Gilbert damping α is often taken as a parameter. We report on a theoretical investigation of α, taking into account crystal symmetries, spin–orbit coupling and thermal reservoirs. The tensor alpha is calculated within the Kamberský breathing Fermi-surface model. The computations are performed within a tight-binding electronic structure approach for the bulk and semi-infinite systems. Slater–Koster parameters are obtained by fitting the electronic structure to first-principles results obtained within the multiple-scattering theory. We address the damping tensor for the bulk and surfaces of the transition metals Fe and Co. The role of various contributions are investigated: intra- and interband transitions, electron and magnetic temperature as well as surface orientation. Our results reveal a complicated non-local, anisotropic damping that depends on all three thermal reservoirs

Exchange interaction and its tuning in magnetic binary chalcogenides

Using a first-principles Green's function approach we study magnetic properties of the magnetic binary chalcogenides Bi2Te3, Bi2Se3, and Sb2Te3. The magnetic coupling between transition-metal impurities is long-range, extends beyond a quintuple layer, and decreases with increasing number of d electrons per 3d atom. We find two main mechanisms for the magnetic interaction in these materials: the indirect exchange interaction mediated by free carriers and the indirect interaction between magnetic moments via chalcogen atoms. The calculated Curie temperatures of these systems are in good agreement with available experimental data. Our results provide deep insight into magnetic interactions in magnetic binary chalcogenides and open a way to design new materials for promising applications

Extensive degeneracy, Coulomb phase and magnetic monopoles in an artificial realization of the square ice model

Artificial spin ice systems have been introduced as a possible mean to investigate frustration effects in a well-controlled manner by fabricating lithographically-patterned two-dimensional arrangements of interacting magnetic nanostructures. This approach offers the opportunity to visualize unconventional states of matter, directly in real space, and triggered a wealth of studies at the frontier between nanomagnetism, statistical thermodynamics and condensed matter physics. Despite the strong efforts made these last ten years to provide an artificial realization of the celebrated square ice model, no simple geometry based on arrays of nanomagnets succeeded to capture the macroscopically degenerate ground state manifold of the corresponding model. Instead, in all works reported so far, square lattices of nanomagnets are characterized by a magnetically ordered ground state consisting of local flux-closure configurations with alternating chirality. Here, we show experimentally and theoretically, that all the characteristics of the square ice model can be observed if the artificial square lattice is properly designed. The spin configurations we image after demagnetizing our arrays reveal unambiguous signatures of an algebraic spin liquid state characterized by the presence of pinch points in the associated magnetic structure factor. Local excitations, i.e. classical analogues of magnetic monopoles, are found to be free to evolve in a massively degenerated, divergence-free vacuum. We thus provide the first lab-on-chip platform allowing the investigation of collective phenomena, including Coulomb phases and ice-like physics.Comment: 26 pages, 10 figure

Mechanisms behind large Gilbert damping anisotropies

A method with which to calculate the Gilbert damping parameter from a real-space electronic structure method is reported here. The anisotropy of the Gilbert damping with respect to the magnetic moment direction and local chemical environment is calculated for bulk and surfaces of Fe$_{50}$Co$_{50}$ alloys from first principles electronic structure in a real space formulation. The size of the damping anisotropy for Fe$_{50}$Co$_{50}$ alloys is demonstrated to be significant. Depending on details of the simulations, it reaches a maximum-minimum damping ratio as high as 200%. Several microscopic origins of the strongly enhanced Gilbert damping anisotropy have been examined, where in particular interface/surface effects stand out, as do local distortions of the crystal structure. Although theory does not reproduce the experimentally reported high ratio of 400% [Phys. Rev. Lett. 122, 117203 (2019)], it nevertheless identifies microscopic mechanisms that can lead to huge damping anisotropies

Heat-conserving three-temperature model for ultrafast demagnetization of 3d ferromagnets

We study the ultrafast magnetization dynamics of bcc Fe and fcc Co using the recently suggested heat-conserving three-temperature model (HC3TM), together with atomistic spin- and lattice dynamics simulations. It is shown that this type of Langevin-based simulation is able to reproduce observed trends of the ultrafast magnetization dynamics of fcc Co and bcc Fe, in agreement with previous findings for fcc Ni. The simulations are performed by using parameters that to as large extent as possible are obtained from electronic structure theory. The one parameter that was not calculated in this way, was the damping term used for the lattice dynamics simulations, and here a range of parameters were investigated. It is found that this term has a large influence on the details of the magnetization dynamics. The dynamics of iron and cobalt is compared with previous results for nickel and similarities and differences in the materials' behavior are analysed following the absorption of a femtosecond laser pulse. Importantly, for all elements investigated so far with this model, we obtain a linear relationship between the value of the maximally demagnetized state and the fluence of the laser pulse, which is in agreement with experiments.Comment: 9 pages, 9 figures, Submitted to Physical Review

Ultrafast demagnetization dynamics of 4f antiferromagnets

We study the ultrafast demagnetization dynamics of LnRh2Si2 (Ln = Pr, Nd, Sm, Gd, Tb, Dy, Ho) antiferromagnets (AFM) after excitation by a laser pulse, using a combination of density functional theory and atomistic spin and spin-lattice dynamics simulations. First, we calculate the Heisenberg interactions using the magnetic force theorem and compare two approaches, where the 4f states of the rare earths are treated as frozen core states or as valence states with added correlation corrections. We find marked quantitative differences in terms of predicted Curie temperature for most of the systems, especially for those with large orbital moment of the rare earth cations. This can be attributed to the importance of indirect interactions of the 4f states through the Si states, which depend on the binding energy of the 4f states and coexists with RKKY-type interactions mediated by the conduction states. However, qualitatively, both approaches agree in terms of the predicted AFM ordering at low temperatures. In the second step, the atomistic dynamics simulations are combined with a heat-conserving two-temperature model, allowing for the calculation of spin and electronic temperatures during the magnetization dynamics simulations. Despite quite different demagnetization times, magnetization dynamics of all studied LnRh2Si2 AFM exhibit similar two-step behavior, in particular, the first fast drop followed by slower demagnetization. We observe that the demagnetization amplitude depends linearly on laser fluence for low fluences, which is in agreement with experimental observations. We also investigate the impact of lattice dynamics on ultrafast demagnetization using coupled atomistic spin-lattice dynamics simulations and a heat-conserving three-temperature model, which confirm linear dependence of magnetisation on laser fluence

Spin-lattice couplings in 3⁢d ferromagnets: Analysis from first principles

Magnetoelasticity plays a crucial role in numerous magnetic phenomena, including magnetocalorics, magnon excitation via acoustic waves, and ultrafast demagnetization, or the Einstein–de Haas effect. Despite a long-standing discussion on anisotropy-mediated magnetoelastic interactions of relativistic origin, the exchange-mediated magnetoelastic parameters within an atomistic framework have only recently begun to be investigated. As a result, many of their behaviors and values for real materials remain poorly understood. Therefore, by using a proposed simple modification of the embedded cluster approach that reduces the computational complexity, we critically analyze the properties of exchange-mediated spin-lattice coupling parameters for elemental 3⁢ ferromagnets (bcc Fe, fcc Ni, and fcc Co), comparing methods used for their extraction and relating their realistic values to symmetry considerations and orbitally decomposed contributions. Additionally, we investigate the effects of noncollinearity (spin temperature) and applied pressure on these parameters. For Fe, we find that single-site rotations, associated with spin temperatures around 100 K, induce significant modifications, particularly in Dzyaloshinskii-Moriya-type couplings; in contrast, such interactions in Co and Ni remain almost configuration independent. Moreover, we demonstrate a notable change in the exchange-mediated magnetoelastic constants for Fe under isotropic contraction. Finally, the conversion between atomistic, quantum-mechanically derived parameters and the phenomenological magnetoelastic theory is discussed, which can be a useful tool towards larger and more realistic dynamics simulations involving coupled subsystems