2.1 什么是Hirshfeld表面？——好比给分子“穿上一件量身定做的衣服”#

想象一下，在一个拥挤的房间里，每个人都占据着一定的空间。如果我们想画出其中一个人所占据的精确区域，最简单的办法是划定一个边界，使得在这个边界内部，这个人的存在感（电子密度）恰好等于周围所有人的存在感之和。这个边界就是Hirshfeld表面。

与传统范德华表面的本质区别在于：范德华表面仅由分子自身原子的大小决定，如同一件按照固定尺码批量生产的成衣；而Hirshfeld表面则充分考虑了晶体的实际堆积环境，如同一件根据您站立时周围人的距离现场裁剪的高级定制礼服。正因为如此，Hirshfeld表面在不同方向上会呈现出不同的凹凸起伏——靠近相邻分子的地方表面会向内凹陷，而远离分子的方向则会向外凸出。

2.2 两个最基本的距离：dᵢ 和 dₑ#

在Hirshfeld表面上的每一个点，我们都可以定义两个距离：

• dᵢ（读作“d inside”）：从该点到表面内部最近原子核的距离
• dₑ（读作“d outside”）：从该点到表面外部最近原子核的距离

您可以这样理解：如果您站在Hirshfeld表面这个“边界”上，dᵢ就是您到分子内部最近邻居的距离，而dₑ则是您到分子外部最近邻居的距离。后续所有指纹图的绘制，本质上就是把晶体中成千上万个表面上点的(dᵢ, dₑ)坐标，像打靶统计一样投影到二维平面上。

2.3 归一化距离 dₙₒᵣₘ——赋予表面色彩的关键指标#

由于不同的原子（如氢、氧、碳）具有不同的范德华半径，直接用dᵢ和dₑ的绝对值来比较不同类型原子间的接触是不公平的。因此，研究者引入了归一化距离dₙₒᵣₘ，其计算公式为：

dₙₒᵣₘ = dᵢ / rᵢᵛᵈʷ + dₑ / rₑᵛᵈʷ

其中rᵢᵛᵈʷ和rₑᵛᵈʷ分别是内部和外部原子的范德华半径。这样处理后，dₙₒᵣₘ的值就具有了统一的物理含义，并通过红-白-蓝三色映射在三维表面上直观呈现：

💡 初学提示：请记住，三维Hirshfeld表面上的红色区域永远是您最先应该关注的焦点，因为它们直接告诉您分子中哪些部位正在与周围分子发生最强的吸引作用。但也要注意，红色区域的面积和深浅需要结合起来看——一个又深又大的红斑往往意味着该处存在贯穿晶体的强氢键网络。

2.4 二维指纹图——一张图看清所有相互作用#

如果说三维Hirshfeld表面告诉您相互作用的空间位置，那么二维指纹图则告诉您每种相互作用的贡献多少。指纹图的横轴是dᵢ，纵轴是dₑ，它将三维表面上所有点的距离信息全部压缩到一张二维散点密度图中。

请想象这样一个场景：如果晶体内全部是松散堆积，分子间没有特定的强相互作用，那么表面上各点的dᵢ和dₑ都会偏大，指纹图呈现为弥散的鹰形分布；如果存在大量O-H···O氢键，则这些接触点会集中在dᵢ和dₑ都很小的左下角区域，形成一对尖峰。正是基于这样的形态差异，经验丰富的研究者仅凭指纹图的整体轮廓，就能快速判断晶体的主要堆积驱动力。

3.1 CrystalExplorer软件的获取与安装#

CrystalExplorer是当前学术界进行Hirshfeld表面分析最主流的软件，由澳大利亚西澳大学的Spackman研究组开发维护。该软件对学术研究免费开放，且支持Windows、macOS和Linux三大平台。

安装步骤：

1. 访问官方网站：https://crystalexplorer.net/
2. 点击页面上的 Download 选项卡
3. 根据您的操作系统选择对应的安装包：
   • Windows用户：选择 .exe 安装包（建议选择64位版本）
   • macOS用户：选择 .dmg 文件
   • Linux用户：选择 .run 或 .AppImage 文件
