如何研究δ，d和多重芳香性

风飞 · 发表于 Post on 2015-12-15 16:50:12

老师：您好！通过看文献了解到芳香性可以分为σ，π，δ，d和多重芳香性。您的帖子里介绍了将核独立化学位移(NICS)分解为sigma和pi轨道的贡献，请问如果要研究δ，d和多重芳香性。您有相关的帖子吗？或者是建议吗？谢谢！

sobereva · 发表于 Post on 2015-12-15 23:37:28

此文很多方法都可以讨论
衡量芳香性的方法以及在Multiwfn中的计算
http://sobereva.com/176

最简单的就是看MO图形，看看形成满足4n+2组合的MO是否是原子的d轨道参与的。MO不好考察时也可以用AdNDP方法考察。
如果你想再定量化一点，比如研究δ芳香性时，用Multiwfn计算多中心键级前或计算ELF前只保留构成δ芳香性的轨道，把其它轨道占据数都清零。
AICD方法也可以只考察某些MO对环电流、NICS的贡献。
多重芳香性就是有多类电子都对芳香性有贡献，比如同时具有sigma芳香性和pi芳香性，如Al4 2-。
这篇文章也可以看看：doi: 10.1002/wcms.1115

风飞 · 发表于 Post on 2015-12-16 20:18:00

老师：您好！我在用高斯算NICS的时候，加入Bq原子以后，高斯报错说“Atomic number out of range in Df2TZV.”我之前的体系是全金属的所以采用的方法和基组是PBE1PBE/def2TZVP,如果是这样，我加入Bq以后我是否需要使用混合基组，加上6-311G之类的，突然想起上次您说计算的NICS的基组可以和结构优化的不一样，如果是这样，对于我的体系只含有（Au/Sb)这两种金属，用什么基组来计算比较合理和划算呢？谢谢！

sobereva · 发表于 Post on 2015-12-17 01:31:28

风飞发表于 2015-12-16 20:18
老师：您好！我在用高斯算NICS的时候，加入Bq原子以后，高斯报错说“Atomic number out of range in Df2TZV ...

我用g09 D.01试了下def2TZVP加上Bq，并没出错，可以把你输入文件传上来看看。
算NICS还用def2TZVP就行。

风飞 · 发表于 Post on 2015-12-17 09:58:44

本帖最后由风飞于 2015-12-17 10:01 编辑

老师：您好！输入文件我是以附录的形式上传的，麻烦您看看。我的Bq原子一共设定了17个，会不会是这个Bq原子的数量超过了呢？谢谢！

sobereva · 发表于 Post on 2015-12-17 11:36:03

风飞发表于 2015-12-17 09:58
老师：您好！输入文件我是以附录的形式上传的，麻烦您看看。我的Bq原子一共设定了17个，会不会是这个Bq原子 ...

你就直接把gjf文件传上来，我好测试

风飞 · 发表于 Post on 2015-12-17 12:07:01

sobereva 发表于 2015-12-17 11:36
你就直接把gjf文件传上来，我好测试

输入文件已上传，谢谢老师

sobereva · 发表于 Post on 2015-12-17 12:34:30

风飞发表于 2015-12-17 12:07
输入文件已上传，谢谢老师

当前确实没法算。你用genecp，去EMSL拷def2-TZVP的定义拿过来用就行了

风飞 · 发表于 Post on 2015-12-29 11:30:16

本帖最后由风飞于 2015-12-31 20:49 编辑

老师：您好！我在用我的体系绘制化学屏蔽表面的时候，在输入文件中用的基组和方法是PBE1PBE/def2TZVP,，然后出现的错误是Atomic number out of range in DF2TZV, 突然想到了您之前说过用def2TZVP 计算全金属芳香性的时候要用混合基组，并且要拷入Def2TZVP 的定义，。我进行修改后，出现的错误是：NSgDBf>NBsAll in quad4. 请问是怎么回事呢？老师修改后的输入文件和对应的出错out文件已上传，谢谢老师！

sobereva · 发表于 Post on 2015-12-31 01:50:45

风飞发表于 2015-12-29 11:30
老师：您好！我在用我的体系绘制化学屏蔽表面的时候，在输入文件中用的基组和方法是PBE1PBE/def2TZVP,，然 ...

D.01容易有这个问题，写上guess=huckel可解决

风飞 · 发表于 Post on 2015-12-31 17:21:43

本帖最后由风飞于 2015-12-31 17:32 编辑

老师：您好！结合您的帖子”将核独立化学位移(NICS)分解为sigma和pi轨道的贡献“+“使用Multiwfn绘制NBO及相关轨道”中的“#P b3lyp/6-31+G* NMR pop=NCSall IOp(6/65=-1)”+“在Gaussian输入文件中加上pop=nboread关键词，在末尾空一行写上比如$NBO plot file=c:\ltwd\NH2COH $END，”的内容我的输入文件如下：”#p pbe1pbe/genecp NMR pop(ncsall,nboread) iop(6/65=-1)；最后一行：+$nbo plot file=F:\quanjinshu\NICS\4huan\SbAuSb(-2)-D4h-fenjie-2 $end“ ，高斯虽然正常结束，但是输出文件中只有一般NICS的输出结果，并没有”各个NLMO对NICS(0)_ZZ的贡献也就是Total那列的负值“并且提示” The NLMO program is not currently set up to handle this.“请问，这是什么原因呢？我的是Sandwich化合物，设置的bq原子比较多是按照（-1，-0.5，0,0.5,1····）这样来说设置的是否是我的Bq原子设置问题，还是？谢谢老师！NICS检验，该物质是有芳香性的，是不是我的体系中没有sigma或者是π芳香性呢？

sobereva · 发表于 Post on 2015-12-31 18:22:00

风飞发表于 2015-12-31 17:21
老师：您好！结合您的帖子”将核独立化学位移(NICS)分解为sigma和pi轨道的贡献“+“使用Multiwfn绘制NBO及 ...

这说明没有进行NLMO分析，看NLMO program is not currently set up to handle this上面的具体提示，多半是由于自带的NBO3.1程序的限制导致的。这和有没有sigma、pi芳香性，以及怎么设Bq没有关系。

风飞 · 发表于 Post on 2015-12-31 20:09:50

老师：我用的是高斯是D.01，版本的请问除了考虑换版本，有其他解决的办法吗？谢谢老师！

风飞 · 发表于 Post on 2015-12-31 22:49:59

sobereva 发表于 2015-12-31 01:50
D.01容易有这个问题，写上guess=huckel可解决

老师：您好！发现加上guess=huckel很费时间，高斯最近出新版本了，老师您觉得有必要换新的算吗？如果有必要，请问咱们群里有新版本的高斯吗？

sobereva · 发表于 Post on 2015-12-31 22:54:42

风飞发表于 2015-12-31 22:49
老师：您好！发现加上guess=huckel很费时间，高斯最近出新版本了，老师您觉得有必要换新的算吗？如果有必 ...

并没觉得费什么额外时间。E.01目前网上没有，应该也不会在这方面有什么改变

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[波函数分析求助] 如何研究δ，d和多重芳香性