24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2950  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wl7171759

新虫 (小有名气)

[求助] VASP做结构优化时用ISIF=3只想对晶格常数中的一个维度放开经行优化已有2人参与

就是用ISIF=3时想固定住X,Y两个方向,而单单只对Z方向的晶格常数经行优化 该如何实现?希望大家不吝赐教
回复此楼

» 收录本帖的淘帖专辑推荐

有用的

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

引用回帖:
3楼: Originally posted by golddoushi at 2014-02-11 15:15:45
修改constr_cell_relax.F, 你把里面的代码复制一段,改个ISIF = 100啥的。。。再把那个小循环里的判定条件改了,再make一下。很简单的,想咋弛豫咋弛豫。我记得以前有个帖子把石墨弛豫成金刚石来着

你好!我现在想固定面内晶格参数,a、b、alpha,弛豫面外晶格参数c、beta、gemma。我的想法是将constr_cell_relax.F中的语句修改成如下:
       FCELL(1,1)=0   
       FCELL(1,2)=0   
       FCELL(2,1)=0   
       FCELL(2,2)=0
或者
       SAVE1=FCELL(1,3)
       SAVE2=FCELL(2,3)
       SAVE3=FCELL(3,1)
       SAVE4=FCELL(3,2)
       SAVE5=FCELL(3,3)
       FCELL=0   ! F90 style: set the whole array to zero
       FCELL(1,3)=SAVE1
       FCELL(2,3)=SAVE2
       FCELL(3,1)=SAVE3
       FCELL(3,2)=SAVE4
       FCELL(3,3)=SAVE5

重新编译vasp之后用ISIF=3,结构还是所有的晶格参数都变了。请问我的方法有什么问题吗?谢谢!
6楼2014-11-13 08:18:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

引用回帖:
3楼: Originally posted by golddoushi at 2014-02-11 15:15:45
修改constr_cell_relax.F, 你把里面的代码复制一段,改个ISIF = 100啥的。。。再把那个小循环里的判定条件改了,再make一下。很简单的,想咋弛豫咋弛豫。我记得以前有个帖子把石墨弛豫成金刚石来着

困扰我的问题被我自己想明白了。我的方法是对的,但是我的结构建的不对。我的结构中三个夹角都不是90度。我之前建的结构是c轴平行于z轴(MS里面的默认设置),其实应该让a、b轴在xOy平面上。按照后一种方法建的模型就能够按照我想要的弛豫方式弛豫了。
感觉vasp里面是不管你建的胞是什么样的点阵,它内部默认的是一个正交点阵。constr_cell_relax.F能够控制的实际上是这个正交点阵。我之前建的胞中a、b轴偏离了正交点阵的xOy平面,所以我实际上固定的只是a、b轴在xOy平面上的投影。
7楼2014-11-13 15:59:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wl7171759 的主题更新
信息提示
请填处理意见