计算化学公社

标题: 请教orca基组和赝势语法格式，在orca里面使用自定义基组和赝势显示无法识别 [打印本页]

作者
Author: morningya 时间: 2023-7-17 19:22
标题: 请教orca基组和赝势语法格式，在orca里面使用自定义基组和赝势显示无法识别
各位老师好，新手求助，我想用PBEO的泛函，给C、N、F、H使用6-31*基组，Ir使用LanLDZ的基组和赝势，输入文件格式如下：

! PBE0 DKH2 DKH-def2-SVP SARC/J RIJCOSX tightSCF miniprint
%pal nprocs 1 end
%maxcore 10000
%basis
NewGTO H "6-31G*" end
NewGTO C "6-31G*" end
NewGTO N "6-31G*" end
NewGTO F "6-31G*" end
NewGTO Ir "LanLDZ" end
NewECP Ir "LanLDZ" end
end

翻了一下手册和之前的帖子，意思是基组和赝势要手动定义？但是运行显示基组无法识别，想问一下老师们应该怎么修改，谢谢！

作者
Author: wzkchem5 时间: 2023-7-17 19:43
因为那叫LANL2DZ，少写一个2程序当然不认

作者
Author: zjxitcc 时间: 2023-7-17 19:47
（1）DKH2哈密顿要搭配相对论全电子基组（如DKH-def2-SVP就是一种），不可以给其他/某些原子使用普通基组或赝势，因此无需探究此问题。
（2）如果你不想用DKH2哈密顿，却写了DKH2 DKH-def2-SVP SARC/J，说明不知道自己在算啥、想算啥，需要说出自己的详细计算目的，让大家告诉你正确的写法。
（3）是PBE0泛函，不是PBEO的泛函。

作者
Author: morningya 时间: 2023-7-17 20:34

wzkchem5 发表于 2023-7-17 19:43
因为那叫LANL2DZ，少写一个2程序当然不认

啊，谢谢老师指出，然后我又改成了

%basis
NewGTO H "6-31G*" end
NewGTO C "6-31G*" end
NewGTO N "6-31G*" end
NewGTO F "6-31G*" end
NewGTO Ir "LANL2DZ" end
NewECP Ir "LANL2DZ" end
end

显示ECP识别不了。。

作者
Author: morningya 时间: 2023-7-17 20:40

zjxitcc 发表于 2023-7-17 19:47
（1）DKH2哈密顿要搭配相对论全电子基组（如DKH-def2-SVP就是一种），不可以给其他/某些原子使用普通基组或 ...

我是想算铱配合物激发态的旋轨耦合，但是之前没用过orca，基态结构是用高斯优化的，用的泛函就是PB0，混合基组6-31g*、LANL2DZ，LANL2DZ是给金属的基组和赝势。想着用orca应该需要基组一致。
现在改成了：

! PBE0 D3BJ def2-SVP def2/J TightSCF Grid4 FinalGrid5
%pal nprocs 1 end
%maxcore 10000
%basis
NewGTO H "6-31G*" end
NewGTO C "6-31G*" end
NewGTO N "6-31G*" end
NewGTO F "6-31G*" end
NewGTO Ir "LANL2DZ" end
NewECP Ir "LANL2DZ" end
end
但是还是不行

作者
Author: zjxitcc 时间: 2023-7-17 21:08
本帖最后由 zjxitcc 于 2023-7-17 21:13 编辑

morningya 发表于 2023-7-17 20:40
我是想算铱配合物激发态的旋轨耦合，但是之前没用过orca，基态结构是用高斯优化的，用的泛函就是PB0，混 ...

（1）关于自定义基组和赝势
无需如此复杂。直接从你之前的高斯fch文件一键生成ORCA inp和mkl文件即可，既能避免繁琐易错的输入文件写法，还能保证基组赝势两个程序完全一致，还能传轨道加速收敛。假设我的fch文件名是test.fch，运行
fch2mkl test.fch
即可生成test_o.inp和test_o.mkl文件。fch2mkl是免费、开源的MOKIT程序中的一个小程序，直接下载就可使用。你想要什么基组赝势、想怎么自定义，在gjf折腾一下，甚至从零快速产生一个fch文件也可以。然后fch2mkl一下就行了，不用把时间花在研究诸多程序的输入文件格式上。
（2）你似乎在使用ORCA 4，这已经老掉牙了，强烈建议立即更新至ORCA 5.0.4版本。
（3）给过渡金属用LANL2DZ精度是比较低的，不太符合现今需求，用LANL2TZ(f)还差不多。

最后，我不知道ORCA算SOC是否支持有赝势，你最后如果卡在这里，再来问。

作者
Author: hebrewsnabla 时间: 2023-7-17 21:13
本帖最后由 hebrewsnabla 于 2023-7-17 21:17 编辑

本来就没有叫做Lanl2dz的赝势，Lanl2是有的。

（别的程序里面可以用Lanl2dz作为赝势名是因为识别比较宽松）

作者
Author: LittlePupil 时间: 2023-7-17 23:46
4.x版本不清楚，但至少ORCA 5.0.4里Lanl2赝势应该写作“HayWadt”。

作者
Author: hebrewsnabla 时间: 2023-7-18 00:48

LittlePupil 发表于 2023-7-17 23:46
4.x版本不清楚，但至少ORCA 5.0.4里Lanl2赝势应该写作“HayWadt”。

对于Ir来说写Lanl2和HayWadt应该是一样的。

作者
Author: LittlePupil 时间: 2023-7-18 21:22
本帖最后由 LittlePupil 于 2023-7-18 21:30 编辑

hebrewsnabla 发表于 2023-7-18 00:48
对于Ir来说写Lanl2和HayWadt应该是一样的。

Yes, indeed.
只是输出的引用信息会有微小的差别。

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