求助：能否让高斯输出偶极算符的积分值？

linzz

如图，想得到这个积分的值，但是使用iop 3/33似乎只能得到重叠积分，核对电子吸引的Hamilton以及动能积分。有没有大神知道如何让高斯打印这个积分的值？或者使用pyscf计算这个积分，学生感觉int1e_r这个积分有点像，但是不确定是不是，希望各位老师指教。

zjxitcc

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}仅仅是这个矩阵的对角元。

linzz

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

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

非常感谢

hebrewsnabla

要用pyscf的话，参考 https://github.com/pyscf/pyscf/blob/master/pyscf/scf/hf.py 的1310行左右。

sobereva

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