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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请教，SIESTA如下Warning该如何消除？谢谢
91/1返回列表
查看: 2653  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

identation

金虫 (著名写手)

[求助] 请教，SIESTA如下Warning该如何消除？谢谢

For Lu, standard SIESTA heuristics set lmxkb to 4
(one more than the basis l, including polarization orbitals).
Use PS.lmax or PS.KBprojectors blocks to override.
Warning: For Lu lmxkb would have been set to 4
Setting it to maximum value of 3 (f projector)

VLOCAL1: 99.0% of the norm of Vloc inside      2.707 Ry
VLOCAL1: 99.9% of the norm of Vloc inside      6.170 Ry
atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge    4.82984
atom: Maximum radius for r*vlocal+2*Zval:    4.31593
GHOST: No ghost state for L =  0
GHOST: No ghost state for L =  1
KBproj: WARNING: Rc( 1)=      4.5943
KBproj: WARNING: Cut of radius for the KB projector too big
KBproj: WARNING: Increasing the tolerance parameter eps
KBproj: WARNING: might be a good idea
GHOST: No ghost state for L =  2
GHOST: No ghost state for L =  3

都是针对Lu原子的 请问这是什么意思呢？ 应该如何消除该warning呢
求指导！谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-11-29 23:49:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emandrei

金虫 (正式写手)

【答案】应助回帖

★ ★
fzx2008(金币+2): 谢谢交流！ 2011-11-30 13:26:26
引用回帖:
1楼: Originally posted by identation at 2011-11-29 23:49:52:
For Lu, standard SIESTA heuristics set lmxkb to 4
(one more than the basis l, including polarization orbitals).
Use PS.lmax or PS.KBprojectors blocks to override.
Warning: For Lu lmxkb would ha ...

1）按照提示，在输入文件里面加上PS.lmax block的参数；
2）重新产生赝势，貌似你的赝势的截断半径太大了。实际上你可以用SIESTSA网站的赝势的，他的很多赝势还是可以使用的
赞一下(1人) 回复此楼
zap65535：用中国字写的。警示作用？版规写得很明白你都不读？我没有义务念给你听。
2楼2011-11-30 09:31:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
2楼: Originally posted by emandrei at 2011-11-30 09:31:52:
1）按照提示，在输入文件里面加上PS.lmax block的参数；
2）重新产生赝势，貌似你的赝势的截断半径太大了。实际上你可以用SIESTSA网站的赝势的，他的很多赝势还是可以使用的

用的就是网站的赝势 我昨天还搜了一下mailing list，貌似稀土的离子有f电子的都会出现这些warning 但是不知道怎么解决
赞一下 回复此楼
3楼2011-11-30 11:31:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vv_china

铜虫 (小有名气)

【答案】应助回帖

引用回帖:
楼: Originally posted by identation at 2011-11-30 11:31:00:
用的就是网站的赝势 我昨天还搜了一下mailing list，貌似稀土的离子有f电子的都会出现这些warning 但是不知道怎么解决

稀土？你是算Lu还是La呀？

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
4楼2011-11-30 11:56:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
4楼: Originally posted by vv_china at 2011-11-30 11:56:28:
稀土？你是算Lu还是La呀？
[ 发自手机版 http://muchong.com/3g ]

Lu
到现在还没解决
SIESTA怎么加快scf收敛呢？跑了300多步都不收敛 好不容易调收敛了 又出现上面的问题了
赞一下 回复此楼
5楼2011-11-30 17:46:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smartly

金虫 (小有名气)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2012-08-14 18:59:44
hi，
我做Si的时候也遇到这个问题。。。生成psf文件时的rc并不大，用lda和gga也没问题，但是vdw计算后，会得到很大的原子rc。
报错信息里有
Check atom: Look for the sentence:
Vna: WARNING: LOCAL NEUTRAL-ATOM PSEUDOPOTENTIAL
Vna: WARNING: Increasing the tolerance parameter EPS
我看到有人说在atom.f里修改eps这个参数就可以算了
但是我没有找到这个atom.f（eps）。。。囧。。。
赞一下(1人) 回复此楼
6楼2012-08-14 15:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoubirong

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by smartly at 2012-08-14 15:22:37
hi，
我做Si的时候也遇到这个问题。。。生成psf文件时的rc并不大，用lda和gga也没问题，但是vdw计算后，会得到很大的原子rc。
报错信息里有
Check atom: Look for the sentence:
Vna: WARNING: LOCAL NEUTRAL-A ...

请问解决了吗？我算Ce也遇到了这样的问题，
赞一下 回复此楼
7楼2018-04-13 11:06:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoubirong

新虫 (初入文坛)
请问楼主，问题解决了吗？我算Ce也遇到了该问题
赞一下 回复此楼
8楼2018-04-13 11:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
8楼: Originally posted by zhoubirong at 2018-04-13 11:08:14
请问楼主，问题解决了吗？我算Ce也遇到了该问题

前两年我有个师弟做过这个Ce的体系我有那么点印象：
1、Ce的基组定义可能有点问题，好像是有轨道未被定义。
2、赝势好像有个空轨道截断半径会很长，不过我觉得应该计算可以不用那个空轨道，测试下来对能量影响不是很大。
3、赝势最好重新生成下，用标量相对论的产生的赝势结果会稍微好点。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
9楼2019-04-02 15:19:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 时间戳又变了8-15 +8 archvillain 2026-08-15 15/750 2026-08-15 16:23 by 医学老男孩
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 15:10 by QXtWNz7PI7wZ
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 7lpolszZVXgi 2026-08-14 3/150 2026-08-15 14:53 by QXtWNz7PI7wZ
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:41 by KxMI1BYBxWX1
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[基金申请] filecode +7 cratir 2026-08-14 11/550 2026-08-15 06:19 by 学员8dgXkO
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