如果out文件名字类似*.out那么打开的时候,会使用ADF的Output窗口打开,窗口自带菜单栏,可以在Properties、Response Properties等下拉菜单选取大部分需要的信息。
但如果是结构优化,我们有时候希望看看优化的每一个结构,电子的能级、占据方式,就需要从out文件里面,手工去翻阅了。
逐渐下翻,可以看到第一次出现SCF字样,表示对第一个结构进行自洽迭代:
S C F CYCLE 1 WARNING: DIIS procedure will be reset orbitals (Q,E): --------------- A :465...484 ( 2.00 -0.2569) ( 2.00 -0.2539) ( 2.00 -0.2467) ( 2.00 -0.2465) ( 2.00 -0.2439) ( 2.00 -0.2400) ( 2.00 -0.2399) ( 2.00 -0.2388) ( 2.00 -0.2385) ( 2.00 -0.2373) ( 0.00 -0.1668) ( 0.00 -0.1642) ( 0.00 -0.1511) ( 0.00 -0.1510) ( 0.00 -0.1478) ( 0.00 -0.1475) ( 0.00 -0.1401) ( 0.00 -0.1356) ( 0.00 -0.1356) ( 0.00 -0.1294) CYCLE 2 d-Pmat mean: 0.49E-02 imax= 1718: -0.85E-01 orbitals (Q,E): --------------- A :465...475 ( 2.00 -0.2427) ( 2.00 -0.2411) ( 2.00 -0.2351) ( 2.00 -0.2349) ( 2.00 -0.2266) ( 2.00 -0.2264) ( 2.00 -0.2264) ( 2.00 -0.2261) ( 2.00 -0.2260) ( 2.00 -0.2247) ( 0.00 -0.1633) CYCLE 3 d-Pmat mean: 0.99E-02 imax= 1718: -0.29E+00 …………………………
上面这些这就是自洽迭代的内容,和*.logfile里面是可以对应起来的(参考:结构优化的*.logfile)。我们拿CYCLE 2来说,就是迭代到第二步的电子能级、占据情况。( 2.00 -0.2247) ( 0.00 -0.1633)表示HOMO为-0.2247 hartree,LUMO为-0.1633Hartree。
往下翻,最后一轮,例如:
………………………… CYCLE 12 *** Setting up for NEW gradients in focky *** Using FIT density in focky orbitals (Q,E): --------------- A :465...484 ( 2.00 -0.2103) ( 2.00 -0.2076) ( 2.00 -0.2047) ( 2.00 -0.2045) ( 2.00 -0.1944) ( 2.00 -0.1938) ( 2.00 -0.1936) ( 2.00 -0.1928) ( 2.00 -0.1925) ( 2.00 -0.1905) ( 0.00 -0.1437) ( 0.00 -0.1374) ( 0.00 -0.1293) ( 0.00 -0.1289) ( 0.00 -0.1209) ( 0.00 -0.1205) ( 0.00 -0.1019) ( 0.00 -0.0971) ( 0.00 -0.0970) ( 0.00 -0.0956)
可以看到HOMO是:-0.1905Hartree,LUMO: -0.1437Hartree。
往下翻,以此出现第二个结构的自洽迭代,内容于此相似。从而可以看到结构优化过程中能级的变化。如果出现HOMO、LUMO简并,也能因此观察到。这在激发态结构优化中,经常需要注意是否出现了这种情况,如果是出现这种情况,激发能往往变得很小,接近0.01eV的量级。
在level图里面,或者out里面,经常可以看见1f、1d、1p、2p、2d这些轨道编号,这些编号和我们常说的3d轨道不完全一致。
其中1f表示第一个f轨道。如果没有冻芯,那么1f实际上就是我们常说的4f轨道,2f就是5f轨道。如果使用了冻芯,加入4f轨道被冻心了,那么1f就是指5f轨道了,类似的,如果3p轨道被冻芯了,1p就是指4p轨道;以此类推。