First-principles study of the structural stability of Mn 3 Z (Z=Ga, Sn and Ge) Heusler compounds Delin Zhang 1, Binghai Yan 1,2,3, Shu-Chun Wu 1, Jürgen Kübler 4, Guido Kreiner 1, Stuart S. P. Parkin 4, Claudia Felser 1,2 1 Max-Planck-Institute für Chemische Physik fester Stoffe, Nöthnitzer Str. 40, 01187 Dresden, Germany 2 Johannes Gutenberg-Universität Mainz, Staudingerweg 9, 55128 Mainz, Germany 3 Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Straße 38 01187 Dresden, Germany 4 Institut für Festkörperphysik, Technische Universität Darmstadt, 64289 Darmstadt, Germany 5 IBM Almaden Research Center San Jose, CA 95120, USA. Abstract We investigate the structural stability and magnetic properties of cubic, tetragonal and hexagonal phases of Mn 3 Z (Z=Ga, Sn and Ge) Heusler compounds using firstprinciples density-functional theory. We propose that the cubic phase plays an important role as an intermediate state in the phase transition from the hexagonal to the tetragonal phases. Consequently, Mn 3 Ga and Mn 3 Ge behave differently from Mn 3 Sn, because the relative energies of the cubic and hexagonal phases are different. This result agrees with experimental observations from these three compounds. The weak ferromagnetism of the hexagonal phase and the perpendicular magnetocrystalline anisotropy of the tetragonal phase obtained in our calculations are also consistent with experiment. 1
I. Introduction The Mn 3 Z (Z=Ga, Sn and Ge) type of Heusler compounds can have three different structural phases, where each phase exhibits different magnetic properties. The hexagonal phase has been known for decades. Mn atoms form a Kagome lattice in a plane with a Z atom in the center of the hexagon. Here, Mn atoms have a triangular antiferromagnetic (AFM) coupling with a weak net magnetic moment [1~8]. The cubic phase is the standard full Heusler structure. In this phase, Mn atoms are present in two unique lattice sites; these different sites have magnetic moments with opposite directions, leading to ferrimagnetic (FiM) order [9, 10]. The cubic phase has a high density of states at the Fermi energy, and hence a Peierls transition could occur [11], giving rise to the third phase, the tetragonal phase. The tetragonal phase can be treated as a cubic phase with a distortion along the z direction. This distortion makes the magnetic moments favor the z axis, meaning the system possesses perpendicular magnetocrystalline anisotropy (PMA) [11], which promises great potential for future high-density spin-transfer torque applications [12,13,14~21]. In experiments, the hexagonal phase of Mn 3 Z (Z=Ga, Sn and Ge) was synthesized by annealing the samples at high temperatures [1~3, 22], while the tetragonal phase of Mn 3 Ga and Mn 3 Ge was realized by low temperature annealing [1, 3, 23, 24]. The cubic phase has not been observed so far, since it may be unstable as mentioned above. On the other hand, the tetragonal phase of Mn 3 Sn has not yet been reported in the literature, although it is expected to behave similarly to the other two compounds. The possible transition between hexagonal, cubic and tetragonal phases catches the attention of researchers [3, 25], but a comprehensive study is yet to be undertaken. In this paper, we investigate the structural stability and magnetic properties of the hexagonal, cubic and tetragonal phases for the three compounds Mn 3 Z (Z=Ga, Sn and Ge) using first-principles calculations. The tetragonal phase of Mn 3 Sn is found to have the lowest total energy among the three phases, similar to Mn 3 Ga and Mn 3 Ge. However, the cubic phase of Mn 3 Sn has higher energy than its hexagonal phase, in contrast to the Mn 3 Ga and Mn 3 Ge compounds. In the case of Mn 3 Sn, we suggest it is difficult for the 2
hexagonal phase to transform from the cubic phase into the tetragonal one. This might explain the lack of the tetragonal Mn 3 Sn phase. II Calculation details First-principles calculations were carried out using Vienna ab initio Simulation Package (VASP) [26]. The ions were described using projector augmented wave (PAW) potentials [27]. The generalized gradient approximation (GGA) in the Perdew-Bruke- Ernzerhof (PBE) form [28] was adopted to describe the exchange-correlation interactions between electrons. An energy cutoff of 500 ev was used for the plane wave basis. Spin-orbit coupling was employed in all calculations to describe the non-collinear spin polarization and magnetocrystalline anisotropy. The volume and shape (c/a) of the cubic, tetragonal and hexagonal structures were fully relaxed to get the stable structural configurations with the lowest energy. In addition, the full potential linear augmented plane wave (FLAPW) method [29] was also used to check the validity of the pseudopotential calculations. III. Results and discussion As shown in Figure 1, the cubic phase belongs to the X 2 YZ full Heusler structure (Fm-3m). Mn II (X position), Mn I (Y position) and Z (Z=Ga, Sn and Ge) atoms occupy the (1/4, 1/4, 1/4), (1/2, 1/2, 1/2) and (0, 0, 0) sites, respectively. The magnetic moments of Mn I and Mn II orient oppositely. The tetragonal phase (I4/mmm) has an elongated c axis and shortened a axis as compared to the cubic lattice. The magnetic order, as in the cubic phase, is ferrimagnetic. Our optimized lattice parameters and magnetic moments (per Mn 3 Z unit) agree well with previous experiments [1, 3, 21, 24] and calculations [17, 23] (see Tables I and II). The PMA energy is defined as the energy difference between the easy direction (001) and the in-plane direction (010). In the tetragonal phase, the PMA is around 1 mev for all three compounds, consistent with previous calculations on Mn 3 Ga and Mn 3 Ge [11, 21]. The Mn atoms in the hexagonal phase exhibit various possible magnetic configurations. On the one hand, in the plane of the triangle lattice, the magnetic moments of Mn atoms may point in-plane or out-of-plane. On the other hand, one finds 3
that the Mn-Mn bonds between neighboring layers in of Mn 3 Z are a little shorter than the in-plane Mn-Mn bonds. Therefore the interlayer magnetic coupling is also important. As a summary, we display the most important configurations in Fig. 2. The AFM type of in-plane ordering in Fig. 2(d) is found to be the most stable. However, the directions of the magnetic moments are not equally separated by 120, therefore they do not fully cancel each other, resulting in a weak ferromagnetic phase [30~32]. The calculated magnetic moments also agree with previous experiments for all three materials. For example, the magnetic moments are 0.03 µ B (exp. 0.045 µ B [1]) for Mn 3 Ga, 0.01 µ B (exp. 0.009 µ B [2]) for Mn 3 Sn and 0.01 µ B (exp. 0.06 µ B [3]) for Mn 3 Ge. In Mn 3 Ge, the variation between the experimental and calculated results originates mostly from the fact that the experimental compositions are off-stoichiometry [3, 22]. It is important to compare the energetic stability of the three phases for Mn 3 Z (Z=Ga, Sn and Ge). Figure 3 shows the total energy dependence on the volume of hexagonal, cubic and tetragonal structures for these three compounds. We relaxed all the structure parameters (volume and shape (c/a) of the lattice) to reach the most stable structures. In the tetragonal and hexagonal phases, for example, the c/a-ratio was fully optimized for any given value of the volume. One unambiguously finds in all three compounds that the tetragonal phase is energetically more stable than the cubic and hexagonal phases. For Mn 3 Ga and Mn 3 Ge, the hexagonal phase has the highest total energy, while the cubic phase exhibits the highest energy for Mn 3 Sn. Although the cubic phase is claimed to be unstable in experimental work [23], its relative energy between these three phases is very important in understanding the stability of the tetragonal and hexagonal phases. We propose that the hexagonal phase does not change into the tetragonal phase directly; rather it transitions through the cubic lattice, which is structurally intermediate. It is simple to see that distortion along the c direction of the cubic phase will lead to the tetragonal phase, while the transition between the cubic and hexagonal phases is not as straightforward. If the cubic lattice is projected along the diagonal direction, one obtains a trianglar lattice. In order to recover a hexagonal phase, four atomic layers including three Mn layers and one Z layer should be compressed into one layer. However, this cannot be realized by a simple projection of ABC sites, for the Z atom overlaps with 4
one Mn atom in this way (Figs. 1(d) and 1(e)). In order to host the additional Mn atoms inside a layer, the honeycomb lattice of Mn (Fig. 1(d)) should change into a Kagome lattice (Figs. 1(f) and 1(g)) to create more available sites for Mn atoms. On the one hand, the transition from the cubic to the hexagonal phase requires compressive pressure along the cubic diagonal direction in order to push Mn and Z atoms inside a plane. On the other hand, the transition from the hexagonal to the cubic phase also needs compressive strain in the ab plane of the hexagonal lattice. For Mn 3 Ga and Mn 3 Ge, the cubic phase is energetically between the hexagonal and tetragonal phases. In this case, the hexagonal phase can easily pass through the cubic phase into the tetragonal phase. Therefore, both phases can be realized in experiments [1, 3, 23, 24]. In contrast to the above two compounds, Mn 3 Sn has a cubic phase whose energy is higher than the hexagonal phase. Consequently, it is difficult to transform from an existing hexagonal phase into a tetragonal one, possibly explaining why tetragonal Mn 3 Sn has not been synthesized so far. We can explain the transition between the hexagonal and tetragonal phases using the cubic phase. However, there still remains a puzzle: the tetragonal structure of Mn 3 Sn has lower energy than the hexagonal one in our calculations, while experiments only observe the hexagonal phase. This may be due to the structural disorder or offstoichiometry compositions that commonly exist in experiments [2]. In these cases, the hexagonal phase may have lower energy than the tetragonal counterpart. We will investigate these effects in a future work. IV. Conclusions We calculated the structures and compared the stabilities of the cubic, hexagonal and tetragonal phases for Mn 3 Ga, Mn 3 Ge and Mn 3 Sn. The cubic phase plays an important role as an intermediate state in the phase transition from the hexagonal to the tetragonal phase. Consequently, the cubic phase is necessary to analyse the stability of the other two phases and understand the experiments. For Mn 3 Ga and Mn 3 Ge, the cubic phase lies between the hexagonal phase with the highest energy and the tetragonal phase with the lowest energy. Consequently, it is possible to transform the hexagonal phase, through the cubic phase, into the tetragonal phase. This is consistent with experimental 5
observations. However, for Mn 3 Sn, the cubic phase has a higher energy than the hexagonal phase. Although the tetragonal phase has the lowest energy, it turns out to be difficult to transform the hexagonal phase into the tetragonal one. This agrees with the fact that only hexagonal Mn 3 Sn has been synthesized in experiments so far. We also propose that external pressure can be utilized to assist the phase transition. In addition, the weak ferromagnetism of hexagonal compounds and the perpendicular magnetocrystalline anisotropy in the tetragonal compounds are consistent with previous experiments. Acknowledgements We acknowledge the funding support by the ERC Advanced Grant (291472) and the kind help from Dr. J. Karel. 6
References 1 E. Krén, and G. Kádár, Solid State Communications 8, 1653 (1970). 2 G. J. Zimmer and E. Krén, Magnetism and Magnetic Materials, 1971 (AIP Conf. Proc. No. 5, American Inst. of Physics, New York, 1972) p.513. 3 G. Kádár and E. Krén, Int. J. Magn. 1, 143 (1971). 4 S. Tomiyoshi, Y. Yamaguchi, and T. Nagamiya, J. Magn. Magn. Mater. 31, 629 (1983). 5 S. Tomiyoshi and Y. Yamaguchi, J. Phys. Soc. Japan 51, 2478 (1982). 6 P. J. Brown, V. Nunez, F. Tasset, J. B. Forsyth and P. Radhakrishna, J. Phys.: Condens. Matter. 2, 9409 (1990). 7 P. Radhakrishna, and J. W. Cable, J. Magn. Magn. Mater. 104-107, 1065 (1992). 8 J. W. Cable, N. Wakabayashi, and P. Radhakrishna, Phys. Rev. B 48, 6159 (1993). 9 J. Kübler, A. R. William, C. B. Sommers, Phys. Rev. B 28, 1745 (1983). 10 S. Wurmehl, H. C Kandpal, G. H. Fecher, and C. Felser, J. Phys.: Condens. Matter 18, 6171 (2006). 11 J. Winterlik, S. Chadov, A. Gupta, V. Alijani, T. Gasi, K. Filsinger, B. Balke, G. H. Fecher, C. A. Jenkins, F. Casper, J. Kübler, G. D. Liu, L. Gao, S. S. Parkin, C. Felser, Advanced Materials, (2012) doi: 10.1002/adma.201201879. 12 S. Ikeda, K. Miura, H. Yamamoto, K. Mizunuma, H. D. Gan, M. Endo, S. Kanai, J. Hayakawa, F. Matsukura, H. Ohno, Nature Mater. 9, 721 (2010). 13 A. D. Kent, Nature Mater. 9, 699 (2010). 14 H. Kurt, K. Rode, M. Venkatesan, P. Stamenov, and J. M. D. Coey, Phys. Status Solidi B 248, 2338 (2011). 15 F. Wu, S. Mizukami, D. Watanabe, H. Naganuma, M. Oogane, Y. Ando, and T. Miyazaki, Appl. Phys. Lett. 94, 122503 (2009). 16 F. Wu, E. P. Sajitha, S. Mizukami, D. Watanabe, T. Miyazaki, H. Naganuma, M. Oogane, and Y. Ando, Appl. Phys. Lett. 96, 042505 (2010). 17 S. Mizukami, F. Wu, A. Sakuma, J. Walowski, D. Watanabe, T. Kubota, X. Zhang, H. Naganuma, M. Oogane, Y. Ando, and T. Miyazaki, Phys. Rev. Lett. 106, 117201 (2011). 18 H. Kurt, K. Rode, M. Venkatesan, P. Stamenov, and J. M. D. Coey, Phys. Rev. B 83, 020405 (2011). 19 T. Kubota, Y. Miura, D. Watanabe, S. Mizukami, F. Wu, H. Naganuma, X. Zhang, M. Oogane, M. Shirai, Y. Ando, and T. Miyazaki, Appl. Phys. Express 4, 043002 (2011). 20 H. Naganuma, M. Oogane, A. Sakuma, Y. Ando, S. Mizukami, T. Kubota, F. Wu, X. Zhang, and T. Miyazaki, Phys. Rev. B 85, 014416 (2012). 7
21 H. Kurt, N. Baadji, K. Rode, M. Venkatesan, P. Stamenov, S. Sanvito, and J. M. D. Coey, Appl. Phys. Lett. 101, 132410 (2012). 22 T. Ohoyama, K. Yasuköchi and K. Kanematsu, J. Phys. Soc. Japan 16, 352 (1961). 23 B. Balke, G. H. Fecher, J. Winterlik, and C. Felser, Appl. Phys. Lett. 90, 152504 (2007). 24 J. Winterlik, B. Balke, G. H. Fecher, C. Felser, M. C. M. Alves, F. Bernardi, and J. Morais, Phys. Rev. B 77, 054406 (2008). 25 T. Ohoyama, J. Phys. Soc. Japan 16, 1995 (1961). 26 G. Kresse, and J. Furthmuller, Phys. Rev. B 54, 11169 (1996). 27 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 28 G. Kresse, and D. Joubert, Phys. Rev. B 59, 1758 (1999). 29 P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, and J. Luitz, 2001 WIEN2k, An Augmented PlaneWave+Local Orbitals Program for Calculating Crystal Properties, Karlheinz Schwarz, Techn. Universitaet Wien, Wien, Austria. 30 J. Sticht, K. H. Höck, and J. Kübler, J. Phys.: Condens. Matter. 1, 8155 (1989). 31 L. M. Sandratskii and J. Kübler, Phys. Rev. Lett. 76, 4963 (1996). 32 R. Laskowski, and G. Santi, arxiv:cond-mat/0408200. 8
Figure 1 (D. L. Zhang et al.) Figure 1 (Colour online) Crystal structure relationships of (a) cubic-to-tetragonal and (b) cubic-tohexagonal phases. The projection from the cubic to hexagonal structures is shown in (c)-(g). Four atomic layers (three Mn layers and one Z layer) will be projected into a single layer, in order to form the hexagonal structure. The orange arrows in (a)-(b) indicate the direction of structural deformation, and the arrows on the Mn atoms in (a) denote the direction of magnetic moment. In addition, Mn atoms are represented by red or blue balls, and Z (Z=Ga, Sn and Ge) atoms by green balls. 9
Figure 2 (D. L. Zhang et al.) Figure 2 (Colour online) A schematic of possible magnetic configurations in the hexagonal lattice of Mn 3 Z (Z=Ga, Sn and Ge). The arrows denote the direction of the magnetic moments. The magnetic moments of Mn atoms lie in the hexagonal basal plane in (a)-(h), and out-of-plane in (i)-(j). The blue and red balls represent the Mn atoms in different planes. The magnetic configuration of (d) is found to be the most stable by our calculations. 10
Figure 3 (D. L. Zhang et al.) E total (ev) E total (ev) E total (ev) -29.4-29.7-30.0-30.3-30.6-29.7-30.0-30.3-30.6-30.9-31.2-31.5-31.8-32.1 (a) Mn 3 Ga cubic tetragonal hexagonal ΔE (tetra-hex) =0.27 ev 85 90 95 100 105 110 115 120 125 130 (b) Mn 3 Sn (c) cubic tetragonal hexagonal ΔE (tetra-hex) =0.11 ev 95 100 105 110 115 120 125 130 Mn 3 Ge cubic tetragonal hexagonal ΔE (tetra-hex) =0.24 ev 85 90 95 100 105 110 115 120 125 Volume (Å 3 ) Figure 3 (Colour online) The dependence of the total energy per Mn 3 Z unit on volume for the cubic, tetragonal and hexagonal structures of (a) Mn 3 Ga, (b) Mn 3 Sn, and (c) Mn 3 Ge. For a given volume, the shape of the lattice (c/a) is also fully optimized. 11
Table I. Optimized lattice parameters of hexagonal, cubic, and tetragonal structures in Mn 3 Z (Z=Ga, Sn and Ge) Heusler compounds. Lattice constant (Å) Mn 3 Ga Mn 3 Sn Mn 3 Ge Hexagonal a=5.26, c=4.26 (a=5.400, c= 4.353 exp.[1] ) a=5.57, c=4.43 (a=5.665, c= 4.531 exp.[2] ) a=5.28, c=4.22 (a=5.360, c=4.320 exp.[3] ) Cubic a=5.82 (a=5.823 calc.[10] ) a=6.04 a=5.75 Tetragonal a=3.77, c=7.16 (a=3.77, c=7.16 calc.[23] a=3.909, c=7.098 exp.[24] ) a=3.93, c=7.47 a=3.75, c=7.12 (a=3.81, c=7.26 exp.[3, 21] ) Table II. Total magnetic moments per Mn 3 Z unit of hexagonal, cubic, and tetragonal structures. For the hexagonal phase, the µ Mn denotes the magnetic moment of each Mn atom. The energy of perpendicular magnetocrystalline anisotropy (per primitive unit cell, Mn 6 Z 2 ) of the tetragonal structure is shown. Magnetic Moment (µ B ) Mn 3 Ga Mn 3 Sn Mn 3 Ge Hexagonal 0.03 (0.045 exp.[1] ) 0.01 (0.009 exp.[2] ) 0.01 (0.06 exp.[3] ) µ Mn =2.5 (2.4±0.2 exp.[1] ) µ Mn =2.9 (3.0 exp.[2] ) µ Mn =2.5 (2.4±0.2 exp.[3] ) Cubic -0.01(-0.01 calc.[10] ) -1.00-1.00 Tetragonal 1.78 (1.77 calc.[23] ) 1.04 0.97(1.00 exp.[21] ) PMA=1.0 mev (1.0 mev calc.[11] ) PMA=1.1 mev PMA=1.0 mev (0.8 mev calc.[21] ) 12