24小时热门版块排行榜    

查看: 1940  |  回复: 10

tanloer

金虫 (小有名气)

[求助] Elk与QE联合使用,是否可行?

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

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

goldenfisher

金虫 (著名写手)

★ ★
liliangfang: 金币+2, 专家考核, 多谢指教 2012-04-17 09:10:51
如果是磁性体系,不加自旋会可能获得不同的结果。因为有了自旋,体系的能量和离子受力都会不一样。当然,第一步你可以不加自旋优化,在此之后,再加自旋。
此外,结合能肯定都要在同一个标准尺度下进行计算,这样才真的有了比较的标准,要不没法比较。如果你的体系不涉及内壳层电子,也没有必要用ELK。我不明白你说的分量不够是什么意思,这些都主要看其中的物理实现是什么。
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 268求调剂 +6 一定有学上- 2026-03-14 7/350 2026-03-17 10:14 by Sammy2
[考研] 材料专硕326求调剂 +5 墨煜姒莘 2026-03-15 5/250 2026-03-16 21:30 by 木瓜膏
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 考研化学学硕调剂,一志愿985 +3 张vvvv 2026-03-15 5/250 2026-03-16 20:25 by 张vvvv
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 材料与化工专硕调剂 +3 heming3743 2026-03-16 3/150 2026-03-16 15:05 by peike
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 298求调剂 +3 Vv呀! 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
信息提示
请填处理意见