当前位置:大红鹰娱乐8055 > 社会科学 >

Gaussian98使用经验谈

发布时间:2017-12-03 阅读:

  Gaussian98谈论使用经验

  经过一段时间的使用高斯,不仅加深了对量子化学的认识,而且巩固了光谱学的理论知识。高斯真的是一件很好的事情,至少在研究光谱理论方面。使用中也有很多问题,并没有包含在高斯的帮助文件中,其中一些问题是由他们自己解决的,有些是与他人讨论的,有些是在网上或文献中找到的。研究方向,涉及到的大部分文章都涉及到分子激发态的问题,如果对分子光谱学(包括红外,拉曼,可见 - 紫外光等)缺乏了解,建议先看对称和光谱学由Harris和Bertroux出版,高等教育出版社,1988),一本简单易读的快速出版书。我对有些问题的理解水平有限可能是错误的,欢迎批评和纠正。 ***************目录******************(一)各种杂项的一些小问题§1-1。内存和速度问题§1-2.GAUSSIAN输入一个注释§1-3.GAUSSIAN缺陷(b)基本集合设置§2-1。基本集合§2-2。基本集合选择(c)GAUSSIAN e xcited状态处理§3-1。输入文件电子状态重复度§3-2。高斯在兴奋状态方法研究§3-3。包含时间的方法(TD)输入方法§3-4.CIS,TDHF,TDDFT方法激发能量比较§3-5.GAUSSIANBUG§3-6处理激发态图3.激发态的势能曲线§3- 7。使用CIS方法的高对称分子(iv)找到分子的光谱常数§4-1。寻求解离能§4-2。寻求电离能§4-3。寻找激发态Te,ωe,Re4-4的谱参数。如何获得谱项的对称性[附录]高斯98的补充1.共振强度f2。 TDDFT做几何优化3.自定义赝势(模型潜能,有效核潜能)(一)各种杂项的一些小问题§1-1。内存和速度问题G98W的默认内存是128兆,我以前用的是64兆,电脑是P233,所以非常慢。我后来发现,在计算之前关闭所有防病毒程序(我的计算机上有三个)可以节省20%的时间。然后加了64兆的内存,时间缩短到了原来的一半。如果将电脑换成P300,用128兆内存,只需要五分之一的时间。大内存支持更大的基本集合。因此,电脑配置可能马虎。 §1-2。高斯输入注意事项输入分子坐标(包括键长和键角)必须是实数,不能是整数,否则将是错误的。例如,密钥长度是2.0埃,可以写成r = 2或r = 2.0,不能写成r = 2。§1-3。高斯的缺点对于含有重原子的分子,光谱的自旋轨道角动量已经失去了意义。此时的光谱项应使用自旋轨道耦合的总角动量中的每一个计算。分量Ω意味着。例如,氯化铊分子的3Π电子态是作为两个非常开放的Ω分量实验观察到的:A0 +态和B1态。但是GAUSSIAN只能处理含有H原子到Cl原子的分子(因为它们不是重原子,所以不需要)。目前,这些问题通常通过免费的MOLFDIR,GAMESS(美国),NWChem,DIRAC,COLUMBUS或商业GAMESS(英国),MOLPRO,ADF等程序来解决。 (二)基本集合§2-1。基本组有两个基本组,一个是全电子基组,另一个是价电子基组。价电子基为内电子逼近含有相对论效应的赝势(缩写PP,也称为模型势,有效核势,缩写为ECP)。在铟原子的情况下,输入以下命令以观察3-21G组铟原子的形式:#3-2​​1G GFPrintIn0,2In输出结果是:3-21G(6D,7F)S 3 1.00 N = 1层,轨道.1221454700E + 05.6124760000E-01.1848913600E + 04.3676754000E + 00.4063683300E + 03.6901359000E + 00SP3 1.00 N = 2层,s和p轨道.5504422550E + 03-1.1107094330E + 00.1523702990E + 00.1197743540 E + 03 .8344349530E-01 .6096507530E + 00.3866926950E + 02 .9696880360E + 00 .3970249500E + 00SP 3 1.00 N = 3层,s和p轨道。 4702931320E + 02 -.2758954250E + 00-1.4808484750E + 00.2249642350E + 02.5977348150E-01.5290866940E + 00.6697116970E + 01.1082147540E + 01.66620681110E + 00D 3 1.00 N = 3层d轨道.1021735600E + 03。 1205559000E + 00.2839463200E + 02 .4884976000E + 00.8924804500E + 01 .5850190000E + 00SP 3 1.00 N = 4层,s和p轨道.6572360380E + 01 .4284830560E + 00 .1091304620E-01.2502157560E + 01 -.4633643610E + 00 .5036758870E + 00.9420245940E + 00 -.8219679320E + 00 .5581808650E + 00D 3 1.00 N = 4层,d轨道。 4535363700E + 01 .2508574000E + 00.1537148100E + 01 .5693113000E + 00.4994922 600E + 00 .3840635000E + 00 \\ n \\这些是内部电子轨道,每个轨道都有三个高斯函数拟合,即在3-21G中3. SP 2 1.00 N = 5层,s和p轨道。 1001221380E + 01-4.364171950E + 00 -.2316333510E-01.1659704190E + 00 .1189893480E + 01 -.9903308880E + 00这是每个轨道的价电子轨道两个高斯函数拟合,即3-21G中的2。 SP 1 1.00 N = 6层,s和p轨道5433974090E-01 .1000000000E + 01 .1000000000E + 01 \\ n \\这是一个空轨道,每个轨道都装有一个高斯函数,即3-21G中的1个。 \\ u0026可见3-21G是全部电子基础组。让我们看看InL原子LanL2DZ基本集合的结果:LANL2DZ(5D,7F)S 2 1.00.4915000000E + 00 -.4241868100E + 01.3404000000E + 00 .4478482600E + 01S 1 1.00.7740000000E-01 .1000000000E + 01P 2 1.00.9755000000E + 00 -1226473000E + 00.1550000000E + 00.1041757100E + 01P 1 1.00.4740000000E-01 .1000000000E + 01这些都是价电子轨道。 ECD中添加了双Zeta(DZ)基本集。所以LanL2DZ是一个价电子基础集合。 \\ u0026在高斯里,CEP也是一个非常重要的价电子基础集合,第98版元素扩展到元素86 Rn。请注意,每个元素的CEP-4G,CEP-31G和CEP-121G与Ar元素相同(但不同于CEP +极化*)。 SDD适用于除Fr和Ra之外的所有元素,但结果的准确性不是很好。 \\ u0026请注意,有些方法不需要输入基本集,如各种半经验方法,G1,G2等。 §2-2。基础设置的选择中国科学B30(5),2000,第419页中的文章给出了选择基础集合的方案。基函数的选择应根据系统中不同原子的性质和实际的化学环境:1.根据元素周期表中的原子位置从左到右描述一般体系,以增加基数功能。 2.负核原子,使用更多的基函数和适当的偏振函数或分散函数。对于核正原子,可以适当减少基函数。 3.对于正常的饱和共价键原子,不需要增加极化或色散功能,氢键,弱相互作用体系,官能团,零价或低价金属原子等敏感体系,增加极化或色散功能。 \\ u0026这使得使用相当完善的基本集合和合理可靠的计算成为可能。程序为HF,MPn,DFT和其他计算。 \\ u0026 GAUSSIAN中的关键字是ExtraBasis和Gen.两者都是基于路径指定的基本集合,加上附加的基础函数。不同之处在于Gen可以定义自己的基本集合。其实这可能不是那么简单。以下是使用ExtraBasis的一个例子:#HF / 3-21G ExtraBasisAgCl分子规格C1 06-311G(d,p)**** \\ u0026 nbsp;这里所有分子中的原子都用3-21G,但是用较大的原子Cl基团6-311G(d,p)。如果计算错误,通常由于以下两个原因:1.不是所有的组都可以使用(哪个组适用于帮助文件所需的原子组),即使每个组被应用于不同的组原子,行不一定在一起,可能有匹配或收敛的问题,这可能是BUG。所以需要反复测试。 2.基本问题的顺序。仍然看上面的例子,我认为在路线部分中设置的碱基应该适用于分子中的所有原子(一般来说,它稍微小一些,尽管基本集合不用于所有原子的出现) (例如,如果上面的输入改变为:#HF / 6-311G(d,那么计算)),而在分子坐标下面的基底集合是可以的(通常是较大的基底集合) p)ExtraBasisAgCl分子说明书Ag 03-21G **** -311G(d,p)的基础设定,但可能会报错。当然,实际情况并非总是如此,另一个分子, (三)高斯处理激发状态§3-1。输入文件中的电子状态多重性输入中的多重性文件是指这种多重性能量最低的电子状态,不是必需的ily状态,但是,一般输入文件与地面状态。激励状态的多样性由输入文件的路径部分(如单倍数,三元组)控制。虽然在高斯的帮助中没有提到,但从所给的例子中可以明显看出。另外,在HyperChem的帮助文件中有一个类似的说明,尽管它是用于HyperChem的,但是我希望GAUSSIAN适用。 §3-2。在高斯高斯研究激发态的方法激发态可以用CIS,CASSCF和ZINDO方法计算。 GAUSSIAN中的CIS实际上是一种零级配置相关的处理分子的激发态。一些从头计算器,如GAMESS,包含用于计算激发态的一级配置相关(FOCI)和二级配置相关(SOCI),我不知道这是否与CIS相对应。从结果来看,FOCI和SOCI略有改善。 ZINDO方法的适用范围很窄,只能用于由周期表的前48个原子组成的分子,这里不再赘述。版本98还添加了一个新的时间依赖方法来处理激发态,包括时间依赖的仇恨 - 伪装(TDHF)和时间相关密度泛函(TDDFT)。我们一般计算激发态称为低激发态,它由原子价电子分子激发态组成。更高的电子态是里德堡状态。从头算法和其他半经验方法对轨道的计算,特别是高海拔轨道的计算都不太好,得到的里德堡态的结果缺乏明确的物理意义。 §3-3。基于时间的方法(TD)输入方法帮助文件对于TDDFT是非常简单的,作为参考CIS的输入方法的一个重要方法,甚至没有一个例子。 CIS的输入法是:#CIS(Triplets,NStates = 10)/ 3-21G Test然后使用B3LYP功能测试:#TBB3LYP(Triplets,NStates = 10)/ 3-21G Tests 实际上,正确的输入法是:#B3LYP / 3-21G TD(Triplets,NStates = 10)测试这个问题浪费了我很长时间,为了说明它的重要性,分成了几个部分。 -4.CIS,TDHF,TDDFT方法激发能量比较一般来说,CIS和TDHF方法结果的准确性相对较差,TDDFT结果最好,但B3LYP,B3P86和MPW1PW91的效果最好各种TDDFTs(取决于所使用的特定分子和碱基集),Physiol.Phys.109(23),1998,10180中;这些文章计算了构成的有机物的垂直激发能我们计算了重元素In的分子激发态,发现如果选择基态,许多TDDFTs的结果与实验值之间的差异甚至小于500波数,而CIS结果可能比5000波数更差。所以现在更多地使用TDDFT方法。另外还有人用TDDFT方法来计算这样一个大分子的DNA。 §3-5.GAUSSIAN处理激动人心的BUG1。对于CIS版本94的计算,50-50选项不支持LanL2DZ和3-21G *基本集合。 98版本不会出现这个问题。这个问题可能与计算的分子有关。 2.同时计算单重激发态和三重激发态(即使用50-50)有时会导致致命的错误:某些状态在某些位置有不规则的能量波动,奇点或奇点三重态(使用单数或三重选项)不显示这个问题。我计算了两次。下面的例子是从互联网上看到的,我觉得德国的兄弟遇到了同样的问题。 \\ u0026他对同一个任务进行三次计算:1)#Bp86 / 6-31G * td =(nstates = 5,50-50)scf =直接密度2)#Bp86 / 6-31G * td = nstates = 1)scf = direct密度 - link1 - #Bp86 / 6-31G * td =(nstates = 10)scf =直接猜测=读取geom = check - link1 - #Bp86 / 6-31G * td =(triplets,nstates = 10)scf = = readgeom = check3)#bp86 / 6-31G * td =(nstates = 10)scf = direct#Bp86 / 6-31G * td =(triplets,nstates = 10)scf =直接密度第一次计算产生10个状态和5T),后两者产生20个。关键字密度对结果没有影响。结果是:3)三元组-B3U 3.1077 3.1077 3.1077三元组B3U 3.9568 3.9568 3.9568三元组AG-4.1612 4.1612 4.1612三元组-B3U 4.1769 4.1769 4.1769三元组-B3G 4.2101 4.2101 4.2101小星B3G 4.3609 4.3609 4.3609小星-B1U 4.3788 4.4864 4.4864小星 - B2U 4.4864 4.6409 4.6409Triplet-AG 4.5752 4.7144 4.7144 Singlet-B3G 4.6409 5.0010通过比较可以发现,对于前六种状态,三种处理的结果是相同的,但在后一种状态中, 50方法与其他两种方法不同。 \\ u0026我个人认为,尽量不要使用这个选项50-50。推荐单线,三线分开计算,虽然花了两倍多的时间,但保险的结果。 §3-6。绘制激发态势能曲线最愚蠢的方法是在完成一个点后,通过改变输入文件的坐标来计算下一个点。这是费时费力的。扫描可以在优化基态几何时使用Scan关键字完成,扫描实际上就像激发态一样工作。所以,打开电脑过夜,到整个早晨数据的兴奋状态的所有能量曲线。例如,计算氢分子从R = 0.6埃至1.2埃的10个单重激发态:#RCIS(NStates = 10)/ 6-311G ScanH2激发态0,1H1H2 H1 rr 0.6 60 0.01需要注意的三件事:结果是垂直激发能量必须转换成能量相对于基态的平衡状态。 2.不要从基态平衡的范围扫描太远,一般从0.75Re到2Re就可以了,太远就会出错。 3.一些理论和一些基本的小组搭配可能不支持扫描,这个时候不得不手动控制。 \\ u0026最后是数据处理,这一步是最无聊的。数据处理需要使用Excel,Oringin,Matlab。这些软件不需要深入学习,知道一些最基本的命令就足够了。 §3-7。使用CIS方法的高对称分子默认的CIS计算找到最低的三个激发态。以此为例,GAUSSIAN通过选择两个状态数(即六个)的向量来迭代初始猜测,直到有三个收敛为止。对于大多数分子系统来说,这个默认的向量数是足够的但是对于分子的高对称性,你需要特殊的程序,否则你不会得到更高的收敛状态。该方法是增加N状态的数量来增加初始向量的数量。建议NStates值是Abel组操作的最大数量,即Standardorientation之前的输出文件的对称部分(可以通过在输入文件中使用#控制打印输出而不是#T来获得)。这是一个计算苯分子的例子:化学计量C6H6Framework group D6H [3C2(HC,CH)] Deg。自由度2Full point group D2H NOp 8最简洁的Abelian子群D2 NOp 4标准方向------------ -----------------------------------------因此在路段中使用CIS(NStates = 8)。 (D)寻找分子光谱常数§4-1。寻求解离能首先我们要注意的是,原子化能和解离能不是一个概念。帮助文件和电子结构方法探索化学的解释不清楚。计算原子能和解离能是两个不同的过程。以三原子分子GaCl2为例,计算解离能的过程是:GaCl2(X2A1)---- Cl(2P)+ GaCl(X1Σ+)过程中的能量变化称为解离能。显然,解离能又可细分为第一步解离能,第二步解离能......,以上是第一步解离能的计算,第二步解离能的过程是:GaCl(X1Σ+ )---- Cl(2P)+ Ga(2P)计算雾化能的过程是:GaCl2(X2A1)---- Cl(2P)+ 2Ga(2P)之间还有一个区别。但是对于双原子分子,这两个概念是相同的。第二点要注意的是,解离能一般可以用符号D0或De来表示(前者不包括零点能,后者包括在内),但在分子光谱中也有一组D0和德。为避免混淆,需要特别指示。第一组D0和De,即解离能,值相对较大,数量级一般为0.1〜1eV。该参数出现在真空紫外振动光谱中。有一个公式:ωe^2ωeχe= ------- \\ u0026 \\ u0026 \\ u0026 \\ u0026 (^是指数符号)4De其中De是解离能。 \\ u0026另一组D0和德称为离心失真常数,数值较小,作为旋转谱的摄动。如: \\ u0026 \\ u0026 \\ u0026 \\ u0026 \\ u0026 \\ u0026 h ^ 3De = ------------------------------ ----- \\ u0026 nbsp; \\ u0026 128(π^ 6)(μ^ 3)(c ^ 3)(re ^ 6)(ωe^ 2)其中De是离心变形常数。 \\ u0026高斯有一个寻求PH2分子原子化能量的例子,可以作为求解离解能量的参考。这是文件e7-01。结果是雾化能= E(P)(Hartree)+ 2E(H)(Hartree)-E(PH2)(Hartree)-ZPE(Hartree)* 2 - ( - 342.50942) - (0.01322)(Hartree)零点ZPE被校正。相差148.3千卡/摩尔和144.7千卡/摩尔的实验值相当小。 §4-2。寻求电离能量电离能量简写为I.P.有两种方法可以找到I.P.在高斯。方法一:参见e7-03文件,本例寻求PH2分子电离能。结果是:ZPE PH2 + -342.14416 0.01347 PH2 -342.50942 0.01322(这是上述原子能的结果)IP =(E(PH2 +)+ ZPE(PH2 +)) - (E(PH2)+ ZPE(PH2 ))\\ u0026 nbsp; =(-342.14416 + 0.01347) - ( - 342.50942 + 0.01322)(Hartree)\\ u0026> = 0.36551(哈特里) = 9.95(eV)方法二:OVGF方法。帮助文件有详细的使用说明。 §4-3。求光谱参数Te,ωe,Re的激发态Re有两种方法。 CIS可以使用结构优化的方法。 Starl的文章讲得非常详细,这里就不详细说明了,只是补充一点,就是如何在设置后强制保留wave:使用Guess = Alter设置wave之后,程序会一直返回此时需要使用SCF关键字中的Symm选项来强制程序在设置后保留波形字符,方法是在路径部分添加SCF = Symm 2.潜能曲线扫描方法,TDDFT方法不能计算激发态的波函数,不能用于研究激发态的平衡构型和分子性质,上述结构优化方法不适用于TDDFT,但是可以利用势能但是多原子分子的能量是一个势能面,很难控制多个输入坐标变量对整个表面进行扫描,因此将这种方法应用于双原子分子。激发态的电势曲线可以在上面的§3-6中找到。利用势能曲线,我们可以找到势能Re的最低点(为了精确起见,在这个区域中的点被认为是少数)。 Re位置处的能量与基态平衡位置处的能量之差就是电子高度Te。势阱的深度是离解能量。 \\ u0026这些都很简单,关键是找到振动参数ωe。求ωe有两种方法。一种是使用莫尔斯势函数拟合的方法。 Morse势函数的形式如下:V(R)= Te + De {1-EXP [-α(R-Re)]} 2(对于基态Te,Te为0)在α之后,可以通过使用公式α=ωe((2(π^ 2)c *μ)/(h * De))0.5(高斯系统)。 \\ u0026这个方法太麻烦了,我们不这样做。只要扫描的莫尔斯势函数为二阶导数R,则导数f2与ωe的导数之间的关系为:4(π^ 2)(ωe^ 2)c *μf2= -h *(10 ^ 18)其中c,h,μ是国际单位制,ωe单位是波数cm ^ -1。另外,在电位扫描中,De的单位是波数cm -1,R和Re的单位是A.否则,上述的f 2需要乘以一个因子。 \\ u0026找到f2 Origin软件就可以完成了。注意采取Re f2的位置。用f2可以计算ωe,不需要计算De和拟合α步骤来避免引入更大的误差。 §4-4。如何获得光谱项的对称性如果用高斯函数来计算较轻的分子,通常给出激发态光谱的正确对称性,但是对于含有重原子的分子,对称性可能不会被获得。这次有必要根据轨道转换信息来获得对称性。 \\ u0026以下是使用CIS方法LanL2DZ基组计算的InCl分子的例子。 (PI)(PI)(PI)(PI)(PI)(PI)(PI)(PI)(PI)(SG)(SG)(SG)输出文件给出了InCl的电子轨道信息。 )这里只给出前4个单重态激发态:CIS波函数对称性不能确定。退出状态1:? Singlet-Sym 4.9226 eV 251.87 nmf = 0.32435 - 6.435875 - 7.59909这个状态用于优化和/或二阶校正。复制这个状态的Cisingles密度为1-particleRhoCI密度.CIS波函数对称性不能确定。预期状态2:单线 - ? Sym 4.9226 eV 251.87 nmf = 0.32435-6.549095-7.43587激发态3:单线态SG 6.7609eV 183.38nmf = 0.28303-6.280384-7.280385-8.57818CIS波函数对称性不能确定。退出状态4:单线 - ? Sym 6.8241 eV 181.68 nmf = 0.00003-7.4496904-6-6.49690因为In原子是重原子,所以大多数电子态不被给予对称性(第三排除态)。我们发现第一个激发态从第五个轨道转换到第六个和第七个轨道,而第六个和第七个轨道是退化的(PI和DELT轨道都是退化的)。换句话说,电子结构的第一激发态是:InCl :(内电子)(1σ)2(2σ)2(1π)4(3σ)1(2π) (当然还有一个三重态)谱项:1Π,这就是1Π态。再往下看,第二激发态和第一激发态的电子结构完全相同,说明这是一个退化态,事实上,Σ态以外的态其次是退化态。因此,第一激发态和第二激发态实际上都是InCl的第一激发态1Π态。 \\ u0026第三激发态给出对称性。如果你不这样做,就像上面所做的那样,找到它的电子结构,它是两种结构的迭层结构:InCl:(1σ)2(2σ)2(1π)3(3σ)2(2π)1InCl :(内电子) (1σ)2(2σ)2(1π)4(3σ)1(4σ)1第一种配置可以耦合出单峰谱:1Δ,1Σ+,1Σ-。因为这不是退化状态,所以只能是1Σ+和1Σ-。这需要进一步的+, - 对称知识。方法是看过渡系数:3 - 6 - .280384 - 7 - .28038 - 由于符号相同,所以是1Σ+状态。那么第四激发态必须是1Σ-。另外你可以看到第二种配置出来的频谱只有1Σ+,所以第三激发态必须是1Σ+。两个结果是一样的。 \\ u0026以上是自旋极限的计算。无限制的计算(每个轨道细分为自旋阿尔法和自旋β)分析是相似的。 ************************************************** ************************************************** ************************************************** ***********************************高斯98的补充1.共振强度f出现共振频率f激发态结果的计算,激发态3:单态 - Sym 5.0236 eV 246.80 nm f = 0.3428 5 - 7 0.68918在早期版本的Gaussian 98中,TD方法计算激发态的共振强度f是错误的(另一种计算激发态的方法没有这个问题) 。 A7版本从一开始就进行了修正。这一点应该注意。 \\ u0026共振强度可以做什么?预测电子激发跃迁的寿命。使用公式t = 1.499 /(f * E ^ 2)\\ u0026其中E是电子状态之间的能隙,单位是波数,f是电子状态的共振强度,t是以秒为单位的电子状态寿命。默认的较低级别是基态。如果我们计算激发态之间的f并在路径段中添加keynsnsDensities,则获得电子状态之间的偶极矩m(x),m(y)和m(z)。 f的解由下式给出:f =(2/3)* E * [m(x)^ 2 + m(y)^ 2 + m(z)^ 2](y)和m单位。因此,激发态之间的能量间隔E必须被转换成au单位来代替上述公式。 2. TDDFT做几何结构优化\\一般认为,由于TDDFT不能计算激发态的波形,所以不能用来研究激发态的平衡态和分子性质。因此,利用Gaussian98提供的TDDFt方法研究激发态,只能得到垂直激发能,不能用来计算平衡键长Re,键角和振动频率ωe。然而,最近的报道说,TDDFT在理论上证明了可行的几何优化。目前已经有消息称高斯计划将TDDFT几何优化添加到新版本的高斯。可以参考:Van Caillie C,Amos RD,使用梯度校正函数的密度泛函理论激发能量的几何导数,Chem。物理学。 Lett。,317(1-2),2000,p.159-164和Hsu,Hirata,Head-Gorden,J.Phys。 A105,2001,451-458。 TDDFT线性多烯烃的激发能和共振强度。 3.自定义赝势(也称为模型电位,有效核电位)高斯伪代码如下:#Becke3LYP / Gen伪= Read Opt Test HF / 6-31G(*)Cr(CO)6 0 1 Cr 0.0 0.0 0.0分子规格继续... CO 0 6-31G(d)**** Cr 0 LANL2DZ **** Cr 0 LANL2DZ这里,Cr原子的基团LanL2DZ进入两次。两个LanL2DZ是不同的。第一个描述价电子,或核电位以外的电子。第二个描述内部电子,即核电位部分。如果您自定义ECP基地,那样输入。当然还有其他的输入法。 \\ u0026这是在CIS计算中定制In原子的基础集合。 In的ECP使用程序内部定义的SDD赝势,Cl原子在路径部分使用LanL2DZ。 #CIS(NStates = 3,50-50)/ LanL2DZ ExtraBasis伪=读取InCI激励状态0,1 In Cl In rr = 2.401 IN 0 S 3 1.00 1.53256200 0.261774000 0.944269000 -0.625406000 0.176733000 0.651304000 S 1 1.00 0.643430000E-01 1.00000000 P 3 1.00 0.641474000 0.211103000 0.342739000 -0.243029000 0.109546000 -0.560191000 P 1 1.00 0.389580000E-01 1.00000000 ****在0 SDD *如果您自定义ECP,则可以将上述SDD更改为:IN 0 IN-ECP 4 46 G POTENTIAL 1 2 1.00000000 0.00000000 SG POTENTIAL 2 2 1.43509100 29.16521900 2 0.69580500 -4.19080600 PG潜力2 2 1.44083200 36.99054200 2 0.70139200 -3.36582000 DG潜力1 2 0.96123600 20.00053100 FG潜力1 2 0.88436900 -6.01909200 \\ u0026 lt; / \\ u0026 gt;价电子群和核势序不能逆转。要下载更多的ECP基础套件,请访问:http://www.emsl.pnl.gov:2080/forms/basisform.html \\ u0026 nbsp;区分核电位和价电子基的方式是:核电位的第一行有一个核电位,而不是内层电子编号,甚至是固定的。像上面的例子一样,核电位取代In的46个内层电子。 \\ u0026如何得到高斯ECP基集的伪势?只需将IOp(3/18 = 1)添加到Route Section行。 IOp的功能非常强大,可以在高斯网页上找到。这只是其中之一。例如,要导出一个Cl原子完整的ECP基本集CEP-121G,您必须使用GFPrint和IOp(3/18 = 1)关键字如下:#HF / CEP-121G GFPrint IOp(3/18 = 1)\\你好吗?其中GFPrint打印CEP-121G的价电子基座部分。当然,如果观察cc-pV [DT] Z,6-31G和其他非ECP基本组,则不需要添加IOp(3/18 = 1),只需使用GFPrint。这是输出中Cl原子的CEP-121G基本部分。第一部分是GFPrint结果,第二部分是IOp(3/18 = 1)的结果:************** *********** *********************************** *********** *原子轨道*高斯函数************ *********** ********************************功能外壳尺寸*编号类型因子*指数S-COEF P-COEF D-COEF F-COEF ************************ *********** ***************************** 1 - 4 SP 1.00 0.222500D + 010.100000D + 010.100000D + 010.000000D + 000.000000D + 00 5-8 SP 1.00 0.117300D + 010.115280D + 000.299520D + 000.000000D + 000.000000D + 00 0.385100D + 000.847170D + 000.583570D + 000.000000D + 000.000000D + 00 9-12 SP 1.00 0.130100D + 000.100000D + 010.100000D + 010.000000D + 000.000000D + 00 ******************************* *********** ************************************************** *********** ======================================= =========== ==================== PSEUDOPOTENTIAL PARAMETERS ================= ========= = =================================== ==============中心原子价电子数量电子数量R指数系数======================= ============== ==================================== === 1 17 7 D及以上1 4.8748300 -3.40738000 S - D 0 17.0036700 6.50966000 2 4.1038000 42.27785000 P - D 0 8.9002900 3.42860000 2 3.5264800 22.15256000 =================== ============== ==================================== ==================以上大部分内容来自www.ccl.net论坛。除了讨论使用GAMESS和其他程序外,其内容大部分是高斯的,少数高斯公司经常参与讨论。建议你有机会看看

关键词: 社会科学