| 查看: 1602 | 回复: 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 ] |
回复此楼
shxincui518
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 6094.7
- 散金: 155
- 红花: 3
- 帖子: 669
- 在线: 214.3小时
- 虫号: 538822
- 注册: 2008-04-03
- 专业: 凝聚态物性I:结构、力学和
5楼2009-04-16 21:16:28
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
20