文件开头,首先是计算该体系的原子,本例中是碳纳米管,因此只有一种原子:C。文件内容如下:
<Sep06-2016> <18:16:36> DIRAC development version RunTime: Sep06-2016 18:16:36 Nodes: 1 Procs: 1 <Sep06-2016> <18:16:36> NORMAL TERMINATION <Sep06-2016> <18:16:36> END <Sep06-2016> <18:16:36> ADF development version RunTime: Sep06-2016 18:16:36 Nodes: 1 Procs: 1 <Sep06-2016> <18:16:36> Carbon (DZ, 1s frozen) <Sep06-2016> <18:16:36> RunType : CREATE <Sep06-2016> <18:16:36> Net Charge: 0 (Nuclei minus Electrons) <Sep06-2016> <18:16:36> Symmetry : ATOM Coordinates Atom X Y Z (Angstrom) 1.C 0.000000 0.000000 0.000000 <Sep06-2016> <18:16:36> >>>> CORORT <Sep06-2016> <18:16:36> >>>> CLSMAT <Sep06-2016> <18:16:36> >>>> ORTHON <Sep06-2016> <18:16:36> >>>> GENPT <Sep06-2016> <18:16:36> Block Length= 68 <Sep06-2016> <18:16:36> >>>> PTBAS <Sep06-2016> <18:16:36> >>>> CYCLE <Sep06-2016> <18:16:36> |Error| MaxErr Wt(A-DIIS) <Sep06-2016> <18:16:36> 1 0.00000000 0.00000000 <Sep06-2016> <18:16:36> 2 0.96064530 0.85132506 <Sep06-2016> <18:16:36> 3 0.70148959 0.61101908 100.0 <Sep06-2016> <18:16:36> 4 0.24621428 0.20460562 100.0 <Sep06-2016> <18:16:36> 5 0.19759987 0.16578230 100.0 <Sep06-2016> <18:16:36> 6 0.02764157 0.02297971 22.2 <Sep06-2016> <18:16:36> 7 0.00500151 0.00417208 3.2 <Sep06-2016> <18:16:36> 8 0.00011217 0.00009337 0.0 <Sep06-2016> <18:16:36> 9 0.00000023 0.00000019 0.0 <Sep06-2016> <18:16:36> 10 0.00000000 0.00000000 0.0 <Sep06-2016> <18:16:36> SCF converged <Sep06-2016> <18:16:36> 11 0.00000000 0.00000000 0.0 <Sep06-2016> <18:16:36> Solutions with partially occupied orbitals may not be <Sep06-2016> <18:16:36> lowest in energy. You might consider lowering the <Sep06-2016> <18:16:36> symmetry in the input and explicitly specifying integer <Sep06-2016> <18:16:36> occupations. In that case always check that you obtain <Sep06-2016> <18:16:36> an aufbau solution. <Sep06-2016> <18:16:36> >>>> POPAN <Sep06-2016> <18:16:36> >>>> DEBYE <Sep06-2016> <18:16:36> NORMAL TERMINATION <Sep06-2016> <18:16:36> END
如果体系没有进行分区,ADF默认以每个原子为一个分区,所以会先计算原子。这一段就是对原子进行DFT计算,求得原子的电子轨道,这些轨道将会作为基组,计算整个分子。如果体系分了片段,那么首先是计算各个片段,各个片段会有自己的*.logfile,那么整体的计算的时候,就不需要计算原子了,就会跳过这一段。 这些内容,左端列出了日期、时间,用户可以对照时间,看到计算的效率,可以估算完成计算需要多少时间。 从CYCLE这一行开始,表示计算进入自洽迭代(一般而言,自洽迭代是所有的计算中,最耗时的部分,当然原子的自洽迭代耗时很少),到SCF converged这一行,表示自洽迭代收敛,再求解一次自洽迭代方程,作为电子态的最终解。所以在SCF converged还有一行。实际上分子的计算也是类似的,下面可以看到。 上述内容一般不需要关注。 下面开始,是分子计算的内容:
<Sep06-2016> <18:16:36> ADF development version RunTime: Sep06-2016 18:16:36 Nodes: 1 Procs: 2
</color>Nodes: 1 Procs: 2表示这里的计算,耗用一个节点,共2个核。</color>
<Sep06-2016> <18:16:36> *** (NO TITLE) *** <Sep06-2016> <18:16:36> RunType : SINGLE POINT <Sep06-2016> <18:16:36> Spin-Orbit: Spin Restricted <Sep06-2016> <18:16:37> Net Charge: 0 (Nuclei minus Electrons) <Sep06-2016> <18:16:37> Symmetry : D(6D)
这一段,列出任务的类型,本例中是单点计算,使用了相对论Spin-Orbit,电荷为0,程序自动检测出点群为D6d。 下面列出分子的坐标(直接ctrl c拷贝坐标的这60行,在ADFinput中ctrl v粘贴,就可以将这60个原子形成的分子粘贴到ADFinput窗口里面):
<Sep06-2016> <18:16:37> >>>> FRAGM Coordinates Atom X Y Z (Angstrom) 1.C 1.661767 -1.661767 2.486510 2.C 2.270016 0.608249 2.486510 3.C 0.608249 2.270016 3.907370 4.C 0.608249 2.270016 -1.776070 5.C -2.270016 -0.608249 -4.617780 6.C -0.608249 2.270016 4.617780 7.C 1.661767 1.661767 0.355210 8.C 2.270016 0.608249 -1.776070 9.C -2.270016 0.608249 -2.486510 10.C -0.608249 -2.270016 -0.355210 11.C -2.270016 -0.608249 -0.355210 12.C 1.661767 -1.661767 -1.776070 13.C 0.608249 -2.270016 4.617780 14.C 0.608249 -2.270016 1.776070 15.C 1.661767 1.661767 -3.907370 16.C 0.608249 2.270016 -0.355210 17.C -2.270016 -0.608249 3.907370 18.C -1.661767 -1.661767 1.776070 19.C 1.661767 -1.661767 3.907370 20.C 1.661767 1.661767 1.776070 21.C -0.608249 -2.270016 -1.776070 22.C -0.608249 -2.270016 3.907370 23.C -1.661767 -1.661767 -3.907370 24.C -2.270016 0.608249 -3.907370 25.C 2.270016 -0.608249 1.776070 26.C 1.661767 1.661767 -2.486510 27.C 1.661767 1.661767 4.617780 28.C 0.608249 2.270016 -4.617780 29.C -1.661767 -1.661767 0.355210 30.C -0.608249 2.270016 1.776070 31.C -0.608249 2.270016 0.355210 32.C -2.270016 0.608249 1.776070 33.C 2.270016 -0.608249 0.355210 34.C -0.608249 -2.270016 2.486510 35.C -1.661767 -1.661767 4.617780 36.C -1.661767 1.661767 -4.617780 37.C 0.608249 -2.270016 -3.907370 38.C 2.270016 0.608249 3.907370 39.C -1.661767 1.661767 2.486510 40.C -1.661767 1.661767 3.907370 41.C -1.661767 -1.661767 -2.486510 42.C 0.608249 -2.270016 -2.486510 43.C -0.608249 2.270016 -2.486510 44.C 2.270016 0.608249 -0.355210 45.C 1.661767 -1.661767 -4.617780 46.C 2.270016 -0.608249 4.617780 47.C -2.270016 -0.608249 -1.776070 48.C 2.270016 -0.608249 -3.907370 49.C -2.270016 0.608249 0.355210 50.C 2.270016 0.608249 -4.617780 51.C 0.608249 -2.270016 0.355210 52.C -1.661767 1.661767 -1.776070 53.C -2.270016 -0.608249 2.486510 54.C -2.270016 0.608249 4.617780 55.C -0.608249 2.270016 -3.907370 56.C 0.608249 2.270016 2.486510 57.C 1.661767 -1.661767 -0.355210 58.C 2.270016 -0.608249 -2.486510 59.C -0.608249 -2.270016 -4.617780 60.C -1.661767 1.661767 -0.355210 <Sep06-2016> <18:16:37> >>>> CORORT <Sep06-2016> <18:16:37> >>>> CLSMAT <Sep06-2016> <18:16:37> >>>> ORTHON <Sep06-2016> <18:16:37> >>>> GENPT <Sep06-2016> <18:16:38> Block Length= 128 <Sep06-2016> <18:16:38> >>>> PTBAS
下面开始自洽迭代:
<Sep06-2016> <18:17:09> >>>> CYCLE <Sep06-2016> <18:17:13> |Error| MaxErr Wt(A-DIIS) <Sep06-2016> <18:17:13> 1 7.86400595 1.00449720 <Sep06-2016> <18:17:15> 2 2.19283096 0.14870767 100.0 <Sep06-2016> <18:17:18> 3 1.44909054 0.10612198 100.0 <Sep06-2016> <18:17:20> 4 1.14359639 0.06934509 69.0 <Sep06-2016> <18:17:23> 5 1.80185993 0.13991257 100.0 <Sep06-2016> <18:17:26> 6 0.68932336 0.05211586 51.6 <Sep06-2016> <18:17:28> 7 0.41756452 0.04369477 43.1 <Sep06-2016> <18:17:31> 8 0.13604083 0.01719497 16.4 <Sep06-2016> <18:17:33> 9 0.03780013 0.00412614 3.2 <Sep06-2016> <18:17:36> 10 0.00849766 0.00092242 0.0 <Sep06-2016> <18:17:39> 11 0.00412811 0.00034462 0.0 <Sep06-2016> <18:17:42> 12 0.00061920 0.00004973 0.0 <Sep06-2016> <18:17:44> 13 0.00007056 0.00000708 0.0 <Sep06-2016> <18:17:47> 14 0.00002504 0.00000241 0.0 <Sep06-2016> <18:17:50> 15 0.00000489 0.00000030 0.0 <Sep06-2016> <18:17:51> SCF converged <Sep06-2016> <18:17:52> 16 0.00000015 0.00000002 0.0
迭代15次后,收敛,再求解一次迭代方程,作为最终的电子态的解。 下面计算总能量、电子布居、结果分析等等,之后以三种单位,分别列出总的Bonding Energy。关于Bonding Energy与Total Energy的区别,参考:Total Bonding Energy与Total Energy的差别与用法。
<Sep06-2016> <18:17:55> >>>> TOTEN <Sep06-2016> <18:18:12> >>>> POPAN <Sep06-2016> <18:18:12> >>>> DEBYE <Sep06-2016> <18:18:12> >>>> AMETS <Sep06-2016> <18:18:12> Bond Energy -18.93003963 a.u. <Sep06-2016> <18:18:12> Bond Energy -515.11258748 eV <Sep06-2016> <18:18:12> Bond Energy -11878.78 kcal/mol <Sep06-2016> <18:18:12> NORMAL TERMINATION <Sep06-2016> <18:18:12> END
到这里为止,计算结束,开始生成结果文件——*.t21(或者TAPE21,二者是一回事,改个名字就一样了):
Job nano has finished
显示这一行的时候,表示*.t21文件生成完毕,计算完全结束。