pyrolysis温度调控的海枣生物炭/环氧复合材料未碰撞中子输运：基于GEANT4的2 MeV与Cf-252中子研究

《Polymers》：Pyrolysis-Dependent Uncollided Neutron Transmission in Date-Palm Biochar/Epoxy Composites: A GEANT4 Study with 2 MeV and Cf-252 Neutrons

　　

热解改变生物炭的组成与物理结构，因此有必要将化学驱动的中子输运效应与密度和几何效应分离。本计算研究考察了在300、500和800 °C下制备的海枣生物炭作为理想化环氧基质填料的情况。主要的2 MeV对比采用每种材料十次独立的10⁶历史GEANT4运行，条件为20 wt%填料、5 cm厚度和相同体密度（1.25 g cm^?3）。从BC300到BC800，总前向径迹长度注量仅增加0.30%，而注量加权的H*(10)响应增加2.86%。另一项十重复分量研究确认，2 MeV条件下BC800–BC300的H*(10)差异（+0.07986 pSv primary^?1）由未碰撞初级透射增加（+0.08398 pSv primary^?1）驱动，而碰撞初级贡献略有下降（?0.00422 pSv primary^?1）。对于代表性²⁵²Cf Watt谱，H*(10)增加2.98%，其中91.1%的净增加来自未碰撞初级。30 × 30 cm²全面几何增加了散射累积，但BC800–BC300的H*(10)增加中仍有84.0%可归因于未碰撞初级。额外的密度、厚度和负载敏感性在所有测试条件下均给出正的BC800–BC300 H*(10)对比（1.18–8.24%），尽管积分注量对比并非普遍为正。因此，该模型支持高能初级存活率存在稳健的化学相关增加，同时表明散射场分配和工程性能仍依赖于几何与工艺。

**研究背景与问题提出**

生物炭是生物质在限氧条件下经热化学转化生成的富碳固体，作为聚合物复合材料填料可兼顾生物质增值与刚度、热稳定性、电响应、流变性等功能性质调控（^[1-4]）。然而这些性质并非单一"生物炭"材料所固有：原料、热解温度、停留时间、粒径、表面化学、孔结构及后处理均会影响填料及其与聚合物基质的相互作用（^[1,3,4]）。环氧/生物炭体系同样表现出配方依赖性，因此辐射屏蔽评估不应将生物炭视为化学或结构上不变的添加剂（^[5]）。

对中子屏蔽而言，这种变异性尤为重要：快中子慢化强烈依赖氢原子密度与弹性散射概率，而吸收与附加散射则取决于整个配方的元素与同位素组成。富氢聚合物因此被广泛用作中子慢化体，常与含硼相配合以增强慢化后的吸收。聚乙烯/B₄C与环氧/UHMWPE/B₄C体系的实验与计算研究说明了氢含量、吸收体负载、密度与微结构在实际中子屏蔽复合材料中的耦合作用（^[6-8]），也表明屏蔽性能不能仅由填料身份推断。

