有限电子温度下铝中质子的电子阻止本领

早期研究在较低温度（<500 K）下发现电子温度对铝的电子阻止本领影响可忽略，经典理论认为，电子温度升高会增加电子热导率和碰撞频率，从而应导致阻止本领增加。然而，初步研究表明可能存在与之相反的量子效应主导的趋势，这需要更深入的理论阐释。已有研究表明，利用高强度激光可以将固体铝等温加热到极高的电子温度（例如300 eV），产生暖密或热密物质。在此极端条件下，电子温度对阻止本领的影响机制尚不明确，缺乏从第一性原理出发的微观理解。

时间为 t =0 时铝晶体和质子的位置（灰色和蓝色的球分别代表质子和铝原子核）

因此，本研究旨在利用PWmat特色功能——含时密度泛函理论（TDDFT）与非绝热分子动力学（NAMD）相结合的第一性原理方法，系统研究有限电子温度（300 K 至 50,000 K） 对质子（速度0.1-1.0 a.u.）在铝中电子阻止本领的影响，借助极高的空间与时间分辨率，提供了阻止过程中所涉及电子态的丰富详细信息，揭示了其背后的微观物理机制，以弥补该领域的知识空白，并为潜在的材料性质调控提供理论依据。

研究结果表明，电子温度会对电子阻止本领产生显著影响。为解释电子阻止本领随靶材电子温度升高而降低的现象，提出了一种基于能级电子跃迁物理图像的量子阻塞机制。文章发表于《Physical Review B》，DOI：10.1103/PhysRevB.111.144303

研究过程与结果

研究团队设计了一套严谨的多尺度计算方法，将量子力学精度的电子结构与经典的离子运动相结合。具体而言，他们采用了基于Ehrenfest动力学的混合框架，其中电子的行为通过含时密度泛函理论（TDDFT）进行描述，而质子（入射离子）和铝原子核的运动则通过经典分子动力学进行模拟。这种方法的优势在于能够自洽地处理质子与铝电子云之间的能量交换，从而实时计算出质子的能量损失。

研究过程始于构建一个代表铝晶体的周期性超胞模型。为了模拟不同的电子温度（从低温0.026 eV到高达10 eV，覆盖了从常温到强激光激发后的状态），研究人员引入了PWmat的特色功能（TDDFT），在此框架下，电子不再全部处于最低能态，而是按照费米狄拉克分布占据各能级，电子温度越高，高能级被占据的概率就越大。随后，让具有不同初始速度（0.1至1.0原子单位）的质子沿铝晶体的特定通道方向入射，并在TDDFT+MD的框架下进行动力学演化，精确追踪每个时间步长内质子动能的减少，最终计算出平均的电子制动能。

不同速度的质子的动能损失

模拟结果揭示了一个与经典直觉相悖但具有清晰量子根源的核心现象：即电子阻止本领会随着电子温度的升高而显著降低。研究团队将这一反常现象归因于他们提出的 “量子阻塞”（Quantum-Blocking）机制——也就是质子损失能量的主要途径是激发铝中的电子，使其从低能态跃迁到未被占据的高能态。

铝晶体对运动质子的电子阻止本领随速度的变化关系

在低电子温度下，低能态被电子占据，而高能态基本空置，为电子跃迁提供了广阔的通道，因此质子能量损失大，制动能值就高。然而，当电子温度升高时，部分电子会热激发到高能态，导致这些目标态被占据，从而阻塞了低能态电子向这些态的跃迁。可供跃迁的通道减少，质子能够激发的电子数量也随之减少，最终表现为阻止本领的下降。“量子阻塞”机制成功地将宏观观测到的阻止本领下降与微观的电子占据统计联系了起来。

电子阻止本领kB τ 的函数

此外，研究还发现了速度依赖的复杂性。对于高速质子（如1.0 a.u.），电子的阻止本领随温度单调下降的趋势最为明显，“量子阻塞”机制占主导地位。而对于低速质子（如0.1 a.u.），电子阻止本领随温度的变化曲线中出现了一个微小的峰值。

进一步的分析表明，这源于在特定中间温度下，质子与铝靶之间的电荷转移效应暂时增强。通过分析投影态密度（PDOS），研究人员发现此时质子的能级与铝的导带之间共振增强，导致质子更有效地捕获周围电子，短暂地增加了相互作用强度，从而引起了制动能的异常升高。这一发现表明，在低能区，阻止本领是“量子阻塞”机制导致的下降趋势与电荷转移效应引起的局部增强共同竞争的结果。

综上所述，该研究通过先进的量子-经典模拟，不仅定量地描绘了电子温度对阻止本领的抑制效应，更重要的是从量子力学层面提出了“量子阻塞”这一创新物理机制作为根本解释，并揭示了在低速情况下电荷转移所带来的附加影响，为理解和主动调控能量粒子在激发态材料中的输运行为提供了深刻的理论基础。

PWmat在该工作中的重要作用

整个研究工作的计算部分均通过PWmat得以实现，作为技术基石，起到了从静态电子结构计算到实时量子动力学模拟的全流程支撑作用，提供了一个完整的、经过验证的模拟环境，使得研究者能够将“有限电子温度下电子阻止本领”这一复杂的物理问题转化为可执行、可计算、可分析的数值实验。没有专业的第一性原理计算软件支撑，这项紧密结合了量子理论和动力学过程的深入研究将难以实现。

本研究最核心的部分是PWmat独特的功能含时密度泛函理论（TDDFT）与非绝热分子动力学（NAMD）的结合模拟。当质子在铝晶格中运动时，它会瞬时扰动周围的电子云，PWmat的TDDFT功能模块能够实时求解时间依赖的Kohn-Sham方程，精确描述电子态在质子场作用下的演化，从而计算出质子与电子之间能量转移的瞬时速率。在PWmat的框架下，可以执行Ehrenfest动力学计算。这意味着软件能够自洽地同时处理两件事：一是根据当前离子（质子和铝核）的位置，计算电子的量子态；二是根据电子产生的势场，通过经典牛顿定律推动离子运动。这种“量子-经典”耦合是模拟能量损失（阻止本领）动态过程的唯一途径。

本研究涉及的是处于有限温度下的金属体系，其电子状态非常复杂。PWmat内置了对有限温度电子占据的支持，能够直接使用费米狄拉克分布来处理不同电子温度下的电子占据数。这是模拟“热”电子气、得到不同电子温度下铝的初始电子基态的关键前提。没有这一功能，整个研究的出发点就不复存在。

研究中提到的100 Ry的高截断能确保了平面波展开的充分性，从而保证了电子波函数描述的能量收敛性。PWmat能够高效地处理如此高的截断能计算，并管理复杂的k点采样，这是获得可靠结果的技术保障。TDDFT+MD模拟计算量极其巨大，PWmat具有优秀的大规模并行计算能力，使得在合理的计算资源下对包含数十个原子的超胞进行长时间动力学模拟成为可能。

在动力学模拟过程中，PWmat不仅输出离子的轨迹和能量，还同时记录下电子结构演化的关键信息，这些数据是后续深入分析物理机制的基础。通过分析PWmat输出的波函数等信息，研究人员可以计算投影态密度（PDOS），从而揭示质子和铝电子能级之间的耦合情况，为解释“电荷转移效应”提供了直接证据。与此同时PWmat还可以输出不同时刻的电荷密度，能够实现可视化质子周围电子的局域和捕获现象。