3.4 一致比较气相与 SMD 溶液模型

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

3.4.1 模型、单位与来源记录

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

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

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

3.4.2.1 输入块 1

%oldchk=case10_water_min.chk
%chk=case12_water_smd_vertical.chk
%mem=4GB
%nprocshared=2
#p B3LYP/6-31G(d) SP SCRF=(SMD,Solvent=Water)
   Geom=AllCheck Guess=Read SCF=Tight Integral=UltraFine

3.4.2.2 输入块 2

%oldchk=case10_water_min.chk
%chk=case12_water_smd_min.chk
%mem=4GB
%nprocshared=2
#p B3LYP/6-31G(d) Opt=(Tight,CalcFC) Freq
   SCRF=(SMD,Solvent=Water) Geom=AllCheck Guess=Read
   SCF=Tight Integral=UltraFine Temperature=298.15 Pressure=1.0

3.4.3 实作研究

3.4.3.1 问题与直观理解

把极性分子放进可极化环境后会发生什么？隐式溶剂以连续介质包围溶质，而不是显式放入成千上万个溶剂分子。SMD 包含静电响应和经验参数化的非静电贡献，适合做受控比较，但不能显示某个具体水合氢键壳层、沿溶剂链的质子转移或溶剂交换动力学。本例以案例 10 的水单体作为小型溶质来学习流程，不把它当成液态水的微观模拟。

前提是保留已验证的 case10_water_min.chk 及对应气相输出。本例区分两个问题：固定几何比较考察核位置不变时模型能量如何变化；独立溶液优化允许结构松弛。两种差值不能混用，它们也不会自动包含完整且与浓度有关的实验自由能。

3.4.3.2 完整检查点输入

先在气相极小点几何上做固定结构的 SMD 单点。

参见上方输入块 1。

再从保留的气相检查点出发，在同一 SMD 环境中独立优化并计算频率。

参见上方输入块 2。

Geom=AllCheck 表示两个输入都不含标题或分子说明区。Guess=Read 只是提供初始轨道猜测，溶液 SCF 仍需在指定溶剂中自洽收敛。各作业使用不同输出检查点，不覆盖气相极小点或固定几何结果。应调用内置 Water 溶剂，而不是手动输入介电常数后就假定所有 SMD 参数仍完全等价。

3.4.3.3 操作步骤

把气相方法、基组、网格、SCF 设置、电荷和多重度记入比较记录；固定几何试验中有意改变的项目应是溶剂模型。
运行固定几何作业，确认输出采用 SMD 和 water，并与气相极小点核坐标比较。分子整体取向改变不等于内部结构改变；发生重定向时可比较键长与键角。
记录最终模型总能量。对于这里常规基态 SMD 计算，应核对输出中关于非静电 CDS 项已计入总能量的说明，不能再加一次 CDS。CDS 诊断项的 kcal mol−1 单位与总能量 hartree 单位应分清。
执行溶液 Opt Freq，并按案例 10 验证极小点。比较气相与溶液的 O–H 键长、H–O–H 角及偶极矩，明确它们是模型计算响应，不预先保证某个数值变化。
分别报告 ΔEfixed = ESMD(Rgas) − Egas(Rgas) 与 ΔErelaxed = ESMD(Rsolv) − Egas(Rgas)。若进一步构建溶液 Gibbs 能估计，应按案例 11 说明所用热校正及标准态转换。

3.4.3.4 检查、错误与解释

SMD 总量包含溶剂化贡献，简单称为“纯电子能差”可能掩盖其含义。反过来，把输出中的 RRHO Gibbs 合计称作完全严格的液相自由能，也会夸大对平移、转动、浓度及溶剂结构的处理能力。应选择一个明确热力学循环，例如将气相 RRHO 校正与 SMD 溶剂化贡献组合的混合循环，但必须注明参考态、几何约定，并避免重复加入溶剂化。

仅在气相优化后的频率步骤加入 SCRF，不能证明得到的是溶液极小点。只对溶液使用更大基组，会把基组效应和溶剂效应混在一起。带电物种还需要特别关注弥散函数、空穴定义和单离子溶剂化的实验约定；这个中性示例不解决这些问题。若具体配位或氢键网络主导化学过程，可能必须加入显式溶剂分子。

中文练习：仅把固定几何作业中的溶剂改成 Solvent=Ethanol，其余设置保持不变。说明哪种比较能隔离溶剂模型变化，以及若要比较完全松弛极小点和一致定义的 Gibbs 能，还需要哪些作业。不要断言能量顺序只由介电常数决定。