24小时热门版块排行榜    

查看: 1127  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[交流] 【求助】生成赝式时出现“jjs wrong” 输入文件如下

想生成一个Au的norm-conserving-TM-fully relativistic 的赝式,但是一直报错,求指教
错误信息

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #         0
     from read_psconfig : error #         1
     jjs wrong
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

输入文件
&input
    atom='Au',
    config='[Xe] 4f14 5d10 6s1 6p0',
    iswitch=3,
    rel=2,
    dft='wc'
    rlderiv=3.0, nld=3, eminld=-3.0, emaxld=+2.0
/
&inputp
   lloc=0,
   pseudotype=2,
   tm=true,
   zval = 11,
   file_pseudopw='Au-wc.UPF',
   author='GaoMin'
/
3
6P  2  1  0.00  -0.10 3.30  3.30
5D  3  2 10.00  0.00  1.80  1.80
6S  1  0  1.00  0.00  2.50  2.50
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
youzhizhe(金币+3): 谢谢指导。 2011-01-17 21:16:21
minmin_0082003(金币+6): 2011-01-18 09:30:47
你说的多个轨道在norm conserving里面同样是需要的,这涉及到一些原子物理的知识。jj这里代表的是LS耦合的参数,每个轨道的总角动量不同,需要把轨道量子数和自旋量子数藕合起来。对于rel=2,至少每个轨道要设置两个nl轨道数,其中有很多经验在里面。尤其是对于重金属原子,没有经验的人不要轻易去造,即使造出来也不一定是对的。
jj=l-1/2先用,再用jj=l+1/2. 比如对于P轨道,l=1,则就要用0.5和1.5两个,这是至少的。
4楼2011-01-17 20:39:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

sunyang1988(金币+1):谢谢指点 2011-01-14 20:29:58
minmin_0082003(金币+4): 2011-01-16 14:06:09
你用的是rel=2,相对论性的,你的轨道给的不对,你应该多看例子,多看input的说明。
2楼2011-01-14 20:17:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by goldenfisher at 2011-01-14 20:17:49:
你用的是rel=2,相对论性的,你的轨道给的不对,你应该多看例子,多看input的说明。

谢谢,但是我看例子怎么都跟我理解的对不上,望指点
比如说config “If rel=2, states with jj=l-1/2 are filled first.”
还有card这块
“   * ELSE IF rel = 2 :

      /////////////////////////////////////////
      // Syntax:                             //
      /////////////////////////////////////////

            nwf
            nl(1)    n(1)    l(1)    oc(1)    jj(1)
            nl(2)    n(2)    l(2)    oc(2)    jj(2)
            . . .
            nl(nwf)  n(nwf)  l(nwf)  oc(nwf)  jj(nwf)”  
这个jj(1)跟jj(2) 的区别,设什么都不对,跟angular momenta 有关系,但是例子里面都是一些超软赝式,比如说,同样的5d轨道就给分成好几个5d, 既然都是5d的话,jj值应该一样,但是例子里面就不一样,困惑ing
3楼2011-01-16 14:05:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by goldenfisher at 2011-01-17 20:39:30:
你说的多个轨道在norm conserving里面同样是需要的,这涉及到一些原子物理的知识。jj这里代表的是LS耦合的参数,每个轨道的总角动量不同,需要把轨道量子数和自旋量子数藕合起来。对于rel=2,至少每个轨道要设置两 ...

那如果我不去造的话,也找不到合适的赝式啊。该怎么办捏
还有一个关于产生赝式的问题
赝式的输出文件中有警告的存在,这个也是金的赝式“
warning: wfc  5 2 not converged
warning: wfc  5 1 not converged
warning: wfc  4 3 not converged
warning: wfc  5 2 not converged
warning: wfc  6 0 not converged
warning: wfc  6 1 not converged

输入文件如下
&input
    atom='Au',
    config='[Xe] 4f14 5d10 6s1 6p0',
    iswitch=3,
    rel=1,
    dft='wc'
    rlderiv=3.0, nld=3, eminld=-3.0, emaxld=+2.0
/
&inputp
   lloc=2,
   rcloc=2.5,
   pseudotype=2,
   tm=.true.,
   zval = 11,
   file_pseudopw='Au-wc.UPF',
   author='GaoMin'
/
3
6S  1  0  1.00  0.00  2.50  2.50
6P  2  1  0.00  -0.10 3.30  3.30
5D  3  2 10.00  0.00  1.80  1.80
/
其实我是想用lloc=-2的,无奈总是说输入有错误,所以我选择2了。
请问lloc该如何选择呢
5楼2011-01-18 09:33:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[考研] 317一志愿华南理工电气工程求调剂 +5 Soliloquy_Q 2026-02-28 9/450 2026-03-01 10:15 by Soliloquy_Q
[考研] 材料化工调剂 +3 今夏不夏 2026-03-01 3/150 2026-03-01 10:13 by 杨杨杨紫
[考研] 0856化工专硕求调剂 +4 董boxing 2026-03-01 4/200 2026-03-01 10:11 by 杨杨杨紫
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 284求调剂 +5 天下熯 2026-02-28 5/250 2026-03-01 08:49 by ms629
[考研] 285求调剂 +6 满头大汗的学生 2026-02-28 6/300 2026-03-01 06:29 by Trying]
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考博] 26申博 +4 想申博! 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 311求调剂 +8 南迦720 2026-02-28 8/400 2026-02-28 21:30 by gaoxiaoniuma
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见