计算化学公社

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

[波函数分析求助] NICS(1)zz计算结果与文献相差很大

[复制链接 Copy URL]

52

帖子

0

威望

1477

eV
积分
1529

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
最近在学习计算NICS(1)zz,就参考了文献On splitting of the NICS(1) magnetic aromaticity index (RSC Adv., 2016, 6, 23900)中corannulene的例子进行练习,但是发现我计算出来的结果与文献里面报道的结果相差很大,请问是我算错了吗?
我先用opt freq b3lyp/6-31g**优化了corannulene的结构,然后用multiwfn给优化后的分子添加gost atom,然后再用nmr=giao b3lyp/6-31g**算它的NMR,最后再用multiwfn计算它的NICS(1)zz.
Ring NICS(-1)zz NICS(0)zz NICS(1)zz
A5(文献结果) -2.37 9.39 2.56
A5(我的计算) 15.22 56.22 11.52
A(文献结果) -12.77-7.17 -6.49
A(我的计算) -23.22 -1.04 -15.16

谢谢大家!

2022-05-19_165724.png (24.95 KB, 下载次数 Times of downloads: 19)

2022-05-19_165724.png

6万

帖子

99

威望

6万

eV
积分
125147

管理员

公社社长

2#
发表于 Post on 2022-5-19 17:06:09 | 只看该作者 Only view this author
1 确保计算级别和文献一致
2 检查gjf文件里的ghost原子位置,看到底对不对
3 检查NICS读取是否正确,说不清楚就把相关文件和相关操作贴出来

以上都确保无误的情况下,和文献对不上大概率是文献里的数据有问题
北京科音自然科学研究中心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

52

帖子

0

威望

1477

eV
积分
1529

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2022-5-19 17:35:05 | 只看该作者 Only view this author
本帖最后由 zhangzb 于 2022-5-19 18:21 编辑
sobereva 发表于 2022-5-19 17:06
1 确保计算级别和文献一致
2 检查gjf文件里的ghost原子位置,看到底对不对
3 检查NICS读取是否正确,说不 ...

谢谢社长!
1. 计算级别与文献一致,都是b3lyp/6-31g**;
2. 检查了gif文件里ghost原子的位置,没有发现错误。
3. 对不起,是我搞错了,我用multiwfn计算出来的是NICS(1)zz, 文献中给的是NICS(1).

cora_nics.gjf.Output.1078531.hpc-mn1.log

59.1 KB, 下载次数 Times of downloads: 4

52

帖子

0

威望

1477

eV
积分
1529

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2022-5-19 17:42:03 | 只看该作者 Only view this author
本帖最后由 zhangzb 于 2022-5-19 18:17 编辑
sobereva 发表于 2022-5-19 17:06
1 确保计算级别和文献一致
2 检查gjf文件里的ghost原子位置,看到底对不对
3 检查NICS读取是否正确,说不 ...

我好像知道问题出在哪里了。文献里的数值取的是Isotropic的负值,而我是根据http://sobereva.com/261这里面的方法计算的。谢谢!

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 14:09 , Processed in 0.531166 second(s), 23 queries , Gzip On.

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