超越DPA：YBa2Cu3O7中辐射损伤定量描述的原子尺度框架

《Small Science》：Beyond DPA: An Atomistic Framework for a Quantitative Description of Radiation Damage in YBa2Cu3O7

　　

高温铜氧化物超导体中的辐射损伤是其在聚变反应堆、加速器设施等严苛环境中部署所面临的主要技术挑战之一。这类材料具有复杂晶体结构，使得其辐照效应建模成为一项尤具挑战性的任务，而现有损伤模型对此仍显不足。本研究通过耦合分子动力学（MD）与二元碰撞近似（BCA）模拟，构建了一种用于描述YBa₂Cu₃O₇中初级辐射损伤的原子尺度方法，并使两者形成互补。当该方法与来自蒙特卡罗（MC）代码的初级碰撞原子（PKA）能谱结合后，研究结果建立了一个用于辐射损伤多尺度建模的框架，从而能够在任意以碰撞级联为主导的特定辐照条件下，对多种损伤描述符进行定量估算，例如缺陷产生、缺陷团簇化以及有效受损体积。该计算方法适用于预测任意复杂功能氧化物中的辐照效应，其应用范围涵盖航空航天、核聚变和高能物理等领域。

本文发表于《Small Science》，围绕高温超导铜氧化物YBa₂Cu₃O_7?δ（YBCO）的辐射损伤定量表征问题，提出了一种超越传统NRT-dpa标量剂量指标的原子尺度多尺度计算框架。研究背景在于，YBCO已成为紧凑型聚变堆和下一代高能加速器高场磁体的重要候选材料，但其服役环境中存在强辐射和极端条件，入射粒子与晶格相互作用后会诱发丰富的点缺陷、反位缺陷、缺陷团簇及局域无序区，进而显著改变临界电流密度J_c和超导临界温度T_c。现有研究通常使用Norgett–Robinson–Torrens displacement-per-atom（NRT-dpa）作为不同辐照实验之间的统一剂量尺度，但这一模型最初是为单质金属建立，应用于具有混合共价键—离子键特征的复杂氧化物时，仅能近似反映沉积到晶格中的能量，而难以刻画缺陷形貌、缺陷空间分布及团簇化程度。论文指出，在相似dpa条件下，电子、质子和裂变中子辐照的YBCO单晶表现出显著不同的T_c下降幅度，说明仅凭dpa无法解释超导功能退化的差异，亟需建立面向复杂功能材料的、更具物理基础的损伤描述体系。

研究人员围绕这一问题，构建了一个以初级碰撞原子（PKA）为核心的层级耦合框架，用以覆盖从数eV到数MeV的完整反冲能区。研究的核心思想是：低能PKA主导的局域级联可通过分子动力学（MD）完整解析，高能PKA引发的大级联则先通过二元碰撞近似（BCA）求得弹道骨架，再把能量降至MD适用窗口内的各个子反冲视作独立子级联，并调用预先建立的MD数据库完成重构。借助这一策略，研究人员不仅能够统计每一类PKA在不同能量下产生的缺陷数目，还能定量给出缺陷团簇数量、小团簇比例以及受损体积分数等多维损伤指标，从而建立起从原子尺度缺陷物理到工程尺度辐照预测的可迁移联系。研究结论表明，与NRT-dpa所暗示的“相同剂量对应相近损伤程度”不同，不同粒子类型即便对应相同名义剂量，也会在YBCO中诱导出截然不同的缺陷景观：电子辐照主要产生高密度、空间孤立的小缺陷团簇，而中子辐照则形成更少但更大、更致密的扩展损伤区。该结果为理解不同辐照条件下超导性能退化差异提供了定量依据，也证明了原子尺度损伤形貌指标比单一dpa值更适合作为复杂氧化物的辐射损伤表征工具。

