单点计算的*.logfile

文件开头,首先是计算该体系的原子,本例中是碳纳米管,因此只有一种原子: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文件生成完毕,计算完全结束。