| 查看: 2117 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
mashuai222银虫 (初入文坛)
|
[求助]
求助：VASP中用ISIF=3优化晶格常数时，如何只让一个方向变化，其他方向不变 已有2人参与
|
||
|
小弟新手上路，ＶＡＳＰ中如何设置参数可以在结构优化时只改变一个方向的晶格常数，其他两个方向的不变？ 不知哪位高手有空指教，小弟不胜感激！ |
回复此楼» 收录本帖的淘帖专辑推荐
 精华网帖收集 |
» 猜你喜欢
0805 316求调剂
已经有6人回复
-
085601求调剂总分293英一数二
已经有3人回复
-
考研调剂
已经有3人回复
-
08工学调剂
已经有17人回复
-
340求调剂
已经有4人回复
-
311求调剂
已经有3人回复
-
食品专硕 一志愿双一流 328
已经有4人回复
-
材料调剂
已经有6人回复
-
300分，材料，求调剂，英一数二
已经有5人回复
-
求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分
已经有5人回复
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6519.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.7小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
3楼2016-03-31 21:25:20
lkccrr
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 3421.3
- 散金: 134
- 帖子: 83
- 在线: 442.7小时
- 虫号: 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
