Magnetic anisotropy and chirality of frustrated Cr nanostructures on Au(111)
Abstract
By using a fully relativistic embedded cluster Green’s function technique we investigated the magnetic anisotropy properties of four different compact Cr trimers (equilateral triangles) and Cr monolayers deposited on Au(111) surface in both fcc and hcp stackings. For all trimers the magnetic ground state was found a frustrated 120° Néel configuration. Applying global spin rotations to the magnetic ground state, the predictions of an appropriate second order spin Hamiltonian were reproduced with high accuracy by the first principles calculations. For the Cr trimers with adjacent Au atoms in similar geometry we obtained similar values for the inplane and outofplane anisotropy parameters, however, the Dzyaloshinskii–Moriya (DM) interactions appeared to differ remarkably. For two kinds of trimers we found an unconventional magnetic ground state showing 90° intheplane rotation with respect to the high symmetry directions. Due to higher symmetry, the inplane anisotropy term was missing for the monolayers and distinctly different DM interactions were obtained for the different stackings. The chiral degeneracy of the Néel configurations was lifted by less then 2 meV for the trimers, while this value raised up to about 15 meV per 3 Cr atoms for the hcp packed monolayer.
pacs:
75.10.Jm, 75.30.Gw, 75.50.Ee, 75.70.Ak, 75.75.Lf1 Introduction
Accurately describing frustrated magnetic systems both experimentally and theoretically is still a challenge. Development in spinpolarized scanning tunnelling microscopy (SPSTM) made it possible to explore noncollinear magnetic structures in atomic resolution. [1, 2] Frustrated noncollinear magnetic structures were reported by Gao et al [3] for Mn islands deposited on Ag(111) surface. Based on topographic measurements which only collect information from the Mn layer and the topmost Ag layer, they concluded that fcc up and hcp down type islands and fcc stacked stripes were present. They demonstrated that the islands exhibit the 120° Néel magnetic structure and the orientation of the Mn moments differs by 30° between fcc and hcp stacked islands most likely due to spinorbit coupling which is different for the two stackings. The latter phenomenon made it possible to distinguish the islands with different stackings.
The most simple system exhibiting geometric frustration is an antiferromagnetically (AF) coupled symmetric trimer. An equilateral compact chromium trimer deposited on the (111) surface of gold is an archetype of such systems. The first noncollinear magnetic calculations of supported metallic 3d triangular trimers were presented fifteen years ago by Uzdin et al. [4] Within the vector Anderson model they showed that a supported equilateral Cr trimer exhibited zero net magnetic moment with the atomic moments enclosing 120° angles. Ab initio investigations of the compact Cr trimer deposited on Au(111) surface were first performed by Gotsits et al. [5] They performed spinpolarized electronic structure calculation using the projector augmentedwave method with the spinorbit coupling (SOC) included for optimizing the geometric and the magnetic configuration. Starting the optimization from fcc sites, they found that the equilateral geometry was preserved in equilibrium and the average of the magnetization at all the sites were laying in the plane enclosing 120° with each other. The 120° Néel state was also found to be the lowest energy magnetic configuration of the Cr trimer supported on a Au(111) surface by Bergman et al [6] using an extension of the real space linear muffin tin orbital method (RSLMTO) within the atomic sphere approximation (ASA) with relativistic effects included within the scalar relativistic approximation.
Using fully relativistic constrained selfconsistent multiple scattering Green’s function electronic structure calculations [7] and adiabatic spin dynamics [8] to search for the ground state, Stocks et al [9] also concluded that the ground state of the compact Cr trimer on Au(111) is the 120° Néel state. The energy of the two magnetic configurations with opposite chirality turned out to be different as it has also been confirmed by Antal et al. [10] They pointed out that the Dzyaloshinskii–Moriya (DM) interaction is responsible for the lifting of the degeneracy of the states with opposite chirality.
The 120° Néel state of Cr monolayer on Pd(111) substrate was observed by Waśniowska et al [11] by SPSTM measurements. From firstprinciples calculations on flat spin spirals, they predicted the ground state to be the 120° Néel state which is in agreement with the measurements. They found the Cr ML to be energetically favourable in fcc than in hcp stacking by about . Palotás et al [12] demonstrated by first principles calculations that the two possible Néel states with opposite chiralities of a Cr monolayer on Ag(111) are energetically inequivalent. They also showed that the magnetic contrast of the simulated SPSTM image was sensitive to the electronic structure of the tip and to the bias voltage.
The aim of the present work is to systematically investigate the magnetic anisotropy of the four Cr trimers and to compare them to those of Cr monolayers deposited on Au(111). After presenting the geometrical structure of the clusters and monolayers we briefly describe the applied methods used to determine the electronic and magnetic structure of the systems. By exploiting the point group symmetry of the systems, a suitable classical Heisenberg model is constructed which is then applied to give analytic forms of the rotational energies. These forms are used to analyze and discuss the results of the abinitio calculations for the four different trimers and for the monolayers with fcc and hcp stacking. The chirality of the trimers and the monolayers is discussed and the values of the accessible model parameters are given.
2 Details of the calculations
The high symmetry adsorption sites of a fcc(111) surface are the fcc hollow, the hcp hollow, the bridge and the ontop sites. The ontop positions of the Cr on Au(111) is energetically unfavourable as it was pointed out by Gotsits et al [5]. For trimers occupying the bridge positions the bond length would be close to the half of the lattice constant of the underlying supporter which is too short compared to the relaxed bond length predicted by abinitio calculations [5]. Considering these facts, in the present paper we investigate fcc and hcp stacked compact trimers located only at the hollow sites of the substrate.
Equilateral trimers can be deposited in the hollow positions of the fcc(111) surface in four different configurations as it is shown in figure 1. By labelling the inequivalent layers of the fcc lattice along the [111] direction with capital letters, the order of the fcc stacking is ABCABCA, while the order for the hcp stacking is ABCABCB, where the last (boldface) symbol corresponds to the deposited trimers or monolayers. Fcc and hcp stacked trimers both can be either up or down triangles, see figure 1, which can be distinguished by the lateral displacement between the supporter atomic rows and the base of the triangle as it is indicated by the blue double arrows in figure 1. Note, that one cannot distinguish between an fcc up and an hcp down island (fcc down and hcp up) if one only sees the island itself and the topmost supporter layer which is the situation in a constant current nonmagnetic STM experiment.
Previous studies on Cr clusters forming equilateral triangle considered the fcc up trimer [5, 6, 9, 10] and, according to our knowledge, magnetic properties of the other compact Cr trimers have not been investigated yet. The four clusters are, however, different: the Cr atoms in an fcc up or an hcp down cluster surround an interstice in the first supporter layer (breezy triangle) while in an fcc down or an hcp up cluster they surround a substrate atom (crammed triangle). The trimers with fcc and hcp stacking together with the substrate exhibit a point group symmetry where the axis intersects the centre of the triangle normal to the substrate and the reflection planes contain the axis and one of the cluster atom. For the monolayers with fcc and hcp stacking the point group symmetry also holds with axes intersecting the centre of an elementary triangle or a Cr atom. Note that the monolayers contain alternating up and down triangles. The two types of triangles are inequivalent, i.e., one type of triangle can not be transformed into the other type by any of the symmetry operations of the system.
The electronic structure of the Cr monolayers and the Cr trimers on top of Au(111) surface was calculated in terms of the fully relativistic screened KKR method [13, 14] within the local spindensity approximation (LSDA), [15] of the density functional theory as parametrized by Perdew and Zunger [16]. In particular, for the clusters we applied the embedded cluster Green’s function technique [17] based on the KKR method. The effective potentials and fields were treated within the atomic sphere approximation (ASA), and a cutoff of for the angular momentum expansion was used. For the case study presented in this work we neglected the relaxation of the layerlayer distances and the 2D lattice constant of gold ( Å) was applied in all calculations. Correspondingly, in case of the trimers the three Cr atoms occupied hollow positions above the topmost Au layer. In order to let the electron density relax around the cluster the first neighbour shell was also included in the embedded cluster. The magnetic ground state configurations for the trimers were determined selfconsistently by using the procedure described in reference 18. The obtained directions of the magnetic moments at the Cr atoms are depicted in figure 1.
The direction of the exchange field on the three Cr sublattices in the monolayers were fixed according to figure 1(a) in the case of fcc stacking, while according to figure 1(d) in the case of hcp stacking. For very low coverage of Cr on Au(111) surface Boeglin et al reported a spin magnetic moment of for single adatoms while the average spin magnetic moment rapidly vanished with increasing island size indicating antiferromagnetic arrangement, however, the details of the AF order has not been explored. [19] The most probable configuration is a 120° Néel structure, however, spin structures like rowwise or doublerowwise antiferromagnetic alignments [20], or even a threedimensional spin structure [21] have been predicted for 2D triangular lattice. Since our interest is focused on the magnetic anisotropy properties of a 120° Néel structure we haven’t attempted to explore the energetics of the different magnetic configurations. In order to confirm that the 120° Néel structure represents a local minimum of the energy, we investigated the spinexcitation spectrum of the monolayer with both stackings. The details of the calculation is given in the Appendix. The spectrum along the special directions of the first Brillouin zone is depicted in figure 2. From our calculations we obtained no imaginary frequencies indicating that the reference magnetic configuration, i.e., the 120° Néel structure is not a metastable state.
3 Spin model
The magnetic properties of thin films and clusters of transition metals can often be successfully described by a classical Heisenberg model
(1) 
where the first summation runs over the interacting pairs of spins, is a unit vector parallel to the magnetization at site , are generalized exchange interaction matrices and represents the secondorder onsite anisotropy matrix. In case of monolayers with Néel spinstructure are sorted out according to the three sublattices. Without loss of generality, the matrices are chosen to be symmetric and traceless. The exchange matrices can be decomposed into an isotropic part of , a traceless symmetric anisotropic part defined as and an antisymmetric part given by . The latter term is usually formulated with the Dzyaloshinskii–Moriya vector, , as .
We start our study of the magnetic anisotropy by symmetry considerations regarded the spinHamiltonian (1) for the trimers. Since the exchange tensors, , and , as well as the onsite anisotropy matrices, , and , are related in terms of appropriate similarity transformations it is sufficient to explore the structure of one of them. Considering the trimers in figure 1 with respect to the local coordinate systems, and are obviously related to each other as
(2) 
where denotes the reflection to the plane. The anisotropy matrix must be invariant under the same operation. The matrices satisfying the above requirements are parametrized as:
(3) 
and
(4) 
For monolayers, the coupling between the sublattices and , , is defined as the sum of the corresponding exchange tensors: for any . If all the sites are located in the same sublattice, will contribute to the anisotropy matrix, which will be the same for all sublattices. Due to the symmetry of the monolayer, it turns out that the terms and in equation (4) disappear.
Following the method proposed by Szunyogh et al [22] to analyze the magnetic anisotropy of and the energies of the systems are calculated during the simultaneous rotation of the antiferromagnetic configuration around the three fold axis and around an axis parallel to the magnetization on one of the sites or sublattices. In the following this energy will be referred to as rotational energy. Similar procedure has been applied to study the magnetic anisotropy for interface [23] and for a cobalt nanocontact [18].
The two type of 120° Néel structures can be distinguished by investigating the chirality of the configurations. According to Antal et al [10] the chirality vector for a trimer is defined as
(5) 
For inplane spinconfigurations, the vector is parallel to the axis and its component, , being either or characterizes the spinconfigurations. The chiralities of the obtained spinconfigurations for the trimers are also given in figure 1.
In the case of an inplane Néel spin structure of a monolayer, the chirality for the up and down triangles alternates between the values of or . In this case, we associate the chirality of the monolayer with the chirality of the up triangles. The energy of the Néel structures with opposite chiralities may be different due to the nonvanishing component of the DM vector. While rotation around the three fold axis does not alter the chirality of an in plane Néel configuration, a 180° rotation around an axis parallel to the magnetization on one of the sites will reverse its value.
The rotational energies were calculated in the spirit of the magnetic force theorem [24]. The effective potentials and exchange fields determined in ground state configurations were kept fixed and the change in energy of the system with respect to the rotational angle is approximated by the change in band energy:
(6) 
where is the Fermi energy, and stand for the density of states (DOS) and for the integrated DOS, respectively, and indicates the dependence of these quantities on the spin configuration of the system. For the trimers, the band energy was calculated by using the Lloyd’s formula [25]:
(7) 
where and denote the singlesite scattering matrix and the scattering path operator (SPO) for the host confined to the sites in the cluster, , respectively, while denotes the singlesite scattering matrices of the embedded atoms [17] and we omitted a constant shift of the energy not affected by the magnetic configuration. Note that using formula (7), the change in the band energy due to the change of the magnetic configuration in cluster is accounted for the whole system, while the direct integration of the local DOS in equation (6) is always restricted to a given environment of only. To calculate the SPO of the layered host systems, the energy integrations were performed by sampling 16 points on a semicircular path in the upper complex semiplane and 3300 points were used in the irreducible wedge of the surface Brillouin zone (SBZ). We have checked the accuracy of the SBZintegrals by performing the same calculations using 1900 points in the irreducible wedge of the SBZ and a deviation of up to 4 % was found in the resulting model parameters. In the case of the monolayers we used a fine adaptive set ranging from 7320 points in the SBZ at the Fermi energy to 840 points at the bottom of the band and we note that the calculations were carried out with three atoms per unit cell, i.e., the magnetic unit cell of the Néel structure.
Using the parametrization given by equations (3) and (4) simple expressions can be derived for the rotational energies based on the Heisenberg model. For the rotations around the three fold axis the energy has the form of
(8)  
(9) 
where is an energy independent on the angle of rotation and the superscripts indicate whether a configuration with positive or negative chirality is rotated rigidly around the axis. Similarly, when the configuration is rotated around the axis parallel to the magnetization at the 2^{nd} Cr atom or at the 2^{nd} sublattice the energy can be given as:
(10) 
The superscript indicates here the chirality at . We note again that the rotation around the axis reverses the chirality but the rotation around the axis does not alter the chirality.
By comparing these rotational energy functions to the results of the first principles calculations the values of the coefficients of the trigonometric functions can be extracted. From equation (3) the energy difference between the positive and the negative chirality configuration can be read off: .
Since the relative angle between the spins does not change during the global spin rotations, the contribution of the isotropic exchange cancel out and the parameter of the model is not accessible through the and functions. Regarding the focus of the recent work, the factual value of the isotropic exchange is, therefore, irrelevant. Furthermore, we note that the rotational invariant fourth order terms introduced by Antal et al [10] neither contribute to the rotational energy.
4 Trimers
In the case of the trimers the spin moments of the Cr atoms scattered between and while the orbital moments between and depending on the geometry and the magnetic structure. These values are in good agreement with the results of previous studies: if geometrical relaxation is included [5]; [6] and [10] if it is neglected. For the orbital moment also small values have been reported: with geometrical relaxation included [5] and [10] without relaxation.
The rotational energies are calculated with a resolution of 3° for all four trimers and the results are shown in figure 3. The parameters in equations (9) and (3) are obtained as the Fourier components of the rotational energies and listed in table 1. Using these parameters the functions given by equations (9) and (3) fit with a high accuracy to the results provided by the abinitio calculations.

