24小时热门版块排行榜    

查看: 3298  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
请填处理意见