计算化学公社

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

[Multiwfn使用咨询] NICS一维扫描分析pai轨道贡献的扫描曲线不平滑

[复制链接 Copy URL]

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

请问一下各位,我正在算一个三元环体系的2pai电子芳香性,用nmr=csgt计算pai轨道的贡献时,扫描的图得到了下面这样一张图,图中的曲线有各种凸起和凹陷,请问这该怎么理解?扫描点数为400(默认为200),范围是从环中心上下10A。

202405240950101363..png (155.41 KB, 下载次数 Times of downloads: 29)

202405240950101363..png

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

2#
发表于 Post on 2024-5-24 11:52:03 | 只看该作者 Only view this author
那叫pi不叫pai

把gview里显示的包含Bq原子位置的结构图贴出来。可能扫描路径都不合理
并且确保SCF收敛情况,保证每个任务SCF都正常收敛了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-5-24 23:46:09 | 只看该作者 Only view this author
非常感谢sob老师的回复与指正。我放上了包含Bq的结构的侧视图与俯视图,还有计算的输入文件以及输出文件。看起来SCF收敛都没有任何问题。

202405242341063292..png (40.39 KB, 下载次数 Times of downloads: 27)

俯视图

俯视图

202405242341288413..png (99.31 KB, 下载次数 Times of downloads: 26)

侧视图

侧视图

PIn2-pai.com

19.62 KB, 下载次数 Times of downloads: 1

PIn2-pai.log

284.06 KB, 下载次数 Times of downloads: 2

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2024-5-31 03:01:28 | 只看该作者 Only view this author
sobereva 发表于 2024-5-24 11:52
那叫pi不叫pai

把gview里显示的包含Bq原子位置的结构图贴出来。可能扫描路径都不合理

谢谢sob老师的回复与指正,我在帖子里贴了图和原始文件,想麻烦Sob老师答疑解惑

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

5#
发表于 Post on 2024-5-31 10:12:42 | 只看该作者 Only view this author
去掉CSGT,先用默认的GIAO看看结果什么样。没准是CSGT数值求解精度问题
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2024-6-1 00:15:12 | 只看该作者 Only view this author
本帖最后由 hcy89115886 于 2024-6-1 00:31 编辑
sobereva 发表于 2024-5-31 10:12
去掉CSGT,先用默认的GIAO看看结果什么样。没准是CSGT数值求解精度问题

谢谢Sob老师的回复,我有nmr=giao的扫描结果,是正常的,我下面附了图,就是当我用csgt拆分sigma与pi作用时出现了曲线不平滑的问题。

202406010014563466..png (58.81 KB, 下载次数 Times of downloads: 24)

202406010014563466..png

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

7#
发表于 Post on 2024-6-1 10:26:47 | 只看该作者 Only view this author
hcy89115886 发表于 2024-6-1 00:15
谢谢Sob老师的回复,我有nmr=giao的扫描结果,是正常的,我下面附了图,就是当我用csgt拆分sigma与pi作用 ...

你看看如果不区分pi和sigma,用NMR=CSGT得到的总NICS曲线什么样
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2024-6-3 23:03:58 | 只看该作者 Only view this author
sobereva 发表于 2024-6-1 10:26
你看看如果不区分pi和sigma,用NMR=CSGT得到的总NICS曲线什么样

谢谢sob老师的回复。我不区分sigma和Pi的时候,整体比较平滑,但是也可以看到一些锯齿型的,和pi贡献的锯齿型是一个位置。如果只看sigma贡献的话曲线是平滑的。
下图为sigma+pi
下图为sigma



5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

9#
发表于 Post on 2024-6-3 23:42:43 | 只看该作者 Only view this author
尝试不同基组,以及尝试带上CPHF=grid=fine看看是否有改进
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2024-6-5 22:49:31 | 只看该作者 Only view this author
sobereva 发表于 2024-6-3 23:42
尝试不同基组,以及尝试带上CPHF=grid=fine看看是否有改进

谢谢sob老师的回复,我试了def2qzvp,因为好像In适配的基组不是特别多,然后方法也换过b3lyp和PBE0,都是一样的结果,没有任何变化。请问Sob老师,还有其他方法能够在计算芳香性的时候区分Pi轨道和sigma轨道的贡献吗?

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

11#
发表于 Post on 2024-6-6 00:46:07 | 只看该作者 Only view this author
hcy89115886 发表于 2024-6-5 22:49
谢谢sob老师的回复,我试了def2qzvp,因为好像In适配的基组不是特别多,然后方法也换过b3lyp和PBE0,都是 ...

不要无视CPHF=grid=fine


用NBO的NCS功能

量子化学波函数分析与Multiwfn程序培训班(http://www.keinsci.com/workshop/WFN_content.html)专门有一部分讲了,诸如可以输出这样的各个轨道的贡献形式:


北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2024-6-6 23:13:27 | 只看该作者 Only view this author
sobereva 发表于 2024-6-6 00:46
不要无视CPHF=grid=fine

谢谢sob老师的回复,我在换基组的时候就加了CPHF=grid=fine,得到的结果是一样的。我尝试一下这个NBO的NCS功能。

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 09:44 , Processed in 0.246997 second(s), 23 queries , Gzip On.

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