就主要技术方法而言，研究人员使用LAMMPS开展MD级联模拟，并采用专门针对YBCO辐射损伤开发的半经验原子间势，结合Buckingham势、库仑相互作用、短程Ziegler–Biersack–Littmark（ZBL）排斥势及基于SRIM的电子阻止项；高能弹道阶段则采用CASWIN进行BCA传播。研究覆盖Y、Ba、Cu、O四类PKA，在20 K和300 K下模拟1 eV至7 keV的MD级联，以及最高至3 MeV的BCA高能级联；缺陷识别采用OVITO中的Wigner–Seitz算法，受损体积通过Gmsh计算。电子与中子辐照应用部分则将MC代码得到的PKA能谱与MD–BCA响应函数卷积，实现具体辐照条件下的损伤预测。本文未涉及生物样本或临床队列。

在结果部分，论文首先给出了“2.1 Method Overview”的系统方法学构建。研究人员指出，要建立YBCO辐射损伤的原子尺度描述，必须量化材料对全谱PKA的响应。纯MD虽能忠实描述多体相互作用、缺陷复合和级联终态，但受限于计算规模，难以处理keV以上反冲；BCA则适合高能弹道阶段的二体碰撞传播。基于此，论文提出分层耦合策略：低能PKA完全由MD处理，高能PKA先经BCA求出子反冲分布，再用MD数据库中对应能量的平均子级联性质完成整体级联重构。该方案在原理上利用了高能级联会碎裂为空间分离子级联这一事实，因此兼具原子尺度精度与高计算效率。

在该小节中，研究人员还重点阐明了缺陷团簇与受损体积的定义方式。论文认为，仅统计总缺陷数不足以预测YBCO这类复杂功能材料的超导退化，因为扩展团簇可能充当有效磁通钉扎中心而提升J_c，孤立小缺陷则更可能增强载流子散射并抑制T_c与J_c。因此，团簇化程度必须被显式纳入描述。研究以低温库珀对相干长度ξ(T = 0) = 1.12 nm作为超导相关的物理长度尺度，将彼此间距小于2ξ的缺陷归并为同一团簇。对于2 keV以下的MD级联，由于缺陷分布较为紧凑，受损区可近似为沿PKA轨迹中心分布的球形区域，其半径取缺陷分布的回转半径R_g(E)。对于高能级联，BCA传播到各子反冲能量低于2 keV后停止，再由相应的虚拟球形损伤区替代，实现完整几何重构。聚类分析据此在“虚拟球”层面完成，簇体积再由Gmsh计算。该处理使得模型能够在不追踪全部显式Frenkel对坐标的前提下，保留空间组织信息。

随后，论文通过与全MD大级联的对比验证了该框架的可靠性。研究人员选取Ba和O两种代表性PKA，在40、60、110 keV等高能点进行纯MD基准测试，并与BCA–MD耦合结果比较。结果显示，在缺陷总数、团簇统计及受损体积三个层面，耦合模型均与纯MD结果高度一致，且当BCA与MD的分界阈值高于2 keV后，结果对阈值选取不再敏感，说明2 keV作为过渡能量是合理的。文中同时指出，由于采用的是非晶格BCA，该框架目前不考虑沟道效应（channeling recoils），因此其适用范围主要是碰撞级联主导、而非强沟道传输主导的辐照情形。

在“2.1 Method Overview”的后半部分，论文系统展示了由该框架提取出的新型损伤描述符。首先是每次级联的缺陷产额与总受损体积，它们均随PKA能量增加而增长，但呈次线性标度，反映出高能级联中损伤效率下降。研究人员从核阻止与电子阻止的能量分配角度解释了这一现象：重原子Y、Ba、Cu在较宽能区内维持较高核阻止比例，因此缺陷产额与能量近似准线性增长，形成更致密的大级联；氧原子则在约10 keV附近达到核阻止峰值，至500 keV后核阻止显著衰减，因此其在更高能量下缺陷产额与受损体积趋于平台。其次是团簇统计：Y、Ba、Cu诱导的级联形貌相近，产生较少但更大、更致密的团簇；氧级联则更易碎裂，形成空间稀疏、聚集程度较低的损伤分布，因此在大部分能区内小团簇数目占优。该结果强调，PKA种类不仅影响缺陷数量，也深刻决定损伤空间形态。

