Arxiv:1506.03735V1 [Cond-Mat.Mtrl-Sci]

arXiv:1506.03735v1 [cond-mat.mtrl-sci] 11 Jun 2015 oantcsaemmr lo ic t icvr [ discovery fer- its understood since best alloy and shape-memory studied romagnetic most the been has which tde aeila ti ietyrltdt Ni to related directly is it as material studied ese ytm,Mn systems, Heusler ancestor and Mn material is pioneering which famous of the pounds, udmna krinrltdrsac ri h field [ alloys the for shape-memory in as magnetic or well of research as skyrmion-related applications, fundamental ac- (STT-RAM) torque–random memory transfer cess spin for anisotropy uniaxial nrltdsses Mn systems. done been related has research on considerable Thus, materials. ing o hti hoeial ikdt Ni to linked theoretically is that loy Mn Mn Mn Recently, material. re- the a a is of As magnetization property of anisotropy. desired orientation strain-induced perpendicular the or sult, shape- rather by properties intrinsic than by controlled best is (ultra-)thin structures in (PMA) anisotropy magnetic perpendicular [ measured Mn of anisotropy talline Mn of [ family systems the familiar are ap- compounds of relatives is relatives These of task of the means class Ad- by structures, proached the Heusler-like practicality. and withing of Heusler the with request proof within or this concepts applica- production dressing device and mass or of functionality devices scope the modern improve of bility to needed are entertclytetd uha Mn as have such or treated, characterized theoretically and synthesized been been have alloys et nadto oMn con- to manganese addition increasing In ferri- the to tent. to ferro- due from ordering transition magnetic a with series, substitution osiuetefudto o antcrctakmem- by racetrack could proposed magnetic these as for ory as foundation such the materials constitute anisotropic that thought h pnrnc omnt ead aeil with materials demands community spintronics The 2 2 ea[ FeGa al- shape-memory ferrimagnetic a is however, NiGa, 7 16 eea ie ndffrn cain.I is It occasions. different on times several ] oet,wihidctsta hs aeil r promisin are materials these applications. that magnetization-switching transition a indicates heavy anisotropy which containing magnetocrystalline moments, compounds the highlighte micro Heusler are on the tetragonal attention temperatures, reveal our changes Curie structural dev focus the electronic they in the force where trends states, driving especially the ou The properties, as single magnetic well of We tetragonal. the set become on detail. large likely effect in a most studied that from is those Starting distortions tetragonal presented. to is theory functional ,adMn and ], 2 opeesv td ftettleeg fmanganese-rich of energy total the of study comprehensive A ntttfu et¨repyi,Tcnsh Universitä für Festkörperphysik, Technische Institut 3 .INTRODUCTION I. a[ Ga ants nttaoa agns-ihHulrcompounds Heusler manganese-rich tetragonal in Magnetism akne al. et Parkin 2 4 ua Wollmann, Lukas Z Y , 2 3 3 5 ia te Mn other NiGa, a[ Ga ahsbe acltd[ calculated been has Ga 2 .Teuixa magnetocrys- uniaxial The ]. eercniee spromis- as reconsidered were , iai atclri well- a is particular in NiGa 4 , 1 a-lnkIsiu ¨rCeicePyi etrStoffe, fester Physik für Chemische Max-Planck-Institut 7 ¨ tnte tas 0 18 rse,Grayand Germany Dresden, 01187 40, Strasse Nöthnitzer .I diint ma- to addition In ]. [ 8 2 .Teindispensable The ]. 1 bsdHulrcom- Heusler -based .Sc materials Such ]. 2 natruha through MnGa 1 2 2 tnsa Chadov, Stanislav bsdHeusler -based oa[ CoGa 2 2 11 MnGa, -based 6 9 and ] 2 – , , 10 15 3 ], ]. ]. 1 naporaecontext. appropriate an w tutr ye,teso-called the types, structure two esrno cuaino ie) Mn sites), of occupation random less Cu verse Mn as nti ae,w hwhwtemgeimand magnetism the Mn how withinthe show change we structure atomic paper, the this In aforemen- the for applications. formation to materials tioned intend the tetragonal We governing of magnetism trends systems. and general of the set comprehen- selected understand i.e., a study, current studying the sively in approach similar a 3 which in terials hssaaoost Cu cubic to were Heusler analogous Fritz of phases name the with associated first that intention same [ The work cre- previous concept our general guided insight. a in- comprehensive into the im- knowledge ates knowledge, an cluster valuable is of in material results clusion single that a task of portant study detailed effects. a bias Although exchange large of context iodrd,Mn disordered), eune fti itrin .. h magnetocrystalline the [ i.e., (MCA) anisotropy distortion, con- this the as of well as sequences distortion, tetragonal the causing tities X cuaino h ykffpstos 8 positions, Wyckoff with the compound, These original of the structures. occupation materials from similar derived Heusler of are of variety structures a family incorporates the that to of compounds definition original Heusler the extending compound”, “Heusler Mn d Jürgen Kübler, 2 ese losaenwdy nomlydvddinto divided informally nowadays are alloys Heusler Y 2 Z Y amtd,629Drsat Germany Darmstadt, 64289 Darmstadt, t eis,euvln ytm ihhairconstituents heavier with systems equivalent series), 2 2 nla h eeec ytm[ system reference the as MnAl ua[ RuGa Y cpcoii fteosre rns We trends. observed the of origin scopic pce aebe netgtd mn hs were these Among investigated. been have species ye,rfrigt h rgnlHulrcompound Heusler original the to referring types, (4 rmte n iia aeil eetu labeled thus were materials similar on, then From . d h ytm htrmi ui from cubic remain that systems the t easacmaidb o magnetic low by accompanied metals dfidatnsigyhg ausfor values high astonishingly find nd ) addtsfrsi-rnfrtorque spin-transfer for candidates g aefo h ltrPuigrl,as rule, Slater–Pauling the from iate .B en fprildniisof densities partial of means By d. a n Mn and Ga, ui aetsses h response the systems, parent cubic ftettaoa itrinadits and distortion tetragonal the of I RSA STRUCTURE CRYSTAL II. 17 ese opud sn density using compounds Heusler wihhsbe on ohv oeor more have to found been has (which ] 2 2 ta[ PtGa Y lui Felser Claudia 3 sa tmfo eidI o the (or IV period from atom an is ,wl ehglgtdadpae into placed and highlighted be will ], 2 Y 19 2 (5 nl ihtestoichiometry the with MnAl, , d 15 ) 20 asre.Tetigrquan- trigger The series. Ga oiae st undertake to us motivated ] ,adMn and ], 1 22 2 ha[ RhGa .Tematerials The ]. 2 c regular tn[ PtIn 4 , b n 4 and , 18 21 2 Y and (cubic, ] ,i the in ], (3 d a ) Ga, in , in- 2

a) b) c)a d) } 4a 4b 4d 4c a/√2}

FIG. 1. (Color online) The conversion from a (a) cubic Heusler structure to a (b) tetragonal derivative phase is displayed in ′ ′ terms of a non-displasive transformation for a system with the general composition of XX Y Z (X, X , Y - transition metals and Z - main-group element, marked as red, orange, blue, and green, respectively) within the fcc lattice. (c) Relationship of the lattice parameters as atet = acub/√2. the spacegroup (SG) 225 through the introduction of va- respectively, depending on the formation of the regular- cancies or slight structural changes. These modiﬁcations or inverse-type Heusler material, as shown in Table I. alter the structure by breaking the inversion symmetry Other Heusler-related structures [23] are the tetragonal when going from a regular Heusler phase to an inverse or derivatives of the cubic parent phases, which have been half-Heusler structure. widely treated in the context of magnetic shape-memory alloys. The relationship and the unit cell transformation Heusler SG 4d 4c 4b 4a between cubic and tetragonal phases is depicted in Fig. 1.

