姜 濤 田 杰 王 寧 彭述明,* 李 梅 韓 偉 張密林
(1中國工程物理研究院核物理與化學研究所,四川綿陽621999;2哈爾濱工程大學材料科學與化學工程學院,教育部超輕材料與表面技術重點實驗室,哈爾濱150001)
LiCl-KCl-LaCl3熔鹽體系中La(III)的電化學行為
姜濤1田杰1王寧1彭述明1,*李梅2韓偉2張密林2
(1中國工程物理研究院核物理與化學研究所,四川綿陽621999;2哈爾濱工程大學材料科學與化學工程學院,教育部超輕材料與表面技術重點實驗室,哈爾濱150001)
在773 K條件下,研究了La(III)在LiCl-KCl熔鹽中W和Ni電極上的電化學行為。La(III)還原反應是一步三電子轉移的準可逆反應;通過在Ni電極上直接電沉積La的方法可以獲得La-Ni金屬間化合物;恒電位電解可以獲得含三種金屬間化合物(LaNi5、La7Ni16和La2Ni3)的La-Ni合金層,并且通過X射線衍射儀(XRD)和掃描電鏡-能譜分析儀(SEM-EDS)確定物相并表征結構。采用開路計時電位法估算了LaNi5金屬間化合物的標準生成吉布斯自由能。揭示了恒電位電解方法是制備La-Ni鍍層合金以及提取熔鹽中La的有效方法。
LiCl-KCl-LaCl3;La-Ni;鍍層合金;吉布斯自由能
另外,含有鑭系元素和第一長周期過渡金屬元素的金屬間化合物因其有趣的物理化學性質,使得它們成為儲氫和永磁材料制備的基礎原料8。LaNi5和La7Ni16是研究較為廣泛的永磁材料9。而La-Ni合金作為儲氫材料的研究開始于20世紀70年代10。LaNi5合金因其易活化、適中的分解氫壓和低平衡壓等特點,成為優(yōu)異的低溫儲氫材料11-13。并且ABx型鑭系-Ni基儲氫合金擁有更大的儲氫性能,其中La7Ni16是典型的此類儲氫合金14。
近年來,許多學者對熔鹽電解法分離Ln并制備Ln-Ni合金進行了研究。Taxil小組15在LiF-CaF2-LnF3熔鹽中的活性Ni電極上通過電化學還原鑭系元素(Nd、Gd和Sm)分別制備了NdNi2、NdNi3、GdNi2、GdNi、SmNi2和SmNi3等金屬間化合物。Sa?la等16研究了LiF-CaF2熔鹽中在W、Ni和Cu電極上Dy(III)的電化學行為,并通過恒電流電解獲得了Dy2Ni17和DyNi4合金鍍層。Chandra等17在LiFCaF2熔鹽中研究了Ce(III)在W和Ni電極上的電化學氧化還原行為,并且制備了CeNi5和Ce2Ni7兩種金屬間化合物。Ito等18,19分別在LiCl-KCl-SmCl3和LiCl-KCl-YbCl3熔鹽中研究了SmNi5、SmNi3和SmNi2,以及Yb2Ni17和YbNi2在Ni電極上的電化學制備。Yasuda等20在NaCl-KCl-NdCl3熔鹽中研究了Nd-Ni合金的電化學制備,得到了NdNi2、NdNi3、 Nd2Ni7和NdNi5等四種金屬間化合物。
綜上所述,研究利用活性Ni電極來實現熔鹽中La的提取,以及實現La-Ni合金制備都具有重大的意義。為此采用多種電化學技術(循環(huán)伏安、方波伏安、計時電位和開路計時電位)研究了La(III)在LiCl-KCl熔鹽中的電化學行為和La-Ni金屬間化合物的制備機理,并估算了LaNi5金屬間化合物的標準生成吉布斯自由能。
實驗中使用的主要試劑有氯化鋰(分析純,98%)、氯化鉀(分析純,99.5%)、氯化鑭(分析純,99.9%)、銀絲(直徑1 mm,99.99%)、鎢絲(直徑1 mm,99.99%)和鎳棒(直徑7 mm,99.99%),以上試劑和材料均為Alfa Aesar中國化學有限公司產品;無水乙醇(分析純,99.7%)為國藥集團化學試劑有限公司產品;石墨棒(直徑6 mm,光譜純)為上海碳素廠產品;氬氣(99.99%)為成都天一科技有限公司產品。
實驗過程中使用的所有水溶液都采用二次蒸餾水配制。整個電化學實驗在裝備井式電爐(SRH-50-10F,哈爾濱丞焱熱處理設備制造有限公司)和Ar氣保護的手套箱(LS1200,成都德力斯實業(yè)有限公司)中進行。電解質LiCl-KCl裝在電爐中的Al2O3坩堝中。LiCl和KCl(50%:50%(w))在973和573 K溫度下分別干燥24 h以除去殘留的水分。La離子以無水LaCl3的形式進入熔鹽中。所有電化學測量采用Autolab AUT302N型電化學工作站(梅特勒公司,瑞士)和Nova1.10軟件。將一根銀絲插入含AgCl(1%(w))和LiCl-KCl(50%:50%(w))的Al2O3管中作為參比電極。所有稱量工作采用電子天平(AB204-N,上海Mettler-Toledo有限公司)進行。實驗中所有電位都相對于該Ag/AgCl電極。W作為工作電極,石墨作為對電極。實驗前W電極末端用SiC砂紙打磨后在乙醇溶液中用超聲波清洗。有效電極面積通過每次插入熔鹽的電極長度來進行計算。長度用游標卡尺(CD-15CPX,Mitutoyo公司,日本)進行測量。通過恒電位的方法來制備La-Ni合金鍍層。電解后,所有樣品在乙醇溶液中用超聲波清洗以去除殘留的鹽,后放置在干燥箱中待測。
樣品的X射線衍射(X′Pert Pro,Philips公司,荷蘭)表征,以Cu-Kα射線作為射線源,工作電壓40 kV,工作電流40 mA,掃描角度(2θ)范圍為10°-90°,掃描步長0.02°,掃描速率0.06(°)·s-1。樣品用掃描電鏡(JSM-6480,JEOL公司,日本)表征,樣品經過打磨和拋光后在乙醇溶液中用超聲波清洗完畢后進行掃描電子顯微鏡測試。
3.1La(III)在LiCl-KCl-LaCl3熔鹽中Ni電極上的電化學行為
3.1.1循環(huán)伏安
采用循環(huán)伏安法研究了W和Ni電極在LiCl-KCl熔鹽中的電化學窗口。之所以選擇W作為Ni電極的對比電極,是因為根據二元W-Li和W-La相圖,在773 K時W與Li和La都不會形成金屬間化合物21,22。圖1為LiCl-KCl熔鹽中W電極上和加入LaCl3后再W和Ni電極上的循環(huán)伏安曲線。在黑色曲線上,僅觀察到一個陰極還原峰(-2.45 V),可以歸屬為Li的沉積,而對應的氧化峰則是Li的溶解。圖1中綠色曲線代表加入LaCl3后在W電極上的循環(huán)伏安曲線。還原峰A(-2.10 V附近)對應于La(III)還原為La(0),該結果與其他學者的觀察相一致1-3。在反向掃描方向,氧化峰A′可以歸屬為La(0)的氧化溶解,這與Fabian等4的觀察結果相一致。黑色曲線代表LiCl-KCl-LaCl3熔鹽中La(III)在Ni電極上的循環(huán)伏安曲線,與在W電極上的循環(huán)伏安曲線不同。除了陰極/陽極峰A/A′對應于La的還原/氧化峰外,其余在La(III)還原之前的三個峰(-1.38、-1.70和-1.97 V)則對應金屬或者金屬間化合物在W電極上的沉積。這三對峰(I/I′、II/II′和III/III′)分別對應于La-Ni金屬間化合物的生成。但是在Ni電極上未觀察到明顯的A/A′峰。
圖1 LiCl-KCl-(3.5%(w))LaCl3熔鹽中W和Ni電極上的循環(huán)伏安曲線圖Fig.1 Cyclic voltammogram curves of LiCl-KCl melts (red curve)on a W electrode,LiCl-KCl-(3.5%(w))LaCl3melts on W(green curve)and Ni(black curve)electrodes
3.1.2方波伏安
方波伏安法是比循環(huán)伏安法更靈敏的方法,并且能計算轉移電子數。圖2(a)顯示在LiCl-KCl-LaCl3熔鹽體系中的不同頻率的方波伏安曲線。在插圖中頻率的1/2次方與電流強度成線性關系。根據半峰寬公式(1)計算23,24:
計算轉移電子數n≈2.85,接近3,因此我們認為La(III)是一步3電子轉移:
圖2 (a)LiCl-KCl-(3.5%(w))LaCl3熔鹽中不同頻率的方波伏安曲線;(b)LiCl-KCl-(3.5%(w))LaCl3熔鹽中W和Ni電極上的方波伏安曲線圖Fig.2 (a)Square-wave voltammogram curves for the reduction of(a)La(III)on a W electrode;(b)square-wave voltammogram curves of the LiCl-KCl-(3.5%(w))LaCl3molten on W and Ni electrodes
我們采用方波伏安法進一步研究La(III)在Ni電極上的電化學行為。圖2(b)所示為773 K、Ni電極、電位步長1 mV,頻率30 Hz時的方波伏安曲線。圖中明顯出現四個峰(A:-2.10 V,I:-1.97 V,II:-1.70 V,III:-1.38 V),分別對應于La還原和3種La-Ni金屬間化合物,這與之前的循環(huán)伏安結果相一致。并且由于方波伏安曲線更加靈敏,循環(huán)伏安曲線中的峰A可以在方波伏安曲線中識別。
3.1.3計時電位
圖3所示為不同電流密度下LiCl-KCl-(3.5%(w)) LaCl3熔鹽中的計時電位曲線。當電流密度分別達到-0.031 A·cm-2(-10 mA),-0.062 A·cm-2(-20 mA),-0.124 A·cm-2(-40 mA)和-0.186 A·cm-2(-60 mA)時,曲線出現三個電位平臺(III、II和I),對應La在Ni電極上的欠電位沉積并生成La-Ni金屬間化合物。當陰極電流密度達到-0.217 A·cm-2(-70 mA)時,新電位平臺A出現,對應La的沉積。
圖3 不同電流下LiCl-KCl-(3.5%(w))LaCl3熔鹽中Ni電極上的方計時電位曲線圖Fig.3 Chronopotentiograms obtained at different currents on a Ni electrode in the LiCl-KCl-(3.5%(w))LaCl3melts S=0.322 cm2,T=773 K
3.2恒電位電解與La-Ni合金鍍層表征
根據La-Ni二元相圖(圖4)25可知,相圖共有8種La-Ni金屬間化合物,我們實驗中出現3種金屬間化合物。因此,我們采用恒電位電解,根據XRD以及SEM/EDS確定是哪種金屬間化合物。
圖4 La-Ni二元合金相圖25Fig.4 Phase diagram of La-Ni binary alloy25A:atomic fraction;w:mass fraction
基于循環(huán)伏安、方波伏安和計時電位的結果,分別在-2.08、-1.85和-1.55 V(vs Ag/AgCl)電位下,在Ni電極上電解制備La-Ni金屬間化合物。
采用XRD表征在-2.08 V下電解3 h后的樣品a,發(fā)現有三種金屬間化合物存在,分別是La2Ni3、La7Ni16和LaNi5(如圖5a所示)。通過SEMEDS分析樣品La和Ni的元素分布,發(fā)現獲得的LaNi5鍍層(圖6a中點A)最接近Ni電極,很致密但非常薄,厚度大約為10 μm。中間層La7Ni16(圖6a中點B)不明顯而且多孔,厚度大約為50 μm。最外層主要是La2Ni3(圖6a中點C),均勻且致密,厚度大約為75 μm。EDS能譜分析結果表明,鍍層主要由Ni和La兩種元素組成。三個選點(A、B和C)的Ni/La原子比分別為6.31、2.29和1.87。當考慮測量誤差后,我們進一步確定分別為LaNi5、La7Ni16和La2Ni3相金屬間化合物。而在-1.85 V下電解3 h后的樣品b中發(fā)現有兩種金屬間化合物存在,分別是La7Ni16和LaNi5(如圖5b所示)。通過SEM-EDS分析樣品La和Ni的元素分布,發(fā)現獲得的LaNi5鍍層(圖6b中點A)最接近Ni電極,很致密、均一且較厚,厚度大約為20 μm。最外層主要是La7Ni16(圖6b中點B),疏松且很厚,厚度大約為30 μm。EDS能譜分析結果表明,鍍層主要由Ni和La兩種元素組成。兩個選點(A和B)的Ni/La原子比分別為8.00和2.06。當考慮測量誤差后,我們進一步確定分別為LaNi5和La7Ni16相金屬間化合物。
圖5 LiCl-KCl-(3.5%(w))LaCl3熔鹽中Ni電極上不同沉積電位下獲得樣品的XRD譜圖Fig.5 XRD patterns of deposit samples obtained by potentiostatic electrolysis at different potentials on a Ni electrode in the LiCl-KCl-(3.5%(w))LaCl3melts
圖6 鎳電極上不同沉積電位(P)下獲得樣品的SEM-EDS譜圖Fig.6 SEM-EDS analysis of the sample obtained by potentiostatic electrolysis at different potentials(P)on a Ni electrode
在-1.55 V下電解3 h后的樣品c中發(fā)現僅有一種金屬間化合物存在,即LaNi5(如圖5c所示)。通過SEM-EDS分析樣品La和Ni的元素分布,發(fā)現獲得的LaNi5鍍層(圖6c中點A)緊挨Ni電極,很致密、均一且很厚,厚度大約為35 μm。EDS能譜分析結果表明,鍍層主要由Ni和La兩種元素組成。一個選點(A)的Ni/La原子比為5.48。當考慮測量誤差后,我們進一步確定為LaNi5相金屬間化合物。
3.3La-Ni金屬間化合物生成的熱力學研究
采用開路計時電位曲線進一步研究了La-Ni金屬間化合物生成的熱力學。根據以上XRD以及SEM-EDS分析可知道,該濃度下的循環(huán)伏安和方波伏安曲線只出現3個I、II和III三種金屬間化合物,并且恒電位電解只獲得3種產物。我們可以確定I、II和III三種金屬間分別是La2Ni3、La7Ni16和LaNi5。根據相圖可知,LaNi5是最富Ni的金屬間化合物,我們可以根據開路計時電位計算LaNi5的生成吉布斯自由能26。
采用開路計時電位曲線進一步研究了La-Ni金屬間化合物生成的熱力學。在LiCl-KCl-(3.5%(w)) LaCl3熔鹽中Ni電極上,-2.50 V沉積50 s,獲得La-Ni合金的一個薄層。圖7為773 K條件下,恒電位的開路計時電位曲線。我們可以清楚地觀察到4個平臺,平臺A表示La(III)/La的平衡。而平臺I、II和III則分別表示La2Ni3、La7Ni16和LaNi5三種金屬間化合物。
電位平臺的平衡電位可以轉化對于La(0)的電動勢,可以估算LaNi5金屬間化合物的生成吉布斯自由能26,數據如表1所示。從表中數據可以看出,本文中獲得的生成吉布斯自由能要大于Dischinger27和Lantelme6等的實驗數據。分析原因認為:(1)實驗溫度不同;(2)平臺對應的電化學反應不同;(3)對平臺的確定和取值不同。(4)測量方法不同。
圖7 電沉積后LiCl-KCl-(3.5%(w))LaCl3熔鹽中Ni電極的開路計時電位曲線圖Fig.7 Open-circuit potential transient curve for a Ni electrode after electrodepositing
表1 LaNi5金屬間化合物的生成吉布斯自由能數據Table 1 Gibbs energies of formation for LaNi5intermetallic compounds
通過多種電化學方法研究了773 K條件下La (III)在LiCl-KCl熔鹽中W和Ni電極上的電化學行為。根據循環(huán)伏安、方波伏安和開路計時電位法的研究可知,通過在Ni電極上直接電沉積La的方法可以獲得La-Ni金屬間化合物。恒電位電解可以獲得含三種金屬間化合物(LaNi5、La7Ni16和La2Ni3)的La-Ni合金層,并且通過XRD和SEM-EDS確定物相并表征結構。采用電動勢方法估算了LaNi5金屬間化合物的生成吉布斯自由能。電解方法可成為制備La-Ni金屬間化合物以及提取熔鹽中La的有效方法。
References
(1)Castrillejo,Y.;Bermejo,M.R.;Martinez,A.M.;DiazArcas,P. J.Min.Metall.Sect.B:Metall.2003,39(102),109.
(2) Masset,P.;Konings,R.J.M.;Malmbeck,R.;Serp,J.;Glatz,J. P.J.Nucl.Mater.2005,344(1-3),173.doi:10.1016/j. jnucmat.2005.04.038
(3) Vandarkuzhali,S.;Gogoi,N.;Ghosh,S.;Prabhakara Reddy,B.; Nagarajan,K.Electrochim.Acta 2012,59,245.doi:10.1016/j. electacta.2011.10.062
(4) Fabian,C.P.;Luca,V.;Chamelot,P.;Massot,L.;Caravaca,C.; Lumpkin,G.R.J.Electrochem.Soc.2012,159(4),F63. doi:10.1149/2.057204jes
(5) Qiao,H.;Nohira,T.;Ito,Y.J.Alloy.Compd.2003,359(1-2), 230.doi:10.1016/S0925-8388(03)00203-2
(6) Lantelme,F.;Cartailler,T.;Berghoute,Y.;Hamdani,M. J.Electrochem.Soc.2001,148(9),C604.doi:10.1149/ 1.1385819
(7)Han,W.;Zhang,Y.;Ye,K.;Yan,Y.;Zhang,M.Metall.Mater. Trans.B 2010,41(5),1123.doi:10.1007/s11663-010-9395-z
(8) Rong,Q.;Schaller,H.J.J.Alloy.Compd.2004,365(1-2),188. doi:10.1016/S0925-8388(03)00667-4
(9) Hsu,S.E.;Beibutian,V.M.;Yeh,M.T.J.Alloy.Compd.2002, 330-332,882.doi:10.1016/S0925-8388(01)01669-3
(10) Srivastava,S.;Srivastava,O.N.J.Alloy.Compd.1999,290(1-2),250.doi:10.1016/S0925-8388(99)00196-6
(11) Prigent,J.;Joubert,J.M.;Gupta,M.J.Solid State Chem.2011, 184(1),123.doi:10.1016/j.jssc.2010.10.037
(12) Dhaou,H.;Askri,F.;Ben Salah,M.;Jemni,A.;Ben Nasrallah, S.;Lamloumi,J.Int.J.Hydrog.Energy 2007,32(5),576. doi:10.1016/j.ijhydene.2006.07.001
(13)Ao,B.Y.;Chen,S.X.;Jiang,G.Q.J.Alloy.Compd.2005,390 (1-2),122.doi:10.1016/j.jallcom.2004.05.092
(14) Kandavel,M.;Bhat,V.V.;Rougier,A.;Aymard,L.;Nazri,G. A.;Tarascon,J.M.Int.J.Hydrog.Energy 2008,33(14),3754. doi:10.1016/j.ijhydene.2008.04.042
(15) Chamelot,P.;Massot,L.;Hamel,C.;Nourry,C.;Taxil,P. J.Nucl.Mater.2007,360(1),64.doi:10.1016/j. jnucmat.2006.08.015
(16) Sa?la,A.;Gibilaro,M.;Massot,L.;Chamelot,P.;Taxil,P.; Affoune,A.M.J.Electroanal.Chem.2010,642(2),150. doi:10.1016/j.jelechem.2010.03.002
(17) Chandra,M.;Vandarkuzhali,S.;Ghosh,S.;Gogoi,N.; Venkatesh,P.;Seenivasan,G.;Reddy,B.P.;Nagarajan,K. Electrochim.Acta 2011,58,150.doi:10.1016/j. electacta.2011.09.012
(18) Iida,T.;Nohira,T.;Ito,Y.Electrochim.Acta 2001,46(16), 2537.doi:10.1016/S0013-4686(01)00470-4
(19) Iida,T.;Nohira,T.;Ito,Y.Electrochim.Acta 2003,48(11), 1531.doi:10.1016/S0013-4686(03)00031-8
(20) Yasuda,K.;Kobayashi,S.;Nohira,T.;Hagiwara,R. Electrochim.Acta 2013,92,349.doi:10.1016/j. electacta.2013.01.049
(21) Sangster,J.;Pelton,A.D.J.Phase Equilib.1991,12(2),203. doi:10.1007/BF02645715
(23) Ramaley,L.;Krause,M.S.Anal.Chem.1969,41(11),1362. doi:10.1021/ac60280a005
(24) Osteryoung,J.G.;Osteryoung,R.A.Anal.Chem.1985,57. doi:10.1021/ac00279a789
(25) Okamoto,H.J.Phase Equilib.2002,23(3),287. doi:10.1361/105497102770331884
(26) von Barner,J.H.;Noyé,P.;Barhoun,A.;Lantelme,F. J.Electrochem.Soc.2005,152.doi:10.1149/1.1833313
(27) Dischinger,J.;Schaller,H.J.J.Alloy.Compd.2000,312(1-2), 201.doi:10.1016/S0925-8388(00)01151-8
Electrochemical Behavior of La(III)in Molten LiCl-KCl-LaCl3
JIANG Tao1TIAN Jie1WANG Ning1PENG Shu-Ming1,*LI Mei2HAN Wei2ZHANG Mi-Lin2
(1Institute of Nuclear Physics and Chemistry,China Academy of Engineering Physics,Mianyang 621999,Sichuan Province,P.R. China;2Key Laboratory of Superlight Materials and Surface Technology,Ministry of Education,College of Materials Science and Chemical Engineering,Harbin Engineering University,Harbin 150001,P.R.China)
Electrochemical and thermodynamic studies on the formation of La-Ni intermetallic compounds in molten LiCl-KCl-(3.5%(w))LaCl3at 773 K were performed.The electrochemical reduction of La(III)ions was investigated on inert W and reactive Ni electrodes by cyclic voltammetry.The reduction potential of La(III)/La on a Ni electrode was observed at more positive potential than that on a W electrode because of the formation of La-Ni intermetallic compounds when La ions reacted with the Ni substrate.Square-wave voltammetry, chronopotentiometry,and open-circuit chronopotentiometry provided further evidence for the formation of La-Ni intermetallic compounds.Potentiostatic electrolysis on a Ni electrode led to the formation of three La-Ni intermetallic compounds,LaNi5,La7Ni16and La2Ni3,according to X-ray diffraction(XRD)and scanning electron microscopy-energy-dispersive X-ray spectroscopy(SEM-EDS)analyses.The standard Gibbs free energies of formation for LaNi5were estimated from open-circuit chronopotentiometric measurements using the Gibbs-Helmholtz equation and Hess law.The formation of the La-Ni alloy layer could be controlled by the applied potential and time.Potentiostatic electrolysis was an effective method for electrochemical extraction of La.
LiCl-KCl-LaCl3;La-Ni;Alloy layer;Gibbs free energy
鑭(La)是鑭系中原子量最小的元素,它的中子吸收截面較大,能夠阻礙可嬗變錒系元素俘獲中子,降低嬗變效率,因此在乏燃料干法后處理中必須考慮它的分離提取問題。對于La(III)在惰性電極上的電化學行為進行了許多研究,包括LiCl-KCl1-4、CaCl2-NaCl1和LiF-NaF-KF-LaF32熔鹽體系。這些實驗結果表明,在LiCl-KCl熔鹽中,La (III)的還原是一步三電子轉移反應3,4。Castrillejo等1利用循環(huán)伏安、計時電位和計時電流技術探討了La(III)的反應機理、成核現象和電活性離子的輸運參數。Vandarkuzhali等3在LiCl-KCl熔鹽中活性Al電極上利用瞬態(tài)電化學技術研究了La(III)的電化學行為,獲得了Al11La3金屬間化合物,并利用開路計時電位曲線計算了其生成熱力學參數(吉布斯自由能)和La(III)的擴散系數。Ito小組5在LiFNaF-KF-LaF3熔鹽中Pd電極上研究了Pd3La合金鍍層的電化學制備。Lantelme等6利用瞬態(tài)電化學技術研究了LiCl-KCl熔鹽中鑭系(La、Ce、Pr、Gd和Y)離子的熱力學和動力學性質,并計算了La-Ni和Gd-Ni金屬間化合物的生成吉布斯自由能。張密林課題組7在LiCl-KCl-MgCl2-KF-Ln2(CO3)3熔鹽中電解制備了Mg-Li-Ce-La四元合金。
March 18,2016;Revised:June 21,2016;Published online:June 22,2016.
.Email:pengshuming@caep.cn;Tel:+86-816-2493835.
O641;O643
10.3866/PKU.WHXB201606223
The project was supported by the Major Research Plan of the National Natural Science Foundation of China(91426302)and National Natural Science Foundation of China(21301163).
國家自然科學基金重大研究計劃-集成項目(91426302)和國家自然科學基金(21301163)資助?Editorial office ofActa Physico-Chimica Sinica
(22) Okamoto,H.J.Phase Equilib.2001,22(6),693.10.1007/ s11669-001-0048-4