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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 手册中计算能带时k点设置的问题
131/2返回列表12下一页
查看: 2918  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiongxiong5712

新虫 (小有名气)

[求助] 手册中计算能带时k点设置的问题

请教各位大侠，siesta3.1手册上第62页关于能带结构设置看不明白：
手册中计算能带时k点设置的问题
这里举的例子是个面心立方晶格的布里渊区k点，但不明白是怎么取的。
我查了fcc的k点坐标：
手册中计算能带时k点设置的问题-1

怎么手册上不是这样取的啊
比如X点的坐标是手册中计算能带时k点设置的问题-2，
手册上怎么是手册中计算能带时k点设置的问题-3？谢谢

[ Last edited by xiongxiong5712 on 2013-7-16 at 15:56 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-07-16 15:48:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

zhangguangping

木虫 (著名写手)
引用回帖:
6楼: Originally posted by xiongxiong5712 at 2013-07-25 10:07:08
手册上的两种单位只是一个2倍关系，这里手册上和我查的坐标不仅仅是2倍关系这么简单啊。您所指的“这两种表示方法采用的单位不同”，是怎么个不同法呢?谢谢您的回复...

下面的图片是我的一个笔记，你可以参考一下。
手册中计算能带时k点设置的问题-4
赞一下(2人) 回复此楼
弘德明志博学笃行
8楼2013-07-25 11:38:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

xiongxiong5712

新虫 (小有名气)
怎么没人理我哪位大侠给点提示啊
赞一下 回复此楼
2楼2013-07-17 05:17:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongxiong5712

新虫 (小有名气)
有哪位好心人给点提示吗？
赞一下 回复此楼
3楼2013-07-17 11:05:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
xiongxiong5712: 金币+5, ★★★很有帮助, 谢谢回复 2013-07-25 10:01:43
手册上的（2.000 0.000 0.000）其实和你的（1,0,0）是等价的点。你说的那个（1/2, 0, 1/2）和（0,1,0）等价。这两种表示方法采用的单位不同。就如同手册上使用的两种单位。
手册上的示例采用的单位是pi/a
赞一下 回复此楼
弘德明志博学笃行
4楼2013-07-24 23:42:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongxiong5712

新虫 (小有名气)
引用回帖:
4楼: Originally posted by zhangguangping at 2013-07-24 23:42:24
手册上的（2.000 0.000 0.000）其实和你的（1,0,0）是等价的点。你说的那个（1/2, 0, 1/2）和（0,1,0）等价。这两种表示方法采用的单位不同。就如同手册上使用的两种单位。
手册上的示例采用的单位是pi/a

非常感谢您的回复。后来我考虑了一下，我说说我的看法，您看对不对。X点等在不同晶胞中表示的点是固定的，比如在面心立方中，X表示四边形的中心，但如果倒易空间的坐标系取法不一样，那这些点的坐标肯定不一样，所以就造成了我看到的手册上的坐标和我查到的坐标不一样。不知道我这样理解对不对？再次感谢您的指教。
赞一下 回复此楼
5楼2013-07-25 10:04:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongxiong5712

新虫 (小有名气)
引用回帖:
4楼: Originally posted by zhangguangping at 2013-07-24 23:42:24
手册上的（2.000 0.000 0.000）其实和你的（1,0,0）是等价的点。你说的那个（1/2, 0, 1/2）和（0,1,0）等价。这两种表示方法采用的单位不同。就如同手册上使用的两种单位。
手册上的示例采用的单位是pi/a

手册上的两种单位只是一个2倍关系，这里手册上和我查的坐标不仅仅是2倍关系这么简单啊。您所指的“这两种表示方法采用的单位不同”，是怎么个不同法呢?谢谢您的回复
赞一下 回复此楼
6楼2013-07-25 10:07:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

引用回帖:
6楼: Originally posted by xiongxiong5712 at 2013-07-25 10:07:08
手册上的两种单位只是一个2倍关系，这里手册上和我查的坐标不仅仅是2倍关系这么简单啊。您所指的“这两种表示方法采用的单位不同”，是怎么个不同法呢?谢谢您的回复...

pi/a是采用Lattice Constant a 来作为单位，而另外一个是采用的倒易格矢作为单位的。另个差别非常大。具体的可以参考一下之前的一个帖子:http://muchong.com/bbs/viewthread.php?tid=4339821. 面心立方的第一布里渊区是具有对称性的，所以(1,0,0)和(0,1,0)以及(0,0,1)都是等价的吧。另外和(-1,0,0)和(0,-1,0)以及(0,0,-1)也等价。如下图所示。
手册中计算能带时k点设置的问题-5
从图中可以看出，X点的坐标是(1/2,0,1/2)这个是用倒易格式作为基矢的。那么具体的坐标可以这样算：1/2*g1+0*g2+1/2*g3=2pi/a(0,1,0),所以用pi/a做单位就是(0,2,0)
赞一下 回复此楼
弘德明志博学笃行
7楼2013-07-25 11:11:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
xiongxiong5712: 金币+5, ★★★★★最佳答案, 非常详细具体，很有耐心，谢谢 2013-07-26 16:13:08
引用回帖:
6楼: Originally posted by xiongxiong5712 at 2013-07-25 10:07:08
手册上的两种单位只是一个2倍关系，这里手册上和我查的坐标不仅仅是2倍关系这么简单啊。您所指的“这两种表示方法采用的单位不同”，是怎么个不同法呢?谢谢您的回复...

SIESTA中计算能带的时候要指定BandLinesScale，也就是说BandLine的单位是什么。

有两个选项：
      pi/a，若是使用这个选项，则需要有LatticeConstant这个标签。其中a是指的LatticeConstant。其中pi/a的单位是1/Bohr（也就说即使你的LattceConstant使用的是Angstrom作单位，这个地方依然会使用1/Bohr作单位）。
     ReciprocalLatticeVectors，若是使用这个选项，不需要有LatticeConstant。同样的，这个选项的单位也是1/Bohr作单位。

下面举一例子来说明这个的用法：
LatticeConstant   5.44 Ang
%block LatticeVectors
0.0   0.5    0.5
0.5   0.0    0.5
0.5   0.5    0.0   
%endblock LatticeVectors

AtomicCoordinatesFormat     ScaledCartesian
%block AtomicCoordinatesAndAtomicSpecies                                                                                                                              
0.00000000   0.00000000    0.00000000   1    Si

0.25               0.25                0.25           1    Si
%endblock AtomicCoordinatesAndAtomicSpecies   

这个是Si晶体的设置，晶格常数使用的是5.44 Angstrom，并且Si晶体是金刚石结构，属面心立方结构。一个原胞中含有两个Si原子。

按照固体物理上，如果设置a为晶体的晶格常数的话，那么面心立方结构的正格子基矢是

a1=a/2(j+k)

a2=a/2(k+i)

a3=a/2(i+j)

上面%block LatticeVectors中就是这么设置的。其中a=5.44 Angstrom.

其倒易格子基矢为：

b1=2π/a(-i+j+k)

b2=2π/a(i-j+k)

b3=2π/a(i+j-k)

因此对于下面倒易空间中O点的计算
%block BandLines
1       0.000   0.000   0.000   \Gamma
100   1.0      0.000    0.000   \O
%endblock BandLines

如果使用的是pi/a选项，那么Gamma点到O点的距离是π*0.529/a=0.305599(1/Bohr),O点的K坐标为：(0.305599,0,0)

如果使用的是ReciprocalLatticeVectors选项，那么那么Gamma点到O点的距离是b1=2*π*1.732*0.529/5.44=1.058626(1/Bohr),O点的K坐标为：0.61068*(-1.0, 1.0,1.0)

按照固体物理上布里渊区高对称点的定义，对于面心立方正格子第一布里渊区常用的高对称点为：

Gamma:2π/a(0,0,0), X:2π/a(1,0,0), K:2π/a(3/4,3/4,3/4), L:2π/a(1/2,1/2,1/2),其中a是正格子的晶格常数。

所以对于上面的X点正确方便的设置是使用pi/a选项，并且设置：

%block BandLines
1       0.000   0.000   0.000   \Gamma
100   2.0      0.000    0.000   \X
%endblock BandLines

如果使用ReciprocalLatticeVectors选项，那么正确的设置应该是：

%block BandLines
1       0.000   0.000   0.000   \Gamma
100   0.0      0.500    0.500   \X
%endblock BandLines
赞一下 回复此楼
弘德明志博学笃行
9楼2013-07-25 11:43:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongxiong5712

新虫 (小有名气)
引用回帖:
9楼: Originally posted by zhangguangping at 2013-07-25 11:43:41
SIESTA中计算能带的时候要指定BandLinesScale，也就是说BandLine的单位是什么。

有两个选项：
      pi/a，若是使用这个选项，则需要有LatticeConstant这个标签。其中a是指的LatticeConstant。其中pi/a的单位是 ...

非常感谢您这么细致的回复，很受益。但是有两个地方不太明白，麻烦再指教一下。您在这里说：
.........
因此对于下面倒易空间中O点的计算
%block BandLines
1       0.000   0.000   0.000   \Gamma
100   1.0      0.000    0.000   \O
%endblock BandLines

如果使用的是pi/a选项，那么Gamma点到O点的距离是π*0.529/a=0.305599(1/Bohr),O点的K坐标为：(0.305599,0,0)

如果使用的是ReciprocalLatticeVectors选项，那么那么Gamma点到O点的距离是b1=2*π*1.732*0.529/5.44=1.058626(1/Bohr),O点的K坐标为：0.61068*(-1.0, 1.0,1.0)
........
这里面不明白为什么如果使用的是pi/a选项，那么Gamma点到O点的距离是π*0.529/a=0.305599(1/Bohr)，这里这个0.529是怎么来的，怎么不直接就是π/a呢？还有后面b1=2*π*1.732*0.529/5.44=1.058626(1/Bohr)，这里的1.732*0.529是怎么来的呢？
不好意思见笑了，还请再次指教。谢谢
赞一下 回复此楼
10楼2013-07-26 16:17:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiongxiong5712 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +6 墨墨漠 2026-03-18 6/300 2026-03-20 00:02 by 23Postgrad
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 288求调剂，一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 271材料工程求调剂 +6 .6lL 2026-03-18 6/300 2026-03-19 15:41 by haoshis
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么！ 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 317求调剂 +3 申子申申 2026-03-19 6/300 2026-03-19 14:16 by 申子申申
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
信息提示
关闭
请填处理意见
关闭