regular L21 225 Mn Mn Y Ga It is seen that a conventional cubic unit cell can be de-

inverse Xa 216 Mn Y Mn Ga scribed in terms of a tetragonal lattice exhibiting a c/a ratio of √2. The cell parameters are interrelated accord- half C1b 216 Mn Mn Ga ing to ctet = ccub , atet = acub/√2. In this study, a set of TABLE I. The structural relationship of the regular and Mn2-based materials, including transition metals of peri- (3d) inverse structure types for Mn2Y Ga is shown, where ods IV, V, and VI, are considered, namely, Mn2Y Ga, ′ (4d) (5d) X=X =Mn and Z=Ga, as compared to Fig. 1. Only one Mn2Y Ga, and Mn2Y Ga. possible conﬁguration of a half-Heusler type is listed exem- plary.

Heusler materials are generally understood to be in- III. COMPUTATIONAL DETAILS termetallic compounds, distinguishing them from general intermetallics which are forming a broad range of The numerical work was done within density func- solid solutions, with no preferred but statistical occupa- tional theory as implemented in the all-electron FP- tion of crystallographic sites. They are also set apart LAPW code, WIEN2k [24], employing the generalized from other ionic or covalent compounds because Heusler gradient approximation (GGA) in the parametrization systems allow the formation of substitutional series of of Perdew, Burke, and Enzerhof as exchange-correlation single sites. However, they maintain the character of an functional [25]. The angular momentum truncation was ordered compound, and thus they are on the borderline set to lmax = 9 and the number of plane-waves deter- between alloys and compounds. Some distinct systems mined by RKmax = 9 to ensure well-converged calcu- exhibit a tendency to form alloys, nevertheless. The gen- lations. The energy convergence criterion for the self- 5 eral composition is given by XX′YZ, where the classical consistent field calculations was set to 10− Ry, whereas 3 definition has been widened to incorporate quaternary the charge convergence was set to 10− . All calculations materials within the family of Heusler compounds. In the where done on a 20 20 20 k-mesh. For a set of given representation of SG 216 (F 43¯ m), the structure contains c/a ratios, the volumes× were× optimized and fitted to the 3 3 3 four highly symmetric Wyckoff positions: 4d ( 4 , 4 , 4 ), Birch–Murnaghan equation of state [26, 27]. From these, 1 1 1 1 1 1 4c ( 4 , 4 , 4 ), 4b ( 2 , 2 , 2 ), and 4a (0, 0, 0). Depending on the optimal ratio and volume were obtained and the lat- the occupation of the crystallographic positions, two or- tice parameters evaluated. On the basis of the crystallo- dering possibilities for ternary alloys (X = X′) are ob- graphic details, complex magnetic properties such as the tained (Table I). In this study, gallium was chosen as exchange parameters and the corresponding Curie tem- the Z element, whereas one manganese atom, X, occu- peratures (TC) were computed by means of the Korringa– pies position 4d. The second manganese atom, X′, and Kohn–Rostoker (KKR) Green’s function method as im- the other transition metal, Y , are located at 4c or 4b, plemented by the Munich SPR-KKR package [28]. The 3

(a) c/a (cubic reference cell) (b) 0.8 1.0 1.2 1.4 1.6 1.8 0.8 1.0 1.2 1.4 1.6 1.8 0.6 0.6

0.4 0.4

0.2 0.2 t-c

[ eV/ f.u. ] 0 0 E ∆

-0.2 -0.2 Mn VGa Mn FeGa Mn 3 Ga Mn 2 ScGa 2 Mn 3 Ga 2 Mn 2 NiGa Mn 2 TiGa Mn 2 CrGa Mn 2 CoGa Mn 2 CuGa -0.4 -0.4 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 c/a (tetragonal reference cell)

FIG. 2. (Color online) Energetic response to volume-conserving elongations and compressions of the crystal structure along (3d) the c-axis for the Mn2Y Ga-series. Panel (a) includes systems with a valence electron count NV (Y ) 7 (early transition (3d) ≤ metals, ETM) and panel (b) shows the remaining combinations of Mn2Y Ga with NV (Y ) 7, with so-called late transition metals (LTM). ≥

angular momentum expansion was truncated for lmax =3 Tetragonally compressed structures are described by which corresponds to the f-wave symmetry. The energy c/a < √2, whereas tetragonally elongated lattices are integration was done on a complex energy mesh with 48 characterized by c/a > √2, as compared to the cubic points along the integration path, using the Lloyd’s for- parent or austenite phase. Elongation occurs with an mula [29] for an improved estimate of the Fermi energy. increase in the length of the c-axis, whereas the ab- The computation of the exchange parameters is based plane is compressed, leaving the volume approximately on the classical Heisenberg model, which was evaluated unchanged. In this study, the volume of the unit cell, by means of the real-space approach [30]. This provides Vcell, was optimized in addition to the c/a-ratio, and it site- and distance-dependent exchange between sites via was found that no significant change occurred for most infinitesimal rotation of the magnetic moments at a par- cases (Table II). ticular site in real space. To account for distance depen- Figs. 2a and 2b show the E (c/a) curves of the dence, an appropriate truncation of the cluster radius, r, around each atomic site had been chosen. This radius was set to 3.5 a lattice spacings to capture even small Material Ct/c Material Ct/c interactions, as the largest contributions to the effective Mn2ScGa 5.08 Mn2OsGa 2.32 exchange constants are found for radii smaller than 1.5 a Mn2CrGa 6.11 Mn2IrGa 2.44 lattice spacings [15, 31]. Mn3Ga 9.08 Mn2PtGa 3.13 − Mn2FeGa 1.93 Mn2AuGa 1.49 − Mn2CoGa 1.82 Mn2OsSn

