用multiwfn算的福井函数与MS dmol3的结果明显不同

vickywu

请各位大神指点：

最近在做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自动生成的，所以绝对没有问题。

Atom	q(N)	q(N+1)	q(N-1)	f-	f+	f0	CDD
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键级是否就可以？

vickywu

sob老师抱歉，我没有注意有【波函数分析求助】的分栏，我选错主题分类了，抱歉抱歉！

卡开发发

一、是。
二、导致差异的因素比较多，这是我觉得有可能出现的情况：
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来做不会有什么问题。

sobereva

用Dmol3算有机物的福井函数在我来看才是异类，见
为什么我非常不建议购买和使用D摩3量子化学程序
http://sobereva.com/508

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

以后分错类直接编辑原帖进行修改就行了，不用在回帖里说明

vickywu

卡开发发 发表于 2020-4-9 06:28
一、是。
二、导致差异的因素比较多，这是我觉得有可能出现的情况：
1、计算简缩fukui函数的时候需要N、N ...

谢谢老师，明白了！受教受教

vickywu

sobereva 发表于 2020-4-9 06:53
用Dmol3算有机物的福井函数在我来看才是异类，见
为什么我非常不建议购买和使用D摩3量子化学程序
http:// ...

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