这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:tssofband [2020/04/06 19:32] – [基本参数设置] liu.jun | adf:tssofband [2023/03/18 13:39] (当前版本) – [修改优化最大步长] liu.jun | ||
---|---|---|---|
行 6: | 行 6: | ||
====参数设置==== | ====参数设置==== | ||
{{ : | {{ : | ||
+ | |||
+ | 关于基组的选择,参考:[[adf: | ||
{{ : | {{ : | ||
- | 关于基组的选择,参考:[[adf: | + | **为需要计算频率的原子,创建一个Region,然后在Properties - IR(Frequencies), |
- | < | + | |
- | Properties | + | |
- | NormalModes Yes | + | |
- | SelectedAtomsForHessian 3 12 13 14 | + | |
- | End | + | |
- | </ | + | |
- | 我们添加了一行SelectedAtomsForHessian 3 12 13 14,表示我们只需要计算编号为3、12、13、14的这4个原子,也就是两个H原子和两个S原子。 | + | |
====注意:==== | ====注意:==== | ||
行 58: | 行 53: | ||
====修改优化最大步长==== | ====修改优化最大步长==== | ||
保存任务后,修改*.run文件: | 保存任务后,修改*.run文件: | ||
+ | < | ||
GeometryOptimization | GeometryOptimization | ||
Method Quasi-Newton | Method Quasi-Newton | ||
行 65: | 行 60: | ||
| | ||
End | End | ||
- | UpdateTSVectorEveryStep True | ||
End | End | ||
InitialHessian | InitialHessian | ||
行 76: | 行 70: | ||
End | End | ||
End | End | ||
+ | </ | ||
+ | 这里我们设置了最大位移量为0.01,这样可以避免优化过程中的振荡行为,这个值越小,越利用收敛,但收敛的也越慢。 | ||
- | 这里我们设置了最大位移量为0.01,这样可以避免优化过程中的振荡行为,这个值越小,越利用收敛,但收敛的也越慢。 | + | 保存*.run文件,并直接运行任务。 |
====结果查看==== | ====结果查看==== | ||
SCM - Movie打开,收敛的最后一个结构(梯度小于阈值),也就是搜索得到的精确过渡态结构(即鞍点)。File - Save Geometry,保存结构。用于下面频率验证。 | SCM - Movie打开,收敛的最后一个结构(梯度小于阈值),也就是搜索得到的精确过渡态结构(即鞍点)。File - Save Geometry,保存结构。用于下面频率验证。 | ||
行 99: | 行 95: | ||
</ | </ | ||
用三种单位显示了其能量。过渡态能量减去反应物能量即得到绝对0度下的能垒。 | 用三种单位显示了其能量。过渡态能量减去反应物能量即得到绝对0度下的能垒。 | ||
- | |||
- | 如何通过计算反应相关的部分原子的频率,近似地得到反应相关的自由能?参考:[[adf: |