materialstudio使用经验总结

1、关于K点1 .应当使用多少个k网格？很难一般地回答，只能给出一般建议。注意：一定要检查k网格，首先用较粗糙的网格计算，接下来用精细的网格计算。通过比较两次的结果，决定选用较粗糙的网格，或是继续进行更精细网格的计算，直到达到收敛。金属体系需要精细的网格，绝缘体使用很少的k点通常就可以。小单胞需要精细格点，大单胞很可能不需要。因此：单位晶胞内原子数很多（比如40-60个）的绝缘体，可能仅需要一个（移动后的）k点。另一方面，面心立方的铝可能需要上万个k点以获得好的DOS对于孤立原子或分子的超晶胞，仅需要在Gamm疝计算。对于表面（层面）的超晶胞计算，仅需要（垂直于表面）z方向上有1个k点。甚至可以增

2、加晶格参数c,这样即使对精细格点，沿z方向上也只产生一个k点（产生k点后，不要忘记再把c改回）。2 .当体系没有出现时间反演对称操作时，是否加入？大多数情况下的回答是“是”，只有包含自旋-轨道耦合的自旋极化（磁性）计算除外。这时，时间反演对称性被破坏（+k和-k的本征值可能不同），因此决不能加入时间反演对称性。3 .是否移动k网格？（只对某些格子类型有效）“移动”k网格意味着把所有产生的k点增加（x,x,x）,把那些位于高对称点（或线）上白kk点移动到权重更大的一般点上。通过这种方法（也即众所周知的“特殊k点方法”）可以产生等密度的，k点较少的网格。通常建议移动。只有一点注意：当对半导体的带隙

3、感兴趣时（通常位于GammaX,或BZ边界上的其它点），使用移动的网格将不会得到这些高对称性的点，因此得到的带隙和预期结果相比或大或小。这个问题的解决：用移动的网格做SCF循环，但对DOS计算，改用精细的未移动网格。关于k空间布点的问题，建议参阅以下文献Phys.Rev.B49,16223（1994）如何构建缺陷晶体结构晶体结构改成P1,然后去掉想抹去的原子就可以了在ms中如何做空穴对于金属缺陷，是直接剪切一个原子？个人经验：就是直接把原子去掉就OK;如果不是正版软件，有可能出现同时去掉其他同位置的原子，如果这种情况，就重新定义，问题就不会出现了.还有，一般考虑孔穴的时候，都要标明哪些原子的迟

4、豫，具体为什么不知道，国外的文献有提到.希望有做空位的一起多讨论.我Q:183876402PDOS选项计算DOS时，选择PDOS,可以画出s,p,d轨道的DOS,但无法画出某一个原子的s,p,d图关于PDOS的Chart中求积分的问题在用Castep计算出PDOS后,如何在Chart中对曲线局部进行积分？将Chart输出为cav格式,然后在excel中求和？简单，把数据导出，在Origin里作图，程序里有积分微分卷积功能，在数据分析下面。作图时选取积分范围。优化结构算能带一般需要优化结构。如果选择实验的参数，全部固定的话就不需要了如何做二维电子密度图MS结果文件夹中*.grd文件内存储的是三维

5、空间各点的电荷密度值，利用这个数据就可以得到二维的电荷密度等值线图，应该有专门的软件能画，不过用matlab编自己编程序也不难，其中关键命令是contourslice,实现在某一平面内绘制等高线对DOS图的分析根据DOS的积分曲线可以计算出，对于表面Nb和C原子，大约有16.8%和14.8%的电子态集中在4.02.0eV的区域，而对体相原子则分别为6.4%和6.0%表面吸附我做H在ZnO上吸附。刚开始时后我构建的吸附构型忘记imposesymmetry了，Groupname是P1。在第二次计算的时候我加上了symmetry。两次计算差别出来了：(1)首先是imposesymmetry后，sup

6、ercell中的原子位置由原来的现面跑到了上面，也就是和真空层换了一下位置！而且吸附原子竟然超出了supercell的白色线框！(2)比较两次计算的DOS,一模一样；但是两次计算的BandStructure却有很大差异！回答：1)只是显示问题(2)BandStructure有很大差异是指那种差异？使用不同对称性计算能带时，默认计算的K点是不同的，所以图像肯定不同。如果你确定是选择计算了同样的K点，能带结构仍然不同，那可能是采用对称后结构变化导致的。如何计算结构中某一元素的分波态密度1 .我最近看了关于氧气锌的论文，上面有锌原子的分波态密度图，可是我怎么都没算出来，一直得到氧和锌原子的分波态密度

