6-关键字分析 Gaussian 03
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1。优化步骤 2。最终结构 3。振动分析 4。电荷分布
Job_Command Line
布居数分析POP和Density
Density控制用什么方法进行布局数分析。缺省 是SCF。 可以用Density=Current指定用当前采 用的方法进行布局数分析。 MP2/6-31G* Density=Current （用MP2方法进行布局数分析） Density可以进行跃迁CIS密度分析。 Density=Transition=（N，M）默认M=0 表示：使用态M和态N之间的CIS跃迁密度。
Polar＝cubic 对解析极化率做数值微分，计算超极化率
Polar＝enonly 进行能量的二阶数值微分，计算极化率
Polar
对HF，DFT和MP2自动计算极化率和超极化率
对CIS，MP3，MP4(SDQ), CID， CISD， CCD，
CCSD， QCISD和CASSCF方法，polar将只计算
IRCMAX计算的输入
# IRCMAX(CCSD/6-311G**:MP2/6-31G*, ZERO, STEPSIZE=N)
在MP2/6-31G*反应路经上寻找点，使CCSD/6-311G**能 量位于其最大值 找到的CCSD能量最大点要进行结构优化。 ZERO选项包含低水平下零点能校正，用于产生零曲率变 化过渡态理论ZC-VTST的数据。 其他选项和IRC计算中相似（Forward，Reverse等）
复合高精度能量方法Gn，CBS和W1u
1. 复合高精度方法是专门为计算高精度的热 力学性质开发的。 2. 每种方法都包含一系列规划好的，不同水 平的计算步骤。 3. 最后的能量是按照事先拟合处的经验公式 计算的。 4. 自动的复合能量方法只对局域极小能量点 进行，不能进行过渡态计算。
CBS-4M, CBS-Lq,CBS-Q, CBS-QB3 和CBS-APNO
单点能计算的表示方法： 高水平方法/高水平基组//低水平方法/低水平基组 例如： CCSD(T)/CC-PVDZ//MP2/6-31G* 表示在MP2/6-31G*水平下进行的结构优化，以优化 的几何结构（坐标），进行CCSD(T)/CC-PVDZ下的 单点能计算
做单点能计算的方法
1. 分步进行优化和单点能计算 a. 拷贝优化的几何构型进行高水平单点能计算 b. 从Checkpoint文件里读取优化的几何构型进行单 点能计算
300.0 1.0 同位素说明
所有复合高精度能量的方法的基本依据是， 合理的结构可以给出合理的能量 依据同样的道理，可以在比标准高的优化 水平下计算几何构型，并用这个几何构型 按标准复合能量方法能量的计算步骤进行 能量计算。
根据G2的计算规则，可以手动计算过渡态 的G2能量或者反应路径上的能量计算
NBO分析（NBO, NPA, SaveNBOs, SaveNLMOs）
POP=Full 打印所有轨道的布局分析 POP=REG打印5个占据轨道，5个虚轨道 的布局分析 NBO分析 NBO 完全自然轨道布局分析 NPA 只进行NBO的自然布局分析阶段
Polar
指示进行偶极电场极化率的计算 Freq和polar可以同时进行计算，但是不能与没有解析梯 度的方法合用。 可以进行静态(或者含频)极化率和超极化率的计算
G1, G2, G2MP2, G3, G3MP2, G3B3和 G3MP2B3
Gaussian-1 G1 Gaussian-2 G2 Gaussian-3 G3 G2MP2进行改良版的G2，也就是G2 （MP2）的计算，它对基组展开修正使用 MP2代替MP4，具有接近完全G2的精度， 但是节省时间。 ReadIsotopes
0.586860 -0.801000 2.739275 3.641534
1 1 2 2
Counterpoise: corrected energy = -2660.083831739527 Counterpoise: BSSE energy = 0.003902746890
用GaussView可视化观察优化结果
POP关键字
打印控制选项（Full, Reg, Minimal）
成键性（Bonding）
自然轨道分析（NO, NOAB, AlphaNatural， BetaNatural, SpinNatural） 静电势和电荷分析（MK, CHelp, CHelpG, Dipole）
根据不同的方案进行选点分析
极化率，polar＝enonly将计算极化率和超极化率。 半经验方法不能计算polar
平衡性计算
对于弱相互作用体系，可能会有基组重叠误差 （BSSE），需要进行平衡性计算
# hf/lanl2dz counterpoise=2 nosymm opt test HBr + HF, optimization with counterpoise correction using ECP basis 010101 H -0.046866 0. Br -0.331864 0. F 0.396755 0. H 0.584835 0.
W1u 和W1BD
这两种方法比CBS-QB3和G3更加昂贵，但 是更加准确。 W1u是在W1方法修改后，对开壳层体系使 用UCCSD代替ROCCSD。 W1BD，它对有关方法用BD代替耦合簇 ReadIsotopes
%Chk=formald # G2 Test G2 on formaldehyde 01 分子说明 --Link1-%Chk=formald %NoSave # G2(Restart,ReadIso) Geom=Check Repeat at 300 K 0,1
高精度能量计算、布居 数分析
用低水平优化的结合结构进行高水 平的能量计算
量化计算的矛盾：
几何结构优化： 1.在较低的水平方法和基组下，优化得到稳定 结构是良 好的，但是能量计算精度不够。 2.高水平的方法和基组，优化的结构更好，能量计算更 精确，但是所需要的计算资源要多许多。
频率计算和热力学校正：
用于进行非常精确的能量计算,使用这些关键 词时不需要指定基组。 这些关键词进行复杂的能量计算， 可以用标准方法指定同位素，ReadIsotopes 现在CBS-4模型已经升级，使用新的局域化 程序并改善了经验参数，建议新的研究使 用新的版本CBS-4M（M表示最小布局局域 化）
CBS-APNO只适合用于第一行原子 CBS-4M,CBS-Lq.CBS-Q,CBS-QB3 只适用于第一和第二行原子
1.高水平方法，基组的频率计算对计算的内存， 硬盘的有很大的需要，而且时间很长。 2. 较低方法和基组水平下频率计算结果，经过 校正因子修正后，也是比较好的，而计算时间 会大大缩短。
为了兼顾计算精度和计算资源的矛盾， 可以 采用单点能计算的方法，即，用较低水平下 进行几何结构的优化和频率计算，再在低水 平下优化的结构上进行高水平的单点能计算
%chk=XXX.chk #p CCSD(T)/CC-PVDZ SP TEST Geom=check Title
01
2.低水平优化后自动进行单点能计算
#P CCSD/6-311G**//MP2/6-31G*
Title 01
这种Hale Waihona Puke Baidu法不能和FREQ结合使用
高水平下过渡态能量IRCMAX
高水平和低水平计算的过渡态位置不同(结 构有差别) 当高水平的过渡态结构和低水平的过渡态 结构相近时，可以用IRCMAX在低水平的 IRC路径上，自动寻找高水平下的能量最高 点。 IRCMAX从过渡态结构开始。
Job_Command Line
布居数分析POP和Density
Density控制用什么方法进行布局数分析。缺省 是SCF。 可以用Density=Current指定用当前采 用的方法进行布局数分析。 MP2/6-31G* Density=Current （用MP2方法进行布局数分析） Density可以进行跃迁CIS密度分析。 Density=Transition=（N，M）默认M=0 表示：使用态M和态N之间的CIS跃迁密度。
Polar＝cubic 对解析极化率做数值微分，计算超极化率
Polar＝enonly 进行能量的二阶数值微分，计算极化率
Polar
对HF，DFT和MP2自动计算极化率和超极化率
对CIS，MP3，MP4(SDQ), CID， CISD， CCD，
CCSD， QCISD和CASSCF方法，polar将只计算
IRCMAX计算的输入
# IRCMAX(CCSD/6-311G**:MP2/6-31G*, ZERO, STEPSIZE=N)
在MP2/6-31G*反应路经上寻找点，使CCSD/6-311G**能 量位于其最大值 找到的CCSD能量最大点要进行结构优化。 ZERO选项包含低水平下零点能校正，用于产生零曲率变 化过渡态理论ZC-VTST的数据。 其他选项和IRC计算中相似（Forward，Reverse等）
复合高精度能量方法Gn，CBS和W1u
1. 复合高精度方法是专门为计算高精度的热 力学性质开发的。 2. 每种方法都包含一系列规划好的，不同水 平的计算步骤。 3. 最后的能量是按照事先拟合处的经验公式 计算的。 4. 自动的复合能量方法只对局域极小能量点 进行，不能进行过渡态计算。
CBS-4M, CBS-Lq,CBS-Q, CBS-QB3 和CBS-APNO
单点能计算的表示方法： 高水平方法/高水平基组//低水平方法/低水平基组 例如： CCSD(T)/CC-PVDZ//MP2/6-31G* 表示在MP2/6-31G*水平下进行的结构优化，以优化 的几何结构（坐标），进行CCSD(T)/CC-PVDZ下的 单点能计算
做单点能计算的方法
1. 分步进行优化和单点能计算 a. 拷贝优化的几何构型进行高水平单点能计算 b. 从Checkpoint文件里读取优化的几何构型进行单 点能计算
300.0 1.0 同位素说明
所有复合高精度能量的方法的基本依据是， 合理的结构可以给出合理的能量 依据同样的道理，可以在比标准高的优化 水平下计算几何构型，并用这个几何构型 按标准复合能量方法能量的计算步骤进行 能量计算。
根据G2的计算规则，可以手动计算过渡态 的G2能量或者反应路径上的能量计算
NBO分析（NBO, NPA, SaveNBOs, SaveNLMOs）
POP=Full 打印所有轨道的布局分析 POP=REG打印5个占据轨道，5个虚轨道 的布局分析 NBO分析 NBO 完全自然轨道布局分析 NPA 只进行NBO的自然布局分析阶段
Polar
指示进行偶极电场极化率的计算 Freq和polar可以同时进行计算，但是不能与没有解析梯 度的方法合用。 可以进行静态(或者含频)极化率和超极化率的计算
G1, G2, G2MP2, G3, G3MP2, G3B3和 G3MP2B3
Gaussian-1 G1 Gaussian-2 G2 Gaussian-3 G3 G2MP2进行改良版的G2，也就是G2 （MP2）的计算，它对基组展开修正使用 MP2代替MP4，具有接近完全G2的精度， 但是节省时间。 ReadIsotopes
0.586860 -0.801000 2.739275 3.641534
1 1 2 2
Counterpoise: corrected energy = -2660.083831739527 Counterpoise: BSSE energy = 0.003902746890
用GaussView可视化观察优化结果
POP关键字
打印控制选项（Full, Reg, Minimal）
成键性（Bonding）
自然轨道分析（NO, NOAB, AlphaNatural， BetaNatural, SpinNatural） 静电势和电荷分析（MK, CHelp, CHelpG, Dipole）
根据不同的方案进行选点分析
极化率，polar＝enonly将计算极化率和超极化率。 半经验方法不能计算polar
平衡性计算
对于弱相互作用体系，可能会有基组重叠误差 （BSSE），需要进行平衡性计算
# hf/lanl2dz counterpoise=2 nosymm opt test HBr + HF, optimization with counterpoise correction using ECP basis 010101 H -0.046866 0. Br -0.331864 0. F 0.396755 0. H 0.584835 0.
W1u 和W1BD
这两种方法比CBS-QB3和G3更加昂贵，但 是更加准确。 W1u是在W1方法修改后，对开壳层体系使 用UCCSD代替ROCCSD。 W1BD，它对有关方法用BD代替耦合簇 ReadIsotopes
%Chk=formald # G2 Test G2 on formaldehyde 01 分子说明 --Link1-%Chk=formald %NoSave # G2(Restart,ReadIso) Geom=Check Repeat at 300 K 0,1
高精度能量计算、布居 数分析
用低水平优化的结合结构进行高水 平的能量计算
量化计算的矛盾：
几何结构优化： 1.在较低的水平方法和基组下，优化得到稳定 结构是良 好的，但是能量计算精度不够。 2.高水平的方法和基组，优化的结构更好，能量计算更 精确，但是所需要的计算资源要多许多。
频率计算和热力学校正：
用于进行非常精确的能量计算,使用这些关键 词时不需要指定基组。 这些关键词进行复杂的能量计算， 可以用标准方法指定同位素，ReadIsotopes 现在CBS-4模型已经升级，使用新的局域化 程序并改善了经验参数，建议新的研究使 用新的版本CBS-4M（M表示最小布局局域 化）
CBS-APNO只适合用于第一行原子 CBS-4M,CBS-Lq.CBS-Q,CBS-QB3 只适用于第一和第二行原子
1.高水平方法，基组的频率计算对计算的内存， 硬盘的有很大的需要，而且时间很长。 2. 较低方法和基组水平下频率计算结果，经过 校正因子修正后，也是比较好的，而计算时间 会大大缩短。
为了兼顾计算精度和计算资源的矛盾， 可以 采用单点能计算的方法，即，用较低水平下 进行几何结构的优化和频率计算，再在低水 平下优化的结构上进行高水平的单点能计算
%chk=XXX.chk #p CCSD(T)/CC-PVDZ SP TEST Geom=check Title
01
2.低水平优化后自动进行单点能计算
#P CCSD/6-311G**//MP2/6-31G*
Title 01
这种Hale Waihona Puke Baidu法不能和FREQ结合使用
高水平下过渡态能量IRCMAX
高水平和低水平计算的过渡态位置不同(结 构有差别) 当高水平的过渡态结构和低水平的过渡态 结构相近时，可以用IRCMAX在低水平的 IRC路径上，自动寻找高水平下的能量最高 点。 IRCMAX从过渡态结构开始。