Sob老师好，我看了您上面发的教程，在实现的过程中遇到了几个问题：
1.“如果wLong或wShort其一为0，本质上就等于是wLong=wShort的情况，即不存在中程部分，此时也相当于只定义了一个w。”我在尝试的过程中发现当我使用只有长短程的范围分离泛函，定义wLong=wShort=w和定义wShort=0，wLong=w时算出来的HOMO,LUMO轨道能和分子单点能不一样，请问这是为什么
我是这么写输入文件的：
A.定义wShort=0，wLong=w
#p wB97X/def2TZVP
scrf=solvent=Dichloromethane
IOp(9/40=4)
IOp(3/107=0178000000,3/108=0178000000)
IOp(3/119=0112001577)
IOp(3/120=0888000457)
IOp(3/131=100000)
B.定义wLong=wShort=w
#p wB97X/def2TZVP
scrf=solvent=Dichloromethane
IOp(9/40=4)
IOp(3/107=0178001780,3/108=0178001780)
IOp(3/119=0112001577)
IOp(3/120=0888000457)
IOp(3/131=100000)

2.“cFull+cLong对应长程系数，cFull+cShort对应短程系数。对于存在中程的情况，cFull相当于中程的系数。”
我想实现附图中的公式，则短程DFT的系数应该为1-1/ε，对应DFx的cFull+cShort
但我发现我设置cFull=0.8423，cShort=0和设置cFull=0，cShort=0.8423的计算结果（单点能）并不相同，这是否说明短程系数并不对应于cFull+cShort？
cFull=0.8423，cShort=0： cFull=0.8423,cShort=0.txt (92 KB, 下载次数 Times of downloads: 0)

2025-9-18 11:50 上传 Uploaded

点击下载
Click to download

IOp(3/107=0178001780,3/108=0178001780)
IOp(3/119=0112001577)
IOp(3/120=0045700000)
IOp(3/131=08423)
cFull=0，cShort=0.8423： cFull=0,cShort=0.8423.txt (88.72 KB, 下载次数 Times of downloads: 0)

2025-9-18 11:49 上传 Uploaded

点击下载
Click to download

IOp(3/107=0178001780,3/108=0178001780)
IOp(3/119=0112001577)
IOp(3/120=0888008423)
IOp(3/131=-1)

sob老师好，我看了上面的教程，尝试用另外一种方式定义LC-PBE0，
#P PBEPBE/6-31G*  IOp(3/76=1000010000,3/77=1000010000) IOp(3/119=1000002500,3/120=1000002500) IOP(3/107=0300000000,3/108=0300000000) IOP(3/130=-1,3/131=-1)
得到的成分cshort=0.25,cfull=0.0,clong=1.00，同样满足cshort+cfull=0.25,clong+cfull=1.00的关系，但是算出来的能量和上面您定义的LC-PBE0不完全一致，差0.05hartree左右（HCHO分子），是哪里我理解错误的？
更新：我想了一下，可能是因为原始定义里只有beta是跟w相关的，这样是不是相当于把跟w相关的部分弄错了应该？

Masshiro 发表于 2021-2-25 15:18
想请问下，对于LC-PBE0中α=0.2、α+β

这种情况没意义，这样就成了长程比短程HF成份还低，违背了初衷

想请问下，对于LC-PBE0中α=0.2、α+β<1.0的情况下，如果α+β很小导致β为负数时，IOp(3/119)和IOp(3/120)要怎么写呢

linq5203 发表于 2020-6-22 18:12
请问sob，文中提到的对DFT相关部分进行范围分离的iop(3/109,3/121,3/132)是不是可以这么用？在用双杂化泛函 ...

相关部分的范围分离我没试过，虽然有这个选项，不代表程序真正的实现了或者被充分测试过
这对于改进弱相互作用我认为没什么意义，双杂化泛函结合色散校正就已经足够理想了，不会有什么额外的油水

请问sob，文中提到的对DFT相关部分进行范围分离的iop(3/109,3/121,3/132)是不是可以这么用？在用双杂化泛函时，如果搭配iop(3/126,3/127,3/128,3/129,3/133,3/134)等对MP2的自旋平行/反平行进行范围分离的iop，理论上是可以在长程部分使用100%MP2（及其变体）成分的。请问这种做法是否对弱相互作用等问题有一定价值？

社长威武！