24小时热门版块排行榜    

查看: 3241  |  回复: 8

wang-lf09

木虫 (小有名气)

[求助] 在vasp计算里,有没有好的只优化原子位置和晶胞体积的方法?

只想优化原子位置和晶胞体积,但ISIF没有对应的设置。在用ISIF=7只优化体积的时候,优化了一段时间不能优化了自己停掉了(估计是原子位置不合理,但我也不知道合理位置是什么样的,没法调整啊,要是vasp有同时优化原子位置和晶胞体积的功能就好了),不知道虫友有没有遇到类似情况?我现在是手动调整晶胞三个方向的长度,然后再优化,如此循环使三个基矢方向的应力小到一定程度为止,但这种方法太费劲了,一两天也做不好一个小例子。虫友们有什么好办法啊?速来救火啊。先谢啦!
回复此楼

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

vasp

» 猜你喜欢

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

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

likai0106

金虫 (小有名气)

【答案】应助回帖

★ ★
fzx2008(金币+2): 谢谢回帖交流! 2011-11-23 14:53:57
wang-lf09(金币+5): 2011-11-25 09:57:34
引用回帖:
1楼: Originally posted by wang-lf09 at 2011-11-23 12:25:26:
只想优化原子位置和晶胞体积,但ISIF没有对应的设置。在用ISIF=7只优化体积的时候,优化了一段时间不能优化了自己停掉了(估计是原子位置不合理,但我也不知道合理位置是什么样的,没法调整啊,要是vasp有同时优化 ...

是这样你可以试试利用脚本来解决你想做的情况,写出一个循环脚本给出一系列的值,之后进行计算。脚本形式类似与loop.sh优化晶格常数的脚本,如果你的a,b,c三个方向有比例关系容易些,如果不成比例的话,可以尝试对其分别进行修改。但是如果你想做固定方向的优化的话,比如a,b不动只优化c方向怕是要改vasp原代码重新编译vasp了。
希望对你有帮助。
2楼2011-11-23 14:17:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


youzhizhe(金币+1): 谢谢交流。 2011-11-24 23:16:55
引用回帖:
1楼: Originally posted by wang-lf09 at 2011-11-23 12:25:26:
只想优化原子位置和晶胞体积,但ISIF没有对应的设置。在用ISIF=7只优化体积的时候,优化了一段时间不能优化了自己停掉了(估计是原子位置不合理,但我也不知道合理位置是什么样的,没法调整啊,要是vasp有同时优化 ...

你把你的计算的错误贴出来吧 不至于什么都没有停掉的,
嘟啦啦嘟啦啦嘟
3楼2011-11-23 14:31:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wang-lf09

木虫 (小有名气)

引用回帖:
3楼: Originally posted by 贺仪 at 2011-11-23 14:31:26:
你把你的计算的错误贴出来吧 不至于什么都没有停掉的,

forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
vasp               000000000079FFE6  Unknown               Unknown  Unknown
vasp               00000000007A1B45  Unknown               Unknown  Unknown
vasp               00000000007A061E  Unknown               Unknown  Unknown
vasp               00000000009CA153  Unknown               Unknown  Unknown
vasp               00000000004570E0  Unknown               Unknown  Unknown
vasp               00000000004382FC  Unknown               Unknown  Unknown
libc.so.6          00002AF5EC478184  Unknown               Unknown  Unknown
vasp               0000000000438229  Unknown               Unknown  Unknown
--------------------------------------------------------------------------
mpirun has exited due to process rank 1 with PID 5940 on
node node8 exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------

报错如上所述,烦请帮忙看看,谢谢!
4楼2011-11-24 11:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likai0106

金虫 (小有名气)

【答案】应助回帖


youzhizhe(金币+1): 谢谢交流。 2011-11-24 23:17:00
引用回帖:
4楼: Originally posted by wang-lf09 at 2011-11-24 11:00:41:
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
vasp               000000000079FFE6  Unknown               U ...

第一你看看你的内存是否够用。
第二看看你的权限有没有问题。
5楼2011-11-24 15:14:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


youzhizhe(金币+1): 谢谢交流。 2011-11-24 23:17:04
wang-lf09(金币+2): 2011-11-25 09:58:12
引用回帖:
4楼: Originally posted by wang-lf09 at 2011-11-24 11:00:41:
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
vasp               000000000079FFE6  Unknown               U ...

上面的问题很多错误都是这样提示的。打开OUTCAR把最后面的东西复制出来,或者你找到提示错误的地方
嘟啦啦嘟啦啦嘟
6楼2011-11-24 15:21:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuvivi

金虫 (小有名气)

【答案】应助回帖

youzhizhe: 非应助回帖 2011-11-24 23:17:29
引用回帖:
2楼: Originally posted by likai0106 at 2011-11-23 14:17:10:
是这样你可以试试利用脚本来解决你想做的情况,写出一个循环脚本给出一系列的值,之后进行计算。脚本形式类似与loop.sh优化晶格常数的脚本,如果你的a,b,c三个方向有比例关系容易些,如果不成比例的话,可以尝 ...

"比如a,b不动只优化c方向怕是要改vasp原代码重新编译vasp了。"
这个是怎么实现的,我知道版上有人说过这个,到时具体是怎么修改constr_cell_relax.F文件??
谢谢了~
7楼2011-11-24 18:29:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likai0106

金虫 (小有名气)

引用回帖:
7楼: Originally posted by fzuvivi at 2011-11-24 18:29:36:
"比如a,b不动只优化c方向怕是要改vasp原代码重新编译vasp了。"
这个是怎么实现的,我知道版上有人说过这个,到时具体是怎么修改constr_cell_relax.F文件??
谢谢了~

是的,他后面有提示,你看看提示,这是一个矩阵的形式,可以编译成任何你想优化的方向。
8楼2011-11-25 08:06:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

峰云陈

银虫 (小有名气)

坚持就是胜利 楼主加油
慢慢来..
9楼2013-07-15 08:23:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wang-lf09 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 9/450 2026-08-04 07:46 by k7OM8YghWbkC
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +9 锦衣卫寒战 2026-07-28 13/650 2026-08-04 07:21 by Ermito
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 06:07 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见