简介这份以JadeMDI Jade5.0为工具的XRD数据分析说明文档面向材料科学、矿物学、化学等专业需要完成物相鉴定与晶粒尺寸、晶胞参数计算的师生和实验人员。文档以编号01N2009534样品的课程作业为例完整演示从原始数据导入、初步物相检索与限定条件检索到判定主晶相CaB6O10·5H2O、扣除背底与Kα2、平滑曲线、全谱或40°前选区拟合再按谢乐公式计算平均晶粒尺寸结果为1888Å约188.8nm并通过标注衍射指数与Calculate Lattice求得晶胞参数与标准PDF卡片对比验证。资源包内仅1个doc文件约2.81MB步骤与配图对应紧密适合作为XRD数据处理的上机参考和作业范本。目前已有706人学习下载。1. 从 01N2009534 这份 XRD 数据文档说起Jade 物相、晶胞参数、晶粒尺寸三件事怎么串起来这份编号 01N2009534 的样品 XRD 数据文档课程作业里要求做三件事物相分析、平均晶粒尺寸、晶胞参数。用 MDI Jade 5.0 打开之后界面会把峰位、峰强、半高宽、背底这些信息摆在同一张谱上很多人会把它当成三个独立按钮先检索物相再点 SizeStrain最后 Calculate Lattice。更稳妥的顺序是先定物相再处理背底和 Kα2之后才谈拟合和参数计算。原因很直接物相选错后续标注的衍射指数和晶面间距全部失去意义。这份数据最终锁定主晶相 CaB6O10·5H2O对应 PDF 卡 12-0528平均晶粒尺寸 1888 Å也就是 188.8 nm晶胞参数算到 a10.9769 Å、b16.5742 Å、c6.5742 Åα90°、β91°、γ90°。适合刚接触 Jade 的材料、矿物、化学方向读者也适合需要把课程作业参数讲清楚的人。后续按物相检索、晶粒尺寸、晶胞参数三段拆开。2. Jade 物相检索的判定逻辑初步检索、限定化学元素与 FOM 取舍2.1 数据导入和文本格式自检XRD 文本数据最常见的格式是两列第一列 2θ第二列强度。拖拽到 Jade 之前先确认文件没有表头混杂、分隔符异常、编码错乱。Jade 对文件并不挑剔但它不会替用户判断扫描范围是否合理。常见做法是先用系统命令行看一眼数据头尾避免把强度列当角度列导入。# 查看前 5 行确认两列数据2θ 和强度 head -n 5 01N2009534.txt # 统计行数XRD 谱通常几百到几千点 wc -l 01N2009534.txt # 查看 2θ 范围确认扫描范围是否覆盖 40° 前的强峰 awk NR1{min$1; max$1} {if($1min)min$1; if($1max)max$1} END{print 2θ range:, min, to, max} 01N2009534.txt这段命令的逻辑很直接head检查列结构wc -l确认数据点数量awk找最小和最大 2θ。参数上NR1表示第一行初始化 min 和 max后续每行更新极值。假如输出不是数字或者 2θ 范围只有几度通常说明文件里有表头、空行或分列符异常。Jade 5.0 里也可以走菜单 File/Patterns/Read 导入但拖拽.txt到软件图标上更快适合课程作业里反复打开同一份数据。初步物相检索时右键检索按钮弹出对话框。常用设置如下检索项建议设置原因数据库类型所有类型避免漏掉无机、矿物或有机金属卡片检索目标Major Phase先找主晶相不急着把微量相全部塞进去Use Chemistry不勾选初步检索不限定元素先看候选范围FOM 排序关注但不要迷信FOM 小不等于峰位一定匹配最好最强峰匹配重点检查最强峰对不上主晶相基本可以排除初步检索结果里出现四个可能性较高的物相CaB5O8(OH)B(OH)3(H2O)3、CaB6O10·5H2O、Ca2.62Al9.8Si26.2O72H4.56 和 C20H20N16O8S4Th。前三个属于无机物第四个是有机金属化合物。Jade 的检索逻辑是先按 d 值和峰强给候选再算 FOM。用户要做的不是看谁排第一而是把标准谱和实验谱叠在一起观察最强峰、次强峰和低角区峰位能不能同时对上。2.2 四个候选相的峰位匹配与排除依据物相检索的关键画面是“标准峰竖线”和“样品峰”的对应关系。图 1-3、图 1-4 给出的候选里Ca2.62Al9.8Si26.2O72H4.56 和 C20H20N16O8S4Th 的标准衍射峰没有与样品最强峰匹配尤其是主峰位置偏得明显这类候选可以先放一边。CaB5O8(OH)B(OH)3(H2O)3 和 CaB6O10·5H2O 都能对上不少峰其中 CaB6O10·5H2O 的更强峰对应更整齐所以主晶相优先锁定它。候选物相类型峰位观察初步判断CaB5O8(OH)B(OH)3(H2O)3无机硼酸盐多数峰能对应最强峰有很小偏离可能共存不能单独定为主相CaB6O10·5H2O无机硼酸盐强弱峰对应良好PDF 12-0528主晶相首选Ca2.62Al9.8Si26.2O72H4.56无机硅铝酸盐标准峰未与样品最强峰匹配初步排除C20H20N16O8S4Th有机金属化合物主峰不匹配含 Th 元素初步排除这里有一个很容易踩的坑FOM 值小不一定比 FOM 稍大的卡片更可信。FOM 受峰数、背景、择优取向和卡片质量影响。我一般会把“最强峰是否重合”“低角区前三个强峰是否重合”“高角区峰位是否同步偏移”作为三个硬条件。CaB6O10·5H2O 在这三点上表现更好所以后续晶粒尺寸和晶胞参数都围绕它做。若样品来源不明比如不知道是合成样还是天然矿物就不能只凭元素稀有性排除第四种物相若明确是人工合成考虑 Th 的稀少性以及元素体系差异才有条件把它排除。2.3 限定条件检索与 PDF 12-0528 的确认初步检索给出大范围候选后需要做限定条件检索。操作上重新打开检索对话框勾选 Use Chemistry把元素限定为 Ca、B、O 这类与主晶相一致的体系或者直接在候选卡片列表里勾选 CaB6O10·5H2O 对应的 PDF 卡。此时检索结果会更集中标准峰与样品峰的对应关系也更容易判断。原文中使用的 CaB6O10·5H2O 标准 PDF 卡号为 12-0528后续晶胞参数计算也以这张卡为参照。限定条件检索之后还要保留一点余地样品可能不是单相。CaB5O8(OH)B(OH)3(H2O)3 的峰与样品峰也能对应一部分只是在最强峰附近存在很小偏离。判断主晶相时优先选择强峰匹配更好的 CaB6O10·5H2O同时记录“可能含有 CaB5O8(OH)B(OH)3(H2O)3”这一备选。物相分析不是非黑即白课程作业里写清“一定含有”和“可能含有”比强行给单相结论更可靠。检索完成后回到主界面准备扣背底、扣 Kα2为后面的半高宽拟合做前处理。3. 平均晶粒尺寸计算谢乐公式、Constant FWHM 与 40° 选区拟合3.1 谢乐公式与半高宽补正选项Jade 计算平均晶粒尺寸用的是谢乐公式D Kλ / (β cosθ)其中 D 是晶粒尺寸K 是形状因子λ 是 X 射线波长β 是衍射峰半高宽θ 是布拉格角。Cu Kα 常用波长取 1.5406 Å。β 和 θ 都要用弧度不是度。Jade 内部会做换算但自己核对时千万注意单位。更麻烦的是实验测到的半高宽不只是晶粒细化和应变带来的还包含仪器宽化。没有标准样品做仪器半高宽补正曲线时原文选择 Constant FWHM 作为半高宽补正选项这是课程作业里最常见的折中方案。参数典型取值说明K0.89 或 0.9球形晶粒常用 0.89不同形状会变λ1.5406 ÅCu Kα1 波长β弧度实测 FWHM 减仪器宽化后再代入θ弧度2θ 的一半补正方式Constant FWHM无标样时使用精度有限但可复现用 Python 算一遍可以检查 Jade 的结果有没有明显数量级错误import math K 0.89 wavelength 1.5406 # Å, Cu Kα1 two_theta_deg 29.5 # 示例峰位实际用拟合峰位 fwhm_deg 0.20 # 示例半高宽实际用 Jade 拟合结果 theta math.radians(two_theta_deg / 2.0) beta math.radians(fwhm_deg) D K * wavelength / (beta * math.cos(theta)) print(D , round(D, 2), Å) print(D , round(D / 10.0, 2), nm)这段代码的逻辑是先把 2θ 和 FWHM 从度换成弧度再套谢乐公式。K控制形状因子wavelength按靶材改two_theta_deg和fwhm_deg来自 Jade 的峰拟合结果。若算出来是几纳米到几百纳米数量级通常合理若算出几万纳米多半是 FWHM 单位没换弧度或者把 2θ 当成 θ 用了。3.2 前处理顺序扣背底、Kα2、平滑一次晶粒尺寸计算前原始谱要先做三步扣背底、扣 Kα2、平滑。顺序不要乱。先显示已有背底再执行扣背底和 Kα2最后做一次平滑。平滑次数不宜多一次通常够用次数多了会把弱峰抹平半高宽也会被改变。Jade 主界面里对应的背底按钮和 Kα2 按钮在物相检索之后使用处理后的谱图峰背比会明显改善。扣除 Kα2 是因为 Cu 靶产生的 Kα1 和 Kα2 双线会让峰形不对称尤其在较高角度更明显。半高宽拟合时双线不剥离会把 FWHM 拉宽晶粒尺寸就偏小。平滑是为了压低高频噪声让拟合稳定。若谱图信噪比本来就高也可以跳过平滑直接在 Report/SizeStrain Plot 里看结果。常见做法是保留处理前后两份截图便于课程作业里说明每一步改了什么。3.3 全谱拟合失败后改选区拟合点击全谱拟合时系统提示衍射峰过多意味着整段谱的峰重叠严重自动拟合很难给出可靠半高宽。这时不要硬跑全谱改选区拟合。从这份 XRD 图看主要衍射峰集中在 40° 以前因此选择 40° 区域做拟合。选区拟合的逻辑是只让 Jade 处理峰形相对清楚、重叠较少的区间半高宽更可信。拟合完成后点击 Report/SizeStrain Plot弹出对话框选择 Constant FWHM得到晶粒平均尺寸 1888 Å即 188.8 nm。这个数值来自多个衍射峰的统计结果不是单峰计算。若某个峰的拟合残差大Jade 在平均时仍可能把它算进去所以要看拟合曲线和实验曲线是否贴合。现象可能原因处理全谱拟合提示峰过多峰重叠、背底高改为 40° 前选区拟合FWHM 偏大Kα2 未剥离、背底扣除不足重新扣 Kα2 和背底晶粒尺寸偏小仪器宽化未补正用标准样建立 FWHM 曲线不同峰算出尺寸差很多应变、择优取向、峰重叠分峰后看残差别只看平均值4. 晶胞参数计算衍射指数标注、选区拟合与 Calculate Lattice4.1 只勾选 CaB6O10·5H2O 并完成前处理晶胞参数计算比晶粒尺寸更依赖物相归属。既然主晶相确定为 CaB6O10·5H2O检索结果里只在它前面的方框打钩其他候选相不参与晶胞计算。操作顺序是导入 01N2009534.txt打开物相检索勾选 CaB6O10·5H2O退回主界面显示背底扣除背底与 Kα2再做一次平滑。每一步都会改变峰位和峰形所以原始数据最好备份一份。# 原始数据只读后续任何操作另存项目 cp 01N2009534.txt 01N2009534_backup.txt # 若有多个样品文件批量列出文件名和数据行数 for f in *.txt; do echo $f $(wc -l $f) donecp用来保留原始谱避免 Jade 里反复扣背底后没法回退。for循环列出当前目录下所有.txt文件及行数适合一次处理多个样品时快速检查数据完整性。参数上wc -l $f只输出行数不把文件名混进去方便后续用脚本筛选。4.2 标注衍射指数和选区拟合扣完背底和 Kα2 后点击主界面右下角的衍射指数标注按钮在图上标出主要衍射峰的 hkl。CaB6O10·5H2O 属于单斜晶系β 约 91°峰位分裂和重叠比立方晶系麻烦。标注时优先选 40° 以前强度较高、轮廓清楚的峰。峰太弱或重叠严重的位置不强求全部标出否则容易把相邻峰指错。标注完衍射指数后再做选区拟合。拟合区间仍以 40° 以前为主因为这一段的强峰多、峰形相对独立。拟合结果会给出每个峰的 2θ 和半高宽Jade 再把这些数据送入 Calculate Lattice。若某些峰位偏差超过 0.05°先检查是不是标错 hkl或者 Kα2 没扣干净。4.3 Calculate Lattice 结果与 PDF 卡对比点击 Option/Calculate Lattice直接得到晶胞参数。原文结果为 a10.9769 Åb16.5742 Åc6.5742 Åα90°β91°γ90°。PDF 卡 12-0528 的标准值为 a11.03 Åb16.4 Åc6.577 Åα90°β91.33°γ90°。两者放在一起看a 和 c 很接近b 偏差约 1%β 差 0.33°。参数Jade 计算值PDF 12-0528 标准值相对偏差a / Å10.976911.03-0.48%b / Å16.574216.41.06%c / Å6.57426.577-0.04%α / °90900β / °9191.33-0.33°γ / °90900用 Python 复算偏差比手工看更稳pdf {a: 11.03, b: 16.4, c: 6.577, beta: 91.33} calc {a: 10.9769, b: 16.5742, c: 6.5742, beta: 91.0} for key in [a, b, c]: delta (calc[key] - pdf[key]) / pdf[key] * 100 print(key, 相对偏差:, round(delta, 2), %) print(beta 差值:, round(calc[beta] - pdf[beta], 2), °)代码里pdf和calc分别存标准卡和计算值循环只处理 a、b、c 三个长度参数。beta差值单独输出因为角度差不能用相对百分比表达。若偏差超过 1%优先查峰位标定、样品高度、零偏移和 Kα2 剥离若只有个别峰偏离查择优取向和峰重叠。4.4 用 Origin 画图和误差来源后续画图可以回到原始数据或扣除背底后的数据用 Origin 绘制 CaB6O10·5H2O 的 XRD 图。40° 以前峰位密集部分峰很低图中只标主要衍射峰的衍射指数即可。不要把每一个小包都标成独立峰那样图会失去可读性。晶胞参数偏差来源主要有四类峰位拟合误差、Kα2 残留、样品偏心或高度偏差、单斜晶系峰重叠。b 轴偏差 1% 很可能与低角区峰重叠和拟合区间选择有关。更严格的做法是用 40° 到 80° 的高角峰做晶胞精修但前提是这些峰能被清楚识别。5. 复现时最容易翻车的参数与交叉验证技巧5.1 用标准样建立仪器半高宽曲线Constant FWHM 只是没有标样时的替代方案。真要把 188.8 nm 这个结果做扎实最好先用 Si 或 LaB6 标准样测一条谱提取不同 2θ 下的仪器半高宽拟合成 FWHM-2θ 曲线。Jade 的 SizeStrain Plot 里可以导入仪器宽化曲线再把实测 FWHM 减去仪器宽化。这样做出来的晶粒尺寸通常比 Constant FWHM 更小一点因为仪器贡献被扣掉了。若课程作业只要求流程复现Constant FWHM 够用但报告里要写明这一限制。5.2 多峰拟合与 Williamson-Hall 交叉验证单看平均晶粒尺寸容易被个别坏峰带偏。更稳的做法是把 40° 前几个强峰分别拟合把 2θ、FWHM 列出来代入谢乐公式看离散程度。若不同峰算出的 D 相差两三倍说明应变或峰重叠在作怪这时可以进一步做 Williamson-Hall 分析把尺寸宽化和应变宽化分开。import math K 0.89 wavelength 1.5406 # Å peaks [(20.1, 0.18), (29.5, 0.20), (35.2, 0.22)] # 示例2θ, FWHM for two_theta, fwhm in peaks: theta math.radians(two_theta / 2.0) beta math.radians(fwhm) D K * wavelength / (beta * math.cos(theta)) / 10.0 print(round(two_theta, 2), ° -, round(D, 1), nm)这段代码对每个峰单独算 Dpeaks里换成 Jade 导出的实际拟合值即可。参数顺序是 2θ 和 FWHM单位都是度代码内部转弧度。若三个峰算出来分别是 180 nm、220 nm、170 nm平均值尚可接受若出现 50 nm 和 500 nm先回查峰拟合不要急着解释成晶粒分布。5.3 晶胞参数和物相检索的交叉检查晶胞参数算完以后回头核对物相检索结果。CaB6O10·5H2O 的 a、c 和标准卡很接近b 稍大β 稍小说明主晶相判断基本站得住。假如 Calculate Lattice 给出的晶胞参数与 PDF 卡差到 5% 以上通常不是晶胞真的变了而是物相选错、衍射指数标错或拟合区间混入了杂相峰。复现时保留每一步的截图和拟合残差比只抄最后三个数字更有用。下一次复现时先把标准样的 FWHM-2θ 曲线拟合出来再回到这份 01N2009534 数据188.8 nm 这个数字才有可追溯的误差棒。本文还有配套的精品资源点击获取 MDI Jade 实战：XRD 物相、晶粒尺寸与晶胞参数计算全流程
简介这份以JadeMDI Jade5.0为工具的XRD数据分析说明文档面向材料科学、矿物学、化学等专业需要完成物相鉴定与晶粒尺寸、晶胞参数计算的师生和实验人员。文档以编号01N2009534样品的课程作业为例完整演示从原始数据导入、初步物相检索与限定条件检索到判定主晶相CaB6O10·5H2O、扣除背底与Kα2、平滑曲线、全谱或40°前选区拟合再按谢乐公式计算平均晶粒尺寸结果为1888Å约188.8nm并通过标注衍射指数与Calculate Lattice求得晶胞参数与标准PDF卡片对比验证。资源包内仅1个doc文件约2.81MB步骤与配图对应紧密适合作为XRD数据处理的上机参考和作业范本。目前已有706人学习下载。1. 从 01N2009534 这份 XRD 数据文档说起Jade 物相、晶胞参数、晶粒尺寸三件事怎么串起来这份编号 01N2009534 的样品 XRD 数据文档课程作业里要求做三件事物相分析、平均晶粒尺寸、晶胞参数。用 MDI Jade 5.0 打开之后界面会把峰位、峰强、半高宽、背底这些信息摆在同一张谱上很多人会把它当成三个独立按钮先检索物相再点 SizeStrain最后 Calculate Lattice。更稳妥的顺序是先定物相再处理背底和 Kα2之后才谈拟合和参数计算。原因很直接物相选错后续标注的衍射指数和晶面间距全部失去意义。这份数据最终锁定主晶相 CaB6O10·5H2O对应 PDF 卡 12-0528平均晶粒尺寸 1888 Å也就是 188.8 nm晶胞参数算到 a10.9769 Å、b16.5742 Å、c6.5742 Åα90°、β91°、γ90°。适合刚接触 Jade 的材料、矿物、化学方向读者也适合需要把课程作业参数讲清楚的人。后续按物相检索、晶粒尺寸、晶胞参数三段拆开。2. Jade 物相检索的判定逻辑初步检索、限定化学元素与 FOM 取舍2.1 数据导入和文本格式自检XRD 文本数据最常见的格式是两列第一列 2θ第二列强度。拖拽到 Jade 之前先确认文件没有表头混杂、分隔符异常、编码错乱。Jade 对文件并不挑剔但它不会替用户判断扫描范围是否合理。常见做法是先用系统命令行看一眼数据头尾避免把强度列当角度列导入。# 查看前 5 行确认两列数据2θ 和强度 head -n 5 01N2009534.txt # 统计行数XRD 谱通常几百到几千点 wc -l 01N2009534.txt # 查看 2θ 范围确认扫描范围是否覆盖 40° 前的强峰 awk NR1{min$1; max$1} {if($1min)min$1; if($1max)max$1} END{print 2θ range:, min, to, max} 01N2009534.txt这段命令的逻辑很直接head检查列结构wc -l确认数据点数量awk找最小和最大 2θ。参数上NR1表示第一行初始化 min 和 max后续每行更新极值。假如输出不是数字或者 2θ 范围只有几度通常说明文件里有表头、空行或分列符异常。Jade 5.0 里也可以走菜单 File/Patterns/Read 导入但拖拽.txt到软件图标上更快适合课程作业里反复打开同一份数据。初步物相检索时右键检索按钮弹出对话框。常用设置如下检索项建议设置原因数据库类型所有类型避免漏掉无机、矿物或有机金属卡片检索目标Major Phase先找主晶相不急着把微量相全部塞进去Use Chemistry不勾选初步检索不限定元素先看候选范围FOM 排序关注但不要迷信FOM 小不等于峰位一定匹配最好最强峰匹配重点检查最强峰对不上主晶相基本可以排除初步检索结果里出现四个可能性较高的物相CaB5O8(OH)B(OH)3(H2O)3、CaB6O10·5H2O、Ca2.62Al9.8Si26.2O72H4.56 和 C20H20N16O8S4Th。前三个属于无机物第四个是有机金属化合物。Jade 的检索逻辑是先按 d 值和峰强给候选再算 FOM。用户要做的不是看谁排第一而是把标准谱和实验谱叠在一起观察最强峰、次强峰和低角区峰位能不能同时对上。2.2 四个候选相的峰位匹配与排除依据物相检索的关键画面是“标准峰竖线”和“样品峰”的对应关系。图 1-3、图 1-4 给出的候选里Ca2.62Al9.8Si26.2O72H4.56 和 C20H20N16O8S4Th 的标准衍射峰没有与样品最强峰匹配尤其是主峰位置偏得明显这类候选可以先放一边。CaB5O8(OH)B(OH)3(H2O)3 和 CaB6O10·5H2O 都能对上不少峰其中 CaB6O10·5H2O 的更强峰对应更整齐所以主晶相优先锁定它。候选物相类型峰位观察初步判断CaB5O8(OH)B(OH)3(H2O)3无机硼酸盐多数峰能对应最强峰有很小偏离可能共存不能单独定为主相CaB6O10·5H2O无机硼酸盐强弱峰对应良好PDF 12-0528主晶相首选Ca2.62Al9.8Si26.2O72H4.56无机硅铝酸盐标准峰未与样品最强峰匹配初步排除C20H20N16O8S4Th有机金属化合物主峰不匹配含 Th 元素初步排除这里有一个很容易踩的坑FOM 值小不一定比 FOM 稍大的卡片更可信。FOM 受峰数、背景、择优取向和卡片质量影响。我一般会把“最强峰是否重合”“低角区前三个强峰是否重合”“高角区峰位是否同步偏移”作为三个硬条件。CaB6O10·5H2O 在这三点上表现更好所以后续晶粒尺寸和晶胞参数都围绕它做。若样品来源不明比如不知道是合成样还是天然矿物就不能只凭元素稀有性排除第四种物相若明确是人工合成考虑 Th 的稀少性以及元素体系差异才有条件把它排除。2.3 限定条件检索与 PDF 12-0528 的确认初步检索给出大范围候选后需要做限定条件检索。操作上重新打开检索对话框勾选 Use Chemistry把元素限定为 Ca、B、O 这类与主晶相一致的体系或者直接在候选卡片列表里勾选 CaB6O10·5H2O 对应的 PDF 卡。此时检索结果会更集中标准峰与样品峰的对应关系也更容易判断。原文中使用的 CaB6O10·5H2O 标准 PDF 卡号为 12-0528后续晶胞参数计算也以这张卡为参照。限定条件检索之后还要保留一点余地样品可能不是单相。CaB5O8(OH)B(OH)3(H2O)3 的峰与样品峰也能对应一部分只是在最强峰附近存在很小偏离。判断主晶相时优先选择强峰匹配更好的 CaB6O10·5H2O同时记录“可能含有 CaB5O8(OH)B(OH)3(H2O)3”这一备选。物相分析不是非黑即白课程作业里写清“一定含有”和“可能含有”比强行给单相结论更可靠。检索完成后回到主界面准备扣背底、扣 Kα2为后面的半高宽拟合做前处理。3. 平均晶粒尺寸计算谢乐公式、Constant FWHM 与 40° 选区拟合3.1 谢乐公式与半高宽补正选项Jade 计算平均晶粒尺寸用的是谢乐公式D Kλ / (β cosθ)其中 D 是晶粒尺寸K 是形状因子λ 是 X 射线波长β 是衍射峰半高宽θ 是布拉格角。Cu Kα 常用波长取 1.5406 Å。β 和 θ 都要用弧度不是度。Jade 内部会做换算但自己核对时千万注意单位。更麻烦的是实验测到的半高宽不只是晶粒细化和应变带来的还包含仪器宽化。没有标准样品做仪器半高宽补正曲线时原文选择 Constant FWHM 作为半高宽补正选项这是课程作业里最常见的折中方案。参数典型取值说明K0.89 或 0.9球形晶粒常用 0.89不同形状会变λ1.5406 ÅCu Kα1 波长β弧度实测 FWHM 减仪器宽化后再代入θ弧度2θ 的一半补正方式Constant FWHM无标样时使用精度有限但可复现用 Python 算一遍可以检查 Jade 的结果有没有明显数量级错误import math K 0.89 wavelength 1.5406 # Å, Cu Kα1 two_theta_deg 29.5 # 示例峰位实际用拟合峰位 fwhm_deg 0.20 # 示例半高宽实际用 Jade 拟合结果 theta math.radians(two_theta_deg / 2.0) beta math.radians(fwhm_deg) D K * wavelength / (beta * math.cos(theta)) print(D , round(D, 2), Å) print(D , round(D / 10.0, 2), nm)这段代码的逻辑是先把 2θ 和 FWHM 从度换成弧度再套谢乐公式。K控制形状因子wavelength按靶材改two_theta_deg和fwhm_deg来自 Jade 的峰拟合结果。若算出来是几纳米到几百纳米数量级通常合理若算出几万纳米多半是 FWHM 单位没换弧度或者把 2θ 当成 θ 用了。3.2 前处理顺序扣背底、Kα2、平滑一次晶粒尺寸计算前原始谱要先做三步扣背底、扣 Kα2、平滑。顺序不要乱。先显示已有背底再执行扣背底和 Kα2最后做一次平滑。平滑次数不宜多一次通常够用次数多了会把弱峰抹平半高宽也会被改变。Jade 主界面里对应的背底按钮和 Kα2 按钮在物相检索之后使用处理后的谱图峰背比会明显改善。扣除 Kα2 是因为 Cu 靶产生的 Kα1 和 Kα2 双线会让峰形不对称尤其在较高角度更明显。半高宽拟合时双线不剥离会把 FWHM 拉宽晶粒尺寸就偏小。平滑是为了压低高频噪声让拟合稳定。若谱图信噪比本来就高也可以跳过平滑直接在 Report/SizeStrain Plot 里看结果。常见做法是保留处理前后两份截图便于课程作业里说明每一步改了什么。3.3 全谱拟合失败后改选区拟合点击全谱拟合时系统提示衍射峰过多意味着整段谱的峰重叠严重自动拟合很难给出可靠半高宽。这时不要硬跑全谱改选区拟合。从这份 XRD 图看主要衍射峰集中在 40° 以前因此选择 40° 区域做拟合。选区拟合的逻辑是只让 Jade 处理峰形相对清楚、重叠较少的区间半高宽更可信。拟合完成后点击 Report/SizeStrain Plot弹出对话框选择 Constant FWHM得到晶粒平均尺寸 1888 Å即 188.8 nm。这个数值来自多个衍射峰的统计结果不是单峰计算。若某个峰的拟合残差大Jade 在平均时仍可能把它算进去所以要看拟合曲线和实验曲线是否贴合。现象可能原因处理全谱拟合提示峰过多峰重叠、背底高改为 40° 前选区拟合FWHM 偏大Kα2 未剥离、背底扣除不足重新扣 Kα2 和背底晶粒尺寸偏小仪器宽化未补正用标准样建立 FWHM 曲线不同峰算出尺寸差很多应变、择优取向、峰重叠分峰后看残差别只看平均值4. 晶胞参数计算衍射指数标注、选区拟合与 Calculate Lattice4.1 只勾选 CaB6O10·5H2O 并完成前处理晶胞参数计算比晶粒尺寸更依赖物相归属。既然主晶相确定为 CaB6O10·5H2O检索结果里只在它前面的方框打钩其他候选相不参与晶胞计算。操作顺序是导入 01N2009534.txt打开物相检索勾选 CaB6O10·5H2O退回主界面显示背底扣除背底与 Kα2再做一次平滑。每一步都会改变峰位和峰形所以原始数据最好备份一份。# 原始数据只读后续任何操作另存项目 cp 01N2009534.txt 01N2009534_backup.txt # 若有多个样品文件批量列出文件名和数据行数 for f in *.txt; do echo $f $(wc -l $f) donecp用来保留原始谱避免 Jade 里反复扣背底后没法回退。for循环列出当前目录下所有.txt文件及行数适合一次处理多个样品时快速检查数据完整性。参数上wc -l $f只输出行数不把文件名混进去方便后续用脚本筛选。4.2 标注衍射指数和选区拟合扣完背底和 Kα2 后点击主界面右下角的衍射指数标注按钮在图上标出主要衍射峰的 hkl。CaB6O10·5H2O 属于单斜晶系β 约 91°峰位分裂和重叠比立方晶系麻烦。标注时优先选 40° 以前强度较高、轮廓清楚的峰。峰太弱或重叠严重的位置不强求全部标出否则容易把相邻峰指错。标注完衍射指数后再做选区拟合。拟合区间仍以 40° 以前为主因为这一段的强峰多、峰形相对独立。拟合结果会给出每个峰的 2θ 和半高宽Jade 再把这些数据送入 Calculate Lattice。若某些峰位偏差超过 0.05°先检查是不是标错 hkl或者 Kα2 没扣干净。4.3 Calculate Lattice 结果与 PDF 卡对比点击 Option/Calculate Lattice直接得到晶胞参数。原文结果为 a10.9769 Åb16.5742 Åc6.5742 Åα90°β91°γ90°。PDF 卡 12-0528 的标准值为 a11.03 Åb16.4 Åc6.577 Åα90°β91.33°γ90°。两者放在一起看a 和 c 很接近b 偏差约 1%β 差 0.33°。参数Jade 计算值PDF 12-0528 标准值相对偏差a / Å10.976911.03-0.48%b / Å16.574216.41.06%c / Å6.57426.577-0.04%α / °90900β / °9191.33-0.33°γ / °90900用 Python 复算偏差比手工看更稳pdf {a: 11.03, b: 16.4, c: 6.577, beta: 91.33} calc {a: 10.9769, b: 16.5742, c: 6.5742, beta: 91.0} for key in [a, b, c]: delta (calc[key] - pdf[key]) / pdf[key] * 100 print(key, 相对偏差:, round(delta, 2), %) print(beta 差值:, round(calc[beta] - pdf[beta], 2), °)代码里pdf和calc分别存标准卡和计算值循环只处理 a、b、c 三个长度参数。beta差值单独输出因为角度差不能用相对百分比表达。若偏差超过 1%优先查峰位标定、样品高度、零偏移和 Kα2 剥离若只有个别峰偏离查择优取向和峰重叠。4.4 用 Origin 画图和误差来源后续画图可以回到原始数据或扣除背底后的数据用 Origin 绘制 CaB6O10·5H2O 的 XRD 图。40° 以前峰位密集部分峰很低图中只标主要衍射峰的衍射指数即可。不要把每一个小包都标成独立峰那样图会失去可读性。晶胞参数偏差来源主要有四类峰位拟合误差、Kα2 残留、样品偏心或高度偏差、单斜晶系峰重叠。b 轴偏差 1% 很可能与低角区峰重叠和拟合区间选择有关。更严格的做法是用 40° 到 80° 的高角峰做晶胞精修但前提是这些峰能被清楚识别。5. 复现时最容易翻车的参数与交叉验证技巧5.1 用标准样建立仪器半高宽曲线Constant FWHM 只是没有标样时的替代方案。真要把 188.8 nm 这个结果做扎实最好先用 Si 或 LaB6 标准样测一条谱提取不同 2θ 下的仪器半高宽拟合成 FWHM-2θ 曲线。Jade 的 SizeStrain Plot 里可以导入仪器宽化曲线再把实测 FWHM 减去仪器宽化。这样做出来的晶粒尺寸通常比 Constant FWHM 更小一点因为仪器贡献被扣掉了。若课程作业只要求流程复现Constant FWHM 够用但报告里要写明这一限制。5.2 多峰拟合与 Williamson-Hall 交叉验证单看平均晶粒尺寸容易被个别坏峰带偏。更稳的做法是把 40° 前几个强峰分别拟合把 2θ、FWHM 列出来代入谢乐公式看离散程度。若不同峰算出的 D 相差两三倍说明应变或峰重叠在作怪这时可以进一步做 Williamson-Hall 分析把尺寸宽化和应变宽化分开。import math K 0.89 wavelength 1.5406 # Å peaks [(20.1, 0.18), (29.5, 0.20), (35.2, 0.22)] # 示例2θ, FWHM for two_theta, fwhm in peaks: theta math.radians(two_theta / 2.0) beta math.radians(fwhm) D K * wavelength / (beta * math.cos(theta)) / 10.0 print(round(two_theta, 2), ° -, round(D, 1), nm)这段代码对每个峰单独算 Dpeaks里换成 Jade 导出的实际拟合值即可。参数顺序是 2θ 和 FWHM单位都是度代码内部转弧度。若三个峰算出来分别是 180 nm、220 nm、170 nm平均值尚可接受若出现 50 nm 和 500 nm先回查峰拟合不要急着解释成晶粒分布。5.3 晶胞参数和物相检索的交叉检查晶胞参数算完以后回头核对物相检索结果。CaB6O10·5H2O 的 a、c 和标准卡很接近b 稍大β 稍小说明主晶相判断基本站得住。假如 Calculate Lattice 给出的晶胞参数与 PDF 卡差到 5% 以上通常不是晶胞真的变了而是物相选错、衍射指数标错或拟合区间混入了杂相峰。复现时保留每一步的截图和拟合残差比只抄最后三个数字更有用。下一次复现时先把标准样的 FWHM-2θ 曲线拟合出来再回到这份 01N2009534 数据188.8 nm 这个数字才有可追溯的误差棒。本文还有配套的精品资源点击获取