版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

zal7891

木虫 (小有名气)

应助: 30 (小学生)
金币: 3485.7
散金: 250
红花: 4
帖子: 252
在线: 326.6小时
虫号: 1252296
注册: 2011-04-01
性别: GG
专业: 凝聚态物性I:结构、力学和

[求助] ELK简单计算失败，求高手解答！已有1人参与

初学ELK软件，测试了一些金属化合物体系，计算还都可以进行。但是最近计算了一下铂系金属氮化物PtN2，发现连简单的自洽都不能进行。化合物的特点是结构中包含N-N这样距离较短的二聚物氮单元，Muffin-tin半径可能对计算有一定影响，所以测试了一下“rgkmax”参数，从7.0到12.0，但是计算仍然不能完成。计算过程中警告很多，如：
Warning(linengy): could not find linearisation energies in s.c. loop
Warning(occupy): minimum eigenvalue less than minimum linearisation energy
Warning(rhonorm): total charge density incorrect for s.c. loop
Warning(occupy): could not find Fermi energy
Error(genlofr): degenerate local-orbital radial functions
在这里求助木虫上的朋友帮忙，大家是否有遇到过类似问题，尤其是对这样的包含N-N二聚物的体系是否有特别要注意的设置，这里先行谢过。
以下是我的ELK输入参数：
tasks
  0

xctype
  20 0 0

rgkmax
  12.0

ngridk
  3  3  3

scale
  1.0

sppath
  '/xxxxx/species/'

avec
  9.197145208    0.000000000    0.000000000
  0.000000000    9.197145208    0.000000000
  0.000000000    0.000000000    9.197145208

atoms
2                                  : nspecies
'Pt.in'                               : spfname
4                                  : natoms; atposl, bfcmt below
0.00000000 0.00000000 0.00000000 0.00000000  0.00000000  0.00000000
0.50000000 0.00000000 0.50000000 0.00000000  0.00000000  0.00000000
0.50000000 0.50000000 0.00000000 0.00000000  0.00000000  0.00000000
0.00000000 0.50000000 0.50000000 0.00000000  0.00000000  0.00000000
'N.in'                               : spfname
8                                  : natoms; atposl, bfcmt below
0.41612397 0.41612397 0.41612397 0.00000000  0.00000000  0.00000000
0.58387603 0.58387603 0.58387603 0.00000000  0.00000000  0.00000000
0.08387603 0.58387603 0.91612397 0.00000000  0.00000000  0.00000000
0.91612397 0.41612397 0.08387603 0.00000000  0.00000000  0.00000000
0.91612397 0.08387603 0.58387603 0.00000000  0.00000000  0.00000000
0.08387603 0.91612397 0.41612397 0.00000000  0.00000000  0.00000000
0.58387603 0.91612397 0.08387603 0.00000000  0.00000000  0.00000000
0.41612397 0.08387603 0.91612397 0.00000000  0.00000000  0.00000000

回复此楼

» 猜你喜欢

建议基金发布提前给出明确的时间点已经有13人回复
2026国自然函评费到账已经有7人回复
今日不放榜？网传国自然预计 8 月 27 日可查结果已经有20人回复
科研孤儿太难了已经有18人回复
看来今天不会放榜了？已经有11人回复
时间戳又变了已经有19人回复
filecode，4个jtjc了已经有14人回复
让我中一个面上吧！已经有13人回复
人气不行了已经有7人回复
感觉是下周放榜了已经有12人回复

» 本主题相关价值贴推荐，对您同样有帮助:

VSAP计算中出现这个提示是什么问题啊，求高手已经有18人回复
求Origin或Matlab高手非线性拟合公式已经有5人回复
求高手相助：origin非线性拟合新建公式怎么编辑已经有8人回复
求matlab高手帮忙画复杂函数的函数图像已经有6人回复
请教如何用VASP计算氧化物的晶格常数？以alpha-Al2O3为例已经有7人回复
求高手解答--高斯软件计算化学反应机理已经有6人回复
向origin高手求助，关于excel和origin做线性拟合已经有7人回复
求Matlab 高手改进程序可付费已经有11人回复
Gaussian 09 计算出来的chk文件无法转化成fchk文件，求高手指教！已经有9人回复
求SWMM软件高手指导已经有15人回复
matlab 气泡边缘提取求助高手！！已经有5人回复
fluent计算发散，仅多了一个烟气回流的结构，求高手指点。已经有5人回复
【求助】求教DPD高手怎样用amorphous cell计算Flory huggins参数已经有8人回复
【求助】CASTEP计算出错，运行失败，请高手帮忙。已经有11人回复
【求助】求高手：怎样用Jade计算晶胞参数已经有4人回复
【原创】DFT计算体会已经有20人回复
【求助成功】如何用VASP计算反铁磁已经有14人回复
【求助】请教各位高手：阻抗值如何计算？已经有5人回复
【求助】计算不收敛怎么办？已经有4人回复
【求助】关于adv mater 08年一篇论文中tddft计算的一些疑问，请高手指正已经有12人回复

我思故我在！

1楼 2014-04-22 16:01:53

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

