计算化学公社

标题: NICS一维扫描分析pai轨道贡献的扫描曲线不平滑 [打印本页]

作者
Author: hcy89115886 时间: 2024-5-24 09:50
标题: NICS一维扫描分析pai轨道贡献的扫描曲线不平滑
请问一下各位，我正在算一个三元环体系的2pai电子芳香性，用nmr=csgt计算pai轨道的贡献时，扫描的图得到了下面这样一张图，图中的曲线有各种凸起和凹陷，请问这该怎么理解？扫描点数为400（默认为200），范围是从环中心上下10A。

作者
Author: sobereva 时间: 2024-5-24 11:52
那叫pi不叫pai

把gview里显示的包含Bq原子位置的结构图贴出来。可能扫描路径都不合理
并且确保SCF收敛情况，保证每个任务SCF都正常收敛了

作者
Author: hcy89115886 时间: 2024-5-24 23:46
非常感谢sob老师的回复与指正。我放上了包含Bq的结构的侧视图与俯视图，还有计算的输入文件以及输出文件。看起来SCF收敛都没有任何问题。

作者
Author: hcy89115886 时间: 2024-5-31 03:01

sobereva 发表于 2024-5-24 11:52
那叫pi不叫pai

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

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

作者
Author: sobereva 时间: 2024-5-31 10:12
去掉CSGT，先用默认的GIAO看看结果什么样。没准是CSGT数值求解精度问题

作者
Author: hcy89115886 时间: 2024-6-1 00:15
本帖最后由 hcy89115886 于 2024-6-1 00:31 编辑

sobereva 发表于 2024-5-31 10:12
去掉CSGT，先用默认的GIAO看看结果什么样。没准是CSGT数值求解精度问题

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

作者
Author: sobereva 时间: 2024-6-1 10:26

hcy89115886 发表于 2024-6-1 00:15
谢谢Sob老师的回复，我有nmr=giao的扫描结果，是正常的，我下面附了图，就是当我用csgt拆分sigma与pi作用 ...

你看看如果不区分pi和sigma，用NMR=CSGT得到的总NICS曲线什么样

作者
Author: hcy89115886 时间: 2024-6-3 23:03

sobereva 发表于 2024-6-1 10:26
你看看如果不区分pi和sigma，用NMR=CSGT得到的总NICS曲线什么样

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

作者
Author: sobereva 时间: 2024-6-3 23:42
尝试不同基组，以及尝试带上CPHF=grid=fine看看是否有改进

作者
Author: hcy89115886 时间: 2024-6-5 22:49

sobereva 发表于 2024-6-3 23:42
尝试不同基组，以及尝试带上CPHF=grid=fine看看是否有改进

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

作者
Author: sobereva 时间: 2024-6-6 00:46

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）专门有一部分讲了，诸如可以输出这样的各个轨道的贡献形式：
(, 下载次数 Times of downloads: 26)

作者
Author: hcy89115886 时间: 2024-6-6 23:13

sobereva 发表于 2024-6-6 00:46
不要无视CPHF=grid=fine

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

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