DFT方法和HF方法地比较

实用文案DFT方法和HF方法的比较计算中非常常用的方法。 DFT方法考虑了电子相关，而且算起来也快，所以很多人喜欢用，其中又以 B3LYP最为常用，所以清楚每种方法的优劣和适用范围有利于我们对问题的理解。HF方法忽略的大部分的电子相关。 相反，在很多时候DFT方法常常是过多的考虑了电子相关，这会使得过渡态的能量偏低，造成算出来的活化能偏低而且计算氢键的键能也会偏低。并且在计算有机分子的芳香性也不好， DFT会过多考虑电子离域，导致计算出来的能量偏低但对于过渡金属、有机生物分子， DFT方法都能很好的处理，这是它比其它方法好的地方。由于HF方法忽略了电子相关，所以在处理弱键的时候是不好的，所以我们一般也不用它来计算活化能（偏高）和解离能（偏低）并且在体系中如果有孤对电子、共扼体系的话，电子比较松散或者离域的比较厉害， 这时候电子相关也很重要，HF方法也不能很好的定量解决。B3LYP与其它泛函相比，对分子基态得到的特性一般相差不大。 所以B3LYP很常用也无可厚非。但是一些人用TDDFT算激发能也用B3LYP，而且只用B3LYP，很差的结果也拿出来发表，根本不考虑原因，这就有问题了。上个世纪末，很多使用TDDFT算激发能的文章都得到一个相同的结论， 就是B3LYP作TDDFT激发能计算的结果是不可靠的：对不同的分子体系，有的时候跟实验值相当接近， 有的时候却差得不得了。因此在做 TDDFT激发能计算的时候，应该多试几种泛函，特别是没标准文档实用文案有实验值，或者 B3LYP的结果比较差的时候。如果图省事，想找一个通用的泛函作TDDFT计算，推荐用 PBE极其几种改进版，Gaussian 98的用户可以用B3PW91。它们比较稳定，虽然不总是最好，但也不会太差。另外， HCTH系列泛函（特别是HCTH407）也值得一试。B3LYP之所以计算 TS能量会偏低,主要在于其交换相关势不够准确 ,特别是在长程区的渐近行为不够好 ,也正是如此,b3lyp是不可能准确计算氢键.上次徐昕老师在作报告时就提到了 b3lyp在计算氢键和vanderwal 力结果是poorandbad,他发展的 X3LYP是可以得到较好的结果的 .用B3LYP计算弱相互作用体系He2是不成键的,而X3LYP计算是成键的,所以X3LYP在远程势的渐近行为是比较好的,如果用来计算过渡态的话也比 b3lyp 要准确的多.当然,现在出现了很多的泛函结果都会比 b3lyp准确,不过要加入程序还尚待时日 ,希望大家在计算前要仔细考虑所用的方法合不合适 .标准文档实用文案量子化学计算方法的演进主要分为(1)分子轨道法（简称 MO法）和(2)价键法。以下只介绍 MO法。分子轨道法的核心是哈特里—福克—罗特汉方法，简称 HFR方程，它是以三个在分子轨道法发展过程中做出卓越贡献的人命名的方程。 1928年D·R·哈特里（Hartree）提出了一个将 N个电子体系中的每一个电子都看成是由其余的N-1个电子所提供的平均势场运动的假设。 这样对于体系中的每一个电子都得到了一个单电子方程（表示这个电子运动状态的量子力学方程） ，称为哈特里方程。使用自洽场迭代方式求解这个方程（自洽场分子轨道法） ，就可以得到体系的电子结构和性质。哈特里方程未考虑由于电子的自旋而需要遵守泡利原理。 1930年，B·A·福克（Fock）和J·C·斯莱特（Slater）分别提出了考虑泡利原理的自洽场迭代方程，称为哈特里—福克方程。它的单电子轨函数 （即分子轨道）取为自旋轨函数（即电子的空间函数与自旋函数的乘积）。泡利原理要求，体系的总电子波函数要满足反对称化要求，即对于体系的任何两个粒子的坐标交换都使总电子波函数改变正负号，而斯莱特行列式波函数正是满足反对称化要求的波函数。将哈特里—福克方程用于计算多原子分子，会遇到计算上的困难。 C·C·J·罗特汉（Roothaan）提出将分子轨道向组成分子的原子轨道（ AO）展开，这样的分子轨道称为原子轨道的线性组合（简称 LCAO）。使用LCAO-MO，原来积分微分形式的哈特里—福克方程就变为易于求解的代数方程，称为哈特里—福克—罗特汉方程，简称 HFR方程。RHF方程闭壳层体系是指体系中所有的电子均按自旋相反的方式配对充满某些壳层（壳层指一个分子能级或能量相同的即简并的两个分子能级） 。这种体标准文档实用文案系的特点，是可用单斯莱特行列式表示多电子波函数（分子的状态），描述这种体系的HFR方程称为限制性的HFR方程，所谓限制性，是要求每一对自旋相反的电子具有相同的空间函数。限制性的HFR方程简称RHF方程。UHF方程开壳层体系是指体系中有未成对的电子（即有的壳层未充满）。描述开壳层体系的波函数一般应取斯莱特行列式的线性组合，这样计算方案就将很复杂，然而对于开壳层体系的对应极大多重度（所谓多重度，指一个分子因总自旋角动量的不同而具有几个能量相重的状态）的状态（即自旋角动量最大的状态）来说，可以保持波函数的单斯莱特行列式形式（近似方法）以描述这些体系。从头计算法原则上讲，有了HFR方程（不论是RHF方程或是UHF方程），就可以计算任何多原子体系的电子结构和性质。真正严格的计算称之为从头计算法。在从头计算法里，分子轨道由组成体系的原子的全部原子轨道线性组合而成。对于原子轨道有不同的选法。 斯莱特型轨道适于描写电子云的分布， 但在计算一些积分时包含对无穷级数的积分， 十分麻烦，所以在从头计算法里，常取高斯型函数做为基函数，取一个高斯型函数或数个高斯型函数的线性组合模拟一个原子轨道。已经有了不少进行多原子体系的从头计算法的标准的计算机程序， 如J·A·波普尔研究集体推出计算机程序系列：高斯系列。该系列第一个公开的版本为高斯70（指1970年），以后差不多每年更新一次，功能逐次增多，算法也日趋完善。赝势价轨道从头计算法 从直观的化学观出发，可以想见，在原子形成分子时，仅仅原子的价层电子发生了较大的变形， 而内层电子分布则改变较小。 为了节省计算时间，而又不失去计算精度，从 20世纪70年代以来，开始只考虑原标准文档实用文案子的价电子，而把内层电子和原子核看成一个凝固的原子实， 用一个核模型势来代替内层电子与价电子的从头计算法， 称为赝势价轨道从头计算法。不同的方案对于赝势的取法不尽相同，且大都能得到与全电子从头计算法相近的结果， 因而大大节省了计算时间，特别对于含有重原子的体系，恰好是全电子从头算难以处理的。赝势价轨道从头计算法将会在过渡金属络合催化的量子化学研究方面发挥重要作用。组态相互作用法 比从头计算法精确度更高的组态相互作用法有两个优点：一是不依赖于试探波函数的形式，就能原则上提非相对论薛丁谔方程的精确解；二是原则上可用于原子或分子体系的任何稳定态。近似计算法 在半经验计算法中，目前最常用的近似是零微分重叠（简称ZDO）；近似程度最高的是全略微分重叠（简称 CNDO）近似，是J·A·波普尔在1965年提出的。此外，在近似计算法中，还有休克尔分子轨道法和推广的分子轨道法。标准文档实用文案Chem3D 使用心得这是一个三维分子结构演示软件，上面有和 Gaussian的链接，可以直接从Chem3D打开Gaussian运算，本身也带有分子力场计算功能，还可以和Mopac链接可以打开的文件主要有以下几种：*.gjf,*.fch,*.cub, 这几种是Gaussian的输入、formcheck 和cube文件*.dat,*.mop,*.mpc,mopac 的输入文件*.zmt,这个是内坐标格式的分子构型文件，和 HyperChem 、mopac通用*.pdb,*.ent,Protein DB 格式，和 HyperChem 、AlChem 等等通用，Gaussian本身也可以用这两种格式来生成内坐标文件*.mol,MDL 的mol 格式，和HyperChem 稍微有点不兼容，但是如果用AlChem 打开，在保存一下（仍然是 *.mol），Chem3D 和HyperChem 都可以打开了。最有用的可能是 *.zmt,*.mol,*.gjf,*.fch 四种格式 ，比如用 Gaussian优化以后，在Ultilities 选项里面选formcheck ，就可以把check文件*.chk转化成*.fch，而*.fch是可以直接用Chem3D 打开的。*.zmt是一内坐标的分子结构文件，我这里有个小软件可以把这种格式的文件转换成 *.gjf，是Gaussian的分子内坐标文件，但没有routesection 等等。Chem3D 除了生成Gaussian输入文件以外，还有许多不错的功能，比如可以手动地转动某一个键， 可以用鼠标显示某个键长，键角和二面角等地大小，等等，还可以显示所以内坐标。显示也有各种模型方式，比如球棍、 stick 等等，可以显示原子编号原子名称，但是没有GaussView和HyperChem 以及AlChem，WinMopac 等等漂亮。标准文档实用文案Gaussian 常见问题分析1．检查是否有初始文件错误在命令行中加入 %kJobL301or%kJobL302如果通过则一般初始文件 ok。常见初级错误：自旋多重度错误变量赋值为整数变量没有赋值或多重赋值键角小于等于0度，大于等于180度分子描述后面没有空行二面角判断错误，造成两个原子距离过近分子描述一行内两次参考同一原子,或参考原子共线2．SCF(自洽场)不收敛则一般是 L502错误省却情况做64个cycle迭代(G03缺省128cycles)修改坐标，使之合理改变初始猜Guess=Huckel或其他的，看Guess关键词。增加叠代次数SCFCYC=N（对小分子作计算时最好不要增加，很可能结构不合理）iop（5/13=1）这样忽略不收敛，继续往下做。3．分子对称性改变标准文档实用文案修改坐标，强制高对称性或放松对称性b.给出精确的、对称性确定的角度和二面角。 如CH4的角度给到109.47122c.放松对称性判据 Symm=loosed.不做对称性检查 iop(2/16=1) （最好加这个选项） iop(2/16=2) 则保持新的对称性来计算4．Opt时收敛的问题修改坐标，使之合理增加叠代次数optcyc=N5．优化过渡态，若势能面太平缓，则不好找到。iop(1/8=10) 默认30（下一个结构和该结构的差别 0.3Å ），可改成10。如果每一步都要用到小的步长，应该加 opt(notrustupdate)6．在CI（组态）方法中如QCISD(T),CCSD(T),CID方法中，省却最大循环50，若出错（L913错误）解决方法：#PQCISD(maxcyc=N) 注：N≤5127．优化过渡态opt=TS(给出过渡态)opt=qst2( 给出反应物和产物)opt=qst3( 给出反应物和产物和过渡态 )标准文档实用文案a. 用G03时的出错 opt=ts 必须加FC(forceconstant) 写法：opt=(TS,calcFc)or opt=(TS,calchffc) 计算HF力常数，对QCISD，CCSD等方法用；or opt=(TS,modRedundant) (最好写这个)如果计算采用QCISD计算（不好计算FC）则写为QCISDopt=(TS,calcHFFC)( 用HF计算FC)8.无法写大的Scratch文件RWFa.劈裂RWF文件%rwf=loc1,size1,loc2,size2, ⋯⋯..,locN,-1b.改变计算方法 MP2=Direct 可以少占硬盘空间c.限制最大硬盘 maxdisk=N GB，****MB,有些系统写 2GB会出错，可以写2000MB9.FOPT出错原因是变量数与分子自由度数不相等。 可用POPT或直接用OPT优化过渡态只能做一个STEP原因是负本征数目不对添加iop（1/11）=1或者noeigentest标准文档实用文案量化小节（转）关于自旋多重度定义:多重度＝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啊什么的），可能会高多重度的能量低，所以需要都算算，看哪个能量更低。关于赝势：标准文档实用文案简单来说，赝势就是不计算内层电子，而是把内层电子的贡献用一个势来描述，放在哈密顿里面。适用于重元素。赝势基组，实际上包括 赝势和基组两个部分，内层电子采用赝势，即 effectivecorepotential(ECP) ，外层价电子采用一般的基组。 比如：LanL2DZ:D95Vonfirstrow,LosAlamosECPplusDZonNa-Bi.就是对第一行原子是 D95V(这个是非赝势基组)，对Na-Bi是使用一个叫做 LosAlamos的有效核势加上一个 DZ基组。所以Lanl2dz就是对前面的原子全电子基组，对后面的原子是赝势基组。 （再次说明，量化里面， C，O那一行，算周期表的第一行）使用赝势的3个原因：1。没有相应的全电子基组。2。减少计算量。3。赝势可以包含重金属相对论效应的修正。在高斯中，lanl2dz 基组，在手册中可以查到其定义为：LanL2DZ:D95Vonfirstrow,LosAlamosECPplusDZonNa-Bi ，也就是说，对于C，O等元素来说（量化中，H大概算第0周期，C，O才是第一周期），Lanl2dz实际上还是全电子基组，而对于Na以后才是对内层电子用 LosAlamosECP赝势，外层电子用 DZ基组。标准文档实用文案使用赝势的输入文件：1.所有原子使用 lanl2dz--------------------------#HF/lanl2dzoptlanl2dzforallatoms0 20.00.00.0C0.00.01.2Cu0.00.03.2--------------------------2.所有原子使用 lanl2dz的另一种输入方法。--------------------------#HF/genecpoptlanl2dzforallatoms标准文档实用文案0 20.00.00.0C0.00.01.2Cu0.00.03.2CO0lanl2dz****Cu0lanl2dz //定义价电子的基组,CO0是碳，氧，零，其中零用作终止符号。****Cu0lanl2dz //定义内层电子的赝势--------------------------3.混和基组，即有的使用全电子，有的使用 Lanl2dz。格式同2。--------------------------#HF/genecpoptlanl2dzforCu,6-31G(d)forCandO标准文档实用文案0 20.00.00.0C0.00.01.2Cu0.00.03.2CO06-31G(d) //另一种全电子基组****Cu0lanl2dz //定义价电子的基组****Cu0lanl2dz //定义内层电子的赝势--------------------------开壳层和闭壳层闭壳层计算就是对于多重度是 1的体系，此时α和β的电子数目相同，可以把α和β配对，成对的α和β使用同一个轨道，一个轨道上填充 2个电子。开壳层计算就是对α和β电子分别计算，一个轨道上只填充 1个电子，一般来说，标准文档实用文案多重度是1时，开壳层计算和闭壳层计算会给出相同的结果。限制开壳层计算是对多重度大于 1的体系，此时α和β的电子数目不同，设有 m个α和n个β电子，m>n，那么让前n个轨道上每个填充一个α和一个β，剩下的m-n个α电子再填充m-n个轨道。即前n个轨道是闭的（每个轨道 2个电子），后m-n个轨道是开的（每个轨道 1个电子）在高斯中，以HF为例，闭、开、限制开壳层计算分别是 RHF，UHF，ROHF。如果只写HF，则按下面的方式取默认方法：对多重度是1的体系，默认为 RHF，对多重度大于 1的体系，默认是UHF。关于收敛问题 (L502,L508,L9999)对于一个优化计算，它的过程是先做一个 SCF计算，得到这个构型下的能量，然后优化构型，再做 SCF，然后再优化构型。。。因此，会有两种不收敛的情况：一是在某一步的 SCF不收敛(L502 错误)，或者构型优化没有找到最后结果L9999错误）。预备知识:计算时保存chk文件,可以在后续计算中使用 guess=read 读初始猜测.对于SCF不收敛，通常有以下的解决方法：使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波标准文档实用文案函数稳定。如果这个还不收敛，会提示 L508错误。改变键长，一般是缩小一点，有时会有用。计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。待补充.对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态）这有几种可能性：看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1).(flyingheart)标准文档实用文案一个建议是，对于大体系，难收敛体系，先用小基组，低精度算法优化一下，以得到较好的初始构型，再用高精度的计算接着算。如果前面的方法保留了 chk文件，重新计算时需要使用 geom=allcheck 读入构型（就不必麻烦地写构型了），guess=read （读入初始波函数，可以加快第一步 SCF收敛）。关于对称性：（原贴 ）刚刚在zixia上看到zhuoliu 的大作。结合自己以前知道的东西，总结一下。1。gaussian 中输入什么对称性，一般优化的结果仍然还是那个对称性，比如CO2，如果初始两个CO键长输入不是完全相等（比如一个1.214,一个1.215），那么程序就会判断为 C∞v对称，那么优化结果虽然键长几乎相等，但仍然认为是C∞v，这个从振动频率或者分子轨道对称性上可以看出来。 --我们知道，CO2实际上是直线的两边对称的构型，其对称性应该是 D∞v。因此，为了得到高的对称性，必须输入的时候，精确地输入数值，比如 sqrt(2)，就要保留很多的小数点，180.0角，就不能写成179.9。2。有时计算过程中对称性会变化，比如做过渡态的时候，这时需要用IOP(2/16=3) ，否则计算会出错退出。3。比如用直角坐标输入一个正三角形构型，其对称性应该是 D3h，但是如果输入的小数点后面的数字不够多，那么常常得到的是 C2v或其它。为了消除输入标准文档实用文案文件中坐标的有效位数的影响， 得到较高的对称性，可以降低对称性判断的严格性。一般可以用 symm=loose ，这等价于 IOP(2/17=4, 2/18=3)。还可以减小这4和3这两个数值，使得更加loose，但不能过小，否则会出错。symm=loose只是在第一步判断输入构型的对称性时用到。此外，也可以用 gaussview 来调整设置初始构型的对称性。4。如果要降低对称性，那么可以用 symm(PG=C3v) 等等来做。使判断出来的对称性为C3v的一个子群。即由 PG来限制最高对称性。附用到的IOP的详细解释。IOp(2/16)action taken if the point group changes during an optimization.IOP(2/16=0)Abort the job. IOP(2/16=1)Keep going. IOP(2/16=2)Keepgoingandleavesymmetryon,usingtheoldsymmetry.IOP(2/16=3)Keepgoingandleavesymmetryon,usingthenewsymmetryIOp(17)Tolerancefordistancecomparisonsinsymmetrydetermination.Default(determinedinthesymmetrypackage,currently1.d-8).N>010-N.N<010N,usethesametolerancefororientation.标准文档实用文案IOp(18)Tolerancefornon-distancecomparisonsinsymmetrydetermination.Default(determinedinthesymmetrypackage,currently1.d-7).N>010-N.N<0 10N,usethesametolerancefororientation.supi）附加一个我的：有时候分子对称性太高经常无法优化收敛，这时可以加一个nosymm，有时挺管用的。关于频率计算 关键字Freq频率计算，需要求能量对坐标的二阶导数，得到力常数，然后除以原子质量，求得振动频率。得到的振动频率，对应于实验上的红外和拉曼光谱（根据对称性，判断是否有红外或拉曼活性）。默认计算红外振动的强度，不计算拉曼强度。频率计算必须在优化好的构型下进行。因此，必须使用和优化一样的方法和基组，进行频率计算。有些方法，使用解析方法求解二次导数，例如 DFT，HF等；而有些方法只能用数值方法，例如 CCSD等。数值方法计算量很大，但是如果中断了，可以用标准文档实用文案freq(restart) 接着计算，而解析方法则不能 restart。频率修正因子，是一个经验的东西，别人也是计算了一系列分子以后，总结（或者说拟合）出来的。所以有多准，是不是适合你在算的体系，就不敢保证了。一些修正因子： （从前一个地方下载回来，便于国内访问）关于虚频（原贴： ）首先，什么是频率。中学的时候我们学过简谐振动，对应的回复力是 f=-kx,对应的能量曲线，是一个开口向上的二次函数 E=kx^2/2. 这样的振动，对应的x=0的点是能量极小值点（简单情况下也就是最小值点）。这时的振动频率我们也会求：ω=2πsqrt(k/m) 。显然它是一个正的频率，也就是通常意义下的振动频率。那么，一维情况下，如果能量曲线是一个开口向下的二次曲线呢？首先， 从能量上看，这是个不稳定的点，中学的物理书上称为“不稳平衡” 。用现在的观点看，就是这一点导数是零（受力为 0），且是能量极大值。如果套用上面的公式，“回复力”f=-k'x（实际上已经不是回复，而是让 x越来越远了），这里k'是个负数，ω=2πsqrt(k'/m) 显然就是一个虚数了，即所谓的虚频。 Gaussian里面给出一标准文档实用文案个负的频率，就是对应这个虚频的。实际情况下，分子的能量是一个高维的势能面， 构型优化的时候，有时得到了极小值点，这样这个点的任意方向上， 都可以近似为开口向上的二次函数， 这样这里对应的振动频率就都是正的。 对于极大值点，在每个方向都是开口向下的二次函数，那么频率就会都是负的——当然一般优化很少会遇到这样的情况。对于频率有正有负的情况，说明找到的点在某些方向上是极大值， 有些方向上是极小值。如果要得到稳定的能量最低构型， 显然需要通过微调分子的构型， 消去所有的虚频。如何微调？要看虚频的振动方向。 想象着虚频对应的就是开口向下的二次函数，显然，把分子坐标按照振动的方向移动一点点， 分子应该就可以顺着势能面找到新的稳定点，但是也不能太小。而所谓的过渡态，则是连接反应物和产物之间的最低能量路径上的能量极大值。好比山谷中的 A，B两点，它们之间的一个小土丘，就是过渡态，从 A到B的反应，需要越过的是这个小土丘，而不是两边的高山。这样，过渡态就是在一个方向上是极大值， 而在其它方向上都是极小值的点。因此，过渡态只有一个虚频。优化得到虚频，消虚频的方法，就是根据虚频对应的振动位移 (输出文件每一个频率下面都有一个类似坐标的 3列数值，给出的就是每个原子在这个振动频率时候的移动方式)，把这些位移，乘以一个适当大小的因子，直接和原来平衡构型相加，就得到新的构型了。举例说明：标准文档实用文案#HF/3-21Goptfreqtest02H0.00.00.0H0.00.01.0H0.00.0-1.0上面的计算，是优化一个 H-H-H 的直线构型然后算频率。我们知道，3个H的稳定构型，应该是一个 H2分子，加上一个H原子。但是由于我们输入的时候，中间H原子两边的H的键长相等，因此，在这个对称性的限制下，结果给出的“稳定构型”是：H0.0000000.0000000.000000H0.0000000.0000000.934091H0.0000000.000000-0.934091频率分析，有一个虚频：1 23SGU PIU标准文档实用文案PIUFrequencies---2291.19091121.14861121.1486Red.masses--1.00781.00781.0078Frcconsts--3.11720.74640.7464IRInten--42.97068.81678.8167RamanActiv--0.00000.00000.0000Depolar(P)--0.00000.00000.0000Depolar(U)--0.00000.00000.0000AtomANXYZXYZXYZ110.000.000.820.820.030.00-0.03 0.82 0.002 1 0.00 0.00 -0.41 -0.41 -0.02 0.000.02 -0.41 0.003 1 0.00 0.00 -0.41 -0.41 -0.02 0.00标准文档实用文案0.02 -0.41 0.00第一个是虚频。可以看出，其振动模式是：中间的 H的z坐标变大，两边的两个H的z坐标变小，（然后是中间的H的z坐标变小，两边的H的z坐标变大，完成一次振动），这个如果还不清楚，画一下图就知道了。现在，我们就把原来的坐标，加上振动模式对应的坐标：H0.0000000.0000000.000000H0.0000000.0000000.934091＋H0.0000000.000000-0.9340910.00 0.00 0.820.00 0.00 -0.41 0.00 0.00 -0.41 直接加的结果，是H0.0000000.0000000.820000H0.0000000.0000000.524091H0.0000000.000000-1.344091画图知道，显然对原来的构型，变化太多了。对这个例子，如果取一个系数 0.1乘以振动坐标，再求和，那么结果就是：标准文档实用文案H0.0000000.0000000.082000H0.0000000.0000000.893091H0.0000000.000000-0.975091显然就合理多了。用这个坐标去重新优化，就可以得到真正的没有虚频的稳定结构了。这个系数0.1，只是个经验的数值，对于不同体系，可以自己设定。设置过小，会得到同样的虚频，设置过大，可能会得到别的构型（当然也可能碰巧得到更好的构型）。关于B