24小时热门版块排行榜    

查看: 1596  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

nkleof

银虫 (正式写手)

[交流] 【求助】自旋轨道耦合计算错误

在计算磁各向异性时,用wien2k做自旋轨道耦合计算,同时考虑hybrid(加上-eece选项)。步骤是先做一般的hybrid计算,完成scf收敛后,initso,在case.inso文件里如下设置:
******************************
WFFIL
4  1  0                      llmax,ipr,kpot
-10.0000   1.50000           emin,emax (output energy window)
   0.  0.  1.                 direction of magnetization (lattice vectors)
******************************
然后运行命令
runsp_lapw -so -dm -eece -it0 -p -i 100 -ec 0.000001 -NI
体系是FeCr2S4,运行过程中在lapw0完成后出错,错误信息如下:
******************************************
**  Error in Parallel LAPW1
**  LAPW1 STOPPED at Sun Apr 12 18:54:28 CST 2009
**  check ERROR FILES!
'SELECT' - no energy limits found for L= 1
'SELECT' - E-bottom   -1.03000   E-top -200.00000
......
******************************************
每个uplapw1_*.error里面都是'SELECT'那两行的错误信息,在wien2k的faq上提到这个问题有两种可能的原因,一个是case.inm 里面的mixing factor 太大,造成太大的电荷波动;另一个是费米能级EF太大,以致case.in1里面的默认E-parameter值0.3不再是一个合适的值。
在我的计算中,grep :FER case.scf的结果显示费米能级最大到1.09,而grep DIS case.scf的结果显示DIS值变化如下
DIS  :  CHARGE DISTANCE       ( 0.0000009 for atom    2 spin 2)      0.0000011
DIS  :  CHARGE DISTANCE       ( 2.7457405 for atom    2 spin 2)      2.7081582
不知道是属于哪种原因引起的。希望高手能给详细讲解一下错误产生的原因,指点个解决方法。
任何建议都欢迎。

[ Last edited by nkleof on 2009-4-15 at 21:13 ]
回复此楼

» 收录本帖的淘帖专辑推荐

自旋轨道耦合操作 自旋轨道耦合

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shxincui518

铜虫 (正式写手)

请问楼主,我没考虑旋轨耦合,只考虑hybrid(加上-eece选项)case.in0文件中加上HYBRID。case.ineece文件也修改了。但在计算是说我的case.in0文件有问题。不知道楼主遇到过这个问题吗?请赐教!
3楼2009-04-16 09:29:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

shirley8295

新虫 (小有名气)


nkleof(金币+1,VIP+0):感谢建议,欢迎讨论 4-16 10:48
这个问题我今天在没有加旋轨耦合的计算中也遇到了,这样的问题最有效的办法是修改RMT,当然这是走投无路的选择。后来我发现建结构的时候原胞单位弄错了,搭出来的模型原子挤了一堆。

[ Last edited by shirley8295 on 2009-4-16 at 08:54 ]
生命在于折腾。
2楼2009-04-16 08:53:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

引用回帖:
Originally posted by shxincui518 at 2009-4-16 09:29:
请问楼主,我没考虑旋轨耦合,只考虑hybrid(加上-eece选项)case.in0文件中加上HYBRID。case.ineece文件也修改了。但在计算是说我的case.in0文件有问题。不知道楼主遇到过这个问题吗?请赐教!

好像没有遇到过in0文件出问题,我不记得做hybrid计算时要修改in0文件,你能详细说一下in0文件里在哪儿加HYBRID吗,并把出错信息列一下吗?
4楼2009-04-16 10:21:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shxincui518

铜虫 (正式写手)

我就是在in0文件里加了一项HYBRID。这个参数在in0说明文件里也有。等我明天去实验室,把in0文件发给你,你看看,你没改这个文件吗?
5楼2009-04-16 21:16:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 求调剂323材料与化工 +4 1124361 2026-03-24 4/200 2026-03-25 11:19 by shulmg
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 11/550 2026-03-25 10:23 by userper
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 招08考数学 +7 laoshidan 2026-03-20 16/800 2026-03-24 17:52 by 乌拉儿山脉
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
信息提示
请填处理意见