Spin-orbit induced non-collinear spin structure in deposited transition metal clusters
The influence of the spin-orbit coupling on the magnetic structure of deposited transition metal nanostructure systems has been studied by fully relativistic electronic structure calculations. The interplay of exchange coupling and magnetic anisotropy was monitored by studying the corresponding magnetic torque calculated within ab-initio and model approaches. It is found that a spin-orbit induced Dzyaloshinski-Moriya interaction can stabilise a non-collinear spin structure even if there is a pronounced isotropic ferromagnetic exchange interaction between the magnetic atoms.
pacs:Valid PACS appear here
As the complexity and technological applicability of nanostructured magnetic materials grows it has become important to develop a reliable quantitative theoretical framework in which to understand them. In principle this is available from relativistic density functional theory (RDFT) (1). The magnetic behaviour of complex magnetic systems is described effectively in terms of ’local moments’, even in metallic systems, which fluctuate as the temperature is increased from 0K. Within the effective single electron picture of DFT, the moments arise from local magnetic fields associated with the different atomic sites which affect the electronic motions and are self-consistently maintained by them. The spin-orbit coupling (SOC) effects on the electronic structure determine the magnetic anisotropy linking the magnetic and spatial structure of a material. As demonstrated by Staunton et al. (2) and others, magnetic torque calculations enable magnetic anisotropy to be studied reliably and magnetic structures to be determined. We have recently developed a method to study ab-initio the magnetic structures of complex nanostructures using this approach (3). It is based on the framework of relativistic density functional theory using the local spin density approximation (LSDA) for exchange and correlation effects (4). The electronic structure was determined in a fully relativistic way from the basis of the Dirac equation for spin-polarized potentials which is solved using the Korringa-Kohn-Rostoker (KKR) multiple scattering formalism (5).
In this letter we set out to understand the relativistic effects on the magnetic structure of nanoclusters in simpler terms of a model ‘moment’ Hamiltonian. In this context the mapping of energetic properties obtained from first principles calculations of complex magnetic systems onto a Heisenberg Hamiltonian has proved over recent years to be a very robust and successful scheme (6); (7). Moreover for many systems the results concur with those from an ab-initio ‘disordered local moment’ theory in which no prior mapping to a Heisenberg system is assumed (8); (9). An extensively used approach to calculate the isotropic exchange interaction parameter for two magnetic moments on sites and for use in the classical Heisenberg Hamiltonian was worked out by Lichtenstein et al. (6) using perturbation theory and the so-called Lloyd formula. A corresponding fully relativistic approach was later introduced by Udvardi et al. (7) that produces an exchange interaction tensor for use in a prescribed extended Heisenberg Hamiltonian. This scheme generates in particular a Dzyaloshinski-Moriya (DM)-type interaction (10); (11), which may explain many interesting phenomena including the magnetic ground state configuration of nanostructures (12); (13) as well as magnetic thin films (14); (15). Here we explore what form the effective moment Hamiltonian should have by a detailed study of its building blocks. With an investigation of Fe, Co and Ni dimers deposited on Pt(111) as most simple cluster examples we identify strong DM-type interactions as well as an additional substrate-generated effect in Ni (for calculational details see (3)). The conclusions of these results can then be applied straightforwardly to larger clusters and nanostructures.
We start by considering the magnetic torque vector acting on an atomic magnetic moment on a site ( or ) and aligned along direction . The torque vector is defined in terms of the change in energy of the system when changing the orientation of the magnetic moment, on site , . The component of with respect to the axis can be determined from first-principles using an expression derived by Staunton et al. (2). Following on from this the derivative describing the change in energy upon changing the orientation of two magnetic moments on sites and can also be obtained (6); (7). By making use of the rigid spin approximation (RSA) (16) this approach leads to a fitting of the magnetic energy landscape at low temperatures obtained from first principles calculations onto a Heisenberg model Hamiltonian. With SOC included an anisotropy in the exchange interaction may occur. We then use these quantities in an extended classical Heisenberg Hamiltonian for the ‘spins’, of the following form, (e.g.(17); (7)):
Here the exchange interaction tensor has been split into its conventional isotropic part , its traceless symmetric part and its anti-symmetric part . We assume that the latter one is represented in terms of the DM vector with ( is the Levi-Civita symbol). Finally, the anisotropy constants account for the so-called on-site anisotropy energy associated with each individual moment oriented along . From our ab-initio calculation of the torques we can test whether such a model Hamiltonian is justified and also find the values of the , , , and parameters.
For the model Eq. (1) the rate of change in energy when a moment on site is rotated about an axis can be partitioned into contributions from different terms of Heisenberg Hamiltonian . The contribution to the torque from the DM coupling is given by:
The last term derives from the single site anisotropy term, i.e. . By focussing on nanoclusters of collinear ‘spin’ arrangements, the magnetically anisotropic terms can be determined. Note that the last term in Eq. (2) does not contribute to the torque in the case of a collinear magnetic structure. Also the sum of the DM contributions from all sites in a nanocluster to the total torque vanishes in this case. On the other hand, the anisotropy of the exchange interaction, represented by the symmetric tensor gives rise to a finite torque hence contributing to the total magnetic anisotropy energy of the nanocluster.
We now see how well our ab-initio calculations of deposited Fe, Co and Ni dimers, using our first-principles multiple scattering formalism, fit such model torques. Firstly we find that the dimers, when constrained to be magnetically collinear exhibit a pronounced out-of-plane magnetic anisotropy (18). The collinear magnetic moments and have been assumed to be orientated at an angle to the surface normal (-axis) and are expressed in terms of polar angles and , as . Fig. 1 shows the atomic configuration together the projection of the moments onto the surface (-plane). In parallel to the fixed frame of reference () we use a second one () rotated by with respect to the fixed one with (see Fig. 1). In all calculations the torque is taken around the axis, i.e. , for fixed to as a function of the azimuthal angle . From the model Eq. (1) an expansion of the anisotropy energy term in terms of spherical harmonics, gives the single site contribution to the torque in present configuration to be whereas deduced from the model has a () variation when lies in the () plane.
In Fig. 2 the ab-initio results for the torque are compared with those deduced from the model Eq. (1). As one notes, the torques for the two atoms are different but are related with respect to their -dependence according to the C symmetry of the system. On comparing the torque of the two Co atoms of Co/Pt(111) obtained directly from the electronic structure calculations with those from the model (see Fig. 2a-c) we find the symmetric part of the exchange interaction tensor to be negligible. The contribution due to the DM interaction Eq. (2) is shown in Fig. 2b. Clearly, and vary with and are opposite in sign. In accordance with the C symmetry, the DM vector lies in the plane (Fig. 1). For the contribution due to the on-site anisotropy one finds to be very small and the dominating terms and to be practically the same for both atomic sites leading to .
The contribution to connected with does not depend on , while that connected with varies with (Fig. 2c). As one can see, the torque derived from the model Hamiltonian reproduces the results calculated directly rather well. The remaining deviation is primarily due to the limitations of the model Hamiltonian with respect to the dependency of magnetic energy on the magnetic moment orientations . From the decomposition of the ’s via the model Hamiltonian it becomes clear that its -dependence is dominated by the DM-contribution while the -contribution gives rise to a minor additional modulation. Owing to the large negative values of and an out-of-plane anisotropy results for the total system. This can be seen also from Fig. 2e where the total torque for Co on Pt(111) is shown together with the individual contributions and the corresponding DM terms .
Once more one sees that the -dependence of the individual torques is determined by while that of the total torque is set by the on-site anisotropy terms. This also holds for Fe (Fig. 2d) for which the DM terms are even more dominant, i.e. the -dependence of the individual torques are nearly exclusively due to .
The situation is noticeably different for Ni on Pt(111) (Fig. 2f) for which the DM terms give only minor contributions to the individual torques and . In contrast to Fe and Co the difference between the torques of different atoms cannot be attributed to DM coupling. The period of oscillations of the torques as functions of are different compared to those of the torques created by DM coupling varying as . Clearly an effect which is not contained in the model Hamiltonian Eq. (1)) is evident and must derive from interactions of the Ni moments with those induced in the Pt substrate. As with a DM-type interaction it produces torques such that but its variation would arise if a term of the form were added to Eq. (1), where lies along and along . These Ni dimer results indicate that an effective Heisenberg model must be used with caution for systems where the magnetic structure of a nanocluster is strongly influenced by the spin polarisability of the substrate.
For the Co and Fe dimers, however, the model works very well. For our chosen geometry and and using the symmetry properties of the elements of the exchange tensor, one finds for Co/Pt(111) with the total torque . Finally, the magnetic anisotropy energy (MAE) of the dimer being the difference in energy when the magnetic moments are both oriented along and , is given by the integral (3). Obviously this has no contribution from the DM interaction. For Co/Pt(111) we find the exchange parameters , and to be nearly the same implying that the total MAE of the dimer is nearly exclusively due to the on-site anisotropy. The values meV and meV for Co/Pt(111) lead as mentioned above to a pronounced out-of plane anisotropy, i.e. in the ground state the total magnetisation points along the surface normal.
Taking the difference between the individual torques one arrives at the relation allowing to be deduced directly from the ab-initio torque calculations. Table 1 shows the corresponding results for dimers Fe and Co on Pt(111) in comparison with data derived from a mapping to our model Hamiltonian, Eq. (1).
The closeness of the results justifies once more the use of the model Hamiltonian for Fe and Co. In addition one notes that has an appreciable value compared to the isotropic exchange constant . Fixing the azimuthal angle to be and performing similar steps one finds to be . This is also in line with the symmetry of the investigated dimer systems (see Fig. 1). , on the other hand, may take a non-zero value and is found to be comparable to (see Table 1).
Thus, the above analysis shows that the torques and may differ even if the total torque is zero, i.e. if the moments are aligned collinearly along the easy axis (surface normal). The difference between these torques is caused exclusively by the terms leading to a rotation around the -axis. Minimizing the magnetic energy of the 2 atom clusters leads to an outward tilting of the moments by an angle given by atan with being the isotropic exchange coupling constant. The corresponding results given in Table 1 show that the DM interaction causes the deposited Fe and Co dimers to have an appreciable deviation from collinear configurations in spite of the pronounced ferromagnetic exchange coupling given by . This effect of SOC is completely in line with the findings of Sandratskii and Kübler (19) for bulk systems.
The substrate clearly plays a crucial role in the DM interaction. Firstly, hybridisation with the substrate breaks the inversion symmetry for the dimer leading to a non-zero DM vector. This symmetry effect is also confirmed by test calculations on free dimers and on dimers embedded in a bulk Pt matrix. Secondly, the hybridisation with the substrate also allows the SOC effects of the substrate to be transferred to the magnetic 3d transition metal dimer. We have confirmed this by further calculations in which the SOC of the substrate and the dimer atoms were manipulated separately. Enhancing the SOC for a Co dimer leads primarily to an increase of the on-site anisotropy . However, enhancing the SOC for the Pt substrate leads to a strong increase in the anisotropy as well as to a larger difference in the individual torques on the two Co atoms, reflecting an increase of the DM interaction (Fig. 2). This behaviour is in line with Levi’s model of the indirect DM interaction between two spin moments (20), which is mediated by nearby atoms. As a consequence, the magnitude of the DM interaction is essentially determined by the SOC strength of the neighbouring atoms.
The effect of anisotropic exchange is even more spectacular in magnetic alloy nanoclusters, where magnetic atoms are separated by non-magnetic atoms with large SOC. This is demonstrated by our calculations for a FePt () 2D alloy cluster deposited on the Pt(111) surface pertinent to an experimental investigation by Honolka et al. (21). Ab-initio torque calculations show the clear trend that the magnetic structure of this cluster is non-collinear. Here we present the results of Monte Carlo simulation based on an effective Heisenberg model with the isotropic and DM exchange coupling parameters calculated in accordance to Ref. (7). Fig. 3 shows the non-collinear magnetic structure obtained for K. The non-collinearity between the Fe chains is essentially caused by the nearest neighbours Fe-Fe interchain DM interaction ( meV) being of similar magnitude when compared to the isotropic exchange interaction ( meV). Within a Fe chain, however, the DM interaction is more than one order of magnitude smaller when compared to the isotropic exchange leading only to a slight screwing of the Fe magnetic moments along the chain. Thus, we see that the FePt cluster’s non-collinear structure is created by an enhanced DM coupling between Fe moments mediated by Pt atoms having large SOC, on the one hand side, while the isotropic exchange between these atoms separated by non-magnetic Pt is small enough to make both these couplings comparable.
In summary our investigation of deposited transition metal dimers as simple but realistic model systems has demonstrated that ab-initio magnetic torque calculations enable the impact of SOC on the magnetic interactions within nanostructures to be monitored in a very detailed way revealing subtle anisotropic effects. The analysis of directly calculated electronic structure quantities within a framework based on a Heisenberg model Hamiltonian gives further insight and identifies the role of the various contributions and also the limitations of such models. For the Fe and Co dimer systems studied here the DM interaction was found to be pronounced owing primarily to the SOC of the substrate. Moreover it leads to non-collinear magnetic configurations of the dimers in spite of the pronounced ferromagnetic coupling and out-of-plane anisotropy. These magnetically anisotropic interactions also have a profound effect in larger clusters and we have demonstrated this explicitly with a study of a deposited FePt cluster. In particular we infer that the magnetic structure around the edges of magnetic nanoparticles is likely to be significantly affected.
Financial support by the Deutsche Forshungsgemeinschaft within the framework of the priority program (DFG-Schwerpunktprogramm 1153) Clusters at Surfaces: Electron Structure and Magnetism is gratefully acknowledged.
- preprint: APS/123-QED
- A. H. MacDonald and S. H. Vosko, J. Phys. C: Solid State Phys. 12, 2977 (1979).
- J. B. Staunton et al., Phys. Rev. B 74, 144411 (2006).
- S. Bornemann et al., Eur. Phys. J. D 45, 529 (2007).
- S. H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 58, 1200 (1980).
- H. Ebert, in Electronic Structure and Physical Properties of Solids, Vol. 535 of Lecture Notes in Physics, edited by H. Dreyssé (Springer, Berlin, 2000), p. 191.
- A. I. Liechtenstein, M. I. Katsnelson, V. P. Antropov, and V. A. Gubanov, J. Magn. Magn. Materials 67, 65 (1987).
- L. Udvardi, L. Szunyogh, K. Palotas, and P. Weinberger, Phys. Rev. B 68, 104436 (2003).
- B. L. Györffy, A. J. P. ans J Staunton, G. M. Stocks, and H. Winter, J. Phys. F: Met. Phys. 15, 1337 (1985).
- J. Staunton, S. Ostanin, S. A. Razee, and B. L. Györffy, Phys. Rev. Letters 93, 257204 (2004).
- I. Dzyaloshinsky, J. Phys. Chem. Solids 4, 241 (1958).
- T. Moriya, Phys. Rev. 120, 91 (1960).
- A. Anta et al., Phys. Rev. B 77, 174429 (2008).
- E. Y. Vedmedenko, L. Udvardi, P. Weinberger, and R. Wiesendanger, Phys. Rev. B 75, 104431 (2007).
- M. Bode et al., 447, 190 (2007).
- P. Ferriani et al., Phys. Rev. Letters 101, 027201 (2008).
- V. P. Antropov et al., Phys. Rev. B 54, 1019 (1996).
- V. P. Antropov, M. I. Katsnelson, and A. I. Liechtenstein, Physica B 237-238, 336 (1997).
- P. Gambardella et al., Science 300, 1130 (2003).
- L. M. Sandratskii, Adv. Phys. 47, 91 (1998).
- P. M. Levi, Solid State Commun. 7, 1813 (1969).
- J. Honolka et al., Phys. Rev. Letters 102, 067207 (2009).