中子衍射探测金属中氢同位素_理论方法报告

B: 理论与方法报告——中子衍射探测金属中氢同位素

1. 中子-H/D/T散射长度基础

1.1 散射长度对比表

中子与原子核的散射由束缚散射长度（bound scattering length） 描述。氢及其同位素的中子散射特性具有极显著的差异，这是中子衍射区分氢同位素的物理基础。以下数据源自 Sears (1992) 的经典编纂，以及 International Tables for Crystallography Vol. C (2006) 中的更新版本。

同位素	b_coh (fm)	σ_coh (barn)	b_inc (fm)	σ_inc (barn)	σ_scat (barn)	σ_abs (barn)
¹H (protium)	−3.7406	1.7583	25.274	80.27	82.03	0.3326
²H (deuterium, D)	+6.671	5.592	4.04	2.05	7.64	0.000519
³H (tritium, T)	+4.792	2.89	−1.04	0.14	3.03	0 (有效)
自然氢 (H)	−3.7390	1.7568	—	80.26	82.02	0.3326

数据来源：Sears, V.F. *Neutron scattering lengths and cross sections*, Neutron News, Vol. 3, No. 3 (1992), pp. 26–37；亦见 International Tables for Crystallography, Vol. C, Chapter 4.4 (2006), Table 4.4.4.1。

1.2 关键物理意义

相干散射长度符号差异： ¹H 的 b_coh = −3.74 fm（负值），而 D 为 +6.67 fm（正值），T 为 +4.79 fm（正值）。负号意味着中子波函数在散射时经历 π 相位反转。这一差异直接导致 H 和 D 在全散射信号中的贡献符号相反，这是同位素替换法（Isotopic Substitution） 的核心原理。

非相干散射截面的巨大差异： ¹H 的 σ_inc = 80.27 barn，比 D 的 2.05 barn 高出约 40 倍，比 T 的 0.14 barn 高出近 600 倍。这意味着：

• 含 H 样品的中子衍射图中存在强烈的、几乎各向同性的非相干散射本底
• 用 D 取代 H 可大幅降低本底（~40 倍），同时增强相干信号（~3 倍）
• 这是中子衍射中常规使用氘代样品而非原始氢化物的根本原因

吸收截面的差异： D 的吸收截面 (σ_abs = 0.000519 barn) 远低于 H (0.3326 barn)，T 则几乎无吸收。这使得 D 样品更适用于需要长时间数据采集的高精度衍射实验。

T 的特殊性： T 的 σ_inc 极低 (0.14 barn) 且 b_coh 为正，理论上是最理想的中子衍射同位素，但其放射性（半衰期 12.32 年，β 辐射）和可得性限制了实际应用。对于涉及氚的核材料研究（如聚变堆中的氚滞留），T 的中子散射特性仍是关键参考数据。

1.3 H/D/T 同位素替换法

同位素替换法通过比较同一化合物在 H、D 或混合同位素条件下的中子衍射数据来提取氢原子的位臵信息。其数学基础是：

总相干散射强度 I_coh ∝ |Σ b_coh,i · exp(iQ·r_i)|²

由于 b_coh(H) 与 b_coh(D) 符号相反，在不同 H/D 比率下的测量可以产生不同的对比度，从而将 H 的贡献从总信号中分离出来。

应用实例： 在 LaNi₅H_x 体系中，研究者通过在 LaNi₅D_x 与 LaNi₅H_x 之间进行对比测量，反演出氢原子在晶格间隙位的精准分布。差傅里叶核密度图（Difference Fourier nuclear density map）技术即基于这一原理（详见第 2 节）。

---

2. 间隙位确定：从 Rietveld 到 MEM

2.1 四面体间隙位（T-site）与八面体间隙位（O-site）

在金属晶格中，氢原子通常占据两类间隙位：

