高斯计算入门讲课文档.ppt

发布时间：2024-11-15 08:54

L7模块输出: 各原子受力情况最大值均方根值现在一百二十四页，总共二百一十二页。根据受力和所采用的优化方法,确定新的构型参数四个构型收敛指标当前值标准阈值(可用loose参数修改) 现在一百二十五页，总共二百一十二页。重复上述过程对新构型计算能量及其梯度,直至满足收敛条件: 现在一百二十六页，总共二百一十二页。 (4) 构型优化说明: 要缩短构型优化时间,需尽可能给出较为准确的初始构型, 例如采用X衍射实验结果等等; 对于较大体系的构型优化, 为了缩短机时, 可采用分步优化的方法, 即首先采用较小的基组进行优化, 再将该步骤得到的构型作为初始构型, 采用较大的基组作进一步的构型优化. 该方法尤其适合于初始构型不太确定的情形; 由于构型优化涉及到多变量的优化过程, 其最终的结果受到初始构型的影响较大, 即不能保证所得的构型对应于能量极小点. 这也是由于化合物势能面的复杂性引起的, 例如: 现在一百二十七页，总共二百一十二页。 local minimum 为了保证得到的构型对应于能量极小点, 通常需在构型优化的基础上进一步进行频率计算, 若计算结果存在明显虚频, 则得到的构型并非对应于能量极小点. 现在一百二十八页，总共二百一十二页。除了初始构型对优化结果有较大影响外, 体系对称性的限制也将影响到最终得到的构型. 默认情况下, 在构型优化过程中, 体系的对称性是保持不变的, 即分子所属的点群是不变的. 例如, 以NH3为例, 如果初始构型将NH3定义为平面型(D3h), 则最终是得不到属于C3V群的NH3. 为了克服该困难, 可用nosymm关键词来取消对称性上的限制, 例如: %mem=32mb %chk=H2O #P B3LYP/6-311G OPT NOSYMM The geometry optimization of water 0,1 O H,1,1.0 H,1,1.0,2,105.4 此时, H2O对称性将放开, 不再为C2v群, 而是C1群现在一百二十九页，总共二百一十二页。 c. 振动光谱(IR)计算: 在构型优化基础上, 通过进一步计算能量的二阶导数, 可求得力常数，进而得到化合物的振荡光谱。与单点能和构型优化相比，IR计算需调用L10和11模块(包括L1002, 1101, 1102, 1110等)。 (1) IR计算的输入：关键词为FREQ，例如： %mem=32mb %chk=H2O #P B3LYP/6-311G OPT FREQ Geometry optimization and the frequency calculation of H2O 0,1 O H,1,1.0 H,1,1.0,2,105.4 在构型优化的基础上作进一步的频率计算现在一百三十页，总共二百一十二页。说明：必须注意，只有对稳定构型进行IR计算才有意义，即所要进行IR分析的构型是必须事先经过优化，而且构型优化与IR 计算必须使用相同的计算方法和基组，否则计算结果也是没有意义的。以H2O为例，要对其进行IR分析须经过下述步骤： a. 首先对H2O的构型进行优化，此时输入文件为： %mem=32mb %chk=H2O #P B3LYP/6-311G OPT Geometry optimization of H2O 0,1 O H,1,1.0 H,1,1.0,2,105.4 现在一百三十一页，总共二百一十二页。 b. 在上一步计算结果的基础上，进一步进行IR计算: %mem=32mb %chk=H2O #P B3LYP/6-311G geom=check guess=check freq Frequency calculation of H2O 0,1 i) 在该步骤所采用的计算方法和基组与前一步骤必须相同； ii) geom=check含义是从checkpoint文件(h2o.chk)中读取构型，注意的是该chk文件是在上一步计算过程中产生的，在opt 结束后，该文件保存了H2O经优化后的构型。由于已确定从chk文件中读取构型，因此在输入文件的后面不要再给出H2O的构型，否则计算将出错； iii) guess=check含义是从checkpoint文件中读取初始波函数，由于该波函数已经SCF迭代优化过，因此可直接作为初始波函数，这样可以减少FREQ计算中SCF所耗费的时间；现在一百三十二页，总共二百一十二页。上述分步计

网址：高斯计算入门讲课文档.ppt https://www.yuejiaxmz.com/news/view/81211

上一篇：11月11日宏利周期混合净值下跌

下一篇：家庭财务规划协议书.docx

高斯计算入门讲课文档.ppt

相关内容

随便看看

最新动态分享

热点动态分享

专题

推荐动态分享