| 查看: 2078 | 回复: 2 | ||||
mashuai222银虫 (初入文坛)
|
[求助]
求助:VASP中用ISIF=3优化晶格常数时,如何只让一个方向变化,其他方向不变 已有2人参与
|
|
小弟新手上路,VASP中如何设置参数可以在结构优化时只改变一个方向的晶格常数,其他两个方向的不变? 不知哪位高手有空指教,小弟不胜感激! |
» 收录本帖的淘帖专辑推荐
精华网帖收集 |
» 猜你喜欢
交叉科学部支持青年基金,对三无青椒是个机会吗?
已经有5人回复
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
招博士
已经有6人回复
限项规定
已经有8人回复
国家级人才课题组招收2026年入学博士
已经有5人回复
Fe3O4@SiO2合成
已经有6人回复
青年基金C终止
已经有4人回复
青椒八年已不青,大家都被折磨成啥样了?
已经有7人回复
26申博求博导推荐-遥感图像处理方向
已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生(考核制)—机械、材料、力学方向
已经有3人回复

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
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6519.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.7小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
3楼2016-03-31 21:25:20













回复此楼