• 八面体间隙位（Octahedral site, O-site）： 配位数为 6，位于六个金属原子构成的八面体中心。在 FCC 晶格（如 Pd, LaNi₅ 的 β 相）中，O-site 通常是主要占据位。
• 四面体间隙位（Tetrahedral site, T-site）： 配位数为 4，位于四个金属原子构成的四面体中心。在 BCC 晶格（如 V, Nb, Ta）和某些 FCC 晶格中（存在纳米尺寸效应时）出现。

典型占据模式的实例：

体系	晶格类型	占据位	参考文献
PdH_x (bulk)	FCC	仅 O-site	Akiba et al. (2016) JACS
PdH_x (纳米颗粒, 8 nm)	FCC	30% T-site + 70% O-site	Akiba et al. (2016) JACS
β-VH_0.5	BCT	O₂-site	Otomo et al. (2020) JPSJ Rev.
Ti₃₈V₃₀Cr₁₄Mn₁₈ (BCC相)	BCC	O-site	J. Power Sources (2013)
Ti₃₈V₃₀Cr₁₄Mn₁₈ (BCT/FCC相)	BCT/FCC	T-site	J. Power Sources (2013)
LaNi₅D_x	六方 P6/mmm	多间隙位 (6m, 12n, 4h, 12o, 3f)	Latroche et al. (2004)

纳米尺寸效应的突破性发现： Akiba 课题组（2016）首次通过中子衍射确定了金属纳米颗粒中的氢位臵。在平均粒径 8.0 nm 的 Pd 纳米颗粒中，30% 的 D 原子占据 T-site，70% 占据 O-site。温度依赖测量表明 T-site 占据仅发生在纳米颗粒的亚表面区域。这是首次在金属纳米颗粒中直接测定氢的占据位。来源：Akiba et al., *J. Am. Chem. Soc.* (2016), "Nanometer-Size Effect on Hydrogen Sites in Palladium Lattice."

2.2 Rietveld 精修方法

Rietveld 方法是对粉末中子衍射数据进行全谱拟合的标准方法，通过最小化计算谱与实验谱的差异来精修结构参数。对于氢化物分析，Rietveld 精修的关键参数包括：

• 晶胞参数（反映氢化/氘化引起的晶格膨胀）
• 原子坐标（确定间隙位臵）
• 占据率（Occupancy）（测定各间隙位的氢填充分数）
• 各向同性/各向异性位移参数（ADP）（描述氢的热运动）
• 各向异性展宽参数（描述氢化引起的微应变，如 LaNi₅ 体系）

LaNi₅ 体系的方法学范例： Latroche 等人对 LaNi₅₊ₓ 氘化物进行了系统 Rietveld 精修，使用 FULLPROF 程序。该体系的空间群为 P6/mmm，具有五个显著占据的 D 位（6m, 12n, 4h, 12o, 3f），其中 12n 和 6m 位容纳了约 68% 的 D 原子。关键挑战包括：

• 氘化/退氘循环引入的严重微应变导致各向异性峰展宽，必须使用各向异性峰型函数
• 基线平面展宽参数 S_AA 和 c 轴方向展宽参数 S_CC 需分别精修

来源：Latroche et al., *J. Solid State Chem.* (2004); Percheron-Guégan et al., *J. Alloys Compd.* (1999).

2.3 差傅里叶核密度图（Difference Fourier Nuclear Density Map）

在 Rietveld 精修过程中，若先从模型中省略氢原子，计算结构因子 F_calc，然后与观测结构因子 F_obs 比较，差值 ΔF = F_obs − F_calc 的反傅里叶变换即得到差傅里叶核密度图。氢原子在该图中以其真实位臵的正（D，b>0）或负（H，b<0）核密度峰出现。

该方法与 X 射线衍射存在本质差异：

• X 射线差图由重原子（金届）的电子散射主导，相位角在大范围正确，轻原子如 O 均可被检出
• 中子差图中，Zr 与 O 的中子散射长度可比拟（Zr: b = 7.16 fm; O: b = 5.803 fm），省略原子可能引入大量错误相位角
• H 和 Ti 在中子图中呈现负散射密度，违反直接法对非负密度假设，增加了位臵确定的难度

