First-principles study of lattice instabilities in the ferromagnetic
martensite Ni$_2$MnGa

Show preview

Hide preview

ar X
iv :c
on d-
m at
/0 30
43 49
v1 [
co nd
-m at.
mt rl-
sc i]
15 A
pr 20
03 First-principles study of lattice instabilities in the ferromagnetic martensite Ni2MnGa
Claudia Bungaro∗ and K. M. Rabe Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854-8019, USA.
A. Dal Corso Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Beirut 2/4, 34014 Trieste, Italy.
The phonon dispersion relations and elastic constants for ferromagnetic Ni2MnGa in the cubic and tetragonally distorted Heusler structures are computed using density-functional and density- functional perturbation theory within the spin-polarized generalized-gradient approximation. For 0.9 < c/a < 1.06, the TA2 tranverse acoustic branch along [110] and symmetry-related directions displays a dynamical instability at a wavevector that depends on c/a. Through examination of the Fermi-surface nesting and electron-phonon coupling, this is identified as a Kohn anomaly. In the parent cubic phase the computed tetragonal shear elastic constant, C′=(C11−C12)/2, is close to zero, indicating a marginal elastic instability towards a uniform tetragonal distortion. We conclude that the cubic Heusler structure is unstable against a family of energy-lowering distortions produced by the coupling between a uniform tetragonal distortion and the corresponding [110] modulation. The computed relation between the c/a ratio and the modulation wavevector is in excellent agree- ment with structural data on the premartensitic (c/a = 1) and martensitic (c/a = 0.94) phases of Ni2MnGa.
1

I. INTRODUCTION
Ferromagnetic shape-memory alloys displaying large magnetic-field-induced strain have recently emerged as a new class of active materials, very promising for actu- ator and sensor applications. The largest known mag- netostrain effects have been observed in Ni2MnGa-based Heusler alloys, where up to 6% and 9.5% strains have been induced by a magnetic field less than 1 Tesla.1,2
The magnetic shape-memory behavior is closely linked to the occurrence of a martensitic transformation in con- juction with a strong magnetocrystalline anisotropy of the low-temperature martensitic phase. For a deeper un- derstanding of the magnetic shape-memory mechanism, a microscopic explanation for the origin of the martensitic transformation and the magnetocrystalline anisotropy is clearly needed. Towards this goal, we have used ab ini- tio techniques to investigate the origin of the martensitic transformation. Ni2MnGa is the most intensively studied of the rel-
atively few known ferromagnetic shape-memory materi- als. A number of thermal and stress-induced martensitic transformations have been observed in Ni2MnGa-based alloys, and both transition temperatures and crystallo- graphic structures are quite sensitive to alloy composi- tion. In the stoichiometric alloy, Ni2MnGa, which is fer-
romagnetic below TC ≈ 380 K, two thermally-induced phase transitions have been observed. From the high- temperature cubic Heusler structure, a premartensitic phase transformation to a modulated cubic structure oc- curs below TPM ≈ 260 K,
3,4 followed by a martensitic transformation to a modulated tetragonal structure be- low TM ≈ 220 K.
5 The low temperature martensitic phase has a tetragonal structure (c/a=0.94, a=5.90A˚) with a superimposed incommensurate modulation along the [110] direction, consisting of a shuffling of (110) planes in the [110] direction with a periodicity of al- most 5 interplanar distances, which corresponds to a wave vector qM ≈
2π a (0.43, 0.43, 0) (Refs. 4, 6). The
low-temperature martensitic phase has a high magne- tocrystalline anisotropy,7 making it useful for applica- tions. At a phenomenological level, the coupling between strain, modulation, and magnetization has recently been described in a Landau-theory framework.8
Some clues to the microscopic origin of the martensitic transitions are provided by the softening of particular phonons and elastic constants. Inelastic neutron scat- tering experiments on the high-temperature phase found a significant, though incomplete, softening in the TA2 phonon branch along the [110] direction, at a wave vec- tor qPM ≈
2π a (0.33, 0.33, 0) (Ref. 3). The phonon soft-
ening has been found to correspond to a premartensitic phase occurring between the high-temperature cubic and low-temperature martensitic structures. The premarten- sitic phase has a cubic structure with a superimposed [110]-transverse shuffling modulation, analogous to the
modulation of the martensitic phase but with a different periodicity of almost 6 interplanar distances.3 While the PM phase is anticipated by a precursor phonon softening at qPM, no phonon softening at qM has been observed above TM. The aim of this study is to provide a unified ex-
planation for the microscopic origin of the rich vari- ety of phase transitions and modulated structures oc- curring in this ferromagnetic shape-memory alloy. To this end we have performed a first-principles study of the phonon dispersions and lattice instabilities of ferromag- netic Ni2MnGa and of their dependence upon uniform tetragonal strain. In our calculations we use the spin- polarized generalized-gradient approximation (σ-GGA) recently developed within the density-functional pertur- bation theory (DFPT) formalism.9 This approach allows us to obtain an accurate description of the structural, magnetic, and vibrational properties of Ni2MnGa. The paper is organized as follows. The details of the
computational method are given in section II. In section III we present the results for the cubic Heusler structure. The crystal structure and magnetization are discussed in Section III.A. Section III.B is devoted to the phonon dis- persion, with particular attention to the phonon anomaly and related dynamical instability. The origin of the phonon anomaly is discussed in Section III.C. The com- puted elastic constants are presented and discussed in Section III.D. In Section IV we investigate the depen- dence of the phonon dispersions and lattice instabilities on uniform tetragonal strain. In Section V the main re- sults are summarized and the microscopic origin of the observed phase transitions explained.
II. COMPUTATIONAL METHOD
Our calculations have been performed within the framework of density-functional and density-functional- perturbation theory (DFT and DFPT). In particular, the vibrational properties have been computed using a recent implementation of ultrasoft pseudopotentials into DFPT.10 We used ultrasoft pseudopotentials11 for Ni and Mn, freezing the 3s and 3p core electrons and treating the 3d and 4s states as valence levels.12,13 The nonlinear core correction is used to account for the overlap between the core and the valence charges.14 For the Ga atom we used a norm-conserving pseudopotential, which also includes a nonlinear core correction, treating the 4s and 4p states as valence levels. To describe the effects of exchange and correlation, we used the Perdew-Burke-Ernzerhof15
functional with the spin-polarized generalized gradient correction (σ-GGA), recently implemented within the DFPT formalism in the PWscf code.9 For comparison, we have performed calculations within the local spin density approximation (LSDA) using the Perdew-Zunger param- eterization of the exchange and correlation energy. The plane-wave basis set had a kinetic-energy cutoff of 25
2

