24小时热门版块排行榜    

查看: 2835  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 0703化学求调剂 +5 很老实人 2026-03-09 5/250 2026-03-14 02:57 by JourneyLucky
[考研] 求调剂 +3 清风问长安 2026-03-09 3/150 2026-03-14 02:15 by JourneyLucky
[考研] 一志愿浙江大学0856材料与化工求调剂 +4 yansheng@211 2026-03-09 5/250 2026-03-14 02:10 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 泣血叩求调剂恩,愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 5/250 2026-03-13 22:06 by 星空星月
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 研究生招生 +3 徐海涛11 2026-03-10 7/350 2026-03-12 14:26 by 徐海涛11
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 083000环境科学与工程调剂 +8 mingmingry 2026-03-09 9/450 2026-03-11 10:23 by 沙漠之狐994
信息提示
请填处理意见