锐利的碎片

木虫 (正式写手)

star watcher

1ST强帖: 1
应助: 136 (高中生)
金币: 3637.1
散金: 252
红花: 22
帖子: 988
在线: 1224.9小时
虫号: 961933
注册: 2010-03-05
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zal7891: 金币+50, ★★★★★最佳答案, 谢谢回复！ 2014-04-29 19:58:46

看看这个，默认的N参数找不到线性化能量，我换成lapw基组了。
elk.in

CODE:

tasks

  0

stype

  0

swidth

  0.01

lradstp

  2

xctype

  20 0 0

trimvg

  .true.

gmaxvr

  14.0

nrmtscf

  4

autolinengy

  .true.

lmaxapw

  10

lmaxvr

  8

mixtype

  2

rgkmax

  8.0

ngridk

  3 3 3

scale

  1.0

avec

  9.197145208       0.000000000       0.000000000

  0.000000000       9.197145208       0.000000000

  0.000000000       0.000000000       9.197145208

atoms

    2                                    : nspecies

'Pt.in'                                 : spfname

    4                                    : natoms; atposl, bfcmt below

    0.00000000    0.00000000    0.00000000    0.00000000  0.00000000  0.00000000

    0.50000000    0.00000000    0.50000000    0.00000000  0.00000000  0.00000000

    0.50000000    0.50000000    0.00000000    0.00000000  0.00000000  0.00000000

    0.00000000    0.50000000    0.50000000    0.00000000  0.00000000  0.00000000

'N.in'                                  : spfname

    8                                    : natoms; atposl, bfcmt below

    0.41612397    0.41612397    0.41612397    0.00000000  0.00000000  0.00000000

    0.58387603    0.58387603    0.58387603    0.00000000  0.00000000  0.00000000

    0.08387603    0.58387603    0.91612397    0.00000000  0.00000000  0.00000000

    0.91612397    0.41612397    0.08387603    0.00000000  0.00000000  0.00000000

    0.91612397    0.08387603    0.58387603    0.00000000  0.00000000  0.00000000

    0.08387603    0.91612397    0.41612397    0.00000000  0.00000000  0.00000000

    0.58387603    0.91612397    0.08387603    0.00000000  0.00000000  0.00000000

    0.41612397    0.08387603    0.91612397    0.00000000  0.00000000  0.00000000

Pt.in

CODE:

 'Pt'                                       : spsymb

 'platinum'                                 : spname

  -78.0000                                  : spzn

   355605.4398                              : spmass

  0.226455E-06    2.500   40.1191   700    : sprmin, rmt, sprmax, nrmt

  22                                        : spnst

   1   0   1   2.00000    T                 : spn, spl, spk, spocc, spcore

   2   0   1   2.00000    T

   2   1   1   2.00000    T

   2   1   2   4.00000    T

   3   0   1   2.00000    T

   3   1   1   2.00000    T

   3   1   2   4.00000    T

   3   2   2   4.00000    T

   3   2   3   6.00000    T

   4   0   1   2.00000    T

   4   1   1   2.00000    T

   4   1   2   4.00000    T

   4   2   2   4.00000    T

   4   2   3   6.00000    T

   4   3   3   6.00000    F

   4   3   4   8.00000    F

   5   0   1   2.00000    T

   5   1   1   2.00000    F

   5   1   2   4.00000    F

   5   2   2   4.00000    F

   5   2   3   5.00000    F

   6   0   1  1.000000    F

   1                                        : apword

    0.1500   0  F                           : apwe0, apwdm, apwve

   0                                        : nlx

   7                                        : nlorb

   0   2                                    : lorbl, lorbord

    0.1500   0  F                           : lorbe0, lorbdm, lorbve

    0.1500   1  F

   1   2                                    : lorbl, lorbord

    0.1500   0  F                           : lorbe0, lorbdm, lorbve

    0.1500   1  F

   2   2                                    : lorbl, lorbord

    0.1500   0  F                           : lorbe0, lorbdm, lorbve

    0.1500   1  F

   3   2                                    : lorbl, lorbord

    0.1500   0  F                           : lorbe0, lorbdm, lorbve

    0.1500   1  F

   3   3                                    : lorbl, lorbord

    0.1500   0  F                           : lorbe0, lorbdm, lorbve

    0.1500   1  F

   -2.5235   0  T

   1   3                                    : lorbl, lorbord

    0.1500   0  F                           : lorbe0, lorbdm, lorbve

    0.1500   1  F

   -1.9151   0  T

   2   3                                    : lorbl, lorbord

    0.1500   0  F                           : lorbe0, lorbdm, lorbve

    0.1500   1  F

   -0.2143   0  T

N.in

CODE:

 'N'                                        : spsymb

 'nitrogen'                                 : spname

  -7.00000                                  : spzn

   25532.65214                              : spmass

  0.755929E-06    1.3000   34.9411   300    : sprmin, rmt, sprmax, nrmt

   4                                        : spnst

   1   0   1   2.00000    T                 : spn, spl, spk, spocc, spcore

   2   0   1   2.00000    F

   2   1   1  1.000000    F

   2   1   2   2.00000    F

   2                                        : apword

    0.1500   0  F                           : apwe0, apwdm, apwve

    0.1500   1  F                           : apwe0, apwdm, apwve

   0                                        : nlx

   0                                        : nlorb

