24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2633  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
信息提示
请填处理意见