The giant negative thermal expansion in the Ge-doped antiperovskite Mn3CuN compound is theoretically studied by using the first principles calculations. We propose that such a negative thermal expansion property is essentially attributed to the magnetic phase transition, rather than to the lattice vibration of the Ge-doped compound. Furthermore, we found that the doped Ge atoms in the compound significantly enhance the antiferromagnetic couplings between the nearest neighboring Mn ions, which effectively stabilizes the magnetic ground states. In addition, the nature of the temperature-dependent changes in the volume of the Ge-doped compound was revealed.
1. Introduction
Controlling the thermal expansion behavior of materials is eagerly required in modern advanced industries [1, 2]. Commonly, such control can be realized through mixing a kind of positive thermal expansion (PTE) materials with a kind of negative thermal expansion (NTE) materials whose volumes contract on heating. As we know, there are many PTE materials in the world, but the NTE ones rarely exist yet. Therefore, a challenging issue in controlling the thermal expansion behavior of the materials is to synthesize materials possessing NTE property.
Mn3CuN, a member of the anti-perovskite manganese nitride Mn3MN (M = Zn, Cu, etc.) [3–5], does not possess the negative thermal expansion property. However, recent experiments [6] reported that, with partial substitution of Ge for Cu in the Mn3CuN compound, the synthesized Mn3(Cu1-xGex)N compounds exhibited isotropic negative thermal expansion. Typically, the Ge-doped compound, Mn3(Cu1-xGex)N(x=0.5), showed the NTE property in a wide temperature window of 316–386 K, together with a giant negative thermal expansion coefficient (The linear thermal expansion coefficient of a material is definited as α=(1/L)(ΔL/ΔT). Here, L is the length of the material at temperature T0, and ΔL is the change of the length when temperature has a change of ΔT=T-T0) of αL=-16×10-6 K−1. More importantly, the working temperature and the width of the temperature window as well as the negative thermal expansion coefficient of the Ge-doped compound were modulated through the introduction of Ge dopant.
Prompted by the pioneering research [6], many other groups have studied the doping effect on the NTE property for some counterparts of Mn3(Cu1-xGex)N, typically like Mn3(Cu1-xSnx)N [7], Mn3(Zn1-xGex)N [8], Mn3(Cu0.6SixGe0.4-x)N [9], and so on [10, 11]. All the efforts reported in these papers showed that it was the doped Ge that brought the NTE property into these anti-perovskite compounds, just like that of the Mn3(Cu1-xGex)N compound. Therefore, the doped Ge was regarded as a “magic element” for the observed NTE property in the class of the anti-perovskite manganese nitride compounds [6].
Recently, revealing the nature of such a “magic element” in the elegant performance of the NTE property has aroused many attractions, Matsuno and coworkers [12] proposed that the Ge dopant caused a local structural distortion in the Mn3(Cu0.7Ge0.3)N compound, leading to the negative thermal expansion of the system in a broad temperature range. This was supported by the neutron powder diffraction and NMR measurements on Mn3(Cu1-xGex)N [13]. On the other hand, through fitting the possible magnetic configurations to their experimental data, Kodama et al. [14] found that the ordering moment of the ground state known as the anti-ferromagnetic Γ5g (seen in Figure 1(a)) of the Mn3(Cu1-xGex)N compound gradually increased with decreasing temperature. This might also be responsible for the NTE property of the concerned compound. Even so, it is still unclear why the doped Ge atoms make the Mn3(Cu1-xGex)N compound be of NTE property.
The magnetic configurations of (a) Γ5g, (b) C0, and (c) C1. The arrows in the figures indicate the orientations of the local magnetic moments at the ions.
Usually, the thermal expansion property of a compound arises from different mechanisms, such as the lattice vibration, the temperature-excited magnetic phase transition, and the diffusion of ions/atoms between interstitial sites. Among these different mechanisms, both the lattice vibration and the phase transition are mainly responsible for the thermal expansion property of many materials. However, for our concerned Mn3(Cu1-xGex)N compound, the relation between its thermal expansion behavior and the lattice vibration is still puzzling.
In this theoretical study, the Mn3(Cu0.5Ge0.5)N compound is selected as a representative, since it exhibited a striking NTE property in experiment. Our calculations indicate that the lattice vibration contributes to the positive thermal expansion, and the magnetic phase transition from anti-ferromagnetic phase to paramagnetic phase contributes to the negative thermal expansion. Both of the PTE and NTE are compensative to each other, resulting in the observed thermal expansion behavior of the Mn3(Cu0.5Ge0.5)N compound. Furthermore, a novel mechanism is proposed to explain the role of the doped Ge for the broad working temperature range of NTE in the Mn3(Cu0.5Ge0.5)N compound.
2. Computational Details
In our calculations, a 2×2×2 supercell consisting of 20 atoms is used to simulate the magnetic structures of Mn3(Cu0.5Ge0.5)N compound, where half of Cu atoms were replaced by Ge atoms. Based on the symmetry of the supercell, there are three kinds of distributions of Ge in it. These configurations are all optimized by using the VASP code [15, 16], in which the projector augmented wave (PAW) pseudopotentials and the Perdew-Burke-Ernzerhof (PBE) parameterizations of the generalized gradient approximation [17] are adopted. The cut-off energy is 400 eV and a gamma-centered k-point mesh of 5×5×3 is used for the geometry optimization. Optimizations are performed until the residual forces acting on the atoms are smaller than 0.01 eV/Å.
By calculations, the energies of the supercell with three kinds of distributions of Ge in the compound are achieved, respectively. We find that the energy differences between the Ge distributions are less than 18 meV per unit cell. Such very small energy differences imply that these different distributions of Ge in Mn3(Cu0.5Ge0.5)N compound may exist at finite temperatures. Even so, we select the lowest-energy one to serve our following calculations.
3. Results and Discussions3.1. Vibrational Effect on the Variation of Volume
Structurally, Mn3(Cu0.5Ge0.5)N is composed of the corner-shared Mn6N octahedron, somewhat like the structural feature in the ZrW2O8 compound [18], where an open-framework structure consisting of corner-linked rigid polyhedra associates with the transverse vibrational modes that cause NTE. Thus, it is necessary to examine the thermal expansion behavior induced by the vibrational modes. Within the quasiharmonic approximation, the thermal expansion coefficient αL is related to the vibrational modes of a material, which is calculated through the formula [19]
(1)αL=kB3BV∑j=1nγj(ℏωjkBT)2exp(ℏωj/kBT)[exp(ℏωj/kBT)-1]2,
where ωj is the jth frequency, B the bulk modulus, V the volume of the system, and γj the Grüneisen parameter relevant to the jth vibrational mode. Here, the Grüneisen parameter is defined as
(2)γj=-(∂lnωj∂lnV)T.
To evaluate the coefficient of the thermal expansion, we first compute the bulk modulus and all vibrational frequencies at Γ point for the supercell. Based on the obtained energy as a function of the volume of the system the bulk modulus of the concerned system is predicted to be about 128.2 GPa (to accurately achieve the vibrational modes, a larger supercell containing 2×2×2 chemical unit cells is selected in our calculations relevant to the vibrational modes.) Meanwhile, the dynamic matrix of the considered system is generated through the forces acting on each atom, with using the frozen phonon approximation. By solving the dynamic matrix, we obtain all of the vibrational frequencies and the related eigenmodes. According to the obtained eigenmodes and the bulk modulus, the thermal expansion coefficient αL is obtained, as shown in Figure 2. Apparently, the value of αL is very small at low temperature and then increases rapidly with increasing temperature. When the temperature reaches 100 K or so, αL almost converges to the value of 11.50×10-6 K−1. This value is in good agreement with that (13×10-6 K−1) from experiment [6]. Such agreement indicates that the widely used frozen phonon approximation is also reliable to evaluate the dynamical property of the compound. More importantly, the calculated values of αL shown in Figure 2 are positive in the whole considered temperature range. This strongly predicts that the transverse motion from the rigid octahedron has no contribution to the observed NTE property in the Mn3(Cu0.5Ge0.5)N compound.
The coefficient of the thermal expansion αL caused by the lattice vibration as a function of the temperature.
3.2. Magnetic Phase Transition versus the Variation of Volume
Experiments found the observed NTE behavior in the Mn3(Cu0.5Ge0.5)N compound accompanied with a magnetic transition from the low temperature Γ5g-type anti-ferromagnetic (AFM) phase to the high temperature paramagnetic (PM) phase [6, 20]. However, such NTE behavior was not found in the Mn3CuN compound yet [6, 21]. Clearly, the observed NTE property of Mn3(Cu0.5Ge0.5)N does tightly couple with the doped Ge atoms. To understand the role of the doped Ge atoms in the concerned NTE behavior, it is valuable to study the Mn3CuN compound firstly.
For the undoped MnCuN compound, where the Mn ions contribute to local magnetic moments, showing non-collinear magnetic configuration. Its ground state was experimentally demonstrated to be ferromagnetic. Furthermore, a magnetic configuration (notated as C0, seen in Figure 1(b)) was suggested to explain the ground state [21]. Our calculations reproduce the magnetic configuration C0, with the ferromagnetic components of local magnetic moments of 0.18 μB and 0.30 μB for the Mn ions locating at two unequal planes, respectively, which are close to the values of 0.20 μB and 0.65 μB reported in experiment [21]. However, such agreement does not directly support C0 configuration as the ground state of the Mn3CuN compound. In fact, a ground state configuration should have the lowest energy among all of the possible magnetic configurations. To examine whether the configuration C0 is the ground state or not, we generate many kinds of magnetic configurations, followed by full relaxations. Among these considered magnetic configurations, a lowest-energy configuration (notated as C1, seen in Figure 1(c)) is found, whose energy is lower by about 0.3 eV per primitive unit cell than that of the C0 configuration. Furthermore, we find that the configuration C1 characterizes the FM feature, and the ferromagnetic component averaged on Mn ions in this configuration is about 0.1 μB, being smaller than the experimental value of 0.35 μB.
Moreover, for this undoped compound, we find that these magnetic configurations have different energies, meaning that these magnetic configurations may be thermally excited at different temperatures. On the other hand, the equilibrium volumes of all considered magnetic configurations are very close to each other. Hence, the thermally excited transformation between different magnetic configurations does not change the volume of the compound nearly. In other words, the Mn3CuN compound does not show NTE on heating. This is consistent with the observation in experiment very well [6].
Then, we turn to the Mn3(Cu0.5Ge0.5)N compound. For this compound, the magnetic state at low temperature was demonstated to be AFM in experiment. Moreover, experiments suggested that a special AFM configuration named as Γ5g just corresponded to the ground state of the compound [13]. On the theoretical side, it is necessary to clarify its magnetic ground state. To do this, we randomly generate about eighty magnetic configurations plus the Γ5g configuration in this compound. After full relaxations of the spins and atomic positions for each concerned case, we only have twenty four inequivalent magnetic configurations. Figure 3(a) displays the energy as a function of the lattice constant for each magnetic configuration. From Figure 3(a), one can find that the lowest-energy magnetic configuration is Γ5g, and its equilibrium lattice constant is evaluated to be 3.864 Å, in agreement with the experimental value of 3.910 Å [13]. Furthermore, our calculations indicate that the local magnetic moments of the Mn ions in the Γ5g configuration are 2.63 μB, being in good consistent with the experimental value of 2.47±0.27μB [21]. Such agreement supports that Γ5g is really the magnetic ground state for Mn3(Cu0.5Ge0.5)N compound.
(a) The energies of the considered magnetic configurations as a function of the lattice constant forMn3(Cu0.5Ge0.5)N. The solid line represents the energy curve of Γ5g and the dashed lines represent the energy curves of the magnetic metastable configurations. (b) The energies as a function of the lattice constant for the distorted Γ5g configurations (dashed lines) and the perfect Γ5g configuration (solid line). The energy at the equilibrium volume of the Γ5g configuration is taken as a reference.
The experiments revealed that with raising temperature, the Mn3(Cu0.5Ge0.5)N compound exhibited a volume contraction. Basically, the temperature effect not only excites more magnetic configurations of the system, but also distorts each magnetic configuration to some extent. To qualitatively evaluate the distorted magnetic configurations, we randomly deviate the orientations of the local moments of Mn ions from their initial orientations for each magnetic configuration, and thus a series of the distorted configurations can be generated through controlling the deviating extent. The volume of the compound at each distorted configuration is optimized with fixing the orientations of the local moments. Our calculations show that the volume of the distorted magnetic configurations almost does not change, except for the Γ5g configuration. As shown in Figure 3(b), with distorting the Γ5g configuration seriously, the energy of the system increases significantly, and the volume of the lattice decreases largely. This predicts that the volume contraction on heating is actually contributed from the distorted Γ5g configuration together with the thermal excited metastable configurations.
On the other hand, when the compound was heated up to the Neel temperature, it transfers from the AFM Γ5g state to the PM state [13]. As we know, the PM state is a magnetic disorder state, which cannot be simply approximated to be either a non-magnetic state or a magnetic ordering state with lowmagnetic moment in a theoretical calculation. Since there are many different magnetic configurations in the compound, the PM state may be approximatively regarded as a mixture of many magnetic structures that can be thermally excited, when the compound is heated up above the Neel temperature. According to this, the volume of the compound at the PM state is expressed as the statistically averaged volume of many different magnetic structures including Γ5g, by using the Boltzmann distribution function. Combining this consideration with the energy curves shown in Figure 3, we statistically compute the lattice constant (a(T)) at each considered temperature. Of these temperature-dependent lattice constants, the lattice constant at 500 K is selected as a reference, and thus we obtain the ratio of the lattice constant contraction, Δa/a=(a(T)-a(500K))/a(500K). As displayed in Figure 4(a), the lattice constant exhibits a gradual decrement with increasing temperature. This implies that the considered magnetic interactions really cause the NTE property in the whole temperature range. To go further, combining this NTE behavior with the PTE behavior from the lattice dynamics shown in Figure 2, we find that the compound exhibits the NTE behavior only in a temperature range of 250–340 K, as seen in Figure 4(b), in consistent with the results reported in experiment [6]. Our calculations above reveal that the volume contraction from the magnetic couplings is actually compensated by the lattice vibrations to some extent. The effect of such compensation is strongly relevant to the temperature. Namely, when the temperature is either below 250 K or over 340 K, the volume change of the compound from the lattice vibration overpasses that from the magnetic interaction, showing the PTE behavior. And in the temperature range between 250 and 340 K, the NTE effect is much stronger than the PTE effect, so that the compound exhibits a net effect of NTE in such a temperature window.
(a) The ratio of the lattice contraction caused by the magnetic interaction. (b) The ratio of the lattice change arising from the magnetic interaction and the lattice vibration. Δa represents the difference between the lattice constant at temperature T and the lattice constant at 500 K.
3.3. The Role of the Doped Ge Atoms in the NTE Property
As mentioned above, it is the doped Ge atoms that cause the compound to be of NTE property. Meanwhile, the observed NTE property is stemmed from the magnetic interaction between ions in the Mn3(Cu0.5Ge0.5)N compound. So, the nature of the observed NTE property is essentially associated with how the doped Ge atoms affect the magnetic couplings between ions. Physically, the magnetic couplings between ions in a compound can be expressed by the exchange parameters between different ions within the framework of the classical isotropic Heisenberg theory [22–25]. In our concerned compound, three typical kinds of exchange parameters, J1, J2, and J2′ are considered. Here, J1 is the coupling between an Mn ion and its nearest neighbor Mn; J2 is the super exchange between an Mn ion and its second neighbor Mn ion, separated by an N atom; J2′ is the coupling between an Mn ion and its second neighbor Mn ion, without any atom between them. Based on the classical isotropic Heisenberg theory [22–25], the energy gain from the magnetic couplings is expressed as 24J1-12J2-24J2′ for Mn3(Cu0.5Ge0.5)N and 8J1-12J2+8J2′ for Mn3CuN.
Table 1 lists the calculated exchange parameters for both compounds. Among the three kinds of exchange parameters, the amplitude of J1 for the Mn3(Cu0.5Ge0.5)N compound is considerably bigger than that for the Mn3CuN compound. This indicates that the doped Ge strongly enhances the magnetic interaction between the nearest-neighboring Mn ions in the Mn3(Cu0.5Ge0.5)N compound. It is worth noting that in our concerned anti-perovskite manganese nitride compound, the couplings between the nearest-neighboring Mn ions form a triangular structure in space. Such a triangular feature is just responsible for the noncollinear anti-ferromagnetic configuration, which is the basic characteristic of the Γ5g configuration. So, the doped Ge atoms effectively enhance the AFM Γ5g configuration much more than the other configurations in the Mn3(Cu0.5Ge0.5)N compound. This is just corresponding to the fact that Mn3(Cu0.5Ge0.5)N has the ground state of Γ5g, whereas Mn3CuN does not have yet [6, 21].
Exchange parameters (meV) for Mn3CuN and Mn3(Cu0.5Ge0.5)N compounds.
Compound
J1
J2
J2′
Mn3CuN
−7.5
22.5
−4.9
Mn3(Cu0.5Ge0.5)N
−19.90
20.3
−8.8
As reported in experiments, the doped Ge atoms not only caused the NTE property, but also made the NTE occurring in a broad temperature range for Mn3(Cu0.5Ge0.5)N. However, these features did not occur in the Mn3(Cu1-xMx)N (M = Zn, Si, and Ga) compounds. Basically, the different behaviors in NTE of the compounds caused by the different dopants tightly couple with the different electronic structures of the dopants. For the isolated atoms of Cu, Ge, Zn, and Mn, the energies of their highest occupied orbitals are in the order of EZn4s(-5.93eV)<EMn4s(-5.47eV)<ECu4s(-4.98eV)<EGe4p(-4.38eV). Clearly, the energy differences of the orbitals between Mn and M (M = Cu, Zn and Ge) are quite big. This means that the 4s orbital of either Cu or Zn and the 4p orbital of Ge do not significantly couple with the valence orbital of Mn in the compound. So, in the compound, most electrons on these atomic orbitals are unpaired and thus locate at the unsaturated states which are almost at the Fermi level. Especially, EGe4p is much higher than the others, and the number of electrons on Ge 4p is larger than the other orbitals considered above. Thus the doped Ge in the compound contribute more states at the Fermi level than the other cases. Meanwhile, the local density of states (LDOS) at the Fermi level from the Mn ions drops down largely when the content of Ge is 50% [26] in the compound. With rising temperature, more electrons from the Ge 4p orbital are thermally excited into the conduction band in such metallic systems. These conduction electrons may polarize some of local electrons in the compound, altering the magnetic couplings between Mn ions and the local magnetic moments accordingly. Namely, some of magnetic states in the Ge-doped compound is easily excited by temperature, even temperature rises slightly. As a result, the Ge-doped compound shows a gradually temperature-dependent change in volume. This is corresponding to the wide temperature range for the NTE working, as reported in experiments [6].
4. Conclusion
We have investigated the lattice vibrations and the magnetic phase transition for the Ge-doped metallic anti-perovskite Mn3(Cu0.5Ge0.5)N compound, based on first-principles calculations. We find that the lattice vibrations contribute to the PTE property rather than the NTE property. Our calculations indicate that the observed NTE is actually resulted from the magnetic phase transition from the Γ5g phase to the PM phase. Furthermore, we suggest that the thermal effect excites more valence electrons from the doped Ge atoms to the conduction bands in the compound. These donated conduction electrons polarize the local electrons and change the local magnetic moments of Mn ions, leading to the contraction of the compound in volume.
Acknowledgments
This work is supported by the National Basic Research Program of China (2009CB939901) and the National Science Foundation of China with Grant no. NSFC10974184. The HP-LHPC of USTC is acknowledged for computational support. The authors gratefully acknowledge the valuable comments from Professor Koshi Takenaka.
SleightA. W.Compounds that contract on heating19983712285428602-s2.0-0001394201BarreraG. D.BrunoJ. A. O.BarronT. H. K.AllanN. L.Negative thermal expansion2005174R217R2522-s2.0-1324429570810.1088/0953-8984/17/4/R03FruchartD.BertautE. F.Magnetic studies of the metallic perovskite-type compounds of manganese19784437817912-s2.0-0001550815KimW. S.ChiE. O.KimJ. C.HurN. H.LeeK. W.ChoiY. N.Cracks induced by magnetic ordering in the antiperovskite ZnNMn32003681742-s2.0-0346880258172402FruchartD.BertautE. F.MadarR.LorthioirG.FruchartR.Structure magnetique et rotation de spin de Mn3NiN1971921179317972-s2.0-0008504878TakenakaK.TakagiH.Giant negative thermal expansion in Ge-doped anti-perovskite manganese nitrides2005872632-s2.0-2974445414310.1063/1.2147726261902TakenakaK.AsanoK.MisawaM.TakagiH.Negative thermal expansion in Ge-free antiperovskite manganese nitrides: tin-doping effect200892132-s2.0-3804904073510.1063/1.2831715011927SunY.WangC.WenY.ZhuK.ZhaoJ.Lattice contraction and magnetic and electronic transport properties of Mn3Zn1-xGexN2007912332-s2.0-3684908849310.1063/1.2822813231913HuangR.LiL.CaiF.XuX.QianL.Low-temperature negative thermal expansion of the antiperovskite manganese nitride Mn3CuN codoped with Ge and Si200893832-s2.0-5134910524110.1063/1.2970998081902WangB. S.TongP.SunY. P.LiL. J.TangW.LuW. J.ZhuX. B.YangZ. R.SongW. H.Enhanced giant magnetoresistance in Ni-doped antipervoskite compounds GaCMn3-xNix (x = 0.05, 0.10)20099522322250910.1063/1.3268786HuangR. J.XuW.XuX. D.LiL. F.PanX. Q.EvansD.Negative thermal expansion and electrical properties of Mn3(Cu0.6NbxGe0.4 - x)N (x=0.05−0.25) compounds20086216238123842-s2.0-4194913510210.1016/j.matlet.2007.12.003MatsunoJ.TakenakaK.TakagiH.MatsumuraD.NishihataY.MizukiJ.Local structure anomaly around Ge dopantsin Mn3Cu0.7Ge0.3N with negative thermal expansion2009941832-s2.0-6554915316410.1063/1.3129169181904IikuboS.KodamaK.TakenakaK.TakagiH.TakigawaM.ShamotoS.Local lattice distortion in the giant negative thermal expansion material Mn3Cu1-xGexN20081012042-s2.0-5684911009510.1103/PhysRevLett.101.205901205901KodamaK.IikuboS.TakenakaK.TakigawaM.TakagiH.ShamotoS.Gradual development of Γ5g antiferromagnetic moment in the giant negative thermal expansion material Mn3Cu1-xGexN (x~0.5)2010812282-s2.0-7795632791610.1103/PhysRevB.81.224419224419KresseG.HafnerJ.Ab initio molecular dynamics for liquid metals19934715585612-s2.0-1284428624110.1103/PhysRevB.47.558KresseG.FurthmüllerJ.Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set1996541611169111862-s2.0-2442537377PerdewJ. P.BurkeK.ErnzerhofM.Generalized gradient approximation made simple19967718386538682-s2.0-4243943295ErnstG.BroholmC.KowachG. R.RamirezA. P.Phonon density of states and negative thermal expansion in ZrW2O8199839667071471492-s2.0-003251205410.1038/24115WangK.ReeberR. R.Mode Grüneisen parameters and negative thermal expansion of cubic ZrW2O8 and ZrMo2O820007616220322042-s2.0-0001730066TakenakaK.TakagiH.Magnetovolume effect and negative thermalexpansion in Mn3(Cu1-xGex)N20064734714742-s2.0-3364677644810.2320/matertrans.47.471IikuboS.KodamaK.TakenakaK.TakagiH.ShamotoS.Magnetovolume effect in Mn3Cu1-xGexN related to the magnetic structure: neutron powder diffraction measurements200877242-s2.0-3854917298310.1103/PhysRevB.77.020409020409RosengaardN. M.JohanssonB.Finite-temperature study of itinerant ferromagnetism in Fe, Co, and Ni1997552214975149862-s2.0-0000480789ZhouY. M.WangD. S.KawazoeY.Effective ab initio exchange integrals and magnetic phase transition in fcc Fe and Mn antiferromagnets19995913838783902-s2.0-0000156599WangJ. T.ZhouL.WangD. S.KawazoeY.Exchange interaction and magnetic phase transition in layered Fe/Au(001) superlattices2000625335433602-s2.0-0001114697GaoX.ZhouY. M.WangD. S.Analysis of the determination of magnetic phase transition temperatures by the first-principles calculations2002251129372-s2.0-003680123710.1016/S0304-8853(02)00352-9QuB. Y.HeH. Y.PanB. C.The elastic properties of Mn3(Cu1-xGex)N compounds201114804212510.1063/1.3658025