化合物TbMn6Sn6的磁性和磁相变
郭光华, 张海贝
(中南大学 物理科学与技术学院, 长沙 410083)
摘 要: 采用分子场理论模型对稀土-过渡族金属间化合物的磁性、 自旋重取向相变和场诱导的一级磁相变进行了研究。 从理论上计算了TbMn6Sn6的易磁化方向随温度的变化以及不同温度时的磁化曲线, 描述了自旋重取向相变和场诱导的一级磁相变。 基于单离子模型计算了Tb离子和化合物的一阶和二阶磁晶各向异性常数及其随温度的变化, 得到了和实验数据基本一致的结果。 结果表明: Mn次晶格和Tb次晶格磁晶各向异性能之间的相互竞争是导致自旋重取向相变的物理机制, 而稀土离子的四阶晶场项对场诱导的一级磁相变产生很大的影响。
关键词: 稀土-过渡族金属间化合物; 自旋重取向; 磁相变; 磁晶各向异性 中图分类号: O482.5
文献标识码: A
Magnetic properties and magnetic phase transition of compound TbMn6Sn6
GUO Guang-hua, ZHANG Hai-bei
(School of Physics Science and Technology,Central South University, Changsha 410083, China)
Abstract: The molecular field theory of exchange interaction was applied to investigate the magnetic properties, spin-reorientation transition and field-induced magnetic phase transitions of intermetallic compound TbMn6Sn6. The temperature dependence of easy magnetization direction and magnetic curves at different temperatures were calculated. The theoretical results have good agreement with the experimental data. In the framework of single ion model the magnetocrystalline anistropic constants and their temperature dependence of compound were calculated. The results show that the competition between the magnetocrtstalline anistropic energy of Tb and Mn sunlattices gives rise to the spin-reorientation transition. The fourth-order crystal field parameter (B04) of Tb ion plays an important role on the field-induced magnetic phase transition.
Key words: rare earth-transition intermetallic compound; spin-reorientation transition; magnetic phase transition; magnetocrystalline anisotropy
近年来, 由于三元系稀土-过渡族金属间化合物RM6X6(其中R为稀土元素, M为Fe, Mn和Cr, X为Sn和Ge)丰富的磁学性质引起了人们的广泛关注[1-15]。 重稀土RMn6Sn6化合物具有HfFe6Ge6型的层状晶体结构, 空间群为P6/mmm, 其晶体结构可描述为由R和Mn原子层沿(001)面按Mn-R-Mn-Mn-R-Mn序列堆栈而成[5], 其中R原子层具有六角平面结构(H), Mn原子层为笼目网格结构(K)。 从磁学的角度看, RMn6Sn6化合物由两种不同的磁次晶格构成, 即R次晶格和Mn次晶格。 中子衍射和磁测量的研究结果表明[5-7], 对于重稀土Tb、 Dy和Ho化合物, 由于R-Mn间具有较强的交换相互作用, 当温度低于居里温度Tc时(Tb、 Dy和Ho化合物的Tc分别为423、 393和376K), Mn次晶和R次晶格同时进入铁磁有序状态, 并且Mn次晶格和R次晶格构成反铁磁耦合。 当温度进一步降低到Tt时, 化合物中发生自旋重取向相变, 其易磁化方向从c面向c轴偏转, Tb, Dy和Ho化合物的自旋重取向温度Tt分别为330、 320和200K(文献[7]测量TbMn6Sn6单晶样品得到的自旋重取向相变温度为310K)。 对DyMn6Sn6和HoMn6Sn6, 其自旋重取向并不充分, 温度约为220和100K时, 易磁化方向与c轴的夹角分别为45°和48°, 并且一直保持到2K; 而对TbMn6Sn6, 当温度低于自旋重取向温度时, 其易磁化方向完全转到c轴。 此外, 在外磁场的作用下, TbMn6Sn6中还发生场诱导的一级磁相变[7-9]。 本文作者在分子场理论模型基础上定量地研究了化合物TbMn6Sn6的磁性以及自旋重取向相变和场诱导的磁相变, 并基于磁晶各向异性的单离子模型计算了该化合物的磁晶各向异性常数及随温度的变化关系。
1 理论模型
TbMn6Sn6由Tb和Mn两种次晶格构成, 根据分子场理论模型, Tb和Mn离子的哈密顿量可表示为
式中 J和S分别为Tb离子和Mn离子的总角动量和自旋角动量; gJ和gS为Tb和Mn离子的朗德因子; μB为玻尔磁子; B02和B04为Tb3+的晶场常数; O02和O04为Stevens算符; K1t为Mn离子的磁晶各向异性能常数; θt为Mn离子磁矩的热平均值与晶轴c的夹角; H为外磁场; Hex, Tb和Hex, Mn为作用在稀土离子和锰离子上的分子场, 其形式为
Hex, Tb=λTb-Tb〈μTb〉+6λTb-Mn〈μMn〉(3)
Hex, Mn=λMn-Mn〈μMn〉+λTb-Mn〈μTb〉(4)
式中 λTb-Tb, λTb-Mn, λMn-Mn分别为Tb-Tb, Tb-Mn和Mn-Mn间交换相互作用分子场系数, 〈μTb〉和〈μMn〉分别代表Tb和Mn离子磁矩的热平均值。
温度为T时, 系统自由能可表示为
式中 kB玻尔兹曼常数; ZTb, ZMn为Tb和Mn离子的配分函数, 可分别表示为
由热平衡条件可以得到〈μMn〉和〈μTb〉的热平均值及它们与c轴间的夹角, 由此确定化合物的易磁化方向、 磁结构和磁化强度随温度和外磁场的变化。
2 结果与讨论
在具体计算中, 取4.2K时Tb和Mn离子的磁矩分别为9μB和2.17μB[7]。 通过拟合TbMn6Sn6化合物的各种实验数据[6-8], 其中包括自旋重取向相变温度、 热磁曲线以及不同温度下磁场沿不同晶轴方向时的磁化曲线, 所确定的其它磁参数为: λMn-Mn=180.1T/μB, λTb-Mn=-7.89T/μB, λTb-Tb=2.19T/μB, B02=-0.5K, B04=-0.008K, K1t0=-1.8K(4.2K时Mn离子的磁各向异性常数)。
图1给出了理论计算的TbMn6Sn6化合物的易磁化方向与c轴夹角随温度的变化曲线。 从图1可知, 高温时, TbMn6Sn6的易磁化方向在c面上, Tb和Mn离子磁矩与c轴夹角为90°, T=310K时发生自旋重取向相变, 易磁化方向偏转到c轴, 并一直保持到最低温度。 在整个温度范围内Tb和Mn离子磁矩始终保持反平行排列。 理论计算结果与中子衍射的实验结果基本一致[6]。
导致稀土-过渡族金属间化合物自旋重取向相变的因素有多种, 两种次晶格磁晶各向异性能之间
图1 化合物的易磁化方向与c轴的夹角随温度的变化
Fig.1 Temperature dependence of easy magnetization direction of compound
的相互竞争以及不同阶次的磁晶各向异性常数间的竞争都有可能导致易磁化方向的改变。 从拟合实验数据所得的磁参数来看, Mn次晶格的磁各向异性常数为负值, 表明其易磁化方向在c面上。 中子衍射的研究结果也证明TbMn6Sn6化合物中Mn次晶格的易磁化方向始终保持在c面[6]。 作用在Tb离子上的晶场参数B02和B04均为负值(相应的A02和A04为正), 这都有利于Tb次晶格的磁矩沿着c轴排列。 因此, Mn次晶格和Tb次晶格磁晶各向异性能之间的相互竞争是导致自旋重取向相变的物理机制, 高温时Mn次晶格磁晶各向异性能起主要作用, 低温时Tb次晶格的磁晶各向异性能占主导地位。
对TbMn6Sn6单晶样品和取向样品的实验研究表明[7, 8], 当温度低于自旋重取向温度TSR且在TSR附近时, 垂直c轴的磁化会出现场诱导的一级磁相变。 图2所示为外磁场垂直c轴方向时不同温度下理论计算的TbMn6Sn6化合物的磁化曲线和相应的实验数据。 很显然, 温度为320K时, 易磁化方向在c面上, 磁化曲线为一直线, 当温度低于自旋重取向温度TSR, 在外场作用下发生了场诱导的一级磁相变。 首先, 随着外场的增大, Mn磁矩逐渐向外场方向偏转, 当外场达到某一临界场时, 磁矩突然转到外场方向, 相变前和相变后, Mn磁矩和Tb磁矩始终保持反平行排列, 磁结构均为亚铁磁结构, 而不是如文献[9]中所述, 这一磁相变伴随着自旋的翻转。 随着温度的降低相变临界磁场逐渐增大, 温度低于160K, 在实验允许的磁场范围
图2 外磁场垂直于c轴时化合物的磁化曲线
Fig.2 Magnetization curves of compound TbMn6Sn6 under magnetic field perpendicular to c-axis at different temperatures
内已观察不到磁相变, 但本文的计算表明, 在更高的磁场也会出现场诱导的一级磁相变。 理论计算表明稀土离子的四阶晶场项对这一磁相变有很大的影响。
基于理论计算所得到的晶场参数, 在磁晶各向异性能的单离子模型基础上, 计算了化合物的磁各向异性常数及随温度的变化。 对Mn次晶格, 可以认为其磁晶各向异性常数满足Akulov-Zener定律, 即Kt(T)∝M(T)3。 对稀土次晶格, 磁晶各向异性系数Kmn与晶场参数Bmn满足如下关系:
式中 J2=J(J-1/2), J4=J2(J-1)(J-3/2), J为总角动量。 根据磁晶各向异性能的单离子模型, 磁晶各向异性系数随温度的变化具有如下关系:
式中 In+1/2为双曲贝塞尔函数; L为朗之万函数。
对RMn6Sn6化合物稀土离子磁晶各向异性常数与各向异性系数之间的关系为
根据上述公式和所得的晶场参数, 可以计算出Tb离子磁晶各向异性常数随温度的变化。 图2(a)给出了理论计算的Tb离子磁晶各向异性常数K1R和K2R随温度的变化曲线, 可以看出, 在整个磁有序温度范围内K1R始终大于零, 说明K1R利于Tb离子磁矩沿c轴方向排列; 而对于K2R, 其值均为负。 TbMn6Sn6总的磁各向异性常数为K1=K1R+6K1t, K2=K2R, 其随温度的变化曲线如图2(b)所示。 可以看出, 理论计算的K1和K2值基本与实验数据相符合。 在整个磁有序温度范围内K1和K2不满足条件0〈-K1〈2K2, 因此TbMn6Sn6不会出现如HoMn6Sn6和DyMn6Sn6的锥面易磁化方向。
3 结论
在分子场理论模型基础上, 定量研究了稀土-金属间化合物TbMn6Sn6的磁性、 自旋重取向相变和场诱导的一级磁相变。 计算了化合物的易磁化方向随温度的变化, 研究表明, Tb次晶格和Mn次晶格磁晶各向异性能之间的相互竞争是导致自旋重取向相变得物理机制。 计算了外磁场垂直于c轴时不同温度下的磁化曲线, 较好地描述了场诱导的一级磁
图3 TbMn6Sn6中Tb次晶格的磁晶各向异性常数(a)和总的磁晶各向异性常数随温度的变化(b)
Fig.3 Temperature dependence of magnetic anistropic constants of Tb sublattice(a) and magnetic anistropic constants of compound TbMn6Sn6 (b)
相变。 计算结果表明, 相变前和相变后, 化合物始终保持亚铁磁结构, 相变伴随着磁化强度向外场方向的突然偏转, 而不是自旋翻转。 基于磁晶各向异性的单离子模型计算了Tb离子和化合物的磁晶各向异性常数及其随温度的变化。 通过对化合物磁晶各向异性常数的分析表明, TbMn6Sn6不会出现如HoMn6Sn6和DyMn6Sn6的锥面易磁化方向。
REFERENCES
[1]Zajkov N K, Mushnikov N V, Gerasimov E G, et al. Magnetic properties of Tb1-xYxMn6Sn6 compounds[J]. J Alloys Comp, 2004, 363: 40-45.
[2]Horia T, Shiraishia H, Indoh K, et al. Magnetic properties of Y1-xDyxMn6Sn6 alloys[J]. J Magn Magn Mater, 2004, 272-276: e433-e434.
[3]WANG Ru-wu, ZHANG Li-gang, ZHANG Shao-ying, et al. Magnetism and magnetoresistance of Dy1-xSmxMn6Ge6 (x=0.2-1.0) compounds[J]. J Phys D: Appl Phys, 2004, 37: 385-388.
[4]Lefèver C, Venturini G A. Magnetic study of HoMn6Sn6-xGax single crystals (0.14≤x≤1.89)[J]. J Magn Magn Mater, 2004, 268: 374-379.
[5]Venturini G, Chafik E l B, Malaman B. Magnetic properties of RMn6Sn6 (R=Sc, Y, Gd-Tm, Lu) compounds with HfFe6Ge6 type structure[J]. J Magn Magn Mater, 1991, 94: 35-42.
[6]Malaman B, Venturini G, Welter R, et al. Magnetic properties of RMn6Sn6 (R=Gd, Er) compounds from neutron diffraction and Mossbauer measurements[J]. J Magn Magn Mater, 1999, 202: 519-534.
[7]Clatterbuck D M, Gschneidner K A Jr. Magnetic properties of RMn6Sn6 (R=Tb, Ho, Er, Tm, Lu) single crystals[J]. J Magn Magn Mater, 1999, 207: 78-94.
[8]Zajkov N K, Mushnikov N V, Bartashevich M I, et al. Magnetization processes in the TbMn6Sn6 compound[J]. J Alloys Comp, 2000, 309: 26-30.
[9]Hu J, Wang K Y, Hu B P, et al. Magnetic transition and coercivity in TbMn6Sn6[J]. J Phys Condens Matter, 1995, 7(5): 889-893.
[10]Schobinger-Papamantellos P, Andre G, Rodriguez-Carvajal J, et al. Magnetic ordering in RMn6-xCrxSn6(x=1, 2) compounds studied by neutron diffraction and magnetic measurements[J]. J Alloys Comp, 2000, 306: 47-55.
[11]Schobinger-Papamantellos P, Andre G, Rodriguez-Carvajal J, et al. Magnetic ordering in DyMn6-xCrxSn6 (x=1, 2) compounds studied by neutron diffraction and magnetic measurements[J]. J Magn Magn Mater, 2000, 219: 22-32.
[12]Li Y B, Zhang S Y, Yao J L, et al. Magnetic properties and magnetoresistance effect of Y1-xGdxMn6Sn6 (x=0-1) compounds[J]. Chin Phys, 2003, 12(3): 328-331.
[13]Zhang Li G, Wang R W, Li L, et al. Magnetic properties and magnetoresistance of HfFe6Ge6-type Dy1-xGdxMn6Ge6 (x=0.1-0.6) compounds[J]. Chin Phys, 2003, 12(5): 562-565.
[14]Canepa F, Napoletano M, Lefevre C, et al. Magnetisation and magnetoresistance of YMn6Sn6-xInx compounds (x=0.10, 0.18, 0.21, 0.38)[J]. J Alloys Comp, 2003, 349: 6-11.
[15]Cadogan J M, Ryan D H, Moze O, et al. Magnetic ordering in ErFe6Sn6[J]. J Phys: Condens Matter, 2003, 15: 1757-1771.
(编辑龙怀中)
基金项目: 湖南省自然科学基金资助项目(03JJY4044); 教育部留学回国人员启动基金资助项目
收稿日期: 2005-04-07; 修订日期: 2005-06-01
作者简介: 郭光华(1964-), 男, 教授, 博士
通讯作者: 郭光华, 教授; 电话: 0731-8836443; E-mail: guogh@mail.csu.edu.cn