LiF-KCl熔盐溶液的MonteCarlo法计算机模拟研究-上海有机化学.PDF

化学学报 AarA OHIMICA 缸NICA 1189, n， ω-5组 LiF-KCl 熔盐溶液的 Monte Carlo 法 计算机模拟研究 II. 熔体结构与物性的关系 徐驰势 江乃雄 陈念胎 〈中国科学院上海冶金研究所，上海) 本文用 Monte Callo 法对互易盐系LiF-KCl 熔盐溶液的局部结构作了计算机模拟. 计算了该 熔盐溶液的总势能和势能分布.在 1200K 模拟温度下，该熔盐溶液中大部分离子组成各种形式的 离子团.根据模拟的互易盐系熔盐溶液模型，讨论了熔体局部结构和物性之间的关系. 在前文(1) 已经报道 LiF-KCl 熔盐洛液的径向分布函数和热力学性质.本文继续报道 Mon切 Carlo 法计算机模拟结果，通过 LiF-KCl 熔盐溶液总势能的计算和熔体局部结构的分 析，讨论了熔体局部结构和物性之间的关系. 凸 能 ，巳品、 势 为研究 KCI-LiF熔盐榕液结构，共模拟计算了温度相同而体积不同的 5 个V， T 点.每 5 个V， T 点一般都进行了 2-4x 10 次 MC 运算. 为便于与已发表的纯态KC1，且F 熔盐的 计算机模拟结果作比较， LiF-KCI 熔盐溶液的模拟温度亦取 1200K. 全部运算中，离子沿 ø , y , Z 方向最大的步长是士0.5Å. 详细算法见文献[1] . 5 个 V， T 点计算的总势能、总内能、维里值及压力等列于表1. 其中总内能为动能与势 能E之和: U = (3/2)NKT+E) 压力 P 根据系统的维里值〈功〉计算: P=NKT/V一〈ψ/3V 而 H，=U+PV 表1 同时列出不同作用力对总势能的贡献， EC， ER及EDD分别代表静电势，近程排 fF势及 Van der Waals 力引起的色散能. 由表1 的 P， V 值可在P， V 曲线上定出当 P=O 时 LiF-KC1 熔盐榕液的摩尔体积V. 向理，亦可定出当 p=o 时各项热力学数据.当 V=32.52o皿3/皿01 时，其P接近 o. 本文以 此状态下模拟结果作为研究对象. 考察 L沪， F- , K飞 Cr-诸离子在 LiF-KCI 熔盐榕被中的势能分布(图 1) ，发现Li飞 F­ 的平均势能较低，而 K飞 cr 的平均势能较高.可见在该熔盐榕液中， K+ 和 cr 处于较不稳 定(或活跃)状态，而半径较小的Li飞 F- 离子由于组成较紧凑的离子集团，相对比较稳寇. 1987 年6 月 1 日收到.前一报见文献[lJ. 上海计算技术研究所. .1530. 化学擎报 ACTA CHIMICA SINICA 198~ 轰1 能量、庄为计算结果停 3 U 〈ψ〉 E 飞E) EJ) ER) EDD) (cmSj皿01) I (10 bar) 3 1.45 I 1.岳6 -789.5