计算化学公社

标题: BDF XIANCI模块一直不收敛 [打印本页]

作者
Author: hhhhhjjdh 时间: 2024-10-18 11:02
标题: BDF XIANCI模块一直不收敛
各位老师好，我正在用BDF的XIANCI模块计算自由基分子的SDSPT2，但程序一直卡在CI循环。把maxiter加到400也不行。请问各位老师有什么办法吗？谢谢！

作者
Author: wzkchem5 时间: 2024-10-18 13:53
正常来说xianci不应该有这么大幅度的振荡。检查有没有基组线性相关等问题

作者
Author: hhhhhjjdh 时间: 2024-10-18 14:41
本帖最后由 hhhhhjjdh 于 2024-10-18 14:45 编辑

wzkchem5 发表于 2024-10-18 13:53
正常来说xianci不应该有这么大幅度的振荡。检查有没有基组线性相关等问题

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

作者
Author: beefly 时间: 2024-10-18 18:13
可能的话，建议用2023年以前的版本，没有这种问题，但是输入不兼容，更繁琐一些。

作者
Author: hhhhhjjdh 时间: 2024-10-18 19:14

beefly 发表于 2024-10-18 18:13
可能的话，建议用2023年以前的版本，没有这种问题，但是输入不兼容，更繁琐一些。

谢谢老师！我去试试

作者
Author: wzkchem5 时间: 2024-10-20 13:24

hhhhhjjdh 发表于 2024-10-18 07:41
谢谢老师回复!基组线性相关是检查MCSCF模块吗，还有一个问题是我用MCSCF模块做casscf计算，但轨道占据数 ...

检查重叠矩阵的最小本征值

作者
Author: beefly 时间: 2024-10-20 17:10
本帖最后由 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的参考态仍然可能不收敛。

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

作者
Author: wzkchem5 时间: 2024-10-20 18:22

beefly 发表于 2024-10-20 10:10
2023年以前的版本中，mcscf模块优化好的CI因子直接被xian-ci读取，这类似于molcas的做法。所以只要mcscf ...

明白了，这么看来的话确实应该在目前的代码里提供2023年以前的做法的选项

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