关注：

1） 在优化过程中保持对称性不变，能否通过固定原子坐标优化实现？

2） 表面的构造

3） 如何判断优化过程是否收敛，如下述情形

1.如何判断优化过程是否收敛，如下述情形

$ grep CPU OUT*4
OUTCAR_0-04:                  Total CPU time used (sec):     1074.722
OUTCAR_100-04:                  Total CPU time used (sec):      695.405
OUTCAR_10-04:                  Total CPU time used (sec):      916.000
OUTCAR_150-04:                  Total CPU time used (sec):      809.826
OUTCAR_200-04:                  Total CPU time used (sec):      925.247
OUTCAR_300-04:                  Total CPU time used (sec):    26966.850   显然没有收敛
OUTCAR_50-04:                  Total CPU time used (sec):      702.251

grep F= OSZICAR_300-04
  1 F= -.59908513E+02 E0= -.59901487E+02  d E =-.599085E+02  mag=    22.2514
  2 F= -.59909170E+02 E0= -.59900929E+02  d E =-.657653E-03  mag=    22.2290
  3 F= -.59909715E+02 E0= -.59898200E+02  d E =-.120211E-02  mag=    22.1501
  4 F= -.59907302E+02 E0= -.59894637E+02  d E =0.121097E-02  mag=    22.0980
  5 F= -.59909552E+02 E0= -.59897659E+02  d E =-.103969E-02  mag=    22.1327

......

2. 如何固定原子坐标优化

POSCAR 里各类原子原子数的下一行写selective dynamics
 
在各个原子坐标后 指定该原子是否固定
 
固定死的 写上 F F F

放开的   写上 T T T

如果第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

diamond
0.529177249
  14.276296295569   0.000000000000   0.000000000000
   2.951447009758   3.516467367877   0.000000000000
   1.475764202501   1.758196839974   3.975852328470
C
4
Selective
Direct
  -0.000000009721   0.000001503872   0.000000126860 T T F
  -0.000000008663   0.333332080250   0.333332843477 T T F
   0.500000008663  -0.333332080250  -0.333332843477 T T F

   0.500000009721  -0.000001503872  -0.000000126860 T T F