| 查看: 2392 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
hanyanli0475铜虫 (小有名气)
|
[求助]
石墨烯纳米带总态度与分波态密度不符，分波态密度之和与总态度不相等
|
|
|
大家好，最近我在重复石墨烯纳米带的电子结构，我已经得到了能带和文献的差不多，但是在计算态密度时，发现PDOS之和不等于DOS，我已经把所有得到的每个原子的s和p态的态密度分别加起来了，但是，还是相差的太多，是不是我输入文件又问题啊，请高手们指点一下啊，真是无语啊，往下的工作就卡在这了，谢谢大家啊。下面是输入文件： 1.scf &CONTROL calculation = 'scf' , restart_mode = 'from_scratch' , outdir = '/' , pseudo_dir = '/' , prefix = 'GNRs' , etot_conv_thr = 1.0D-5 , forc_conv_thr = 1.0D-4 , tstress = .true. , tprnfor = .true. , / &SYSTEM ibrav = 8, celldm(1) = 12.85444,celldm(2) = 2.64706,celldm(3) =0.62131, nat = 18,ntyp = 2, ecutwfc = 125.0 ,ecutrho = 500 ,nosym = .false., occupations = 'smearing' ,degauss = 0.02 ,smearing = 'mp' , / &ELECTRONS conv_thr = 1.0d-8 , mixing_beta = 0.7 , mixing_mode = 'plain' , diagonalization = 'david', / ATOMIC_SPECIES C 12.00000 C.pz-vbc.UPF H 1.00000 H.pz-vbc.UPF ATOMIC_POSITIONS crystal C 0.866030546 0.704408575 0.695292364 C 0.866030398 0.300318030 0.695284939 C 0.866028084 0.434371782 0.702120006 C 0.866028204 0.570355384 0.702121677 C 0.866028670 0.637974084 0.202842982 …………………………………… H 0.866033838 0.247278886 0.250636660 K_POINTS automatic 1 1 21 0 0 0 2.nscf &CONTROL calculation = 'nscf' , restart_mode = 'from_scratch' , outdir = './' , pseudo_dir = '。/' , prefix = 'GNRs' , tprnfor= .true. , tstress= .true. , etot_conv_thr= 1.0e-5, forc_conv_thr= 1.0e-4, / &SYSTEM ibrav = 8,celldm(1) = 12.85444,celldm(2) =2.64706,celldm(3) =0.62131, nat = 18,ntyp = 2, ecutwfc = 125.0 ,ecutrho = 500 ,nosym = .false., nbnd=80, / &ELECTRONS conv_thr = 1.0d-8 , mixing_beta = 0.7 , mixing_mode = 'plain' , diagonalization = 'david', / ATOMIC_SPECIES C 12.00000 C.pz-vbc.UPF H 1.00000 H.pz-vbc.UPF ATOMIC_POSITIONS crystal C 0.866030546 0.704408575 0.695292364 C 0.866030398 0.300318030 0.695284939 C 0.866028084 0.434371782 0.702120006 C 0.866028204 0.570355384 0.702121677 ………… H 0.866033669 0.247279466 0.819441547 H 0.866033841 0.757447327 0.250639886 H 0.866033838 0.247278886 0.250636660 K_POINTS automatic 1 1 300 0 0 0 3.dos.x &inputpp outdir='./' prefix='GNRs' fildos='GNRs.dos', Emin=-15.0, Emax=25.0, DeltaE=0.01 / 4.pdos.x &inputpp outdir='./' prefix='GNRs' Emin=-20.0, Emax=25.0, DeltaE=0.01 ngauss=1, degauss=0.02 / 图已给出，请打家帮帮忙呗，谢谢哈！  DOS.jpg  pdos.jpg |
回复此楼» 猜你喜欢
梦想
已经有5人回复
-
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 态密度图分析求助
已经有3人回复
- 总的态密度为什么和分波态密度不相等？
已经有19人回复
- 石墨烯的总态密度与分波态密度
已经有25人回复
- CASTEP几何优化，能带结构和态密度分析
已经有11人回复
- CASTEP对优化后的晶体进行分态密度分析，应该怎么操作
已经有14人回复
- 分波态密度的计算问题？
已经有3人回复
- 分波态密度是什么？如何在CASTEP中做出分波态密度的图？
已经有7人回复
- [求助]关于态密度图的单位和分析
已经有10人回复
- 【求助】总态密度与分波态密度
已经有9人回复
- 【求助】怎么看分波态密度呢
已经有5人回复
- 【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析
已经有30人回复
- 【求助】请帮忙分析能带和态密度图！（红包赠送）
已经有35人回复
- 【分享】关于态密度的一些好的分享
已经有14人回复
- 【求助】MS５.０ Castep 计算不出部分态密度（PDOS）怎么办？（站内有人提问MS4.3，�
已经有19人回复
- 【求助】关于锂原子vasp计算后的1s2s2p分波态密度处理
已经有12人回复
- 【求助】关于castep中的分波态密度
已经有7人回复
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
★ ★
fzx2008: 金币+2, 谢谢指教！ 2012-12-04 20:16:25
fzx2008: 金币+2, 谢谢指教！ 2012-12-04 20:16:25
|
nbnd=30,电子数是60，那么是不可能看到带隙的，所以数据只能显示到你填满的最高能级（-0.6）。带隙的上方是空带，但是你取的nbnd都被占满了，你画出来的只是满带，也就是带隙下方的能带。nbnd再取多一点，哪怕是一条，带隙立马出来。 |
hanyanli04757楼2012-12-04 20:07:18
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
2楼2012-12-03 22:13:59
hanyanli0475
铜虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 173.1
- 散金: 32
- 红花: 3
- 帖子: 296
- 在线: 130.5小时
- 虫号: 1160645
- 注册: 2010-12-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2012-12-03 22:19:14
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
4楼2012-12-04 08:44:08
30