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

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

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

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

非常感谢

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