6-关键字分析 Gaussian 03解析

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

复合高精度能量方法Gn,CBS和W1u
1. 复合高精度方法是专门为计算高精度的热 力学性质开发的。 2. 每种方法都包含一系列规划好的,不同水 平的计算步骤。 3. 最后的能量是按照事先拟合处的经验公式 计算的。 4. 自动的复合能量方法只对局域极小能量点 进行,不能进行过渡态计算。
CBS-4M, CBS-Lq,CBS-Q, CBS-QB3 和CBS-APNO
1.高水平方法,基组的频率计算对计算的内存, 硬盘的有很大的需要,而且时间很长。 2. 较低方法和基组水平下频率计算结果,经过 校正因子修正后,也是比较好的,而计算时间 会大大缩短。
为了兼顾计算精度和计算资源的矛盾, 可以 采用单点能计算的方法,即,用较低水平下 进行几何结构的优化和频率计算,再在低水 平下优化的结构上进行高水平的单点能计算
用于进行非常精确的能量计算,使用这些关键 词时不需要指定基组。 这些关键词进行复杂的能量计算, 可以用标准方法指定同位素,ReadIsotopes 现在CBS-4模型已经升级,使用新的局域化 程序并改善了经验参数,建议新的研究使 用新的版本CBS-4M(M表示最小布局局域 Baidu Nhomakorabea)
CBS-APNO只适合用于第一行原子 CBS-4M,CBS-Lq.CBS-Q,CBS-QB3 只适用于第一和第二行原子
IRCMAX计算的输入
# IRCMAX(CCSD/6-311G**:MP2/6-31G*, ZERO, STEPSIZE=N)
在MP2/6-31G*反应路经上寻找点,使CCSD/6-311G**能 量位于其最大值 找到的CCSD能量最大点要进行结构优化。 ZERO选项包含低水平下零点能校正,用于产生零曲率变 化过渡态理论ZC-VTST的数据。 其他选项和IRC计算中相似(Forward,Reverse等)
%chk=XXX.chk #p CCSD(T)/CC-PVDZ SP TEST Geom=check Title
01
2.低水平优化后自动进行单点能计算
#P CCSD/6-311G**//MP2/6-31G*
Title 01
这种方法不能和FREQ结合使用
高水平下过渡态能量IRCMAX
高水平和低水平计算的过渡态位置不同(结 构有差别) 当高水平的过渡态结构和低水平的过渡态 结构相近时,可以用IRCMAX在低水平的 IRC路径上,自动寻找高水平下的能量最高 点。 IRCMAX从过渡态结构开始。
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
布居数分析POP和Density
Density控制用什么方法进行布局数分析。缺省 是SCF。 可以用Density=Current指定用当前采 用的方法进行布局数分析。 MP2/6-31G* Density=Current (用MP2方法进行布局数分析) Density可以进行跃迁CIS密度分析。 Density=Transition=(N,M)默认M=0 表示:使用态M和态N之间的CIS跃迁密度。
POP关键字
打印控制选项(Full, Reg, Minimal)
成键性(Bonding)
自然轨道分析(NO, NOAB, AlphaNatural, BetaNatural, SpinNatural) 静电势和电荷分析(MK, CHelp, CHelpG, Dipole)
根据不同的方案进行选点分析
300.0 1.0 同位素说明
所有复合高精度能量的方法的基本依据是, 合理的结构可以给出合理的能量 依据同样的道理,可以在比标准高的优化 水平下计算几何构型,并用这个几何构型 按标准复合能量方法能量的计算步骤进行 能量计算。
根据G2的计算规则,可以手动计算过渡态 的G2能量或者反应路径上的能量计算
单点能计算的表示方法: 高水平方法/高水平基组//低水平方法/低水平基组 例如: CCSD(T)/CC-PVDZ//MP2/6-31G* 表示在MP2/6-31G*水平下进行的结构优化,以优化 的几何结构(坐标),进行CCSD(T)/CC-PVDZ下的 单点能计算
做单点能计算的方法
1. 分步进行优化和单点能计算 a. 拷贝优化的几何构型进行高水平单点能计算 b. 从Checkpoint文件里读取优化的几何构型进行单 点能计算
G1, G2, G2MP2, G3, G3MP2, G3B3和 G3MP2B3
Gaussian-1 G1 Gaussian-2 G2 Gaussian-3 G3 G2MP2进行改良版的G2,也就是G2 (MP2)的计算,它对基组展开修正使用 MP2代替MP4,具有接近完全G2的精度, 但是节省时间。 ReadIsotopes
高精度能量计算、布居 数分析
用低水平优化的结合结构进行高水 平的能量计算
量化计算的矛盾:
几何结构优化: 1.在较低的水平方法和基组下,优化得到稳定 结构是良 好的,但是能量计算精度不够。 2.高水平的方法和基组,优化的结构更好,能量计算更 精确,但是所需要的计算资源要多许多。
频率计算和热力学校正:
Polar=cubic 对解析极化率做数值微分,计算超极化率
NBO分析(NBO, NPA, SaveNBOs, SaveNLMOs)
POP=Full 打印所有轨道的布局分析 POP=REG打印5个占据轨道,5个虚轨道 的布局分析 NBO分析 NBO 完全自然轨道布局分析 NPA 只进行NBO的自然布局分析阶段
Polar
指示进行偶极电场极化率的计算 Freq和polar可以同时进行计算,但是不能与没有解析梯 度的方法合用。 可以进行静态(或者含频)极化率和超极化率的计算
相关文档
最新文档