IV. RESULTS Mn2NiGa 1.09 Mn2IrSn 1.94 − − Mn2CuGa 2.37 Mn2PtSn 0.35 − A. Lattice Relaxation Mn2MoGa 3.46 Mn2OsIn 0.36 − Mn2RuGa 1.63 Mn2IrIn 1.01 − The total energy E (c/a) as a function of the c/a ratio Mn2RhGa 1.48 Mn2PtIn 0.86 − was calculated and the results are shown in Figs. 2–5. Mn2PdGa 2.78 Mn3Ge 3.98 − The energy zero are deﬁned with respect to the cubic Mn2FeGe 3.80 parent compound, and consequently the energy diﬀer- ence for all phases can be compared easily. The case of TABLE II. Relative volume change between the austenite and Mn3Ga is used as a benchmark and is repeatedly plotted martensite phase in percent, C =(Vtet Vcub) 100/Vcub in Figs. 2 and 3. Table III contains the numerically opti- t/c − · mized lattice parameters. The study reveals that a large (3d) number of the herein-treated materials are most stable Mn2Y Ga series for early transition metals (ETM, in their respective tetragonal structures, with c/a > √2. Y = Sc, Ti, V, Cr) and late transition metals (LTM, 4

(a) c/a (cubic reference cell) (b) 0.8 1.0 1.2 1.4 1.6 1.8 0.8 1.0 1.2 1.4 1.6 1.8 0.6 0.6

0.4 0.4

0.2 0.2

t-c 0 0 [ eV/ f.u. ] E ∆

-0.2 -0.2 Mn 3 Ga Mn 2 NbGa Mn 2 RuGa Mn 2 PtGa Mn 3 Ga Mn 2 WGa Mn 2 PtGa Mn 2 ZrGa Mn 2 MoGa Mn 2 RhGa Mn 2 AgGa Mn TaGa Mn IrGa Mn OsGa Mn AuGa -0.4 2 2 2 2 -0.4 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 c/a (tetragonal reference cell)

FIG. 3. (Color online) Energetic response to volume-conserving elongations and compressions of the crystal structure along (4d) (5d) the c-axis for the (a) Mn2Y Ga and (b) Mn2Y Ga series. Systems involving late transition metals and systems exhibiting at least a local tetragonal energy minimum are shown.

Y = Mn, Fe, Co, Ni, Cu), respectively. A preference only be formed in the series of Mn2Y Ga that includes for tetragonal structures is seen for materials including LTMs, making them derivatives of inverse Heusler sys- (3d) LTM for Mn2Y Ga, with Y being Mn, Fe, or Ni. tems. The information obtained from the lattice opti- In the group of ETM systems, Mn2ScGa was found to mization in terms of the relative positions of the energy be tetragonal, which is an exception. The situation for minima is condensed in Fig. 4, which gives the preferred (4d) (5d) Mn2Y Ga and Mn2Y Ga is seen to be similar. crystal structure visualized in the manner of the periodic In contrast to Mn2CoGa, for which a cubic structure table. (c/a = √2) is preferred, the systems with the same va- Each compound in the family of Mn2Y Ga materials lence electron count (NV = 26) involving heavier species shown in this figure is depicted by one cell that symbol- (ruthenium and osmium) exhibit a global energy mini- izes a transition metal Y of the 3d, 4d, or5d series. The mum for the tetragonal structure. corresponding cell is built up by three subcells, which Discontinuities in the E (c/a) curves are observed for the represent the two variants of chemical coordination (the cases of Mn3Ga and Mn2FeGa, resembling a first-order first or second subcell) and the existence of a global en- transition, whereas E (c/a) for Mn2NiGa is continuous ergy minimum for c/a = √2 (third subcell). The color (Fig. 2b). Strongly composition-dependent modulated code symbolizes the energy6 level of a configuration on martensitic phases [32, 33] and premartensitic phases the energy landscape relative to one another: dark green [34, 35] have been experimentally produced; in combi- - global minimum, light green - local minimum, red - no nation with first-principles calculations [36], the onset of minimum. In cases in which the investigated materials the martensitic transition has been thought to be initi- did not exhibit a cubic minimum, the first two subcells ated by a displacement of atomic planes orthogonal to the are understood to be the type of coordination around crystallographic c-axis of the tetragonal cell. Experimen- the 4d crystallographic site: the symmetry of the coordi- tally [35] and theoretically [36], it has been shown that nating shell is either centrosymmetric (in relation to SG phonon softening along [ζ, ζ, 0] exists in shape-memory 225) or non-centrosymmetric (in relation to SG 216). As materials, and thus the transition has been related to shown in a preceding publication on the cubic variants the occurrence of the tetragonal distortion. [15], systems involving ETMs adopt the L21-type struc- The approach undertaken in the current study resulted ture, whereas compounds containing LTMs are found to in elongated structures in those cases, which were unsta- have the inverse Heusler structure (Xa-type). ble towards a tetragonal distortion. Additionally, global Inspection of Fig. 4 in combination with Figs. 2 and 3 energy minima for distorted phases were found mostly for reveals interesting details, such as the fact that tetrag- the inverse structures, as depicted in Fig. 4. In contrast, onal derivative phases of cubic Heusler alloys, implying compressed variants possessing a global energy minimum that a global energy minimum is present, are observed were not observed in the Mn–Ga Heusler family. The cal- only for a valence electron count of NV 24. It is also culated lattice data and the figures shown (Figs. 2–3) lead clearly seen that the onset of the formation≥ of tetrago- to the conclusion that stable tetragonal structures may nally elongated structures evolves over the periods from 5

225 216 139/119 L2 1 Xa Tet Tet TetCub Tet Sc Ti V Cr Mn Fe Co Ni Cu 225 216 139 225 216 139 225 216 139 225 216 139 225 216 119 225 216 119 225 216 119 225 216 119 225 216 119 Y Zr Nb Mo Tc Ru Rh Pd Ag 225 216 139 225 216 139 225 216 139 225 216 139 225 216 119 225 216 119 225 216 119 225 216 119 225 216 119 Hf Ta W Re Os Ir Pt Au 225 216 139 225 216 139 225 216 139 225 216 119 225 216 119 225 216 119 225 216 119 225 216 119

L2 1 X a

FIG. 4. (Color online) Schematic overview of the preferred site occupancy and crystal structure of Mn2Y Ga Heusler compounds. Stable, metastable, and instable lattices are marked by dark-green, light-green, and red subcells, respectively.

