24小时热门版块排行榜    

查看: 1964  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 297求调剂 +9 戏精丹丹丹 2026-03-17 9/450 2026-03-21 01:49 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 一志愿西安交通大学 学硕 354求调剂 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见