计算化学公社

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

[ORCA] 求助CASSCF计算Ho3+自由离子的步骤

[复制链接 Copy URL]

20

帖子

0

威望

345

eV
积分
365

Level 3 能力者

我之前在老师们的提点下完成了用ORCA CASSCF计算Yb3+自由离子的任务,由于Yb2+和Yb原子都是闭壳的,因此可以先计算它们的轨道再运用到Yb3+的计算中。
但是,现在我需要计算Ho3+自由离子的话,直接计算Ho3+离子的轨道能量得出的结果并不正确,请问该如何着手这个问题?

4289

帖子

4

威望

9550

eV
积分
13919

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2020-9-9 10:11:14 | 只看该作者 Only view this author
Ce3+, Yb3+,Ho3+这些不都是同一套技巧么。。。咋一个原子问一次

你可以试试从Ho-或Ho3-负离子出发
自动做多参考态计算的程序MOKIT

20

帖子

0

威望

345

eV
积分
365

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2020-9-9 10:12:55 | 只看该作者 Only view this author
zjxitcc 发表于 2020-9-9 10:11
Ce3+, Yb3+,Ho3+这些不都是同一套技巧么。。。咋一个原子问一次

你可以试试从Ho-或Ho3-负离子出发

可以可以,老师之前讲的方法很好,我现在还做不到举一反三

20

帖子

0

威望

345

eV
积分
365

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2020-9-12 19:06:22 | 只看该作者 Only view this author
本帖最后由 yinhaoyu98 于 2020-9-12 19:15 编辑
zjxitcc 发表于 2020-9-9 10:11
Ce3+, Yb3+,Ho3+这些不都是同一套技巧么。。。咋一个原子问一次

你可以试试从Ho-或Ho3-负离子出发

可是Ho-和Ho3-离子不存在啊,应该不能这么算吧?
算所有镧系金属3+离子都是要这样调用闭壳的轨道嘛? 求指点。

4289

帖子

4

威望

9550

eV
积分
13919

Level 6 (一方通行)

MOKIT开发者

5#
发表于 Post on 2020-9-12 19:25:32 | 只看该作者 Only view this author
yinhaoyu98 发表于 2020-9-12 19:06
可是Ho-和Ho3-离子不存在啊,应该不能这么算吧?
算所有镧系金属3+离子都是要这样调用闭壳的轨道嘛? 求 ...

你说的不存在是哪种不存在?

这只是计算过程中的一步而已,造出个合理的初始轨道,不是最终的结果。也有其他办法,但是说起来都比较复杂,门槛更高。
自动做多参考态计算的程序MOKIT

20

帖子

0

威望

345

eV
积分
365

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2020-9-12 19:30:52 | 只看该作者 Only view this author
zjxitcc 发表于 2020-9-12 19:25
你说的不存在是哪种不存在?

这只是计算过程中的一步而已,造出个合理的初始轨道,不是最终的结果。也 ...

那么像Pr,Nd这种元素也能造出来嘛?

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

7#
发表于 Post on 2020-9-12 21:23:27 | 只看该作者 Only view this author
可以用FOD轨道作为初猜(参见ORCA说明书的FOD部分),这样即便你的初猜用的不是f0、f7或f14的离子,也能保持球对称性。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

20

帖子

0

威望

345

eV
积分
365

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2020-9-12 21:31:17 | 只看该作者 Only view this author
wzkchem5 发表于 2020-9-12 21:23
可以用FOD轨道作为初猜(参见ORCA说明书的FOD部分),这样即便你的初猜用的不是f0、f7或f14的离子,也能保 ...

好的,我试试,谢谢!

20

帖子

0

威望

345

eV
积分
365

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2020-9-15 18:34:45 | 只看该作者 Only view this author
本帖最后由 yinhaoyu98 于 2020-9-16 09:14 编辑

现在遇到的问题是算得的数据和nist上的标准数据完全不匹配,是我的input有问题吗?Pr3+离子

  1. ! Sapporo-DKH3-QZP-2012 TightSCF conv

  2. %pal nprocs 10 end

  3. %casscf  
  4.         nel 2      
  5.         norb 7
  6.         mult 3,1
  7.         nroots 21,28

  8.         PTMethod SC_NEVPT2

  9.         rel
  10.         dosoc true
  11.         end
  12. end
  13. *xyz 3 3
  14. Pr 0 0 0
  15. *   

复制代码


得到的SOC数据:                                                                     


NIST的数据:
                                                                                                            


本版积分规则 Credits rule

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

GMT+8, 2026-2-24 06:35 , Processed in 0.193090 second(s), 23 queries , Gzip On.

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