left to right and from lower to higher NV . (NV = 25) behaving analogously to Mn2FeGa (NV = 25) This leads to the question of the mechanism behind this and of Mn2FeGe (NV = 26) being similar to Mn2CoGa distortion, which we are going to approach in Sec. IVA1. (NV = 26). A comparison of the corresponding c/a The lattice parameters for all cubic compounds are found curves demonstrates that this is approximately true and ˚ within a range of ∆ac/√2 =0.35 A for the Ga series. The thus that isoelectronicity is an appropriate concept in the values increase from the borders of the series towards chemistry and physics of Heusler materials. the middle of the range. The same behavior is found Although our calculations agree well with those of previ- for tetragonal compounds, where the range spans from ous works [37], comparison with the experimental data of atet = 3.68 to 4.11 (∆atet = 0.43 A)˚ and ctet =6.91 to Mn–Ni–Ga systems exhibits an interesting discrepancy. 7.48 (∆ctet =0.54 A),˚ whereas Mn2ScGa is an exception. The deviation between experiment and theory was traced The c/a coordinate exhibits inverted behavior, decreas- back to the deviating structural model in terms of the oc- ing from the middle to the left and right borders of the cupation of the involved sites. Neutron diffraction studies series. The overall similarity of the lattice data opens the on Mn2NiGa highlighted the fact that the order is differ- possibility of intermixability with each other and thus ent from the expected Xa-type in the austenite phase. tunability of the whole class of materials. Therefore, the Thus, the chemical formula reads as (Mn,Ni)2MnGa and magnetization may be adjusted over a large range, al- is called L21b-type because the point symmetry includes lowing for the formation of tetragonal compensated fer- inversion symmetry through random occupation of 4d rimagnets via adequate substitution. and 4c sites with Ni and Mn [38]. Similarly in Mn2FeGa, Having discussed the structural trends as a function a deviation from theory was found in an experimental of the valence electron change in the d-electron sys- study because the site occupation was expected to differ tem, the effects of the variation of the main group el- from the perfect MnFeMnGa ordering [16]. Similar issues (4d) ement, Z, from Ga to Ge, In, and Sn in a small sam- are present for the Mn2Y Ga series. From available ple of compounds [Mn2(Mn,Fe)Ge, Mn2(Os,Ir,Pt)In, and data, including neutron diffraction studies, the site occu- Mn2(Os,Ir,Pt)Sn] are briefly discussed. The correspond- pancies have been clarified. Several authors found that ing data, including those for Mn3Ga, are graphed in members of the Mn2RuZ series (Z = Ga, Ge, Si, Sn) Fig. 5. and Mn2RhZ series (Z = Ga) [17, 18, 39] not to be Fig. 5 reveals that the placement of In and Sn at the Z po- tetragonal under the respective reaction conditions. In sition leads to the emergence of a deep energy minimum contrast, they have been realized as cubic Heusler alloys for Mn2IrSn with c/a < √2 and even deeper minima exhibiting a strong degree of anti-site disorder, which has for elongated phases with large c/a-ratios of Mn2OsIn, been characterized as an alloying tendency [18]. Mn2IrIn and Mn2OsSn whose lattice parameters resem- Orthorhombic deformations of the unit cell have been ble those of layered structures. Heusler alloys are of- observed in some systems such as Mn2NiGa [37]. Thus, ten interpreted in terms of a rigid, band model-like ap- the restriction to tetragonal distortions and ordered com- proach, and an electron-filling scheme is employed for pounds in conducted studies leads to a simplified descrip- the prediction of magnetic moments in Co2-based al- tion of these materials. Nevertheless, a general under- loys. Interpretation of the tetragonal instabilities us- standing can be obtained in approaching the Mn2-based ing such an approach leads to the assumption of Mn3Ge Heusler systems through this ansatz. In future studies 6 that aim to predict ground-state structures and magnetic of Y (4c) and Mn(4d), which are found to be more widely configurations, the parameter space for the atomic sites dispersed even though the majority and minority states and relaxations paths for the electronic and spin degrees are separated owing to exchange splitting. The PDOSs of freedom has to be enlarged and restrictions that are of Mn(4d) and Y (4c) are strongly dispersed, with con- widely used have to be dropped. As a consequence of siderable overlap of the spectral weight between the t2g the applied restrictions, disorder effects were not incorpo- and eg states in the majority channel, whereas the minor- rated into this study, thus leaving open any explanations ity channel is gapped, with t2g characterizing the lower of the deviations from experimental results. boundary and eg comprising the upper boundary of the gap. As the DOS is gapped in the minority spin channel c/a (cubic reference cell) the study of the tetragonal system is strongly facilitated 0.8 1.0 1.2 1.4 1.6 1.8 as the states at the Fermi edge in the majority channel 0.6 mostly constitute the origin of the tetragonal distortion. Mn 3 Ga Mn 2 OsIn Mn 2 OsSn The majority PDOSs in the range of 5 eV up to εF Mn FeGe Mn IrIn Mn IrSn − 0.4 2 2 2 exhibit a characteristically shaped peak structure. From Mn Ge Mn PtIn 3 2 Mn3Ga to Mn2CoGa, the majority spin channel (lower 0.2 panels of Figs. 6a–8c) is continuously filled. The Fermi energy is consequently shifted to higher band energies and thus εF sweeps over a range of the majority DOSs, t-c [ eV/ f.u. ] 0 E whereas the minority spin channel remains unchanged. It ∆ is clearly seen from Figs. 6 and 7 that the tetragonal tran- -0.2 sition correlates with the peak structure of the majority DOSs. For εF being centered on a peak of the majority -0.4 DOS, the tetragonal distortion can be triggered. These 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 local maxima are mainly composed of states of the Y (4c) c/a (tetragonal reference cell) (Mn, Fe, Co, Ni) and Mn (4d) atoms. These energy lev- els are of t2g symmetry mainly. In simple interpretation, FIG. 5. (Color online) The energy landscapes, E(c/a), the DOS can be understood in a rigid-band-like fashion. of Mn3Ge, Mn2FeGe, Mn2(Ir,Pt)Sn, and Mn2(Ir,Pt)In are The limit of this interpretation is reached with Mn2NiGa, shown. where the d-PDOS is rearranged and the Slater–Pauling rule is no longer valid for the cubic phase [15]. Compar- ing this to the PDOS of the tetragonally distorted systems, it is observed that the resulting PDOSs are widely 1. Analysis of the Densities of States dispersed and significantly less structured. Mn2NiGa, however, behaves differently. Further filling of the majority states, as intuitively expected, does not occur. In Various attempts to explain the instability of the cu- contrast, the gap in the minority channel closes as the bic phase have been given in the literature using mod- states of e symmetry are pulled towards ε . Thus, the els such as the band Jahn–Teller effect [5, 16], anoma- g F tetragonal distortion in Mn2NiGa is formed by another lous phonon vibrations [49, 50], and Fermi-surface nest- mechanism, which may explain why Mn NiGa is found ing [51]. These different approaches describe the same 2 to be a magnetic shape memory alloy, whereas Mn Ga behavior, i.e., the instability of the cubic phase, from 3 and Mn2FeGa are found in their respective tetragonal different perspectives and extracting different types of crystal structures, although the total energy differences information. Commonly, the densities of states (DOSs) are comparable to that of Mn2NiGa. of related austenite and martensite are compared and We emphasize that these findings differ from other mod- contrasted marking the starting and the endpoint of the (3d) els in which the instability is thought to depend solely transition. As example for the Mn2Y Z series, the on states of Mn(4b) that is found in a tetrahedral envi- DOSs of Mn3Ga, Mn2FeGa, and Mn2NiGa are shown ronment. The DOS at εF in mainly composed of states in Figs. 6–8. For the cubic variants, the partial DOSs of the 4c and 4d positions with minor contributions from (PDOSs) are shown in their corresponding projections on Mn(4b). Here the instability removes the strong peaks of the sites and in terms of projections on the irreducible Mn, Fe, and Co at the 4c position, whereas the states of representations. The peaks in the PDOS of Mn(4b) are Mn(4b) are not rearranged. well separated on the energy scale. On the one hand, this separation is due to the strong crystal field splitting of the Mn(4b) d-states, where the occupied eg states are located in a range between 4 and 3 eV, whereas 2. Spin-Polarization the t PDOS is found between− 1 eV− and the Fermi 2g − edge, εF . On the other hand, the separation of occu- Half-metallicity [52, 53] (complete or nearly complete pied and empty states follows from the strong exchange spin-polarization P (ε ) 100%) is generally observed F ≈ split of the Mn(4b) d-states, in contrast to the d-states in cubic Co2- and Mn2-based Heusler compounds. The 7

TABLE III. Calculated lattice parameters of the cubic parent phases and corresponding tetragonal structures, in comparison with existing literature data. The lattice parameters are given in Angstroms.˚ The saturation magnetization data, MS, are given in µB/ f.u. The listed literature data refer to theoretical (xt) or experimental (xe) investigation of the mentioned cubic (cy) or tetragonal (ty) phases, respectively.

NV SG c/a ac at MS SG c/a at ct MS SG c/a at ct ordering Lit.

Mn2ScGa 20 225 √2 6.15 4.35 4.00 139 1.94 3.98 7.70 5.36 − − Mn2TiGa 21 225 √2 5.95 4.21 2.97 139 225 √2 5.95 ct MnMnTiGa [40] − Mn2VGa 22 225 √2 5.82 4.12 1.98 139 225 √2 5.91 ce MnMnVGa [41],[42] − Mn2CrGa 23 225 √2 5.76 4.07 1.00 139 1.82 3.82 6.95 2.75 225 √2 5.77 ct MnMnCrGa [43] − − Mn3Ga 24 225 √2 5.82 4.12 0.01 139 1.82 3.90 7.08 1.89 139 1.77 3.77 7.16 tt MnMnMnGa [4] − 139 1.82 3.90 7.09 te MnMnMnGa [4]

Mn2FeGa 25 216 √2 5.79 4.09 1.03 119 1.98 3.68 7.29 0.78 119 1.90 3.79 7.19 te Mn(Fe,Mn)2Ga [16] − Mn2CoGa 26 216 √2 5.78 4.09 2.00 119 1.93 3.71 7.14 0.17 216 √2 5.86 ce MnCoMnGa [44]

Mn2NiGa 27 216 √2 5.85 4.14 1.18 119 1.82 3.79 6.91 1.00 216 √2 5.91 ce [45] 119 1.72 3.91 6.70 te [45]

225 √2 5.94 ce (Mn,Ni)2MnGa [38]

139 1.74 3.89 6.77 te (Mn,Ni)2MnGa [38] 216 √2 5.84 ct MnNiMnGa [37]

Mn2CuGa 28 216 √2 5.94 4.20 0.33 216 √2 5.94 ct MnCuMnGa [37]

Mn2ZrGa 21 225 √2 6.14 4.34 3.00 119 − Mn2NbGa 22 225 √2 6.00 4.24 2.00 119 − Mn2MoGa 23 225 √2 5.91 4.18 1.01 119 1.81 3.89 7.04 2.99 − − Mn2RuGa 25 216 √2 5.96 4.22 1.03 119 1.96 3.80 7.45 0.24 216 √2 6.00 ce (Mn 2 ,Ru 1 )3Ga [17] − 3 3 Mn2RhGa 26 216 √2 5.98 4.23 1.64 119 1.94 3.82 7.43 0.10 225 √2 6.03 ce (Mn,Rh)2MnGa [18]