实例应用： 在 BaSn₀.₅In₀.₅O₂.₇₅₊α 的质子导体研究中，差 MEM 法通过比较 H 溶入样和 D 溶入样的结构因子差，成功分离出氢的同位素分布，揭示 O–D 键倾向于偏转向第二近邻氧原子。来源：Nagasaki et al. (2009), *J. Solid State Chem.*

2.4 最大熵法（Maximum Entropy Method, MEM）

MEM 是一种从有限衍射数据重建核密度分布的非参数方法。其核心思想是：在所有与实验数据一致的密度分布中，选择信息熵最大的一个（即最无偏的分布）。

MEM 在中子晶体学中的优势：

• 可在仅有有限数量衍射数据的条件下产生高分辨率核密度分布
• 不预设原子位臵的高斯分布形状，能捕捉非简谐振动和离域分布
• 特别适用于粉末衍射数据（通常衍射峰重叠，信息量有限）

代表性工作：

1. 氘代甲烷水合物中的氢分布（2006）： Hoshikawa 等人结合 Rietveld 和 MEM，在 7.7–185 K 范围内分析了氘代甲烷水合物的中子粉末衍射数据，清晰揭示了 D₂O 笼子中的低密度 D 位和 CD₄ 分子的取向无序。

1. KDP（磷酸二氢钾）中的氢键非简谐振动： Kumazawa 等人在室温下对 KDP 进行了中子粉末衍射 MEM 分析，发现氢原子沿 O···O 方向拉长，表现出与氢键相关的非简谐热振动特征。

1. BaSn₀.₅In₀.₅O₂.₇₅₊α 中的质子位臵（2009）： 前文所述的差 MEM 法解决了常规 MEM 与 Rietveld 精修之间的不一致性：Rietveld 精修暗示 O–D 键偏向第二近邻氧，而常规 MEM 倾向于 12h 等分位，差 MEM 明确了偏向事实。

1. 中子蛋白质晶体学（2014）： Yamamura 和 Sugawara 将 MEM 应用于分辨率 1.1 Å 的部分氘代蛋白质中子数据，成功观测到大多数 H 和 D 原子，并探索了 MEM 在低分辨率情况下减少全氘代需求的潜力。

2.5 Switendick-Westlake 准则（最小 H–H 距离）

中子衍射实验还验证和修正了一个重要的经验法则：

Switendick-Westlake 准则： 在金届氢化物中，氢原子之间的最小近邻距离约为 2.1 Å。这一距离源自 H 之间的库仑排斥，从根本上限制了金属氢化物的体积储氢密度。

违反该准则的罕见案例（中子衍射直接观测）：

• RNiInH_x 体系（R = Ce, La）中存在 Ni–D···D–Ni–D···D–Ni 链，其 D···D 距离仅约 1.6 Å
• LaNiInD₁.₂₂: D–D = 1.635(8) Å
• CeNiInD₁.₂₄: D–D = 1.606(6) Å

DFT 计算表明，此类短距离通过 Ni 与 H 之间的电子密度局域化（形成 H–Ni–H 共价碎片）屏蔽库仑排斥得以实现。来源：Klein et al. (2022), *J. Alloys Compd.* 905, 162381; 亦见 Switendick (1978), *Z. Phys. Chem. NF*.

---

3. 德拜-沃勒因子（ADP）与氢振动

3.1 各向异性位移参数（ADP）

在中子衍射精修中，原子的热运动通过德拜-沃勒因子描述。对于简谐振动，原子位移用 3×3 张量 U^ij（各向异性位移参数, ADP）表示：

T_k(h) = exp[ −2π² ΣᵢΣⱼ hᵢ aⁱ U^ij aʲ hⱼ ]

中子衍射是研究氢 ADP 的首选方法，因为：

