求助：vasp优化二维晶体的constr_cell_relax.F

nextnest

各位老师：       大家好！目前尝试用vasp5.4.1优化二维晶体，限制c轴，优化a和b轴，参考了网上的constr_cell_relax.F文件进行了编译，但是编译报错如下：
make[2]: *** [constr_cell_relax.o] Error 1。
       想请教一下大家如何修改constr_cell_relax.F文件实现a，b的同时优化？谢谢！


附：修改后的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

ggdh

为啥不直接isif=4？

nextnest

ggdh 发表于 2017-11-3 12:14
为啥不直接isif=4？

在计算中，需要优化a，b，固定c，同时还要优化原子坐标，可能会改变晶胞体积，ISIF=4限制了不改变晶胞体积。

ggdh

nextnest 发表于 2017-11-3 12:19
在计算中，需要优化a，b，固定c，同时还要优化原子坐标，可能会改变晶胞体积，ISIF=4限制了不改变晶胞体 ...

计算二维材料，C方向大小无所谓，固定体积，让ab任意变动，让c去代偿ab的变动。

nextnest

明白了，我试一试，多谢！