关于LC-ωPBE计算的几个问题

小范范1989

sobereva 发表于 2015-11-28 09:20

谢谢sob老师的指点。谢谢老师

卡开发发

1、只是人为划分为远近而已，可以认为是erf函数性质本身导致的，选取这个erf的目的也是如此，至于远或近对r12贡献的大小是通过omega调整的，你自己可以使用Excel看看，也内置了erf函数；
2、（Gaussian的iop我不熟悉，就不班门弄斧了）；
3、范围分离泛函其实也就是把长程项和短程项用不同的泛函描述，如果使用两种局域或半局域泛函，一般对计算通常没改善，所以基本上也就是一种杂化泛函，杂化泛函这样就能够分成全局和范围分离的。全局的杂化泛函长程短程HF比例都是一样的，这个你的理解没问题；
4、LC-GGA一般LR的部分是100%的HF，能修改的只有omega项。

当然，长程短程的形式也不是绝对纯粹erf，比如CAM-B3YP，短程项就是[1-α+β·erf(μr12)]/r12，长程则是[α+β·erf(μr12)]/r12。其实2#李大师分享的那份ppt当中Scuseria对三类如LC-wpbe、HISS、HSE泛函还是做了些详尽的介绍。

小范范1989

卡开发发 发表于 2015-11-28 11:53
1、只是人为划分为远近而已，可以认为是erf函数性质本身导致的，选取这个erf的目的也是如此，至于远或近对r ...

谢谢老师的指点。

小范范1989

咨询一下两位老师， @sobereva  @卡开发发  这个LC-wpbe泛函怎么在高斯09d01中实现，
我使用：#p lc-wpbe/6-31g(d) iop(3/107=0203600000,3/108=0203600000)  opt TDA scrf=(solvent=toluene)
结果：
The selected state is a singlet
CISAX will form     1 AO SS matrices at one time.
NMat=     1 NSing=     1 JSym2X= 0.
No func 3rd derivs with HSE.
Error termination via Lnk1e in /home/software/g09D01/g09/l914.exe at Wed Jan 27 15:17:00 2016.
Job cpu time:       0 days  0 hours 30 minutes 24.2 seconds.
No func 3rd derivs with HSE.
Error termination via Lnk1e in /home/software/g09D01/g09/l914.exe at Wed Jan 27 15:17:00 2016.
No func 3rd derivs with HSE.
Error termination via Lnk1e in /home/software/g09D01/g09/l914.exe at Wed Jan 27 15:17:00 2016.
No func 3rd derivs with HSE.
Error termination via Lnk1e in /home/software/g09D01/g09/l914.exe at Wed Jan 27 15:17:00 2016.
File lengths (MBytes):  RWF=    249 Int=      0 D2E=      0 Chk=      6 Scr=      1

也就是这个方法不能来做s1的优化？那文章中的s1-t1的能差是咋找的嫩？谢谢老师

sobereva

小范范1989 发表于 2016-1-27 15:12
咨询一下两位老师， @sobereva  @卡开发发  这个LC-wpbe泛函怎么在高斯09d01中实现，
我使用：#p lc-wpbe/ ...

都有lc-wpbe关键词，还写那些IOp干嘛
有那个提示说明那个泛函在g09里有局限性，没法得到解析梯度。要优化只能用数值梯度了，会很慢。

小范范1989

sobereva 发表于 2016-1-27 17:17
都有lc-wpbe关键词，还写那些IOp干嘛
有那个提示说明那个泛函在g09里有局限性，没法得到解析梯度。要优 ...

恩，明白了，谢谢sob老师。