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 ]

回复此楼

» 猜你喜欢

国自科送审了吗已经有9人回复
博士招生已经有5人回复
夜，静悄悄的已经有8人回复
研究生做的很差，你们会让毕业吗？已经有6人回复
26年博士申请自荐-电催化已经有5人回复
2026博士或科研助理转27年博士已经有5人回复
26年申博自荐-计算机视觉已经有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的回帖