可能不少人对CASSCF下的ICSS和NICS感兴趣,我最近在MOKIT里实现了此功能,可以自动调用Dalton做CASSCF下的核磁计算、添加格点、汇总格点生成cub文件,cub文件可直接用Multiwfn/VMD/GaussView观看。由于ICSS计算量很大(Dalton程序内部有限制,导致这个像高斯DFT下ICSS那样要分很多文件算),目前只能算(约)20原子以下的小活性空间体系,可以用于有多组态/多参考态特征的体系做点研究。计算实例见PDF手册5.4 ICSS and NICS of Cyclobutadiene章节,示例输入文件见16-C4H4.gjf。使用该功能发文章时请记得在正文里引用MOKIT(另外Dalton程序也是要引用的)。如果觉得ICSS算不动,可以改算NICS(1)_ZZ之类的NICS值。
贴两张环丁二烯和C9H9+的图,前者CASSCF(4,4)/6-31+G(d)下算的,后者CASSCF(8,8)/6-31G(d)下算的,用的都是比较粗糙的格点。图是GaussView显示的。 (, 下载次数 Times of downloads: 18)