计算化学公社

标题: DLPNO-CCSD(T)计算不返回输出文件 [打印本页]

作者
Author:
化学势    时间: 2017-2-27 09:50
标题: DLPNO-CCSD(T)计算不返回输出文件
我最近使用DLPNO-CCSD(T)/aug-cc-pVTZ计算单点能。
但是算到如下图的地方就一直不返回输出文件了,也不报错,差不多有2天了,机器显示ORCA还在运行,我试着交了3个类似的SP任务都卡在这里。之前,有其他方法没出现过类似事情,请问这种情况正常么?希望高人指点,谢谢!


我的体系是有一个S,去其余都是C、H、O、N的40~50个原子的体系。
下边是我的输入文件的关键词,附件中是具体的输出文件。


作者
Author:
yflchx    时间: 2017-2-27 09:57
再耐心等等?抑或用cc-pVTZ试算一次。

对于DLPNO-CCSD(T),我觉得def2-TZVPP性价比更高。
作者
Author:
yflchx    时间: 2017-2-27 10:00
本帖最后由 yflchx 于 2017-2-27 10:02 编辑

对于ORCA的几乎所有计算,我建议def2-TZVPP。速度、稳定性、准确度三者通盘考虑def2-TZVPP表现最佳。

当然RI一定要用。

作者
Author:
sobereva    时间: 2017-2-27 10:18
aug-cc-pVTZ是相当大的,弥散会令耗时暴增,尤其还是DLPNO-CCSD(T)。在确认def2-TZVPP这个档能正常跑完之前不要用aug-cc-pVTZ
作者
Author:
化学势    时间: 2017-2-27 10:19
yflchx 发表于 2017-2-27 10:00
对于ORCA的几乎所有计算,我建议def2-TZVPP。速度、稳定性、准确度三者通盘考虑def2-TZVPP表现最佳。

当 ...

def2-TZVPP确实很不错
老师让用aug-cc-pVTZ 也只能算了。
同样的方法,20个原子的体系20分钟不到就算完了,40多个就跑了这么久还不动了
作者
Author:
yflchx    时间: 2017-2-27 10:23
化学势 发表于 2017-2-27 10:19
def2-TZVPP确实很不错
老师让用aug-cc-pVTZ 也只能算了。
同样的方法,20个原子的体系20分 ...

可不能用线性思维来考虑计算时间。
作者
Author:
化学势    时间: 2017-2-27 11:11
sobereva 发表于 2017-2-27 10:18
aug-cc-pVTZ是相当大的,弥散会令耗时暴增,尤其还是DLPNO-CCSD(T)。在确认def2-TZVPP这个档能正常跑完之前 ...

好哒,谢谢我去试试Def2-TZVPP
作者
Author:
冰释之川    时间: 2017-2-27 11:41
化学势 发表于 2017-2-27 11:11
好哒,谢谢我去试试Def2-TZVPP

(, 下载次数 Times of downloads: 81)

作者
Author:
化学势    时间: 2017-2-27 19:00
sobereva 发表于 2017-2-27 10:18
aug-cc-pVTZ是相当大的,弥散会令耗时暴增,尤其还是DLPNO-CCSD(T)。在确认def2-TZVPP这个档能正常跑完之前 ...

测试了在DLPNO-CCSD(T)/cc-pVTZ方法下,1 hour跑完了。
感觉可能用aug-遇到bug了,因为又算了一次之前能跑的小体系也卡在第一张图的位置。
作者
Author:
sobereva    时间: 2017-2-28 02:56
化学势 发表于 2017-2-27 19:00
测试了在DLPNO-CCSD(T)/cc-pVTZ方法下,1 hour跑完了。
感觉可能用aug-遇到bug了,因为又算了一次之前能 ...

如果没有特殊必要就别用aug-
作者
Author:
bd1225an    时间: 2020-3-7 12:21
化学势 发表于 2017-2-27 19:00
测试了在DLPNO-CCSD(T)/cc-pVTZ方法下,1 hour跑完了。
感觉可能用aug-遇到bug了,因为又算了一次之前能 ...

