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

登录注册

应《网络安全法》要求，自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

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有41人回复
回收溶剂求助已经有7人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

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

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 ：电子结构

引用回帖:

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 个回答

锐利的碎片

木虫 (正式写手)

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的回帖