7、，希望哪位大侠指教2 .在计算性质选择了计算densityofstate是，对话框下面有一个口calculatePDOS,把它勾选上，计算成功后，在分析就能看到总态密度和分态密度，随你选择。3 .按照楼上的指导计算成功后，需要哪个原子的PDOS,就选中哪个原子，然后在analysis里DOS项前打勾，partial项打勾，再view就ok了。有关能带分析能带图分析能带图的横坐标是在模型对称性基础上取的K点。为什么要取K点呢？因为晶体的周期性使得薛定谓方程的解也具有了周期性。按照对称性取K点，可以保证以最小的计算量获得最全的能量特征解。能带图横坐标是K点，其实就是倒格空间中的几何点。其中最重要也

8、最简单的就是gamma那个点，因为这个点在任何几何结构中都具有对称性，所以在castep里，有个最简单的K点选择，就是那个gamma选项。纵坐标是能量。那么能带图应该就是表示了研究体系中，各个具有对称性位置的点的能量。我们所得到的体系总能量，应该就是整个体系各个点能量的加和。记得氢原子的能量线吧？能带图中的能量带就像是氢原子中的每条能量线都拉宽为一个带。通过能带图，能把价带和导带看出来。在castep里，分析能带结构的时候给定scissors这个选项某个值，就可以加大价带和导带之间的空隙，把绝缘体的价带和导带清楚地区分出来。DOS叫态密度，也就是体系各个状态的密度，各个能量状态的密度。从DOS

9、图也可以清晰地看出带隙、价带、导带的位置。要理解DOS,需要将能带图和DOS结合起来。分析的时候，如果选择了full,就会把体系的总态密度显示出来，如果选择了PDOS,就可以分别把体系的s、p、d、f状态的态密度分别显示出来。还有一点要注意的是，如果在分析的时候你选择了单个原子，那么显示出来的就是这个原子的态密度。否则显示的就是整个体系原子的态密度。要把周期性结构能量由于微扰裂分成各个能带这个概念印在脑袋里。最后还有一点，这里所有的能带图和DOS的讨论都是针对体系中的所有电子展开的。研究的是体系中所有电子的能量状态。根据量子力学假设，由于原子核的质量远远大于电子，因此奥本海默假设原子核是静止不

10、动的，电子围绕原子核以某一概率在某个时刻出现。我们经常提到的总能量，就是体系电子的总能量。这些是我看书的体会，不一定准确，大家多多批评啊！Energy(eV)10icissors工源1口£凰鹿冲-3001234567891011121314DensitycfStates(efertrontj/eV)七改CflSTEPBandSlrurtursCASTEPDensityof9At：e$如何分析第一原理的计算结果转自：转自量化网摘要：本文总结了对于第一原理计算工作的结果分析的三个重要方面，以及各自的若干要点用第一原理计算软件开展的工作，分析结果主要是从以下三个方面进行定性/定量的讨论：1

11、、电荷密度图(chargedensity);2、能带结构(EnergyBandStructure);3、态密度(DensityofStates,简称DOS)。电荷密度图是以图的形式出现在文章中，非常直观，因此对于一般的入门级研究人员来讲不会有任何的疑问。唯一需要注意的就是这种分析的种种衍生形式，比如差分电荷密图(def-ormationchargedensity)和二次差分图(differencechargedensity)等等，力口自旋极化的工作还可能有自旋极化电荷密度图(spin-polarizedchargedensity)。所谓差分”是指原子组成体系(团簇)之后电荷的重新分布，七次”是

12、指同一个体系化学成分或者几何构型改变之后电荷的重新分布，因此通过这种差分图可以很直观地看出体系中个原子的成键情况。通过电荷聚集(accumulation)/损失(depletion)的具体空间分布，看成键的极性强弱；通过某格点附近的电荷分布形状判断成键的轨道(这个主要是对d轨道的分析，对于s或者p轨道的形状分析我还没有见过)。分析总电荷密度图的方法类似，不过相对而言，这种图所携带的信息量较小。能带结构分析现在在各个领域的第一原理计算工作中用得非常普遍了。但是因为能带这个概念本身的抽象性，对于能带的分析是让初学者最感头痛的地方。关于能带理论本身，我在这篇文章中不想涉及，这里只考虑已得到的能带，如

