6-关键字分析-Gaussian-03解析

高精度能量计算、布居数分析用低水平优化的结合构造进展高水平的能量计算量化计算的冲突：几何构造优化：1.在较低的水平方法和基组下，优化得到稳定构造是良好的，但是能量计算精度不够。2.高水平的方法和基组，优化的构造更好，能量计算更准确，但是所需要的计算资源要多很多。频率计算和热力学校正：1.高水平方法，基组的频率计算对计算的内存，硬盘的有很大的需要，而且时间很长。2.较低方法和基组水平下频率计算结果，经过校正因子修正后，也是比较好的，而计算时间会大大缩短。为了兼顾计算精度和计算资源的冲突，可以承受单点能计算的方法，即，用较低水平下进展几何构造的优化和频率计算，再在低水平下优化的构造上进展高水平的单点能计算单点能计算的表示方法：高水平方法/高水平基组//低水平方法/低水平基组例如：CCSD(T)/CC-PVDZ//MP2/6-31G*表示在MP2/6-31G*水平下进展的构造优化，以优化的几何构造〔坐标〕，进展CCSD(T)/CC-PVDZ下的单点能计算做单点能计算的方法1.分步进展优化和单点能计算a.拷贝优化的几何构型进展高水平单点能计算b.从Checkpoint文件里读取优化的几何构型进展单点能计算%chk=XXX.chk#pCCSD(T)/CC-PVDZSPTESTGeom=check

Title012.低水平优化后自动进展单点能计算#PCCSD/6-311G**//MP2/6-31G*Title01这种方法不能和FREQ结合使用高水平下过渡态能量IRCMAX高水平和低水平计算的过渡态位置不同(构造有差异)当高水平的过渡态构造和低水平的过渡态构造相近时，可以用IRCMAX在低水平的IRC路径上，自动查找高水平下的能量最高点。IRCMAX从过渡态构造开头。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复合高精度方法是特地为计算高精度的热力学性质开发的。每种方法都包含一系列规划好的，不同水平的计算步骤。最终的能量是依据事先拟合处的阅历公式计算的。自动的复合能量方法只对局域微小能量点进展，不能进展过渡态计算。CBS-4M,CBS-Lq,CBS-Q,CBS-QB3和CBS-APNO用于进展特别准确的能量计算,使用这些关键词时不需要指定基组。这些关键词进展简单的能量计算，可以用标准方法指定同位素，ReadIsotopes现在CBS-4模型已经升级，使用新的局域化程序并改善了阅历参数，建议新的争论使用新的版本CBS-4M〔M表示最小布局局域化〕CBS-APNO只适合用于第一行原子CBS-4M,CBS-Lq.CBS-Q,CBS-QB3只适用于第一和其次行原子

G1,G2,G2MP2,G3,G3MP2,G3B3和G3MP2B3Gaussian-1G1Gaussian-2G2Gaussian-3G3G2MP2进展改进版的G2，也就是G2〔MP2〕的计算，它对基组开放修正使用MP2代替MP4，具有接近完全G2的精度，但是节省时间。ReadIsotopesW1u和W1BD这两种方法比CBS-QB3和G3更加昂贵，但是更加准确。W1u是在W1方法修改后，对开壳层体系使用UCCSD代替ROCCSD。W1BD，它对有关方法用BD代替耦合簇ReadIsotopes%Chk=formald#G2TestG2onformaldehyde01分子说明--Link1--%Chk=formald%NoSave#G2(Restart,ReadIso)Geom=CheckRepeatat300K0,1300.01.0同位素说明全部复合高精度能量的方法的根本依据是，合理的构造可以给出合理的能量依据同样的道理，可以在比标准高的优化水平下计算几何构型，并用这个几何构型按标准复合能量方法能量的计算步骤进展能量计算。依据G2的计算规章，可以手动计算过渡态的G2能量或者反响路径上的能量计算布居数分析POP和DensityDensity掌握用什么方法进展布局数分析。缺省是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〕依据不同的方案进展选点分析NBO分析〔NBO,NPA,SaveNBOs,SaveNLMOs〕POP=Full打印全部轨道的布局分析POP=REG打印5个占据轨道，5个虚轨道的布局分析NBO分析NBO完全自然轨道布局分析NPA只进展NBO的自然布局分析阶段Polar指示进展偶极电场极化率的计算Freq和polar可以同时进展计算，但是不能与没有解析梯度的方法合用。可以进展静态(或者含频)极化率和超极化率的计算Polar＝cubic对解析极化率做数值微分，计算超极化率Polar＝enonly进展能量的二阶数值微分，计算极化率Polar对HF，DFT和MP2自动计算极化率和超极化率对CIS，MP3，MP4(SDQ),CID，CISD，CCD，CCSD，QCISD和CASSCF方法，polar将只计算极化率，polar＝enonly将计算极化率和超极化率。半阅历方法不能计算polar平衡性计算对于弱相互作用体系，可能会有基组重叠误差〔BSSE〕，需要进展平衡性计算#hf/lanl2dzcounterpoise=2nosymmopttestHBr+HF,optimizationwithcounterpoisecorrectionusingECPbasis010101H-0.0468660.0.5868601Br-0.3318640.-0.8010001F0.3967550.2.7392752H0.5848350.3.6415342Counterp