统计自洽场INFERNO模型中电子共振态-计算物理.pdfVIP

第32卷第3期 计 算 物 理 V01．32，No．3 CHINESE OFCOMPUTATIONALPHYSICS 2015年5月 JOURNAL May，2015 文章编号：1001—246X(2015)03-0361-08 统计自洽场INFERNO模型中电子共振态 马桂存， 齐 进， 王 敏 (北京应用物理与计算数学研究所，北京 100094) 电子共振态，给出共振能量区间内电子的态密度，波函数，相移因子随能量的变化及位相一振幅形式的波函数中的 位相随能量的变化，得到与非相对论情况下完全一致的结果． 关键词：INFERNO模型；电子共振态；态密度；散射相移因子 中图分类号：052；0413．1；0482文献标志码：A 0 引言 物质在高温高压下的状态方程是天体物理中恒星构造，地球物理中地核运动，激光惯性约束聚变问题及 核武器设计等领域中不可缺少的重要参数…．在惯性约束聚变数值模拟中所使用的物态方程其温度、密度 范围都很宽广，温度从室温到几千万度，压强从常压到几千万大气压． 在极高温度和极高密度下(压强P10TPa)，托马斯费米理论是生成高温高压物态方程的重要理论模 型．在较低的温度和密度下(压强P1TPa)，经验物态方程的理论模型比较多，如德拜模型或德拜模型加非 谐效应修正，其次还可以采用第一原理方法计算物态方程．在千万大气压致托马斯费米理论适用的范围内， 物质中的原子实际上是处于部分离化的状态，这种状态在物态方程的研究中被称之为过渡区．为了能够计算 这个区域的物态方程，人们提出了很多的物态方程理论模型，其中最有代表性的是统计自洽场INFERNO模 型旧1．该模型考虑了束缚态电子的压致离化和温致离化效应，能够比较好地适用于过渡区． 中加入一个原子而形成的平均原子模型．在这个模型中，电子分为束缚态，散射态，共振态．电子的状态由求 解Dirac方程获得；高能量和高角动量的电子态采用托马斯费米近似计算旧J，其计算结果不能很好地保证物 电子态，得到了很好的结果．但对计算束缚态并不适用．2009年法国原子能科学研究中心的科学家提出了采 用大规模并行计算的方法H1来克服高能量积分、高角动量求和的困难．2011年俄罗斯科学家也用类似的方 子要从束缚态变为巡游态，在这个过程中电子要发生共振，即其态密度有尖峰出现．在某些温度密度点处，电 子态密度能谱的峰值非常尖锐，峰高达到107以上，峰半宽度达到10—7以下，这给计算带来了一定的困难，有 时会影响迭代过程的收敛性．因此，有必要开展INFERNO模型中电子共振态的研究． 电子受到势场的散射作用而发生共振．Sidky等人M1通过解薛定鄂方程，只对非相对论情况下的电子共 振态特点作了简要分析．在相对论情况下电子共振态是否仍然具有这些特点，及如何描述他们，都是一个值 得研究的问题．在本文中，我们要求解的方程就是相对论形式的狄拉克方程，且必须采用自洽场形式的求解 方法，即势场和波函数在初始时都是未知的，通过自洽迭代的求解方法才能获得有用的结果，每一个迭代循 环过程所需的计算量都相当大．因此，开展INFERNO模型中电子共振能态特征研究，对提高计算的效率，改 善计算的精度都具有十分重要的意义． 收稿日期：2014—02—23；修回日期：2014—06—05 362 计 算 物 理 第32卷 1 INFERNO理论模型 在通常的固体结构中，带正电的离子按照一定的方 式排列，在每个离子周围又有很多电子，它们代表核芯 一 电子与价电子(或巡游电子)．核芯电子只能在束缚它 B 的原子核周围运动，而巡游电子属于整个系统，可以在 物质内部自由运动．在量子自洽场统计模型中，认为正 扫一害Qa R Radius 离子是均匀分布在整个空间的，同时，为了保持系统的