应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP做结构优化时用ISIF=3只想对晶格常数中的一个维度放开经行优化
81/1返回列表
查看: 2591  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wl7171759

新虫 (小有名气)

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

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

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

有用的

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-02-11 09:22:34
已阅   回复此楼   关注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一下。很简单的，想咋弛豫咋弛豫。我记得以前有个帖子把石墨弛豫成金刚石来着
赞一下(1人) 回复此楼
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 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 饶议：什么制度能保障大学普通教师不用为领导拎包，不用看领导脸色 +9 zju2000 2024-06-12 15/750 2024-06-16 20:45 by 北溟鱼1318
[基金申请] F口信息学部拿面上，大概需要什么样的成果 +3 _奋黎_ 2024-06-16 4/200 2024-06-16 20:21 by _奋黎_
[考博] 2025考博 +6 自强不息a?a 2024-06-15 8/400 2024-06-16 19:39 by 鱼翔浅底1
[基金申请] 2024国社科通讯评审 +9 qsd10086 2024-06-13 17/850 2024-06-16 19:14 by onroad24
[教师之家] 请问事业编制和年薪制冲突吗？ +8 ZHONGWU_U 2024-06-14 8/400 2024-06-16 17:37 by spj860iuy
[教师之家] 每次骚扰女学生的都是院系领导，而不是普通教师，小编们要注意措辞正确 +9 zju2000 2024-06-15 11/550 2024-06-16 14:49 by appleapple2
[找工作] 成都产品质量检测研究院 200+3 鲸鱼663 2024-06-11 9/450 2024-06-16 10:08 by SNaiL1995
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[硕博家园] 关于硕博连读的一些疑问？ +4 Lwenter 2024-06-14 4/200 2024-06-14 14:32 by ou0551
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 化学会年会改成一年开一次了？ +5 babu2015 2024-06-12 13/650 2024-06-14 10:46 by babu2015
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[论文投稿] 最近写了一篇控制优化领域的文章，可以投哪里啊？有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
[硕博家园] 机械研究生如何拿到年薪40+w +13 阿巴阿巴哦哦 2024-06-11 15/750 2024-06-13 15:40 by 113745685
[考博] 博导选择 +3 bing85977 2024-06-12 3/150 2024-06-13 15:34 by 我是邱尧
[硕博家园] 考博应该从哪里开始准备 +6 亮姐2019 2024-06-11 8/400 2024-06-13 12:53 by 帕秋莉GO
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
信息提示
关闭
请填处理意见
关闭