计算化学公社

标题: 量子化学中的冻芯近似是怎么实现的？ [打印本页]

作者
Author: 清微 时间: 2018-1-4 16:46
标题: 量子化学中的冻芯近似是怎么实现的？
量子化学中的冻芯近似是怎么实现的？是对内层电子不定义基组还是自洽场迭代时把内层电子的影响归结到核电荷上?

作者
Author: Warm_Cloud 时间: 2018-1-4 20:24
本帖最后由 Warm_Cloud 于 2018-1-4 20:26 编辑

冻芯近似是给post-hf方法用的，正常做完HF后，在做posh-hf的时候，以前求和是从1到Nocc（占据轨道数），使用了冻芯近似后，改成FCstar+1到Nocc，FCstar是你冻结的轨道数，举个例子，对于CO2,一共22个电子，Nocc就是11，对于C和O，分别冻结2个电子，总共就是冻结6个电子，3个占据轨道，因此，求和由原来的1到11，变成4到11，你可以看这段代码：

Tempc = 0.0d0
Tempk = 0.0d0
do k = FCstar+1,Nocc
   do l = FCstar+1,Nocc

      do a = Nocc+1,Nbas
      do b = Nocc+1,Nbas
         Temp = (En(k)+En(l)-En(a)-En(b))
         Tempc = Tempc + g_Mo(EI(b,l),EI(a,k))*g_Mo(EI(l,b),EI(k,a))/Temp
         Tempk = Tempk + g_Mo(EI(b,l),EI(a,k))*g_Mo(EI(l,a),EI(k,b))/Temp
      end do
      end do
   end do
end do

这是做MP2的代码，不使用冻芯近似的时候，FCstar=0，使用的时候，FCstar就是每个原子的冻芯轨道数目加起来。上面的例子就是1（C）+1（O）+1（O）=3
另外，你说的在SCF中的有类似思想的那叫赝势。

作者
Author: lindlar 时间: 2018-1-5 02:35
学习啦！！！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)