赞一下

回复此楼

2楼2014-04-27 20:17:34

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

zal7891

木虫 (小有名气)

应助: 30 (小学生)
金币: 3485.7
散金: 250
红花: 4
帖子: 252
在线: 326.6小时
虫号: 1252296
注册: 2011-04-01
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

2楼: Originally posted by 锐利的碎片 at 2014-04-27 20:17:34
看看这个，默认的N参数找不到线性化能量，我换成lapw基组了。
elk.in

tasks
  0

stype
  0

swidth
  0.01

lradstp
  2

xctype
  20 0 0

trimvg
  .true.

gmaxvr
  14.0

nrmtscf
  ...

谢谢锐利的碎片的回复，用你的参数自洽计算可以进行了。非常感谢！
但是不知道具体哪些参数控制lapw基组的使用，是elk.in中的参数autolinengy吗？还是需要修改species文件呢？
目前在尝试ELK的几何优化功能，不知道对这样的体系优化的结果是否会可靠。

赞一下

回复此楼

我思故我在！

3楼2014-04-29 20:10:35

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

锐利的碎片

木虫 (正式写手)

star watcher

1ST强帖: 1
应助: 136 (高中生)
金币: 3637.1
散金: 252
红花: 22
帖子: 988
在线: 1224.9小时
虫号: 961933
注册: 2010-03-05
专业: 凝聚态物性 II ：电子结构

引用回帖:

3楼: Originally posted by zal7891 at 2014-04-29 20:10:35
谢谢锐利的碎片的回复，用你的参数自洽计算可以进行了。非常感谢！
但是不知道具体哪些参数控制lapw基组的使用，是elk.in中的参数autolinengy吗？还是需要修改species文件呢？
目前在尝试ELK的几何优化功能，不知 ...

rgkmax, gmaxvr, lmaxapw, lmaxvr是控制平面波的。species里面设的是local orbital和apw，一般不用改，要看明白可以看http://www.springer.com/west/hom ... 10041-22-86705782-0 这本书，论坛里有下载。

赞一下

回复此楼

4楼2014-05-01 17:14:00

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

大懒猫王浩

金虫 (著名写手)

学渣

应助: 5 (幼儿园)
金币: 315.4
散金: 456
红花: 7
帖子: 1198
在线: 285.2小时
虫号: 3536036
注册: 2014-11-14
性别: GG
专业: 核技术及其应用

想算个二聚体的能量与距离变化曲线，求大佬带一下ELK的使用...

赞一下

回复此楼

一步一步往上走，，，，

5楼2019-09-18 15:46:03

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

相关版块跳转我要订阅楼主 zal7891 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 建议基金发布提前给出明确的时间点 +10	kulium 2026-08-21	13/650	2026-08-21 21:58 by alongwaytogo
[基金申请] 2026国自然函评费到账 +7	羊腰板 2026-08-21	7/350	2026-08-21 21:28 by iwuli
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16	医学老男孩 2026-08-20	20/1000	2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17	我4大白菜 2026-08-20	18/900	2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8	chengyan1220 2026-08-21	11/550	2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13	wuchongjun 2026-08-20	19/950	2026-08-21 17:21 by 紫杉醇
[基金申请] filecode，4个jtjc了 +11	ziyangfang 2026-08-19	14/700	2026-08-21 16:51 by applepetty
[基金申请] 人气不行了 +7	fansofjerry 2026-08-21	7/350	2026-08-21 14:13 by Vivilian
[基金申请] 我面上完蛋了 +7	且听虎啸 2026-08-20	8/400	2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23	yaoyewhu2008 2026-08-20	25/1250	2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖？ +7	CrisMessi 2026-08-18	7/350	2026-08-21 12:15 by weiyin
[基金申请] 今天放榜没戏了吧 +9	yuleib84 2026-08-19	11/550	2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3	archvillain 2026-08-18	3/150	2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护，明天很可能放榜！ +10	zju2000 2026-08-16	11/550	2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18	基诺咪客 2026-08-17	23/1150	2026-08-20 17:19 by Godzela
[基金申请] 重要消息，中午系统在维护 +11	yuleib84 2026-08-18	12/600	2026-08-20 11:09 by xskun
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5	瞬息宇宙 2026-08-15	5/250	2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗？ +14	布布和一二 2026-08-19	15/750	2026-08-19 18:07 by gltch
[基金申请] 今天维护系统维护祝所有人高中 +8	gjjjzhong 2026-08-18	9/450	2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4	且听虎啸 2026-08-17	5/250	2026-08-18 00:49 by 蔡棒棒菂

24小时热门版块排行榜

zal7891

[求助] ELK简单计算失败，求高手解答！ 已有1人参与

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

锐利的碎片

【答案】应助回帖

zal7891

锐利的碎片

大懒猫王浩

[求助] ELK简单计算失败，求高手解答！已有1人参与