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

登录注册

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

24小时热门版块排行榜

返回列表

minmin_0082003

金虫 (正式写手)

[交流] 【求助】生成赝式时出现“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

回复此楼

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

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

请问如下期刊，哪些期刊对英语语言要求不是很高已经有17人回复
反应溶剂中有水，反应过程有无机盐生成，怎样后处理去掉无机盐和水已经有24人回复
格氏副反应生成烯烃已经有9人回复
in-situ formation 原位生成已经有9人回复
[文件夹加密软件].Railfoxsoft.Folder.Armor.v7.0【转载】已经有15人回复
急问各位虫友，青年基金专家看到的是不是最后生成的pdf 文件，单位生成后乱了。已经有25人回复
VASP中POSCAR文件的生成已经有8人回复
PVDF膜在水中浸泡后，水中会生成一些白色絮状的膜，请教下已经有13人回复
求教NETWORK 4.6 输入文件(RDF)怎么生成的已经有3人回复
配置人工唾液出现絮状沉淀还能用吗？已经有13人回复
【讨论】讨论使用vasp超软赝势参数生成pwscf超软赝势文件的可行性已经有6人回复

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

查看全部散金贴

1楼 2011-01-14 13:34:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★
sunyang1988(金币+1):谢谢指点 2011-01-14 20:29:58
minmin_0082003(金币+4): 2011-01-16 14:06:09

你用的是rel=2，相对论性的，你的轨道给的不对，你应该多看例子，多看input的说明。

赞一下(1人)

回复此楼

2楼2011-01-14 20:17:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

引用回帖:

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

goldenfisher

金虫 (著名写手)

★ ★ ★
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两个，这是至少的。

赞一下(1人)

回复此楼

4楼2011-01-17 20:39:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

引用回帖:

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

goldenfisher

金虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe(金币+3): 谢谢交流。 2011-01-19 20:25:04

lloc指的是你最后用的那个局域轨道的轨道角动量，你这里用的是D，肯定是lloc=2。出现warning最大的可能就是你的轨道截断不对。可以精细地调整。
不过我建议你这里不要用5D轨道做最后一个，或者调整电子占据，使5D轨道有空余轨道出来。

赞一下(2人)

回复此楼

6楼2011-01-19 17:24:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

引用回帖:

Originally posted by goldenfisher at 2011-01-19 17:24:37:
lloc指的是你最后用的那个局域轨道的轨道角动量，你这里用的是D，肯定是lloc=2。出现warning最大的可能就是你的轨道截断不对。可以精细地调整。
不过我建议你这里不要用5D轨道做最后一个，或者调整电子占据，使5 ...

把5d空出来？金原子的最外层本来就是5d106s1.
把d空出来，那只有5d96s2,但是我以前测试过这个，不准确，而且没见过这么用滴。如果理解错误，请多多指教，谢谢

[ Last edited by minmin_0082003 on 2011-1-20 at 11:25 ]

赞一下

回复此楼

7楼2011-01-20 11:24:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