24小时热门版块排行榜    

查看: 6261  |  回复: 15

condensed

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
Illusionist: 金币+5, 有帮助, 108个原子,这是折算的晶格常数 2012-05-17 08:35:43
liliangfang: 金币+1, 谢谢交流 2012-05-17 08:54:31
引用回帖:
3楼: Originally posted by Illusionist at 2012-05-14 14:15:57:
打破对称性之后是否需要提高ENCUT?因为我打破之后,发现晶格常数变大了许多·······3.18变到了3.21左右

发现晶格常数变大了许多·······3.18变到了3.21左右
'''''''''''''''''''''''''''''''''''''''''''''''''''''

胞(原子数)这么小(少)啊,可靠么
11楼2012-05-17 00:34:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
fzx2008: 金币+3, 谢谢指教 2012-05-17 14:57:51
引用回帖:
9楼: Originally posted by dxf_dxf_dxf at 2012-05-17 00:07:44:
掺杂的问题,建议把K格点设成只有Gamma点,再用更大的晶格(例如4x4x4,5x5x5,6x6x6等)。
至于收敛性测试,要注意收敛条件是,而不是
在优化几何结构的时候,建议不要使用ISIF=3选项!速度会非常慢。在同样的计算 ...

先在1x1x1体系下用ISIF=7优化得到粗略的晶格常数
然后做ENCUT和KPOINTS的收敛性测试
按满足收敛条件的ENCUT和KPOINTS值,计算不同晶格常数下体系的能量
用murn.f拟合得到平衡态晶格常数。
把1x1x1体系按比例放大,得到2x2x2/3x3x3/4x4x4等体系。
对3x3x3体系:
KPOINTS设成
CODE:
K-Points
0
Monkhorst Pack
1 1 1
0  0  0

然后进行ENCUT收敛性测试,记录不同ENCUT情况下能量和cpu time
然后在4x4x4x体系下做同样的事情,记录不同ENCUT情况下能量和cpu time
再重复至5x5x5,6x6x6,7x7x7
计算时间超过48小时的任务就可以忽略了。
综合考虑性能和精度的影响,确定最后选用的掺杂体系,例如我就选用4x4x4体系。
在此体系下,加入掺杂原子,按照之前收敛性测试的结果设置ENCUT
接下来就开始几何优化了。这里一般要求保持原胞的形状和体积,不过你也可以先优化原胞的体积,一般不要优化原胞的形状。不过这也说不定,例如极高压情况中体系的晶格会改变(例如从体心立方变成六角密堆,所以原胞形状是否改变跟你所做的具体工作有关。)
优化好原胞以后再做离子迟豫。你可以自己写脚本计算,也可以让vasp自己算。
12楼2012-05-17 13:44:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)

【答案】应助回帖


fzx2008: 金币+1, 谢谢指教 2012-05-17 14:58:02
引用回帖:
11楼: Originally posted by condensed at 2012-05-17 00:34:06:
发现晶格常数变大了许多·······3.18变到了3.21左右
'''''''''''''''''''''''''''''''''''''''''''''''''''''

胞(原子数)这么小(少)啊,可靠么

3.18变到3.21,变化达到9%。一般来说,LDA近似和GGA近似计算得到的晶格常数和实验结果的差别在5%左右。
13楼2012-05-17 13:49:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanjiayu0116

金虫 (初入文坛)

引用回帖:
5楼: Originally posted by 贺仪 at 2012-05-14 14:32:59
首先对你的数据:
2)ENCUT
280   4093.144030
300   4093.144030
320   4093.144030
340   4093.144030
表示怀疑

其次对称性没必要打破,设置isif=3就可以了

楼主可不可以把你测试ENCUT的脚本给我一个
14楼2013-03-26 16:11:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

随便注册

银虫 (初入文坛)

LZ,我也想要你的ENCUT测试脚本。当伸手党了,不好意思,非常感谢您!
tom4587@hotmail.com
15楼2013-03-29 08:16:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)

引用回帖:
2楼: Originally posted by xueht987 at 2012-05-14 14:13:29
首先打破对称性是可以的,只不过计算可能会多费点时间;
至于替位掺杂,要考虑一下掺杂原子与被掺杂原子是否是等价的,若果是非等价掺杂要考虑电荷补偿的问题,相关参数是NELECT,具体关于这方面之前本版有很多的相 ...

这个电荷补偿一定要设置吗?本身不就是电中性的吗
希望有一天成为大神
16楼2013-10-21 15:58:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Illusionist 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] 生生命学部会评 +3 wuke100666 2026-07-15 3/150 2026-07-15 17:11 by kingkocxr
[基金申请] 小木虫没落了,除了祈祷帖子,几乎看不到有价值的帖子 +13 botar 2026-07-14 13/650 2026-07-15 17:08 by kingkocxr
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗? 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[基金申请] 生命口C13面上会评时间 +6 luckinging 2026-07-14 9/450 2026-07-15 14:39 by 雨冰共舞
[教师之家] 内心匮乏 +13 水冰月月野兔 2026-07-12 13/650 2026-07-15 14:21 by jiduquegai
[文学芳草园] 你可曾有过这种感受? +3 Jeanya 2026-07-09 6/300 2026-07-15 11:07 by mbtwt
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[论文投稿] 跨出版社商投稿 20+4 倪好520 2026-07-13 4/200 2026-07-15 08:35 by bobvan
[基金申请] 没收到消息,再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +13 sadpencil 2026-07-12 20/1000 2026-07-14 13:12 by 杠就是你对
[基金申请] 明天E口面上会评 +8 布布和一二 2026-07-12 9/450 2026-07-13 22:59 by QTSorange
[论文投稿] MSER送审了还被拒稿 +4 S778950906 2026-07-08 6/300 2026-07-13 16:22 by tfang
[基金申请] 祈祷青基必中 50+6 hadengry 2026-07-09 16/800 2026-07-12 20:31 by 登山虫虫
[考博] 27届辽宁大学应届毕业生申博 +3 可乐微风大海 2026-07-09 3/150 2026-07-10 21:23 by 贪玩的研究僧
[基金申请] 关于如何从代码看上不上会 +17 且听虎啸 2026-07-09 23/1150 2026-07-10 19:51 by sdfapple719
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
信息提示
请填处理意见