24小时热门版块排行榜    

查看: 2113  |  回复: 8
【悬赏金币】回答本帖问题,作者wang-lf09将赠送您 3 个金币

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 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后网上办公系统维护 +7 取名字烦人 2024-06-25 8/400 2024-06-26 09:44 by 没入门滴
[有机交流] 过柱子,产品,杂质在是 pe:ea=100:1 也一起出来? +5 w256 2024-06-25 5/250 2024-06-26 09:31 by 小木木cc
[硕博家园] 博士该不该读? +6 L1009225316 2024-06-25 6/300 2024-06-26 08:41 by 陶渊明2024
[硕博家园] 考博英语 +5 198新用户 2024-06-25 10/500 2024-06-26 08:09 by 地狱伞兵
[考博] 关于读博感觉自己很抓马 +15 小九月 2024-06-19 21/1050 2024-06-25 23:57 by version211
[找工作] 老博士选高校 5+3 woodhead8386 2024-06-24 4/200 2024-06-25 22:16 by wangtaotall
[有机交流] 对苯乙烯磺酰氯的合成机理 25+3 该死的科研 2024-06-24 5/250 2024-06-25 17:30 by 王学士
[基金申请] 焦虑没有毛线用,默默前行是王道! +4 漠上藜梭 2024-06-24 8/400 2024-06-25 14:32 by 漠上藜梭
[硕博家园] 夏至,要不要硕博联谊 +4 我是王小帅 2024-06-21 6/300 2024-06-25 14:25 by 庭前花开月儿明
[硕博家园] 数据不好 +5 Hetai 2024-06-23 7/350 2024-06-25 12:37 by 1591099
[基金申请] 演员杨幂在核心期刊发表论文,中国广播电视学期刊 +9 babu2015 2024-06-22 12/600 2024-06-25 11:05 by 我是勇敢牛牛
[教师之家] 有没有今年的影响因子? +6 jurkat.1640 2024-06-22 10/500 2024-06-24 22:48 by 潇湘之迷
[考博] 申请25博士,可以提前进组做科研助理 +6 逐梦途中w 2024-06-22 8/400 2024-06-24 21:49 by keyaner23
[基金申请] 今年什么时候会评啊 +7 lancet0903 2024-06-24 7/350 2024-06-24 18:48 by Pickfoot
[基金申请] 自然基金变化?? +13 今晚推荐22 2024-06-21 23/1150 2024-06-24 17:20 by 舟言川
[基金申请] b口yq上会情况 +5 XUAAAAA 2024-06-19 6/300 2024-06-24 01:13 by sunrise3148
[基金申请] 请问三个字的中文名字,在文中到底怎么表示? 5+3 zhongyuan764 2024-06-22 6/300 2024-06-23 17:42 by 鱼翔浅底1
[基金申请] 国产期刊影响因子大于12的有20多个个了 +17 babu2015 2024-06-20 22/1100 2024-06-23 13:12 by 德尚中行
[论文投稿] 第一篇论文投稿进程记录 +4 慎独的小花卷 2024-06-20 9/450 2024-06-20 20:37 by 慎独的小花卷
[论文投稿] 投稿求助 +4 平凡的日子 2024-06-19 5/250 2024-06-20 16:24 by yueyueyue@
信息提示
请填处理意见