Mn2PdGa 27 216 √2 6.12 4.33 0.55 119 1.84 3.93 7.23 0.93

Mn2AgGa 28 216 √2 6.22 4.40 0.34 119

Mn2HfGa 21 225 √2 6.12 4.33 2.99 119 − Mn2TaGa 22 225 √2 6.00 4.24 1.99 119 − Mn2WGa 23 225 √2 5.92 4.19 0.94 119 − Mn2OsGa 25 216 √2 5.95 4.21 1.02 119 1.97 3.80 7.48 0.28 − Mn2IrGa 26 216 √2 5.97 4.22 2.00 119 1.95 3.83 7.44 0.11

Mn2PtGa 27 216 √2 6.13 4.33 0.44 119 1.87 3.91 7.31 0.75 119 1.38 4.37 6.05 te MnPtMnGa [19]

Mn2AuGa 28 216 √2 6.26 4.42 0.19 119 1.73 4.11 7.13 0.14

Mn2OsSn 25 216 √2 6.21 4.39 1.50 119 1.95 3.97 7.75 0.02 − Mn2IrSn 26 216 √2 6.31 4.46 0.41 119 1.91 4.01 7.67 0.45 119 1.54 4.29 6.59 te MnIrMnSn [46]

Mn2PtSn 27 216 √2 6.39 4.52 0.19 119 1.81 4.15 7.52 0.02 119 1.35 4.51 6.08 te MnPtMnSn [1] − Mn2OsIn 25 216 √2 6.26 4.43 0.62 119 2.02 3.93 7.93 0.27 − Mn2IrIn 26 216 √2 6.30 4.45 0.68 119 1.98 3.97 7.85 0.07

Mn2PtIn 27 216 √2 6.37 4.51 0.31 119 1.84 4.12 7.57 0.38 119 1.57 4.32 6.77 te MnPtMnIn [21]

Mn3Ge 25 216 √2 5.76 4.07 1.01 119 1.90 3.74 7.10 0.98 225 1.91 3.81 7.26 te MnMnMnGe [47] − 225 1.90 3.75 7.12 tt MnMnMnGe [6]

Mn2FeGe 26 216 √2 5.73 4.05 2.01 119 2.05 3.63 7.42 0.06 216 √2 5.80 ct MnFeMnFe [48] − highly symmetric structure and the peculiar electronic of a pseudo-gap in one spin-channel is observed instead. properties due to covalent bonding lead to the appear- ance of a gap in the minority density of states (DOS). The emergence of the tetragonal distortion reduces the spin-polarization of the half-metallic cubic parent phases. The degeneracy of the t2g and eg states is lifted due to the change in local coordination caused by the distortion, which is seen in Figs. 6–8 and Fig. 16, and the emergence 8

4 a) Mn(4b) 225 4 a) Mn(4b) 216

0 0 total total e eg g 4 4 t2g t2g

4 b) Mn(4d) 4 b) Y(4c)

0 0

4 4

139 4 c) Mn(2b) 4 c) Mn(4d)

DOS [1/eV] 0 0 x 2 - y 2 z2 4 xy 4 xz+yz

4 d) Mn(2d) 4 d) Mn(2b)

DOS [1/eV] 0 0 x 2 - y 2 z2 4 xy 4 xz+yz -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 4 e) Fe(2c) 119 E - ε F [eV] 0 FIG. 6. (Color online) DOS of cubic (SG 225) and tetragonal (SG 139) Mn3Ga. 4

4 f) Mn(2d) B. Magnetic Ground State 0 1. The Slater–Pauling Rule 4 Referring to a previous publication [15] on the cubic parent compounds of the investigated materials, the -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 results are presented by means of the Slater–Pauling E - ε [eV] curves. As can be seen in Fig. 10, the Slater–Pauling F rule experiences strong changes, so that calling it the Slater–Pauling rule is done for reasons of convenience. FIG. 7. (Color online) DOS of cubic (SG 216) and tetragonal Fig. 10 visualizes the fact that the magnetization of (SG 119) Mn FeGa. (d) 2 all tetragonal alloys of the Mn2Y Ga family experiences a shift to smaller values. In the case of the (3d) Mn2Y Ga group, this shift is found to be constant throughout the set of compounds, which results in a result in nearly vanishing net moments for Mn2RuGa, linear dependence of the net moment on the valence Mn2RhGa, Mn2PdGa, Mn2OsGa, Mn2IrGa, Mn2PtGa, electron count, thus giving rise to pseudo-Slater–Pauling and Mn2AuGa, as is seen in Fig. 10. Nevertheless, com- behavior, even though half-metallicity and thus integer pensation of spin moments may be achieved for an elec- net moments are not observed (Fig. 9). A decrease of tron count close to NV = 25.7, which can be realized by (4d) the net moment is also found for the Mn2Y Ga and intermixtures of stoichiometric phases of Mn2Y Ga, such Mn Y (5d)Ga compounds. Unlike the lighter compounds, as m = 1.89 µ and a corresponding proportion 2 Mn3Ga − B the changes are not constant over the series and therefore of Mn2NiGa with mMn2NiGa = 1.00 µB or Mn2CoGa. 9

a) Mn(4b) 216 Group Number (Y) 4 34 5 6 7 8 9 10 11 100 0

total / %

)

80 eg F ε ( 4 t2g P 60 4 b) Ni(4c) 40 0 Y (3d ) Mn2 Ga 20 (4d ) Spin- Polarization Y Mn2 Ga 4 5d Mn Y ( )Ga 0 2 4 c) Mn(4d) 2021 22 23 24 25 26 27 28 Number of valence electrons NV 0 FIG. 9. (Color online) The spin-polarization at the Fermi energy for both cubic and tetragonal phases. For a system with 4 NV < 27, the spin-polarization is reduced by the tetragonal distortion. 4 a) Mn(2b) 119

DOS [1/eV]

0 x 2 - y 2 z2 16 20 24 28 16 20 24 28 32 4 xy xz+yz b) c) 4 Rh Ir 4 b) Ni(2c) Pd Pt Cr Ru Cr Os Au Ag 0 0 Mo W Nb Ta Zr Hf

B (4d) (5d) 4 µ -4

[ ] Y Y Mn2 Ga Mn2 Ga

m 12 4 c) Mn(2d) a) Fe

8 Co 0

4 4 Ni Magnetic moment Co Ni Mn Ga Fe -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 3 Cu 0 E - ε Cr F [eV] V Ti Sc Y (3d) -4 Mn2 Ga FIG. 8. (Color online) DOS of cubic (SG 216) and tetragonal (SG 119) Mn2NiGa. 16 20 24 28 32 36 40 44 48 Number of valence electrons NV

FIG. 10. (Color online) Slater–Pauling curves of (3d) (4d) Thus, fractions of 0.435 mMn2NiGa and 0.565 mMn2FeGa (a) Mn2Y Ga, (b) Mn2Y Ga compounds, and (c) (5d) could ideally lead to complete compensation of the mag- Mn2Y Ga compounds. netization. A similar approach was undertaken by Nayak et al. that obtained a compensated ferrimagnet by varying the Mn/Pt-ratio in Mn3 xPtxGa leading to complete compensation of magnetization− for x 0.59 theo- ≈ retically [20]. In Mn2RuxGa thin ﬁlms the compensation of the spin-moment has been achieved through variation of the ruthenium concentration [54]. 10

Y (3d) Y (4d) Y (5d) (a) Mn2 Ga (b)Mn2 Ga (c) Mn2 Ga ScTi V Cr Mn Fe Co Ni Cu Mo Tc Ru Rh Pd Ag Os Ir Pt Au

4 L2 1 Xa L2 1 Xa Xa 4

[µ / f.u.] m Mn(4b) m Mn(4b) m Mn(4b)

m 2 2 m m Y(4b) m Y(4b) Y(4c) m Y(4c) m Y(4c) 0 0

