Gaussian如何进行势能面扫描

Gaussian软件是量子化学计算的强大工具，用于探索分子系统的势能面（Potential Energy Surface, PES）。PES扫描是通过系统地改变一个或多个内部坐标（如键长、键角或二面角），计算相应能量变化，从而映射出反应路径、构象转变或振动模式等。PES扫描有助于识别最小值、过渡态和反应机理，尤其在研究化学反应、分子动力学和光谱时不可或缺。Gaussian支持两种主要PES扫描类型：刚性扫描（Rigid Scan）和柔性扫描（Relaxed Scan）。刚性扫描仅固定扫描坐标进行单点能量计算，效率高但忽略结构松弛；柔性扫描则在每个步骤优化其他坐标，更准确但计算量大。以下详细介绍如何在Gaussian中进行PES扫描，包括输入准备、执行和结果分析。建议使用GaussView辅助构建输入和可视化结果。

1. 势能面扫描的两大类型

在 Gaussian 中，扫描主要分为以下两种模式：

刚性扫描 (Rigid Scan)

• 关键词： Scan
• 特点： 只改变指定的几何参数，分子的其他所有参数保持固定。
• 适用场景： 粗略观察能量趋势，或者当分子非常大、几何优化极其耗时时使用。由于不准许其他自由度松弛，得到的势能曲线通常比实际能量偏高。

柔性扫描 (Relaxed Scan)

• 关键词： Opt=ModRedundant
• 特点： 每改变一次指定的参数，程序都会对分子的其余自由度进行几何优化。
• 适用场景： 最常用的方法。它能模拟真实的物理过程，得到的势能面更平滑，且能更准确地寻找过渡态的大致位置。

2. 准备输入文件 (Input File)

刚性扫描示例

在输入文件的 Route Section 写入 Scan。在分子坐标部分，通过指定起始值、步数和步长来定义扫描。

输入格式示例： B 1 2 S 10 0.1

• B: 代表键长 (Bond)。
• 1 2: 原子编号。
• S: 代表 Scan。
• 10: 扫描的步数。
• 1: 每步增加的长度（单位为 Angstrom）。

柔性扫描示例（推荐）

这是最常用的方式，使用 Opt=ModRedundant 并在坐标末尾定义扫描参数。

输入文件结构：

PBE1PBE/6-31G(d) Opt=ModRedundant

(此处省略标题和电荷多重度)

C 0.000 0.000 0.000 H 0.000 0.000 1.089 ...

D 1 2 3 4 S 36 10.0

• D: 代表二面角 (Dihedral)。
• 1 2 3 4: 原子编号。
• S: 代表 Scan。
• 36: 扫描步数。
• 0: 步长（单位为度 deg）。

3. 关键希腊字母与单位说明

在处理扫描结果和撰写报告时，请注意以下物理量的表述：

• 能量单位： 通常使用 Hartree（原子单位）或转换后的 kcal/mol、kJ/mol。
• 角度参数： 二面角常标注为 phi 或 psi。
• 波函数对称性： 有时涉及不可约表示，如 alpha 或 beta 自旋轨道。
• 振动分析： 若在扫描点进行频率计算，会涉及到频率字母 nu 或波长 lambda。

4. 结果处理与分析

完成计算后，您需要提取每个点对应的能量。

1. 使用 GaussView 查看： 打开 .log 或 .chk 文件，点击 "Results" -> "Scan"，可以直接看到平滑的能量曲线。
2. 寻找波峰与波谷：

• 能量最低点（Local Minimum）： 对应分子的稳定构象。
• 能量最高点（Maximum）： 往往是过渡态（Transition State）的良好初始结构。

     3.提取数据： 可以使用 grep 命令在 Linux 终端提取 "SCF Done" 的能量值，或者在 GaussView 中将数据导出为文本，方便在 Excel 或 Origin 中绘图。

    5.注意事项

• 步长选择： 步长过大可能导致丢失势能面上的重要极值点；步长过小则浪费计算资源。一般键长扫描建议05 到 0.1 Angstrom，二面角建议 5 到 15 度。
• 对称性： 如果扫描破坏了分子原本的对称性，建议在 Route Section 加入 Nosymm 关键词，防止计算中途报错。
• 多重度： 确保整个扫描过程中体系的电荷和多重度保持一致。