24小时热门版块排行榜

返回列表

fyuewen

金虫 (正式写手)

应助: 30 (小学生)
金币: 1620.4
帖子: 480
在线: 125.2小时
虫号: 1506790

[交流] 虚晶近似方法得到的赝势怎么用？如何修改相应的in文件中的原子位置和原子种类？

最近我用虚晶赝势方法生成了一个N和O的混合赝势（比如叫做new.pbe.UPF)，混合比例1:9现在想利用这个赝势计算我之前算过的某体系的band，DOS等。但是问题是输入文件的原子位置，种类等这个地方不知道怎么写笔记恰当。
我的in文件里面原先的数据中包含部分以下内容
ATOMIC_SPECIES
##  #####  ####.UPF    （我用#代替了一些信息）
##  #####  ####.UPF
O  15.9994  O.pbe-van_ak.UPF       这一行，我应该用什么原子来替代呢？

ATOMIC_POSITIONS {crystal}
#    0.000000000 0.000000000 0.0
#    0.500000000 0.500000000 0.5
O       0.500000000 0.000000000 0.5    氧原子的位置也要被虚拟原子代替，如何改？
O       0.000000000 0.500000000 0.5
O       0.500000000 0.500000000 0.0

我想这几个关键地方改好了，之后的计算就可以了。希望高手解答。谢谢

[ Last edited by fyuewen on 2012-8-25 at 10:03 ]

回复此楼

» 猜你喜欢

什么是人一生最重要的？已经有10人回复
面上可以超过30页吧？已经有13人回复
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？已经有10人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有13人回复
版面费该交吗已经有17人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有19人回复
【博士招生】太原理工大学2026化工博士已经有8人回复
280求调剂已经有4人回复

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

分数占据的Cu2Se晶体的热电性质模拟（Nat. Mat+JACS+CM）已经有10人回复
原子间距小于赝势的截断半径时的情况应如何处理呢？已经有6人回复
请问截止能是不是只跟体系中的原子种类和赝势类型有关，与体系原子数无关？已经有5人回复

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

查看全部散金贴

1楼 2012-08-24 16:49:07

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

souledge

专家顾问 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fyuewen: 金币+10, 谢谢souledge，我明天根据您的建议试一下。谢谢您的指导。 2012-08-26 23:27:16
xueht987: 金币+1, 鼓励交流！ 2012-08-27 09:50:34

程序在意的只有赝势中的势，而不在意你给它取个啥名字。只要species中定义的跟position中的对应就OK。
赝势文件中的元素名只在+U计算中用到。
元素原子量则是在分子动力学和声子频率计算中用到。
混合后的赝势容易出现各种计算错误，最典型的是对角化失败，这对赝势质量及两赝势的生成方法要求很多。所以如果没有同作者生成的，最好还是自己生成。
顺便，N和O的化学差别比较大，不敢说这么搞靠不靠谱。至少Vanderbilt建议只用VCA于同族的元素之间。

赞一下

回复此楼

2楼2012-08-26 00:46:57

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

fyuewen

金虫 (正式写手)

应助: 30 (小学生)
金币: 1620.4
帖子: 480
在线: 125.2小时
虫号: 1506790

引用回帖:

2楼: Originally posted by souledge at 2012-08-26 00:46:57
程序在意的只有赝势中的势，而不在意你给它取个啥名字。只要species中定义的跟position中的对应就OK。
赝势文件中的元素名只在+U计算中用到。
元素原子量则是在分子动力学和声子频率计算中用到。
混合后的赝势容 ...

in文件里面的ATOMIC_SPECIES这个地方：
O 15.9994 O.pbe-van_ak.UPF
应该换成哪个原子量呢？
我是按照N:O=1:9这个赝势混合比来写这个原子量吗？我试过15.7979

赞一下

回复此楼

3楼2012-08-26 23:32:10

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