24小时热门版块排行榜    

查看: 2141  |  回复: 12

甲斐之虎

银虫 (正式写手)

[求助] 关于晶格常数弛豫,一直有个疑问 已有3人参与

比如,对于石墨烯,假设PBE泛函能量最小那点为晶格常数=2.45A
按照,侯博士那个测试脚本可以进行测试,比如设置i = 2.44 2.45 2.46 2.47 2.48……
但是,如果不那么做,而是一开始建立模型取实验值a=2.46A,然后对石墨烯进行全弛豫(ISIF=3),是否一样可以得到a=2.45A这个结果呢?
另外,如果这个方法可行,那么这个方法对于其他二维材料是否也同样可行呢?
不知道有没有表达清楚,请各位前辈不吝赐教,先行谢过了

附侯博士测试脚本:
#!/bin/sh
rm WAVECAR
for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30
do
cat > POSCAR <<!
Al-fcc
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1
Direct
0.0 0.0 0.0
!
echo / a = $i angstrom /; time vasp
E=‘grep /TOTEN0 OUTCAR | tail -1 | awk -{printf /%12.6f \n0, $5 }.‘
V=‘grep /volume0 OUTCAR | tail -1 | awk -{printf /%12.4f \n0 , $5}.‘
echo $V $E >>EtVo.dat
done
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

别人笑我太疯颠,我笑他人看不穿。不见五陵豪杰墓,无花无酒锄作田。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

小虫虫LG

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-04-29 07:37:19
甲斐之虎: 金币+5, ★★★★★最佳答案, 多谢指导 2015-04-29 10:20:20
我感觉石墨烯不能用3吧,不然真空层也跟着变
4楼2015-04-28 14:05:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢交流 2015-04-29 07:37:27
甲斐之虎: 金币+13, ★★★★★最佳答案, 确实如此,神了…… 2015-04-29 10:22:32
ISIF=3 ???? 4.4 ?r???????


С???LG   ?f????^???_
????????ī??? ISIF=3
??????????????? (????????)
????^???ī??????????С???С (??????С?? 3A ??? 6A?????c????)
??K??????   "?ī?bulk structure"
????????ī? thin film ??

??????S?Y??(????????Y??)???  ??????? ISIF=3 ?@???????
5楼2015-04-28 16:15:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

future_wl

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
ljw4010: 金币+2, 谢谢交流! 2015-04-28 10:10:08
甲斐之虎: 金币+2, 多谢交流 2015-04-29 10:20:59
我觉得侯老师那个脚本是4.6版本时候用的吧,那时候ISIF=3还没有,不能进行全结构弛豫,所以要进行不同体积和能量的拟合,新版本的话一般的体系用ISIF=3就能优化。
个人见解!
未来就是现在
2楼2015-04-28 08:45:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

q965240010

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by future_wl at 2015-04-28 08:45:24
我觉得侯老师那个脚本是4.6版本时候用的吧,那时候ISIF=3还没有,不能进行全结构弛豫,所以要进行不同体积和能量的拟合,新版本的话一般的体系用ISIF=3就能优化。
个人见解!

4.6的时候ISIF=3不能使用吗?
3楼2015-04-28 09:46:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)

引用回帖:
5楼: Originally posted by jpchou at 2015-04-28 16:15:36
ISIF=3 早在 4.4 時就已存在


小虫虫LG   說的比較正確
如果你做石墨烯用 ISIF=3
真空如果取的不夠厚的話 (要很厚很厚很厚)
計算過程石墨烯真空層就會一直變小一直變小 (好像會變小到 3A 還是 6A,有點忘了) ...

请问ISIF=3和侯博士@valenhou001 的方法是否有优劣区别?
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
6楼2015-04-28 18:12:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

引用回帖:
6楼: Originally posted by magicmonk at 2015-04-28 18:12:39
请问ISIF=3和侯博士valenhou001 的方法是否有优劣区别?...

以石墨烯為例的話
不是優劣之分的問題
是 ISIF=3 "不適用"
7楼2015-04-29 04:59:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)

引用回帖:
7楼: Originally posted by jpchou at 2015-04-29 04:59:30
以石墨烯為例的話
不是優劣之分的問題
是 ISIF=3 "不適用"...

哦,那bulk结构的呢?
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
8楼2015-04-29 08:03:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫虫LG

木虫 (正式写手)

引用回帖:
8楼: Originally posted by magicmonk at 2015-04-29 08:03:33
哦,那bulk结构的呢?...

Bulk结构当然是可以用的。
9楼2015-04-29 10:44:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

引用回帖:
5楼: Originally posted by jpchou at 2015-04-28 16:15:36
ISIF=3 早在 4.4 時就已存在


小虫虫LG   說的比較正確
如果你做石墨烯用 ISIF=3
真空如果取的不夠厚的話 (要很厚很厚很厚)
計算過程石墨烯真空層就會一直變小一直變小 (好像會變小到 3A 還是 6A,有點忘了) ...

其实取决于你设置的参数,尤其是收敛条件。
参数设置的好的话,一般情况下10A~15A就够了。虽然会变小一点,但是不会太大。
10楼2015-04-29 16:30:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 甲斐之虎 的主题更新
信息提示
请填处理意见