计算化学公社

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

[ORCA] 求助:socme计算与文献差很大

[复制链接 Copy URL]

11

帖子

0

威望

234

eV
积分
245

Level 3 能力者

各位大佬,
下面是我的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


202303221640528046..png (43.18 KB, 下载次数 Times of downloads: 6)

图2

图2

202303221639001012..png (143.38 KB, 下载次数 Times of downloads: 5)

图1

图1

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

2#
发表于 Post on 2023-3-22 17:05:57 | 只看该作者 Only view this author
gridx4报错是因为gridx4是ORCA 4.x.x及以前才支持的关键字,新版orca用的格点变了,需要用defgrid3等关键字增加格点
是不是你用的构象和文献不一样?
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
本团队长期招收研究生,有意者可私信联系

245

帖子

0

威望

780

eV
积分
1025

Level 4 (黑子)

3#
发表于 Post on 2023-3-22 18:45:13 | 只看该作者 Only view this author
要不把TDA开了看看?完整复现一下,TDA和TDDFT方法可能也会引入一些不同吧。

11

帖子

0

威望

234

eV
积分
245

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-3-22 19:29:28 | 只看该作者 Only view this author
本帖最后由 风风光风光 于 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

paper-HOMOandLUMO.png (37.81 KB, 下载次数 Times of downloads: 4)

paperHOMO\LUMO

paperHOMO\LUMO

LUMO.png (53.62 KB, 下载次数 Times of downloads: 5)

LUMO

LUMO

HOMO.png (55.18 KB, 下载次数 Times of downloads: 5)

HOMO

HOMO

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

5#
发表于 Post on 2023-3-22 20:14:08 | 只看该作者 Only view this author
风风光风光 发表于 2023-3-22 12:29
感谢感谢,我用文献的泛函和基组优化完的构象不太一样,LUMO、HOMO图对不上。

检查你的isovalue值是否和文献一致
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

6#
发表于 Post on 2023-3-23 06:50:46 | 只看该作者 Only view this author
风风光风光 发表于 2023-3-22 19:29
感谢感谢,我用文献的泛函和基组优化完的构象不太一样,LUMO、HOMO图对不上。

gview默认0.05等值面数值,Multiwfn默认0.02。统一了才能比较。文献如果没明确提,一般用的是程序默认的等值面数值
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

144

帖子

0

威望

2533

eV
积分
2677

Level 5 (御坂)

7#
发表于 Post on 2023-3-23 19:55:59 | 只看该作者 Only view this author
'我用文献的泛函和基组优化完的构象不太一样,'
这个十分影响SOC的计算值,惊不惊喜,意不意外

11

帖子

0

威望

234

eV
积分
245

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2023-3-23 21:16:50 | 只看该作者 Only view this author
本帖最后由 风风光风光 于 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老师有什么好的建议或意见吗?

11

帖子

0

威望

234

eV
积分
245

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-3-23 21:18:45 | 只看该作者 Only view this author
wzkchem5 发表于 2023-3-22 20:14
检查你的isovalue值是否和文献一致

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

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

10#
发表于 Post on 2023-3-23 22:22:04 | 只看该作者 Only view this author
风风光风光 发表于 2023-3-23 14:16
谢谢sob大佬的回答,我调整了等值面数值,现在能把HOMO和LUMO图对上了。我还想请教您一下,
1. 为什么我 ...

1. 改的时候没写错就行,本来orca的输入格式就是为用户直接手写设计的,而不是为Multiwfn自动产生来设计的。主要注意M06-2X在orca里不要加横杠
2. 这种问题自己看orca说明书detailed documentation里的basis set相关章节,有一个大列表,每个基组支持哪些元素
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

11#
发表于 Post on 2023-3-24 00:31:03 | 只看该作者 Only view this author
风风光风光 发表于 2023-3-23 21:16
谢谢sob大佬的回答,我调整了等值面数值,现在能把HOMO和LUMO图对上了。我还想请教您一下,
1. 为什么我 ...

对于ORCA用户来说用M06-2X的意义本来就不大,所以下文的功能里没有给M06-2X留位置。M06-2X擅长计算有机体系热力学和弱相互作用,但这方面都比不过ORCA支持的wB97M-V
详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://sobereva.com/490http://bbs.keinsci.com/thread-13560-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 10:03 , Processed in 1.361113 second(s), 23 queries , Gzip On.

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