4. 下载完成后，按照常规方式安装（Windows系统双击exe文件按提示操作即可）

硬件与系统要求提示：

• 操作系统建议Windows 10以上、macOS 10.14以上或主流Linux发行版
• 内存至少4 GB（推荐8 GB以上）
• 硬盘空间至少2 GB
• 若处理较大分子（如含50个以上非氢原子），建议配备独立显卡以加速三维渲染

3.2 输入文件的准备——CIF文件从哪来？#

CrystalExplorer的输入文件格式为 CIF（Crystallographic Information File） ，这是晶体学领域通用的标准格式文件。获取CIF文件的途径主要有三种：

途径一：如果您已经完成了单晶X射线衍射实验

您可以通过晶体解析软件（如Olex2、SHELXTL、CrysAlisPro等）完成结构精修后，在软件中导出CIF格式文件。通常，精修完成后软件会自动生成CIF文件，您只需检查该文件是否包含完整的原子坐标、温度因子和晶胞参数即可。

途径二：从剑桥结构数据库（CSD）中获取已有晶体数据

CSD是世界上最大的小分子晶体结构数据库，由剑桥晶体学数据中心（CCDC）维护。您可以通过以下方式获取：

• 访问CSD网络门户：http://www.ccdc.cam.ac.uk/structures/
• 在搜索框中输入化合物名称、分子式或CCDC编号
• 找到目标结构后，点击下载CIF文件（下载可能需要注册或机构订阅权限）

途径三：使用CrystalExplorer自带的“教学子集”练习数据

如果您只是想练习操作而手头暂时没有自己的CIF文件，CrystalExplorer官网提供了名为“Teaching Subset”的示例数据集，其中包含多个经典化合物的CIF文件，例如：

• ACETAC07：乙酸（展示强氢键尖峰的经典案例）
• BENZEN：苯（展示CH-π相互作用特征）
• HCLBNZ11：六氯苯（展示卤素参与相互作用的案例）
• BARBAD01：巴比妥酸二水合物（展示结晶水参与氢键网络的案例）

💡 建议：初次上手时，强烈建议先使用上述教学子集中的BENZEN和ACETAC07这两个文件进行练习。它们分子小、计算快、特征明显，能帮助您快速建立对不同指纹图特征的直观感觉。

3.3 ⚠️ 最容易被忽视的关键步骤：氢原子位置的标准化处理（新版本CrystalExplorer已经自动化处理）#

这是整个分析流程中最容易导致结果偏差却最常被初学者忽略的环节。常规的X射线衍射数据在测定氢原子位置时存在一个系统性缺陷——X射线主要被电子云散射，而氢原子只有一个电子，其电子云偏向于与之成键的重原子（如碳、氧、氮），导致测定出的氢原子位置会向重原子方向偏移。如果直接使用这些“偏短”的氢原子坐标进行Hirshfeld表面分析，计算出的氢键距离会系统性偏小，导致dₙₒᵣₘ值异常偏低，可能错误地将一些原本较弱的接触误判为强氢键。

正确的处理方法是：在将CIF文件导入CrystalExplorer后，不要急着生成表面，而应当首先进行氢原子位置的标准化修正。具体操作路径为：

在CrystalExplorer主界面上方菜单栏中，依次点击：

Tools → Neutron Normalization

这个功能会将所有X-H键长修正为更准确的中子衍射标准值（如C-H修正至1.083 Å，N-H修正至1.009 Å，O-H修正至0.983 Å）。完成修正后，您才能在后续分析中得出具有化学合理性和统计可比性的结果。

📌 核心要点：无论您使用何种方式获得CIF文件，请在生成Hirshfeld表面之前务必执行氢原子标准化操作。如果您的论文最终要发表，也应当在方法学部分明确说明已进行了此步处理。