生物炭近期被视为低成本、生物质来源的屏蔽添加剂。Martellucci与Torsello用PHITS研究生物炭增强混凝土，表明衰减取决于生物炭分数与复合材料密度（^[9），但该工作未分离热解依赖的组成变化如何改变聚合物基质中未碰撞中子与散射中子的平衡。这一区分很重要，因为即使总透射注量变化很小，能谱分布乃至注量-剂量加权量（如H*(10)）也可能显著变化。

海枣残渣提供了受控研究案例：Usman等报道300至800 °C间组成发生显著变化（^[10]），碳随炭化温度升高而增加，氢和氧减少，灰分富集。同时，真实热解还改变颗粒密度、体密度、孔隙率、比表面积、颗粒结构、吸湿性与界面行为（^[3,4,11]）。因此等密度计算可分离化学贡献，但本身无法预测实际制备复合材料的行为。

本研究的创新点并非泛泛证明富氢材料慢化中子，而是量化实测的热解依赖海枣组成变化在环氧基质中如何分配于未碰撞初级、碰撞初级与次级中子贡献，检验该机制在代表性²⁵²Cf裂变谱下是否持续，并评估其对密度、厚度、填料负载及更大全面几何的敏感性。目标应用是在制备辐射防护部件前对聚合物配方进行计算筛选。这些计算仍属模型预测，不能替代实验表征与透射测量。

**主要关键技术方法**

研究采用GEANT4 10.6.2（Shielding参考物理列表、G4NDL4.6中子数据）进行中子输运模拟；填料元素组成取自Usman等报道的300、500、800 °C海枣生物炭（BC300、BC500、BC800），灰分中Ca、Mg、K、Na按氧化物记账处理，余下灰分以SiO₂代理；基质为理想化环氧C₂₁H₂₇O₅，经倒数混合法则得1.25 g cm^?3；采用薄计分体径迹长度注量估计，并以ICRP 74/ICRU 57注量-环境剂量当量转换系数计算H*(10)；设2 MeV单能源与²⁵²Cf Watt谱两种源条件；主比较用十次独立10⁶历史重复，另设分量分类、30 × 30 cm²全面几何及密度/厚度/负载单因素敏感性计算，并以Beer–Lambert解析重建作一致性校验。

**研究结果**

3.1 主要2 MeV比较：十重复实验中，BC800总前向注量仅比BC300高约0.30%，而H*(10)加权响应由2.6957升至2.7729 pSv primary^?1（+2.86%）。BC500–BC300积分注量差异仅+0.078%且95%区间含零，故不支持每个中间热解步骤总注量单调增加。能量分布变化比积分注量更明显：1.9–2.0 MeV份额由BC300的某值升至BC800，差2.578个百分点（95% CI 2.529–2.627）；低于0.1 MeV份额相应下降，差异远大于重复级蒙特卡罗不确定度。

3.2 十重复输运分量分析：独立分量分类实验再现主要排序并揭示其来源。2 MeV下BC300与BC800总H*(10)响应分别为某值与某值（+2.96%），BC800–BC300差异为某值（95% CI 0.07772–0.08200）；未碰撞初级分量增加某值（95% CI 0.08224–0.08571），碰撞初级分量下降某值。即净增加源于未碰撞初级存活增强，碰撞场在原始窄束配置中略作反向抵消。对²⁵²Cf，总H*(10)增加2.98%，其中91.1%来自未碰撞初级、8.0%来自碰撞初级、0.9%来自次级中子，三者Welch 95%区间均不含零。

3.3 评价数据一致性检验：Beer–Lambert重建给出2 MeV宏观总截面EP0、BC300、BC500、BC800分别为0.254214、0.238498、0.231016、0.225753 cm^?1；蒙特卡罗未碰撞透射系统性低于重建值0.43–0.51%，因符号偏移超出抽样不确定度而视为小的系统性实现差异，且不改变材料排序。氢对重建宏观截面的贡献分别约64.6%、62.5%、60.8%、59.8%。

3.4 几何敏感性（30 × 30 cm全面照射）：大几何使俘获散射注量增加，真空归一化未碰撞注量与原始几何几乎相同（BC300 0.30220、BC800 0.32177），而碰撞注量因侧向逃逸减少而升高（BC300 0.5830、BC800 0.5687）。尽管如此，BC800总H*(10)仍比BC300高3.47%，其中84.0%来自未碰撞初级、16.0%来自碰撞初级、<0.1%来自次级。碰撞对比在窄束与全面几何间变号，确认散射场分配依赖几何，但未碰撞初级存活增强在两种情形下均为主要机制。

3.5 密度、厚度与填料负载敏感性：单因素计算显示绝对响应强烈依赖参数（尤其厚度），但BC800–BC300的H*(10)对比在所有条件下均为正：1.15 g cm^?3时+2.53%，1.35 g cm^?3时+3.35%，3 cm时+1.18%，10 cm时+8.24%，10 wt%时+1.51%，30 wt%时+4.19%，六种情形绝对差异的Welch 95%区间均不含零。总注量对比则不稳健：3 cm厚度时BC800总前向注量比BC300低0.295%，但H*(10)仍高1.18%，说明主要热解依赖信号是能量依赖的透射与剂量加权，而非积分注量的普遍单调增加。

**讨论与结论**

讨论指出应区分两个结论：积分前向注量随热解化学仅温和变化且并非在所有敏感性条件下单调；而能谱与H*(10)加权响应呈现更持久的BC300–BC800差异。这与氢数密度由BC300至BC800递减在物理上一致——相互作用概率降低使更高能量初级更易穿透板体，而该部分在注量-剂量转换中权重更大。分量分析直接支持这一解释：原始几何中未碰撞分量大于总差异，因碰撞场部分反向补偿；全面几何中碰撞贡献转正但仍仅占16%，84%净差异来自未碰撞初级。解析G4NDL重建在同一核数据族内独立再现材料排序，支持氢耗竭为一阶主导因素。与聚合物屏蔽文献相比，既往体系多将富氢相与B₄C等吸收体结合以同时实现慢化与吸收，而本研究提出的是不同的筛选问题：在引入专门吸收体之前，以热解依赖的生物质填料替代部分含氢基质将如何改变输运。因此结果宜解释为多功能复合材料设计的组成筛选约束，而非主张以纯生物炭替代常规中子屏蔽填料。局限方面，环氧以理想化有效化学式表示，生物炭按均匀元素混合物建模，残余灰分矿物学以氧化物记账近似，孔隙率、颗粒形貌、水分、填料分散与界面效应均未显式建模；评价数据衰减计算与GEANT4在同一核数据族内吻合约0.5%，不构成独立实验验证，仍需制备并表征的试样在相应源-样-探测器配置下进行中子测量加以确认。

结论部分译文：对理想化海枣生物炭/环氧复合材料开展了组成分辨的GEANT4研究，以确定哪些热解依赖中子输运效应稳健、哪些依赖几何或材料假设。在主要20 wt%、5 cm、等密度比较中，BC800总前向注量仅比BC300高0.30%，但H*(10)加权响应高2.86%。十重复分量计算表明该差异由未碰撞初级中子存活增强驱动；2 MeV下未碰撞贡献增加0.08398 pSv primary^?1，而碰撞分量略降。对代表性²⁵²Cf谱，H*(10)增加2.98%，其中91.1%的净增加来自未碰撞初级。30 × 30 cm²全面敏感性增加散射累积并使碰撞BC800–BC300对比变号，表明散射场依赖几何；即便如此，未碰撞初级仍占总H*(10)增加的84.0%。密度、厚度与填料负载敏感性同样在所有测试条件下保持正的BC800–BC300 H*(10)对比（1.18–8.24%），而总注量差异并非普遍为正。因此稳健结论并非"提高热解温度总提高积分中子透射"，而是氢耗竭增加高能初级存活并可在所建模条件下提高下游剂量加权响应。这些结果为环氧/生物炭配方提供计算筛选框架而非制备处方；真实复合材料性能将取决于实测密度、孔隙率、水分、颗粒形貌、分散、灰分矿物学、源谱与屏蔽几何，在将计算解释为工程屏蔽性能前仍需以制备并表征组成的试样进行实验验证。