结构搜索的三个阶段提供不同信息
生成、排序和实验确认分别回答“有哪些”“哪个优先”“是否成立”。
读图时保留这项条件：第一名只在已生成集合中领先，正确结构可能根本不在该集合。
加入一项可靠清单，观察候选集合怎样缩小
从三碳八氢一氧的三个连接候选开始，再选择已核实的逐碳类型清单。
8.5.1 把解谱问题写成可检查的对象
计算机辅助结构解析通常简称 CASE。它接收分子式、谱学指认和连接限制，生成符合条件的分子图，再对候选作比较。前一节的全谱拟合从一个既定自旋模型出发调整参数；这里则要回答哪些原子连接方式值得建立自旋模型。两项工作可以相互配合，但一条拟合很好的曲线不会自动补足分子图中尚未确定的键。
展开说明与推导
分子图把原子看成顶点，把化学键看成带键级的边。图中每个碳都需要一个稳定编号，这个编号可以对应峰表里的 C1、C2 等记录，不能随着候选图的绘制位置改变。实验侧还要分别保存信号组、谱线、积分范围和二维交叉峰。一个三重峰含三条谱线，却通常只代表一组化学环境；把三条线建立成三个原子，后续所有约束都会失真。
实际输入最好先经过人工复核。分子式是否来自足够可靠的质量测量，碳峰是否有重合或未检出的季碳，积分是否包含交换氢，编辑 HSQC 的相位是否与实验约定一致，都决定图生成的前提。数据缺口要作为缺口保留。例如两个碳可能共用一个位移值，可写成一个包含两个原子的信号记录，不能为了让计数相等而凭空拆出不同位移。
CASE 也包含不同实现路线。LSD 强调寻找满足输入限制的连接结构；SENECA 的原始工作采用随机搜索和模拟退火；Sherlock 将谱处理、约束编辑、结构生成和排序连在一个界面中。[1–3] 这些名称指向不同算法和工作流。以下的小例子只演示透明的图枚举，不代表已经在这些软件中完成过一轮实测。
候选是满足规定化学空间和约束的一张分子图；指认是实验信号到图中原子或原子组的映射。同一张图可能有多种暂时可行的指认，必须连同映射一起核查。
8.5.2 约束传播怎样减少搜索
先区分必须满足的条件与允许偏差的证据。已确认的分子式、正常价态和可靠的直接连接常可作为硬约束；预测位移与观察位移的接近程度通常适合作为软评分。若把一个经验位移范围设成不可越过的边界，罕见环境可能在生成阶段就被排除，之后再精细排序也找不回来。对自动识别的官能团和数据库片段尤其要保留来源与置信程度。
展开说明与推导
约束传播的意思是把一条已知限制所必然带来的后果继续用于其他变量。一个碳被可靠地确认为 CH3，在通常的中性四价碳模型里，它剩余的键级和只有一；若该位置再被指定同时连接两个重原子，就出现矛盾。一个碳被认作 CH2，则还要容纳两个单位的重原子键级，可以是两个单键，也可以在其他条件允许时是一个双键。碳类型本身还没有给出邻居的身份。
二维数据提供不同强度的限制。HSQC 的可靠单键相关把氢组与承载碳配对；COSY 通常说明氢间存在可观测的标量耦合通路，常用于拼接相邻片段。HMBC 则把氢与较远碳联系起来，常见两键或三键相关，实际体系也可能出现更长程峰或漏峰。软件输入若把每个 HMBC 都锁死成恰好三键，所得“唯一结构”可能只是错误距离假设的结果。
对重叠峰可使用析取约束：该相关来自 Ha 或 Hb，至少一个映射成立。它比任意选一个峰归属更忠实，却会增加搜索分支。记录约束时还应写明实验名称、谱中坐标、处理版本以及人为解释。日后发现一个峰属于杂质，只需撤回它对应的限制，便能重新计算，而不必猜测旧答案为什么发生变化。
未观察到 HMBC 峰，通常不能直接写成“这两个位置之间不存在两键或三键通路”。只有实验灵敏度、耦合选择与峰重叠都经过针对性评估，负证据才可能承担明确的排除作用。
8.5.3 一个能够穷尽的小型枚举
取教学分子式 C3H8O，并规定单一连通、中性闭壳层、碳四价、氧二价的常见结构，不允许非通常价态。其不饱和度为零，所以在这个模型内没有环和多键。四个重原子组成一棵树，需要三条单键。先给三个碳编号，氧保持独立编号；从四点之间六条可能边中选三条，检查连通性、氧的度数以及氢数，最后将只差碳编号置换的图归为同一个连接结构。
展开说明与推导
随书脚本逐一检查所有边集合，保留十五张带碳编号的合法图，经同构去重后得到三种连接方式：丙烷-1-醇、丙烷-2-醇和甲氧基乙烷。这里的“十五”包含重复编号，“三”才是不同构造异构体数。脚本的输出与图 8.5-4 可直接核对；该结论的完备性只覆盖刚刚列出的化学模型，不能推到带电、混合物或其他分子式的任务。
连接图枚举练习包包含仅依赖 Python 标准库的脚本及预期报告。解压后运行脚本，先核对每张图的氢总数，再检查碳编号的六种置换是否归到同一规范表示。更改氧的价态上限前，应先明确更改后的化学含义，否则得到更多图也不代表多发现了合理异构体。
展开说明与推导
现在加入一组理想化教学证据：独立碳类型信息确认一个 CH3、两个 CH2、没有 CH；氢谱可归并为三组非交换氢，积分为 3∶2∶2。B 需要一个 CH，C 需要两个 CH3，均与已确认类型冲突，所以只剩 A。一个 CH3 和两个 CH2 合计七个碳键合氢，分子式中的第八个氢可分配给 OH。即使 OH 因交换未检出，也没有必要伪造一个积分为一的峰。
随后用 HSQC 建立这三组氢各自的承载碳，用 COSY 检查甲基与中间亚甲基、两个亚甲基之间的联系。OH 的 COSY 在可观测时能补充信息，但不能预先要求它一定存在。这个顺序把碳类型、氢计数和耦合网络放在不同证据层次上，读者可以逐条检查候选为何被排除，而无需接受一个不透明的综合得分。
8.5.4 图与峰映射需要共同求解
较大的体系通常无法像上例一样直接完成指认。设有两个重叠氢组，共同落入一个 HMBC 拾峰区域；若先把该区域全部归给某个位置，再枚举结构，程序会继承这个选择。更稳妥的做法是把两个映射都保留，分别检查每个映射与直接连接、积分和其他相关是否相容。输出不仅列出候选名称，还要列出支撑它的映射与尚有歧义的原子组。
展开说明与推导
对称性又带来另一种重复。丙烷-2-醇中两个甲基在通常的非手性环境下等效，交换其编号不产生第二个物理指认。程序应合并这种由图的自同构导致的重复，而不是把两个编号方案当成不同化合物。不过，两个不同位置恰好具有接近的位移属于偶然重叠，不能仅凭数值相等就视为对称等效；它们在二维谱中仍可能给出不同联系。
指认表还应允许一个原子有多个氢位移。非对映异位亚甲基的两氢可能分开出现，但依然由同一个碳承载。若软件把它们误合为两个碳，就会改变图的规模；若误合成一个等效氢组，又会抹去可用于检验候选的耦合信息。
位移预测参与指认时还应避免重复计算证据。若映射已经是通过最小化位移误差选出来的，再把同一组误差当成一份独立的验证结果，会高估支持程度。可以保留部分相关作为后续检查，或比较多种可行映射，并如实说明最优映射的选择过程。NMReDATA 这类记录格式可帮助保存原子编号、指认与谱学参数，[5] 但文件通过格式检查并不等于其中的结构解释已经成立。
8.5.5 搜索覆盖范围与排名各回答什么
候选列表长度同时受化学空间、约束和搜索方法影响。穷尽枚举要求每个允许的连接图都受到检查；随机搜索可能在有限时间里找到若干低评分候选，却不能仅因多次找到同一个答案就宣称没有其他结构。大体系常需要启发式搜索来节省计算。报告应写明算法、终止条件、是否完成搜索，以及是否限制某些键型或强制插入某个片段。
展开说明与推导
排名解决的是已生成候选之间的比较。位移残差、官能团合理性和相关满足程度可以组合为评分，但权重、误差分布与训练资料都会改变次序。一个候选排第一，既不保证真正结构进入了列表，也不自动给出有校准的概率。数据库片段能大幅减少候选数；若片段来自同一套待解释位移，其匹配误差和后续排名之间也可能有关联，不能把两次筛选当成两项独立确认。
若无结构满足全部限制，应先定位互相冲突的最小约束组合。逐项撤回待定指认、可疑长程峰或经验片段，再检查候选变化；不要默默放宽全部条件，只保留程序最终给出的名称。
8.5.6 用候选差异安排下一项实验
有两个以上候选时，先找到它们真正不同的原子联系，再选择能检验该差异的实验。若差别在某个碳是 CH 还是 CH2，可靠的碳类型信息比进一步细调位移评分更直接；若差别是两个片段的连接方向，可针对相应 HMBC 或选择性实验作验证；若连接相同而立体关系不同，则要转向耦合与空间相关。新增数据必须具有候选区分能力，而不是仅增加大家都能解释的峰。
展开说明与推导
最终记录应包括原始谱入口、处理参数、峰表、分子式依据、硬约束、歧义映射、生成限制、候选全集或搜索状态、排序规则，以及后续排除证据。把得到唯一候选的条件写进结论，例如“在给定分子式与正常价态模型内，完整枚举后只剩这一连接结构”。若立体关系未研究，也应在原子编号对应的结构图旁明确保留这个问题。下一节将讨论公共数据库检索，它可以提供候选与对照数据，仍须经过这些检查。
习题与解答
习题 8.5-1 为什么小例子只需三条重原子边？
在规定模型内，C3H8O 的不饱和度为零，因此没有环和多键。四个重原子必须连通，无环连通图是树，边数为顶点数减一，即三。若改变价态或允许带电结构，要重新设定模型，不能继续沿用此证明。
习题 8.5-2 仅有 3∶2∶3 的积分，可否排除甲氧基乙烷？
不能。甲氧基乙烷的碳键合氢本来可以分为两个甲基和一个亚甲基，即 3∶2∶3。还要检查积分是否完整、HSQC 碳类型及耦合网络。若确认两个 CH2 和一个 CH3，才与这个候选的原子计数冲突。
习题 8.5-3 随机搜索十次都找到候选 A，可以称 A 唯一吗？
不能据此证明唯一。重复收敛说明算法在这些初值和设置下倾向找到 A，未证明允许的全部连接图已检查。应报告搜索方式与终止条件，保留其他候选的可能性，或在可行时用完备枚举与区分实验补足证据。
术语
- CASE
- 计算机辅助结构解析，将谱学限制转为候选结构。
- 分子图
- 以原子和键表示连接关系的图对象。
- 硬约束
- 候选必须满足的明确条件。
- 软评分
- 允许偏差并用于候选比较的评价量。
- 约束传播
- 将一个限制的必然后果应用到其他变量。
- 同构去重
- 合并只因顶点编号不同而重复的连接图。
- 析取约束
- 在若干可能解释中至少一种需要成立的限制。
- 峰指认
- 实验信号到原子或原子组的映射。
- 穷尽枚举
- 检查给定化学模型内的全部允许候选。
参考文献
- Nuzillard J M, Massiot G. Logic for Structure Determination. Tetrahedron, 1991, 47: 3655–3664. DOI；作者维护的 LSD 文档。
- Steinbeck C. SENECA: A Platform-Independent, Distributed, and Parallel System for Computer-Assisted Structure Elucidation in Organic Chemistry. J. Chem. Inf. Comput. Sci., 2001, 41: 1500–1507. DOI。
- Wenk M, Nuzillard J M, Steinbeck C. Sherlock—A Free and Open-Source System for the Computer-Assisted Structure Elucidation of Organic Compounds from NMR Data. Molecules, 2023, 28: 1448. 开放原文。
- Han Y, Steinbeck C. Evolutionary-Algorithm-Based Strategy for Computer-Assisted Structure Elucidation. J. Chem. Inf. Comput. Sci., 2004, 44: 489–498. DOI。
- Pupier M, Nuzillard J M, Wist J, et al. NMReDATA, a standard to report the NMR assignment and parameters of organic compounds. Magn. Reson. Chem., 2018, 56: 703–715. DOI。