计算化学公社

标题: NICS(1)zz计算结果与文献相差很大 [打印本页]

作者
Author: zhangzb 时间: 2022-5-19 17:00
标题: NICS(1)zz计算结果与文献相差很大
最近在学习计算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

谢谢大家！

作者
Author: sobereva 时间: 2022-5-19 17:06
1 确保计算级别和文献一致
2 检查gjf文件里的ghost原子位置，看到底对不对
3 检查NICS读取是否正确，说不清楚就把相关文件和相关操作贴出来

以上都确保无误的情况下，和文献对不上大概率是文献里的数据有问题

作者
Author: zhangzb 时间: 2022-5-19 17:35
本帖最后由 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).

作者
Author: zhangzb 时间: 2022-5-19 17:42
本帖最后由 zhangzb 于 2022-5-19 18:17 编辑

sobereva 发表于 2022-5-19 17:06
1 确保计算级别和文献一致
2 检查gjf文件里的ghost原子位置，看到底对不对
3 检查NICS读取是否正确，说不 ...

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

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