4.1 第一步：打开CIF文件并补全分子#

1. 启动CrystalExplorer程序
2. 点击菜单栏 File → Open，在弹出的文件选择对话框中找到您的CIF文件（此处以BENZEN.cif为例）
3. 单击选中后点击 Open 或 打开，晶体结构将呈现在三维图形窗口中

此时您可能会发现，窗口中显示的不一定是完整的苯分子，而可能只是不对称单元中的一部分原子。这是完全正常的，因为CIF文件为了节约存储空间，往往只记录独特部分，其余原子可以通过对称操作生成。

为了得到完整的分子以进行后续分析，请执行以下操作： 4. 在三维窗口中，使用鼠标左键框选或逐个点击选择属于您目标分子的所有原子 5. 在菜单栏中选择 Edit → Complete Molecule（或使用工具栏中的相应图标） 6. 软件将自动根据晶体的对称性补全所选原子的所有对称等价原子，形成一个完整的分子

💡 操作技巧：您也可以通过右键点击任意一个原子，在弹出的菜单中选择 Select Molecule，软件会自动识别并选中与该原子属于同一独立分子的全部原子，然后再执行Complete Molecule操作。这种方法在处理复杂分子时更为高效。

4.2 第二步：生成Hirshfeld表面#

确保已选中完整的目标分子（分子周围应有高亮边框表示被选中），然后执行以下操作：

1. 点击主工具栏中的 Surface Generation 按钮（通常是一个带有球体表面图案的图标）
2. 在弹出的参数设置对话框中，按以下选项配置：
   • Surface Type：选择 Hirshfeld
   • Property：建议首先选择 dₙₒᵣₘ（这是我们最常使用的默认属性）
   • Resolution：选择 High（以获得较精细的表面细节）
   • Wavefunction：选择 New wavefunction（软件将调用内置Tonto程序计算波函数）
3. 点击 OK 确认

此时，软件将开始计算波函数并生成Hirshfeld表面。计算过程中，程序窗口下方会显示进度条。对于苯这样的小分子，通常只需要几秒钟即可完成。计算完成后，您会看到目标分子被一层红-白-蓝三色相间的半透明表面所包裹。

⚠️ 计算时间提醒：如果您处理的是含有20个以上非氢原子或包含过渡金属的分子，计算时间可能会显著延长（从几分钟到十几分钟不等）。建议在计算前关闭其他占用CPU的大型程序。如果计算时间过长，可以在上述对话框中将Resolution从High降低为Medium，这会牺牲部分细节但大幅提升计算速度。

4.3 第三步：基本的三维表面观察与导航#

成功生成表面后，您可以通过鼠标操作从任意角度观察分子及其表面：

请首先观察苯分子的Hirshfeld表面颜色分布。您会发现，整体表面以白色和浅蓝色为主，几乎没有明显的红色区域。这个观察结果就是第一个重要的科学结论——苯晶体中不存在强氢键或强偶极相互作用，其分子间主要依靠较弱的范德华力和CH-π相互作用来维持堆积。请记住这个“没有红色”的苯表面，后面我们将其与乙酸（含强氢键）的表面进行对比。

4.4 第四步：生成二维指纹图#

1. 在图形窗口右侧或下侧的 Surface Controller 面板中，确认您需要分析的那个Hirshfeld表面处于选中状态（表面列表中被高亮显示）
2. 点击该面板中的 Display Fingerprint Plot 按钮（通常是一个带有双轴坐标图样的图标）
3. 一个独立的指纹图窗口将弹出，显示整张指纹散点密度图

此时您看到的是所有类型原子间接触的总和。苯的指纹图整体呈弥散状，dᵢ和dₑ主要分布在1.0到2.0 Å之间。这种形态正是典型范德华力占主导的特征。

4.5 第五步：使用过滤功能查看特定类型的相互作用#

这一步是整个分析中最能体现Hirshfeld表面方法价值的地方——我们可以“剥离”出特定原子对之间的接触，单独查看它们的贡献。

