计算化学公社

标题: 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里显示笛卡尔坐标轴就知道了

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

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