13、何能从里面看出有用的信息。首先当然可以看出这个体系是金属、半导体还是绝缘体。判断的标准是看费米能级和导带(也即在高对称点附近近似成开口向上的抛物线形状的能带)是否相交，若相交，则为金属，否则为半导体或者绝缘体。对于本征半导体，还可以看出是直接能隙还是间接能隙：如果导带的最低点和价带的最高点在同一个k点处，则为直接能隙，否则为间接能隙。在具体工作中，情况要复杂得多，而且各种领域中感兴趣的方面彼此相差很大，分析不可能像上述分析一样直观和普适。不过仍然可以总结出一些经验性的规律来。主要有以下几点：1)因为目前的计算大多采用超单胞(supercell)的形式，在一个单胞里有几十个原子以及上百个电子，所

14、以得到的能带图往往在远低于费米能级处非常平坦，也非常密集。原则上讲，这个区域的能带并不具备多大的解说/阅读价值。因此，不要被这种现象吓住，一般的工作中，我们主要关心的还是费米能级附近的能带形状。2)能带的宽窄在能带的分析中占据很重要的位置。能带越宽，也即在能带图中的起伏越大，说明处于这个带中的电子有效质量越小、非局域(non-local)的程度越大、组成这条能带的原子轨道扩展性越强。如果形状近似于抛物线形状，一般而言会被冠以类sp带(sp-likeband)之名。反之，一条比较窄的能带表明对应于这条能带的本征态主要是由局域于某个格点的原子轨道组成，这条带上的电子局域性非常强，有效质量相对较大。

15、3)如果体系为掺杂的非本征半导体，注意与本征半导体的能带结构图进行对比，一般而言在能隙处会出现一条新的、比较窄的能带。这就是通常所谓的杂质态(dopingstate),或者按照掺杂半导体的类型称为受主态或者施主态。4)关于自旋极化的能带，一般是画出两幅图：majorityspin和minorityspin。经典的说，分别代表自旋向上和自旋向下的轨道所组成的能带结构。注意它们在费米能级处的差异。如果费米能级与majorityspin的能带图相交而处于minorityspin的能隙中，则此体系具有明显的自旋极化现象，而该体系也可称之为半金属(halfmetal)。因为majorityspin与费米

16、能级相交的能带主要由杂质原子轨道组成，所以也可以此为出发点讨论杂质的磁性特征。5)做界面问题时，衬底材料的能带图显得非常重要，各高对称点之间有可能出现不同的情况。具体地说，在某两点之间，费米能级与能带相交；而在另外的k的区间上，费米能级正好处在导带和价带之间。这样，衬底材料就呈现出各项异性：对于前者，呈现金属性，而对于后者，呈现绝缘性。因此，有的工作是通过某种材料的能带图而选择不同的面作为生长面。具体的分析应该结合试验结果给出。(如果我没记错的话，物理所薛其坤研究员曾经分析过$beta$-Fe的(100)和(111)面对应的能带。有兴趣的读者可进一步查阅资料。)原则上讲，态密度可以作为能带结构

17、的一个可视化结果。很多分析和能带的分析结果可以一一对应，很多术语也和能带分析相通。但是因为它更直观，因此在结果讨论中用得比能带分析更广泛一些。简要总结分析要点如下：1)在整个能量区间之内分布较为平均、没有局域尖峰的DOS,对应的是类sp带，表明电子的非局域化性质很强。相反，对于一般的过渡金属而言，d轨道的DOS一般是一个很大的尖峰，说明d电子相对比较局域，相应的能带也比较窄。2)从DOS图也可分析能隙特性：若费米能级处于DOS值为零的区间中，说明该体系是半导体或绝缘体；若有分波DOS跨过费米能级，则该体系是金属。此外，可以画出分波(PDOS)和局域(LDOS)两种态密度，更加细致的研究在各点处

18、的分波成键情况。3)从DOS图中还可引入魔能隙"(pseudogap)的概念。也即在费米能级两侧分别有两个尖峰。而两个尖峰之间的DOS并不为零。震能隙直接反映了该体系成键的共价性的强弱：越宽，说明共价性越强。如果分析的是局域态密度（LDOS）,那么震能隙反映的则是相邻两个原子成键的强弱：震能隙越宽，说明两个原子成键越强。上述分析的理论基础可从紧束缚理论出发得到解释：实际上，可以认为震能隙的宽度直接和Hamiltonian矩阵的非对角元相关，彼此间成单调递增的函数关系。4）对于自旋极化的体系，与能带分析类似，也应该将majorityspin和minorityspin分别画出，若费米能级

