计算化学公社

标题: 用multiwfn算的福井函数与MS dmol3的结果明显不同 [打印本页]
作者
Author: vickywu    时间: 2020-4-8 23:29
标题: 用multiwfn算的福井函数与MS dmol3的结果明显不同
请各位大神指点：

最近在做quantitative structure-reactivity relationship，然后看的文献上有机物的福井函数都是用dmol3算的，但是我还是想用multiwfn，所以做了一些功课，但是算出来的结果和文献上MS算的结果差挺大的，因此请教各位，指点一二，看看我的计算步骤或者概念上是不是有什么问题？在此先感谢各位！

首先，先放上文献上所述的MS计算参数设定：
[url=][/url]
[url=][/url]

我用Multiwfn算的时候是完全按照 http://sobereva.com/484 进行的计算，得到CDFT.txt文档。该结果如下所示：其中N, N+1,N-1的gjf输入文件是multiwfn自动生成的，所以绝对没有问题。
Atomq(N)q(N+1)q(N-1)f-f+f0CDD
19(H
0.1263
0.1115
0.1474
0.0211
0.0148
0.0179
-0.0063
14(S
0.5085
0.4903
0.5331
0.0247
0.0182
0.0214
-0.0065
18(H
0.1282
0.1056
0.1557
0.0275
0.0226
0.025
-0.0048
10(H
0.0506
-0.0046
0.0813
0.0307
0.0552
0.0429
0.0246
7(H
0.0542
0.0042
0.0852
0.031
0.05
0.0405
0.0191
8(H
0.0436
-0.0145
0.0803
0.0367
0.0581
0.0474
0.0214
9(H
0.043
-0.0119
0.0811
0.0381
0.0548
0.0465
0.0168
17(N
-0.2079
-0.2214
-0.1647
0.0432
0.0135
0.0283
-0.0296
5(C
-0.0308
-0.1577
0.0136
0.0444
0.127
0.0857
0.0826
1(C
-0.0259
-0.13
0.0189
0.0448
0.1041
0.0745
0.0593
12(H
0.1275
0.0982
0.1783
0.0508
0.0293
0.04
-0.0214
13(H
0.1275
0.0963
0.1786
0.0511
0.0312
0.0411
-0.0199
15(O
-0.3333
-0.3694
-0.2784
0.0549
0.0361
0.0455
-0.0187
16(O
-0.333
-0.3675
-0.2763
0.0567
0.0345
0.0456
-0.0222
3(C
0.0579
0.0133
0.119
0.061
0.0446
0.0528
-0.0164
2(C
-0.0602
-0.1826
0.0084
0.0685
0.1224
0.0955
0.0539
4(C
-0.0632
-0.1733
0.0128
0.076
0.1101
0.0931
0.0341
6(C
-0.0487
-0.0845
0.041
0.0897
0.0358
0.0628
-0.0539
11(N
-0.1647
-0.2022
-0.0154
0.1493
0.0375
0.0934
-0.1118

第一个问题是：表中的f+,f-,f0是否就是各原子的简缩福井函数？
第二个问题是：multiwfn中计算N,N-1,N+1的wfn文件是选择的默认的B3LYP/6-31G*，但是文献里的MS计算似乎并没有说明基组，是否是由于基组不同导致的数据差异？

下表是文献值和我计算出来的值，进行了一个对比：
Organic compounds
BOx
BOn
f(+)x
f(+)n
f(-)x
f(-)n
f(0)x
f(0)n
Sulfanilamide/文献值
1.963
0.927
0.091
0.028
0.123
0.025
0.096
0.026
测试值
1.668

0.127
0.014
0.149
0.021
0.096
0.018


可见，尤其f(+)的最大值、最小值；f(-)的最大值以及f(0)的最小值都有很明显的差异，因此请大神指点一下，我做的福井函数计算是否有问题。

PS.在做键级计算时，我直接用multiwfn的mayer键级分析的，可以看到和文献值还是对不上，而文献值也是从dmol3算的，具体哪种键级分析文献没有提及，请问我直接算mayer键级是否就可以？


作者
Author: vickywu    时间: 2020-4-8 23:33
sob老师抱歉，我没有注意有【波函数分析求助】的分栏，我选错主题分类了，抱歉抱歉！
作者
Author: 卡开发发    时间: 2020-4-9 06:28
一、是。
二、导致差异的因素比较多，这是我觉得有可能出现的情况：
1、计算简缩fukui函数的时候需要N、N±△N的原子电荷，原子电荷的选择并不是唯一的，若是用Mulliken还存在基组依赖的问题，数值基组和GTO结果不见得一样，文献对此也没做什么说明。
2、DMol3计算fukui函数的时候△N的默认选择并不是1。
3、方法上文献也没具体说明，到底是用的何泛函。
4、DMol3能做的键级有Mulliken和Mayer，文献还是没说明清楚。

综合上面的来说：
1、你找的这篇文献给的信息太少，该写到文章的没写，不该写的如那个SCF设置到500圈写的很6啊，让人一看就知道是外行做的233333.
2、DMol3算小分子体系没啥优势，老老实实按http://sobereva.com/484来做不会有什么问题。
作者
Author: sobereva    时间: 2020-4-9 06:53
用Dmol3算有机物的福井函数在我来看才是异类，见
为什么我非常不建议购买和使用D摩3量子化学程序
http://sobereva.com/508

1 一个原子一个数值，显然是简缩福井函数。
2 计算时候必须说明计算级别，如果这么重要信息都不说的话，作者水平堪忧
不同计算级别下的波函数分析结果没有可比性。
纯泛函、杂化泛函差异是明显的。Dmol3的用户往往都在用过时的纯泛函（因为杂化泛函在Dmol3里巨慢），结果自然和B3LYP会相差很大。更何况Dmol3的基函数类型和主流量子化学程序用的Gaussian型基函数相差也很大。

以后分错类直接编辑原帖进行修改就行了，不用在回帖里说明
作者
Author: vickywu    时间: 2020-4-9 12:04
卡开发发 发表于 2020-4-9 06:28
一、是。
二、导致差异的因素比较多，这是我觉得有可能出现的情况：
1、计算简缩fukui函数的时候需要N、N ...

谢谢老师，明白了！受教受教
作者
Author: vickywu    时间: 2020-4-9 12:05
sobereva 发表于 2020-4-9 06:53
用Dmol3算有机物的福井函数在我来看才是异类，见
为什么我非常不建议购买和使用D摩3量子化学程序
http:// ...

谢谢老师，那我就安心按您的帖子做了。十分感谢！



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3