20计算化学_10.计算类型及高斯实现方法_3.pptVIP

袁淑萍 国家重点实验室培育基地408室 E-mail: spyuanqdu@163.com 电话算作业类型;几何优化的目的;OPT的更多选项;GEOM的选项;Handling Difficult Optimization Cases 复杂体系的优化; 对于大体系，先用简单方法作一粗略优化、后用高级方法和高级基组优化，有时可能更快地达到目的。 后续作业可以从前步作业的结果中取得几何构型和波函数作为继续优化的起点，此时，后续作业再不需要输入分子几何（但仍需要分子电荷和自旋多重度)。 也可先用分子力学优化，作为SP的输入构型或进一步作量子化学优化的初始构型。分子力学优化相当快，机时远少于量子化学优化。;优化中途变量出错。最常见的现象是键角超越了最大限度180o： ;步数超出 检查非常容易改变的坐标和/或强烈耦合的坐标 增加几何优化的最大步数 OPT (Restart, Maxcyc=N) 在优化过程中点群改变了 检查结构和/或不使用对称性(NoSymm);批处理任务;计算作业类型;搜索过渡态的目的;过渡态优化算法;Optimization of Transition States;寻找过渡态用OPT=QST2 或 OPT=QST3;Opt=QST2 的输入;为过渡态估算Hessian矩阵;过渡结构算例;过渡结构的确认;寻找过渡态失败怎么办;需要注意的问题;计算作业类型;振动频率分析计算可以作为一个单独的作业，但更常见的是附在分子构型优化之后成为多步作业，以验证优化结果是否稳定构型或过渡态。 频率计算与构型优化用同样的计算方法和基组进行，关键词是Freq。下面是频率计算作为后继作业的输入实例： % chk=H2O.chk # HF/6-31G(d) Freq Geom=Checkpoint Test H2O Frequencies 0 1;频率计算的目的 预测分子的红外和拉曼光谱(频率和强度) 在HF水平计算时，由于忽略相关能，计算的频率高于实验值。通常在计算前或计算后利用校正因子来校正。 影响频率计算准确性的因素还有：确定力常数矩阵所需的微分技术的不准确性和平衡几何构型的不确定性等。;几何构型优化时为难进行的优化计算力常数矩阵 力常数是能量对坐标的导数。从数学的角度，导数显示着函数的变化趋势（可以根据导数值确定函数在某个范围内是变大还是变小）。 优化的目的是寻找对应能量最低的构型。所以对难进行的优化可以通过提供力常数值来指引优化的方向，使优化能够顺利完成。 ; 确定势能面上驻点性质（极小点？鞍点？） 虚频率的有无及其数目可以帮助我们判断优化结果是鞍点或能量极小构型，分子稳定构型不应有虚频率，而过渡态则应该且只有一个虚频率。 计算零点能和热化学数据，对体系总能量、 焓、熵进行校正 极化度与超极化度。;Hessian矩阵本征值与振动模式的关系;振动强度;振动分析算例;;;振动频率的校正因子;;;;;;;;;;热化学数据;Freq(ReadIsotopes)选项： Specify alternate temperature, pressure, and/or isotopes This information appears in a separate input section having the format: temp pressure [scale]????? Must be real numbers isotope mass for atom 1 isotope mass for atom 2 ... isotope mass for atom n where temp, pressure, and scale are the desired temperature, pressure, and an optional scale factor for frequency data when used for thermochemical analysis (the default is unscaled). ;%Chk=freq # HF/6-31G(d,p) Freq Test ????? Frequencies at STP ????? molecule specification ?? --Link1– %Chk=freq %NoSave # HF/6-31G(d,p) Freq(ReadIso, ReadFC) Geom=Check Test Repeat at 300 K 0,1 ????? 300.0 1