24小时热门版块排行榜    

查看: 2333  |  回复: 8

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呢
求指导!谢谢
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注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网站的赝势的,他的很多赝势还是可以使用的
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)。。。囧。。。
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +3 Micky11223 2026-03-25 5/250 2026-03-25 22:32 by barlinike
[考研] 329求调剂 +3 钮恩雪 2026-03-25 3/150 2026-03-25 14:43 by 糖加冰
[考研] 285求调剂 +3 AZMK 2026-03-24 3/150 2026-03-25 12:23 by userper
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 11/550 2026-03-25 10:23 by userper
[考研] 07化学280分求调剂 +7 722865 2026-03-23 7/350 2026-03-25 09:29 by aa331100
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
信息提示
请填处理意见