| 查看: 2197 | 回复: 2 | |||
mashuai222银虫 (初入文坛)
|
[求助]
求助：VASP中用ISIF=3优化晶格常数时，如何只让一个方向变化，其他方向不变 已有2人参与
|
|
小弟新手上路，ＶＡＳＰ中如何设置参数可以在结构优化时只改变一个方向的晶格常数，其他两个方向的不变？ 不知哪位高手有空指教，小弟不胜感激！ |
回复此楼» 收录本帖的淘帖专辑推荐
 精华网帖收集 |
» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有10人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
青基如何破局
已经有4人回复
-
基金不中，共勉
已经有8人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
lkccrr
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 3277
- 散金: 379
- 帖子: 98
- 在线: 444.4小时
- 虫号: 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
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6523.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.8小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
3楼2016-03-31 21:25:20
投票:
10