| 查看: 1548 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 自旋轨道耦合操作 | 自旋轨道耦合 |
» 猜你喜欢
求助一下有机合成大神
已经有3人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有13人回复
-
存款400万可以在学校里躺平吗
已经有30人回复
-
求推荐英文EI期刊
已经有5人回复
-
请教限项目规定
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
shirley8295
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 60.5
- 帖子: 52
- 在线: 33.2小时
- 虫号: 690730
- 注册: 2009-01-10
- 专业: 凝聚态物性 II ：电子结构
2楼2009-04-16 08:53:35
shxincui518
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 6094.7
- 散金: 155
- 红花: 3
- 帖子: 669
- 在线: 214.3小时
- 虫号: 538822
- 注册: 2008-04-03
- 专业: 凝聚态物性I:结构、力学和
3楼2009-04-16 09:29:40
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
4楼2009-04-16 10:21:01
shxincui518
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 6094.7
- 散金: 155
- 红花: 3
- 帖子: 669
- 在线: 214.3小时
- 虫号: 538822
- 注册: 2008-04-03
- 专业: 凝聚态物性I:结构、力学和
5楼2009-04-16 21:16:28
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
结贴
★ ★
zdhlover(金币+2,VIP+0):请反馈重新计算后的结果后返还金币，多谢反馈和交流 4-18 00:57
zdhlover(金币+2,VIP+0):请反馈重新计算后的结果后返还金币，多谢反馈和交流 4-18 00:57
|
版主帮我把剩下的金币返还给我吧，这个帖子可以结了。 在mailing list 上，Peter Blaha给的建议是调整计算的方法，因为自旋耦合计算可能会造成晶格对称性变化，因此第一步的hybrid设置可能与initso后的要求相悖，而case.ineece文件不会被程序自动修改。所以如果我没有手动update case.ineece文件，就有可能出错。 我准备按照他的建议，按下面步骤重新进行计算: runsp initso runsp -so runsp -so -eece |
6楼2009-04-17 20:31:35