• 中子被原子核散射，而非电子云；H/D 的散射长度与重元素可比
• 因此可以不加约束地精修氢的各向异性 ADP，而 X 射线几乎无法做到
• 金属-氢键键长可达到极高精度（优于 0.01 Å）

来源：Madsen (2012), *Structure and Bonding* Vol. 146, "Modeling and Analysis of Hydrogen Atoms."

3.2 H 与金属原子 ADP 的差异

同位素效应在 ADP 中显著体现：

• 氢原子的振幅通常比金属原子大一个数量级（如 Pd 中 H 的均方根位移 ~0.3–0.5 Å，而 Pd 约 0.07–0.1 Å）
• D 的质量是 H 的两倍，其零点和振动振幅显著小于 H（约 1/√2 倍），这导致在相同温度下 D 的中子衍射峰强度更高、热漫散射本底更低
• 在同一化合物中比较 H 与 D 的 ADP 可获取同位素效应和键合信息

3.3 各向异性与局域振动模式

氢原子的热运动具有高度各向异性：

• 键弯曲运动（bond-bending mode） 的振幅远大于键伸缩运动（bond-stretching mode）
• 平动和转动分子运动进一步加强了各向异性
• 在 Mg₂NiH₄ 中，中子衍射显示氢原子具有高各向异性温度参数，在 Ni 周围的八面体位臵以约 10⁻¹¹ s 的跳变弛豫时间高度局域化运动

3.4 非简谐振动效应

中子 MEM 核密度图可捕捉非简谐振动的特征。KDP 中的 MEM 分析是经典范例：H 核密度沿 O···O 氢键方向显著拉长，不能用简谐高斯模型描述。通过 MEM 重建的有效单粒子势（OPP）可以拟合非简谐参数。

使用 MEM 核密度分布进行非简谐分析的优势：

• 减少参数间的相关性（传统精修中，非简谐参数与高阶多极参数之间存在强耦合）
• 即使对粉末衍射数据也能精确确定势参数
• 可区分静态无序和动态非简谐振动

来源：Kumazawa et al., *J. Phys. Soc. Jpn.*; Nagasaki et al. (2009).

---

4. 非相干散射挑战与处理

4.1 问题来源的定量分析

H 的 σ_inc = 80.27 barn 比 D 的 σ_inc = 2.05 barn 高约 40 倍。这意味着，在典型的粉末中子衍射实验中：

• 含 H 样品：非相干散射贡献的能量宽度覆盖整个探测范围，产生几乎平坦但强度较高的本底（典型值 ~10³ 计数/通道）
• 同一样品经 D 取代后：本底降低约 40 倍，信噪比显著提高
• 因此，对于定量结构分析，常规做法是制备氘代样品

然而，D 取代 H 并非完美方案：

• D 和 H 具有不同的零点能和振动频率（动力学同位素效应），D 占据位可能与 H 存在细微差异
• 在某些体系（如质子导体）中，D 和 H 的扩散路径和占据偏好可能不同
• 研究 H 的原始行为时，需要用方法学手段减去非相干散射本底而非简单地用 D 替代

4.2 非相干散射扣除/分离方法

方法一：空白参比样减法

• 在 SANS 和衍射实验中，用零氢浓度的同种空白样品测量本底
• 将含 H 样品扣除空白本底后得到净散射信号
• 示例：Zircaloy-4 中 δ-氘化物析出相的表征

方法二：非弹性扫描定本底法

• 在中子三轴谱仪上，在固定的 2θ 角度（该角度下无任何相干散射峰）进行能量扫描（如 -1.5 至 +1.5 meV）
• 中心弹性通道的计数即为纯非相干散射贡献
• 用已知 H 含量的聚丙烯标准片校准，建立线性关系 (I = I₀ + Ac) 用于定量
• 示例：Garlea et al. (2007) 在 HFIR HB1A 三轴谱仪上使用 λ = 2.367 Å，2θ = 62°，实现了 Zircaloy-4 中 H 含量的无损测量，检测限约 5–10 μg/g (5–10 ppm)

