请问Multiwfn是否可以计算基态波函数对位矢的体积分？

rugals

需要做一个如下图公式的计算：



(公式出处：https://www.sciencedirect.com/to ... localization-length)

其中ψ_gs是基态波函数，R是位矢，积分是对三维全域的体积分。

sobereva

这个公式相当于假定当前体系只是单电子或双电子体系（后者用一个双占据轨道描述）
这种情况可以实现，需要自己改写Multiwfn源代码定义一个用户自定义函数（波函数与坐标矢量乘积），然后用主功能100的子功能4对此函数进行全空间积分
用户自定义函数的说明看Multiwfn手册2.7节

rugals

sobereva 发表于 2023-12-31 10:38
这个公式相当于假定当前体系只是单电子或双电子体系（后者用一个双占据轨道描述）
这种情况可以实现，需要 ...

谢谢sob老师指点。
原公式出处在1L已经补充。
既然是“ground state eigenfunction”，您也说了单电子体系。我理解这里的波函数Φ_gs应当就是单电子近似下，能量本征值对应的波函数。按您Multiwfn培训班中修改代码的PPT应当对应的是MOene，请您费心看下我的理解是否正确？

sobereva

rugals 发表于 2023-12-31 18:04
谢谢sob老师指点。
原公式出处在1L已经补充。
既然是“ground state eigenfunction”，您也说了单电子 ...

单电子体系和单电子近似不是一码事
单电子体系的波函数变量本身就只有一个r（三维空间矢量）
而N电子体系，不管用不用单电子近似求解，多电子波函数都变量都是r1, r2, r3... rN

MOene是轨道能量数组，fmo(x,y,z,iorbsel)才是iorbsel轨道的波函数。参考function.f90源代码这部分

real*8 function calcfuncall(ifunc,x,y,z)
integer ifunc
real*8 x,y,z
if (ifunc==1) then
        if (allocated(b)) then
                calcfuncall=fdens(x,y,z)
        else
                calcfuncall=calcprodens(x,y,z,0) !Use promolecular density
        end if
else if (ifunc==2) then
        calcfuncall=fgrad(x,y,z,'t')
else if (ifunc==3) then
        calcfuncall=flapl(x,y,z,'t')
else if (ifunc==4) then
        calcfuncall=fmo(x,y,z,iorbsel)
else if (ifunc==5) then
        calcfuncall=fspindens(x,y,z,'s')

...略