19、与majority的DOS相交而处于minority的DOS的能隙之中，可以说明该体系的自旋极化。5）考虑LDOS,如果相邻原子的LDOS在同一个能量上同时出现了尖峰，则我们将其称之为杂化峰（hybridizedpeak）,这个概念直观地向我们展示了相邻原子之间的作用强弱。以上是本人基于文献调研所总结的一些关于第一原理工作的结果分析要点。期冀能对刚进入这个领域内的科研工作者有所启发。受本人的水平所限，文章的内容可能会有理论上的不足甚至错误之处，希望大家指出，共同发展第一原理计算物理的方法和研究内容。希望大家多读几遍。http:/cmt.dur.ac.uk/sjc/Castep_Publicat

20、ions/pubblications.html看看人家用castep怎样选择鹰势和K点的关于出错信息的解读（对自己解决问题有帮助）GeomrtryOptimization时，到77步（共100步）就没有反应，我等了一天终于把它停止了，文件中有下列出错信息，请大侠帮忙解读一下。精度是Fine,之前在coarse下可以24步就完成。1、Warning:Therearenoemptybandsforatleastonekpointandspin;thismayslowtheconvergenceand/orleadtoaninaccurategroundstate.Ifthiswarningpers

21、ists,youshouldconsiderincreasingnextra_bandsand/orreducingsmearing_widthintheparamfile.Recommendusingnextra_bandsof6to12.2、BFGS:Warning-RepeatedconsecutiveresetofinverseHessianBFGS:withoutsatisfyingconvergencecriteriawhichBFGS:lookslikeBFGShasrunoutofsearchdirections.BFGS:Warning-Letstryallowingsome

22、uphillstepsandseeifBFGS:wecangetaroundthisbarrier.BFGS:Warning-ItispossiblethatthesystemmaynowconvergetoBFGS:astationarypointOTHERthanthedesiredminimum.BFGS:Hint-thismaybeanindicationthateither:BFGS:a）youareusingapoorguessatgeom_frequency_estBFGS:and/orgeom_modulus_est,orBFGS:b)youareusingunrealisti

23、cconvergencecriteria.BFGS:Suggestthereforethatyouconsiderchangingthem!3、energiesnotcorrectedforfinitebasisset什么意思？要修改那些参数？回答：你可以在property中的SCF增力口emptyband和减小smearing来试试关于结合能形成能的问题我不太明白形成能结合能和总能之间的关系？是否可以指教回答：举例来说吧。对于Si块体材料，通过CASTEP计算得到的能量就叫总能。如果你想把它切开，得到(100)面，通过添加真空层构建一个spercell,计算得到的也叫总能。刚才对于Si块体材

24、料的总能乘上系数再减去Si(100)面总能这就是Si(100)面的形成能。乘上的系数保证他们的原子总数一致。结合能：如果在Si(100)面上吸附了一个Mn原子(假如你关心磁性)，和原来一样计算出总能后，减去Si(100)面的总能再减去Mn原子的总能，得到的就是bindingenergy那位高手能够介绍一下MS中CASTEP的势能模型的适用范围？比如说LDA的(PWC,VWN),GGA的(PW91,BP,PBE,BLYPBOP,VWN-BP,PRBE,HCTH),能够简述一下其大致差别更好LDA局域密度近似(LDA):局域密度近似(LDA)是第一阶梯。它仅仅采用空间点r处的电子密度n(r)来决定

25、那点交换-相关能密度的形式。交换-相关能密度由密度相同的均匀电子气完全确定。泛函的交换部分就准确的用均匀电子气的微分表达。各种不同的局域密度近似(LDA)仅仅是相关部分表示方法不同，所有现代应用的局域密度泛函都基于Ceperly和Alder's在80年代对均匀电子气总能量的MonteCarlo模拟。广义梯度近似(GGA):GGA是Jacob阶梯的第二个台阶，将电子密度的梯度也作为一个独立的变量(|?n(r)|),在描述交换相关能方面，梯度引入了非定域性。GGA泛函包含了两个主要的方向：一个称为无参数”，泛函中新的参数通过已知形式中参数或在其它准确理论帮助下得到。另外一个就是经验方法，未

26、知参数来自于对实验数据的拟和或通过对原Emzerhof(PBE)以及Perdew-WangfromGGA,比如Becke,Lee,ParrandLaplace算符，因此严格上讲属于Jacob子和分子性质准确的计算。Perdew,Burkeand1991(PW91)是无参数的，在量子化学中广泛采用的Yang(BLYP)是经验性。LYP校正采用了密度的二阶阶梯的第三阶，但通常仍然归类为GGA.(见后面英文介绍)请问：在用CASTEP计算材料空位缺陷时，如何才能实现只取一个空位呢* 在build下的symtrey有makeP1* 不是降低对称性，是构建超胞makesuppercell然后就可以删去一

