计算化学公社

标题: 计算占据分子轨道上的NICS(0)_zz数值 [打印本页]

作者
Author: xxx 时间: 2016-6-28 12:33
标题: 计算占据分子轨道上的NICS(0)_zz数值
老师您好，我想用Multiwfn分别计算一下各个轨道上的NICS(0)_zz数值（如下图所示），但是没有看到相关计算的帖子。老师指点一下计算原理和步骤吧。
另外，我用Multiwfn计算了一下文献上的带两个负电荷的Al4团簇的NICS(0)_zz数值，我用的基组和方法以及结构优化的结果和文献上是一样的，但发现NICS(0)_zz数值与文献差别较大（文献上是-30.9,我算出来是-54.8730），计算步骤我确定是正确的，但是不知道哪里出错了。请老师指点

作者
Author: sobereva 时间: 2016-6-28 12:44
1 看此文第一节第二段
将核独立化学位移(NICS)分解为sigma和pi轨道的贡献
http://sobereva.com/145

另外利用AICD貌似也可以实现只考虑某些轨道，因为计算时可以在输入文件里指定CSGT计算时只考虑哪些轨道。

2 本身就是个平面体系，计算NICS(0)_ZZ用不着Multiwfn。你再仔细检验下磁屏蔽张量值分量。如果就是和文献对不上就只能问作者了。

作者
Author: xxx 时间: 2016-6-28 20:37

sobereva 发表于 2016-6-28 12:44
1 看此文第一节第二段
将核独立化学位移(NICS)分解为sigma和pi轨道的贡献
http://sobereva.com/145

谢谢Sob老师的回答。因为不太会所以还想再确认一下，对于平面的结构用高斯计算过后会直接在结果中得到NICS(0)_zz或者NICS(1)_zz（如下图的结果所示）。如果对于非平面的结构，要想确定NICS(0)_zz或者NICS(1)_zz数值还需要用Multiwfn进行计算。请问Sob老师，是这样的吗？

作者
Author: xxx 时间: 2016-6-28 21:15

xxx 发表于 2016-6-28 20:37
谢谢Sob老师的回答。因为不太会所以还想再确认一下，对于平面的结构用高斯计算过后会直接在结果中得到NIC ...

刚刚忘了把图片传上，现在补传一下

作者
Author: sobereva 时间: 2016-6-29 03:19

xxx 发表于 2016-6-28 20:37
谢谢Sob老师的回答。因为不太会所以还想再确认一下，对于平面的结构用高斯计算过后会直接在结果中得到NIC ...

是。因为Z具体是指垂直于分子平面方向。计算平面结构时，高斯自动会把分子放到某个笛卡尔平面上，故平面的法向量和某个笛卡尔轴平行，直接从输出文件中读输相应分量就行。而非平面分子，垂直于分子平面的方向肯定不是正好顺着某个笛卡尔轴的，只能用Multiwfn帮助获得相应分量的值。

作者
Author: xxx 时间: 2016-6-29 08:27

sobereva 发表于 2016-6-29 03:19
是。因为Z具体是指垂直于分子平面方向。计算平面结构时，高斯自动会把分子放到某个笛卡尔平面上，故平 ...

谢谢Sob老师的详细解答

作者
Author: janihyhubo 时间: 2017-3-22 18:35

sobereva 发表于 2016-6-29 03:19
是。因为Z具体是指垂直于分子平面方向。计算平面结构时，高斯自动会把分子放到某个笛卡尔平面上，故平 ...

社长，您的意思是，如果分子平面是yz平面，那么垂直于分子平面的就是X轴方向，所以提取所谓的nics(1)_zz应该是提取xx方向的分量吧？

作者
Author: sobereva 时间: 2017-3-22 19:06

janihyhubo 发表于 2017-3-22 18:35
社长，您的意思是，如果分子平面是yz平面，那么垂直于分子平面的就是X轴方向，所以提取所谓的nics(1)_zz ...

对

作者
Author: janihyhubo 时间: 2017-3-22 19:18

sobereva 发表于 2017-3-22 19:06
对

好的，谢谢社长

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