计算化学公社

标题: 求助DFT计算超精细耦合常数遇到原子序数大的元素怎么解决 [打印本页]

作者
Author: Lucienmlj 时间: 2021-3-13 10:01
标题: 求助DFT计算超精细耦合常数遇到原子序数大的元素怎么解决
求助：我在计算超精细耦合常数时，使用EPR-II基组，但是我的分子中含有Br原子，程序一直报错“Atomic number out of range in EPR2”，请问各位大神，我该怎么处理？
分子式如下：

作者
Author: biogon 时间: 2021-3-13 10:39
换个基组

作者
Author: 北大-陶豫 时间: 2021-3-13 10:44
基组并不是对所有原子都有定义的，你选用基组要让它对当前所有原子都有定义才行。

详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

谈谈量子化学中基组的选择
http://sobereva.com/336

作者
Author: Lucienmlj 时间: 2021-3-13 17:20

biogon 发表于 2021-3-13 10:39
换个基组

谢谢，我试试

作者
Author: Lucienmlj 时间: 2021-3-13 17:20

北大-陶豫发表于 2021-3-13 10:44
基组并不是对所有原子都有定义的，你选用基组要让它对当前所有原子都有定义才行。

详解Gaussian中混合基 ...

感谢，我试下

作者
Author: wzkchem5 时间: 2021-3-13 22:30
另外注意相对论效应的处理，尤其假如你关心Br的超精细耦合的话

作者
Author: Lucienmlj 时间: 2021-3-13 22:34

wzkchem5 发表于 2021-3-13 22:30
另外注意相对论效应的处理，尤其假如你关心Br的超精细耦合的话

请问我该怎么算呢，还是想算Br的

作者
Author: wzkchem5 时间: 2021-3-14 00:31

Lucienmlj 发表于 2021-3-13 22:34
请问我该怎么算呢，还是想算Br的

你用的是ORCA吗？如果用ORCA的话，建议用ZORA显式考虑相对论效应，同时用一个相对论基组。适合EPR计算的相对论基组我还真不太清楚，但是我猜ZORA-def2-TZVPP应该不会太差

作者
Author: sobereva 时间: 2021-3-14 02:37
可以用UGBS，完全未收缩，同时满足相对论计算和算超精细耦合需要的紧核s函数。
但这基组缺点就是收敛难，对积分格点精度要求极高。

要么你在给相对论全电子计算优化的基组基础之上手动加一些紧核s函数增加对原子核处自旋密度的描述。可以从比如pcH等专门给超精细耦合计算的基组那里借来

以上只是假定只考虑Fermi-Contact部分。

作者
Author: Lucienmlj 时间: 2021-3-15 09:58

wzkchem5 发表于 2021-3-14 00:31
你用的是ORCA吗？如果用ORCA的话，建议用ZORA显式考虑相对论效应，同时用一个相对论基组。适合EPR计算的 ...

我暂时按文献，用的高斯

作者
Author: Lucienmlj 时间: 2021-3-15 09:58

sobereva 发表于 2021-3-14 02:37
可以用UGBS，完全未收缩，同时满足相对论计算和算超精细耦合需要的紧核s函数。
但这基组缺点就是收敛难， ...

感谢老师，我这就去试试

作者
Author: 小李子赛 时间: 2022-4-27 16:16
求问输入文件怎么写得，刚接触这个，不太知道

作者
Author: QuantumicGuy 时间: 2022-4-27 18:19

小李子赛发表于 2022-4-27 16:16
求问输入文件怎么写得，刚接触这个，不太知道

https://easyspin.org/easyspin/documentation/orca.html

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