理论计算的赝势基组解析.docx

赝势基组【分享】?? 2009-06-01 15:37:08|??分类： HYPERLINK /ares_119/blog/ \l m=0t=1c=fks_087071093084080075092095082095085087084068092083087066087 \o 计算学习总结 计算学习总结 |??标签： |字号大中小?订阅 简单来说，赝势就是不计算内层电子，而是把内层电子的贡献用一个势来描述，放在哈密顿里面。适用于重元素。 赝势基组，实际上包括 赝势和基组两个部分，内层电子采用赝势，即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 opt lanl2dz for all atoms 0 2 O 0.0 0.0 0.0 C 0.0 0.0 1.2 Cu 0.0 0.0 3.2 -------------------------- 2.所有原子使用lanl2dz的另一种输入方法。 -------------------------- #HF/genecp opt lanl2dz for all atoms 0 2 O 0.0 0.0 0.0 C 0.0 0.0 1.2 Cu 0.0 0.0 3.2 C O 0 lanl2dz **** Cu 0 lanl2dz //定义价电子的基组, C O 0 是碳，氧，零，其中零用作终止符号。 **** Cu 0 lanl2dz //定义内层电子的赝势 -------------------------- 3.混和基组，即有的使用全电子，有的使用Lanl2dz。格式同2。 #HF/genecp opt lanl2dz for Cu, 6-31G(d) for C and O 0 2 O 0.0 0.0 0.0 C 0.0 0.0 1.2 Cu 0.0 0.0 3.2 C O 0 6-31G(d) //另一种全电子基组 **** Cu 0 lanl2dz //定义价电子的基组 **** Cu 0 lanl2dz //定义内层电子的赝势 -------------------------- ? ? ? ? ? ? PS:我的算例 ? ? %chk=XbridgeSDD%mem=50MW%nproc=1#p opt ccsd genecp C H F use 6-311++G(3df,3pd)?X use SDD 0 1?C???????????????? -2? -2?? 0?H???????????????? -2? -3?? 0?H???????????????? -3? -2?? 0?H???????????????? -2? -2? -0?X??????????????? -4? -2?? 0?F???????????????? -6? -2?? 0C H F 0?????? 6-311++G(3df,3pd)?????? ****?????? X?0?????? SDD ?????? **** X?0?????? SDD 赝势的一些东西 (2011-01-04