计算化学公社

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

[新手求助] orca计算报错求助

[复制链接 Copy URL]

2

帖子

0

威望

15

eV
积分
17

Level 1 能力者

求助各位老师!近期在使用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

威望

5274

eV
积分
6895

Level 6 (一方通行)

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

你试试CP2K
又菜又爱玩

2

帖子

0

威望

15

eV
积分
17

Level 1 能力者

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

你试试CP2K

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

4337

帖子

4

威望

9943

eV
积分
14360

Level 6 (一方通行)

MOKIT开发者

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

115

帖子

0

威望

1508

eV
积分
1623

Level 5 (御坂)

5#
发表于 Post on 前天 16:44 | 只看该作者 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有具体怎么做计算的说明。

本版积分规则 Credits rule

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

GMT+8, 2026-6-14 08:34 , Processed in 0.302582 second(s), 24 queries , Gzip On.

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