方法三：极化中子分析（详见第 6 节）

• 纵向极化分析可将自旋翻转（SF，含 2/3 非相干信号）与非自旋翻转（NSF，含全部相干 + 1/3 非相干信号）分离
• 数学公式：σ_inc = (3/2) × I_SF, σ_coh = I_NSF − (1/2) × I_SF
• 示例：Stunault et al. (2016) 在 ILL D3 衍射仪上对液态水实现了"前所未有的相干信号提取精度"

方法四：大 Q 区渐进拟合法

• 在粉末衍射数据中，高 Q 区的非相干本底通常视为缓慢变化的平坦基底
• 通过拟合高 Q 区的基线并扣除，可作为粗略处理

4.3 H 含量定量和成像

利用 H 的高非相干散射截面可以反用于定量检测：

• 在 Zircaloy-4 中，通过测量非相干散射强度校准 H 含量，非线性标定曲线为 I = I₀ + A[1 − exp(−c/c₀)]。低浓度区 (<1000 ppm) 为线性；高浓度区 H 自吸收导致饱和效应。
• 通过小计量体积（1 × 20 × 1 mm³）实现厚度方向空间分辨 mapping，观察到了 δ-ZrH₂ 氢化物层的近表面富集（~400 μm 厚度）。

来源：Garlea et al. (2007), *Mater. Sci. Forum* 539–543, 1443; Choi et al. (2002), *Appl. Phys. A* 76, 797.

---

5. 全散射与 PDF 方法

5.1 方法学基础

中子全散射（Total Neutron Scattering）同时测量相干布拉格散射和漫散射，在绝对标度上得到结构因子 S(Q)。通过对 S(Q) 做傅里叶变换得到对分布函数（Pair Distribution Function, PDF）G(r)：

G(r) = (2/π) ∫ Q[S(Q) − 1] sin(Qr) dQ

G(r) 代表在距离 r 处找到任意一对原子的概率，是"冻结快照"中全部原子间距的加权分布。

PDF 是各偏对分布函数（partial PDF）g_{u-v}(r) 的加权和，每个 g_{u-v}(r) 表示以 u 原子为中心、v 原子在距离 r 处的密度。对于多组分体系（如金属氢化物），需要计算机建模（如逆蒙特卡罗 RMC）来解析重叠的偏 PDF。

来源：Sørby (2008), *Z. Kristallogr.* 223, 617; Sørby (2016), in *Neutron Scattering and Other Nuclear Techniques for Hydrogen in Materials*, Springer.

5.2 与 Rietveld 方法的互补

特征	Rietveld (长程有序)	PDF (局域结构)
信号来源	仅布拉格峰	布拉格 + 漫散射
描述范围	无限周期平均结构	≤ 几十 Å 的局域关联
无序容忍度	有限（需周期模型拟合）	高（可描述无序和缺陷）
对 H 的灵敏度	通过 D 替代较好	同样依赖 D，但可提取 M–H 间距分布

关键互补场景：

• 当体系存在短程有序但长程无序（如非晶氢化物）时，Rietveld 失效而 PDF 奏效
• 当局部结构偏离平均结构（如 H 从理想四面体中心偏移）时，PDF 可捕捉偏移而 Rietveld 只能得到平均位臵
• 原位时间分辨测量可实现充氢/放氢过程中的结构动态追踪

5.3 典型应用

α-VD_0.77（钒氘化物）的短程有序： Sørby 等人使用 RMC 建模全散射数据，约束 D 于理想四面体和八面体间隙位时，确认了 Switendick 准则（D–D ≥ 2 Å）。允许 D 偏移理想位臵时，约 12% 的 D–D 间距短至 1.56 Å——但这些是否为真实物理特征或系统误差仍存争议。

