第8章 · 8.1

8.1 数据格式与可复现记录

Data Formats and Reproducible NMR Records
本页范围与学习目标
区分原始FID、数值谱、峰表与结构指认;理解厂商目录、JCAMP-DX、nmrML和NMReDATA,并检验轴单位、转换与复现记录。
先读图解和核心解释;需要细节时再展开。
先看图 · 阅读入口

一份实验记录有三个不同层次

保存这三个层次,才能分清数据变化与解释变化。

一份实验记录有三个不同层次 步骤 1 原始数据 FID 与采集参数 保留频率轴和复数信息 步骤 2 处理结果 谱与全部处理步骤 能够从输入重算 步骤 3 解释记录 峰表、指认和候选 分别标明证据来源 一份实验记录有三个不同层次 步骤 1 原始数据 FID 与采集参数 保留频率轴和复数信息 步骤 2 处理结果 谱与全部处理步骤 能够从输入重算 步骤 3 解释记录 峰表、指认和候选 分别标明证据来源
图 8.1-1 本页概念关系图。箭头表示本页采用的检查顺序。

读图时保留这项条件:同一 FID 的图片和峰表属于同源信息,不是两次独立实验。

8.1.1 一次实验留下多层数据

计算解析首先需要确定程序读入了什么。仪器采集的复数时间信号、处理后的数值谱、导出的图片、人工峰表和结构指认,记录的是不同层次的信息。它们可以相互关联,却不能任意互换。只收到一张谱图时,无法重新检验接收器数字滤波、相位和原始采集窗口;只收到峰表时,还失去了基线、肩峰及峰间噪声。

展开说明与推导

原始数据通常还需要配套的采集参数,才能解释二进制数组。处理谱则需要频率轴、参照及处理历史。峰表把连续谱压缩成局部位置和强度;信号表进一步把同一多重峰的谱线组织起来;指认把这些观察对应到分子原子。后两个步骤包含人的判断或程序模型,应保留不确定关系。

定义 · 数据与解释

FID 和谱数组保存测量及处理结果;峰、信号和原子指认保存提取与解释。可复现记录需要连接这些层次,使读者能追查一个原子标签对应哪组信号、哪张谱以及哪份原始实验。

展开说明与推导

一条“δ 1.20,t,3H”的报告没有列出三条谱线的各自位置,也未说明积分尺度和重叠程度。程序若将它当成一个窄单峰,会建立错误的前向模型。把原子编号与信号编号分开,能避免把三个等价氢误读成三个独立共振,或把一个亚甲基的两个非等价氢合并成已证明等价的对象。

从测量到原子指认的可追查数据链复数 FID采集参数数值谱轴与处理记录峰/信号表提取条件原子指认结构与置信度谱图图片便于阅读每个产物引用其来源;保留独立的原始副本图片未保留完整数值和处理能力峰位变化、删除峰与指认修订都应能回到相应实验
图 8.1-2 原创数据层次图。横线表示有记录的处理或解释关系,图片分支主要承担展示。

8.1.2 厂商目录必须作为整体理解

常见 Bruker 实验目录包含一维时间域文件 fid、采集参数文件,以及 pdata 下的处理结果;二维原始数组常使用 ser。Varian/Agilent 数据目录通常也有名为 fid 的数据文件,并用 procpar 保存参数。相同文件名并不代表二进制布局相同,读取器必须按实际格式解释。

展开说明与推导

这些名称是辨识线索,具体目录会随实验类型和软件版本改变。不要只复制一个大二进制文件就认为已经备份实验。数据类型、字节序、实虚交错、块头、点数及维度解释来自格式和元数据;二维还涉及间接维采样与正交检测方式。误读后出现的镜像或条纹可能是导入错误,不能马上解释成化学现象。

nmrglue 官方教程展示了厂商读取与通用字典的区别:读取结果中的数组可供数值处理,字典保留有关参数。转换到另一格式时,必须检查轴和处理所需信息是否仍在。某个读取器列出支持 Bruker 或 Varian,并不保证任何年代、任何脉冲实验都已经实测。

先保持源目录只读,在工作副本中读取。记录程序版本、格式识别结果和警告;检查数组长度、复数类型、最大幅度与首尾数据。若文件长度与声明点数不符,先追查填充、块结构和采集状态,不要删除若干点直至谱“看起来正常”。数字滤波修正也应作为明确处理步骤保存。

8.1.3 交换格式各有保存对象

JCAMP-DX 是用于交换谱数据的文本格式族。NMR 规范由 Davies 与 Lampen 在1993年发表,后来又有扩展。文本中不仅有数字,也有标签、单位、缩放和压缩规则。读取软件必须正确解码;用普通文本搜索发现几个数字,不能证明已经还原完整数组。

展开说明与推导

