| 查看: 2177 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
mashuai222银虫 (初入文坛)
|
[求助]
求助：VASP中用ISIF=3优化晶格常数时，如何只让一个方向变化，其他方向不变 已有2人参与
|
||
|
小弟新手上路，ＶＡＳＰ中如何设置参数可以在结构优化时只改变一个方向的晶格常数，其他两个方向的不变？ 不知哪位高手有空指教，小弟不胜感激！ |
回复此楼» 收录本帖的淘帖专辑推荐
 精华网帖收集 |
» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6519.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.8小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
3楼2016-03-31 21:25:20
lkccrr
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 3254.5
- 散金: 237
- 帖子: 89
- 在线: 443.6小时
- 虫号: 2027965
- 注册: 2012-09-25
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
mashuai222: 金币+10, ★★★很有帮助 2016-04-01 08:44:26
感谢参与，应助指数 +1
mashuai222: 金币+10, ★★★很有帮助 2016-04-01 08:44:26
|
这个得VASP编译的时候把constr_cell_relax.F这个文件里面的一些参数做一下更改，比如：!----------------------------------------------------------------------- ! ! At present, VASP does not allow to relax the cellshape selectively ! i.e. for instance only cell relaxation in x direction. ! To be more precisse, this behaviour can not be achived via the INCAR ! or POSCAR file. ! However, it is possible to set selected components of the stress tensor ! to zero. ! The most conveninent position to do this is the routines ! CONSTR_CELL_RELAX (constraint cell relaxation). ! FCELL contains the forces on the basis vectors. ! These forces are used to modify the basis vectors according ! to the following equations: ! ! A_OLD(1:3,1:3)=A(1:3,1:3) ! F90 style ! DO J=1,3 ! DO I=1,3 ! DO K=1,3 ! A(I,J)=A(I,J) + FCELL(I,K)*A_OLD(K,J)*STEP_SIZE ! ENDDO ! ENDDO ! ENDDO ! where A holds the basis vectors (in cartesian coordinates). ! !----------------------------------------------------------------------- SUBROUTINE CONSTR_CELL_RELAX(FCELL) USE prec REAL(q) FCELL(3,3) ！ orthorhombic system relaxation in a and b Saa=FCELL(1,1) Sbb=FCELL(2,2) FCELL=0 FCELL(1,1)=Saa FCELL(2,2)=Sbb ! 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 ! relaxation in z direction only ! SAVE=FCELL(3,3) ! FCELL=0 ! F90 style: set the whole array to zero ! FCELL(3,3)=SAVE RETURN END SUBROUTINE |
2楼2016-03-31 16:13:54
20