计算化学公社

标题: 求助：socme计算与文献差很大 [打印本页]

作者
Author: 风风光风光 时间: 2023-3-22 16:43
标题: 求助：socme计算与文献差很大
各位大佬，
下面是我的inp文件写法
! M062X ZORA ZORA-def2-TZVP SARC/J RI-SOMF(1X) cpcm(Acetonitrile) tightSCF miniprint nopop pal4
%basis
NewGTO I "SARC-ZORA-TZVP" end
end
%maxcore 1000
%pal nprocs 4 end
%tddft
nroots 3
tda false
dosoc true
printlevel 3
end
* xyz 0 1

图1是文献中的计算方法和计算结果，图2是我用上方inp文件计算的计算结果

另外inp中写gridx4运行之后会报错，删掉之后不报错了，这是为什么，版本是ORCA5.0.3

作者
Author: wzkchem5 时间: 2023-3-22 17:05
gridx4报错是因为gridx4是ORCA 4.x.x及以前才支持的关键字，新版orca用的格点变了，需要用defgrid3等关键字增加格点
是不是你用的构象和文献不一样？

作者
Author: KSeGaSn 时间: 2023-3-22 18:45
要不把TDA开了看看？完整复现一下，TDA和TDDFT方法可能也会引入一些不同吧。

作者
Author: 风风光风光 时间: 2023-3-22 19:29
本帖最后由风风光风光于 2023-3-22 19:30 编辑

wzkchem5 发表于 2023-3-22 17:05
gridx4报错是因为gridx4是ORCA 4.x.x及以前才支持的关键字，新版orca用的格点变了，需要用defgrid3等关键字 ...

感谢感谢，我用文献的泛函和基组优化完的构象不太一样，LUMO、HOMO图对不上。C:\Users\12971\Desktop\LUMO.pngC:\Users\12971\Desktop\HOMO.pngC:\Users\12971\Desktop\paper-HOMOandLUMO.png

作者
Author: wzkchem5 时间: 2023-3-22 20:14

风风光风光发表于 2023-3-22 12:29
感谢感谢，我用文献的泛函和基组优化完的构象不太一样，LUMO、HOMO图对不上。

检查你的isovalue值是否和文献一致

作者
Author: sobereva 时间: 2023-3-23 06:50

风风光风光发表于 2023-3-22 19:29
感谢感谢，我用文献的泛函和基组优化完的构象不太一样，LUMO、HOMO图对不上。

gview默认0.05等值面数值，Multiwfn默认0.02。统一了才能比较。文献如果没明确提，一般用的是程序默认的等值面数值

作者
Author: 风起~ 时间: 2023-3-23 19:55
'我用文献的泛函和基组优化完的构象不太一样，'
这个十分影响SOC的计算值，惊不惊喜，意不意外

作者
Author: 风风光风光 时间: 2023-3-23 21:16
本帖最后由风风光风光于 2023-3-23 21:30 编辑

sobereva 发表于 2023-3-23 06:50
gview默认0.05等值面数值，Multiwfn默认0.02。统一了才能比较。文献如果没明确提，一般用的是程序默认的 ...

谢谢sob大佬的回答，我调整了等值面数值，现在能把HOMO和LUMO图对上了。我还想请教您一下，
1. 为什么我在使用multiwfn转ORCA的输入文件时，没有找到M06-2X的选项，所以这个泛函是我在文件里手动改的，我算不对会是这个原因造成的吗？
2. 我对碘使用了全电子基组，书写方法如下
%basis
NewGTO I "SARC-ZORA-TZVP" end
end
我想问这个基组对碘有效吗？我怀疑可能是这个基组选择的问题，请问sob老师有什么好的建议或意见吗？

作者
Author: 风风光风光 时间: 2023-3-23 21:18

wzkchem5 发表于 2023-3-22 20:14
检查你的isovalue值是否和文献一致

谢谢大佬的指导，我查看了isovalue值，能把HOMO和LUMO构型图对上。

作者
Author: wzkchem5 时间: 2023-3-23 22:22

风风光风光发表于 2023-3-23 14:16
谢谢sob大佬的回答，我调整了等值面数值，现在能把HOMO和LUMO图对上了。我还想请教您一下，
1. 为什么我 ...

1. 改的时候没写错就行，本来orca的输入格式就是为用户直接手写设计的，而不是为Multiwfn自动产生来设计的。主要注意M06-2X在orca里不要加横杠
2. 这种问题自己看orca说明书detailed documentation里的basis set相关章节，有一个大列表，每个基组支持哪些元素

作者
Author: sobereva 时间: 2023-3-24 00:31

风风光风光发表于 2023-3-23 21:16
谢谢sob大佬的回答，我调整了等值面数值，现在能把HOMO和LUMO图对上了。我还想请教您一下，
1. 为什么我 ...

对于ORCA用户来说用M06-2X的意义本来就不大，所以下文的功能里没有给M06-2X留位置。M06-2X擅长计算有机体系热力学和弱相互作用，但这方面都比不过ORCA支持的wB97M-V
详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://sobereva.com/490（http://bbs.keinsci.com/thread-13560-1-1.html）

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