-
The lanthanide and actinide metals and alloys are of great interest in experimental and theoretical high-pressure research, because of the unique behavior of the f electrons under pressure and their delocalization and participation in bonding. Cerium (Ce) metal is the first lanthanide element with a 4f electron. It has a very complex phase diagram and displays intriguing physical and chemical properties. In addition, it is expected to be an excellent surrogate candidate for plutonium (Pu), one of the radioactive transuranic actinides with a 5f electron. The bulk properties and phase transformation characteristics of Ce-based alloys are similar to those of Pu and its compounds. Thus, the investigations of Ce-based alloys are necessary and can potentially advance the understanding of the behavior of Pu. In this work, the equation of state, phase transition, elastic and thermodynamic properties of Ce0.8La0.1Th0.1 alloy at high pressure are investigated by using first-principles calculations based on the density-functional theory. The structural properties of the Ce0.8La0.1Th0.1 alloy are in good agreement with the available experimental and theoretical data. The lattice constant a decreases with pressure increasing, while c shows an opposite variation. It is found that the lattice parameter c shows abnormal jump. And the critical volume is located at 20.1 Å3. The axial ratio jumps from a value of about
√2 (corresponding to the fcc structure) to a higher value, which indicates that the fcc-bct transition occurs. And the corresponding transition pressure is located at ~31.6 GPa. When the pressure rises to 34.9 GPa, the bct structure displays a saturated c/a axial ratio close to about 1.67. The Young's modulus E, shear modulus G and the Debye temperature of the fcc phase tend to be " softened” around the phase transition pressure. The vibrational free energy is obtained by using the quasi-harmonic Debye model. And then the thermodynamic properties including the thermal equation of state, heat capacity and entropy under high pressure and high temperature are also predicted successfully. The results show that the heat capacity and entropy increase rapidly with temperature increasing, and decrease with pressure increasing. The high pressure can suppress part of the anharmonicity caused by temperature.-
Keywords:
- Ce-based alloy /
- density functional theory /
- phase transformation /
- high pressure
1. 引 言
稀土元素指镧(La)、铈(Ce)、镨(Pr)等17种元素, 能实现光、电、磁范畴内的特殊能量转换、传导输送及储存等功能. 随着科技的不断发展, 稀土资源已广泛应用于电子信息、污染治理、清洁能源、航空航天、生命技术、医疗卫生等各个方面, 与稀土元素相关的研究成为了现代世界前沿科学的重要组成部分. 储量丰富的稀土元素Ce目前在钢铁、有色金属及合金和发火合金等产业中已有广泛应用, 在电子工业、原子能工业、精密仪器工业的应用也有了新的突破, 为军事、科技、生产助力. Ce作为多形相变材料的一种, 电子态极容易出现跃迁, 具有丰富的相变信息. 常温常压下的Ce存在稳定的γ相, 该相具有fcc面心立方晶体结构. Bridgman[1-3]在测量Ce电导率和等温压缩线时首先发现了Ce在0.8 GPa左右的γ—α相变, 这一现象引发了大量学者的关注. 1949年, Lanson和Tang[4]在X射线衍射实验中证实了Ce的γ—α相变其实是一种同构相变, 相变过程伴随着17%左右的体积塌缩. 2012年, 潘昊等[5]低压冲击金属Ce的γ—α相变, 进行了数值模拟研究. 我们前期运用第一性原理计算结合晶格动力学的方法证明, 随着压强的升高, Ce的相变顺序为α-Ce→α''-Ce→bct-Ce, 相变压强分别为5.36 GPa和14.37 GPa[6].
Ce基合金由于其良好的固溶性, 是核武器关键材料钚(Pu)潜在的替代品, 由于目前的国际规则, 并不能对Pu进行直接的实验研究, 因此Ce基合金作为其潜在的替代品成为了研究的热点. 目前已有学者对Ce-Th, Ce-La, Ce-La-Th合金的物理性质进行了一系列的理论预测与实验测量, 主要关注同构相变的信息. Lawson等[7]用实验测得了Ce0.9Th0.1存在γ—α同构相变. Lawrence等[8]测量了温度在4.2—300 K区间、压强在0—12 kPa区间时不同组分Ce0.9-xLaxTh0.1体系(x = 0.10, 0.11, 0.14, 0.17)的电阻率变化情况, 以此确定了Ce-La-Th合金的γ—α同构相变的压强温度边界. Drymiotis等[9]实验测量了Ce0.8La0.1Th0.1中γ—α相变中伴随的外加磁场变化, Ruff等[10]通过X射线衍射实验证实了Ce0.8La0.1Th0.1也存在体积塌缩的同构相变. 目前对于Ce基合金高压相变的研究较为匮乏. 我们前期用第一性原理计算了CexTh1–x合金体系的结构相变随组分的变化规律[11], 发现fcc-bct相变压强随Ce组分的增加而减小, 也即随着Th的加入, 提高了合金的相变压强. 我们采用密度泛函理论结合准谐近似对两种不同组分的Ce-La合金的相变及状态方程(EOS)进行了研究[12], 发现在压力的作用下, Ce-La合金从fcc结构转变到bct结构, 而在fcc结构中, 又有γ—α同构相变发生. 本文分别计算了两种不同组分的Ce-La合金(Ce0.875La0.125和Ce0.963La0.037)的相变压强, 结果表明随着La含量的增加, 相变压强升高. 高压环境下对Ce-La-Th体系的科学研究还较为匮乏, 本文对Ce0.8La0.1Th0.1在压力下fcc-bct的相变情况和热力学性质进行了研究讨论.
2. 计算细节
采用基于密度泛函理论的虚晶近似(VCA)对无序合金Ce0.8La0.1Th0.1进行计算模拟. VCA是通过将元素的赝势进行混合来产生一种新的势函数. 我们前期也采用这种方法成功预测了Ce-Th合金随组分变化的相图[11]. Ce存在f电子, 在计算中, 需要慎重处理. 我们前期采用密度泛函理论(DFT)+U的方法对纯Ce做了较为深入的探讨[12]. 发现f电子在较低的压强区间, 对Ce的同构相变起决定性的作用. 但是在较高压强区间(同构相变之后), 对fcc-bct的相变影响并不明显. 在前期计算纯Ce和Ce-La合金时, 高压对f电子的关联效应有较强的抑制作用, DFT+U方法计算的体系能量与直接的DFT方法没有明显的差异. 因此, 本文采用DFT直接计算较高压强下的fcc和bct结构的物理性质. 在电子结构计算时, 使用了CASTEP程序包的平面波赝势方法, 交换关联相互函数采用的是广义梯度近似(GGA)框架下的Perdew-Burke-Ernzerhof (PBE)形式, 为获得准确的计算结果, 我们的平面波截断能均取为500 eV, Ce, La和Th的价电子分别取5s25p64f15d16s2, 5s25p65d16s2和6s26p66d27s2. 计算中对fcc和bct两种结构分别采用了15 × 15 × 15 和13 × 13 × 17 Monkhorst pack网格, 为了保证计算的精度使得总能的收敛精度达到10–6 eV/atm. 本文利用准谐德拜模型计算热力学性质, 非平衡的Gibbs自由能(G*)包含了静态能量、晶格振动能以及体积变化带来的能量改变, 可以写成如下形式:
G∗(x;P,T)=E(x)+PV(x)+AVib(x;T), (1) 其中E(x)为静态能量, 可以直接从电子结构计算中得到; PV项代表压强作用下焓值的改变; AVib是Helmholtz自由能, 包含了晶格振动对内能的贡献和熵的改变. Blanco等[13]发展了考虑热效应的准谐Debye模型. 根据这个模型, AVib可以用如下近似求得[14-17]:
AVib(Θ;T)=nkT[98ΘT+3ln(1−e−Θ/T)−D(Θ/T)], (2) 其中Θ是Debye温度, n是每个单位分子式中原子的数目, D(Θ/T)是德拜函数, 定义如下:
D(y)=3y3y∫0x3ex−1dx. (3) 当获得了高温高压的自由能之后, 通过适当的热力学表达式就能求其他热力学性质. 采用这种方法, 邓世杰等[18]对Ti2AlX (X = C, N)的热力学性质做了成功的预测.
3. 结果与分析
3.1 结构相变
我们计算了Ce0.8La0.1Th0.1合金fcc相及bct相的能量-体积关系, 获得了其各自的状态方程. 得到的fcc相及bct相在0 GPa, 0 K下的平衡体积V0, 体积模量B0的数据如表1所列, 同时表中也列出了Ce0.875La0.125及纯Ce的理论值和实验值作为参照. Ce0.8La0.1Th0.1合金bct相的平衡体积V0 = 28.91 Å3, 体积模量B0 = 35.96 GPa. 随着La, Th的加入, 相对于纯Ce而言, 平衡体积V0增加, 而体模量B0减小. 随着Th的加入和La成分的减少, Ce0.8La0.1Th0.1合金的体模量略大于Ce0.875La0.125合金.
表 1 零温零压下fcc相Ce-La-Th合金的平衡体积(V0)及体积模量(B0)Table 1. Equilibrium volume (V0) and bulk modulus (B0) of Ce-La-Th of fcc phase at 0 GPa and 0 K.对Ce0.8La0.1Th0.1合金的fcc—bct相变进行计算和分析, 将得到的fcc相和bct相的体积-压强关系, 与纯Th[19]、纯Ce[20]、Ce0.76Th0.24[21]及Ce0.875La0.125[12]的数据进行比对, 如图1. 本文的数据与已有的Ce基合金及纯Ce的体积-压强变化的规律相吻合.
图 1 体积随压强变化的规律(黑色实点为直接加压结构优化后的结果, 黑色实线为状态方程拟合结果), 并与已有的Ce[20], Th[19], Ce0.875La0.125[12]的计算值及Ce0.76Th0.24[21]实验值进行比较Fig. 1. The EOS of fcc and bct Ce-La-Th together with the experimental data (the black solid point is the result of the structure optimization, the black solid line is the fitting result of the EOS), together with the experimental data for Ce0.76Th0.24[21] and the calculated results for Ce[20], Th[19], Ce0.875La0.125[12].通过对总能的计算结果表明, fcc和bct两种相结构之间焓的差异非常小, 几乎接近计算误差值, 这意味着不能再通过比较焓变的差异来判断相变压强.
对bct相固定体积进行结构优化, 晶格常数的计算结果如图2所示. 随着体积的减小, 当体积减小至20.14 Å3时, 晶格常数a, c的线性规律发生变化, 这意味着晶体结构出现了变化, 相变开始发生. 计算得到的轴向比c/a与已有的Ce0.76Th0.24[21], Ce[6], Th[11], Ce0.875La0.125[12]数据进行了对比. 当压强小于31.6 GPa时, c/a的值约等于, 随着压强的增加, fcc相变得不稳定, 开始向bct相转变. 在34.9 GPa附近时, bct相趋于稳定 c/a的值约等于1.65, 与我们前期计算得到的CeTh, CeLa合金的fcc相, bct相的c/a的值在误差范围内吻合[11, 12].
3.2 弹性性质
本文计算了Ce0.8La0.1Th0.1合金fcc相及bct相结构的弹性性质, 在不同压强下的弹性常数如图3所示. 对于fcc相而言, 随压强增加, 所有的弹性常数都是线性增加的. 弹性常数C11, C44在相变压强附近开始出现非线性的变化. 对bct相而言, 在压强增加的变化趋势下, 越过相变压强后弹性常数线性增加的趋势变得更加明显. 根据弹性常数, 能够计算出多晶的弹性模量, 对于Ce0.8La0.1Th0.1合金, 多晶的弹性模量(剪切模量G、体模量B和杨氏模量E )能够根据Voigt-Reuss-Hill近似得出.
如图4所示, 随着压强的增加, 弹性模量呈线性增加的趋势, 而杨氏模量E和剪切模量G则在相变压强附近出现“变软”的趋势; 随着相变的完成, 新相趋于稳定, 杨氏模量E和剪切模量G随压强增加的关系再次趋于线性. 通过弹性模量,能够计算德拜温度. 如图5所示, 在0 GPa压强下, Ce0.8La0.1Th0.1合金的德拜温度为228.85 K, 高于纯Ce已有的德拜温度研究值[26, 27]. 根据已有的研究[12], 德拜温度随着La组分的增加而增加, Ce0.875La0.125在零温零压下德拜温度为140.9 K, 遂认为在Ce0.8La0.1Th0.1合金中, La和Th元素的掺杂导致了德拜温度的升高, 原子间作用力也会因此升高, 本文的计算结果是合理的. 在大约34.4 GPa时, fcc相德拜温度随压强增加呈减小的趋势, 而bct结构与fcc结构的德拜温度十分接近, 并随压强增加呈线性增加的趋势. 意味着结构相变开始发生, 在41.6 GPa附近, fcc结构转变为bct结构, 此时德拜温度为330 K.
3.3 热力学性质
利用准谐德拜模型获得了Ce0.8La0.1Th0.1的热力学性质和不同温度下的等温压缩曲线. 如图6所示, 当温度为300 K, 压强从0 GPa上升到40 GPa时, fcc相体积缩小了约34.6%. 而在零压下, 当温度从300 K上升到1000 K时, fcc的体积膨胀了约11.1%. 高压状态下温度对体积的影响逐渐减小, 高温的非谐效应在压力的作用下被抑制. 当压强大于40 GPa时bct相稳定存在, 当温度为300 K, 压强从40 GPa上升到80 GPa时, bct 相体积缩小了约13.9%, 而在零压下, 当温度从300 K上升到1000 K时, bct相体积膨胀了约8.2%.
分别计算了fcc相和bct相定容热容CV随温度和压强变化的关系, 以及熵S随温度和压强的变化关系. 如图7所示, CV随着温度升高迅速增加, 在高温下接近25 J/(mol·K)的极限, 热容在不同温度下随压强的变化情况几乎是单调递减的. 当温度超过600 K后, 不同等温线之间的差距变小. 计算得到在常温常压下的熵值约为49.44 J/(mol·K), 随温度的升高, 熵值几乎迅速单增, 在不同的温度下, 熵值随压强的增加而减小.
4. 结 论
本文采用第一性原理计算对Ce-La-Th合金在高压下fcc—bct的相变及热力学性质进行了研究讨论. 获得了Ce-La-Th合金体积随压强变化的规律, 随着压强的增加, 体积减小至20.1 Å3时, 晶格常数a, c的线性变化规律发生突变, 这意味着晶体结构出现了改变, fcc相开始变得不稳定. 通过分析轴向比c/a的值随压强变化的关系发现, 压强在大约31.6 GPa时fcc相开始向bct相转变, 在大约34.9 GPa时, bct相趋于稳定. 弹性性质的计算结果表明, 随着压强的增加, 在相变压强附近, 某些弹性常量出现非线性增加的情况, 当新相结构变得稳定再次趋于线性. 根据弹性常量计算得出的体模量随压强呈线性增加的趋势, 而杨氏模量E和剪切模量G则在相变压强附近出现“变软”的趋势, 当新相bct相趋于稳定时弹性模量再次趋于线性. 弹性德拜温度的计算结果与结构变化的规律类似, 从另一角度说明了相变的信息. 热力学性质的研究结果表明定容热容CV随着温度升高而迅速增加, 在不同温度下随压强的变化呈单调递减. 随温度升高, 熵值迅速单增, 而在不同的温度下, 熵值随压强的增加而减小. 高压对温度带来的非谐效应有明显的抑制作用.
[1] Bridgman P W 1927 Proc. Am. Acad.Arts Sci. 62 207
Google Scholar
[2] Bridgman P W 1951 Proc. Am. Acad. Arts Sci. 79 149
Google Scholar
[3] Bridgman P W 1954 Proc. Am. Acad. Arts Sci. 83 1
[4] Lanson A W, Tang T Y 1949 Phys. Rev. 76 301
Google Scholar
[5] 潘昊, 胡晓棉, 吴子辉, 戴诚达, 吴强 2012 物理学报 61 206401
Google Scholar
Pan H, Hu X M, Wu Z H, Dai C D, Wu Q 2012 Acta Phys. Sin. 61 206401
Google Scholar
[6] Hu C E, Zeng Z Y, Zhang L, Chen X T, Cai L C 2011 Physica B 406 669
Google Scholar
[7] Lawson A C, Williams A, Wire M S 1988 J. Less-common Met. 142 177
Google Scholar
[8] Lawrence J M, Thompson J D, Fisk Z, Smith J L, Batlogg B 1984 Phys. Rev. B 29 4017
Google Scholar
[9] Drymiotis F, Singleton J, Harrison N, Lashley J C, Bangura A, Mielke C H, Balicas L, Fisk Z, Migliori A, Smith J L 2005 J. Phys.: Condens. Matter 17 L77
Google Scholar
[10] Ruff J P C, Islam Z, Das R K, Nojiri H, Cooley J C, Mielke C H 2012 Phys. Rev. B 85 024104
Google Scholar
[11] Hu C E, Zeng Z Y, Zhang L, Chen X R, Cai L C 2010 Solid State Commun. 150 2362
Google Scholar
[12] Zeng Z Y, Hu C E, Li Z G, Zhang W, Cai L C 2015 J. Alloys Compd. 640 201
Google Scholar
[13] Blanco M A, Francisco E, Luaña V 2004 Comput. Phys. Commun. 158 57
Google Scholar
[14] Blanco M A, MartínPendás A, Francisco E, Recio J M, Franco R 1996 J. Mol. Struct.: Theochem 368 245
Google Scholar
[15] Francisco E, Recio J M, Blanco M A, Pendás A M 1998 J. Phys. Chem. 102 1595
Google Scholar
[16] Francisco E, Sanjurjo G, Blanco M A 2001 Phys. Rev. B 63 094107
Google Scholar
[17] Flórez M, Recio J M, Francisco E, Blanco M A, Pendás A M 2002 Phys. Rev. B 66 144112
Google Scholar
[18] 邓世杰, 赵宇宏, 侯华, 文志勤, 韩培德 2017 物理学报 66 146101
Google Scholar
Deng S J, Zhao H Y, Hou H, Wen Z Q, Han P D 2017 Acta Phys. Sin. 66 146101
Google Scholar
[19] Vohra Y K, Holzapfel W B 1993 High Pressure Res. 11 223
Google Scholar
[20] Olsen J S, Gerward L, Benedict U, Itié J P 1985 Physica B+C (Amsterdam)
133 129 Google Scholar
[21] Gu G, Vohra Y K, Winand J M, Spirlet J C 1995 Scr. Metall. Mater. 32 2081
Google Scholar
[22] Svane A 1996 Phys. Rev. B 53 4275
[23] Soderlind P, Eriksson O, Wills J M, Boring A M 1993 Phys. Rev. B 48 9306
Google Scholar
[24] Koskenmaki D C, Gschneidner K A 1978 Handb. Phys. Chem. Rare Earths 1 337
Google Scholar
[25] Gerward L, Olsen J S, Diffr P 1993 Powder Diffr. 8 127
Google Scholar
[26] Decremps F, Antonangeli D, Amadon B, Schmerber G 2009 Phys. Rev. B 80 132103
Google Scholar
[27] Lipp M J, Kono Y, Jenei Z, Cynn H, Aracne-Ruddle C, Park C, Kenney-Benson C, Evans W J 2013 J. Phys: Condens. Matter 25 34
期刊类型引用(0)
其他类型引用(1)
-
图 1 体积随压强变化的规律(黑色实点为直接加压结构优化后的结果, 黑色实线为状态方程拟合结果), 并与已有的Ce[20], Th[19], Ce0.875La0.125[12]的计算值及Ce0.76Th0.24[21]实验值进行比较
Fig. 1. The EOS of fcc and bct Ce-La-Th together with the experimental data (the black solid point is the result of the structure optimization, the black solid line is the fitting result of the EOS), together with the experimental data for Ce0.76Th0.24[21] and the calculated results for Ce[20], Th[19], Ce0.875La0.125[12].
表 1 零温零压下fcc相Ce-La-Th合金的平衡体积(V0)及体积模量(B0)
Table 1. Equilibrium volume (V0) and bulk modulus (B0) of Ce-La-Th of fcc phase at 0 GPa and 0 K.
-
[1] Bridgman P W 1927 Proc. Am. Acad.Arts Sci. 62 207
Google Scholar
[2] Bridgman P W 1951 Proc. Am. Acad. Arts Sci. 79 149
Google Scholar
[3] Bridgman P W 1954 Proc. Am. Acad. Arts Sci. 83 1
[4] Lanson A W, Tang T Y 1949 Phys. Rev. 76 301
Google Scholar
[5] 潘昊, 胡晓棉, 吴子辉, 戴诚达, 吴强 2012 物理学报 61 206401
Google Scholar
Pan H, Hu X M, Wu Z H, Dai C D, Wu Q 2012 Acta Phys. Sin. 61 206401
Google Scholar
[6] Hu C E, Zeng Z Y, Zhang L, Chen X T, Cai L C 2011 Physica B 406 669
Google Scholar
[7] Lawson A C, Williams A, Wire M S 1988 J. Less-common Met. 142 177
Google Scholar
[8] Lawrence J M, Thompson J D, Fisk Z, Smith J L, Batlogg B 1984 Phys. Rev. B 29 4017
Google Scholar
[9] Drymiotis F, Singleton J, Harrison N, Lashley J C, Bangura A, Mielke C H, Balicas L, Fisk Z, Migliori A, Smith J L 2005 J. Phys.: Condens. Matter 17 L77
Google Scholar
[10] Ruff J P C, Islam Z, Das R K, Nojiri H, Cooley J C, Mielke C H 2012 Phys. Rev. B 85 024104
Google Scholar
[11] Hu C E, Zeng Z Y, Zhang L, Chen X R, Cai L C 2010 Solid State Commun. 150 2362
Google Scholar
[12] Zeng Z Y, Hu C E, Li Z G, Zhang W, Cai L C 2015 J. Alloys Compd. 640 201
Google Scholar
[13] Blanco M A, Francisco E, Luaña V 2004 Comput. Phys. Commun. 158 57
Google Scholar
[14] Blanco M A, MartínPendás A, Francisco E, Recio J M, Franco R 1996 J. Mol. Struct.: Theochem 368 245
Google Scholar
[15] Francisco E, Recio J M, Blanco M A, Pendás A M 1998 J. Phys. Chem. 102 1595
Google Scholar
[16] Francisco E, Sanjurjo G, Blanco M A 2001 Phys. Rev. B 63 094107
Google Scholar
[17] Flórez M, Recio J M, Francisco E, Blanco M A, Pendás A M 2002 Phys. Rev. B 66 144112
Google Scholar
[18] 邓世杰, 赵宇宏, 侯华, 文志勤, 韩培德 2017 物理学报 66 146101
Google Scholar
Deng S J, Zhao H Y, Hou H, Wen Z Q, Han P D 2017 Acta Phys. Sin. 66 146101
Google Scholar
[19] Vohra Y K, Holzapfel W B 1993 High Pressure Res. 11 223
Google Scholar
[20] Olsen J S, Gerward L, Benedict U, Itié J P 1985 Physica B+C (Amsterdam)
133 129 Google Scholar
[21] Gu G, Vohra Y K, Winand J M, Spirlet J C 1995 Scr. Metall. Mater. 32 2081
Google Scholar
[22] Svane A 1996 Phys. Rev. B 53 4275
[23] Soderlind P, Eriksson O, Wills J M, Boring A M 1993 Phys. Rev. B 48 9306
Google Scholar
[24] Koskenmaki D C, Gschneidner K A 1978 Handb. Phys. Chem. Rare Earths 1 337
Google Scholar
[25] Gerward L, Olsen J S, Diffr P 1993 Powder Diffr. 8 127
Google Scholar
[26] Decremps F, Antonangeli D, Amadon B, Schmerber G 2009 Phys. Rev. B 80 132103
Google Scholar
[27] Lipp M J, Kono Y, Jenei Z, Cynn H, Aracne-Ruddle C, Park C, Kenney-Benson C, Evans W J 2013 J. Phys: Condens. Matter 25 34
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 11846
- PDF下载量: 139
- 被引次数: 1