[求助] 刚接触QE。做了石墨烯的能带和态密度，不是是否正确，求各位鉴定下

» 猜你喜欢

心脉受损已经有5人回复
博士读完未来一定会好吗已经有15人回复
Springer期刊投稿求助已经有4人回复
读博已经有3人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复

» 本主题相关商家推荐: (我也要在这里推广)

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

石墨烯纳米带总态度与分波态密度不符，分波态密度之和与总态度不相等已经有9人回复
石墨烯的总态密度与分波态密度已经有25人回复

1楼2013-04-17 14:18:39

回帖支持 ( 显示支持度最高的前 50 名 )

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

对了，再把input发上来吧

&control
calculation= ‘scf‘
restart_mode= ‘from_scratch‘
pseudo_dir= ‘/home/wxzhang/soft/espresso-4.2.1/pseudo/
outdir= ‘/home/zyy/tmp/’
prefix= ‘graphene’
tprnfor = .true.
tstress = .true.
&system
ibrav= 4,
celldm(1)=4.01, (*1.8897269), celldm(3)= 10.00,
nat= 2,
ntyp= 1,
ecutwfc= 40.0,
ecutrho= 480,
occupations=‘smearing’ , smearing=‘gaussian’, degauss=0.05，
nbnd= 16
/
&electrons
diagonalization=‘david’
mixing_beta= 0.7
mixing_mode= ‘plain’
conv_thr= 1.0d-10
/
ATOMIC_SPECIES
C 12.0107          C.UPF
ATOMIC_POSITIONS（crystal)
C 2/3 1/3 1/2
C 1/3 2/3 1/2
K_POINTS（automatic)
8 8 8 0 0 0

赞一下(1人)

7楼2013-04-17 15:09:41

普通回帖

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

能看到能带和态密度图么，由于第一次发，不知成功否

2楼2013-04-17 14:57:53

newton3915

铁杆木虫 (正式写手)

应助: 4 (幼儿园)
金币: 4742.3
散金: 40
红花: 4
帖子: 418
在线: 391.3小时
虫号: 609189
注册: 2008-09-22
专业: 凝聚态物性 II ：电子结构

引用回帖:

2楼: Originally posted by 异域神迹 at 2013-04-17 14:57:53
能看到能带和态密度图么，由于第一次发，不知成功否

三个图片都不到哦

3楼2013-04-17 14:58:53

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

不好意思，现在来上传试试

graphene.bands.png

Graphene(1).png

Graphene(2).png

4楼2013-04-17 15:07:37

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

有点儿混乱，因为是第一次用QE，请教下，谢谢了

5楼2013-04-17 15:08:22

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

引用回帖:

3楼: Originally posted by newton3915 at 2013-04-17 14:58:53
三个图片都不到哦...

已上传，请看下，谢谢了。

6楼2013-04-17 15:08:57

newton3915

铁杆木虫 (正式写手)

应助: 4 (幼儿园)
金币: 4742.3
散金: 40
红花: 4
帖子: 418
在线: 391.3小时
虫号: 609189
注册: 2008-09-22
专业: 凝聚态物性 II ：电子结构

引用回帖:

7楼: Originally posted by 异域神迹 at 2013-04-17 15:09:41
对了，再把input发上来吧

&control
calculation= ‘scf‘
restart_mode= ‘from_scratch‘
pseudo_dir= ‘/home/wxzhang/soft/espresso-4.2.1/pseudo/
outdir= ‘/home/zyy/tmp/’ ...

你这个“celldm(1)=4.01, (*1.8897269)” 是什么意思？另外k点真空方向是1就行了，8*8的格子太小了。还有计算能带之前还要做nscf计算的，你的band是不是直接用的scf的结果啊，nscf中要选k点的路径的。

8楼2013-04-17 15:35:32

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

引用回帖:

8楼: Originally posted by newton3915 at 2013-04-17 15:35:32
你这个“celldm(1)=4.01, (*1.8897269)” 是什么意思？另外k点真空方向是1就行了，8*8的格子太小了。还有计算能带之前还要做nscf计算的，你的band是不是直接用的scf的结果啊，nscf中要选k点的路径的。...

不好意思，“celldm(1)=4.01,”，后面括号里是我自己加的换算单位，不在input里的。由于是想先试试而已所以k点的选取比较少。我是直接用的scf的结果的，请问nscf中要选k点的路径的是什么意思啊？

9楼2013-04-17 16:04:28

newton3915

铁杆木虫 (正式写手)

应助: 4 (幼儿园)
金币: 4742.3
散金: 40
红花: 4
帖子: 418
在线: 391.3小时
虫号: 609189
注册: 2008-09-22
专业: 凝聚态物性 II ：电子结构

引用回帖:

9楼: Originally posted by 异域神迹 at 2013-04-17 16:04:28
不好意思，“celldm(1)=4.01,”，后面括号里是我自己加的换算单位，不在input里的。由于是想先试试而已所以k点的选取比较少。我是直接用的scf的结果的，请问nscf中要选k点的路径的是什么意思啊？...

就是你MS算band不是也要用GMKG这些布里渊区的高对称点吗，算nscf的话，KPOINTS就是要写出这些点的坐标，你还是先看看qe的example吧，里面有算能带的例子。

10楼2013-04-17 16:30:06