Rwp越小越好吗？Rietveld精修的6个影响因素与判据

R_wp 这个数是怎样算出来的？

Rietveld 精修把整条粉末衍射谱当作拟合对象：晶体结构模型给出每个 Bragg 反射的位置和积分强度，峰形函数把这些强度展宽为有一定宽度的峰形，背景函数补上峰之外的散射，三者相加得到每个 2θ 数据点的计算强度 y_i,calc。

最小二乘法反复调整晶胞参数、原子坐标、占位率、位移参数、峰形参数、背景系数和尺度因子，让计算强度与观测强度 y_i,obs 之间的加权残差平方和 Σw_i(y_i,obs − y_i,calc)² 降到最小。

R_wp 就是把这个被最小化的量按观测强度归一化后再开方：R_wp = [Σw_i(y_i,obs − y_i,calc)² / Σw_iy_i,obs²]^1/2。

权重 w_i 通常取计数方差的倒数，计数服从泊松统计时 w_i = 1/y_i,obs，于是分母 Σw_iy_i,obs² 等于整条谱的总计数，分子则是每个点的残差平方除以该点计数后的累加。分子记录模型没有解释掉的强度，分母记录数据一共有多少强度，R_wp 度量的是二者之比。

分母里的观测强度包含背景。峰之外的每一个数据点都把自己的计数加进分母，只要背景函数把这些点拟合得不差，它们对分子就几乎没有贡献；同一份残差放在高背景数据上，R_wp 会比放在低背景数据上小得多。

与之对应的 R_p 不加权重，直接用强度差绝对值之和除以强度之和，强峰在其中占的份额比在 R_wp 里更高；R_B 则不逐点比较，只用各反射的积分强度计算 Σ|I_obs − I_calc|/ΣI_obs，重叠峰的观测积分强度按计算强度比例分配，背景和峰形都不进入这个比值。

精修软件还会给出期望值 R_exp = [(N − P)/Σw_iy_i,obs²]^1/2，其中 N 是数据点数，P 是自由参数个数。模型完美、残差只剩计数噪声时，R_wp 的期望值就是 R_exp；两者之比或比值的平方被记为拟合优度 χ²，接近 1 表示残差与统计噪声同量级。

同一份数据的 R_exp 由计数决定，与模型无关，所以 R_wp 与 R_exp 之比才把数据本身的噪声和模型未拟合的残差区分开。

数据条件对 Rwp 的影响

分子和分母都来自数据，所以采集条件会在结构模型不动的情况下改写 R_wp。

凝固后的锌样品在同步辐射平板探测器上给出了一组直接对照：样品不转动时，二维图上只有零星亮斑，一维谱中若干反射的强度偏离随机取向粉末的预期，R_wp 为 13.6%，R_B 高达 38%；

单轴旋转后德拜环变得连续一些，R_B 降到 18%，R_wp 反而升到 17.8%；双轴旋转让更多晶粒满足衍射条件，R_wp 降到 7.6%，R_B 降到 5.4%。

三组数据用的是同一块样品和同一个锌结构模型，变化的只有参与衍射的晶粒数量。晶粒太少时，落在探测器上的只是少数取向满足 Bragg 条件的颗粒，各反射的积分强度不再服从粉末平均，R_B 记录的正是这种积分强度偏差。

单轴旋转使部分反射的强度回到粉末平均值附近，但每张一维谱的总强度和峰形也随之改变，分母和分子同时变化，R_wp 与 R_B 的走向可以不一致。

这种不一致来自两个指标的对象差别：R_B 只比较各反射的积分强度，R_wp 把每个数据点的残差都计入。取向不均匀、颗粒粗大、样品量不足、透射与反射几何的差别，都会在结构模型完全正确时改变 R_wp 的数值；跨仪器、跨样品比较 R_wp 前，这些采集条件必须先对齐。

背景高的数据把大量计数计入分母，却几乎不增加分子中的残差。锰矿石样品加入刚玉内标后的精修谱中，背景约 1000 计数，最强峰约 5000 计数，两组矿石样品的 R_wp 平均值分别为 2.9 和 3.8；这样低的数值来自背景和非晶散射对分母的抬高，而非更精确的结构模型。

