应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于晶格常数弛豫，一直有个疑问
131/2返回列表12下一页
查看: 2435  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

甲斐之虎

银虫 (正式写手)

[求助] 关于晶格常数弛豫，一直有个疑问 已有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
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。
1楼 2015-04-28 06:27:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

小虫虫LG

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-04-29 07:37:19
甲斐之虎: 金币+5, ★★★★★最佳答案, 多谢指导 2015-04-29 10:20:20
我感觉石墨烯不能用3吧，不然真空层也跟着变
赞一下(1人) 回复此楼
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 ?@???????
赞一下(1人) 回复此楼
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的回帖
相关版块跳转 我要订阅楼主 甲斐之虎 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 听说今天filecode变了 +25 布布和一二 2026-08-06 48/2400 2026-08-11 11:37 by since—2010
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +5 布布和一二 2026-08-10 6/300 2026-08-11 10:15 by 医学老男孩
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +15 Lanmanbaby 2026-08-09 23/1150 2026-08-11 09:51 by demonstreets
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +6 Tide man 2026-08-10 6/300 2026-08-11 09:45 by tangmnt
[基金申请] 为什么网上很多人说本周 12号出结果 +5 瞬息宇宙 2026-08-10 5/250 2026-08-11 09:29 by 铁帽子农民
[基金申请] 帮忙看看fileCode +6 wwncly 2026-08-10 11/550 2026-08-11 09:14 by xiaruohan
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
关闭
请填处理意见
关闭