27、个原子了正常的话如果晶体是原胞，直接删除一个原子可能造成很大影响可以建大一点的超晶胞（比如2X2X2超晶胞）然后删掉一个原子，这样比较合理，毕竟缺陷浓度是很低的。* 都可以啊！先看一下showsymmetry是P1就可以直接删一个，如果不是就换成P1,一般makesupercell后自动变成P1了，你可以根据自己所需的浓度选择其中一种方式MS无法收敛请大家帮忙看看了，掺杂Eu的体系，一开始提示增加emptyband,减小smear,都做了，把精度设置到coarse,可最后还是在100轮后提示无法收敛.<-SCF* Warning*max.SCFcyclesperformedbutsyst

28、emhasnotreachedthegroundstate.Currenttotalenergy,E=-10141.84970726eVCurrentfreeenergy(E-TS)=-10141.84970726eV(energiesnotcorrectedforfinitebasisset)NBest.0Kenergy(E-0.5TS)=-10141.84970726eV*Warning:electronicminimisationdidnotconvergeduringfinitebasissetcorrection.*Checkpointfilecannotbewritten.Erro

29、rcalculate_finite_basis:Convergencefailedwhendoingfinitebasissetcorrection.Errorcalculate_finite_basis:Convergencefailedwhendoingfinitebasissetcorrection.*Warning:electronicminimisationdidnotconvergeduringfinitebasissetcorrection.*Checkpointfilecannotbewritten.1MPIAbortbyuserAbortingprogram!1Abortin

30、gprogram!0MPIAbortbyuserAbortingprogram!0Abortingprogram!回答：1.增加smearing容易收敛，但是太大可能结果会有些问题。2 .增加循环的圈数，3 .增加空带数目，4 .体系有自旋单电子的，要设置自旋极化，如果设置非自旋极化，可能不收敛5 .如果是金属体系，比较难收敛，可以试试修改k点和emptyband试试看。1 .几何优化可以在多大范围内改变原子的位置啊？感觉每次优化后的位置调整不是很多啊？2 .大家在算能带和态密度的时候的时候是用单点能计算给出能带还是几何优化给出能带啊？回答：问题1几何优化的目的是寻找压力最小的几何结构，原子位

31、置改变不多是因为你建的构形比较合理，比如直接从软件数据库中导入的结构在0压力下改变很小，如果加个压力就会变化大些。问题2一般都是优化的时候算，如果你想算特定构形的能带或其它性质，就可以用单点能；【表面模拟】castep计算中遇到的几个问题CODE:Copytocipboard)CODE:Copytoclpboard第一个问题：大家回答的差不多。你的结果肯定是有问题的，必须要调试到没这个警告出现.噌加bands的做法不是很可靠，最好是调整smearing,怎么调2最好的方法是测试，就是大小都调.不过你这里提示是要你减小这个数值，可以第二个向题，L每个位置都要算，通过能量高低判断最隹吸附位苴.对于

32、分子可能有多种组合的位苴，部要计算，不过最好参考下又献，一般是在高对称点咯！推荐下，计算的时候固定乂挫标，保持高度变化就。上！2,根绝原子半径判断，呵呵，在表面就不要考虑，如果你考虑他可能吸附到体内去，比如那里是fee位苴，有一个洞，那是可以考虑原子半径.因为半径大了，那个洞容不下他肯定不会进去咯.不过你把他放到洞里算下就知道哒.一般是用吸附能判断，吸附能为正就是可以吸附，为负不吸附.不过吸附能定义的计算方法要搞正确就是.吸附能定义为吸附前后物质总能量的变化，是结合能的反数，也就是吸附体系的能豌减去清洁表面的能量，再够去被吸附物质的能量.3,是否选择optimizedcell一般的体系都可以不

33、选择，选择上了计算量噌大太多.如果不选择，也不是位置固定.只是那个合子，或者说整体的体积和外形不变化而已.我是这样的理解的！4,怎么模拟表而一般采说计算表面吸附是不会固定下面的基地的，固定了可信度小，发出去自己也没底气人_人一般采用的2种办法，一种是楼上说的，固定的最下面层，然后上面和吸附原子接触的就“JLax，这个是相当于只考虑一个表面.固定的层数对半开.第二种方法是不固定层数，选择很大的层数，或者选择小一点的层数，固定中间层数，两面的原子层数不固定，这个就是有两个表面的情况还有少数是全都固定原子层数.对干楼主来说，caste的话，肯定在你计算的时候你选择的任务是只算能量，我估计你应该是选择