Toby 用模拟数据给出了同样的对照：模型偏差与 Bragg 强度完全相同的两组拟合，无背景时 R_wp = 23%、χ² = 2.54，加入显著背景后 R_wp = 3.5%、χ² = 1.31。

玻璃陶瓷、含非晶相的复合材料、Cu Kα 下荧光背景严重的含 Fe 样品，以及毛细管、Kapton 膜和低角空气散射抬高的基线，都属于高背景数据。

α-堇青石玻璃陶瓷的精修谱里，晶相峰叠在起伏的非晶散射之上，R_p 和 R_wp 都低于 10%，其中相当一部分来自非晶基底对分母的贡献。

部分软件另行给出扣除背景后的 R_wp，把分母换成扣背景强度，这个数值通常大得多，也更贴近峰区的拟合情况。

权重 w_i = 1/y_i,obs 时，R_exp = [(N − P)/Σy_i,obs]^1/2，总计数越高，R_exp 越小。短时间扫描或低功率光源得到的数据总计数少，R_exp 处在百分之几到百分之十几，R_wp 再怎么优化也停在这个噪声下限附近；

同步辐射或长时间累积把总计数提高两个数量级后，R_exp 远低于实验室数据，同样的峰形失配此时会高出噪声幅度，χ² 随之升高。数据越好，R_wp 反而可能更高，这正是 Toby 所说的计数统计悖论。

步长和角度范围也进入 N 和总计数。Rowles 对定量物相分析数据的系统考察表明，步长在 0.01–0.04° 2θ 之间时结果稳定，超过 0.08° 后偏差迅速增大，步长接近或超过峰半高宽时峰顶强度损失；最大计数在 2 万到 5 万之间才能同时兼顾 R_wp 和 χ²，含少量相的样品需要提高到 20 万。

高角上限低于 70° 2θ 时，弱峰密集的高角区被裁掉，N 减小，相含量的偏差同时变大。

模型偏差对 Rwp 的贡献

数据条件确定后，R_wp 剩下的部分才归模型。残差平方在强峰附近最大，因为那里的强度和强度差都最大，权重 1/y 只能部分抵消这种放大；一条结构信息丰富但强度很弱的超结构峰，即使完全没有被模型解释，对分子的贡献也只有强峰残差的几十分之一。

哪些模型偏差能改变 R_wp，取决于它们出现在强峰区还是弱峰区。

峰形函数描述每个反射的宽度、形状和不对称性，它的失配出现在每一个强峰的两翼。

α-Al₂O₃ 纳米粉体的精修用对称赝 Voigt 函数时，(104)、(113) 等强峰底部留下 S 形残差；改用带轴向发散不对称修正的赝 Voigt 后，S 形残差被压平，χ² 从 2.6 降到 2.48，R_wp 只从 21.9% 降到 21.5%。

差谱形状的改善清楚可辨，R_wp 的变化却不到半个百分点。

反过来也成立：峰形参数 U、V、W、Lorentz 分数和不对称参数一旦被放开，R_wp 会先于任何结构参数下降。

晶胞参数、零点和样品位移决定峰位，它们的偏差在强峰两侧形成成对的正负残差，对 R_wp 的贡献同样集中在强峰。

原子坐标、占位率和位移参数影响的是各反射的相对强度，它们造成的残差分散在整条谱上，单个结构参数的错误对 R_wp 的贡献常小于一个峰宽参数。

片状或针状晶粒在样品台上定向排列后，某些晶面族的反射强度被抬高，另一些被压低，结构模型本身没有错，R_wp 却会升高。

Ni(OH)₂ 纳米片的精修引入 March 择优取向参数 r₀₀₁ = 0.81 后，R_wp 为 8.08%，拟合优度 S 为 0.99；把 r₀₀₁ 从 1（随机取向）逐步减小到 0.6，(001) 反射的计算强度随之增长，其他反射相对减弱。

