24小时热门版块排行榜    

查看: 3204  |  回复: 22

souledge

专家顾问 (著名写手)

【答案】应助回帖


zzy870720z(金币+1): 鼓励交流 2011-09-14 17:47:01
改王之王(金币+10): 感谢回复………… 2011-09-14 19:19:13
  分析对了,果然是fhi的赝势本身提供信息太少,所以转化起来很……这封P.Giannozzi的回信估计会比较有用吧:
引用回帖:
shudong wang wrote:

> I want to use fhi2upf.x to convert 27-Co.LDA.fhi to upf format, but I
> have some questions:
> 1.What the wavefunction1-4 should I write?

you should ask this to whoever generated that pseudopotential.
The FHI format does not contain all the info needed to be used
with other codes.

You should in any case know how many valence electrons there
are. Occupancies are used only for starting the calculation
and are not critical. The correct angular momentum is important.
I think it is s,p,d,f, in this order, but no warranty.

> If my converted pp is not precise, how can I solve this? Which
> parameters should I change to enhance the precision??

this is a big question that doesn't have a quick answer.
Please have a look at the atomic_doc/pseudo-gen.pdf file

P.
--
Paolo Giannozzi, Democritos and University of Udine, Italy

思想重于技巧,内涵重于表象
11楼2011-09-14 00:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

引用回帖:
11楼: Originally posted by souledge at 2011-09-14 00:38:10:
  分析对了,果然是fhi的赝势本身提供信息太少,所以转化起来很……这封P.Giannozzi的回信估计会比较有用吧:

谢谢,看来赝势的转换不比生成容易多少,只能先看看atomic_doc/pseudo-gen.pdf慢慢学习了,就像maillist里面说的It is incomplete but it better than nothing.…………呵呵
12楼2011-09-14 19:27:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


uuv2010(金币+1): 欢迎参加交流 2011-09-14 22:07:18
引用回帖:
12楼: Originally posted by 改王之王 at 2011-09-14 19:27:43:
谢谢,看来赝势的转换不比生成容易多少,只能先看看atomic_doc/pseudo-gen.pdf慢慢学习了,就像maillist里面说的It is incomplete but it better than nothing.…………呵呵

实在急用,可以试试这几个模守赝势,分别是LDA和PBE的。质量上实在不敢保证,嗯,尽量高的截断能和K点吧~如果是做一下拉曼之类的定性看看趋势,估计没啥问题。
思想重于技巧,内涵重于表象
13楼2011-09-14 20:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

引用回帖:
13楼: Originally posted by souledge at 2011-09-14 20:20:44:
实在急用,可以试试这几个模守赝势,分别是LDA和PBE的。质量上实在不敢保证,嗯,尽量高的截断能和K点吧~如果是做一下拉曼之类的定性看看趋势,估计没啥问题。

实在太感谢了……我这就试试去…………呵呵
14楼2011-09-14 20:48:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

引用回帖:
13楼: Originally posted by souledge at 2011-09-14 20:20:44:
实在急用,可以试试这几个模守赝势,分别是LDA和PBE的。质量上实在不敢保证,嗯,尽量高的截断能和K点吧~如果是做一下拉曼之类的定性看看趋势,估计没啥问题。

额…………试了试…………算到ph.x的时候提示no elec. filed with metal,不知道是不是赝势的问题,SCF到是能通过…………
15楼2011-09-14 21:55:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


zzy870720z(金币+1): 谢谢提示 2011-10-09 11:23:40
引用回帖:
15楼: Originally posted by 改王之王 at 2011-09-14 21:55:54:
额…………试了试…………算到ph.x的时候提示no elec. filed with metal,不知道是不是赝势的问题,SCF到是能通过…………

你要算有效电荷呢?记得那个是针对电场的线性响应,所以你需要先加个电场~
思想重于技巧,内涵重于表象
16楼2011-09-14 23:02:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

★ ★
zzy870720z(金币+2): 谢谢提供 2011-10-09 11:30:34
做了新的Zn和O的PBE的模守赝势,要不要来试试,是不是比你自己转化的那个好啊,话说,100Ry的截断才能勉强的收敛,实在有点过吧~~~~~
思想重于技巧,内涵重于表象
17楼2011-10-08 17:32:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

感谢提供…………我这就去试试…………呵呵
18楼2011-10-08 19:59:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

★ ★
zzy870720z(金币+2): 鼓励讨论 2011-10-09 11:30:45
引用回帖:
17楼: Originally posted by souledge at 2011-10-08 17:32:39:
做了新的Zn和O的PBE的模守赝势,要不要来试试,是不是比你自己转化的那个好啊,话说,100Ry的截断才能勉强的收敛,实在有点过吧~~~~~

你的赝势截断能设到47ry就收敛的比较好了,和我算的87ry的结果相当,可能是我的赝势转化的有错误,我是按下面的方法做的
Zn
label1, occupancy > 3d 10
label1, occupancy > 4s 2
label1, occupancy > 4p 0
label1, occupancy > 4f 0

O
label1, occupancy > 2s2
label1, occupancy > 2p 4
label1, occupancy > 3s 0
label1, occupancy > 3p 0

然后我去maillist搜了一下看到paolo说的这句话
On Tuesday 14 April 2009 05:19, shangy wrote:

> When I run the fhi2upf.x, it asked the "Wavefunction label and occupancy"
> [...] I tried it like this
> ---------------------
> Wavefunction # 1: label, occupancy > 4s 2
> Wavefunction # 2: label, occupancy > 4p 2
> Wavefunction # 3: label, occupancy > 4d 0
> Wavefunction # 4: label, occupancy > 4f 0
> ---------------------
> but I don't know if it is right?

the FHI format doesn't contain - to the best of my very limited
knowledge - any information on the electronic configuration
used to generate the pseudopotential. Atomic orbitals are
actually used only to produce the starting wavefunctions
and charge density for self-consistency. You should specify
the labels of l=0 to lmax valence orbitals (in the order:s,p,d,f,..)
and a reasonable value for their occupancy. You have Ge with
4 valence electrons, so your choice seems to be correct.


我感觉是不是我Zn的设错了?4s层电子是在3d的内部,
所以Zn应该这样排
label1, occupancy > 4s 2
label1, occupancy > 3d 10
label1, occupancy > 4p 0
label1, occupancy > 4f 0
19楼2011-10-08 22:11:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

昨天试了试这个…………还是不太好…………
label1, occupancy > 4s 2
label1, occupancy > 3d 10
label1, occupancy > 4p 0
label1, occupancy > 4d 0
20楼2011-10-09 17:25:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见