mMn(4d) m Mn(4d) m Mn(4d) -2 -2 mMn(8c) mMn(8c)

Total magnetic moment Total 2021 22 232425 26 27 28 232425 26 27 28 2526 27 28 N Number of valence electrons V

(3d) (4d) (5d) FIG. 11. (Color online) Atomic magnetic moments in (a) Mn2Y Ga, (b) Mn2Y Ga, and (c) Mn2Y Ga compounds. Open symbols denote the L21-type coordination. Filled symbols denote the Xa-type coordination. The squares stand for the cubic materials, whereas the tetragonal systems are represented by rectangles.

2. Local Magnetic Moments systems, the nearest-neighbor distance basically does not change, whereas the direction of the nearest-neighbor interaction does (Fig. 16). In the case of Mn3Ga, Ct/c The change of the total magnetic moment is to be un- is approximately 9%, and thus may suppress the shape derstood in terms of the site moments. Inspection of memory effect in this material. Fig. 11 clearly reveals the dependencies. The change in the local moments indicates the specific impact of the distortion on the single sites. The local moments responding strongest to elongation/compression of the crystal axes 3. Magnetocrystalline Anisotropy are members of the magnetic sublattices formed by the Mn(8c), Mn(4d), and Y (4c) sites (the former so-called Inherent in noncubic crystals is a directional preference tetrahedral sites as compared to SG 225). Depending of magnetization that is absent in cubic materials, which on the characters of the local moments, which are either is related to the tetragonal modification of the crystal itinerant or localized in nature [12, 52], we find major dif- axes. The magnetocrystalline anisotropy (MCA) energy ferences in the influence of the tetragonal distortion on is here defined as the energy difference between states these sites. In Fig. 11, the localized moment of Mn(4b) is with magnetization pointing along the z-axis and the found in the upper part of the plot for a positive value of x- or y-axis, that is, E = E E , whereas MCA (100) − (001) the magnetic moment of approximately 3 µB. The effect other crystallographic directions are not considered. The of elongation along the c-axis and compression of the ab- anisotropy energy is phenomenologically thought to de- plane has a stronger influence on moments of itinerant pend on the value of the c/a ratio, which is more or character found in the lower part of the plot, referring less equal for most of the compounds investigated in the to the 4d and 4c sites. These are located in the same present study. Therefore, the underlying mechanism is lattice plane (compare the atoms depicted as red and understood as a band-filling effect, which affects the spin blue spheres in Fig. 1). Manganese on site 4b [Mn(4b)] orbit coupling (SOC) symmetry. This interpretation can exhibits a large localized moment of 3 µB, is thus is gener- be directly taken from Figs. 12 and 13. Increasing the ally much less affected. Apart from the magnitude of the SOC strength by varying the Y element through the 3d, local moments, Fe(4c) in Mn2FeGa exhibits a spin-flip 4d, and 5d series increases the MCA energy by a factor of from parallel to anti-parallel alignment of the Mn(4b)- approximately 3 for Y = Fe, Ru, Os. The effect of band Fe(4c) interaction upon the tetragonal distortion. Apart filling is deciphered by sweeping the Y elements along a from the changes in magnitude of the local atomic mo- series. Going from left to right in any set of compounds (xd) ments, the effective anti-parallel coupling of the nearest Mn2Y Ga, the MCA is altered from preferred out-of- neighbor manganese atoms does not suffer from the struc- plane to in-plane orientation. The same situation holds (5d) tural transformation, as will be quantified in terms of the for Mn2Y Z compounds, whose preferred orientation exchange interaction constants in Sec. IVC. Since the is graphed in Fig. 13. Two compounds in this figure de- volume change Ct/c is in the order of 1-2% for the most serve special attention, which are Mn2PtIn and Mn2IrSn. 11

Previously they were predicted to possess a non-collinear magnetic order does not change. A decrease or change magnetic order just as Mn2RhSn, which was described of sign of the coupling constant JMn(4d) Mn(4d) is typical in great detail in Ref. [46]. Although we did not con- owing to preferential anti-parallel coupling− as caused for sider non-collinear order in the present case, it might be short Mn–Mn distances. Instead, the cubic to tetragonal worth while to point out that it is the Mn-atom on site 4d transition even results in increased values, thus leading that reacts by canting to the incipient spin-reorientation to an increase of TC. Prima facie, the overall trend of that is seen clearly in Figs. 12. Similar physics is ob- increased TCs cannot be traced back to a common mech- served in the famous Rare-Earth magnets Nd2Fe14B and anism. In contrast, there is no exchange interaction that Er2Fe14B. Manganese thus shares properties with Rare exhibits a similar behavior over the series. For instance, Earths [55–57]. the findings in the case of Mn2FeGa are related to the re- duction of magnetic frustration, which is due to the com- Group Number (Y ) peting anti-parallel interaction Mn(4b)–Mn(4d), and par- 3 4 5 6 7 8 9 10 11 allel interactions of Fe(4c) with both Mn(4d) and Mn(4b) neighbors. Upon the tetragonal distortion a spin-flip of L2 1 Xa 3 34.81 Mn 2 OsGa Mn 2 IrGa the Fe(4c) local moment is observed, due to the change of the sign of Mn(4b)–Fe(4c) interaction. Upon the dis- 2 23.21 tortion, the strength of Mn(4b)–Mn(4d) interaction is altered by approximately ∆JMn(4b) Mn(4b) = 56 meV in Mn MoGa − 2 Mn RuGa Mn RhGa Mn PdGa meV / f.u. 1 2 2 2 11.60 Mn2FeGa. [ ] Mn ScGa 2 Mn 2 NiGa A similar, but smaller effect is found in Mn2CoGa, indi- Mn CrGa Mn Ga MCA out of plane (z) 2 3 E 0 [ ] cating the magnetic frustration that had been present in 0 in plane (x,y) Mn FeGa Mn CoGa Mn AuGa 2 2 2 TK the cubic phase and the weakening of the exchange inter- Mn 2 PtGa action Y (4c)–Mn(4b) (Y = Fe, Co) due to the tetragonal -1 Mn Y (3d) Ga Mn Y (4d) Ga Mn Y (5d) Ga -11.60 2 2 2 distortion. This might be one of the contributions pro- 20 21 22 23 24 25 26 27 28 hibiting the shape memory effect in Mn2FeGa. The main N Number of valence electrons V contribution to TC, the Mn(4d)–Mn(4b) exchange, does not suffer from the structural transition. Similarly, the FIG. 12. (Color online) Calculated MCA energy of Mn2Y Ga. influence of the distortion on the exchange causes the smaller contributions to be altered and some frustration to be diminished. For example, the Mn(4d)–Mn(4d) interaction vanishes, with preferred anti-parallel alignment L2 1 Xa in the cubic case, whereas the major Mn(4d)–Mn(4b) in- 3 34.81 Y (3d) teraction remains unchanged. Mn2 Ge Mn Y (5d) In Exceptions to the general observation of increased T in 2 2 23.21 C Mn 2 IrIn Y (5d) Mn2 Sn tetragonally distorted phases are Mn2NiGa, Mn2PdGa,

Mn 3 Ge and Mn2PtGa systems, in which the TCs are reduced

meV / f.u. 1 11.60 [ ] upon the tetragonal distortion. The signiﬁcant reduc-

Mn 2 RhSn MCA out of plane (z) tion (Fig. 14) is caused by a weakened Mn(4d)–Mn(4b) E 0 0 [ ] in plane (x,y) interaction (Fig. 15) that may indicate an unstable mag- Mn 2 IrSn TK Mn 2 PtSn netic groundstate. A relation to the Heusler compound -1 -11.60 Mn RhSn [46] can theoretically be established as these Mn 2 PtIn 2 materials possess the same number of valence electrons, 22 23 24 25 26 27 28 29 N . Mn RhSn has been shown to exhibit a non-collinear Number of valence electrons NV V 2 magnetic groundstate due to canting of the diﬀerent man- FIG. 13. (Color online) Calculated MCA energy of some cho- ganese moments [46]. Non-collinear spin conﬁgurations have not been considered in this work, so that canting of sen Mn2Y Z alloys: Mn3Ge, Mn2IrIn, Mn2IrSn, Mn2PtIn, and Mn2PtSn. the spins or spin spiral groundstates in these materials may be present.

