计算化学公社

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

[新手求助] oniom计算激发态及背景电荷求助

[复制链接 Copy URL]

53

帖子

0

威望

739

eV
积分
792

Level 4 (黑子)

初次使用oniom算晶体激发,重复别人的工作,因为文献的计算方法说的太简单,重复的数据一直对不上,所以打算自己改进做。有几个问题望老师们解答:1. 优化时,使用关键词# opt freq=modelmodes oniom(b3lyp/6-311g(d):uff=qeq)=embedcharge geom=connectivity,两个相同的输入文件放进G09 A01和G16 A03的计算结果不一样。一个是计算过程中,G16的MM层的Br,马利肯电荷值都为0,总能量也不一样。因为之后想用G16在High层加D3校正,所以想弄清楚原因。
2. 计算激发态时,关键词部分为 # oniom(td=(50-50,nstates=7)/b3lyp/6-311g(d):uff=qeq)=embedcharge geom=connectivity。oniom里面没有QM的轨道波函数信息,如果想看相关的激发轨道,是否必须做NTO等轨道分析。
3. 2中的关键词是没法做NTO,之前看社长回答的帖子,要用背景电荷来代替MM区域的电子后照常做TDDFT就可以了,这一过程是否还是用oniom来做。可不可以将1中优化好的High层拿出来做单点,然后转fchk后放到Multiwfn里做RESP,之后再算激发态。另外,加背景电荷时,要在关键词部分加charge,具体加在哪一块呢?说明书的oniom没有看到这个关键词。
4. 如果从最开始的优化我就想MM层用背景电荷,是否可以将单个分子拿出来,冻结除H以外的原子,优化后转fchk再算RESP呢?之后算激发态时,再把优化好的QM层单拿出来算个单点再算RESP,是否是这样的过程?
out文件太大无法上传,gjf上传,另附部分截图。十分感谢。

G09.png (9.44 KB, 下载次数 Times of downloads: 70)

G09 总能量

G09 总能量

G16.png (4.07 KB, 下载次数 Times of downloads: 50)

G16 总能量

G16 总能量

屏幕截图 2020-10-15 203619.png (267.28 KB, 下载次数 Times of downloads: 49)

计算的体系

计算的体系

G16 Br电荷.png (4.41 KB, 下载次数 Times of downloads: 53)

G16 Br的马利肯电荷

G16 Br的马利肯电荷

G09 Br电荷.png (3.77 KB, 下载次数 Times of downloads: 35)

G09 Br的电荷

G09 Br的电荷

g16.gjf

81.52 KB, 下载次数 Times of downloads: 22

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

2#
发表于 Post on 2020-10-16 02:10:57 | 只看该作者 Only view this author
如果就是想计算晶体环境下的某个分子的激发能、光谱,根本没必要用ONIOM,纯粹自找麻烦
直接抠出个团簇,冻结重原子,优化氢原子位置,然后把周围分子都用拟合静电势电荷当成背景电荷描述就行了。
QEq这个对静电势重现很垃圾的电荷根本没法恰当描述晶体环境对当前分子的静电作用。
而且用ONIOM之后fch文件里没有波函数信息,之后都没法用Multiwfn做NTO、空穴-电子分析之类的研究
北京科音自然科学研究中心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

53

帖子

0

威望

739

eV
积分
792

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2020-10-16 11:52:14 | 只看该作者 Only view this author
sobereva 发表于 2020-10-16 02:10
如果就是想计算晶体环境下的某个分子的激发能、光谱,根本没必要用ONIOM,纯粹自找麻烦
直接抠出个团簇, ...

明白了,谢谢sob老师

53

帖子

0

威望

739

eV
积分
792

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2020-10-16 15:40:13 | 只看该作者 Only view this author
sobereva 发表于 2020-10-16 02:10
如果就是想计算晶体环境下的某个分子的激发能、光谱,根本没必要用ONIOM,纯粹自找麻烦
直接抠出个团簇, ...

Sob老师,我详细重复一下这个工作您看是否正确
用VMD抠出团簇,整个团簇冻结重原子,优化氢原子,将这个过程产生的chk转fchk算拟合静电势电荷,然后将周围分子用背景电荷描述,用TDDFT算激发态,之后可以做NTO等分析。

53

帖子

0

威望

739

