VASP并行可执行软件包,可对晶胞参数进行部分优化
下载地址:http://ishare.iask.sina.com.cn/f/17331535.html
在优化表面,二维或一维体系时,必须固定某个或某几个晶胞方向不优化。
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 返回小木虫查看更多
你好,我优化一个2D构型,(011),x轴不优化,放宽y,z轴优化,但最终x,y,z都在优化,似乎这于这种情况不适合,你检查下代码是哪有问题?当然如果是(110)这样的话,程序就没有问题
,
你好,我优化一个2D构型,(011),x轴不优化,放宽y,z轴优化,但最终x,y,z都在优化,似乎这于这种情况不适合,你检查下代码是哪有问题?当然如果是(110)这样的话,程序就没有问题。
太牛了,谢谢
非常感谢!
同问,如果非正交呢
顶顶更健康
我也发现了 011肯定是有问题的 代码有问题 暂时还没学程序,看不懂哪里有问题
你好~我也遇到了这个问题,请问你是如何解决的,谢谢大神虫友~