Ry. The augmentation charges, required by the use of ultrasoft pseudopotentials, were expanded with an en- ergy cutoff of 450 Ry. The Brillouin zone integration was performed using the smearing technique,16 suitable for metallic systems. The structural properties and most of the phonon frequencies are well converged using a first-order smearing function with a smearing parameter σ=0.03 Ry and an fcc (6, 6, 6) Monkhorst-Pack grid,17
yielding 28 k-points in the irreducible wedge of the Bril- louin zone (IBZ) for the cubic structure. To obtain an ac- curate description of the anomalous TA2 phonon branch, a smaller smearing parameter is necessary and, conse- quently, a finer k-point sampling. An accuracy to within a few cm−1 has been obtained using σ=0.01 Ry and an fcc (10, 10, 10) k-point grid, yielding 110 k-points in the IBZ, for the cubic structure. For the tetragonal struc- ture, an (8, 8, 8) k-point grid yielding 144 k-points in the IBZ of the face-centered orthorhombic unit cell was used. To compute the full phonon dispersions of the cu-
bic structure, we computed the interatomic force con- stants by Fourier transformation of the dynamical matri- ces computed on a (6, 6, 6) q-point grid in the fcc BZ. The phonon dispersions along the [110] ([011]) direction were obtained by interpolating the dynamical matrices computed on a finer mesh of 24 q-points between Γ and the shortest reciprocal-space vector parallel to the [110] ([011]) direction. An even denser mesh of q-points has been used in the proximity of the anomaly for a more accurate interpolation of the TA2 branch.
III. CUBIC Ni2MnGa
In this section, we report the results of calculations for Ni2MnGa in the fcc L21 Heusler structure (see figure 1).
A. Crystal structure and magnetization
In table I are given the minimum-energy lattice param- eter, a0, the bulk modulus, B0, and the magnetic moment per unit cell, µ0, computed within the spin-polarized σ- GGA. Note that to obtain convergence to within 0.01µB, µ0 has been computed on a (10, 10, 10) k-point grid with σ=0.01 Ry. These results agree well with those of all-electron calculations (FLAPW) performed using the σ-GGA.18 The theoretical results are also in very good agreement with the experimental data; the theoret- ical lattice parameter obtained with σ-GGA is equal, to within theoretical precision, to the experimental value. For comparison, the LSDA results are included in ta- ble I. The LSDA lattice parameter is 2.5% smaller than the experimental value. This underestimate of the lat- tice parameter is correlated with an increase in the bulk modulus to 38% higher than the experimental value.
B. Phonon dispersion relation
The phonon dispersions of the ferromagnetic cubic structure (a0 = 11.03 a.u.) have been computed, using σ-GGA, for ~q along high-symmetry lines in the first Bril- louin zone, as shown in figure 2. The solid lines indicate the computed phonon dispersion converged to within a few wave numbers. The theoretical phonon dispersion curves are in excellent agreement with the available in- elastic neutron scattering data.4,19
To evaluate the accuracy of the theoretical method, the phonon frequencies computed using different approxima- tions for the exchange and correlation energy, at the zone boundary X point, are given in table II. Since the experi- mental frequency of the anomalous (ζζ0)-acoustic branch is not strongly temperature dependent at this q-point, it is a good reference for evaluating the accuracy of the theoretical method. The experimental frequency for the transverse acoustic mode X′5 is in very good agreement with the phonon frequency computed using the σ-GGA, with the computed value being about 6% softer. For com- parison, we have also done the computation with LSDA. Calculations for the cubic structure with the LSDA lat- tice constant yield frequencies about 10% harder than those computed using the σ-GGA, except for the lowest X′5 mode which is 7% softer, increasing the discrepancy with experiment. In addition, the ordering of the modes in LSDA is different. Specifically, the X′4 mode is softer than the nearby X′5 and X5 modes. We have also com- puted the phonon frequencies within the LSDA but fixing the lattice parameter to the equilibrium value computed in the σ-GGA. The result is that all the modes soften, so that the frequencies are about 10% smaller than with σ-GGA. This softening results in an even greater discrep- ancy with the experimental X′5 mode (20%). We conclude that the better accuracy of σ-GGA is not merely an effect of the more accurate value of the equilibrium lattice pa- rameter. All calculations reported below were performed with σ-GGA. The most striking feature of the phonon dispersion
relation in figure 2 is the anomalous dip in the low- est branch of the transverse acoustic modes (TA2) along [110]. The minimum at the incommensurate wave vec- tor q0=
2π a (ζ0, ζ0, 0), with ζ0 = 0.34, is at imaginary
frequency. Thus, the crystal is dynamically unstable to the lattice distortion corresponding to the eigenvector of this mode, which consists of a nearly rigid displacement of the (110) atomic planes along the [110] direction with a modulation period of slightly less than 6 interplanar distances along the [110] direction. This energy-lowering distortion can be specified by an amplitude, u, and a phase, φ,
um = u cos(mζ0π + φ), (1)
where um is the displacement of the m th (110) atomic
plane along the [110] direction. Due to the cubic symme- try, there are 12 equivalent anomalies along the <110>
3

