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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Wien2k&FLAPW&ELK » Elk与QE联合使用，是否可行？
111/2返回列表12下一页
查看: 1851  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tanloer

金虫 (小有名气)

[求助] Elk与QE联合使用，是否可行？

最近利用ELK计算一系列TM-Al体系，涉及空位、元素掺杂等。遇到了计算量的问题。
（1）首先，含缺陷体系属于超晶胞，用ELK来进行优化，几乎是不可能完成的任务。于是利用QE，进行VC-relax，得到结构用于ELK计算。这样是否会被杂志接受？为了减少计算量，QE VC-relax时，可否不考虑自旋激化？
（2）其次，在ELK计算中，为了得到结合能，需要计算sublattice的能量，但烦人的是，当对总体系计算时，空带已经足够了，但对TM的sublattice总提示空带不够，如果增加空带，这不就不在一个精度或水平下来进行比较，得到的结合能就不可靠了，不知理解是否正确？
（3）想要得到弹性模量，如果用ELK计算不同体积下的基态能量，然后拟合得到，这个计算量也相当要命。如果只用QE计算的话，感觉分量不够！
各位有何见解？请指教！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-04-15 00:53:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

goldenfisher

金虫 (著名写手)
★ ★
liliangfang: 金币+2, 专家考核, 多谢指教 2012-04-17 09:10:51
如果是磁性体系，不加自旋会可能获得不同的结果。因为有了自旋，体系的能量和离子受力都会不一样。当然，第一步你可以不加自旋优化，在此之后，再加自旋。
此外，结合能肯定都要在同一个标准尺度下进行计算，这样才真的有了比较的标准，要不没法比较。如果你的体系不涉及内壳层电子，也没有必要用ELK。我不明白你说的分量不够是什么意思，这些都主要看其中的物理实现是什么。
赞一下(2人) 回复此楼
6楼2012-04-17 02:54:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 专家考核, 谢谢指教 2012-04-15 08:08:09
tanloer: 金币+5, ★★★很有帮助, 谢谢支持！ 2012-04-15 11:15:18
即使是同一个软件，使用不同的XC泛函或者赝势，都可能收敛的模型尺寸有出入，更何况一个是平面波赝势方法，一个是全势线性缀加波方法的呢……最好是用QE做个大体的结构出来，然后再用ELK仔细查找，并且可以把QE的工作在此省略。
至于只用QE计算，是完全可以的，别说QE这种功能强悍的，就是VASP那样功能比较少的软件，在不挂其他辅助程序的前提下做一个有好的idea的计算，照样也能发nature级的文章啊～
赞一下 回复此楼
思想重于技巧，内涵重于表象
2楼2012-04-15 03:53:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher
★
uuv2010: 金币+1, 多谢指导~ 2012-04-17 15:25:26
没什么问题,elk作者也建议这么干.当然赝势要严格检查一致性,最好是用PAW, abinit有个相关的教程.
赞一下 回复此楼
3楼2012-04-15 08:43:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tanloer

金虫 (小有名气)
引用回帖:
2楼: Originally posted by souledge at 2012-04-15 03:53:11:
即使是同一个软件，使用不同的XC泛函或者赝势，都可能收敛的模型尺寸有出入，更何况一个是平面波赝势方法，一个是全势线性缀加波方法的呢……最好是用QE做个大体的结构出来，然后再用ELK仔细查找，并且可以把QE的 ...

那么，用QE relax时，不考虑自旋计划，妥否？
赞一下 回复此楼
4楼2012-04-15 11:19:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★
uuv2010: 金币+1, 专家考核, 多谢指导~ 2012-04-17 15:25:44
引用回帖:
4楼: Originally posted by tanloer at 2012-04-15 11:19:10:
那么，用QE relax时，不考虑自旋计划，妥否？

没做过磁性体系，不敢乱说。但是个人认为，如果是需要考虑磁性的体系，那还是加上nspin和starting_magnetization(i)稳妥点～
赞一下 回复此楼
思想重于技巧，内涵重于表象
5楼2012-04-15 11:22:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tanloer

金虫 (小有名气)
引用回帖:
6楼: Originally posted by goldenfisher at 2012-04-17 02:54:14:
如果是磁性体系，不加自旋会可能获得不同的结果。因为有了自旋，体系的能量和离子受力都会不一样。当然，第一步你可以不加自旋优化，在此之后，再加自旋。
此外，结合能肯定都要在同一个标准尺度下进行计算，这样 ...

谢谢！的确如你所言，结果刚出来，我就是先没加自旋优化，然后在此基础上加上自旋，得到的结果还是有比较明显的差异。因为FeAl金属间化合物中，B2结构的FeAl实验基态是非磁态，但无论是GGA还是LDA都只能得到铁磁基态，当采用LDA+U后，可以得到非磁基态。文献认为是Fe的3d电子的t2g和eg劈裂时在U的引入后所占比重不同。
我也面临这样一个问题， 当对FeAl单胞选择U=6.33 eV可得到非磁基态，但对3x3x3FeAl时，相同的计算参数，磁性有所增加，约为0.04。难道还要针对这样的超晶胞重新选择U值？在ELK里就会计算量很大。
另外，当计算3x3x3FeAl中心一个Fe空位时，ELK的输入文件是否需要把vkloff做一定的偏移？
赞一下 回复此楼
7楼2012-04-17 11:45:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★
uuv2010: 金币+1, 专家考核, 多谢指导~ 2012-04-17 15:30:28
引用回帖:
7楼: Originally posted by tanloer at 2012-04-17 11:45:58:
谢谢！的确如你所言，结果刚出来，我就是先没加自旋优化，然后在此基础上加上自旋，得到的结果还是有比较明显的差异。因为FeAl金属间化合物中，B2结构的FeAl实验基态是非磁态，但无论是GGA还是LDA都只能得到铁磁 ...

准确的说，确定U值的最初方法就是扩大超胞，达到U值收敛的～
赞一下 回复此楼
思想重于技巧，内涵重于表象
8楼2012-04-17 13:03:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
uuv2010: 金币+1, 专家考核, 欢迎参加讨论 2012-04-17 15:30:45
引用回帖:
8楼: Originally posted by souledge at 2012-04-17 13:03:58:
准确的说，确定U值的最初方法就是扩大超胞，达到U值收敛的～
eb/68/941486_1334639027_102.png

这种办法很早就不适用了，现在有自洽的方法获得U值；不过LDA+U一向是比较复杂的，很多因素需要考虑。合金的体系，更是麻烦。
赞一下 回复此楼
9楼2012-04-17 13:56:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★
uuv2010: 金币+1, 专家考核, 欢迎参加讨论 2012-04-17 15:30:58
引用回帖:
9楼: Originally posted by goldenfisher at 2012-04-17 13:56:24:
这种办法很早就不适用了，现在有自洽的方法获得U值；不过LDA+U一向是比较复杂的，很多因素需要考虑。合金的体系，更是麻烦。

是说Marzari的那个自洽的方法么？那个方法不是很直观啦，并且在LZ的问题的说明上，不如老的超胞法给力啊～
赞一下 回复此楼
思想重于技巧，内涵重于表象
10楼2012-04-17 15:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tanloer 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