LiAlD₄ 的放氢前兆： 使用 NOVA（J-PARC）和 NOMAD（SNS）的中子全散射 PDF 研究表明，在 40–300 K 范围内，Al–D 最近邻距离几乎不变，但 Li–D 对分布在 200–250 K 以上开始展宽和缩短，这可能标志着 LiD 局域生成作为放氢前兆现象。

V–Ti–Cr 固溶体二氢化物中的 H 偏移： Sakaki 等人（2022, *Acta Materialia*）利用同步辐射 X 射线和中子全散射 PDF，结合特殊准随机结构（SQS）模型和第一性原理计算，发现 H 并不占据四面体间隙的中心，而是向 Cr 原子方向偏移、远离 Ti 原子，形成 H–Cr 键（更短）和 H–Ti 键（更长）。

高熵合金氢化物 TiVCrNbD_x： 中子全散射 RMC 建模显示，在部分脱氢态 (x=2.2) 下，D 以低占据率同时占据四面体和八面体间隙位。占据率与近邻金属原子的价电子浓度之间存在强关联。

5.4 实验装置：NOVA 衍射仪（J-PARC BL21）

NOVA 是专为储氢材料研究设计的高强度中子全散射衍射仪：

• Q 范围：1–40 Å⁻¹，实空间分辨率 0.16 Å
• 配备 Paris-Edinburgh 压机，可实现 20 GPa 以上高压中子衍射
• 配备氢气压-成分-温度 (PCT) 仪表系统：最高 10 MPa H₂，50–473 K
• 样品量小至 ~100 mg 仍可获得可用的结构分析数据
• 时间分辨原位测量可在充/放氢过程中进行瞬态结构追踪

来源：Otomo et al. (2010), *J. Phys. Conf. Ser.* 215, 012024; J-Stage (2012) *HESS* 37, 328.

---

6. 新兴方法前沿

6.1 极化中子衍射（Polarized Neutron Diffraction）

使用极化中子束结合纵向极化分析可在实验上分离相干与非相干散射：

测量通道	包含成分
非自旋翻转 (NSF)	全部相干 + 1/3 非相干
自旋翻转 (SF)	2/3 非相干

→ σ_inc = (3/2) × I_SF→ σ_coh = I_NSF − (1/2) × I_SF

关键进展：

• Stunault 等人（2016）在 ILL D3 极化热中子衍射仪上对氢化液体实现了"前所未有精度的相干信号提取"
• J-PARC BL21 NOVA 衍射仪配备了 ³He 自旋过滤器，结合 Geant4 蒙特卡罗模拟校正多重散射效应
• 该方法的局限：极化分析的强度损失（通常丢失 ~80% 通量），使得实验时间显著延长

动态核极化（DNP）增强的中子晶体学：

• 在亚开尔文温度（~1 K）和高磁场（2.5–5 T）下对样品进行 DNP，使 H 核极化
• 极化后，H 的相干散射截面从 1.8 barn 增强至 14.8 barn（~8 倍），同时非相干本底被完全抑制至零
• 可实现 10–100 倍的性能增益，将所需晶体体积从 ~1 mm³ 降至 ~0.01 mm³
• 已在中子蛋白质晶体学中演示验证

来源：Stunault et al. (2016), *J. Phys. Conf. Ser.* 711, 012003; *Methods in Enzymology* (2020).

6.2 中子全息术（Neutron Holography）

中子全息术是一种新型原子尺度成像技术，可直接重建 H 原子周围的局域原子分布。

原理： H 原子（质子）对热中子有巨大的非相干弹性散射截面（~80 barn）。当 H 嵌入单晶金属样品并在单色中子束中照射时，H 原子作为点状球面中子波源。入射球面波与周围原子散射波的干涉图案在多个样品取向后记录，形成全息图。

两种对偶实现方式：

• 内部源法（Inside-source）： H 作为点源，记录与周围原子的干涉
• 内探测法（Inside-detector）： 强中子吸收性原子核作为点状探测器，入射平面波为参考束，周围原子的散射波为物束

