计算化学公社

标题: 求助：能否让高斯输出偶极算符的积分值？ [打印本页]

作者
Author: linzz 时间: 2024-4-7 11:00
标题: 求助：能否让高斯输出偶极算符的积分值？
如图，想得到这个积分的值，但是使用iop 3/33似乎只能得到重叠积分，核对电子吸引的Hamilton以及动能积分。有没有大神知道如何让高斯打印这个积分的值？或者使用pyscf计算这个积分，学生感觉int1e_r这个积分有点像，但是不确定是不是，希望各位老师指教。 (, 下载次数 Times of downloads: 5) (, 下载次数 Times of downloads: 4)

作者
Author: zjxitcc 时间: 2024-4-7 14:17
本帖最后由 zjxitcc 于 2024-4-7 14:18 编辑

iop(3/33=1)有在输出文件里打印偶极矩积分矩阵，搜关键字Multipole matrices

你贴的图片里公式不是很严谨，没有正确表明偶极矩积分的维度。偶极矩算符是矢量算符，有三个分量，例如x分量的矩阵元是R_{ij,x}=<psi_i|R_x|psi_j>，对于y, z分量依此类推。R_{ij,x}是一个三维数组，有两个维度的大小均是分子轨道数目，第三个维度大小是3（即x,y,z三个分量）。
左右矢两个分子轨道可以相同，也可以不同，因此R_{ij,x}对于i,j是一个矩阵，而R_{ii,x}仅仅是这个矩阵的对角元。

作者
Author: linzz 时间: 2024-4-7 16:13

zjxitcc 发表于 2024-4-7 14:17
iop(3/33=1)有在输出文件里打印偶极矩积分矩阵，搜关键字Multipole matrices

你贴的图片里公式不是很严 ...

非常感谢

作者
Author: hebrewsnabla 时间: 2024-4-7 16:13
要用pyscf的话，参考 https://github.com/pyscf/pyscf/blob/master/pyscf/scf/hf.py 的1310行左右。

作者
Author: sobereva 时间: 2024-4-7 19:19
Multiwfn载入fch文件，主功能6里面选7 Print various kinds of integral matrix between basis functions，再选2 Electric dipole moment integral，这是最方便的做法，所有基函数之间的偶极算符的各个分量的积分全都立马给出了。

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