应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 我就纳闷了，siesta计算graphene为什么总是有gap
121/1返回列表
查看: 1877  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wring9005

金虫 (著名写手)

[交流] 我就纳闷了，siesta计算graphene为什么总是有gap

SystemName          graphene
SystemLabel         graphene

NumberOfSpecies     1
NumberOfAtoms       2
%block ChemicalSpeciesLabel
1   6   C
%endblock ChemicalSpeciesLabel

PAO.BasisType   split
PAO.BasisSize     DZP
LatticeConstant    2.46000000 Ang
%block LatticeVectors
1.0000000000  0.0000000000   0.0000000000
0.5000000000  0.8660254038   0.0000000000
0.0000000000  0.0000000000  10.0000000000
%endblock LatticeVectors
AtomicCoordinatesFormat  Fractional
%block AtomicCoordinatesAndAtomicSpecies
0.0000000000 0.0000000000 0.00000000 1
0.6666666667 0.6666666667 0.00000000 1
%endblock AtomicCoordinatesAndAtomicSpecies

%block kgrid_Monkhorst_Pack
  8  0   0   0
  0  8   0   0
  0  0   1   0
%endblock kgrid_Monkhorst_Pack

XC.functional       LDA
XC.authors          CA
SpinPolarized        F

MaxSCFIterations        501
UseSaveData          F
DM.MixingWeight         0.1
DM.MixSCF1              F
DM.NumberPulay          15
DM.Tolerance            1.d-6

MeshCutoff           300.0 Ry

SolutionMethod         diagon
OccupationFunction     FD
ElectronicTemperature  10 meV

MD.TypeOfRun              CG     # Default value
MD.NumCGsteps             50     # Default value
MD.MaxCGDispl              0.2000000000     Bohr     # Default value
MD.MaxForceTol             0.1555740000E-03 Ry/Bohr     # Default value

LongOutput              F     # Default value
WriteBands              T
BandLineScale ReciprocalLatticeVectors
%block BandLines
   1  0.0000000000 0.0000000000 0.000 \Gamma
300 -0.3333333333 0.6666666667 0.000 K
300  0.0000000000 0.5000000000 0.000 M
300  0.0000000000 0.0000000000 0.000 \Gamma
%endblock BandLines

WriteMullikenPop        2
MullikenInSCF           T

能带计算结果，如图

未命名.JPG
回复此楼

» 收录本帖的淘帖专辑推荐

siesta问题

» 猜你喜欢

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

查看全部散金贴
1楼 2013-04-08 11:27:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)
★ ★
wring9005(franch代发): 金币+2, 谢谢回帖交流，， 2013-04-08 21:29:05
注意K点的坐标！
回复此楼
3楼2013-04-08 15:36:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sakuraxlp

银虫 (正式写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 谢谢回帖交流，， 2013-04-08 21:29:15
gap energy 在 0- 0.25 eV 之间都算正常。具体看你graphene 的结构啊。。
给你个wikipedia 的链接，希望你那能打开。。。。这里计算了 单排到6排 蜂格结构的band structure, band gap  多数为非零
http://en.wikipedia.org/wiki/File:Cnt_gnrarm_v3.gif
赞一下 回复此楼
4楼2013-04-08 15:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wring9005

金虫 (著名写手)
引用回帖:
3楼: Originally posted by fzx2008 at 2013-04-08 15:36:24
注意K点的坐标！

K点坐标，没有错啊
赞一下 回复此楼
5楼2013-04-08 16:06:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wring9005

金虫 (著名写手)
引用回帖:
4楼: Originally posted by sakuraxlp at 2013-04-08 15:39:18
gap energy 在 0- 0.25 eV 之间都算正常。具体看你graphene 的结构啊。。
给你个wikipedia 的链接，希望你那能打开。。。。这里计算了 单排到6排 蜂格结构的band structure, band gap  多数为非零
http://en.wiki ...

那个属于纳米带，以前也算过，一般都没有这么大gap
不知道这次怎么就算不对了
赞一下 回复此楼
6楼2013-04-08 16:08:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫于

金虫 (正式写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 谢谢回帖交流，， 2013-04-08 21:29:27
你的K点坐标对么，是-0.333 0.6666 0   还是0.333  0.666 0啊
赞一下 回复此楼
7楼2013-04-08 18:24:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
BandLineScale ReciprocalLatticeVectors
这一行写错了，应该是
BandLinesScale ReciprocalLatticeVectors

同时，按照这个晶胞（ab轴夹角60）的话，k点x轴坐标应该是 1/3 2/3。 夹角120才是-1/3 2/3
赞一下 回复此楼
8楼2013-04-08 23:07:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

to_Landau

铁杆木虫 (正式写手)
K点设置有问题。
回复此楼
9楼2013-04-09 09:04:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wring9005

金虫 (著名写手)
引用回帖:
7楼: Originally posted by 小木虫于 at 2013-04-08 18:24:59
你的K点坐标对么，是-0.333 0.6666 0   还是0.333  0.666 0啊

这两个K点分别对应k点和k'点，都是dirac锥点
赞一下 回复此楼
10楼2013-04-09 10:33:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wring9005

金虫 (著名写手)
引用回帖:
8楼: Originally posted by KalaShayminS at 2013-04-08 23:07:26
BandLineScale ReciprocalLatticeVectors
这一行写错了，应该是
BandLinesScale ReciprocalLatticeVectors

同时，按照这个晶胞（ab轴夹角60）的话，k点x轴坐标应该是 1/3 2/3。 夹角120才是-1/3 2/3

非常感谢！
错在BandLinesScale，这样参数就从ReciprocalLatticeVectors变为pi/a
K点值没错
赞一下 回复此楼
11楼2013-04-09 13:17:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wring9005

金虫 (著名写手)
引用回帖:
9楼: Originally posted by to_Landau at 2013-04-09 09:04:36
K点设置有问题。

确实是K点的问题
但K点设置没问题
赞一下 回复此楼
12楼2013-04-09 13:18:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
mlanqiang2楼
2013-04-08 11:51   回复  
相关版块跳转 我要订阅楼主 wring9005 的主题更新
121/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