一个 JCAMP-DX 文件可能保存处理谱,也可能涉及时间域或多维数据,必须检查该记录的类型。不能把扩展名作为“含原始 FID”的保证。导出者是否保留完整元数据、读入者是否支持相应表示,以及转换后是否改变精度,都是需要检验的环节。

nmrML 用 XML 模式和受控词表组织 NMR 数据及元数据,覆盖原始或处理信息的交换。它有助于把不同厂商的参数映射到共同表达,但不会替代相位校正、自动指认或结构求解。一个符合模式的文件,仍可能含有错写的温度、轴单位或样品标识。

NMReDATA 则着重关联化学结构与提取的位移、耦合、积分和二维相关等数据,采用扩展的 SDF 记录。Pupier 等提出的 NMR record 还把有关谱和原始记录结合起来,便于检查指认依据。它与 nmrML 的重点不同,实际保存可同时包含厂商原始目录、开放数值格式和结构指认文件。

辨析 · 格式通过不等于化学正确

模式校验检查字段与语法;数值校验检查轴、点数及变换;化学校验检查原子计数和相关证据。三者分别回答不同问题。结构图能打开,并不意味着其中的全部指认已被验证。

8.1.4 用一个数组追查轴的含义

设教学实验保存4096个复数采样点,复数采样间隔为0.25 ms。这里的“复数点”各含一个实部和一个虚部;它不等于某些厂商参数中的实数数据字数。先把元数据定义转换到这一明确约定,再计算采样频率和数字频率间隔。

展开说明与推导
\[f_s=1/\Delta t=4000\ {\rm Hz},\qquad \Delta f=f_s/N=0.9765625\ {\rm Hz}.\]
(8.1-1)

采集时间按该离散采样约定约为1.024 s,最后采样时刻为1.02375 s。在400 MHz的氢观察频率下,0.9765625 Hz 对应0.00244140625 ppm。这个数是频率格点间距,不保证能分辨这么近的两条实际谱线;峰宽、噪声和加窗仍然限制分辨能力。

若导出程序把数组补零至8192点,显示格点间距变为0.48828125 Hz,而原始采集仍只有4096点。补零改善数值插值,并未增加观察时间。记录中应分别保存采集点数和变换点数,避免另一位使用者误以为数据是重新采集的高分辨率谱。

要获得绝对位移轴,还需载频、谱宽、FFT 排列及参照。频率正负方向与图上左右方向不是同一个问题;图可以选择高位移在左,数组则可能递增或递减保存。对已知参照峰核验位置,再检查一个非参照信号及谱端,能发现翻转、平移和单位误解。

补零改变显示采样网格而非原始采集原始 N = 4096,复数间隔 0.25 ms,采集约 1.024 s4096 点变换格点间距 0.9765625 Hz补零至 8192 点格点间距 0.48828125 Hz上下横轴表示同一小段频率范围,刻度为网格示意两者都来自相同采集时间,未凭补零增加真实分辨能力
图 8.1-3 原创采集与显示网格图。刻度数量示意补零前后的间隔比例,不代表完整谱的全部点。
习题 8.1-1 · 换场强后的格点

保持上例的时间域采样,在600 MHz下,4096点变换的格点间距是多少 ppm?

参考解答

频率间距仍为4000/4096=0.9765625 Hz;位移间距为0.9765625/600=0.00162760417 ppm。数字采样设置未改变,但同一 Hz 在更高观察频率下对应更小的 ppm。

8.1.5 建立可复算的处理记录

一个最小记录应明确样品标识、实验标识、核种、溶剂、温度、观察频率、采集方式、原始文件路径与校验和。未知浓度或 pH 应标为未知,不填写一个惯用数字。相同样品的重复实验仍使用不同实验编号,才能区分重复制样、重复采集与同一采集的重复处理。

展开说明与推导

处理记录按实际次序保存滤波修正、截取、加窗、补零、傅里叶变换、相位、基线和参照参数。峰提取还需阈值、峰合并条件及积分边界。用图形界面手工调整也可以复现,前提是保存最终参数和操作记录;只保留一张屏幕截图无法还原全部调整。

对教学计算,除了输入与脚本,还要保留依赖版本和输出。设置随机种子的步骤需记录种子与算法;不使用随机过程时无需为了形式而添加种子。结果中写清数值精度,便于区分浮点微小差异与峰轴发生明显变化。

建议把实验观察、提取结果和候选解释分别编号。一个待定亚甲基可保留“两个氢属于同一碳,具体立体指认未定”的信息,而不是任意分配 Ha/Hb 后当成定论。修订指认时保留旧版本和原因,重新处理时则更新关联的谱产物标识。

方法 · 先复算一个小结果

另一环境读入同一原始目录,核对复数点数和采样间隔,复算频率轴及参照峰,再比较一段数值谱和一个积分。小结果一致后,才把全部峰表交给后续模型。

习题 8.1-2 · 文件名相同

