求助查找每一个分子S1态和T1态的交换能和库仑积分

York

老师们好，我想把每一个分子S1态和T1态的交换能和库仑积分找到，请问高斯输出文件能查到吗？

北大-陶豫

如果想把能量分解为交换、相关以及其他贡献，应该使用能量分解，参考FAQ：

http://gaussian.com/faq1/
可以通过读取checkpoint文件里的density，再通过extralinks=l608实现。例如：

%chk=h2co.chk
#p b3lyp/6-31g(d)

Formaldehyde

0 1
C                -0.545946     0.000000    0.000000
O                 0.698692     0.000000    0.000000
H                -1.156931     0.000000    0.938173
H                -1.156931     0.000000   -0.938173

--link1--
%chk=h2co.chk
# Geom=AllCheck ChkBas Guess=(Read,Only) Density=Check extralinks=l608

<数值1> <数值2>

注意<数值1> <数值2>独立为一行，之前有一空行。<数值1>需要查看泛函对应的IOp(3/74)值实现，而<数值2>取决于DFT积分网格。例如，B3LYP 对应的<数值1>是"-5"；“FineGrid” (G09默认) 和"UltraFineGrid"(G16默认)对应的<数值2>分别是4和5。输出格式如下：

Results using SCF density:
FoFJK:  IHMeth= 0 ICntrl=     500 DoSepK=F KAlg= 0 I1Cent=   0 FoldK=F
IRaf=         0 NMat=   1 IRICut=       1 DoRegI=T DoRafI=F ISym2E= 0.
FoFCou: FMM=F IPFlag=           0 FMFlag=           0 FMFlg1=           0
         NFxFlg=           0 DoJE=F BraDBF=F KetDBF=F FulRan=T
         wScrn=  0.000000 ICntrl=     500 IOpCl=  1 I1Cent=           0 NGrid=           0
         NMat0=    1 NMatS0=      1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
Symmetry not used in FoFCou.
ET=  113.323772 EV= -329.178547 EJ=   85.580333 EK=  -14.005407 ENuc=   30.433784

A=0 I=   1 IExCor=  402 IRadAn=       4 ScaHFX=  0.200000
            ScaDFX=  0.800000  0.720000  1.000000  0.810000
HFx=  0.000000  0.000000  1.000000  0.000000  0.000000
DFx=  0.000000  0.000000  1.000000  0.000000  0.000000
DFc=  0.000000  0.000000  1.000000  0.000000  0.000000
            ENTVJ=  -99.840658 Ex=  -13.961347 Ec=   -0.695994 ETotM2e= -197.2772497241  ETot= -114.4979980220


其中 ETot = ET + EV + EJ + ENuc + Ex + Ec

ET：electron kinetic energy
EV：nuclear attraction
EJ： Coulomb exchange
ENuc: nuclear repulsion
Ex：<数值> 对应方法给出的 exchange energy
Ec：<数值> 对应方法给出的 correlation energy
ETot: <数值> 对应方法的总能量

York

北大-陶豫 发表于 2021-12-2 12:44
如果想把能量分解为交换、相关以及其他贡献，应该使用能量分解，参考FAQ：

http://gaussian.com/faq1/

好的，谢谢老师

sobereva

你先明确你到底要做什么。库仑作用能和库仑积分、交换能和交换积分，完全是两码事

York

sobereva 发表于 2021-12-3 03:51
你先明确你到底要做什么。库仑作用能和库仑积分、交换能和交换积分，完全是两码事

好的老师，谢谢您的提醒

York

sobereva 发表于 2021-12-3 03:51
你先明确你到底要做什么。库仑作用能和库仑积分、交换能和交换积分，完全是两码事

老师，我想问您一下怎么样能够单独算出下图中公式里面说的两个未成对电子的交换能，我用楼上老师说的方法算的好像是整个体系的交换能

sobereva

York 发表于 2021-12-19 23:01
老师，我想问您一下怎么样能够单独算出下图中公式里面说的两个未成对电子的交换能，我用楼上老师说的方法 ...

Multiwfn目前可以基于数值格点计算特定两个轨道之间的库仑和交换积分，看手册3.200.17节。以后有可能会基于更快也更精确的解析积分计算。

hebrewsnabla

PySCF可以算。fcidump功能可以输出所有分子轨道之间的积分。也可以指定任意若干个轨道之间的。

York

sobereva 发表于 2021-12-20 18:05
Multiwfn目前可以基于数值格点计算特定两个轨道之间的库仑和交换积分，看手册3.200.17节。以后有可能会基 ...

好的，谢谢老师

York

hebrewsnabla 发表于 2021-12-20 18:38
PySCF可以算。fcidump功能可以输出所有分子轨道之间的积分。也可以指定任意若干个轨道之间的。

好的，谢谢老师

York

sobereva 发表于 2021-12-20 18:05
Multiwfn目前可以基于数值格点计算特定两个轨道之间的库仑和交换积分，看手册3.200.17节。以后有可能会基 ...

老师，我在文献中看到有人求轨道间的质心距离用了第一张图片的公式并引用了您的Multiwfn，但是我没有在Multiwfn手册找到这样的公式，我只看到第二张图片求质心距离的公式，请问老师这两个公式有什么区别吗？

北大-陶豫

York 发表于 2021-12-26 20:44
老师，我在文献中看到有人求轨道间的质心距离用了第一张图片的公式并引用了您的Multiwfn，但是我没有在Mu ...

基本算是大一普通物理的水平
想想两个质点的质心怎么求，三个质点的质心怎么求，很容易就可以推广到N个质点。
然后，分立的质点是用求和的方法计算质心，连续的质点自然是用积分的方法计算质心。
不明白的话可以找一本大一的普通物理教材看看。第一张图应该是求整个分子的质心以及求某种原子的质心，第二张图才是求轨道的质心

York

北大-陶豫 发表于 2021-12-26 21:01
基本算是大一普通物理的水平
想想两个质点的质心怎么求，三个质点的质心怎么求，很容易就可以推广到N个 ...

好的好的老师，谢谢老师