1. 在指纹图窗口中，点击菜单栏 Filter → By Element
2. 在弹出的过滤对话框中：
   • Inside（表面内部原子）：从下拉菜单中选择 C（碳）
   • Outside（表面外部原子）：从下拉菜单中选择 H（氢）
   • 勾选 Include reciprocal contacts（这将把内部H外部C的情况也计入，即CH-π相互作用的双向贡献）
3. 点击 Apply 或 OK

此时，指纹图将只显示涉及C（内部）和H（外部）的接触（以及互反的H（内部）和C（外部）接触）。对于苯分子，您会观察到在dᵢ和dₑ约1.8-2.0 Å的区域出现了延展的“翅膀”状分布——这正是CH-π相互作用的特征指纹。同时，窗口中还会显示该类接触占总所有接触的百分比数值。

💡 实用技巧：依次对H-H、O-H（若分子含氧）、C-H、C-C等不同原子对进行过滤查看，并记录各类型接触的百分比，这些数据是您在论文结果与讨论部分的重要素材。通常建议至少统计H-H、O-H/H-O、C-H/H-C、C-C以及N-H/H-N（若含氮）这五类接触的数据。

4.6 第六步：生成静电势表面（进阶操作）#

除dₙₒᵣₘ之外，静电势表面是另一类常用的显示模式。它可以帮助您理解分子表面哪些区域带正电荷（亲电，容易吸引负电基团）或负电荷（亲核，容易吸引正电基团）。

切换至静电势表面的方法：

1. 在 Surface Controller 面板中，找到您已生成的Hirshfeld表面条目
2. 点击其旁边的 Properties 下拉菜单，选择 Electrostatic Potential
3. 软件将重新绘制表面，此时颜色含义变为：
   • 🔴 红色：静电势为负（富电子区域）
   • 🔵 蓝色：静电势为正（缺电子区域）

提示：为了更好地突出静电势的分布细节，建议在Surface Controller中将显示范围手动调整为 ±0.025 au（原子单位），这个范围能清晰显示强电负性原子（如氧、卤素）周围的负电势区域以及氢原子附近的正电势区域。

5.1 三维表面的判读逻辑：三步定位法#

第一步：寻找红色区域

观察整个Hirshfeld表面的颜色分布，用肉眼找出所有呈现红色的区域。这些地方就是分子间距离最短、相互作用最强的部位。

第二步：确认红色区域的原子归属

将鼠标悬停在红色区域上，或者右键点击该区域选择显示邻近原子，确认这个红色凹陷是由目标分子的哪个原子与周围分子的哪个原子相互靠近所造成的。例如，如果红色区域出现在羰基氧附近，且外侧邻近的是一个氢原子，那么这就很可能是经典的C=O···H氢键。

第三步：判断相互作用的性质

结合dₙₒᵣₘ的绝对值大小、红色区域的面积以及原子对类型，综合判断该相互作用的性质：

• 红色区域深且面积小：为强方向性相互作用（如氢键），可能对晶体堆积方向起决定性作用
• 红色区域浅且面积大：为较弱的面-面接触（如π-π堆积），可能主导层状堆积方向
• 多个分散的浅红色区域：为多重弱相互作用协同作用

5.2 指纹图的判读逻辑：从整体到局部#

判读指纹图的正确顺序是先整体，后局部。

第一步：观察整体形状

• 若整体呈鹰形弥散分布，无任何尖峰或翅膀 → 范德华力占绝对主导
• 若整体形状左下角出现明显的尖峰 → 存在氢键网络
• 若整体形状出现左右不对称的翼状延伸 → 存在CH-π或π-π堆积

第二步：查看百分比数据

过滤后的各类接触百分比是高度可量化的数据。一个实用的经验数据参考如下（以常见有机小分子晶体为例）：

第三步：关键对比分析

