24小时热门版块排行榜    

查看: 2049  |  回复: 5

fyuewen

金虫 (正式写手)


[交流] 虚晶近似方法得到的赝势怎么用?如何修改相应的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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+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

金虫 (正式写手)


引用回帖:
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的回帖

yanfang2741

禁虫 (文坛精英)

本帖内容被屏蔽

4楼2013-06-25 19:41:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jfgood

木虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
问楼主:虚晶赝势方法生成了一个N和O的混合赝势,这个是怎么实现的呢?是用fhi98PP么?里边的各个量子数怎么设定?
5楼2014-04-29 22:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyuewen

金虫 (正式写手)


★ ★
franch: 金币+2, 鼓励交流。。 2014-05-05 22:51:35
引用回帖:
5楼: Originally posted by jfgood at 2014-04-29 22:03:55
问楼主:虚晶赝势方法生成了一个N和O的混合赝势,这个是怎么实现的呢?是用fhi98PP么?里边的各个量子数怎么设定?

N-O混合,需要使用pwscf里的工具。很久没做这方面工作了。有时间再看看。
fhi98pp这个我没有用过。
6楼2014-04-30 09:26:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fyuewen 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 05:13 by winsaint
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +9 Kittylucky 2026-08-27 9/450 2026-08-28 01:08 by jurkat.1640
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见