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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 我就纳闷了，siesta计算graphene为什么总是有gap
121/1返回列表
查看: 2228  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 02:59 by l48hFd8Nwae2
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 21:29 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