24小时热门版块排行榜    

查看: 516  |  回复: 2

mansamsung8113

木虫 (正式写手)

[求助] 关于结构参数的精修

看了侯老师的教程,有一个是关于Mg的参数精修的例子,最后一步教程上说“在该体积下,重复上面第2和3步,得到平衡状态下的其他晶胞参数。这一步就是:在得到了E(V)曲线后,通过拟合状态方程得到平衡状态下的体积,计算出上面脚本中的变量$i的值,并改变$i的循环值,再运行run_cell计算一次,得到其他的结构参数c和位置u.。
我的理解是在第一次运行完run_cell,用murn.x拟合E-V曲线,得到平衡状态的体积V0和晶格常数a0后,然后在把run_cell的 $i 值就取为a0后,再运行run_cell 吧?
第二次运行的run_cell时,对应的脚本开头是不是应该这样:
#!/bin/sh
rm WAVECAR
for i in a0
do
cat > POSCAR < Mg-hex
$i
0.0 -1.0 0.0
0.8660254037844 0.5 0.0
0.0 0.0 1.6230529595
.....................
还请各位大侠赐教
回复此楼

» 猜你喜欢

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

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

mansamsung8113

木虫 (正式写手)

题目说错了,应该是计算Mg的结构参数
坚持不懈
2楼2012-07-28 16:51:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流。 2012-07-30 12:02:12
如果不需要计算弹性模量的话,直接在INCAR里设置ISIF=3和你想要的NSW、IBRION、EDIFFG就行了。
3楼2012-07-29 23:41:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mansamsung8113 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 07:55 by GBJlw5zzHWpZ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 20/1000 2026-11-06 07:55 by GBJlw5zzHWpZ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 ncTVBM6AOJHF 2026-11-04 16/800 2026-11-06 07:55 by GBJlw5zzHWpZ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 05:43 by qL1cU46YIdK0
[基金申请] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 05:43 by qL1cU46YIdK0
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-06 05:22 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 05:21 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 05:15 by Equinoxhua
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ncTVBM6AOJHF 2026-11-04 16/800 2026-11-06 04:58 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 13/650 2026-11-06 04:49 by qL1cU46YIdK0
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 4/200 2026-11-06 02:37 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
信息提示
请填处理意见