34、了结构优化,子位置的.第三个问题8选择自旋极化，推荐采用默认设置.正常情况下，郃会收敛到比较合理的值.除非自己了解很透彻，再自定义.PW91Perdew-Wanggeneralized-gradientapproximationPerdewandWang(1992)BPBeckeexchangeplusPerdewcorrelationBecke(1988),PerdewandWang(1992)PBEPerdew-Burke-ErnzerhofcorrelationPerdewetal.(1996)RPBERevisedPBEfunctionalbyHammeretal.Hammeretal

35、.(1999)HCTHHamprecht,Cohen,TozerandHandyfunctionalBoeseandHandy(2001)BLYPBeckeexchangeplusLee-Yang-ParrcorrelationBecke(1988),Leeetal.(1988)BOPBeckeOneParameterfunctionalTsunedaetal.(1999)VWN-BPBPfunctionalwiththelocalcorrelationreplacedbytheVWNfunctional.Voskoetal.(1980),Becke(1988),PerdewandWang(1

36、992)VWN:TheVosko-Wilk-Nusair(VWN)functionalisthemostpopularLSDcorrelationpotential.Itusesafittoaccuratenumericalresults(byCeperlyandAlder)ofauniformelectrongas.CeperleyandAlderperformedquantumMonteCarlocalculationsonauniformelectrongasatlowandhighspinlimitsforseveralelectrondensities.VWNusesthePadei

37、nterpolationproceduretofittheCAresultsforboththeparaandferrostatesandforlowandhighdensities.DMol3usesthebestVWN(socalled"Fit")parameters.PWC:ThePerdew-Wang(PWC)functionalisarecentparameterizationoftheCeperleyandAlderdata,whichcorrectssomeVWNproblemswithfitting.PWCisthedefaultfunctionalforD

38、Mol3calculations.Thelocalspin-density(LSD)approximationaccuratelypredictsstructures,vibrations,andrelativeenergiesofcovalentsystems;however,bondenergiesareseriouslyoverestimated.ThelocalDFTshouldnotbeusedforsystemswithweakbonds,suchashydrogenbonds.TheseproblemswiththeLSDmethodcanbecorrectedtoalargee

39、xtentbyusingtheso-calledgradient-corrected(ornonlocal)functionals.P91,BP,BLYP,BOP:DMol3supportsseveralnonlocalexchangeandcorrelationfunctionals.Themostpopular,theBeckeexchangefunctional(B88)isusedinconjunctionwiththePerdew-Wangcorrelationfunctional(BP)ortheLee-Yang-Parrcorrelationfunctional(BLYP).Th

40、eso-calledgeneralizedgradientcorrected(GGA)functional,byPerdewandWang(P91)wasderivedbyconsideringlowandhighdensityregimesandbyenforcingvarioussummationrules.PBE:ThePBE(Perdew,BurkeandEnzerhof)functional(1996,1997)isanotherGGAfunctionalinwhichalltheparameters(otherthenthoseinitsLDAcomponent)arefundam

41、entalconstants.TheexchangepartofthisfunctionalissimilartotheBeckeformula(1986),andthecorrelationpartisclosetothePerdew-Wangfunctional(1986).Thisfunctionalhasastrongphysicalbackground,reliablenumericalperformanceanditisfrequentlyusedinDFTcalculations.RPBE:Morerecently,Hammer,HansenandNorskov(1999)pro

42、posedamodifiedversionofthePBEformulathatimprovesconsiderablythermochemicalresults.Sofarthisfunctional,calledRPBE,hasbeenmainlyusedforsolids.HCTH:TheHCTHfunctional,namedfortheauthors'initials(Hamprechtetal.,1998),representsa"pragmatic"philosophyindesigningaDFTfunctional.Assumingthatthee

43、xactfunctionalwillneverbefound,theyproposeaflexibleformofgradientcorrectedfunctionalthatisfittedtothetrainingsetofmolecules.Thedefaultforthecurrentimplementation,thesocalledHCTH/407functional,thatwasobtainedbyfittingtothesetofthe407atomicandmolecularsystems(BoeseandHandy,2001).Thisfunctionalwasfoundtopredictamuchimprovedthermochemistryforinorganicandh