实验演示： 在 PdH_0.78 单晶中成功实现了 Pd 原子相对于 O-site H 的全息重建，直接"看到"了 H 周围的金届原子位臵。

优势： 中子和 X 射线/电子相比，穿透深度大，对轻元素和重元素的灵敏度差异小，特别适用于金属-氢体系。

来源：Cser et al. (2001), *Europhys. Lett.* 54, 747; *J. Alloys Compd.* (2006) "Neutron holography of metal–hydrogen systems."

6.3 逆蒙特卡罗（RMC）建模

RMC 是分析全散射数据的关键工具，尤其适用于无序氢化物：

• 可同时建模布拉格和漫散射信号
• 可包含多种原子种类、多种间隙位类型、原子位移约束
• 适合研究 H 短程聚集和元素偏聚

重要发现（RMC 揭示的聚集信息）：

• Ti₀.₆₃V₀.₂₇Fe₀.₁₀D₁.₇₃ 中的 Fe 形成纳米尺度团簇（~nm 级），这些团簇内几乎不含 D——因为 Fe 本身不形成稳定氢化物——解释了 Fe 替代 V 后容量损失的微观起源
• ZrCr₂D₄ 中，室温（无序态）下的短程有序与低温（完全有序相）结构在 ~4 Å 尺度一致

来源：Sørby (2008), *Z. Kristallogr.*; *J. Alloys Compd.* 448 (2007).

6.4 非弹性中子散射（INS）的补充作用

虽然主题聚焦于衍射，非弹性中子散射与衍射存在重要联系：

• 氢的局域振动模式： INS 直接测量 H 的振动能谱，结合 MEM 核密度和 ADP 可以互相验证
• 德拜-沃勒因子的非弹性基础： DW 因子的温度依赖本质上源于非弹性散射过程的积分
• 在 RNiInH_x 体系中，INS 在 ~90 meV 处观测到 "配对 H 原子" 的特征振动峰，与中子衍射发现的 D···D ~1.6 Å 短距结构完美对应（Klein et al., 2022）

6.5 高压方法

NOVA 的 Paris-Edinburgh 压机（>20 GPa）与中子衍射的组合开辟了极端条件下氢化物结构的探索：

• 镧氢化物（LaH/D）在 GPa 压力下自发电离成两个不同 H 含量的相
• 氢在保持金届亚晶格宏观晶体性的同时输运扩散

---

7. 总结与方法选择流程图

在选择金属氢化物中氢位臵确定方法时，应考虑以下因素：

1. 是否为晶体样品？
2. 是 → Rietveld + 差Fourier / MEM（长程有序）
3. 否（无序/纳米晶/非晶）→ 全散射 PDF + RMC

1. 氢含量足够高（>10 at%）？
2. 是 → 直接 D 替代后衍射
3. 否 → 可能需要极化分析或 DNP 增强

1. 需要区分 H 和 D 吗？
2. 是 → 同位素替换 + 差MEM 或差Fourier

1. H 的局域位臵信息是否足够？
2. 以确定占位为主 → 粉末 Rietveld 即可
3. 需精确键长和振动信息 → 单晶中子衍射 + 非简谐MEM
4. 需可视化 H 周围环境 → 中子全息术

---

8. 参考文献

1. Sears, V.F. (1992). Neutron scattering lengths and cross sections. *Neutron News*, 3(3), 26–37.

1. Sears, V.F. (2006). Scattering lengths for neutrons. In E. Prince (Ed.), *International Tables for Crystallography*, Vol. C, Chapter 4.4, pp. 444–454. IUCr/Wiley.

1. NIST Center for Neutron Research. Neutron scattering lengths and cross sections. https://www.ncnr.nist.gov/resources/n-lengths/

1. Otomo, T. et al. (2020). Structural Studies of Hydrogen Storage Materials with Neutron Diffraction: A Review. *J. Phys. Soc. Jpn.*, 89, 051001.

