计算化学公社

标题: 求助：Dalton CC2 激发态计算报错 [打印本页]

作者
Author: zizi 时间: 2024-6-23 00:03
标题: 求助：Dalton CC2 激发态计算报错
各位老师好，我尝试使用Dalton对过渡金属配合物的激发态进行CC2计算。前面的波函数计算已经收敛，但计算到CC2的部分时出现如下错误。我在作业运行过程中监控了内存使用情况，没有发现使用超过设定内存的情况。我在手册中搜索相关关键词，没有找到对于“CCSD_SORTAO NUMBAT”、"CCRDAO"等关键词的说明，所以不清楚应该如何纠错。初学Dalton，恳请各位老师指导，非常感谢！

作者
Author: zjxitcc 时间: 2024-6-23 09:01
信息严重缺失，例如
（1）提交Dalton任务时，dalton -gb 后面写了几？（使用几GB内存）
（2）Dalton编译时有没加--int64？
（3）体系多少原子，多少基函数？
（4）.dal文件关键词是什么，是否合理？

作者
Author: zizi 时间: 2024-6-24 12:40

zjxitcc 发表于 2024-6-23 09:01
信息严重缺失，例如
（1）提交Dalton任务时，dalton -gb 后面写了几？（使用几GB内存）
（2）Dalton编译 ...

谢谢老师的回复！抱歉没有提供足够的信息，下面是一些详细信息：
（1）我在提交Dalton任务时候，写了“dalton -gb 960”，指定使用960Gb。
（2）我使用的Dalton版本为2020.1。是按照Sobereva老师的博文“量子化学程序Dalton的编译方法和运行方式简介”（http://sobereva.com/463）中的“1.4 编译Dalton的串行+MKL并行版”进行编译，编译时没有加“--int64”。请问老师我是否应该加入“--int64”重新进行编译？
（3）我的体系有91个原子，含有重金属，共896个基函数。
（4）这是我的.dal文件：

**DALTON
.RUN WAVE FUNCTIONS
**WAVE FUNCTION
.CC
*SCF INP
.MAX DIIS
1000
*CC INP
.CC2
*CCEXCI
.NCCEXCI
2
**END OF DALTON

复制代码

感谢老师的耐心指导！

作者
Author: zjxitcc 时间: 2024-6-24 12:48
本帖最后由 zjxitcc 于 2024-6-24 12:53 编辑

zizi 发表于 2024-6-24 12:40
谢谢老师的回复！抱歉没有提供足够的信息，下面是一些详细信息：
（1）我在提交Dalton任务时候，写了“d ...

没有加--int64编译的Dalton最多只能用16GB内存，因此-gb 960容易造成报错，也没有意义（暂且不论你的机器是否真有960GB内存可用）。从你的体系大小看有可能16GB内存不够，需要加--int64重新编译和安装Dalton。我有点怀疑Dalton的CC2是否真能算得动这么大的体系。

另外，我记得Dalton的CC默认不冻结芯轨道，需要用户自己指定。注意通常情况下一定要指定，否则徒增计算量结果还未必好。要自己在输出文件里检查被冻结芯轨道的数目。

作者
Author: zizi 时间: 2024-7-2 12:45

zjxitcc 发表于 2024-6-24 12:48
没有加--int64编译的Dalton最多只能用16GB内存，因此-gb 960容易造成报错，也没有意义（暂且不论你的机器 ...

谢谢老师，我按照您的建议重新编译了Dalton，现在可以使用更大的内存了。我参考gaussian结果，判断了考虑电子相关的轨道范围，设置冻结了46个轨道。但是由于体系较大，设置基组为def2-SVP，CC2计算进行了7天，coupled cluster CC2 iteration只完成了8圈，没有收敛。

我希望使用比TDDFT精度更高的方法计算激发态，对于这么大的体系，请问我应该如何改进计算的设置，才能算得动？或者，请问老师是否有其它更好的方法值得尝试？非常感谢您的耐心指导！

作者
Author: wjc404 时间: 2024-7-2 13:31
本帖最后由 wjc404 于 2024-7-2 13:32 编辑

可以试试MRCC程序，其CC2方法有一些加速的策略，详见http://bbs.keinsci.com/thread-29156-1-1.html

也可以考虑双杂化泛函的TDDFT，使用ORCA软件开RI去算。

作者
Author: zizi 时间: 2024-7-2 16:08

wjc404 发表于 2024-7-2 13:31
可以试试MRCC程序，其CC2方法有一些加速的策略，详见http://bbs.keinsci.com/thread-29156-1-1.html

也 ...

谢谢老师！我试试您建议的方法。

作者
Author: Yaoyl 时间: 2024-12-10 14:49

zjxitcc 发表于 2024-6-23 09:01
信息严重缺失，例如
（1）提交Dalton任务时，dalton -gb 后面写了几？（使用几GB内存）
（2）Dalton编译 ...

老师您好，我想问一下您通过Multiwfn软件做出双光子光谱图后，它的纵坐标是摩尔吸收系数，怎么才能转换成双光子吸收系数呢？

作者
Author: zjxitcc 时间: 2024-12-10 14:55

Yaoyl 发表于 2024-12-10 14:49
老师您好，我想问一下您通过Multiwfn软件做出双光子光谱图后，它的纵坐标是摩尔吸收系数，怎么才能转换成 ...

在本论坛上一个问题短时间内多处提问是不允许的，你很快就会受到警告。我没有“通过Multiwfn软件做出双光子光谱图”，不要问我。

作者
Author: Yaoyl 时间: 2024-12-10 14:59

zjxitcc 发表于 2024-12-10 14:55
在本论坛上一个问题短时间内多处提问是不允许的，你很快就会受到警告。我没有“通过Multiwfn软件做出双光 ...

好的，谢谢老师提醒

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