directions which give rise to six different but crystallo- grafically equivalent lattice modulations um. The phonon anomaly is shown in more detail in fig-
ure 3, where the theoretical dispersion is compared with inelastic neutron scattering data taken at two different temperatures: T=370 K (squares) and T=250 K (trian- gles).3 The experimental data display an anomaly in the TA2 branch that corresponds to the anomaly predicted by our theoretical dispersion. The wave vector of the experimentally-observed anomaly is ζPM ≈ 0.33,
3 in ex- cellent agreement with the theoretical value. For those modes that do not have a strong dependence upon tem- perature, such as the LA branch and the TA2 modes away from the anomaly (ζ ≥ 0.45), the agreement between theory and experiment is very good. To compare theory and experiment for modes that are strongly temperature dependent, we need to extrapolate the experimental val- ues to T = 0 K. In accordance with the soft mode the- ory, (h¯ω)2 is experimentally observed to decrease linearly with temperature above TPM=260 K. Extrapolation to T = 0 K gives an imaginary frequency of 30i cm−1, which is in reasonable agreement with the computed value. Examination of the phonon dispersion throughout the
entire BZ shows that imaginary frequencies leading to dy- namical instabilities occur only in a very localized region in q-space. This is in agreement with the experimen- tal dispersion measured in the direction perpendicular to q0, q=q0+(ζ,−ζ, 0).
4 The contour plot in figure 4 shows that the lattice instability is confined to a “drop- shaped” small region in reciprocal space, located along the [110] direction and centered at the critical wave vec- tor q0, where the dominant instability occurs. The lo- calized nature of the anomaly in q-space is a signature of its electronic origin, as will be discussed in detail in the next section.
C. Origin of the anomaly
Screening due to electron-phonon coupling involving electronic states near the Fermi level in metals can give rise to anomalous dips in the phonon dispersion, called Kohn anomalies. The occurrence of these anomalies de- pends mainly on the geometry of the Fermi surface, as well as on the q-dependence of electron-phonon matrix elements. If the Fermi surface has flat portions with nest- ing vector q0, there generally will be very strong screen- ing of the potential perturbation due to atomic displace- ments at that wavevector, leading to a pronounced soft- ening highly localized in q-space. The effect of Fermi surface geometry can be quantified by calculation of the generalized susceptibility, as in Ref. 20. In cubic Ni2MnGa, the Fermi level crosses both minor-
ity and majority spin bands. Of the two Fermi surfaces, only the one for the minority spin bands, plotted in fig- ure 5, shows obvious nesting features. In particular, a large fraction of the opposite sides of the “pipes,” running
along the faces of the cube and crossing at their center, are nested by wavevectors of the form (ζ0, k, 0)
2π a with ζ0
= 0.34 and 0 ≤ k ≤ 0.34 and the cubic-symmetry-related equivalent wavectors. Along [110], this corresponds to the critical wave vector of the anomaly. It should be noted that this differs from the position of the nesting wavevector and peak in generalized susceptibility along [110] found in Refs. 20 and 21. This difference might be attributable to the differences in method, in particular, the difference in choice of density functional. [110] ap- pears to be a direction for which the anomaly is strongest. For example, although there is nesting along [100] as well, there is no sign of an anomaly in the phonon dispersion. This can be attributed to the vanishing of the relevant electron-phonon matrix elements. The strength of a Kohn anomaly is expected to be very
sensitive to electronic temperature, with an increasing temperature reducing the sharpness of the Fermi surface and thus weakening the anomaly. In our calculations, the sensitivity to electronic temperature can be directly investigated, as the smearing parameter σ plays the role of a fictitious electronic temperature. In table III we show the dependence on σ of those modes at the critical wave vector, q0, that have the same symmetry as the soft TA2 branch (i.e. the three modes of Σ3 symmetry). The anomalous TA2 mode, Σ3(1), is much more sensitive to the fictitious electronic temperature than the other two modes. All the frequencies are well converged for σ=0.01 Ry. The dependence of the phonon anomaly upon the electronic temperature is also shown in figure 2, where the dispersion of the anomalous TA2 branch computed with σ=0.03 Ry (dashed line) is compared with the fully converged calculation for σ=0.01 Ry (solid line). Only the modes with wavevectors close to the anomaly, (ζ0 − 0.16) < ζ < (ζ0 + 0.16), are affected by the change in the electronic temperature. The anomaly is smoothed out as the fictitious electronic temperature is increased. Thus, we conclude that the computed dip is indeed a Kohn anomaly.
D. Elastic constants
The computed elastic constants for the ferromagnetic cubic structure are reported in table IV. In a first- principles framework, the standard approach for obtain- ing elastic constants is to compute the stress tensor for a selected set of small strains. To obtain convergence at the level of a few GPa, we have used very dense k- point meshes (up to (14 14 14), corresponding to 10976 points in the BZ). To compute C11 and C12 the cubic crystal is distorted by the tetragonal deformation: ǫzz=ǫ, ǫxx=ǫyy=0. For small deformations the stress, σij , is lin- ear with ǫ: σxx=σyy=C12ǫ, and σzz=C11ǫ. The values so obtained for the elastic constants C11, C12, and the shear modulus C′=(C11−C12)/2 are indicated as “theory 1” in table IV.
4