C. Exchange Coupling and Curie Temperatures V. SUMMARY

The details of the calculations pertaining to Curie tem- Using total energy calculations within density func- peratures (TC) have been described by us previously [15]. tional theory, we investigated in detail the response to Even though the in-plane next-nearest-neighbor distance tetragonal distortions for a large set of cubic Heusler com- (3d,4d,5d) decreases due to the tetragonal distortion, the next- pounds, Mn2Y Ga, and some other chosen mate- nearest neighbor Mn(8c)-Mn(8c) coupling is still positive rials. We were able to single out the systems that remain for the L21-derived tetragonal phases; thus, the overall cubic from those that favor a tetragonal structure. The 12

ScTi V Cr Mn Fe Co Ni Cu Pc Pt TC,c TC,t ∆TC,t−c EMCA Ku L2 Mn2ScGa 87 35 464 0.616 1.62 50 1 Xa Mn2TiGa 83 557 0 Mn2VGa 94 587 Mn2CrGa 97 50 578 970 392 0.779 2.46 -50

Mn3Ga 96 56 221 610 389 0.906 2.7 [meV]

eff -100

Mn2FeGa 95 56 601 848 247 0.359 1.16 J Mn2CoGa 93 60 928 1124 196 0.236 0.77 -150 Mn2NiGa 35 42 1005 750 255 0.193 0.62 Mn(4d) - Mn(4d) Mn(4d) - Mn(4d) − -200 Mn(4d) - Mn(4c) Mn(4d) - Y (4c) Mn2CuGa 53 1491 Mn(4d) - Y (4b) Mn(4d) - Mn(4b) Mn2ZrGa 82 207 -250 Y (4b) -Y (4b) Mn(4b) - Mn(4b) Mn NbGa 98 289 2 20 21 22 23 24 25 26 27 28 Mn2MoGa 85 65 140 335 196 0.636 1.91 Number of valence electrons NV Mn2RuGa 95 1 619 1315 696 0.564 1.68 Mn2RhGa 15 59 576 1351 776 0.322 0.95 FIG. 15. (Color online) The evaluated eﬀective exchange in- Mn2PdGa 7 46 809 335 473 0.040 0.11 teraction parameters are shown and compared for all cubic − Mn2AgGa 24 1240 and tetragonal cases. Thus, the underlying mechanism of the increase of TC due to the distortion is visualized. Circles Mn2HfGa 89 correspond to cubic compounds and rectangles to tetragonal Mn2TaGa 96 compounds. Mn2WGa 83

Mn2OsGa 96 273 1075 802 3.270 9.72

Mn2IrGa 74 5 411 1122 711 2.388 7.02 Mn Mn2PtGa 23 51 799 326 472 0.293 0.84 Z − − − Mn2AuGa 8 26 1027 897 130 0.731 1.94 − − −

TABLE IV. The calculated Curie temperatures in Kelvin of Mn tetragonal TC,t and cubic TC,c parent compounds (taken from Ref. 15). The changes due to the tetragonal transformation are listed as ∆TC,t−c. EMCA represents the magnetocrystalline anisotropy energies in meV per formula unit, whereas MJ the anisotropy constant, K , is given in 3 . u m Mn Z

Group Number (Y ) 4 5 6 7 8 9 10 11 FIG. 16. (Color online) A two-dimensional projection of Cu the nearest-neighbor coordination of Mn(4b) is shown. The Mn Y (3d) Ga symmetry of the coordination changes while undergoing the 1400 2 Ru Rh Y (4d) Mn2 Ga Ag martensitic transition, whereas the nearest-neighbor distance [K] Y (5d) Ir remains unchanged as the volume change is on the order of c 1200 Mn2 Ga Os T Ni only 1–9%. Co 1000 Fe 800 Mn Au Ti V Cr details of the total energy as a function of the distortion 600 were found to be similar for materials exhibiting the same 400 Pt number of valence electrons. The magnetizations of the Curie Temperature tetragonal alloys were found to be shifted to smaller val- 200 Nb Mo Pd Zr ues, which we could attribute to changes of the itinerant local moments. This led to characteristic modiﬁcations 21 22 23 24 25 26 27 28 of the Slater–Pauling curve. By means of partial den- N Number of valence electrons V sities of states, the changes to the electronic structures revealed the microscopic origin of the observed trends. FIG. 14. (Color online) The calculated Curie temperatures As compared to the cubic parent phases, a strengthening of the Heusler compounds containing Ga. The results shown of the exchange interaction between neighboring sites was here are obtained in the Mean-Field Approximation (MFA) observed, which resulted in an increase of the Curie tem- and highlight the consequence of the structural relaxation. perature. Focusing our attention on the magnetocrys- Squares correspond to cubic compounds and rectangles to tetragonal compounds. talline anisotropy we observed an interesting trend that describes a spin reorientation over our series of com- 13 pounds, furthermore very large anisotropies are found ACKNOWLEDGMENTS for tetragonal Heusler compounds containing heavy transition metals accompanied by low magnetic moments, Financial support from the Deutsche Forschungsge- which indicates that these materials may be promising meinschaft (DfG) (research unit FOR No. 1464 ASPI- candidates for spin-transfer torque applications. MATT, project P 1.2-A) and European Research Council Advanced Grant (ERC-AG) No. 291472 IDEA Heusler! is gratefully acknowledged.

