24小时热门版块排行榜

>论坛更新日志 (4056)
>虫友互识 (422)
>文献求助 (396)
>导师招生 (373)
>硕博家园 (123)
>考博 (111)
>论文投稿 (109)
>休闲灌水 (104)
>博后之家 (99)
>招聘信息布告栏 (97)
>催化 (59)
>公派出国 (55)
>基金申请 (51)
>考研 (51)
>教师之家 (45)
>论文道贺祈福 (44)

返回列表

fyuewen

金虫 (正式写手)

应助: 30 (小学生)
金币: 1618.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 ]

回复此楼

» 猜你喜欢

2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有4人回复
磺酰氟产物，毕不了业了！已经有6人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有9人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
论文终于录用啦！满足毕业条件了已经有25人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

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

分数占据的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 (小学生)
金币: 1618.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的回帖