eV
积分
792

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2020-10-16 16:15:27 | 只看该作者 Only view this author
sobereva 发表于 2020-10-16 02:10
如果就是想计算晶体环境下的某个分子的激发能、光谱,根本没必要用ONIOM,纯粹自找麻烦
直接抠出个团簇, ...

还有一个问题就是如果要计算的团簇很大,是不是可以将中心分子和中心分子周围小范围的分子冻结重原子,最外围全部冻结来优化,然后取出中心分子来算拟合静电势电荷,之后将除了中心分子外全部用背景电荷来描述,正常算TDDFT呢?

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

6#
发表于 Post on 2020-10-17 14:25:10 | 只看该作者 Only view this author
王采岑 发表于 2020-10-16 15:40
Sob老师,我详细重复一下这个工作您看是否正确
用VMD抠出团簇,整个团簇冻结重原子,优化氢原子,将这个 ...

拟合静电势电荷要对单个分子计算,而不是对团簇计算
北京科音自然科学研究中心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

53

帖子

0

威望

739

eV
积分
792

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2020-10-17 15:48:11 | 只看该作者 Only view this author
sobereva 发表于 2020-10-17 14:25
拟合静电势电荷要对单个分子计算,而不是对团簇计算

再请问老师,分子优化的时候,是单个分子冻结重原子优化还是放在团簇里冻结外层分子优化?

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

8#
发表于 Post on 2020-10-18 06:47:03 | 只看该作者 Only view this author
王采岑 发表于 2020-10-17 15:48
再请问老师,分子优化的时候,是单个分子冻结重原子优化还是放在团簇里冻结外层分子优化?

前者
北京科音自然科学研究中心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

53

帖子

0

威望

739

eV
积分
792

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2020-10-18 10:13:55 | 只看该作者 Only view this author

好的,十分感谢您

2

帖子

0

威望

438

eV
积分
440

Level 3 能力者

10#
发表于 Post on 2020-10-19 10:01:03 | 只看该作者 Only view this author

请问社长,要得到激发态的构型是不是可以在优化团簇的氢之后,把周围的分子当作背景电荷,然后优化中心分子的激发态构型

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

11#
发表于 Post on 2020-10-21 12:12:13 | 只看该作者 Only view this author
dsss 发表于 2020-10-19 10:01
请问社长,要得到激发态的构型是不是可以在优化团簇的氢之后,把周围的分子当作背景电荷,然后优化中心分 ...

往往可以接受,但不算严格,毕竟这假定了周围原子的位置都是固定的。

如果用ONIOM,把周围几层分子都用分子力场表示,并且把最外层MM分子坐标冻住,这样会好很多
北京科音自然科学研究中心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

7

帖子

0

威望

175

eV
积分
182

Level 3 能力者

12#
发表于 Post on 2023-2-15 11:02:45 | 只看该作者 Only view this author
sobereva 发表于 2020-10-21 12:12
往往可以接受,但不算严格,毕竟这假定了周围原子的位置都是固定的。

如果用ONIOM,把周围几层分子都 ...

请问老师,那这样用oniom分层做激发之后,在multiwfn中就没法进行NTO和e-h等相关激发态分析了吗?可以直接在输入文件中加入pop(saveNTO,NTO)关键词之后,在高斯中进行NTO分析吗?
谢谢老师!

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

13#
发表于 Post on 2023-2-16 05:14:00 | 只看该作者 Only view this author
guodx 发表于 2023-2-15 11:02
请问老师,那这样用oniom分层做激发之后,在multiwfn中就没法进行NTO和e-h等相关激发态分析了吗?可以直 ...



上面的方式产生fch。至于基于这样的fch和Gaussian输出文件用Multiwfn的空穴-电子分析是否兼容,没测试过,自行测试
北京科音自然科学研究中心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

7

帖子

0

威望

175

eV
积分
182

Level 3 能力者

14#
发表于 Post on 2023-2-16 10:14:53 | 只看该作者 Only view this author
sobereva 发表于 2023-2-16 05:14
上面的方式产生fch。至于基于这样的fch和Gaussian输出文件用Multiwfn的空穴-电子分析是否兼容,没测 ...

好的,谢谢sob老师,我试试

24

帖子

0

威望

169

eV
积分
193

Level 3 能力者

15#
发表于 Post on 2025-3-14 15:13:15 | 只看该作者 Only view this author
guodx 发表于 2023-2-16 10:14
好的,谢谢sob老师,我试试

兄弟 兼容吗

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 19:47 , Processed in 0.285335 second(s), 23 queries , Gzip On.

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