| 查看: 11584 | 回复: 17 | |||||||||||||
| 【奖励】 本帖被评价14次，作者chemvasp增加金币 10.4 个 | |||||||||||||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | |||||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||||
[资源]
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 |
» 猜你喜欢
287求调剂
已经有5人回复
-
274求调剂
已经有4人回复
-
308求调剂
已经有14人回复
-
本科郑州大学，一志愿华东师范大学282求调剂
已经有24人回复
-
化学308分求调剂
已经有18人回复
-
求材料调剂，一志愿郑州大学289分
已经有19人回复
-
277求调剂
已经有5人回复
-
一志愿211，化学310分，本科重点双非，求调剂
已经有10人回复
-
266求调剂
已经有14人回复
-
申博
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于vasp收敛和优化
已经有24人回复
- ansys中接触单元的参数如何定义
已经有10人回复
- vasp.5.2.12 说明书
已经有4人回复
- 正交设计试验后，是否以及怎么样进行后续的进一步优化？
已经有6人回复
- 在vasp计算里，有没有好的只优化原子位置和晶胞体积的方法？
已经有8人回复
- vasp 那些参数设置可以影响到磁性计算的磁矩
已经有12人回复
- dmol在对晶胞参数进行优化时是否应该固定晶胞角度（fix angles）
已经有12人回复
- 利用CASTEP对建立的晶胞进行几何优化的问题 （急！！！）
已经有13人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
- 【求助】能否发一份用vasp计算超晶胞某个表面的输入文件给我
已经有7人回复
- 【求助】Castep中结构优化后晶胞参数差别较大 vasp高手也帮我看一下
已经有22人回复
- 【求助】vasp编译时，报错：cannot find -lmkl_em64t
已经有11人回复
5楼2011-07-30 15:47:27
2楼2011-07-22 17:08:16
4楼2011-07-23 12:21:24
6楼2012-05-22 11:17:08
