24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1603  |  回复: 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

铜虫 (正式写手)

我就是在in0文件里加了一项HYBRID。这个参数在in0说明文件里也有。等我明天去实验室,把in0文件发给你,你看看,你没改这个文件吗?
5楼2009-04-16 21:16:28
已阅   回复此楼   关注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的回帖

shxincui518

铜虫 (正式写手)

请问楼主,我没考虑旋轨耦合,只考虑hybrid(加上-eece选项)case.in0文件中加上HYBRID。case.ineece文件也修改了。但在计算是说我的case.in0文件有问题。不知道楼主遇到过这个问题吗?请赐教!
3楼2009-04-16 09:29:40
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿武汉理工,总分321,英一数二,求老师收留。 +5 nnnnnnn5 2026-03-25 5/250 2026-03-27 04:42 by wxiongid
[考研] 284求调剂 +11 junqihahaha 2026-03-26 12/600 2026-03-27 04:37 by wxiongid
[考研] 复试调剂,一志愿南农083200食品科学与工程 +4 XQTJZ 2026-03-26 4/200 2026-03-27 01:23 by dgnhs
[考研] 0703化学338求调剂! +5 Zuhui0306 2026-03-26 6/300 2026-03-26 23:42 by 381988266
[考研] 一志愿华理,数一英一285求A区调剂 +8 AZMK 2026-03-25 10/500 2026-03-26 22:37 by 学员8dgXkO
[考研] 求调剂 +6 白QF 2026-03-21 6/300 2026-03-26 20:37 by fmesaito
[考研] 07化学280分求调剂 +8 722865 2026-03-23 8/400 2026-03-26 20:00 by 不吃魚的貓
[考研] 286求调剂 +13 Faune 2026-03-21 13/650 2026-03-26 19:52 by peike
[考研] 085602化学工程求调剂。 +4 平乐乐乐 2026-03-26 4/200 2026-03-26 17:57 by fmesaito
[考研] 材料考研求调剂 +3 Dendel 2026-03-23 6/300 2026-03-26 17:51 by fmesaito
[考研] 机械学硕310分,数一英一,一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 求调剂 +3 QiMing7 2026-03-25 3/150 2026-03-25 21:13 by 给你你注意休息
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 招08考数学 +8 laoshidan 2026-03-20 17/850 2026-03-25 17:52 by 一个红太阳
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 276求调剂。有半年电池和半年高分子实习经历 +9 材料学257求调剂 2026-03-23 10/500 2026-03-24 07:36 by wangy0907
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
信息提示
请填处理意见