1. Zhang, T. et al. (2021). Light-weight solid-state hydrogen storage materials characterized by neutron scattering. *J. Alloys Compd.*, 899, 163254.

1. Klein, R.A. et al. (2023). Neutron Scattering Studies of Materials for Hydrogen Storage. In *Comprehensive Inorganic Chemistry III*, Vol. 10, Chapter 10.02, pp. 3–50. Elsevier.

1. Akiba, H. et al. (2016). Nanometer-Size Effect on Hydrogen Sites in Palladium Lattice. *J. Am. Chem. Soc.*, 138, 10238.

1. Latroche, M. et al. (2004). Neutron diffraction study of the deuterides of over-stoichiometric LaNi₅₊ₓ compounds. *J. Solid State Chem.*, 177, 2382.

1. Nagasaki, T. et al. (2009). Neutron powder diffraction and difference maximum entropy method analysis of protium- and deuterium-dissolved BaSn₀.₅In₀.₅O₂.₇₅₊α. *J. Solid State Chem.*, 182, 2632.

1. Hoshikawa, A. et al. (2006). Observation of hydrogen in deuterated methane hydrate by maximum entropy method with neutron powder diffraction. *Physica B*, 385–386, 139.

1. Kumazawa, S. et al. A new analysis technique of anharmonic thermal vibration by using nuclear density distribution of the Maximum Entropy Method. *J. Phys. Soc. Jpn.* (相关出版物).

1. Yamamura, S. & Sugawara, Y. (2014). MEM neutron protein crystallography. *Acta Cryst. A*, 70, a51407.

1. Madsen, A.Ø. (2012). Modeling and Analysis of Hydrogen Atoms. *Structure and Bonding*, 146, 21–52.

1. Sørby, M.H. et al. (2008). Pair Distribution Functions and Reverse Monte Carlo modelling of disordered metal hydrides. *Z. Kristallogr.*, 223, 617.

1. Sørby, M.H. (2016). Total Neutron Scattering. In *Neutron Scattering and Other Nuclear Techniques for Hydrogen in Materials*, Springer.

1. Garlea, V.O. et al. (2007). Neutron Incoherent Scattering Measurements on Hydrogen-Charged Zircaloy-4. *Mater. Sci. Forum*, 539–543, 1443.

1. Choi, Y.N. et al. (2002). Measurement of very small hydrogen content in zirconium alloys by measuring thermal neutron incoherent scattering. *Appl. Phys. A*, 76, 797.

1. Stunault, A. et al. (2016). Structure of hydrogenous liquids: separation of coherent and incoherent cross sections using polarised neutrons. *J. Phys. Conf. Ser.*, 711, 012003.

1. Cser, L. et al. (2001). Atomic-resolution neutron holography. *Europhys. Lett.*, 54, 747.

1. Sakaki, K. et al. (2022). Displacement of hydrogen position in di-hydride of V–Ti–Cr solid solution alloys. *Acta Materialia*, 226, 117620.

1. Sørby, M.H. et al. (2007). Reverse Monte Carlo investigation of the short-range order in α-VD₀.₈. *J. Alloys Compd.*, 446, 267.

1. Sørby, M.H. et al. (2007). Reverse Monte Carlo investigation of temperature dependent deuterium distribution in ZrCr₂D₄. *J. Alloys Compd.*, 448, 249.

1. Klein, R.A. et al. (2022). Neutron vibrational spectroscopic evidence for short H···H contacts in the RNiInH₁.₄;₁.₆ (R = Ce, La) metal hydride. *J. Alloys Compd.*, 905, 162381.

1. Otomo, T. et al. (2010). For high-pressure experiments using total scattering spectrometer NOVA at J-PARC. *J. Phys. Conf. Ser.*, 215, 012024.

1. Switendick, A.C. (1978). Band structure calculations for metal hydrogen systems. *Z. Phys. Chem. NF*, 117, 89.

1. Westlake, D.G. (1983). A model for hydrogen absorption in transition metals. *J. Less-Common Met.*, 91, 1.