计算化学公社

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

[ORCA] orca复现文献里CASSCF计算(BiO8)^14-报错求助

[复制链接 Copy URL]

2

帖子

0

威望

23

eV
积分
25

Level 2 能力者

求助各位老师!近期在使用orca复现一篇论文时出现报错,原文献中使用计算方法CASSCF,ANO基组,考虑相对论效应,计算模型为(BiO8)^14-,文献中初始Bi-O键长为2.55A,模型图为 下面是我的orca输入文件:
! X2C ANO-RCC-TZP   RIJCOSX TightSCF AutoAux

%pal
      nprocs 24
end
%MaxCore 6000

%casscf
  nel 1      
  norb 3      

  mult 2      
  nroots 3   

   maxiter 200
  ShiftUp 2.0

  rel
       DoSOC true
  end

  PTMethod SC_NEVPT2
  TrafoStep RI
  PTSettings
      F12 true
      D4Step efficient
  end

end

%rel
  PictureChange 2
end

*xyz -14 2
Bi      0.000000    0.000000    0.000000
O       1.472243    1.472243    1.472243
O      -1.472243    1.472243    1.472243
O      -1.472243   -1.472243    1.472243
O       1.472243   -1.472243    1.472243
O       1.472243    1.472243   -1.472243
O      -1.472243    1.472243   -1.472243
O      -1.472243   -1.472243   -1.472243
O       1.472243   -1.472243   -1.472243
*
请教各位老师,我是应该修改活性空间来解决报错吗?


202606121035138765..png (18.34 KB, 下载次数 Times of downloads: 0)

202606121035138765..png

1621

帖子

0

威望

5294

eV
积分
6915

Level 6 (一方通行)

2#
发表于 Post on 2026-6-12 11:21:21 | 只看该作者 Only view this author
用orca计算周期性体系,本身就不是很合适吧?

你试试CP2K
又菜又爱玩

2

帖子

0

威望

23

eV
积分
25

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2026-6-12 14:08:04 | 只看该作者 Only view this author
牧生 发表于 2026-6-12 11:21
用orca计算周期性体系,本身就不是很合适吧?

你试试CP2K

好的谢谢 我去试一下CP2K 另外请教一下这个问题在orca中是没有办法解决的吗 我们课题组要求是使用orca进行研究

4347

帖子

4

威望

9993

eV
积分
14420

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2026-6-12 16:42:22 | 只看该作者 Only view this author
如果原文献CASSCF计算用的结构模型是有限原子的团簇,并且采用Molpro/OpenMolcas做的计算,那么你用ORCA肯定可以复现。当然,即使你可以复现理论计算结果,但与实际体系的实验结果比较 是另一码事,要考虑到体系是否可以简化为团簇模型进行计算,还是说不能如此粗糙 而应当考虑周期性。这是两个独立的问题。
自动做多参考态计算的程序MOKIT

116

帖子

0

威望

1528

eV
积分
1644

Level 5 (御坂)

5#
发表于 Post on 2026-6-12 16:44:26 | 只看该作者 Only view this author
你这个标题不是很合适,建议改一下,这个标题完全不知道你在问什么

如果使用CASSCF的画,这种一般是用EMB算的,也就是中心团簇加基组,团簇外用点电荷势场。点电荷可以用Ewald程序生成,而且做CASSCF一般是基于之前计算的参考轨道的,本身EMB模型就有可能有收敛问题,你还不读参考轨道。

个人更推荐用Molcas做,ORCA也可以实现,点电荷的关键词我记得是point charge。靠近团簇的点电荷是需要用赝势包裹的,但是ORCA没有对应的赝势,偷懒的话可以最直接用SMD。
可以参考:https://molcas.gitlab.io/OpenMol ... xamples/ex-ecp.html以及Sergentu 21年的几篇工作,他的SI有具体怎么做计算的说明。

6万

帖子

99

威望

6万

eV
积分
127652

管理员

公社社长

6#
发表于 Post on 2026-6-15 18:55:15 | 只看该作者 Only view this author
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题清楚、准确反映出帖子具体内容,避免有任何歧义和含糊性,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “orca计算报错求助” 改了,以后务必注意!
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-6-24 19:24 , Processed in 0.236828 second(s), 25 queries , Gzip On.

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