从第一性原理赝势构建稠密等离子体中内禀屏蔽重粒子势

《Atoms》：Construction of Intrinsically Screened Heavy-Particle Potentials in Dense Plasmas from First-Principles Pseudopotentials

　　

研究人员提出一个计算框架，用于从第一性原理赝势构建稠密等离子体中有效的、内禀屏蔽的重粒子势。从统一赝势格式（Unified Pseudopotential Format, UPF）中可用的投影缀加波（projector-augmented-wave, PAW）数据出发，重建角动量相关的通道势（s、p、d、f），并根据等离子体温度下的Boltzmann占据概率进行混合。通过为每个电荷态计算自洽Hartree势，纳入热占据价电子的内禀屏蔽，从而产生元素特异、温度依赖的有效散射势，用于中性原子和离子。经Stewart–Pyatt模型修正电离势降低（ionisation-potential depression, IPD）的Saha电离平衡决定物种组成。数密度加权的混合物平均产生多组分等离子体的单一有效重粒子势。该程序在D_plas_V中实现，这是一个零外部依赖的C++17代码，且所有方法均针对解析测试用例进行了验证。

在温稠密等离子体（如惯性约束聚变、行星内部和恒星大气）中，重粒子（中性原子和离子）间的有效相互作用势既区别于裸Coulomb势，也不同于孤立原子赝势。自由电子会屏蔽离子芯，而在高温高密度条件下，常规离子间多极展开失效，必须使用完整离子势。第一性原理量子分子动力学（quantum molecular dynamics, QMD）虽能自洽计算这些势，但计算成本极高，难以进行宽范围温度、密度和化学组成的参数扫描。基于截断Coulomb势的简化模型虽已用于天体等离子体不透明度计算，却缺乏元素特异性和角动量分辨。因此，研究人员旨在利用现成的第一性原理赝势，发展一种介于全QMD与解析屏蔽Coulomb势之间的高效方法，以构建元素特异、角动量分辨、温度依赖的重粒子有效散射势。该研究发表在《Atoms》。

研究人员利用DFT库中已有的PAW赝势，从JTH v1.1验证表获取UPF格式的赝势数据，重建s、p、d、f角动量通道势；在等离子体温度下按Boltzmann占据概率混合通道势；对每个电荷态求解径向Schr?dinger方程和球Poisson方程，得到价电子自洽Hartree势，纳入内禀屏蔽；采用Saha电离平衡并施以Stewart–Pyatt电离势降低（IPD）修正确定物种组成；最后进行数密度加权混合平均，得到单一有效重粒子势。该方法在C++17程序D_plas_V中实现，并已通过解析测试用例验证。该研究为稠密等离子体输运和光学计算提供了一个计算高效、物理动机清晰的势构建框架。

主要关键技术方法包括：从UPF文件解析PAW赝势数据，重建局域势和非局域投影子的角动量通道势；基于Saha方程并耦合Stewart–Pyatt IPD修正求解多组分电离平衡；采用Boltzmann占据概率对通道势进行温度热平均；利用Numerov方法求解放射状Schr?dinger方程，获得束缚态本征值和波函数，构造价电子密度，再通过球Poisson方程求解Hartree势；最后按数密度加权平均所有带电物种的势。测试样本为90% Ar–10% H混合等离子体，总重粒子数密度和自由电子数密度由Saha解自洽确定。

4.1 Saha组成与物种平均势：对90% Ar–10% H等离子体求解Saha电离平衡，得到六种主要重粒子物种，其中Ar⁺占据绝对主导，其次是H⁺和Ar²⁺，中性Ar和H贡献很小，Ar³⁺可忽略。每种物种的势经通道重建、Boltzmann热平均和Hartree屏蔽独立计算，再按数密度加权得到平均重粒子势。结果表明，平均势主要由Ar⁺的贡献决定，其他物种在中距离上产生一定偏移。各物种势随电离度增加而系统性加深，因为剩余价电子减少、屏蔽变弱，且在大距离处趋于各自的Coulomb渐近尾部。

4.2 与解析Debye屏蔽Coulomb势的比较：将计算得到的重粒子平均势与以平均离子电荷为中心的Debye屏蔽Coulomb势进行对比。在大距离处两者趋于一致；在中间距离处，重粒子平均势明显更深，体现了物种特异性Hartree屏蔽信息；在短距离处，核排斥芯主导。该偏差表明，输运计算需纳入物种特异性屏蔽效应，并可通过冲击波实验的稠密等离子体电导率数据加以检验。

讨论部分指出，上述与Debye模型的偏差应对应于电导率的修正，且可用Fortov–Mintsev小组的冲击波电导率测量作为实验基准。研究人员还计划在强耦合条件下引入Ebeling–F?rster–Fortov离子关联修正IPD形式，并基于Vitel团队与L. G. D’yachkov合作发展的理论框架，将结果与动态介电函数和频率相关输运性质联系起来。

翻译研究结论部分：研究人员已开发了一个计算框架，从现成的第一性原理PAW赝势出发构建稠密等离子体中有效的内禀屏蔽离子势。该过程包括从UPF数据重建角动量相关通道势，应用Boltzmann加权热平均以捕捉温度依赖的电子构型，以及通过Hartree势计算来自价电子的内禀屏蔽修正。Saha电离平衡结合Stewart–Pyatt IPD修正确定物种组成，数密度加权平均所有电荷态得到单个重粒子平均势。该方法提供了一种计算高效的框架，以物理动机明确的方式构建包含内禀屏蔽效应的有效离子势，其结构随着底层物理要素（通道分辨赝势、Boltzmann混合、Saha平衡和Hartree屏蔽）自洽地纳入而趋近于真实散射势。对于90% Ar–10% H测试案例，重粒子平均势在中间距离与解析Debye屏蔽Coulomb势显著不同，表明解析物种特异性屏蔽贡献的重要性。未来需系统研究角动量通道数对平均势的敏感性，并与QMD模拟和实验数据全面验证；同时，强耦合条件下的离子关联修正IPD形式以及与动态介电函数和频率相关输运的联系，将作为后续扩展方向。