24小时热门版块排行榜    

查看: 1961  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +4 shxchem 2026-03-20 5/250 2026-03-21 00:24 by 星空星月
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 一志愿中南化学(0703)总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
请填处理意见