这次练习的主题编号是 786核心就两个字拆孪晶。更准确地说是把单晶 X 射线衍射数据处理里最常见的两个难点串起来——数据还原和孪晶拆分。很多同学做结构解析练习时指标化不干净、Rint 偏高、结构精修就是不收敛问题往往不在最后的结构解析阶段而是在数据还原阶段就已经把孪晶数据当成单晶数据去处理了。本文按照“练习 786”的目标拆成四步走先判断晶体是不是孪晶再选择正确的还原策略然后完成孪晶拆分最后用拆分后的 HKLF 5 数据做结构精修。整套流程覆盖 Bruker 和 Rigaku 两大常见数据采集体系的处理思路也给出 Olex2 SHELX 精修环节的配置示例。文章适合这几类读者刚开始接触单晶结构解析想系统走一遍数据处理流程的学生手头拿到一组衍射数据指标化不理想怀疑是孪晶但不知道怎么验证已经在用 APEX3、CrysAlisPro 或 Olex2但对 HKLF 5、BASF、孪晶矩阵这些概念还比较模糊。1. 核心能力速览先直观看一下这个练习需要掌握的能力边界。能力项说明项目类型单晶 X 射线衍射数据处理教学练习聚焦数据还原与孪晶拆分数据输入衍射图谱帧文件Bruker .sfrm、Rigaku .img 等数据输出拆分后的 .hkl 强度数据、.ins/.res 精修文件、CIF 文件核心操作指标化、积分、吸收校正、孪晶检测、孪晶拆分、结构精修常用软件APEX3/APEX4 SAINT SADABS、CrysAlisPro、Olex2 SHELXSuite、PLATON是否支持批量支持软件内可批处理多颗晶体数据也可用命令行脚本GPU 需求一般不需要 GPU瓶颈在 CPU 多核利用率与磁盘 IO操作系统Windows 为主SHELX 精修也支持 Linux 环境难度等级中等偏上难点在孪晶判定和 HKLF 5 格式理解不要以为“数据还原”只是点几个按钮真正决定结构解析成败的恰恰是这个阶段有没有把孪晶信息保留下来。2. 数据还原在单晶结构解析中的位置单晶结构解析的完整链路一般是晶体挑选 → 数据收集(衍射实验) → 数据还原 → 结构解析 → 结构精修 → CIF/CCDC 归档其中“数据还原”是承上启下的关键阶段它的任务是把探测器上记录的一堆衍射斑点变成一张表每个衍射点的指标 hkl、强度 I 和标准偏差 sigma(I)。数据还原通常包含五个子步骤指标化Indexing从衍射图中确定晶胞参数和取向矩阵。积分Integration把每一帧衍射图上的衍射斑强度积分出来。吸收校正Absorption Correction校正晶体形状、X 射线路径带来的吸收差异。标度Scaling把不同帧之间的数据统一到同一相对尺度。合并Merging按 Laue 群对称性合并等效衍射点输出最终 hkl 文件。在练习 786 里第 1 步和第 5 步最容易出问题。如果样品本身是孪晶指标化阶段可能拿到两套取向矩阵如果这两套矩阵对应的衍射点在探测器上重叠合并时又不做拆分最后输出的 hkl 文件就是“污染的强度数据”。换句话说孪晶拆分不是精修阶段的补救手段而是数据还原阶段就应该考虑的操作。这也是练习 786 想强调的核心观念。3. 孪晶是怎么产生的如何快速判定孪晶是指一块晶体由两个或多个取向不同的晶域组成这些晶域共享部分晶格或呈对称关联。常见原因包括晶体生长过程中的相变、温度变化、机械应力等。从数据处理角度孪晶主要分三类类型特征对数据处理的影响缺测孪晶merohedral两套晶域倒易格子完全重叠衍射点叠加表观对称性升高单晶指标化也能成功但 Rint 偏高赝缺测孪晶pseudo-merohedral两套格子几乎重叠有微小偏离衍射点部分分离或重叠积分后强度互相污染非缺测孪晶non-merohedral两套格子明显不重合衍射图上有“多余”衍射点单晶指标化时大量点无法定位下面给出判断孪晶的五个信号任何一个命中都应该做进一步检查。3.1 指标化拟合优度差在 APEX3 或 CrysAlisPro 中做单晶指标化后查看拟合优度Fit和未指标化衍射点比例。一般来说未指标化峰占比超过 1% 到 2%就要考虑第二套取向矩阵。3.2 Rint 异常偏高Rint 是等效衍射点合并时的一致度指标。晶胞体积、点群判断正确的情况下Rint 偏高经常是孪晶叠加或点群判定过低的信号。3.3 倒易空间截面出现“双点”在倒易空间查看器中观察 hk 平面可以看到本该是单个的衍射点出现成对排列说明存在两个取向域。3.4 精修时 R1 降不下来结构解析阶段如果流程正确但 R1 始终停在 0.08 以上且差值电子密度有明显正负峰回去检查数据是否是孪晶。3.5 PLATON / XPREP 的 Twin 检测把初步合并好的 hkl 文件放入 PLATON运行 Twin Law 检测或者在 XPREP 中执行 Twin 选项程序会给出可能的孪晶法则。如果确认是孪晶接下来就不是简单重跑一遍而是要在数据还原阶段把两套衍射强度拆开。4. 数据还原阶段的孪晶拆分Bruker 体系实操Bruker 体系的常见软件组合是 APEX3/APEX4 采集、SAINT 积分、SADABS 吸收校正、XPREP 数据转换。下面以这套流程说明孪晶拆分如何操作。4.1 先做单晶指标化再做 Twin 指标化流程上不要上来就切 Twin 模式先执行 Standard 指标化在 APEX3 主界面进入 Index 模块选择自动指标化查看“Solution”列表中的拟合优度如果拟合优度差点击“Twin”按钮让程序搜索第二套取向矩阵确认两套矩阵都能解释可见衍射点保存两个 p4p 文件。这一步的关键是两套矩阵必须物理上可解释。不要盲目接受程序给出的所有候选矩阵要结合晶胞参数和显示方向的衍射点分布判断。4.2 双矩阵积分SAINT 积分时在“Crystal”设置中同时导入两个取向矩阵文件Matrix 1: crystal_0.p4p Matrix 2: crystal_1.p4p积分策略建议选择“Global”积分让程序按两套矩阵分别给衍射点分配强度查看每个矩阵对应的积分统计如果某个衍射点同时被两个矩阵解释需要看程序输出的归属信息。积分得到的 hkl 文件仍然包含两套晶域的贡献需要通过吸收校正和 XPREP 进一步拆分。4.3 SADABS 吸收校正对双矩阵积分结果执行 SADABS。这个步骤不会自己把孪晶拆开但能校正吸收效应并输出更干净的强度统计。SADABS 运行完成后观察 Rint 和合并统计。如果 Rint 仍不理想说明需要 XPREP 做孪晶判定并生成 HKLF 5 数据。4.4 XPREP 输出 HKLF 5 数据XPREP 是 Bruker 体系数据还原和格式转换的中心环节读入 SAINT 输出的 .raw 文件执行“Twin”检测得到可能的孪晶法则在输出选项中选择“HKLF 5”格式程序会生成一个包含两套孪晶域的 hkl 文件并同步写出对应 .ins 文件样板。这里需要特别注意XPREP 判定孪晶法则时会尝试给出一个 3×3 变换矩阵。这个矩阵一定要结合晶胞参数和倒易点验证不能直接默认接受一个表面拟合好但物理意义不明确的矩阵。5. 数据还原阶段的孪晶拆分CrysAlisPro 体系实操Rigaku 仪器配套的 CrysAlisPro 处理思路类似但操作入口不同。5.1 自动指标化与未指标化峰检查打开 CrysAlisPro 后先运行 AutoChemist 或手动执行 Index。查看“Unindexed Peaks”数量如果未指标化峰比例偏高说明存在第二套晶域在数据列表中选中这些峰尝试“Add Crystal”或“Add Domain”程序会再次指标化给出第二套取向关系。5.2 双晶域积分CrysAlisPro 支持多晶域积分。在“Crystal”列表中确认有两个晶域Domain 0 和 Domain 1然后执行 Integrate。积分时可以设置每个衍射点的归属阈值重叠区域的拆分策略输出格式选择 HKLF 5 或保留双矩阵信息。CrysAlisPro 会输出拆分后的 hkl 文件并附带孪晶法则信息。这个文件可以直接拖进 Olex2 做结构解析。5.3 输出前检查输出前重点看三个统计两套晶域各自的拟合优度未指标化峰比例是否明显下降拆分后 Rint 是否比拆分前改善。如果拆分后 Rint 反而更差可能是指标化矩阵选错也可能是两个晶域的衍射点本身没有明显分离需要重新评估孪晶类型。6. 拆分后的结构解析与 HKLF 5 精修完成数据还原和孪晶拆分后hkl 文件以 HKLF 5 格式进入 SHELX 精修流程。这一步是练习 786 最容易出细节错误的地方。6.1 HKLF 5 与 HKLF 4 的区别普通单晶数据用 HKLF 4 指令。孪晶数据要用 HKLF 5并在指令后给出每个孪晶域的变换矩阵和初始比例。HKLF 5 的典型 .ins 文件示例TITL twin_example_786 CELL 7.50 8.00 9.00 90.00 95.00 90.00 ZERR 2 0.01 0.01 0.01 0.01 0.01 0.01 LATT -1 SFAC C H O UNIT 20 30 10 HKLF 5 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 -1.0 0.000 0.000 1.000 0.000 -1.000 0.000 1.000 0.000 0.000 -1.0 BASF 0.1注意HKLF 5 后面每行前 9 个数字是孪晶域的 3×3 变换矩阵最后一个是该晶域的初始体积分数。这里写的是两域示例实际矩阵必须来自数据还原时的指标化结果不能照抄。BASF 是孪晶比例参数初始值一般取 0.1精修后自动收敛。6.2 TWIN 指令和 HKLF 5 的区别缺测孪晶场景下有时可以直接写 TWIN 指令配合 BASFTWIN BASF 0.1TWIN 指令适合“缺测孪晶”且孪晶操作是一个 180 度旋转、程序能自动识别的情况。HKLF 5 则更通用适合非缺测孪晶和需要手动给出变换矩阵的场景。练习 786 建议优先掌握 HKLF 5因为它覆盖的场景更完整。6.3 Olex2 中加载孪晶数据如果使用 Olex2打开 Olex2使用 File Load选择 XPREP 输出的 .insOlex2 自动读取对应 .hkl在 Work 面板中确认是否识别为 HKLF 5用 SHELXT 直接解结构再用 SHELXL 修。Olex2 中可以在“Info”或“Matrix”面板看到孪晶比例。也可以使用 Tools Twinning 相关功能做进一步检查。6.4 精修结果如何判断成功孪晶数据精修成功后重点看这几个指标指标判定参考BASF 值收敛在 0.05~0.5 之间且精修后误差较小R1一般可接受 0.05~0.08特殊情况可放宽wR20.15 以下较理想孪晶数据可适度放宽GooF接近 1.0差值电子密度无明显大于 0.5 e/ų 的残留峰如果 BASF 跑到接近 0 或 1说明“孪晶比例”实际是单晶可能根本不存在孪晶如果 BASF 不收敛或者精修完还有大量残余峰先回去查变换矩阵是否正确。7. 资源占用与处理时间观察单晶数据还原不是深度学习任务基本不依赖 GPU但硬件配置直接影响效率。练习 786 这类数据通常包含几千帧衍射图每帧几 MB 到几十 MB 不等实际处理时可以从下面几个维度观察性能。7.1 CPU 利用率指标化和积分阶段软件会在部分算法中使用多核。观察 CPU 利用率是否跑满如果利用率低且处理慢可能是软件单线程限制而不是电脑性能不够。7.2 内存占用积分大尺寸 CCD/CMOS 探测器数据时内存占用会明显上升。一般而言16 GB 内存能覆盖大多数教学级数据如果数据集特别大关闭其他大内存软件再跑。7.3 磁盘 IO 与存储数据还原会反复读写原始帧文件SSD 能明显缩短积分时间。同时注意保留原始数据目录不建议在积分过程中移动或删除帧文件。7.4 时间估算方法单颗晶体的数据还原时间可以从单帧积分时间的测试外推先抽 10 帧做积分测试得到每帧平均耗时再乘以总帧数就能粗略估算全量积分时间。这个数字因 CPU 和探测器型号差异很大不要拿别人报的时间做严格参照。8. 常见问题与排查方法问题现象可能原因排查方式解决方案指标化拟合优度差样品是孪晶或晶体质量差查看未指标化峰比例尝试 Twin 指标化重新选晶体收集数据指标化后解出超大晶胞使用了错误取向矩阵或存在双晶查看晶胞体积与衍射点密度换矩阵重新指标化单晶指标化成功但 Rint 高可能缺测孪晶XPREP/PLATON 做孪晶检测改用 HKLF 5 BASF 精修积分时间过长帧数多、探测器大、磁盘慢测试单帧积分耗时换 SSD在软件中降低帧分辨率或范围SADABS 报错晶面文件缺失或形状错误检查 .fac 文件重新生成晶体形状描述XPREP 找不到孪晶法则孪晶比例太低或点群判定错误尝试不同 Laue 群手动输入候选矩阵验证HKLF 5 精修 BASF 不收敛变换矩阵错误或初始值差检查 .ins 中矩阵从数据还原中确认正确的孪晶法则R1 降不下来点群、空间群判断错误检查吸收校正和合并统计回到 XPREP 重新定空间群Olex2 加载 hkl 后衍射点异常.ins 和 .hkl 不匹配检查文件格式确认是否 HKLF 5必要时手动指定CIF 输出时缺孪晶信息精修后未添加孪晶法则说明检查 CIF 的 _diffrn 相关条目手动补充孪晶矩阵和体积分数9. 最佳实践与合规提醒9.1 数据处理最佳实践收集数据时多测几颗晶体优先选择衍射点锐利、无拖尾的晶体指标化后先确认未指标化峰比例再决定是否需要 Twin 处理孪晶数据不要过早做高强度合并保持 hkl 文件中的孪晶分量每次数据处理保留原始帧、积分日志、XPREP 输出的副本精修时给 BASF 合理初值不用 0 或 0.5 以外的极端值完成后至少用 PLATON 做一次空间群检查和结构验证。9.2 学术与合规提醒单晶结构解析练习中涉及晶体数据的管理需要注意实验数据必须来自自己或授权使用的晶体样品不得把他人未发表数据用于非授权练习如果结构将发表或提交到 CCDCCIF 中必须如实报告孪晶信息包括孪晶法则、孪晶比例和 BASF 精修结果不要通过隐藏孪晶、强制合并或“挑数据”来得到表面漂亮的精修结果使用商业软件APEX3、CrysAlisPro 等时遵守软件许可协议教学练习优先使用本校或单位已授权的工具涉及第三方晶体结构数据库时注意引用来源和版权规定。10. 总结与下一步练习 786 最值得花时间的地方不是结构精修那一步而是前期的数据还原阶段。建议练习时先用一套已知孪晶比例的数据跑通“指标化 → 双矩阵积分 → XPREP 输出 HKLF 5 → Olex2 精修”这条主链路再换一套未知数据做判断。最容易踩的坑有三个一是指标化时直接忽略未指标化峰二是在 XPREP 中不检查孪晶法则就默认接受三是把 HKLF 5 的变换矩阵写错或漏掉 BASF。后续可以继续扩展的方向用 PLATON 做孪晶法则的自动验证、处理赝缺测孪晶数据、学习 CrysAlisPro 的批处理脚本以及把多颗孪晶数据的还原流程固化成自己的标准操作模板。建议收藏备用下次拿到复杂衍射数据时按这个流程跑一遍比临时翻手册要快得多。