• 同类结构间的对比：例如，您合成了两个结构相似的配合物，但一个含结晶水而另一个不含。对比二者的O-H接触百分比，可以定量描述水分子对氢键网络的增强作用。
• 多晶型之间的对比：同一化合物的不同晶型，即使指纹图整体形态相似，在dᵢ和dₑ的峰值位置上也往往存在微小偏移，这些偏移对应着晶胞参数差异导致的分子间距变化。

5.3 一个需要特别警惕的误区#

在对含卤素（如Cl、Br、I）或含硫的分子进行Hirshfeld表面分析时，初学者常常预期会看到C-X···X或C-X···H等卤键或硫键的明显信号。然而，实际结果往往令人失望——这些含有重原子的接触在指纹图中的百分比贡献通常很低（往往低于5%），甚至完全无法通过肉眼在三维表面上找到对应的大面积红色区域。

这并不意味着晶体中不存在卤键，而是反映了Hirshfeld分析方法本身的一个固有特性：该方法基于电子密度分配，而重原子的范德华半径较大，导致即使存在较强的卤键，其dₙₒᵣₘ值也往往因分母偏大而被“稀释”，难以在红-白-蓝着色中形成醒目的对比。

正确的应对策略是：不要因为指纹图中卤素相关接触的百分比低就轻易排除卤键的存在。相反，您应当结合静电势表面来检测卤素原子的σ-hole区域，如果卤素原子的外侧（沿C-X键延长方向）呈现蓝色正电势区域，且该方向指向邻近的富电子原子（如O或N），则可以合理推断存在卤键相互作用。这一结论需要通过独立实验证据（如拉曼光谱或太赫兹光谱）或DFT理论计算来进一步验证。

Q1：计算波函数时提示错误或卡住不动怎么办？#

可能原因：CrystalExplorer内置的Tonto程序在某些Windows系统上可能由于路径包含中文字符而运行异常。

解决方案：

1. 检查您的CIF文件路径，确保不包含任何中文字符或空格（例如放在 D:\data\ 而非 D:\我的数据\）
2. 检查文件命名是否包含特殊字符，建议仅使用英文字母、数字和下划线
3. 如仍无法解决，可尝试重新安装软件并重启电脑

Q2：生成表面后分子显示不完整，有的原子悬在半空中？#

原因：您可能没有正确补全分子就直接生成了表面。

解决方案：

1. 删除已生成的不完整表面（在Surface Controller中选中后右键删除）
2. 重新选择目标分子的所有原子
3. 执行 Edit → Complete Molecule
4. 确认补全后的分子结构合理（所有原子成键完整）后，再重新生成Hirshfeld表面

Q3：如何将分析结果导出用于论文插图？#

导出三维表面图像：

• 在三维窗口中调整到您满意的视角和表面显示模式
• 点击菜单栏 File → Export → Image
• 选择图片格式（建议使用PNG或TIFF以获得高清质量）
• 设定分辨率（建议300 dpi用于出版）

导出指纹图数据：

• 在指纹图窗口中点击 File → Export Data
• 可选择导出为CSV格式（可用Excel或Origin打开）或矢量图格式（SVG）

Q4：如何快速找到与红色区域相互作用的邻近分子？#

操作方法：

1. 在三维窗口中，右键点击目标红色区域
2. 在弹出的菜单中选择 Generate External Fragment
3. 软件将自动生成并显示与该区域发生接触的邻近分子片段
4. 您可以通过鼠标旋转来观察这两个分子之间相互作用的精确几何构型

Q5：对于含有多个独立分子的体系（如Z’ > 1），需要注意什么？#

当不对称单元中含有两个或以上晶体学独立的分子时（例如文献[53]中的案例），必须为每一个独立分子分别生成Hirshfeld表面。然后，您可以比较不同分子的指纹图——如果它们的指纹图几乎完全重叠，说明分子在晶体中的环境相似；如果存在明显差异（如一个分子的指纹图出现尖峰而另一个没有），则说明这两个分子在晶体中扮演着不同的堆积角色，这是一个值得在论文中详细讨论的有趣现象。