24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2395  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试324 中药学 一志愿天中医 求调剂 +3 李的Lucky 2026-04-10 3/150 2026-04-15 23:54 by 陈皮皮
[考研] 求助调剂,跨调 +19 X十甫寸Y 2026-04-11 20/1000 2026-04-15 21:18 by cuisz
[考研] 291求调剂 +7 关忆北. 2026-04-14 7/350 2026-04-15 14:39 by 黑科技花岗岩
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 0831生医工第一轮调剂失败求助 +12 小熊睿睿_s 2026-04-11 16/800 2026-04-12 16:28 by 钰璞
[考研] 326求调剂 +6 Shansyn 2026-04-10 6/300 2026-04-12 09:46 by hammer3
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 085600材料与化工329分求调剂 +16 叶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_新宇
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
信息提示
请填处理意见