Trimer/ML fcc up fcc down hcp up hcp down fcc ML hcp ML
From figure 3 it can be inferred that the energy minimum corresponds to for the fcc up and hcp down trimers, while to for the fcc down and hcp up trimers, see also figure 1. Note, however, that in case of the fcc up trimer the energy difference between the two chiral states is found to be eV which is near the computational accuracy of our method.
In the case of the , the relative orientation of the magnetization vector and the easy direction set by the onsite anisotropy term () is the same for the three Cr atoms and this situation is preserved during the global inplane rotation, therefore, the anisotropy energies of the single atoms are simply summed up. The same argument holds for the twosite anisotropies (), thus, as indicated by equation (9), we expect a angular dependence for the inplane rotational energy. This is clearly confirmed by the first principles calculations, see the lower graph of figure 3.
For the inplane rotational energy of the trimers we expect the anisotropy terms to cancel since the second order inplane anisotropy energies are sampled at angles , , and , for which , i.e., independent of the angle of rotation. The magnitude of the rotational energies of the positive chirality trimers was indeed found below 7 eV, indicating a very small deviation between the spin model (1) and the ab initio calculation. Similarly, Szunyogh et al [22] found a angular dependence of the rotational energy of with an amplitude of in the socalled state with negative chirality, while for the states with positive chirality the rotational energy had no angular dependence up to an absolute error of .
It can be read off from table 1 that the trimers with similar environment, i.e., the breezy and the crammed triangles, exhibit similar parameter values. This is, in particular, valid for the outofplane and inplane anisotropy parameters, and , respectively. The component of the DM vector turned out to be similar for the fcc down and the hcp up (crammed) trimers, but for the fcc up and the hcp down (breezy) trimers are rather different. We notice that for an fcc up trimer Antal et al [10] reported a value of which is a remarkable difference compared to our present value of . There are, however, distinct differences between the two calculations. On the one hand, here we included one shell of environment around the atoms forming the trimer, whereas in reference 10 only the Cr atoms were taken into account in the selfconsistent calculations. On the the other hand, we calculated the rotational energies in terms of the Lloyd’s formula, equation (7), while Antal et al [10] used equation (6) to evaluate the bandenergy.
In case of the fcc up trimer, we repeated the magnetic force theorem calculation of the rotational energies by using the selfconsistent effective potentials and fields from the negative chirality configuration and found that even in this case the positive chirality state was lower in energy. Remarkably, however, Stocks et al [9] found considerable difference for if they used the (negative chirality) Néel state () or the outofplane ferromagnetic state () for the selfconsistent reference potential and field calculation.
Regarding the inplane anisotropy (see the lower graph of figure 3), for the fcc up and hcp down trimers we found a value of which is about 50 % larger in magnitude as compared to the fcc up trimer calculations of Stocks et al [9]. The reason for this difference is the same as mentioned above in context to . In the case of fcc down and hcp up trimers, the reversed (positive) value of means that the ground state of these trimers is rotated by 90° with respect to the conventional, high symmetry directions of the Néel state, see figs. 1(b) and 1(c). Similarly, Gao et al [3] found that triangular Mn islands of different stackings exhibit different easy directions inside the 120° Néel structure.
5 Monolayers
For both fcc and hcp stacked monolayers we obtained a spin magnetic moment of and an orbital magnetic moment of for the Cr atoms. The band energies while rotating the magnetic configuration around the axis lying in the plane of the monolayer are shown in figure 4. Using the parameters in table 1 the results of the first principles calculations can be fitted with a high accuracy by the function given in equation (3). Due to the fact that each site in the monolayer forms a symmetry centre, the energy of the system contains at best fourth order terms () in case of inplane global rotations. Correspondingly, the band energy turned out to be practically independent on the angle of rotation around the axis, . Note that Szunyogh et al found nonvanishing inplane anisotropy parameters for the (111) layers in bulk [22] or at the interface [23] since for these systems the above symmetry does not apply.
Fitting the outofplane rotational energy in figure 4 to equation (3), we obtain nearly the same parameters for the fcc and hcp monolayers. These parameters are somewhat reduced in magnitude as compared to those for the trimers. As obvious from the nearly like dependence of the bandenergy curves in figure 4, the outofplane rotational energies are dominated by the normaltoplane components of the DM interactions and is opposite in sign for the fcc and the hcp monolayers. Note that for the monolayer case is the only interaction which distinguishes between the two Néel states with opposite chiralities.
The parameters are almost an order larger in magnitude for the monolayers than for the trimers, see table 1. This can be understood due to the following reasoning. The main contribution to the DM interactions is due to the nearest Cr neighbours. Since in the case of monolayers the number of nearest neighbours is three times larger then those in the trimers, a corresponding enhancement of is expected. More quantitatively, the magnetic unit cell of the monolayers is composed of 3 up and 3 down elementary trimers. Hence, a first estimation of the energy difference between the two chirality states of the monolayer could be 3 times the sum of the energy differences of the trimers. (It should be recalled that the chirality index of the up and down trimers are opposite in a monolayer, therefore, the chiral energy of the down trimer should be subtracted from that of the up trimer.) From the data of table 1 we calculate for the fcc monolayer and 5.118 for the hcp monolayer. The large deviation of these values from 3 indicates that the spininteractions in a monolayer are rather different from those in the trimers and/or interactions between more distant pairs have important contributions.
6 Conclusions
With regard to the 120° Néel structure, we systematically investigated the inplane, outofplane and chiral magnetic anisotropy energies of Cr trimers in four different geometries and also of fcc and hcp stacked Cr monolayers deposited on the Au(111) surface. Only one out of the four geometric positions for the trimer was considered earlier and here we showed that the DM interactions depend intriguingly on the geometry. The magnetic ground state of the systems turned out to be formed due to an interplay between the DM interactions and the twosite anisotropy. The actual values of the corresponding parameters determine the energy barrier between the local energy minima related to the magnetic states with different chiral indices. Moreover, we revealed an unconventional inplane easy axis in fcc down and in hcp up trimers. It should be noted that, in terms of SPSTM experiments, the different magnetic states of otherwise indistinguishable fcc and hcp stacked Mn overlayers on Ag(111) became identifiable [3]. The theoretical investigation of the underlying phenomena might thus gain considerable attention.
Appendix A Magnetic excitation in a 3sublattice antiferromagnet
In order to derive the spinwave Hamiltonian for the Néel antiferromagnetic configuration we will follow the method detailed in reference 26. The time evolution of the magnetic moments is given by the LandauLifshitz equation. [27] In a local coordinate system where the direction of the magnetization in sublattice , , and two transverse unit vectors, and , form a righthand system it has the following form [28]:
(11)  
(12) 
where is the band energy of the system, is the magnetization and and are the angle of rotation around the transverse directions and at site of the sublattice . Introducing the variables
(13) 
and expanding the energy up to second order around the Néel state the LandauLifshitz equation has the form of [26]:
(14)  
(15) 
where the , and matrices are related to the second derivatives of the band energy with respect to the transverse change of the exchange field:
(16) 
The analytic formulas for the second derivatives can be found in reference 18. After applying a lattice Fourier transform in space and a continuous Fourier transform in time the magnon frequencies will be the solutions of the following eigenvalue equation:
(17) 
where the , , and matrices are the Fourier transform of the corresponding second derivatives, e.g.,
(18) 
with running over the lattice sites of the th sublattice.
References
 [1] M. Bode, M. Heide, K. von Bergmann, P. Ferriani, S. Heinze, G. Bihlmayer, A. Kubetzka, O. Pietzsch, S. Blügel, and R. Wiesendanger. Chiral magnetic order at surfaces driven by inversion asymmetry. Nature, 447:190–193, 2007.
 [2] Roland Wiesendanger. Spin mapping at the nanoscale and atomic scale. Rev. Mod. Phys., 81:1495–1550, Nov 2009.
 [3] C. L. Gao, W. Wulfhekel, and J. Kirschner. Revealing the 120° Antiferromagnetic Néel Structure in Real Space: One Monolayer Mn on Ag(111). Phys. Rev. Lett., 101:267205, Dec 2008.
 [4] S. Uzdin, V. Uzdin, and C. Demangeat. Magnetic trimer on nonmagnetic substrate: From frustration towards noncollinearity. EPL (Europhysics Letters), 47(5):556, 1999.
 [5] H. J. Gotsis, Nicholas Kioussis, and D. A. Papaconstantopoulos. Evolution of magnetism of Cr nanoclusters on Au(111): Firstprinciples electronic structure calculations. Phys. Rev. B, 73:014436, Jan 2006.
 [6] Anders Bergman, Lars Nordström, Angela Burlamaqui Klautau, Sonia FrotaPessôa, and Olle Eriksson. A firstprinciples study of the magnetism and electronic structure of Cr clusters supported on a Au(111) surface. Journal of Physics: Condensed Matter, 19(15):156226, 2007.
 [7] B. Újfalussy, XinDong Wang, D. M. C. Nicholson, W. A. Shelton, G. M. Stocks, Yang Wang, and B. L. Gyorffy. Constrained density functional theory for first principles spin dynamics. Journal of Applied Physics, 85(8):4824–4826, 1999.
 [8] V. P. Antropov, M. I. Katsnelson, M. van Schilfgaarde, and B. N. Harmon. Ab Initio Spin Dynamics in Magnets. Phys. Rev. Lett., 75:729–732, Jul 1995.
 [9] G. Malcolm Stocks, M. Eisenbach, B. Újfalussy, B. Lazarovits, L. Szunyogh, and P. Weinberger. On calculating the magnetic state of nanostructures. Progress in Materials Science, 52(2–3):371–387, 2007.
 [10] A. Antal, B. Lazarovits, L. Udvardi, L. Szunyogh, B. Újfalussy, and P. Weinberger. Firstprinciples calculations of spin interactions and the magnetic ground states of Cr trimers on Au(111). Phys. Rev. B, 77:174429, May 2008.
 [11] M. Waśniowska, S. Schröder, P. Ferriani, and S. Heinze. Real space observation of spin frustration in Cr on a triangular lattice. Phys. Rev. B, 82:012402, Jul 2010.
 [12] Krisztián Palotás, Werner A. Hofer, and László Szunyogh. Simulation of spinpolarized scanning tunneling microscopy on complex magnetic surfaces: Case of a Cr monolayer on Ag(111). Phys. Rev. B, 84:174428, Nov 2011.
 [13] L. Szunyogh, B. Újfalussy, P. Weinberger, and J. Kollár. Selfconsistent localized KKR scheme for surfaces and interfaces. Phys. Rev. B, 49:2721–2729, Jan 1994.
 [14] L Szunyogh, B Újfalussy, P Weinberger, and J Kollár. The selfconsistent fully relativistic SKKR Green function method: applications to the (100), (110) and (111) surfaces of Au and Pt. Journal of Physics: Condensed Matter, 6(18):3301, 1994.
 [15] D. M. Ceperley and B. J. Alder. Ground State of the Electron Gas by a Stochastic Method. Phys. Rev. Lett., 45:566–569, Aug 1980.
 [16] J. P. Perdew and Alex Zunger. Selfinteraction correction to densityfunctional approximations for manyelectron systems. Phys. Rev. B, 23:5048–5079, May 1981.
 [17] B. Lazarovits, L. Szunyogh, and P. Weinberger. Fully relativistic calculation of magnetic properties of Fe, Co, and Ni adclusters on Ag(100). Phys. Rev. B, 65:104441, Mar 2002.
 [18] László Balogh, Krisztián Palotás, László Udvardi, László Szunyogh, and Ulrich Nowak. Theoretical study of magnetic domain walls through a cobalt nanocontact. Phys. Rev. B, 86:024406, Jul 2012.
 [19] C. Boeglin, P. Ohresser, R. Decker, H. Bulou, F. Scheurer, I. Chado, S. S. Dhesi, E. Gaudry, and B. Lazarovits. Spin magnetic moments from single atoms to small Cr clusters. physica status solidi (b), 242(9):1775–1778, 2005.
 [20] P. Krüger, M. Taguchi, and S. MezaAguilar. Magnetism of transitionmetal monolayers on Cu(111) and Ag(111). Phys. Rev. B, 61:15277–15283, Jun 2000.
 [21] Ph. Kurz, G. Bihlmayer, K. Hirai, and S. Blügel. ThreeDimensional Spin Structure on a TwoDimensional Lattice: MnCu(111). Phys. Rev. Lett., 86:1106–1109, Feb 2001.
 [22] L. Szunyogh, B. Lazarovits, L. Udvardi, J. Jackson, and U. Nowak. Giant magnetic anisotropy of the bulk antiferromagnets and from first principles. Phys. Rev. B, 79:020403, Jan 2009.
 [23] L. Szunyogh, L. Udvardi, J. Jackson, U. Nowak, and R. Chantrell. Atomistic spin model based on a spincluster expansion technique: Application to the interface. Phys. Rev. B, 83:024401, Jan 2011.
 [24] H. J. F. Jansen. Magnetic anisotropy in densityfunctional theory. Phys. Rev. B, 59:4699–4707, Feb 1999.
 [25] P Lloyd. Wave propagation through an assembly of spheres: II. The density of singleparticle eigenstates. Proceedings of the Physical Society, 90(1):207, 1967.
 [26] L. Udvardi, L. Szunyogh, K. Palotás, and P. Weinberger. Firstprinciples relativistic study of spin waves in thin magnetic films. Phys. Rev. B, 68:104436, Sep 2003.
 [27] L. Landau and E. Lifshitz. On the theory of the dispersion of magnetic permeability in ferromagnetic bodies. Physik. Zeits. Sowjetunion, 8:153–169, 1935.
 [28] L. Rózsa, L. Udvardi, and L. Szunyogh. Langevin spin dynamics based on ab initio calculations: numerical schemes and applications. ArXiv eprints, January 2014. arXiv:1401.7885v1 [condmat.mtrlsci].