计算化学公社

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

[NBO相关] 请教在计算NBO和WBI键级遇到输出文件中无法找到结果的问题。

[复制链接 Copy URL]

9

帖子

0

威望

1149

eV
积分
1158

Level 4 (黑子)

      请问各位老师,我在得到物质单晶结构后对物质进行了几何优化,使用的版本是Gaussian 03,优化时用的关键词是# opt freq b3lyp/gen pseudo=read,所得到的结果在GaussView上看结构有一些键没连上或者键型不对。我在得到优化结果后用优化的结果进行了NBO计算,计算的关键词是# b3lyp/gen pop=(full,nboread,savenbos) pseudo=read,得到的结果有一些键的NBO数值没有显示,请问这个结果可信吗?在之后想得到WBI键级,在输入文件关键词# b3lyp/gen pop=nbo pseudo=read,并在输入文件最后加入了$nbo bndidx $end,在输出文件中没有找到Wiberg bond index matrix in the NAO basis: ,请问是在计算过程中有什么错误我没有找到结果吗?所有的计算都得到了Normal termination of Gaussian 03。
谢谢各位老师解答与指正。

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2018-7-4 02:19:27 | 只看该作者 Only view this author
仔细看
谈谈原子间是否成键的判断问题
http://sobereva.com/414

你不写nboread就代表不会从末尾读取$NBO信息,显然不会输出Wiberg键级矩阵

建议用Multiwfn计算Mayer键级。如果是有机体系,同时推荐用Multiwfn计算我提出的拉普拉斯键级,看J. Phys. Chem. A, 117, 3100−3108 (2013)
北京科音自然科学研究中心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

9

帖子

0

威望

1149

eV
积分
1158

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2018-7-4 17:39:28 | 只看该作者 Only view this author
sobereva 发表于 2018-7-4 02:19
仔细看
谈谈原子间是否成键的判断问题
http://sobereva.com/414

Sob老师,我用fchk文件里的结构进行了计算,使用的# b3lyp/gen pop=(full,nboread,savenbo) pseudo=read,计算结尾加入了$nbo bndidx $end,在输出文件中出现了 NBStor is confused about NOcc. Error termination via Lnk1e in /home/luyq/g03/l607.exe我进行了几次更改是加入文件末尾的$nbo bndidx $end造成的,请问如过需要计算WBI,应该怎么调整,还有如果计算拉普拉斯键级在高斯计算中需要计算那些值?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

4#
发表于 Post on 2018-7-4 23:49:17 | 只看该作者 Only view this author
q_7997 发表于 2018-7-4 17:39
Sob老师,我用fchk文件里的结构进行了计算,使用的# b3lyp/gen pop=(full,nboread,savenbo) pseudo=read ...


Gaussian自带的NBO 3.1 bug很多,当前不成功就甭用NBO了

直接用fch就能用Multiwfn计算各种键级
详见
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379
Multiwfn入门tips
http://sobereva.com/167
北京科音自然科学研究中心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, 2025-8-15 16:59 , Processed in 0.901729 second(s), 20 queries , Gzip On.

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