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（及其变体）成分的。请问这种做法是否对弱相互作用等问题有一定价值？

社长威武！