计算化学公社

 找回密码 Forget password
 注册 Register
Views: 118|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[其它量化程序] BDF XIANCI模块一直不收敛

[复制链接 Copy URL]

4

帖子

0

威望

273

eV
积分
277

Level 3 能力者

各位老师好,我正在用BDF的XIANCI模块计算自由基分子的SDSPT2,但程序一直卡在CI循环。把maxiter加到400也不行。请问各位老师有什么办法吗?谢谢!

8683ce34-fb05-4771-a468-d3058c6b5465.png (62.47 KB, 下载次数 Times of downloads: 2)

报错

报错

1cz-sds-1.inp

2.67 KB, 下载次数 Times of downloads: 1

输入文件

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

2#
发表于 Post on 2024-10-18 13:53:24 | 只看该作者 Only view this author
正常来说xianci不应该有这么大幅度的振荡。检查有没有基组线性相关等问题
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

4

帖子

0

威望

273

eV
积分
277

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-10-18 14:41:03 | 只看该作者 Only view this author
本帖最后由 hhhhhjjdh 于 2024-10-18 14:45 编辑
wzkchem5 发表于 2024-10-18 13:53
正常来说xianci不应该有这么大幅度的振荡。检查有没有基组线性相关等问题

谢谢老师回复!基组线性相关是检查MCSCF模块吗,还有一个问题是我用MCSCF模块做casscf计算,但轨道占据数并不是由大到小的(这个问题能说明存在线性相关的问题吗)

5049df67-e862-487d-b840-d3b2bd7504cc.png (453.87 KB, 下载次数 Times of downloads: 1)

5049df67-e862-487d-b840-d3b2bd7504cc.png

689

帖子

21

威望

5019

eV
积分
6128

Level 6 (一方通行)

4#
发表于 Post on 2024-10-18 18:13:31 | 只看该作者 Only view this author
可能的话,建议用2023年以前的版本,没有这种问题,但是输入不兼容,更繁琐一些。

4

帖子

0

威望

273

eV
积分
277

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-10-18 19:14:57 | 只看该作者 Only view this author
beefly 发表于 2024-10-18 18:13
可能的话,建议用2023年以前的版本,没有这种问题,但是输入不兼容,更繁琐一些。

谢谢老师!我去试试

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

6#
发表于 Post on 2024-10-20 13:24:09 | 只看该作者 Only view this author
hhhhhjjdh 发表于 2024-10-18 07:41
谢谢老师回复!基组线性相关是检查MCSCF模块吗,还有一个问题是我用MCSCF模块做casscf计算,但轨道占据数 ...

检查重叠矩阵的最小本征值
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

689

帖子

21

威望

5019

eV
积分
6128

Level 6 (一方通行)

7#
发表于 Post on 2024-10-20 17:10:28 | 只看该作者 Only view this author
本帖最后由 beefly 于 2024-10-20 17:14 编辑
wzkchem5 发表于 2024-10-20 13:24
检查重叠矩阵的最小本征值

2023年以前的版本中,mcscf模块优化好的CI因子直接被xian-ci读取,这类似于molcas的做法。所以只要mcscf收敛了,xian-ci计算就没问题。
从2023年开始,xian-ci会重新算一遍mcscf参考态的能量和CI因子,这是molpro的做法。由于算法存在差异,不如老代码稳定,导致即便前面的mcscf已经收敛了,xian-ci的参考态仍然可能不收敛。

我认为最保险的做法应该是二者都支持,慢慢地过渡到后者。

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

8#
发表于 Post on 2024-10-20 18:22:42 | 只看该作者 Only view this author
beefly 发表于 2024-10-20 10:10
2023年以前的版本中,mcscf模块优化好的CI因子直接被xian-ci读取,这类似于molcas的做法。所以只要mcscf ...

明白了,这么看来的话确实应该在目前的代码里提供2023年以前的做法的选项
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 11:03 , Processed in 0.170866 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list