稀土金属氯化物熔体结构的MD计算.docVIP

稀土金属氯化物熔体结构的MD计算 单纯氯化稀土的熔点较高，粘度和挥发性较大，本身不稳定、容易同水和氧作用，尤其是在熔融状态下对稀土金属溶解度很大，因此不能仅用它作为电解质进行电解。实际生活中，熔盐电解质都是两元或多元熔体，但在稀土熔盐以及含稀土的混合熔盐中，稀土金属氯化物（RECL3）是最简单的离子熔体，研究稀土金属氯化物系的微观结构和其内在规律，对研究实际中含RECL3的两元或多元熔体有不容忽视的作用。 目前，对稀土金属氯化物熔体结构的研究不多，Tatlipinar等用超网链方程计算了若干镧系金属三价氯化物的偏径向分布函数等结构信息。Saboungi等和Mochinaga等分别用中子衍射和X射线方法测定了YCL3和一系列镧系金属三价氯化物的结构。上海大学黄世萍等用分子动力学方法模拟了YCL3熔盐结构，昆明理工大学侯怀宇等用分子动力学方法模拟了LaCL3、GdCL3和YCL3熔盐结构。 1 模拟设置 模拟中采用的势函数形式如下： （3—1） 式中Zi为离子电荷数，为离子的有效硬度参数，Ri是离子有效半径，f是常数。根据文献，计算所取势函数参数值为：ZRE=+3，Zcl=-1，f=1.153×10J/m，ρ=0.007nm, ρ=0.171nm,稀土元素三价离子半径列于表3—1.其它模拟条件根据文献选取，列于表3—2. 本模拟中，元胞中包含864个离子其中有216个RE离子和648个Cl离子。元胞中尺寸由熔体密度和元胞中所含离子数算出。本模拟中，粒子初始位置和速度均随机选取，计算步长2fs，程序经运行10000步以后，体系的动能，压力，温度等值已趋于稳定，即认为达到平衡，在弛豫2000步，统计各结构参数。 2 偏径向分布函数 计算所得几种熔盐的偏径向分布函数绘于图3—1.从偏径向分布函数曲线中可以得到如图表3—3所示的特征参数，其中r，r分别表示第一峰位置，第一谷底位置，h，d分别表示第一峰高度和离子间最近距离。 图3—1 各熔盐的偏径向分布函数 表3—3 径向分布函数特征参数 从图3一1及表3我们可以较清楚的看出: 氯化稀土熔盐的RE{l径向分布函数主峰高度随RE卜离子半径的减 小有增大的趋势，LaCI:熔盐的主峰最低，这表明随RE，+离子半径的减小， 氛化稀土熔盐中RE卜离子附近在平均距离处出现Cl一离子的机率随之增 大.换言之，随着RE3+离子半径的减小，Cl一离子在REa+离子周围分布越 来越集中，即熔盐中正负离子分布越来越紧凑。Ycl3熔盐的偏径向分布 函数主峰要比Dyc13的低一些，这可能是因为YCI，和DyC13本身的离子半 径就很接近而Ycl，的熔点又比Dycl，的高很多的缘故。 各熔盐的偏径向分布函数第一峰位置、第一谷底位置以及离子间最 近距离(d)随RE3+离子半径的减小有减小的趋势，这可能是由于REs+离 子半径的减小使得离子间距离也随之减小.径向分布函数g+一(r)的第一 谷底位置(r.，。)中YCI:的数值略高一些，这可能是离子半径和温度、密 度共同作用的结果，这是因为正负离子间的作用主要是库仑作用，温度 相对高的熔盐，相应的距离较远(处于第一配位层边缘)的正负离子间 库仑作用相对弱一些，导致Yc13熔盐第一谷底的位置要比Gdcl3和Dyc13 的还大。另外，径向分布函数g--一(r)的氛氛离子间最近距离都相同，这 是由于在本来结构和性质就很接近的氯化稀土熔盐中，作为同一种离子 的抓离子间的作用性质和形式相同的缘故。 为了更清楚地比较各熔盐的偏径向分布函数，我们把径向分布函数 g++(r)、g+一(r)、冬一(r)分别绘于图3一2、图3一3、图3一4。 由图3一2可以看出，LaC13熔盐中正离子间的分布较为混乱松散，径 向分布函数肠(r)第二峰相对不明显，随着REa+离子半径的减小，这种浪 乱程度逐渐减小，g++(r)第二峰相对明显而且峰值增大，这就得出和前面 相似的结果，即随着REa+离子半径的减小，各熔盐中正离子间分布越来 越集中、紧凑。 各熔盐的径向分布函数g+一(r)、g-一(r)分布较为相似，这可能是由 于稀土元素性质较为相似，故和CI一的径向分布相似，而负离子同为一 种离子，故g--一(对分布也较为相似。 3 配位数曲线 计算所得配位数曲线如图3一5所示。表3一4中的配位数来自相应于偏径向分布函数第一谷底(几1。)处的积分值，即粒子第一配位层内的配位数。 从图3一5及表3一4我们可以看出: 随着RE卜离子半径的减小，配位数曲线的平台部分位置逐渐降低而且渐趋平缓。平缓意味着在第一配位层内，Cl一离子在RE卜离子径向分布的几率变化小，Cl一离子离R扩离子原子核距离比较一致，也即Cl一离子在平均位置(‘)出现的几率大一些，分布比较紧凑一些。 4 健角分布 计算方法是根据