Dmol3计算过程中频率计算会占大部分计算时间么

石头圭

由于之前计算时间都比较快，没仔细看过输出文件，最近计算速度特别慢，计算过程中看输出文件时，发现大部分时间都是在计算频率，一个65原子的模型，每个原子都会进行频率优化，基本一个原子的优化要用0.45-1小时.24核计算下来要花60h左右，用的是Linux2019版MS，请问这样的速度正常么，有什么办法可以缩短计算时间么。输入文件如下

# Task parameters
Calculate                     optimize
Opt_energy_convergence        1.0000e-05
Opt_gradient_convergence      2.0000e-03 A
Opt_displacement_convergence  5.0000e-03 A
Opt_iterations                50
Opt_max_displacement          0.3000 A

# Cartesian constraints
Opt_fixed
1   XYZ
2   XYZ
4   XYZ
5   XYZ
8   XYZ
9   XYZ
12   XYZ
13   XYZ
16   XYZ
17   XYZ
20   XYZ
21   XYZ
24   XYZ
25   XYZ
28   XYZ
29   XYZ
32   XYZ
33   XYZ
36   XYZ
37   XYZ
40   XYZ
41   XYZ
44   XYZ
45   XYZ
48   XYZ
49   XYZ
52   XYZ
53   XYZ
56   XYZ
57   XYZ
60   XYZ
61   XYZ

Symmetry                      off
Max_memory                    24000
File_usage                    smart
Scf_density_convergence       1.000000e-06
Scf_charge_mixing             1.000000e-01
Scf_spin_mixing               5.000000e-01
Scf_diis                      6 pulay
Scf_iterations                2000
Preconditioner                4.000000

# Electronic parameters
Spin_polarization             unrestricted
Charge                        0
Basis                         dnp
Pseudopotential               ecp
Functional                    pbe
Aux_density                   hexadecapole
Dftd                          G06
Integration_grid              fine
Occupation                    thermal 0.0050
Cutoff_Global                 4.6000 angstrom
# Kpoint definition file (intervals/offset):
Kpoints                       file     2 2 1 0.0000 0.0000 0.0000
BSi.kpoints
# Calculated properties
Frequency_analysis            old

卡开发发

有限差分的方法65原子大致相当于算65*6个单点能，假定10分钟算一次SCF，大体上这个结果还是合理的。当然，freq做有限差分的时候，个别原子的移动可能确实导致会SCF收敛偏慢。周期体系的频率计算确实是十分昂贵的。

石头圭

卡开发发 发表于 2020-8-14 17:32
有限差分的方法65原子大致相当于算65*6个单点能，假定10分钟算一次SCF，大体上这个结果还是合理的。当然，f ...

卡开发发老师，确实是用的有限差分计算的，有什么办法可以缩短计算时间么，现在这个计算速度有点过于慢了...

卡开发发

石头圭 发表于 2020-8-14 17:37
卡开发发老师，确实是用的有限差分计算的，有什么办法可以缩短计算时间么，现在这个计算速度有点过于慢了 ...

有人是计算partial Hessian来代替Hessian的。一般来说，找过渡态使用初猜Hessian，用partial Hessian是合理的。若是计算热化学/热力学，确实有人这么做，但我个人觉得可靠性存在争议。

石头圭

卡开发发 发表于 2020-8-14 17:32
有限差分的方法65原子大致相当于算65*6个单点能，假定10分钟算一次SCF，大体上这个结果还是合理的。当然，f ...

卡开发发

石头圭 发表于 2020-8-14 17:58
别人给我了个方法，是使用部分hessian的方法计算非锁定的关键原子的频率，来提高计算速度。这样可行么

上面提到了，看做何种任务。

石头圭

卡开发发 发表于 2020-8-14 17:59
上面提到了，看做何种任务。

不好意思，刚刚没有刷新。我是用来找过渡态的，所以应该是可以的，那么比如我是四层原子，我只计算最上层的原子以及吸附分子频率是否足够说明问题。

卡开发发

石头圭 发表于 2020-8-14 18:03
不好意思，刚刚没有刷新。我是用来找过渡态的，所以应该是可以的，那么比如我是四层原子，我只计算最上层 ...

如果你是作为TS opt的Hessian初猜，原则上是可以，但选取原子的范围就非常依赖经验，最好是反应前后显著参与反应的原子尽可能多地考虑进来，否则TS opt可能会失败。

石头圭

卡开发发 发表于 2020-8-14 19:59
如果你是作为TS opt的Hessian初猜，原则上是可以，但选取原子的范围就非常依赖经验，最好是反应前后显著 ...

多谢大佬，还有个问题要麻烦一下您，我找了一下partial hessian的设置或者关键词都没有找到，请问您那里有具体步骤么

卡开发发

石头圭 发表于 2020-8-14 20:49
多谢大佬，还有个问题要麻烦一下您，我找了一下partial hessian的设置或者关键词都没有找到，请问您那里 ...

选中要算Hessian的原子，DMol3 Calculations->properties->frequency里面->add，然后再勾选use partial Hessian。

石头圭

卡开发发 发表于 2020-8-14 21:07
选中要算Hessian的原子，DMol3 Calculations->properties->frequency里面->add，然后再勾选use partial H ...

好的，万分感谢，解决了头疼了很久的问题了，谢谢！

石头圭

卡开发发 发表于 2020-8-14 21:07
选中要算Hessian的原子，DMol3 Calculations->properties->frequency里面->add，然后再勾选use partial H ...

我弄了partial hessian以后输出文件提示了这个，这是没有使用成功的意思么

Warning:
Hessian in Internals can not be calculated for
Periodic System
DMol automatically will use Cartesian option

卡开发发

石头圭 发表于 2020-8-14 23:31
我弄了partial hessian以后输出文件提示了这个，这是没有使用成功的意思么

Warning:

不是，这个就是字面翻译：周期体系Hessian没办法使用内坐标来做，只能用笛卡尔坐标，程序将会自动使用笛卡尔坐标。

石头圭

卡开发发 发表于 2020-8-14 23:42
不是，这个就是字面翻译：周期体系Hessian没办法使用内坐标来做，只能用笛卡尔坐标，程序将会自动使用笛 ...

好的，谢谢

Reminder

老师您好，请问用Linux系统计算Dmol3时候在.bashrc文件写入什么脚本计算？之前一直是算高斯的，算这个困扰了很久，万分感谢!