VASP中原子位置的指定（POSCAR文件详解）

POSCAR = position + CAR

  第1行：任意文字注释

  第2行：晶格常数，单位A，后面所有的长度值得自原值除以此值
  a=b=c时取a即可，否则个人习惯取三者最大
  若取负值，则为晶胞体积，单位A3

  第3-5行：定义晶矢 参见《固体量子化学——材料化学的理论基础》赵成大

  如 对于正交晶体 a=20.022 b=19.899 c=13.383 α=β=γ=90
  可以这样定义
  20.022
  1.00000 0.00000 0.00000
  0.00000 0.99386 0.00000
  0.00000 0.00000 0.66841
  又如 对于面心立方晶体 a=b=c=3.57 α=β=γ=90
  可以定义如下
  3.57
  0.0 0.5 0.5 (1/2(b+c))
  0.5 0.0 0.5 (1/2(a+c))
  0.5 0.5 0.0 (1/2(a+b))

  第6行：每种元素的原子数，特别注意顺序,要与下面的坐标顺序以及POTCAR中
  的顺序一致

  第7行：可省略，无需空行。
  做动力学时，是否需要固定部分离子的坐标。若是，此行以'S'或者's'首字即可。

  第8行开始为离子的坐标，格式为
  option line
  coordinate1 of element1
  coordinate2 of element1
  ...
  coordinateN of element1
  option line
  coordinate1 of element2
  coordinate2 of element2
  ...
  coordinateM of element2
  ...
  其中，option line指定输入坐标的格式，除了第一个以外，如果后面的输入格式同前，

  则都可以无空行省略。
  option line可指定的输入坐标格式有两种
  'D'or'd' for direct mode
  'C'or'c'or'K'or'k' for cartesian mode
  顾名思义，前者是定义在三个晶矢方向上的坐标
  R=R1×x+R2×y+R3×z R1,R2,R3为前面的晶矢，x,y,z为输入的三个坐标，R为坐标位矢
  而后者只是简单的将直角坐标除以前面第二行定义的晶胞常数
  两者可以混用，但不推荐。

  如果第7行设定了S(Selective Dynamic),则可以用以下形式定义各坐标是否可以移动
  Selective dynamics
  Cartesian
   0.00 0.00 0.00 T T F
   0.25 0.25 0.25 F F F