占位率、位移参数和尺度因子也改变强度分配，并且相互补偿：占位率降低、位移参数增大和尺度因子缩小都会减弱某些反射，最小二乘可以在三者之间换算，得到几乎相同的 R_wp；反射数有限时，择优取向参数与占位率之间也可能出现这类替代。

R_wp 只报告残差总量，替代关系之下哪一组参数是物理真实的，R_wp 本身无法分辨。

Le Bail 或 Pawley 拟合只用晶胞、峰形和背景，把每个反射的积分强度当作自由变量，得到的 R_wp 是当前数据在该峰形模型下所能达到的下限。

Sm_0.8La_0.2NiC₂ 按两组晶格参数略有差别的同结构相做 Le Bail 拟合，R_wp = 11.7、R_exp = 8.65、χ² = 1.84，这组数字在没有任何原子坐标的情况下就已经出现；在同一数据上再加入结构模型做 Rietveld 精修，R_wp 不会低于这个下限。

结构模型对 R_wp 的贡献等于 Rietveld 值与 Le Bail 值之差。差值很小，意味着强度模型与数据吻合；差值很大，就要在原子坐标、占位、择优取向或漏掉的相之间找原因。

额外增加的自由参数，例如更高阶的背景多项式或各向异性位移参数，会压低 R_wp，但 N − P 同时减小，R_exp 升高，χ² 未必改善；参数数量的惩罚只在 χ² 里出现，R_wp 里没有。

Rwp 之外，精修结果还要对照哪些量？

一份精修结果通常包含 R_wp、R_exp、χ²、R_B 或 R_F、差谱以及一组结构参数，它们各自记录不同的偏差。

R_wp 与 R_exp 之比给出残差与计数噪声的倍数；R_B 由各反射的积分强度算出，对峰形失配不敏感，对占位率、位移参数和择优取向敏感；差谱给出偏差所在的 2θ 区间和符号；键长、键角、占位率和相含量则要与离子半径、化学式、合成配比和其他测量逐项对照。

峰形函数的选择会改变精修得到的原子坐标。α-Al₂O₃ 的两组精修 R_wp 只差 0.4 个百分点，AlO₆ 八面体的两组 Al–O 键长却从 1.990/1.845 Å 变为 1.992/1.843 Å，O–Al–O 键角也在 80°–101° 区间内重新分配。

键长要与离子半径和文献值对照，键价和要接近阳离子的形式价态，位移参数要保持正值并处于合理数量级，占位率之和要满足化学式和元素分析结果。

相含量要与合成配比、内标加入量和 SEM/EDS 的估计一致；晶粒尺寸和微应变要与电镜观察的尺度对应；内标法给出的非晶含量要与 DSC 或密度测量相互印证。

参数超出物理限值，或者两组参数只靠互相补偿来维持低 R_wp，这份结果就应当在固定弱约束变量后重新精修。

χ² 接近 1 表示残差与计数噪声同量级；远大于 1 时，或者模型仍有系统偏差，或者权重低估了数据的不确定度；小于 1 时通常意味着不确定度被高估，或者自由参数多到开始拟合噪声。

Hill 和 Flack 把 Durbin–Watson 的 d 统计量作为序列相关统计量引入 Rietveld 分析：相邻残差互不相关时 d 接近 2，残差连续同号时 d 降低，这种序列相关会让最小二乘给出的标准不确定度偏小，而 R_wp 对它没有反应。

按标准差归一化的差谱 (y_obs − y_calc)/σ 把强峰和弱峰放在同一尺度上。Ti–6Al–4V 合金原位加热的 437 组同步辐射数据逐一精修后，α-Ti、β-Ti 相分数和各相晶格应变随温度连续变化。某一组精修中跳变的参数或突然放大的残差会在序列里直接暴露出来。

四个代表性时刻的 Δ/σ 曲线幅度随相组成改变：室温 α 相的残差被压在 ±0.05 以内，高温 β 相单相区和 α–β 两相共存区的峰区残差都高出两倍以上，对应峰重叠和峰形变化最剧烈的温度段；同一温度下的晶格应变则由 α-Ti 的 a、c 轴和 β-Ti 分别给出，三者随温度的走势与相分数的转折点一一对应。