The small magnitude of C′ requires additional atten- tion. We performed an independent computation by con- sidering the tetragonal deformation ǫxx=ǫyy=ǫ, ǫzz=−2ǫ. For small ǫ we have that σxx=σyy=2C
′ǫ and σzz=−4C ′ǫ.
This direct calculation gives C′=(2±2) GPa, in good agreement with the value obtained as a difference of C11 and C12 in theory 1. The near-zero value of C
′ implies that there is almost no energy cost for a small tetragonal shear distortion of the type 2ǫxx=2ǫyy=−ǫzz. Evidence of this marginal elastic instability towards small volume- preserving tetragonal distortions was also found in Ref. 22, where it was shown that the energy surface as a func- tion of c/a is remarkably flat; specifically, the change in energy associated with varying the c/a ratio in the range of values between 0.97 and 1.01 is almost zero to within numerical accuracy. In table IV (theory 2) we also show the values of the
elastic constants estimated from the slope of the long wavelength acoustic modes along the [110] direction. The elastic constants C44, C
′, and CL=(C11+C12+2C44)/2 correspond to the TA1, TA2, and LA modes, respectively. These determine the values of C11 and C12 given in the table. The evaluation of the elastic constants from the phonon dispersion is less accurate, especially for C′ where the corresponding TA2 branch deviates from a linear be- havior already at very small q. The elastic constants in the high-temperature cu-
bic phase have been experimentally determined between room temperature and the pre-martensitic structural phase transition at T = 260 K by measuring the veloc- ity of ultrasonic waves.23,24 These experiments show that the elastic constants CL and C11 are almost temperature independent. In contrast, the transverse elastic constants C44 and C
′ exhibit an anomalous behavior with cooling, softening as the pre-martensitic phase transition temper- ature is approached. While the softening is small for C44, it is more dramatic for C′, which decreases by 60% from room temperature to the transition. We find good agree- ment between theory and room temperature measure- ments for CL and C11, which are not strongly modified by finite temperature. The pronounced softening of C′
is also consistent with the very small zero-temperature theoretical value, C′ ≈ 0.
IV. TETRAGONAL Ni2MnGa
The identification of unstable phonons at q0= 2π a (ζ0,
ζ0, 0), with ζ0 = 0.34, in the cubic Heusler structure leads naturally to an understanding of the transition to the premartensitic cubic-modulated phase with decreas- ing temperature. To explore the subsequent transition to the low-temperature martensitic phase with c/a = 0.94 and modulation wavevector ζ0 = 0.43, we extended the calculations of phonon dispersion to Heusler structures with uniform volume-preserving tetragonal strains rang- ing from c/a = 0.88 to c/a = 1.06, with particular atten-
tion to the TA modes along the <110> directions.
A. Tetragonal structure with c/a=0.94
We first consider in detail the volume-preserving tetragonal distortion with c/a=0.94, corresponding to the low temperature martensitic phase. The phonon dis- persions computed along the inequivalent [110] and [011] directions of the tetragonal BZ are shown in figure 6. The only mode substantially affected by the tetragonal distortion is the anomalous TA2 mode. We find an over- all softening (hardening) of the TA2 branch in the [110] ([011]) direction. The anomaly is particularly affected by the tetragonal distortion. Along the <110> direc- tions, which are perpendicular to the c axis, the phonon anomaly is more pronounced and is shifted to a larger q-vector (ζ = 0.43) than in the cubic structure. This wavevector is in excellent agreement with the long pe- riod modulation observed experimentally in the tetrag- onal martensitic phase. In contrast, along the <011> directions, the anomaly has almost completely disap- peared; the entire branch is stable and there is only a very small wiggle in the TA2 branch. Thus, the entropy from low-frequency phonons should be roughly a factor of three less than in the premartensitic phase, which may explain why the premartensitic phase is more favorable at higher temperatures. As in the cubic phase, the anomaly can be associated
with features of the Fermi surface. In figure 7 we show the minority-spin Fermi surface computed for Ni2MnGa in the tetragonal structure. It is related to the Fermi surface of the cubic structure with two major differences: (i) there are no flat surfaces perpendicular to the (001) direction, and (ii) the “pipes” running along the faces of the tetragonal BZ are wider than the corresponding features in the cubic Fermi surface. As a consequence of these changes, induced by the tetragonal distortion, there is significant nesting only in the <110> directions per- pendicular to the c axis and the edges of the “pipes” are nested by the larger (ζ ζ 0) critical wave-vector, ζ = 0.43. This explains why in the tetragonal structure there are pronounced Kohn anomalies only in the <110> direc- tions, at a larger wave vector than in the cubic structure.
B. Evolution of phonon anomaly with tetragonal
deformation
To understand better the lattice instabilities and re- lated structural energetics of Ni2MnGa, we studied the evolution of the TA2 phonon anomaly with uniform volume-preserving tetragonal strain. Studying the phonon dispersions of the cubic Heusler
structure (c/a = 1) and its tetragonal distortion with c/a = 0.94, we have identified the existence of the un- stable mode at ζ = 0.34 and ζ = 0.43, respectively, that
5

