计算化学公社

标题: 请问Multiwfn是否可以计算基态波函数对位矢的体积分？ [打印本页]

---

作者
Author: rugals    时间: 2023-12-30 15:00
标题: 请问Multiwfn是否可以计算基态波函数对位矢的体积分？


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

(, 下载次数 Times of downloads: 22)

上传 Uploaded

点击下载
Click to download

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

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

---

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

---

作者
Author: rugals    时间: 2023-12-31 18:04

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

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

---

作者
Author: sobereva    时间: 2024-1-1 03:43

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')

...略

---

欢迎光临 计算化学公社 (http://bbs.keinsci.com/)

Powered by Discuz! X3.3