fch文件与wfx文件做拓补分析结果不一致，该如何选择？

wuzhao541

大家好，我对某一聚合物体系做AIM拓补分析的时候发现，fch文件和wfx文件在相同的两个原子之间找到的键临界点稍有差异，
例如某一点用fch得到的结果是：
Position (Bohr):  -11.46304627561961    0.41468070588894   -0.00382641801727
Density of all electrons:  0.1864601015E+00
Density of Alpha electrons:  0.9323005074E-01
Density of Beta electrons:  0.9323005074E-01
Spin density of electrons:  0.0000000000E+00
Lagrangian kinetic energy G(r):  0.2270374465E+00
Hamiltonian kinetic energy K(r):  0.1904218924E+00
Potential energy density V(r): -0.4174593389E+00
Energy density E(r) or H(r): -0.1904218924E+00
Laplacian of electron density:  0.1464622161E+00
而用wfx得到的结果是：
Position (Bohr):  -11.51042994622739    0.39560042599133   -0.00352607742299
Density of all electrons:  0.1870492451E+00
Density of Alpha electrons:  0.9352462256E-01
Density of Beta electrons:  0.9352462256E-01
Spin density of electrons:  0.0000000000E+00
Lagrangian kinetic energy G(r):  0.1790811160E+00
Hamiltonian kinetic energy K(r):  0.1988960882E+00
Potential energy density V(r): -0.3779772042E+00
Energy density E(r) or H(r): -0.1988960882E+00
Laplacian of electron density: -0.7925988867E-01
两者的电子密度拉普拉斯值差别较大。是否是因为两者的坐标精度不同导致的？wfx文件中原子坐标精确到了小数点后12位，而用来转化成fch文件的chk文件中原子坐标精确到小数点后8位。
那么实际计算该以哪个为准呢？因为有些功能wfx文件用不了，是否该按fch文件为准？

冰释之川

wfx 这种波函数文件一般更加精确一点（比wfn精确）

sobereva

你第一个输出对应的点是
-11.46304627561961    0.41468070588894   -0.00382641801727
第二个是
-11.51042994622739    0.39560042599133   -0.00352607742299
坐标都不一样，结果肯定不一样

fch就已经足够精确了

beefly

如果用了ECP，用fch做aim分析的结果可能不可靠

wuzhao541

sobereva 发表于 2016-5-6 12:04
你第一个输出对应的点是
-11.46304627561961    0.41468070588894   -0.00382641801727
第二个是

谢谢sob老师，这两个键临界点的坐标都是由Multiwfn计算出来的，对应成键的两个原子在wfx文件中坐标为：
-1.033140474471e+01  8.546260006671e-01 -1.079054049250e-02
-1.331889190292e+01 -4.484265514901e-01  9.733535689518e-03
对应成键的两个原子在fch文件中坐标为：
-1.03314047E+01  8.54626001E-01 -1.07905405E-02
-1.33188919E+01 -4.48426551E-01 9.73353569E-03
可能这点结构的细微差异导致找到的键临界点位置有偏差吧。

wuzhao541

beefly 发表于 2016-5-6 14:12
如果用了ECP，用fch做aim分析的结果可能不可靠

谢谢讨论，我的体系用的MP2方法和aug-cc-pVDZ基组，没有用ECP

sobereva

wuzhao541 发表于 2016-5-6 15:56
谢谢sob老师，这两个键临界点的坐标都是由Multiwfn计算出来的，对应成键的两个原子在wfx文件中坐标为：
...

MP2计算时，用了density关键词时，out=wfn或wfx得到的波函数文件里记录的是MP2自然轨道
但是如果你直接用fch文件，里面还是fch轨道，显然结果会和wfn/wfx不一样。要得到与wfn/wfx一致的结果，必须把自然轨道转存到fch文件里，看Multiwfn手册第四章开头的说明。