No Arabic abstract
We present a computationally efficient general first-principles based method for spin-lattice simulations for solids. Our method is based on a combination of atomistic spin dynamics and molecular dynamics, expressed through a spin-lattice Hamiltonian where the bilinear magnetic term is expanded to second order in displacement, and all parameters are computed using density functional theory. The effect of first-order spin-lattice coupling on the magnon and phonon dispersion in bcc Fe is reported as an example, and is seen to be in good agreement with previous simulations performed with an empirical potential approach. In addition, we also illustrate the abilities of our method on a more conceptual level, by exploring dissipation-free spin and lattice motion in small magnetic clusters (a dimer, trimer and quadmer). Our method opens the door for quantitative description and understanding of the microscopic origin of many fundamental phenomena of contemporary interest, such as ultrafast demagnetization, magnetocalorics, and spincaloritronics.
We present a method for performing atomistic spin dynamic simulations. A comprehensive summary of all pertinent details for performing the simulations such as equations of motions, models for including temperature, methods of extracting data and numerical schemes for performing the simulations is given. The method can be applied in a first principles mode, where all interatomic exchange is calculated self-consistently, or it can be applied with frozen parameters estimated from experiments or calculated for a fixed spin-configuration. Areas of potential applications to different magnetic questions are also discussed. The method is finally applied to one situation where the macrospin model breaks down; magnetic switching in ultra strong fields.
An accurate and easily extendable method to deal with lattice dynamics of solids is offered. It is based on first-principles molecular dynamics simulations and provides a consistent way to extract the best possible harmonic - or higher order - potential energy surface at finite temperatures. It is designed to work even for strongly anharmonic systems where the traditional quasiharmonic approximation fails. The accuracy and convergence of the method are controlled in a straightforward way. Excellent agreement of the calculated phonon dispersion relations at finite temperature with experimental results for bcc Li and bcc Zr is demonstrated.
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 ...
A description of non-collinear magnetism in the framework of spin-density functional theory is presented for the exact exchange energy functional which depends explicitly on two-component spinor orbitals. The equations for the effective Kohn-Sham scalar potential and magnetic field are derived within the optimized effective potential (OEP) framework. With the example of a magnetically frustrated Cr monolayer it is shown that the resulting magnetization density exhibits much more non-collinear structure than standard calculations. Furthermore, a time-dependent generalization of the non-collinear OEP method is well suited for an ab-initio description of spin dynamics. We also show that the magnetic moments of solids Fe, Co and Ni are well reproduced.
Details are presented of an efficient formalism for calculating transmission and reflection matrices from first principles in layered materials. Within the framework of spin density functional theory and using tight-binding muffin-tin orbitals, scattering matrices are determined by matching the wave-functions at the boundaries between leads which support well-defined scattering states and the scattering region. The calculation scales linearly with the number of principal layers N in the scattering region and as the cube of the number of atoms H in the lateral supercell. For metallic systems for which the required Brillouin zone sampling decreases as H increases, the final scaling goes as H^2*N. In practice, the efficient basis set allows scattering regions for which H^{2}*N ~ 10^6 to be handled. The method is illustrated for Co/Cu multilayers and single interfaces using large lateral supercells (up to 20x20) to model interface disorder. Because the scattering states are explicitly found, ``channel decomposition of the interface scattering for clean and disordered interfaces can be performed.