总结：高斯量化计算总结.docVIP

总结：高斯量化计算总结（必看） 丫85丫(金币+2,VIP+0):谢谢关于自旋多重度定义: 多重度＝2S+1, S=n*1/2,n为单电子数。所以，关键是单电子的数目是多少。当有偶数个电子时，例如 O2，共有16个电子，那么单电子数目可能是0，即8个alpha和8个beta电子配对，对应单重态，但是也可能是有9个α电子和7个β电子，那么能成对的是7对，还剩2个α没有配对，于是n=2,对应的是多重度3。同理还可以有多重度5，7，9， ...一般而言,是多重度低的能量低,最稳定,所以,一般来说,偶数电子的体系多重度就是1。但是也有例外，例如O2就是一个大家都知道的例子，它的基态是三重态，其单重态反而是激发态。所以对于未知的体系，还是算几个保险一点，看哪个能量更低。所以，总结一下，就是电子数目是偶数，未成对电子数目n=0,2,4,6,...自旋多重度是1,3,5,7,...电子数目是奇数，未成对电子数目n=1,3,5,7,...自旋多重度是2,4,6,8,...多数情况是多重度低的能量低，有时（特别是有“磁”性的时候，例如顺磁的O2，以及Fe啊什么的），可能会高多重度的能量低，所以需要都算算，看哪个能量更低。关于赝势：简单来说，赝势就是不计算内层电子，而是把内层电子的贡献用一个势来描述，放在哈密顿里面。适用于重元素。赝势基组，实际上包括 赝势和基组两个部分，内层电子采用赝势，即effective core potential (ECP)，外层价电子采用一般的基组。 比如：LanL2DZ: D95V on first row, Los Alamos ECP plus DZ on Na-Bi.就是对第一行原子是 D95V (这个是非赝势基组)，对Na-Bi是使用一个叫做Los Alamos的有效核势加上一个DZ基组。所以Lanl2dz就是对前面的原子全电子基组，对后面的原子是赝势基组。 （再次说明，量化里面，C，O那一行，算周期表的第一行）使用赝势的3个原因：1。没有相应的全电子基组。2。减少计算量。3。赝势可以包含重金属相对论效应的修正。在高斯中，lanl2dz基组，在手册中可以查到其定义为：LanL2DZ: D95V on first row , Los Alamos ECP plus DZ on Na-Bi，也就是说，对于C，O等元素来说（量化中，H大概算第0周期，C，O才是第一周期），Lanl2dz实际上还是全电子基组，而对于Na以后才是对内层电子用Los Alamos ECP赝势，外层电子用DZ 基组。使用赝势的输入文件：1.所有原子使用lanl2dz--------------------------#HF/lanl2dz optlanl2dz for all atoms0??2O??0.0 0.0 0.0C??0.0 0.0 1.2Cu 0.0 0.0 3.2--------------------------2.所有原子使用lanl2dz的另一种输入方法。--------------------------#HF/genecp optlanl2dz for all atoms0??2O??0.0 0.0 0.0C??0.0 0.0 1.2Cu 0.0 0.0 3.2C O 0lanl2dz****Cu 0lanl2dz? ? //定义价电子的基组, C O 0 是碳，氧，零，其中零用作终止符号。****Cu 0lanl2dz? ?//定义内层电子的赝势--------------------------3.混和基组，即有的使用全电子，有的使用Lanl2dz。格式同2。--------------------------#HF/genecp optlanl2dz for Cu, 6-31G(d) for C and O0??2O??0.0 0.0 0.0C??0.0 0.0 1.2Cu 0.0 0.0 3.2C O 06-31G(d)? ?//另一种全电子基组****Cu 0lanl2dz? ? //定义价电子的基组****Cu 0lanl2dz? ?//定义内层电子的赝势--------------------------开壳层和闭壳层闭壳层计算就是对于多重度是1的体系，此时α和β的电子数目相同，可以把α和β配对，成对的α和β使用同一个轨道，一个轨道上填充2个电子。开壳层计算就是对α和β电子分别计算，一个轨道上只填充1个电子，一般来说，多重度是1时，开壳层计算和闭壳层计算会给出相同的结果。限制开壳层计算是对多重度大于1的体系，此时α和