TDDFT计算软件Octopus学习笔记(二)：基态计算(Si)

本文以硅为例，来进行基态计算。

inp 文件

CalculationMode = gs
PeriodicDimensions = 3
Spacing = 0.5
%LatticeVectors
  0.0 | 0.5 | 0.5 
  0.5 | 0.0 | 0.5
  0.5 | 0.5 | 0.0
%
a = 10.18
%LatticeParameters
 a | a | a
%
%ReducedCoordinates
 "Si" | 0.0 | 0.0 | 0.0 
 "Si" | 1/4 | 1/4 | 1/4 
%
nk = 2
%KPointsGrid
  nk |  nk |  nk
 0.5 | 0.5 | 0.5
 0.5 | 0.0 | 0.0
 0.0 | 0.5 | 0.0
 0.0 | 0.0 | 0.5
%
KPointsUseSymmetries = yes
ExtraStates = 1
Output = dos

命令详解

CalculationMode = gs

此命令用于指定计算模式，这里设定为基态计算ground state(gs)，还有许多其选项，例如 unocc(unoccupied)非占据态计算， td(time dependent)时间依赖计算，等等，共14种选项。默认为gs。

PeriodicDimensions = 3

这个变量指定了你所考虑的应用周期性的维度数量，这里设定为三维，这样就可以在所有方向设定周期性边界条件，意味着我们所考虑的是一个有完整周期性的无限晶体。如果不设置，则默认为0，即没有周期性的孤立体系，例如分子。

Spacing = 0.5

Octopus是一款实空间，实时间的计算软件，所谓实空间即体现在这里。通过设定网格将实空间离散化，而间距(spacing)则用来控制网格之间的距离，自然，间距越小，则离散化的精度越高，相应的计算成本增加。单位默认为原子单位，即波尔(Bohr).

$1Bohr=$