24小时热门版块排行榜    

CyRhmU.jpeg
查看: 10595  |  回复: 17
【奖励】 本帖被评价14次,作者chemvasp增加金币 10.4
本帖产生 1 个 1ST强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chemvasp

新虫 (初入文坛)


[资源] 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
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理计算-tg vasp计算资料 第一性原理计算经验 材料计算模拟软件
VASP资料 量化软件 有效质量 网贴收集
vasp soft study vasp

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chemvasp 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见