计算化学公社

标题: 请教orca加速问题 [打印本页]

作者
Author: dingliping1108 时间: 2016-5-23 11:16
标题: 请教orca加速问题
请教各位老师，我分别用两个orca命令计算同一个物质的基态结构优化，前面一个命令计算用了大约24h，后面一个命令现在已经计算了50h了还没有算完，而且在后面一个命令里，我用了加速命令，理论上应该比前面的命令用时更短，是我后面这个命令有错误吗？

%pal nprocs 4
end
! RHF 6-31G Opt TightSCF SlowConv
! NumFreq
! COSMO(THF)
%geom Calc_Hess true; NumHess true
end

%pal nprocs 4
end
! B3LYP RIJCOSX def2-TZVP def2-TZVP/J TightSCF Opt Numfreq
! COSMO(THF)
%geom Calc_Hess true; NumHess true
end

作者
Author: amani73986 时间: 2016-5-23 13:06
上面的有SlowConv，下面没有啊，而且上面的的HF，下面是DFT，上面是6-31G，下面TZVP基组大多了啊.

作者
Author: fhh2626 时间: 2016-5-23 14:04
"理论上应该比前面的命令用时更短" 没看出来
还有你理解"Numfreq"的意思吗？为何不用"freq"

虽然orca使用很简单，但是还是建议你把每个关键词代表的意思理解了再开始算

作者
Author: dingliping1108 时间: 2016-5-23 15:29

fhh2626 发表于 2016-5-23 14:04
"理论上应该比前面的命令用时更短" 没看出来
还有你理解"Numfreq"的意思吗？为何不用"freq"

我只知道Mum是算数值频率，freq是算解析频率，数值频率慢

作者
Author: dingliping1108 时间: 2016-5-23 15:30
本帖最后由 dingliping1108 于 2016-5-23 15:41 编辑

amani73986 发表于 2016-5-23 13:06
上面的有SlowConv，下面没有啊，而且上面的的HF，下面是DFT，上面是6-31G，下面TZVP基组大多了啊.

请问slowconv是不是在有过渡金属的时候才用啊？

作者
Author: sobereva 时间: 2016-5-23 15:50
数值频率不仅很慢精度也比解析的差。
def2-TZVP大体在6-311G(2df,p)程度，比6-31G慢得不是一星半点。
slowconv仅当你计算出现SCF不收敛再用，平时没事没必要用。

作者
Author: amani73986 时间: 2016-5-23 16:25

dingliping1108 发表于 2016-5-23 15:30
请问slowconv是不是在有过渡金属的时候才用啊？

看看sob的回答吧，数值频率就是NUMFREQ那个，一般都是方法没有解析频率的时候采用算的特别慢。

作者
Author: dingliping1108 时间: 2016-5-23 16:27
本帖最后由 dingliping1108 于 2016-5-23 16:40 编辑

sobereva 发表于 2016-5-23 15:50
数值频率不仅很慢精度也比解析的差。
def2-TZVP大体在6-311G(2df,p)程度，比6-31G慢得不是一星半点。
slo ...

谢谢sob老师指导，再请问sob老师，在寻找过度态和结构优化的时候的时候，是不是用Hybrid_Hess这个命令比用Calc_Hess true命令计算要快一些？ Hybrid_Hess是不是只计算某几个原子的hess？ Calc_Hess true命令是不是每一步都计算精确地hess矩阵？

%geom Calc_Hess true # calculation of the exact Hessian before the first optimization step

Hybrid_Hess # requests the calculation of a Hybrid Hessian for the atoms 0, 1, 5 and 6 instead of the exact Hessian.
[0 1 5 6]
end

作者
Author: dingliping1108 时间: 2016-5-23 16:41

amani73986 发表于 2016-5-23 16:25
看看sob的回答吧，数值频率就是NUMFREQ那个，一般都是方法没有解析频率的时候采用算的特别慢。

非常感谢

作者
Author: sobereva 时间: 2016-5-23 16:54

dingliping1108 发表于 2016-5-23 16:27
谢谢sob老师指导，再请问sob老师，在寻找过度态和结构优化的时候的时候，是不是用Hybrid_Hess这个命令比 ...

Hessian肯定每一步都要计算，否则优化没法做。但是Hessian可以精确计算，也可以近似估算。
Calc_Hess true是指在第一步的时候精确计算Hessian，相当于高斯的calcfc。如果以后打算每步或隔几步精确算一次Hessian，用recalc_hess
Hybrid_Hess是指计算Hessian的时候只对其中几个原子精确计算而且它的近似估算，所以比所有原子都精确计算更省时。

作者
Author: dingliping1108 时间: 2016-5-24 10:12

sobereva 发表于 2016-5-23 16:54
Hessian肯定每一步都要计算，否则优化没法做。但是Hessian可以精确计算，也可以近似估算。
Calc_Hess tr ...

谢谢sob老师指导，再请教您，您之前给我推荐的M062X这个方法，orca如果用M062X只能计算数值频率，不能计算解析频率，对于60个原子的有机物（主要原子C H O Ｎ）除了起步的B3LYP，orca还有没有别的更精确地方法？

作者
Author: sobereva 时间: 2016-5-24 15:54

dingliping1108 发表于 2016-5-24 10:12
谢谢sob老师指导，再请教您，您之前给我推荐的M062X这个方法，orca如果用M062X只能计算数值频率，不能计 ...

如果弱相互作用对体系构象影响不明显，就用B3LYP就够了，或者用纯泛函，这样在ORCA里快得多。弱相互作用若影响明显加个D3就完了。

作者
Author: dingliping1108 时间: 2016-5-24 16:50

sobereva 发表于 2016-5-24 15:54
如果弱相互作用对体系构象影响不明显，就用B3LYP就够了，或者用纯泛函，这样在ORCA里快得多。弱相互作 ...

谢谢sob老师指导

欢迎光临计算化学公社 (http://bbs.keinsci.com/)