计算化学公社

标题: 求助用CASSCF计算Yb3+自由离子时roots能量异常 [打印本页]

作者
Author: yinhaoyu98 时间: 2020-8-31 17:19
标题: 求助用CASSCF计算Yb3+自由离子时roots能量异常
本帖最后由 yinhaoyu98 于 2020-8-31 17:21 编辑

我的输入文本如图所示：
(, 下载次数 Times of downloads: 22)
遇到的问题是这样的：（输出文本的"INITIAL CI STATE CHECK" 和 "CASSCF STATES" 部分得到的能量出奇的大）
(, 下载次数 Times of downloads: 24) (, 下载次数 Times of downloads: 26)
检查后，我发现是"orbital energies"部分有异常：
(, 下载次数 Times of downloads: 28)
从27开始应该是4f轨道，能量应该简并，但是图上看起来大小不一。我还尝试了其它基组，"SARC"系列的基组都有这种问题，但是如果我把基组改为def2-TZVP 或者 def2-QZVP,所得到的结果理论上是说的通的。
请问我该怎么做，使用SARC系列的基组来计算该自由离子？因为我是新手，所以希望各位老师给予详细的解答，谢谢大家！

作者
Author: beefly 时间: 2020-8-31 18:18
casscf的初猜不好，例如hf轨道就是不简并的。建议用闭壳层离子产生hf轨道

作者
Author: yinhaoyu98 时间: 2020-8-31 19:13

beefly 发表于 2020-8-31 18:18
casscf的初猜不好，例如hf轨道就是不简并的。建议用闭壳层离子产生hf轨道

请问我该怎么做? 用！RHF 先做出.gbw再调用吗？

作者
Author: zjxitcc 时间: 2020-8-31 20:43
本帖最后由 zjxitcc 于 2020-8-31 20:46 编辑

我记得你在群里问过Ce3+和Yb3+怎么算吧。。我都回复得挺详细的了，咋就不看咧。。。

随意举例，请参考计算步骤，至于哪个步骤用哪个软件算，爱用啥用啥，并不唯一。
Step 1. Yb原子的RHF计算。
Step 2. Yb3+的SA-CASSCF计算。
这两步的高斯输入文件（内存核数自己设）如下

%chk=Yb_all.chk
#p RHF/gen nosymm int(nobasistransform,DKH2)
title
0 1
Yb 0.0 0.0 0.0
[空行]
去BSE上复制SARC-DKH2基组贴于此
[空行]
[空行]
--Link1--
%chk=Yb_all.chk
#p CAS(13,7,StateAverage,FullDiag,nroot=7)/chkbasis nosymm int(nobasistransform,DKH2) guess=read geom=check
title
3 2
0.14285714
0.14285714
0.14285714
0.14285714
0.14285714
0.14285714
0.14285714
[空行]

复制代码

Step 3. 上述两步完成后，将高斯chk文件转化为fchk文件，用MOKIT中的fch2mkl小程序（https://gitlab.com/jxzou/mokit）从fchk文件生成ORCA的inp和mkl文件，用orca_2mkl将mkl文件转化为gbw文件。修改.inp文件前几行如下（内存核数自己设）

! SARC-DKH-TZVP
%rel
method DKH
order 2
end
%casscf
nel 13
norb 7
mult 2
nroots 7
...

复制代码

SA-CASSCF将在5圈以内迅速收敛，7个根简并。你的帖子里，忘记加上DKH2部分关键词了。

作者
Author: yinhaoyu98 时间: 2020-8-31 20:53

zjxitcc 发表于 2020-8-31 20:43
我记得你在群里问过Ce3+和Yb3+怎么算吧。。我都回复得挺详细的了，咋就不看咧。。。

随意举例，请参考计 ...

我刚刚开始这方面的学习，我试过群里的方法，但是自己得不到正常的结果

我试试老师您说的方法，谢谢哈！

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