计算化学公社

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

[PSI4] 求助 PSI4能量分解SCFConvergenceError

[复制链接 Copy URL]

31

帖子

0

威望

1123

eV
积分
1154

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 hmlyz 于 2021-4-12 15:06 编辑

我在做能量分解时使用了guess GW,Hsoscf true和soscf_max_iter 30后仍然显示: Could not converge SCF iterations in 100 iterations.请问还有什么办法吗?(对不起我在官网网页中手册没有找到解决办法)Did not take a SOSCF step, using normal convergence methods
   @DF-UHF iter  90: -1967.65319450610241    2.62976e+00   5.79728e-03 DIIS
   @DF-UHF iter  91: -1986.37494718918242   -1.87218e+01   2.40559e-04 SOSCF, nmicro=30
   @DF-UHF iter  92: -1986.17125761674106    2.03690e-01   1.01378e-03 SOSCF, nmicro=30
   @DF-UHF iter  93: -1984.95188761410827    1.21937e+00   9.23616e-04 SOSCF, nmicro=30
   @DF-UHF iter  94: -1985.00052609940667   -4.86385e-02   7.37083e-04 SOSCF, nmicro=30
   @DF-UHF iter  95: -1984.80110780120413    1.99418e-01   9.02371e-04 SOSCF, nmicro=30
   @DF-UHF iter  96: -1984.90624767214831   -1.05140e-01   7.96589e-04 SOSCF, nmicro=30
   @DF-UHF iter  97: -1983.82907595182496    1.07717e+00   1.48462e-03 SOSCF, nmicro=30
   @DF-UHF iter  98: -1981.03686986099819    2.79221e+00   2.88602e-03 SOSCF, nmicro=30
   @DF-UHF iter  99: -1976.68592011800229    4.35095e+00   4.45237e-03 SOSCF, nmicro=30
   @DF-UHF iter 100: -1974.19661636077672    2.48930e+00   4.17423e-03 SOSCF, nmicro=30

PsiException: Could not converge SCF iterations in 100 iterations.

  Failed to converge.

Traceback (most recent call last):
  File "/home/haorui/psi4conda/bin/psi4", line 287, in <module>
    exec(content)
  File "<string>", line 70, in <module>
  File "/home/haorui/psi4conda/lib//python3.6/site-packages/psi4/driver/driver.py", line 556, in energy
    wfn = procedures['energy'][lowername](lowername, molecule=molecule, **kwargs)
  File "/home/haorui/psi4conda/lib//python3.6/site-packages/psi4/driver/procrouting/proc.py", line 3325, in run_sapt
    monomerA_wfn = scf_helper('RHF', molecule=monomerA, **kwargs)
  File "/home/haorui/psi4conda/lib//python3.6/site-packages/psi4/driver/procrouting/proc.py", line 1363, in scf_helper
    e_scf = scf_wfn.compute_energy()
  File "/home/haorui/psi4conda/lib//python3.6/site-packages/psi4/driver/procrouting/scf_proc/scf_iterator.py", line 94, in scf_compute_energy
    raise e
  File "/home/haorui/psi4conda/lib//python3.6/site-packages/psi4/driver/procrouting/scf_proc/scf_iterator.py", line 87, in scf_compute_energy
    self.iterations()
  File "/home/haorui/psi4conda/lib//python3.6/site-packages/psi4/driver/procrouting/scf_proc/scf_iterator.py", line 447, in scf_iterate
    raise SCFConvergenceError("""SCF iterations""", self.iteration_, self, Ediff, Dnorm)

psi4.driver.p4util.exceptions.SCFConvergenceError: Could not converge SCF iterations in 100 iterations.






6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2021-4-12 08:13:39 | 只看该作者 Only view this author
说清楚什么体系、用的什么关键词、之前的迭代收敛情况。没这些信息没人能帮你
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

31

帖子

0

威望

1123

eV
积分
1154

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-4-12 09:12:39 | 只看该作者 Only view this author
sobereva 发表于 2021-4-12 08:13
说清楚什么体系、用的什么关键词、之前的迭代收敛情况。没这些信息没人能帮你

好的卢老师,非常抱歉。已经上传了输入输出文件

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

4#
发表于 Post on 2021-4-12 09:15:35 | 只看该作者 Only view this author
C-Fe之间是化学键作用,原理上就不适合用SAPT
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

31

帖子

0

威望

1123

eV
积分
1154

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2021-4-12 09:20:37 | 只看该作者 Only view this author
sobereva 发表于 2021-4-12 09:15
C-Fe之间是化学键作用,原理上就不适合用SAPT

是的卢老师,我看您的文章发现了这个问题,只是在这里尝试做了一下。请问还有其他办法做这个能量分解吗?

1242

帖子

1

威望

3672

eV
积分
4934

Level 6 (一方通行)

6#
发表于 Post on 2021-4-12 11:56:07 | 只看该作者 Only view this author
hmlyz 发表于 2021-4-12 09:20
是的卢老师,我看您的文章发现了这个问题,只是在这里尝试做了一下。请问还有其他办法做这个能量分解吗?

用gamess的EDA

4104

帖子

4

威望

8871

eV
积分
13055

Level 6 (一方通行)

MOKIT开发者

7#
发表于 Post on 2021-4-12 11:58:58 | 只看该作者 Only view this author
hmlyz 发表于 2021-4-12 09:20
是的卢老师,我看您的文章发现了这个问题,只是在这里尝试做了一下。请问还有其他办法做这个能量分解吗?

向厦大苏培峰老师发邮件要GKS-EDA
自动做多参考态计算的程序MOKIT

31

帖子

0

威望

1123

eV
积分
1154

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-4-12 15:06:57 | 只看该作者 Only view this author

好的,非常感谢

31

帖子

0

威望

1123

eV
积分
1154

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2021-4-12 15:07:47 | 只看该作者 Only view this author
zjxitcc 发表于 2021-4-12 11:58
向厦大苏培峰老师发邮件要GKS-EDA

好的,非常感谢 我先查查

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 16:16 , Processed in 0.181694 second(s), 21 queries , Gzip On.

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