计算化学公社

标题: NICS(1)计算问题 [打印本页]

作者
Author: 落花飘雨679 时间: 2017-6-17 11:15
标题: NICS(1)计算问题
请问各位老师，当计算体系的NICS（1）时，重复了文献中的操作，得到的值和文献不同，我比较困惑的是：NICS(1)默认的是在z方向加1后，看结果中的zz值，还是NICS(1)是指整个Bq的xyz全部加1，而NICS(1)-zz是指在z方向加1？

作者
Author: 冰释之川 时间: 2017-6-17 11:52
本帖最后由冰释之川于 2017-6-17 11:55 编辑

（1）关于NICS(1)_ZZ：NICS(1)_ZZ就是指的在环平面上方1埃处（具体来说，是以环中心为起点向垂直于环平面方向挪1埃），垂直于环平面方向的屏蔽张量分量值的负值。假设分子平面是XY、YZ、XZ平面，那么这个分量值就是分别指的ZZ=、XX=、YY=后面的值 (这个假设很重要！)。对于被研究的环不处在笛卡尔平面上，或者环是扭曲的情况，计算NICS(1)_ZZ是很麻烦的，不过Multiwfn有专门的功能来解决这个问题，见《利用Multiwfn计算倾斜、扭曲环的NICS_ZZ》（http://sobereva.com/261）。

（2）关于NICS(1)和NICS(-1)：在RSC Adv.,6,23900(2016)当中作者对于非平面体系分别定义了NICS(1)和NICS(-1)两个量，分别是凸面和凹面上方1埃处的NICS，还定义了几个相关的量，比如NICS(1)av=[NICS(-1)+NICS(1)]/2用来衡量两个面NICS的平均值、NICS(1)diff=NICS(1)-NICS(-1)用来衡量两个面的NICS差值。文中指出体系的非平面性，以及手性，可能会造成两个面NICS值的分裂，即NICS(1)diff不为零。文中还强调，NICS(1)和NICS(-1)两个量衡量的是两个面各自的芳香性，而不是环的芳香性，用这两个值可以讨论两个面特征的差异性，比如反应性的不同。

作者
Author: sobereva 时间: 2017-6-17 16:33
NICS(1)是平面上方1埃处的各项平均磁屏蔽值的负值，显然不能看ZZ，你算的又不是NICS(1)_ZZ。

仔细看此文对NICS的说明
衡量芳香性的方法以及在Multiwfn中的计算
http://sobereva.com/176

作者
Author: 落花飘雨679 时间: 2017-6-19 09:26

sobereva 发表于 2017-6-17 16:33
NICS(1)是平面上方1埃处的各项平均磁屏蔽值的负值，显然不能看ZZ，你算的又不是NICS(1)_ZZ。

仔细看此文 ...

读了老师的文章受益匪浅，我理解的是NICS(1)为平面上方1埃处，也就是在z方向加1 ，但是看得是Isotropic的值，而不是zz的值，不知道这样的理解是不是正确的？

作者
Author: sobereva 时间: 2017-6-19 16:21

落花飘雨679 发表于 2017-6-19 09:26
读了老师的文章受益匪浅，我理解的是NICS(1)为平面上方1埃处，也就是在z方向加1 ，但是看得是Isotropic的 ...

如果你的体系平面正在XY平面上，是这样

作者
Author: 落花飘雨679 时间: 2017-6-19 16:40

sobereva 发表于 2017-6-19 16:21
如果你的体系平面正在XY平面上，是这样

还要请问老师，我如何确定分子在哪个平面上？我计算的分子本身是平面的，是默认在xy吗？

作者
Author: sobereva 时间: 2017-6-19 17:51

落花飘雨679 发表于 2017-6-19 16:40
还要请问老师，我如何确定分子在哪个平面上？我计算的分子本身是平面的，是默认在xy吗？

在gview里显示笛卡尔坐标轴就知道了

作者
Author: 落花飘雨679 时间: 2017-6-19 20:18

sobereva 发表于 2017-6-19 17:51
在gview里显示笛卡尔坐标轴就知道了

真是太感谢老师了，终于明白了

作者
Author: lycheeho 时间: 2025-9-1 16:08

sobereva 发表于 2017-6-19 17:51
在gview里显示笛卡尔坐标轴就知道了

卢老师您好！我有个分子基态是开壳层单重态，请问算NICS和AICD时直接nmr=csgt uB3LYP/6-31G* iop(10/93=1)这样算呢，还是nmr=csgt uB3LYP/6-31G* iop(10/93=1) guess=mix亦或者先uB3LYP/6-31G* guess=mix stable=opt获得稳定开壳层单重态波函数，再nmr=csgt uB3LYP/6-31G* iop(10/93=1) guess=read读取开壳层波函数计算呢？谢谢！

作者
Author: sobereva 时间: 2025-9-2 09:23

lycheeho 发表于 2025-9-1 16:08
卢老师您好！我有个分子基态是开壳层单重态，请问算NICS和AICD时直接nmr=csgt uB3LYP/6-31G* iop(10/93=1 ...

第一种没意义，得到的还是闭壳层状态的结果
第二种大多数情况可以，但不如第三种稳妥

作者
Author: lycheeho 时间: 2025-9-2 10:20

sobereva 发表于 2025-9-2 09:23
第一种没意义，得到的还是闭壳层状态的结果
第二种大多数情况可以，但不如第三种稳妥

明白了，非常感谢！

作者
Author: lycheeho 时间: 2025-9-2 13:18

sobereva 发表于 2025-9-2 09:23
第一种没意义，得到的还是闭壳层状态的结果
第二种大多数情况可以，但不如第三种稳妥

卢老师，还想请教您一个问题：我想单独考虑π轨道对AICD的贡献，用Multiwfn很方便地输出对应占据共轭π轨道序号后，填入输入文件，为何输出文件里面的轨道数与输入对应不上？就是最后一行轨道的计算没有写入对应test.log文件。

作者
Author: sobereva 时间: 2025-9-3 10:32

lycheeho 发表于 2025-9-2 13:18
卢老师，还想请教您一个问题：我想单独考虑π轨道对AICD的贡献，用Multiwfn很方便地输出对应占据共轭π轨 ...

弄清楚那到底是不是占据轨道

作者
Author: lycheeho 时间: 2025-9-3 15:07

sobereva 发表于 2025-9-3 10:32
弄清楚那到底是不是占据轨道

你好卢老师，解决了，是占据轨道，需要在.gjf文件最好多打几个空行，谢谢~

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