两份名为 fid 的文件,一份来自 Bruker,一份来自 Varian。直接用同一二进制读取参数导入有何问题?

参考解答

两者文件布局、块头及配套元数据规则可能不同。应选择对应读取器,保留完整目录并核对复数类型、点数与采样轴。文件名相同只提供线索,不能确定布局。

8.1.6 校验转换与公开记录

转换前后可以先比较点数、单位、轴端点和已知参照位置,再比较峰高、面积与虚部处理。允许误差应由保存精度和转换目的确定,不能看到结果不一致才放宽标准。若目标格式只保存实谱,应明确说明虚谱与原始时间信息没有包含在该产物中。

展开说明与推导

文件校验和用于检查字节是否变化,不能证明内容来自正确样品。压缩包重新打包可能改变压缩文件的校验和,而内部数据字节保持相同;因此可同时保存档案和关键内部文件的校验和。不同处理版本可以有不同校验和,只要关系和原因清楚。

还应检验转换是否保留缺失或待定信息。原始表中的空白可能表示未测量、未检出或无法分离,这些含义不同。若转换软件把空白一律填成零,后续模型可能将其解释为已证明不存在的相关。保存状态字段和文字理由,比用一个看似完整的数字矩阵更可靠。

积分单位也需要单独记录。仪器任意面积、以某信号定为一的相对面积、按氢数归一的面积,以及由标准物换算的浓度,都可以出现在表里。假设某甲基定为3,另一峰面积2,只说明当前响应尺度下相对计数相容;若甲基被重叠污染,两个数字的化学含义一起改变。

二维坐标要标明直接与间接维、核种和单位。一个 HSQC 点的两个数字若顺序交换,会使碳位移进入氢轴。转换后的数组即使形状仍是矩形,也可能把维度排列弄反。以已知的氢碳相关和两个轴范围检查,比单纯确认图能显示更有诊断作用。

采集元数据与化学条件要区分来源。观察频率可来自仪器参数,溶剂或温度也可能由操作者手工填写。自动读取不意味着自动核实了样品状态。记录中可以注明原始参数、人工补充和后续校正,使复核者知道哪些字段可从原目录重建,哪些依赖实验记录。

最后检验一项实际问题:用导出的峰表能否定位到原谱中的区域,并说明积分边界。若峰表的一个峰包含两个重叠信号,保留这一关系;若结构文件重排了原子编号,给出编号映射。没有这些连接,文件数量再多也难以复核具体指认。

nmrXiv 的公开提交流程提供数据与元数据组织的实例。使用平台时还要检查具体记录的访问状态、许可和版本,不把软件仓库许可套到所有上传谱。公开的预测谱、实验处理谱和原始 FID 应保持类型标签;一个可访问链接也不能替代本地可复算的输入说明。

当后续程序给出一个候选结构时,读者应能够追问:输入是原始谱还是峰表,是否移除了溶剂,是否改变积分,结构编号如何映射,哪些相关仍待定。把这些关系保存完整,才能判断错误发生在数据转换、参数提取还是结构模型。下一节将从可下载的原始记录执行一次具体处理。

习题 8.1-3 · 校验和能证明什么

原始文件校验和一致,但两个处理谱的参照峰相差0.20 ppm。可以排除哪些问题,还应检查什么?

参考解答

一致的原始字节支持两次使用了同一文件副本,不能证明处理相同。应检查载频与参照、单位、轴方向和处理参数,也需确认原始文件本身属于目标样品。不能凭校验和判定某一套指认正确。

关键术语

原始数据
采集信号与解释它所需的参数。
数值谱
带轴信息的处理后数组。
峰表
局部谱线位置、强度等的提取记录。
信号表
按共振信号组织的多重峰及积分信息。
指认
观察信号到分子原子的对应关系。
JCAMP-DX
文本形式的谱数据交换格式族。
nmrML
组织 NMR 数据和元数据的 XML 标准。
NMReDATA
关联结构与提取 NMR 参数的记录格式。
校验和
用于检查文件字节变化的摘要。

参考文献与资源

  1. Davies AN, Lampen P. JCAMP-DX for NMR. Applied Spectroscopy, 1993, 47:1093–1099. 原始规范。
  2. Schober D, Jacob D, Wilson M, et al. nmrML: A Community Supported Open Data Standard for the Description, Storage, and Exchange of NMR Data. Analytical Chemistry, 2018, 90:649–656. 论文。
  3. Pupier M, Nuzillard JM, Wist J, et al. NMReDATA, a standard to report the NMR assignment and parameters of organic compounds. Magnetic Resonance in Chemistry, 2018, 56:703–715. 论文。
  4. nmrglue 项目. 官方读取与转换教程;文档支持范围不等于全部实验的实测兼容。
  5. nmrXiv 项目. Publication Deposition;保存数据、元数据和公开记录的操作背景。