你好,可以分享一下您的输入文件的格式吗?我在高斯中进行计算的,用方法DLPNO-CCSD(T)/cc-pVTZ 计算,出错信息显示是输入文件的格式不正确。谢谢您。
作者
Author:
卡开发发    时间: 2020-3-7 12:23
bd1225an 发表于 2020-3-7 12:21
你好,可以分享一下您的输入文件的格式吗?我在高斯中进行计算的,用方法DLPNO-CCSD(T)/cc-pVTZ 计算,出 ...

Gaussian不能做DLPNO-CCSD(T),这是免费的ORCA做的。
作者
Author:
bd1225an    时间: 2020-3-7 14:11
卡开发发 发表于 2020-3-7 12:23
Gaussian不能做DLPNO-CCSD(T),这是免费的ORCA做的。

谢谢您

作者
Author:
sobereva    时间: 2020-3-7 14:25
bd1225an 发表于 2020-3-7 12:21
你好,可以分享一下您的输入文件的格式吗?我在高斯中进行计算的,用方法DLPNO-CCSD(T)/cc-pVTZ 计算,出 ...

仔细看
详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://sobereva.com/490http://bbs.keinsci.com/thread-13560-1-1.html
作者
Author:
xxzj    时间: 2021-3-9 13:52
sobereva 发表于 2017-2-28 02:56
如果没有特殊必要就别用aug-

老师,您好,我想研究金属间的弱相互作用,主要包括开壳层和闭壳层的cu, zn, pt ,pd 等需要比较精确的计算单体或者二聚体的构型,能级,相互作用等等,应该选择怎样的泛函和基组,还是使用DLPNO-CCSD(T)理论进行计算与分析?谢谢老师,辛苦啦
作者
Author:
wzkchem5    时间: 2021-3-9 19:38
xxzj 发表于 2021-3-9 13:52
老师,您好,我想研究金属间的弱相互作用,主要包括开壳层和闭壳层的cu, zn, pt ,pd 等需要比较精确的计 ...

分子多大?如果分子不是特别大(比如小于100个原子且只有一两个金属),用DLPNO-CCSD(T)/cc-pVTZ差不多了,金属用一些带弥散的相对论基组,显式描述相对论效应。当然前提还是不能有特别明显的强相关
作者
Author:
xxzj    时间: 2021-3-9 21:26
wzkchem5 发表于 2021-3-9 19:38
分子多大?如果分子不是特别大(比如小于100个原子且只有一两个金属),用DLPNO-CCSD(T)/cc-pVTZ差不多了 ...

谢谢您,单分子是大约140多个原子,其中一个金属,二聚体后就280多个原子,两个金属了,想就针对两个金属间作用进行分析
作者
Author:
wzkchem5    时间: 2021-3-9 22:11
xxzj 发表于 2021-3-9 21:26
谢谢您,单分子是大约140多个原子,其中一个金属,二聚体后就280多个原子,两个金属了,想就针对两个金属 ...

如果你是用orca做的计算,那就用PBE0-D3BJ/ZORA-def2-TZVPP//PBE0-D3BJ/ZORA-def2-SV(P)吧,用ZORA显式处理相对论效应。其中金属原子要选一个大一些、可能需要带弥散的相对论基组,但是三个过渡周期都有定义的相对论基组不太好选,单点能用aug-cc-pVTZ-DK、结构优化用aug-cc-pVDZ-DK或许可以,但是这个也看看别人的意见。如果PBE0-D3BJ/ZORA-def2-SV(P)跑不动,可以考虑把PBE0换成BP86。
作者
Author:
xxzj    时间: 2021-3-10 09:47
wzkchem5 发表于 2021-3-9 22:11
如果你是用orca做的计算,那就用PBE0-D3BJ/ZORA-def2-TZVPP//PBE0-D3BJ/ZORA-def2-SV(P)吧,用ZORA显式处 ...

好滴,谢谢您




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3