explain the transitions to the premartensitic and marten- sitic phases with decreasing temperature. The fact that for the two tetragonal strains the soft mode is always the TA2 mode but with a different wave vector suggests that the wavevector of the soft mode and the tetragonal strain are coupled. We therefore studied in more detail the dependence of ζ upon tetragonal strain. The dispersion of the TA2 branch along the [110] direc-
tion is shown in figure 8 for several values of the c/a ratio. An overall softening of the TA2 branch is observed with decreasing c/a ratio. With compressive strain (c/a < 1) the anomaly occurs at a larger ζ0 than in the cubic struc- ture, and becomes broader and more pronounced. For c/a <∼ 0.91 the entire branch is unstable. With tensile strain (c/a > 1) the anomaly shifts to a smaller ζ0 and becomes less pronounced. The lattice instability is com- pletely eliminated for c/a >∼ 1.06. This is in agreement with the experimental observation of no superimposed modulations in a stress induced tetragonal phase with c/a = 1.18.25
In figure 8, we have shown the TA2 branch for several different values of c/a. For each of these values of c/a, we can find the value of the wavevector ζ0 for which the imaginary frequency of the TA2 branch has its minimum. This gives us the dependence of the modulation wavevec- tor upon strain and is shown in figure 8. Three regimes can be identified: (i) c/a <∼ 0.91, where the entire TA2 branch is unstable and there is no unique minimum; (ii) 0.9 < c/a < 1.06, where the TA2 branch displays a well- defined dynamical instability localized at the wavevector ζ0 that depends on c/a; and (iii) c/a >∼ 1.06, where the lattice instability is completely eliminated and the TA2 branch is stable over the whole BZ, so that no minimum is defined. Only in regime (ii) can a modulated structure be expected to occur. In order to quantify these trends, in figure 10 we de-
compose the phonon frequency ω2 into a short-range “normal” part ω2n and a long-range “anomalous” part ω
2 a:
ω2(ζ)=ω2n(ζ) + ω 2 a(ζ). The “normal” contribution, de-
fined as ω2n(ζ)=A[1− cos(ζπ)], depends only upon short- range interatomic force constants and is connected to the local chemistry of the crystal. It corresponds to the sim- ple model for the transverse [110] vibrations of the [110] planes when they interact only with a first-neighbor in- terplanar force constant: K=-Am/2, wherem is the mass associated with each plane. The long-range “anomalous” contribution ω2a depends on electronic screening effects that can be strongly q dependent. To fit the constant A in the definition of ω2n(ζ), we assume that for ζ=1 the anomalous contribution is zero. The short-range contribution softens with decreasing
c/a ratio, causing the overall softening observed for the TA2 branch. Decreasing c/a corresponds to a volume- preserving tetragonal distortion with a shorter c and larger a. Therefore smaller c/a values correspond to a larger [110] interplanar distance which leads to a weaker interplanar force constant and therefore a softer ω2n. The anomalous contribution ω2a is negative, large, and
peaked at a critical wave vector ζ0, indicating a strong screening due to the electon-phonon coupling and Fermi surface nesting. Both ζ0 and the intensity of the anomaly are strongly dependent upon strain. The value of ζ0 in- creases with decreasing c/a, the intensity reaching its maximum at c/a=0.97.
V. SUMMARY AND CONCLUSIONS
Our ab initio study of the phonon dispersion and lat- tice instabilities in cubic and tetragonal Ni2MnGa can be summarized as follows. First, the parent cubic phase exhibits a marginal elastic
instability, C′ ≈ 0 , meaning that there is almost no energy cost for a small uniform tetragonal distortion. Second, for 0.91 < c/a < 1.06 a Kohn anomaly is
present due to electron-phonon coupling and Fermi sur- face nesting. It develops into a deep minimum corre- sponding to a localized dynamical instability whose wave vector ζ0 is related in a one-to-one fashion to the c/a ratio. Therefore, tetragonal structures with c/a in this range are unstable towards a particular transverse (110)- shuffling modulation with a specific wavevector ζ0, the one-to-one relationship of the c/a ratio and the modula- tion period being due to the behavior of the topology of the Fermi surface under tetragonal distortions. Consequently, the energy of the cubic structure may be
lowered by any one of a family of deformations, {ǫ,u(ζǫ)}, consisting of a tetragonal strain, ǫ, and a superimposed (110)-shuffling modulation, u(ζǫ), whose periodicity, ζ, is directly related in a one-to-one fashion to the c/a ratio. The near-vanishing of C′ together with the fact that
the energy surface as a function of c/a is remarkably flat permits a modulated ground state with c/a significantly different from 1. Indeed, the computed relation between c/a and the modulation wavevector is in excellent agree- ment with structural data on the premartensitic (c/a = 1) and martensitic (c/a = 0.94) phases of Ni2MnGa, and also with the unmodulated stress-induced phase (c/a = 1.18). The next step is to carry out total energy calculations
of the complex modulated structures in order to deter- mine the equilibrium amplitude and phase of the lattice modulation and the associated energy gain. The energy could then be expanded in terms of symmetry invariants of the relevant structural degrees of freedom (strain and modulation) and determine the coefficients of each term of the expansion using ab initio calculations. In this way one could develop an ab-initio based finite-temperature statistical model of the phase transition in Ni2MnGa. The three phases of Ni2MnGa observed with increas-
ing temperature can be understood as follows. At low temperatures, a minimum-energy modulated tetragonal structure is expected, consistent with the observed c/a = 0.94 phase. As temperature increases, the favorable entropy arising from the larger phase space for the low-
6

