24小时热门版块排行榜    

查看: 2094  |  回复: 8

fanny0007

金虫 (小有名气)


[交流] 【求助成功】abinit 计算2X2X2supercell MgO 的DOS (16个原子)


如图我用abinit计算了MgO的DOS图,其中MgO是2X2X2的超晶包(8个Mg8个O原子),但是,次价带和价带怎么就展宽了????求高人解释

下图是2个原子的MgO的DOS,,请看两组图的差距


[ Last edited by fanny0007 on 2010-12-8 at 11:05 ]
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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

fanny0007

金虫 (小有名气)


这是我的输入文件:
ndtset 2
#dataset 1

kptopt 1
nshiftk1 4
shiftk1
0.5 0.5 0.5
0.5 0.0 0.0
0.0 0.5 0.0
0.0 0.0 0.5
ngkpt1 8 8 8
prtden1 1
toldfe1 1.0d-6
iscf1 7

#dataset2DOS
iscf2 -3   


getden2 1
getwfk2 1
prtdos2 2    #2OS 3DOS
dosdeltae 0.00005
natsph2 16
ratsph2 1.448 1.448 1.448 1.448 1.448 1.448 1.448 1.448 2.592 2.592 2.592 2.592 2.592 2.592 2.592 2.592
ngkpt2 16 16 16
tolwfr2 1.0d-16


#Definition of the unit cell

chkprim 0

diemac 3
acell 3*5.955536 Angstrom
angdeg 60 60 60


rprim
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
ntypat 2
znucl 12 8
natom 16
typat 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1
xred
0.2499999999999991   0.2499999999999999   0.2500000000000006
0.7499999999999991   0.2499999999999999   0.2500000000000006
0.2499999999999991   0.7499999999999997   0.2500000000000006
0.7499999999999988   0.7499999999999999   0.2500000000000007
0.2499999999999992   0.2499999999999999   0.7500000000000006
0.7499999999999990   0.2500000000000001   0.7500000000000004
0.2499999999999990   0.7500000000000000   0.7500000000000004
0.7499999999999988   0.7499999999999999   0.7500000000000007
0.0000000000000000   0.0000000000000000   0.0000000000000000
0.4999999999999999   0.0000000000000000   0.0000000000000000
0.0000000000000000   0.4999999999999999   0.0000000000000000
0.4999999999999999   0.4999999999999999   0.0000000000000001
0.0000000000000000   0.0000000000000001   0.4999999999999999
0.4999999999999999   0.0000000000000001   0.4999999999999998
-0.0000000000000001   0.5000000000000001   0.4999999999999999
0.4999999999999998   0.5000000000000000   0.5000000000000001


ecut 24.0
nstep 45
2楼2010-12-07 09:36:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988(金币+1):谢谢交流 2010-12-07 19:36:42
fanny0007(金币+1): 2010-12-07 20:59:26
基矢有问题?
3楼2010-12-07 14:29:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanny0007

金虫 (小有名气)


引用回帖:
Originally posted by fanny0007 at 2010-12-07 09:36:33:
这是我的输入文件:
ndtset 2
#dataset 1

kptopt 1
nshiftk1 4
shiftk1
0.5 0.5 0.5
0.5 0.0 0.0
0.0 0.5 0.0
0.0 0.0 0.5
ngkpt1 8 8 8
prtden1 1
toldfe1 1.0d-6
iscf1 7

#dataset2DOS
i ...

记矢错在哪里?比如说::::
4楼2010-12-07 18:41:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanny0007

金虫 (小有名气)


引用回帖:
Originally posted by fzx2008 at 2010-12-07 14:29:03:
基矢有问题?

比如说……怎么改
5楼2010-12-07 18:42:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z(金币+1):谢谢提示 2010-12-07 21:14:47
ABINIT中angdeg和rprim不同时使用
6楼2010-12-07 21:13:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
fanny0007(金币+1): 2010-12-07 21:33:50
ellsaking(金币+2):谢谢详细解释~ 2010-12-08 09:50:13
fanny0007(金币+1): 2010-12-08 11:03:19
fanny0007(金币+2): 2010-12-08 11:04:32
你用的是MgO原胞然后扩成2x2x2的超单胞。
从分数坐标可以看出来,超单胞你是借助MS构造的,分数坐标是正确的。
基矢你不妨也用MS导出来的。注意基矢的长度。

用了rprim后就不用那个angdeg了。我记得以前两个参数in文件中同时出现的话会有提示。

仅供参考。
7楼2010-12-07 21:19:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanny0007

金虫 (小有名气)



zzy870720z(金币+1):谢谢总结 2010-12-08 12:01:17
引用回帖:
Originally posted by fzx2008 at 2010-12-07 21:19:32:
你用的是MgO原胞然后扩成2x2x2的超单胞。
从分数坐标可以看出来,超单胞你是借助MS构造的,分数坐标是正确的。
基矢你不妨也用MS导出来的。注意基矢的长度。

用了rprim后就不用那个angdeg了。我记得以前两个 ...

谢谢啦,采用了你的建议,我吧rprim去掉,问题迎刃而解。谢谢
8楼2010-12-08 11:04:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)



小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by fanny0007 at 2010-12-07 18:41:52:



记矢错在哪里?比如说::::

fzx2008   说的有道理。rprim和angdeg不能同时用。
9楼2010-12-13 09:44:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fanny0007 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 aOn84t7OD1JN 2026-08-19 4/200 2026-08-19 22:31 by BqkHsOAGISWU
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 PnssoZglMj3v 2026-08-19 4/200 2026-08-19 22:28 by BqkHsOAGISWU
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 22:19 by BqkHsOAGISWU
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 PnssoZglMj3v 2026-08-19 4/200 2026-08-19 22:16 by BqkHsOAGISWU
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 17:42 by oCMhZAsgNFoO
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
信息提示
请填处理意见