5.1 从 Hessian 矩阵到红外与拉曼光谱

这些输入和起始模型为未执行的教学示例。原创配图为示意图，不是计算结果。使用前请验证版本相关语法、授权数据、数值收敛及科学模型。

5.1.1 模型、单位与来源记录

结构长度使用 Å，电子能使用 hartree，振动波数使用 cm⁻¹。核对各输出字段，区分热修正与标准态项。

5.1.2 未执行输入与明确修改片段

按照下方说明确认每个片段的母计算与替换位置；片段不自动构成独立输入。保留空行与文件来源要求。

5.1.2.1 输入块 1

%chk=water17.chk
%mem=2GB
%nprocshared=2
#p B3LYP/6-31+G(d,p) Opt=Tight Freq=Raman
   SCF=(Tight,XQC) Int=UltraFine

Water: optimized harmonic IR and Raman teaching example

0 1
O    0.000000    0.000000    0.000000
H    0.758000    0.000000    0.586000
H   -0.758000    0.000000    0.586000

5.1.3 实作研究

5.1.3.1 问题、直觉与前置知识

为什么同一种振动在红外与拉曼光谱中会呈现不同的强弱？把振动理解为所有原子核协调运动，而不是某一根化学键单独变化。红外活性反映振动过程中偶极矩的变化；拉曼活性反映极化率的变化。强红外峰不一定也是强拉曼峰。水分子足够小，可以逐一观察并解释一个弯曲和两个伸缩正常模式。前置知识包括电荷/多重度、结构优化、能量二阶导数，以及用分子查看器打开日志和检查点文件。

5.1.3.2 完整输入：water17.gjf

参见上方输入块 1。

优化与频率计算采用相同泛函、基组和积分网格。Freq=Raman 在谐振频率/红外分析之外明确要求拉曼性质。本文件还生成案例18需要的力常数。这里描述的是孤立气相分子，不能直接代表宽峰明显的液态水光谱。

5.1.3.3 编号操作流程

提交前检查三个原子、零电荷及单重态设置，并将输入、完整日志及 water17.chk 一起保存。
确认优化收敛并正常结束。读取最终频率段，不要误用中间结构或追加在文件中的旧任务。非线性三原子分子应有 \(3N-6=3\) 个振动模式。
对每个模式记录频率、约化质量、红外强度、拉曼活性及可用的去偏振信息。通过振动动画识别弯曲、对称伸缩或反对称伸缩，不要仅凭频率顺序命名。
分别绘制红外强度与拉曼活性的棒状谱，并标明各自单位。若需要平滑谱线，可统一采用明确的线型与展宽，例如半高全宽为20 cm⁻¹的洛伦兹线型；展宽属于绘图设定。
应先保存未缩放频率，再决定是否采用文献缩放因子。使用时注明来源、适用计算方法及用途。频率缩放因子不是红外强度的乘数。

5.1.3.4 输出解读与数值检查

Frequencies -- 给出以 cm⁻¹ 为单位的谐振波数；负值表示虚频。IR Inten -- 通常采用 km mol⁻¹；Raman Activ -- 通常采用 Å⁴ amu⁻¹，表示活性。模拟拉曼散射强度还涉及激光波长、温度及散射约定。未经转换的活性图不能标为实测拉曼强度。极小值应有三个实振动模式。仅作为排错范围而非本次计算答案，孤立水分子的弯曲振动通常位于约1500–1800 cm⁻¹，普通谐振计算下的 O–H 伸缩约为3500–4100 cm⁻¹。如果把200 cm⁻¹的模式称为 O–H 伸缩，应立即检查原子、结构和单位。

很小的虚频应促使检查，而不是直接删除。应检查位移方向、提高优化精度并测试积分网格的影响。实验基频还受到非谐性、温度、环境及仪器展宽影响。任意移动峰位后获得的“吻合”并不能验证计算方法。

5.1.3.5 练习

中文练习： 哪个模式的“拉曼活性/红外强度”比值最大？解释为什么该比值有量纲，且不能作为通用的可探测性指标。仅将绘图展宽改为10和40 cm⁻¹，说明哪些分子结论不随之改变。