frequency phonons in the cubic structures leads to a phase transition to the cubic-modulated structure, with a modulation wavevector shifted down from that in the tetragonal phase in accordance with our results. Finally, the frequencies of the unstable phonons associated with the modulation are anharmonically renormalized at the highest temperatures and the observed structure is the cubic Heusler structure, with large fluctuating local dis- tortions. In conclusion, we have provided a first-principles mi-
croscopic explanation for the origin of the physically interesting and technologically important martensitic phase transitions in the ferromagnetic shape-memory Heusler alloy Ni2MnGa.
Calculations in this work have been done using the PWscf package.26 We thank Morrel Cohen, Richard James and Xiangyang Huang for valuable discussions, and Vitaly Godlevsky for his help in the early stages of this work.
∗ Current email: bungaro@physics.rutgers.edu 1 S. J. Murray, M. A. Marioni, A. M. Kukla, J. Robinson, R. C. O’Handley, and S. M. Allen, J. Appl. Phys. 87, 5774 (2000).
2 A. Sozinov, A. A. Likhachev, N. Lanska, and K. Ullakko, Appl. Phys. Lett. 80, 1746 (2002).
3 A. Zheludev, S. M. Shapiro, P. Wochner, A. Shwartz, M. Wall, and L. E. Tanner, Phys. Rev. B 51, 11310 (1995).
4 A. Zheludev, S. M. Shapiro, P. Wochner, and L. E. Tanner, Phys. Rev. B 54, 15045 (1996).
5 P. J. Webster, K. R.A. Ziebeck, S. L. Town, and M. S. Peak, Philos. Mag. 49, 295 (1984).
6 V. V. Martynov and V. V. Kokorin, J. Phys. (France) III 2, 739 (1992).
7 R. Tickle and R. D. James, J. Magn. Mater 195, 627 (1999).
8 V. Buchelnikov, A. Zayak, A. Vasilev and T. Takagi, Int. J. Appl. Electromag. and Mech. 12, 19 (2000).
9 A. Dal Corso and S. de Gironcoli, Phys. Rev. B 62, 273 (2000).
10 A. Dal Corso, A. Pasquarello, and A. Baldereschi, Phys. Rev. B 56, 11 369 (1997).
11 D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 12 The details for the Ni pseudopotential can be found in Ref. 9.
13 We generated the pseudopotential for Mn, using as refer- ence the all-electron configuration 3d5, 4s2. The core radii are (in a.u.): 4s (2.0, 2.2), 4p (2.2, 2.4), 3d (1.6, 2.2). The all electron potential pseudized inside rc=1.8 has been taken as a local potential. Two values of the core radii indi- cate a channel which has been pseudized with the ultrasoft scheme. In such a case the first value is the norm-conserving core radius and the second is the ultrasoft one.
14 S. Louie, S. Froyen and M. L. Cohen, Phys. Rev. B 26, 1738 (1982).
15 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
16 M. Methfessel and A. T. Paxton, Phys. Rev. B 40, 3616 (1989).
17 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976).
18 A. Ayuela, J. Enkovaara, K. Ullakko, and R. M. Nieminen, J. Phys. Condens. Matter 11, 2017 (1999).
19 A. Zheludev, S.M. Shapiro, P. Wochner, A. Shwartz, M. Wall, and L. E. Tanner, J. Phys. IV Colloq. 5, C8-1139 (1995).
20 Y. Lee, J. Y. Rhee and B. N. Harmon, Phys. Rev B 66, 054424 (2002).
21 O. I. Velikokhatny and I. I. Naumov, Phys. Solid State 41, 617 (1999).
22 See Figure 1 in: V. V. Godlevsky and K.M. Rabe, Phys. Rev. B 63, 134407 (2001).
23 J. Worgull, E. Petti, and J. Trivisonno, Phys. Rev. B 54, 15695 (1996).
24 L. Man˜osa, A. Gonza`lez-Comas, E. Obrado´, A. Planes, V. A. Chernenko, V. V. Kokorin, and E. Cesari, Phys. Rev. B 55, 11068 (1997).
25 V. V. Martynov, J. Phys. (France) IV C8, 91 (1995). 26 S. Baroni, A. Dal Corso, S. de Gironcoli, and P. Giannozzi, http://www.pwscf.org.
27 P. J. Webster, Contemp. Phys. 10, 559 (1969).
7

TABLE I. Equilibrium lattice parameter, a0, bulk modu- lus, B0, and magnetic moment, µ0, for Ni2MnGa in the L21 Heusler structure. Calculations using the spin-polarized GGA (σ-GGA), are compared with experimental data and with the results obtained using LSDA.
a0(a.u.) B0(Mbar) µ0(µB)
σ-GGA 11.03 155 4.27 LSDA 10.74 202 3.92 FLAPW-σ-GGAc 10.98 156 4.09
Exp. 11.01a 146b 4.17a
aRef. 27 bRef. 23 cRef. 18
TABLE II. Phonon frequencies (in cm−1) computed at the zone-boundary point X. Results obtained using different ap- proximations for the exchange and correlation energy, σ-GGA and LSDA, are compared with the experiment.
σ-GGA LSDA LSDA Exp. a0 a0(σ-GGA) a0(LSDA) a0(σ-GGA)
X′5 82 76 70 87 X′1 160 179 151 X5 174 191 161 X′5 176 195 166 X′4 178 185 159 X1 191 203 168 X′1 227 249 212 X′5 259 286 233
TABLE III. Dependence upon the smearing parameter, σ, of the frequency of the three Σ3-modes at the critical wavevec- tor, q0. The Σ3(1) mode is the anomalous TA2 mode, the other two, Σ3(2) and Σ3(3), are higher frequency optical modes. The k-point mesh needed to achieve convergence for each of the values of σ is also reported. Frequencies are in cm−1.
σ(Ry) k-mesh Σ3(1) Σ3(2) Σ3(3)
0.03 (6 6 6) 8 195 244 0.01 (10 10 10) 39i 195 242 0.005 (12 12 12) 39i 195 242
TABLE IV. Elastic constants, in GPa, for ferromagnetic Ni2MnGa in the cubic Heusler structure. Theory 1, elas- tic constants from stress calculation under strain. Theory 2, elastic constants estimated from the slope of the [ζζ0] acous- tic branches close to Γ: C′, C44, and CL correspond to the TA2, TA1, and LA modes, respectively. Room temperature measurements are shown for comparison.
C11 C12 C ′ C44 CL
theory 1 153±2 148±2 2.5±2 theory 2 138±9 143±9 −2.5±5 100±5 240±5 Exp.a 152 143 4.5 103 250
Exp.b 136±3 – 22±2 102±3 222±9
aRef. 23 bRef. 24
FIG. 1. The fcc L21 Heusler structure of Ni2MnGa. The small circles, large open circles and large filled circles repre- sent Ni, Mn and Ga, respectively.
8

0 0.25 0.5 (ζ ζ ζ)
10 i
0
10
20
30
E( me
V)
00.51 (ζ ζ 0)
0 0.25 0.5 0.75 1 (2ζ ζ 0)
00.51 (ζ 0 0)
100 i
0
100
200
300 fre
qu en
cy (c
m− 1 )
Γ W X Γ LKX
FIG. 2. Full phonon dispersion of ferromagnetic Ni2MnGa in the fcc Heusler structure, along high symmetry lines of the fcc BZ. Solid lines are fully converged ab initio calculations obtained using an effective electronic temperature σ=0.01Ry. The dashed line is the dispersion obtained using a higher effective electronic temperature, σ=0.03Ry, showing the dependence of the anomaly in the lowest acoustic branch along the (ζζ0) direction on the effective electronic temperature. The circles indicate the neutron scattering data from Refs. [ 19, 4]. The imaginary values of phonon frequencies are plotted along the negative frequency axis. The wave-vector coordinate ζ is in units of 2pi
a .
9

0 0.2 0.4 0.6 0.8 1 q=ζ(110) 2pi/a
50i
25i
0
25
50
75
100
125
fre qu
en cy
(c m−
1 )
theory 370oK 250oK
Γ XK
LA
TA1
TA2
FIG. 3. Partial phonon dispersion of Ni2MnGa in the fcc Heusler structure, along the Γ-K-X line in the [110] direc- tion. The theoretical data are the same as in figure 2. The experimental data taken at 250 K and 270 K are shown for comparison.
[100]
[01 0]
Γ X
X XW
W
K
[00 1]
[110] Γ q0 XK
X K Γ
FIG. 4. Contour plots of the imaginary phonon frequencies in the (001) plane (a) and (110) plane (b) of the fcc BZ. These contour plots show the region in q-space where the crystal is dynamically unstable. This occurs in a narrow valley situated along the [110] direction with its minimum at the critical wave vector q0=(0.34, 0.34, 0).
FIG. 5. Fermi surface of minority spin bands for cubic Ni2MnGa.
0 0.2 0.4 0.6 0.8 1 q=ζ(110) 2pi/a
[110]
100i
0
100
200
300
fre qu
en cy
(c m−
1 )
c/a=0.94
0 0.2 0.4 0.6 0.8 1 q=ζ(0 1/a 1/c) 2pi
[011]
c/a=0.94
FIG. 6. Ab initio [110] and [011] phonon dispersions for fer- romagnetic Ni2MnGa in the tetragonal structure (c/a=0.94, a=11.26 a.u.)
FIG. 7. Fermi surface of minority spin bands for tetragonal Ni2MnGa (c/a=0.94).
10

0 0.2 0.4 0.6 0.8 1
q=ζ[110] 2pi/a −5000
0
5000
10000 ω
2 (c m−
2 )
c/a 1.06
1.03
1
0.97 0.94
0.91 0.88
FIG. 8. (color online only). Dispersion of the squared fre- quencies, for the anomalous TA2 phonon branch, computed for different volume-preserving tetragonal distortions, with c/a ranging from 0.88 to 1.06.
0
0.2
0.4
0.6
0.8
1
1.2
ζ 0
0.8 0.9 1 1.1 1.2 c/a
2
4
6
8
10
m o
du la
tio n
pe rio
d
I II III
dynamical instability
no modulation
whole branch unstable
instability localized
at ζ 0
no
FIG. 9. Dependence upon c/a of the modulation wavevec- tor ζ0 and the related modulation period (in unit of (110) interplanar distances).
0
5000
10000
15000
ω n
2 (c m−
2 )
1 0.97 0.94 0.91
1.03
1.06 c/a
0 0.2 0.4 0.6 0.8 1 q=ζ[110] 2pi/a
−5000
−2500
0
ω a2
(cm
− 2 )
1
0.97 0.94
0.91 1.03
1.06
FIG. 10. Decomposition of ω2 into “normal” (ω2n) and anomalous (ω2a) parts for different volume-preserving tetrag- onal distortions, with c/a varying from 0.91 to 1.06. Dashed lines are for c/a > 1, dot-dashed lines are for c/a < 1 and solid lines are for c/a = 1.
11