24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3143  |  回复: 7

wl7171759

新虫 (小有名气)

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

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

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

有用的

» 猜你喜欢

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

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

甲斐之虎

银虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
我想,这个是不是要用脚本来实现,固定xy的晶格常数,已z方向为变量来优化……
ISIF=3貌似都变吧……
别人笑我太疯颠,我笑他人看不穿。不见五陵豪杰墓,无花无酒锄作田。
2楼2014-02-11 14:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

golddoushi

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
fzx2008: 金币+2, 谢谢指导 2014-02-13 12:59:30
wl7171759: 金币+10, ★★★很有帮助, 谢谢你的热心帮助! 2014-02-14 20:46:22
修改constr_cell_relax.F, 你把里面的代码复制一段,改个ISIF = 100啥的。。。再把那个小循环里的判定条件改了,再make一下。很简单的,想咋弛豫咋弛豫。我记得以前有个帖子把石墨弛豫成金刚石来着
3楼2014-02-11 15:15:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wl7171759

新虫 (小有名气)

引用回帖:
2楼: Originally posted by 甲斐之虎 at 2014-02-11 14:59:52
我想,这个是不是要用脚本来实现,固定xy的晶格常数,已z方向为变量来优化……
ISIF=3貌似都变吧……

ISIF=3是都变..所以我来想大家请教的
4楼2014-02-11 16:23:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wl7171759

新虫 (小有名气)

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

嗯  我在之前的帖子也看到过这个,不过我基本不懂程序...总之谢谢你哈
5楼2014-02-11 16:27:23
已阅   回复此楼   关注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的回帖

gaozhibin

铜虫 (正式写手)

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

您说的把石墨弛豫成金刚石来着 的帖子能 给个链接吗?

谢谢大神
加油努力,坚持爬行的小乌龟
8楼2015-12-22 19:31:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wl7171759 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[基金申请] 今年审到国自然15份,谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +7 水和泥不是水泥 2026-05-13 10/500 2026-05-18 07:50 by 水和泥不是水泥
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 7/350 2026-05-17 19:42 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-17 19:39 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-17 19:36 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:06 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:11 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:57 by ue3ir18jc3
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见