VASP单一轴优化

在优化表面，二维或一维体系时，必须固定某个或某几个晶胞方向不优化。
VASP本身没有这个功能，但是提供了constr_cell_relax.F文件来设置。
根据文件中的提示，添加几行代码后，就可以任意固定晶胞方向了。

使用方法：
在目录下建立文件OPTCELL，第一行输入三个数字（比如110），
数字1表示优化，数字0表示不优化；三个数字依次对应于
xyz方向。比如110表示z轴不优化；001表示x和y轴都不优化。
注1. 无OPTCELL文件则完全是原版软件的功能；
注2. 只适用与正交晶胞，所有角都是直角。

附：修改后的constr_cell_relax.F文件
SUBROUTINE CONSTR_CELL_RELAX(FCELL)
USE PREC
REAL(Q) FCELL(3,3), SAVE(3)
LOGICAL FILFLG
INTEGER ICELL(3)

  INQUIRE(FILE='OPTCELL',EXIST=FILFLG)
  IF (FILFLG) THEN
    OPEN(67,FILE='OPTCELL',FORM='FORMATTED',STATUS='OLD')
      READ(67,"(3I1)") (ICELL(I),I=1,3)
    CLOSE(67)
    DO I=1,3
      SAVE(I)=FCELL(I,I)
    ENDDO
    FCELL=0.0d0
    DO I=1,3
      IF (ICELL(I)==1) FCELL(I,I)=SAVE(I)
    ENDDO
  ENDIF

! just one simple example
! relaxation in x directions only
! SAVE=FCELL(1,1)
! FCELL=0 ! F90 style: set the whole array to zero
! FCELL(1,1)=SAVE

  RETURN
  END SUBROUTINE