[1] J. Winterlik, S. Chadov, A. Gupta, V. Alijani, T. Gasi, C. Felser, Nat. Mat. , 3 (2015). K. Filsinger, B. Balke, G. H. Fecher, C. A. Jenkins, [21] A. K. Nayak, C. Shekhar, J. Winterlik, A. Gupta, and F. Casper, J. Kuebler, G.-D. Liu, L. Gao, S. S. P. Parkin, C. Felser, Appl. Phys. Lett. 100 (2012). and C. Felser, Adv. Materials 24, 6283 (2012). [22] F. Heusler, Verh. d. DPG 5, 219 (1903). [2] L. Wollmann, G. H. Fecher, S. Chadov, and C. Felser, [23] J. C. Suits, Phys. Rev. B 14 (1976). J. Phys. D: Appl. Phys. 48, 164004 (2015). [24] P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, [3] S. Chadov, S. W. D’Souza, L. Wollmann, J. Kiss, G. H. and J. Luitz, WIEN2k, An Augmented Plane Wave + Fecher, and C. Felser, Phys. Rev. B 91, 094203 (2015). Local Orbitals Program for Calculating Crystal Properties [4] B. Balke, G. H. Fecher, J. Kübler, and C. Felser, Appl. (Techn. Universitat Wien, Austria, 2001). Phys. Lett. 90, 152504 (2007). [25] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. [5] J. Winterlik, B. Balke, G. H. Fecher, C. Felser, M. C. M. Lett. 77, 3865 (1996). Alves, F. Bernardi, and J. Morais, Phys. Rev. B 77, [26] F. Birch, Phys. Rev. 71, 809 (1947). 054406 (2008). [27] F. D. Murnaghan, Proc. Natl. Acad. Sci. USA 30, 244 [6] D. Zhang, B. Yan, S.-C. Wu, J. Kübler, G. Kreiner, (1944). S. S. P. Parkin, and C. Felser, J. Phys.: Condens. Matter [28] H. Ebert, D. Ködderitzsch, and J. Minár, Rep. Prog. 25, 206006 (2013). Phys. 74, 096501 (2011). [7] H. Kurt, K. Rode, M. Venkatesan, P. S. Stamenov, and [29] P. Lloyd and P. V. Smith, Adv. Phys. 21, 69 (1972). J. M. D. Coey, Phys. Rev. B 83, 020405 (2011). [30] A. I. Liechtenstein, M. I. Katsnelson, P. V. Antropov, [8] S. S. P. Parkin, M. Hayashi, and L. Thomas, Science and A. V. Gubanov, J. Magn. Magn. Materials 67, 65 320, 190 (2008). (1987). [9] V. Chernenko, E. Cesari, V. Kokorin, and I. Vitenko, [31] M. Meinert, J.-M. Schmalhorst, and G. Reiss, J. Phys.: Scripta Metallurgica et Materialia 33, 1239 (1995), shape Condens. Matter 23, 116005 (2011). memory, SMSA, first Ni2MnGa. [32] G. Fritsch, V. V. Kokorin, V. A. Chernenko, A. Kempf, [10] K. Ullakko, J. K. Huang, C. Kantner, R. C. O’Handley, and I. K. Zasimchuk, Phase Transitions 57, 233 (1996). and V. V. Kokorin, Appl. Phys. Lett. 69, 1966 (1996). [33] V. A. Chernenko, C. Segu´ı, E. Cesari, J. Pons, and V. V. [11] V. Alijani, J. Winterlik, G. H. Fecher, and C. Felser, Kokorin, Phys. Rev. B 57, 2659 (1998). Appl. Phys. Lett. 99, 222510 (2011). [34] A. Zheludev, S. M. Shapiro, P. Wochner, A. Schwartz, [12] M. Meinert, J.-M. Schmalhorst, C. Klewe, G. Reiss, M. Wall, and L. E. Tanner, Phys. Rev. B 51, 11310 E. Arenholz, T. Böhnert, and K. Nielsch, Phys. Rev. (1995). B 84, 132405 (2011). [35] U. Stuhr, P. Vorderwisch, V. V. Kokorin, and P.-A. [13] S. Ouardi, T. Kubota, G. H. Fecher, R. Stinshoff, Lindg˚ard, Phys. Rev. B 56, 14360 (1997). S. Mizukami, T. Miyazaki, E. Ikenaga, and C. Felser, [36] A. T. Zayak, P. Entel, J. Enkovaara, A. Ayuela, and Appl. Phys. Lett. 101, 242406 (2012). R. M. Nieminen, Phys. Rev. B 68, 132402 (2003). [14] T. Kubota, S. Ouardi, S. Mizukami, G. H. Fecher, [37] A. Chakrabarti, M. Siewert, T. Roy, K. Mondal, C. Felser, Y. Ando, and T. Miyazaki, J. Appl. Physics A. Banerjee, M. E. Gruner, and P. Entel, Phys. Rev. 113 (2013). B 88, 174116 (2013), checked. [15] L. Wollmann, S. Chadov, J. Kübler, and C. Felser, Phys. [38] P. J. Brown, T. Kanomata, K. U. Neumann, B. Oulad- Rev. B 90, 214420 (2014). diaf, A. Sheikh, and K. R. A. Ziebeck, J. Phys.: Condens. [16] T. Gasi, A. K. Nayak, J. Winterlik, V. Ksenofontov, Matter 22, 506001 (2010). P. Adler, M. Nicklas, and C. Felser, Appl. Phys. Lett. [39] K. Endo, T. Kanomata, H. Nishihara, and K. R. A. 102 (2013), checked. Ziebeck, J. Alloys. Compd. 510, 1 (2012). [17] T. Hori, M. Akimitsu, H. Miki, K. Ohoyoama, and [40] M. Meinert, J.-M. Schmalhorst, and G. Reiss, J. Phys.: Y. Yamaguchi, Appl. Phys. A 74, 737 (2002). Condens. Matter 23, 036001 (2011). [18] G. Kreiner, A. Kalache, S. Hausdorf, V. Alijani, J.-F. [41] K. Ramesh Kumar, N. Harish Kumar, G. Markandeyulu, Qian, G. Shan, U. Burkhardt, S. Ouardi, and C. Felser, J. Arout Chelvane, V. Neu, and P. D. Babu, J. Magn. Z. Anorg. Allg. Chem. 640, 738 (2014). Magn. Materials 320, 2737 (2008). [19] A. K. Nayak, M. Nicklas, S. Chadov, C. Shekhar, Y. Sk- [42] K. H. J. Buschow and P. G. van Engen, J. Magn. Magn. ourski, J. Winterlik, and C. Felser, Phys. Rev. Lett. 110, Materials 25, 90 (1981). 127204 (2013). [43] H. Luo, Z. Zhu, G. Liu, S. Xu, G. Wu, H. Liu, J. Qu, [20] A. K. Nayak, M. Nicklas, S. Chadov, P. Khuntia, and Y. Li, J. Magn. Magn. Materials 320, 421 (2008). C. Shekhar, A. Kalache, M. Baenitz, Y. Skourski, V. K. [44] G. D. Liu, X. F. Dai, H. Y. Liu, J. L. Chen, Y. X. Li, Guduru, A. Puri, U. Zeitler, J. M. D. Coey, and G. Xiao, and G. H. Wu, Phys. Rev. B 77 (2008). 14

[45] G. D. Liu, J. L. Chen, Z. H. Liu, X. F. Dai, G. H. Wu, [51] S. R. Barman, S. Banik, A. Shukla, C. Kamal, and B. Zhang, and X. X. Zhang, Appl. Phys. Lett. 87, 262504 A. Chakrabarti, Europhys. Lett. 80, 57002 (2007). (2005). [52] J. Kübler, A. R. Williams, and C. B. Sommers, Phys. [46] O. Meshcheriakova, S. Chadov, A. K. Nayak, U. K. Rev. B 28, 1745 (1983). Rössler, J. Kübler, G. André, A. A. Tsirlin, J. Kiss, [53] R. A. de Groot, F. M. Mueller, P. G. van Engen, and S. Hausdorf, A. Kalache, W. Schnelle, M. Nicklas, and K. H. J. Buschow, Phys. Rev. Lett. 50, 2024 (1983). C. Felser, Phys. Rev. Lett. 113, 87203 (2014). [54] H. Kurt, K. Rode, P. Stamenov, M. Venkatesan, Y.-C. [47] H. Kurt, N. Baadji, K. Rode, M. Venkatesan, P. S. Sta- Lau, E. Fonda, and J. M. D. Coey, Phys. Rev. Lett. 112, menov, S. Sanvito, and J. M. D. Coey, Appl. Phys. Lett. 027201 (2014). 101, 132410 (2012). [55] D. Haskel, J. C. Lang, Z. Islam, A. Cady, G. Srajer, [48] H. Z. Luo, H. W. Zhang, Z. Y. Zhu, L. Ma, S. F. Xu, M. van Veenendaal, and P. C. Canfield, Phys. Rev. Lett. G. H. Wu, X. X. Zhu, C. B. Jiang, and H. B. Xu, J. 95, 217207 (2005). Appl. Physics 103, 083908 (2008). [56] W. B. Yelon and J. F. Herbst, J. Appl. Physics 59, 93 [49] A. T. Zayak, P. Entel, K. M. Rabe, W. A. Adeagbo, and (1986). M. Acet, Phys. Rev. B 72, 054113 (2005). [57] J. Chaboy, L. M. Garc´ıa, F. Bartolomé, A. Marcelli, [50] P. Entel, V. D. Buchelnikov, V. Khovailo, A. T. Zayak, G. Cibin, H. Maruyama, S. Pizzini, A. Rogalev, J. B. W. A. Adeagbo, M. E. Gruner, H. C. Herper, E. F. Goedkoop, and J. Goulon, Phys. Rev. B 57, 8424 (1998). Wassermann, and V, J. Phys. D: Appl. Phys. 39, 865 (2006).