计算化学140917例析.ppt

opt 关键词 GView中 opt关键词的使用： 在计算seadditional keywords中输入(例如opt=rfo), 并点击update更tup界面中新命令行 每一个关键词即代表具体计算方法中的一个选项！ opt 冗余内坐标优化 opt=Cartesian 笛卡尔内坐标优化， opt=Z-matrix 内坐标优化 冗余内坐标优化次数最短，内坐标次数最多。 Redundant 在冗余内坐标内用Berny 算法进行优化计算。这对具有解析梯度的方法是默认的。 Z-matrix 用内坐标进行Berny 算法的优化计算。在这种 情况下，用关键词Fopt 而不是Opt 要求程序确认进行了完全优化计算（即，含非活性变量在内的变量之间是线性无关的，并且张开的自由度符合分子对称性）。Popt 形式要求进行内坐标的部分优化。在优化的最后它还禁止频率分析，即每一点的二阶导数（通过CalcAll 选项）。 Cartesian 在笛卡尔坐标中使用Berny 算法进行优化计算。注意初始结构的输入可以用任何一种坐标系。在使用纯笛卡尔坐标系进行优化计算时，不能做部分优化或冻结变量的计算；如果需要的话，可以使用Opt=Z-matrix 选项，并在Z矩阵中把所有要优化的原子以直角坐标定义，进行混合坐标的优化计算。如果分子说明部分的Z 矩阵不含变量，并且指定了Opt=Z-matrix，那么优化计算实际上用直角坐标执行。注意，其它各种坐标系，例如距离矩阵坐标，可以用ModRedundant 选项构造。 选择坐标系关键词 Opt=Restart 从检查点文件重新开始几何优化。在这种情况下，整个执行路径部分由Opt关键词以及和原计算任务相同的选项（以及Restart）构成，不需要其它输入。 MaxCycle=N 指定优化计算最大迭代次数为N。默认是所用冗余内坐标数目的两倍（默认的流程），或是进行优化的变量个数两倍（其它流程），与20 之间的最大值。 MaxStep=N 设定优化走步的最大尺寸（也就是初始可信半径）为0.01N Bohr 或弧度。N的默认值是30。 流程关键词 初始力常数有关的关键词 除非另行指定，否则Berny 几何优化从二阶导数矩阵（也即Hessian）的初始猜测开始，Hessian 矩阵通过从原子半径和简单的价电子力场得到的原子连接关系确定。在每一点使用计算的一阶导数改善这个近似的矩阵。这一方案通常工作得很好，但是在有些情况下初始猜测可能会非常差，使得优化无法正常开始，或者在优化初期花费较多的步骤改善Hessian 而无法接近优化的结构。此外，对于过渡态优化，需要知道鞍点附近的曲率，因此默认的近似Hessian 必须在优化过程中总是加以改善。 以下选项选择的方法用于改善力常数(Gview输入前需要加opt=)： ReadFC 从检查点文件提取力常数。这通常来自较低理论级别优化计算中的最终近似力常数，或是在较低理论级别的频率计算正确得到的（更准确的）力常数（后者远远好于前者）。 CalcFC 指定在第一个点用当前的方法计算力常数（只能用于HF，MP2，CASSCF，DFT 和半经验方法）。 CalcHFFC 指定在第一个点进行HF 解析力常数的计算。CalcHFFC 和MP2 优化一起使用，相当于DFT，AM1，PM3，PM3MM，PM6，和PDDG 等方法的CalcFC。 RCFC 指定从检查点文件读入来自于频率计算的笛卡尔坐标力常数。通常情况下，最好是提取已经转化为内坐标的力常数，说明见上（CalcFC）。但是，频率计算偶尔会发生分子需要变形为低对称性的情况。这种情况下，用旧的内坐标计算的力常数无法使用，因此应当用Opt=RCFC 读取笛卡尔力常数并进行转换。注意，只有在完成频率计算后，检查点文件中才有笛卡尔力常数。在优化失败后不能使用该选项，因为近似二阶导数矩阵中的负本征值个数是错的。在后一种情况下，你可能想要从最近的结构开始，并用数值方法计算导数（见下）。ReadCartesianFC 是RCFC 的同义词。 CalcAll 指定在每一个点使用当前的方法计算力常数（只能用于HF，MP2，CASSCF，DFT 和半经验方法）。注意在收敛的结构上自动进行振动频率分析，且计算的结果做为频率任务保存。 VCD 在Hartree-Fock 或DFT 的Opt=CalcAll 优化每一个点，指定计算VCD 强度。 NoRaman 在Hartree-Fock 的Opt=CalcAll 任务每一点，指定不计算拉曼强度（因为它包含用最后一个优化点的结果进行的频率分析）。拉曼强度对每一个中间二阶导数点增加10-20%的计算量 收敛有关的关键词