请下载ms excel文件 [ Last edited by cadick on 2009-12-14 at 01:24 ]
评价为good/medium的有这几个程序:gamess,molcas,molpro,pcgamess,orca,utchem。后三个程序都采用了gamess的LS方法,因此可归为一类。 本人一直在用前四种程序做MRPT2计算,迄今为止算过大约30种分子,如果连测试也算上(比如,用不同的基组,不同的活性空间,不同的参考波函来源,不同的LS值,不同的构型,...),估计要好几百个样本。我的经验是,molpro的LS对一半的样本解决不了is,故评价为medium。换成gamess/pcgamess/molcas后,95%的情况下不会有is,故评价为good。但是在扫描势能面的时候,个别构型仍可能有解决不了的is,这时候可以把它们跳过去,通过插值方法进行估算,因此不影响结果。 另外,所有已发表的LS文献几乎都是molcas和gamess程序作者做的,反映出这两个程序在解决is问题上有优势,
在Molcas中,可以打印基态和激发态的reference weight, 如果相差很大可能是intruder state的问题,通过shift使基态和激发态的weight接近,从而解决intruder state的问题。我不知道其他程序包是什么情况的,评估其他程序包也是通过reference weight判断的? 还是说发现计算出的能量异常的低,所以可能是分母过小,微扰计算趋于发散?或者其他判断方法?
molpro:有入侵态自动停止计算 gamess/pcgamess/utchem:打印参考权重
了解,谢谢!
谢谢诸位高手。
xie xie
搞量子化学计算的好多呀 支持楼主 哪天也出来个关于凝聚态物理计算的呵呵
评价为good/medium的有这几个程序:gamess,molcas,molpro,pcgamess,orca,utchem。后三个程序都采用了gamess的LS方法,因此可归为一类。
本人一直在用前四种程序做MRPT2计算,迄今为止算过大约30种分子,如果连测试也算上(比如,用不同的基组,不同的活性空间,不同的参考波函来源,不同的LS值,不同的构型,...),估计要好几百个样本。我的经验是,molpro的LS对一半的样本解决不了is,故评价为medium。换成gamess/pcgamess/molcas后,95%的情况下不会有is,故评价为good。但是在扫描势能面的时候,个别构型仍可能有解决不了的is,这时候可以把它们跳过去,通过插值方法进行估算,因此不影响结果。
另外,所有已发表的LS文献几乎都是molcas和gamess程序作者做的,反映出这两个程序在解决is问题上有优势,
在Molcas中,可以打印基态和激发态的reference weight, 如果相差很大可能是intruder state的问题,通过shift使基态和激发态的weight接近,从而解决intruder state的问题。我不知道其他程序包是什么情况的,评估其他程序包也是通过reference weight判断的? 还是说发现计算出的能量异常的低,所以可能是分母过小,微扰计算趋于发散?或者其他判断方法?
molpro:有入侵态自动停止计算
gamess/pcgamess/utchem:打印参考权重
了解,谢谢!
谢谢诸位高手。
xie xie
搞量子化学计算的好多呀 支持楼主
哪天也出来个关于凝聚态物理计算的呵呵