在“2.2 Application to Representative Irradiation Cases”中，研究人员将上述原子尺度响应函数与实际辐照条件下的PKA能谱卷积，用于比较相同名义剂量1 mdpa下电子辐照与裂变中子辐照对YBCO造成的损伤差异。这一部分是对NRT-dpa局限性的直接检验。结果表明，虽然两种情形具有相同名义NRT-mdpa，且Frenkel对密度处于同一数量级，但缺陷景观存在显著分化。电子辐照以低能氧反冲为主，因此生成大量小而孤立的缺陷团簇；中子辐照可产生高能PKA，形成包含数百个Frenkel对的扩展致密损伤区。定量上，中子辐照的Frenkel对密度为5.48 × 10²⁰ cm^?3，电子辐照为2.45 × 10²⁰ cm^?3；团簇密度则分别为1.44 × 10¹⁸ cm^?3和3.45 × 10¹⁸ cm^?3；小团簇比例分别为53%和99%；平均每个团簇缺陷数分别为405和7；受损体积分数分别为0.98%和1.22%。这说明，电子辐照对应的是“更多、更小、更分散”的损伤组织，而中子辐照对应的是“更少、更大、更致密”的损伤组织。论文据此指出，即便总缺陷数接近，超导性能退化机制也可能因缺陷形貌和空间关联不同而显著分离。

讨论部分的核心在于：对于YBCO这类复杂功能氧化物，辐射损伤绝不能再被压缩为单一标量剂量。传统NRT-dpa更接近于沉积损伤能量的经验性量度，而不足以表征与功能退化直接相关的微观缺陷景观。本文建立的MD–BCA多尺度框架通过引入缺陷产额、团簇密度、小团簇比例和受损体积分数等互补指标，为不同粒子辐照造成的损伤组织差异提供了统一且具物理含义的定量描述。研究的意义不仅在于解释既有YBCO辐照实验之间看似矛盾的结果，更在于为未来辐照实验设计、材料筛选及复杂氧化物在核聚变、航天和高能物理中的应用评估提供预测工具。论文同时审慎指出，当前框架虽已能细致刻画辐照诱导结构无序，但如何将这些原子尺度缺陷特征进一步与宏观超导性质演化直接建立联系，仍是开放问题；未来需要在受控单晶辐照实验、缺陷—超导行为关联理论，以及更长时间尺度下的缺陷迁移与回复过程模拟方面继续推进。

结论部分可译为：总之，本研究证明了当NRT-dpa等简单剂量指标应用于复杂功能材料时所具有的内在局限，并以YBCO为例，提出了一种耦合分子动力学（MD）与二元碰撞近似（BCA）的多尺度计算框架，用于提供具有物理基础的辐射损伤描述。该方法借助BCA进入高能区，克服了纯MD受限于keV能区和小超胞的固有约束。由此建立的MD–BCA数据集，在原子尺度物理与工程尺度预测之间搭建了可迁移的联系，能够以与常规蒙特卡罗（MC）计算相当的代价，提供远超标准NRT-dpa模型描述能力的结构细节。当该框架应用于真实辐照条件时，可捕捉不同粒子类型诱导的差异性缺陷形貌，并通过缺陷产额、团簇密度和受损体积分数等互补指标对受损景观进行定量表征。这一能力既为既有实验结果提供了强有力的解释工具，也为未来辐照研究设计提供了预测基础。尽管当前框架已能详细呈现辐照诱导无序的原子尺度图景，但如何将这些结构变化与超导性质演化直接关联仍有待解决。推动这一方向的发展，需要进一步扩展实验基础，尤其是开展在严格控制辐照条件下的单晶系统研究，并持续推进将微观缺陷景观与宏观功能行为联系起来的理论工作；与此同时，还需将计算分析扩展至更长时间尺度，以捕捉缺陷迁移、回复过程及辐照微结构的长期演化。